program_03_07.m

数字信号处理和小波分析的源码

% Program 3_7
% Determination of the Factored Form
% of a Rational z-Transform
%
colordef black;
num = input('Type in the numerator coefficients = ');
den = input('Type in the denominator coefficients = ');
[z,p,k] = tf2zp(num,den);
m = abs(p);
disp('Zeros are at');disp(z);
disp('Poles are at');disp(p);
disp('Gain constant');disp(k);
disp('Radius of poles');disp(m);
sos = zp2sos(z,p,k);
disp('Second-order sections');disp(real(sos));
zplane(num,den,'g')