t2.m

来自「数字信号matlab程序,是科学出版社编的书的上机作业」· M 代码 · 共 12 行

wp=pi/3;N=61;
k=0:1:60;
a=-60*pi*k/61;
Hk=[ones(1,11) zeros(1,39) ones(1,11)];
Hd=Hk.*exp(j*a); 
hn=ifft(Hd,N)
figure(1)
n=0:1:60;
stem(n,hn,'.');grid on
figure(2)
[H,w]=freqz(hn);            
plot(w,20*log(abs(H))); grid on