p6_1.m

数字信号处理源代码

% Program P6_1
% Conversion of a rational transfer function
% to its factored form
num = input('Numerator coefficient vector = ');
den = input('Denominator coefficient vector = ');
[z,p,k] = tf2zp(num,den);
sos = zp2sos(z,p,k)