c122alt.m

来自「matlab例程MATLAB 6.x版的内建数据类型（Built-in data」· M 代码 · 共 24 行

a=0.83

num=[1 0 -1]; den=[1 -sqrt(2)*a a^2];
W=-pi:pi/100:pi;
H=freqz(num,den,W);
mag1=abs(H);
phase1=180/pi*unwrap(angle(H));

a=0.9
num=[1 0 -1]; den=[1 -sqrt(2)*a a^2];
H=freqz(num,den,W);
mag2=abs(H);
phase2=180/pi*unwrap(angle(H));

a=0.96
num=[1 0 -1]; den=[1 -sqrt(2)*a a^2];
H=freqz(num,den,W);
mag3=abs(H);
phase3=180/pi*unwrap(angle(H));

subplot(2,1,1);
plot(W,mag1,W,mag2,W,mag3);
subplot(2,1,2);
plot(W,phase1,W,phase2,W,phase3);