examp66.m

来自「ADSP TOOLBOX: Version 2.0 and gui m-file」· M 代码 · 共 13 行

clc,echo on 
%EXAMPLE 66 
n=[1 2 0.75];d=[1 27.5 261.5 1039 1668 864]; % H(s)
[nz,dz]=s2zni(n,d,100,'trap')                % H(z) 
dt=numdig(dz,7);                             % Denominator of HT(z)
F=0:0.001:0.1;w=2*pi*F;                      % Frequency array
H=freqz(nz,dz,w);                            % Response of H(z)
HT=freqz(nz,dt,w);                           % Response of HT(z)
plot(F,abs([H;HT]))                          % Overplot
abs([roots(nz) roots(dz) roots(dt)])         % Check for stability
echo off %end of example