ex041101.m

来自「stft(含:短时傅立叶变换 Gabor变换等)」· M 代码 · 共 11 行

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('frequency in units of pi'); ylabel('|H|')
title('Magnitude Part')
subplot(2,1,2); plot(w/pi,phaH)/pi;grid
%axis([-2,2,-1,1])
xlabel('frequency in units of pi'); ylabel('radians/pi')
title('Angle Part')