3-12.m

来自「是《MATLAB数字信号处理与应用》一书的源代码,该书由李正周编著 清华大学出版」· M 代码 · 共 16 行

%例程3-12 求解单位冲击响应的幅频

 N=20;
 a=[0.4 0.3 -1];
 b=[1 -0.2 0.6 0.3];
 x=[1 zeros(1,N-1)];
 k=0:N-1;
 y=filter(a,b,x);
 Fs=1000;
 w=0:pi/Fs:pi;
 h=freqz(a,b,w);
 subplot(2,1,1);
 stem(k,y);
 subplot(2,1,2);
 plot(w/pi,abs(h));