examp64.m

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

clc,echo on 
%EXAMPLE 64 
S=2;                                   % Sampling frequency  
[N1, D1]=s2zmatch(1,[1 1],S,'0',0);    % Matched z-transform
[N2, D2]=s2zinvar(1,[1 1],S,'i',0);    % Impulse invariance
[N3, D3]=s2zni(1,[1 1],S,'bili',0);    % Bilinear transform
F=0:0.005:0.5;w=2*pi*F;                % Frequency array
H1=freqz(N1,D1,w);                     % Frequency response of H1(z)
H2=freqz(N2,D2,w);                     % Frequency response of H2(z)
H3=freqz(N3,D3,w);                     % Frequency response of H3(z)
HA=freqs([0 1],[1 1],w*S);             % Frequency response of H(s)
f=F*S;plot(f,abs([H1;H2;H3;HA])),grid  % Overplot against analog frequency
echo off %end of example