ex3011.m

来自「the file contain many matlab signal sour」· M 代码 · 共 17 行

b=[0.0181,0.0543,0.0543,0.0181];
a=[1.0000,-1.7600,1.1829,-0.2781];
m=0:length(b)-1;
l=0:length(a)-1;
K=500;
k=1:1:K;
w=pi*k/K;
num=b*exp(-j*m'*w);
den=a*exp(-j*l'*w);
H=num./den;
magH=abs(H);
angH=angle(H);
subplot(211)
plot(w/pi,magH)
subplot(212)
plot(w/pi,angH/pi*180)