signal.m

很多matlab的源代码

b=[1,0];a=[1,-0.9];
[H,w]=freqz(b,a,100);
magH=abs(H);phaH=angle(H);
subplot(2,1,1);
plot(w/pi,magH);
grid;
xlabel('频率 单位：pi')
ylabel('幅度')
title('幅度响应');
subplot(2,1,2);
plot(w/pi,phaH);grid
xlabel('频率 单位：pi')
ylabel('相位 单位：pi')
title('相位响应');