s42.m

来自「数字信号处理的MATLAB仿真加实验指导」· M 代码 · 共 19 行

%a25.m
M=21;alpha=(M-1)/2;l=0:M-1;wl=(2*pi/M)*l;
Hrs=[1,1,1,1,zeros(1,14),1,1,1];
Hdr=[1,1,0,0];wdl=[0,0.35,0.35,1];
k1=0:floor((M-1)/2);k2=floor((M-1)/2)+1:M-1;
angH=[-alpha*(2*pi)/M*k1,alpha*(2*pi)/M*(M-k2)];
H=Hrs.*exp(j*angH);
h=real(ifft(H,M));
[db,mag,pha,w]=freqz_m(h,[1]);
[Hr,ww,a,L]=hr_type1(h);
subplot(2,2,1);plot(wl(1:11)/pi,Hrs(1:11),'.',wdl,Hdr);
axis([0 1 -0.1 1.1]);title('频率样本H（k）：M=21')
subplot(2,2,2);stem(l,h);
axis([-1 M -0.1 0.4]);title('脉冲响应h(n)');xlabel('');
subplot(2,2,3);plot(ww/pi,Hr,wl(1:11)/pi,Hrs(1:11),'.');
axis([0 1 -0.2 1.2]);title('振幅响应H(w)');xlabel('频率单位:pi');
subplot(2,2,4);plot(w/pi,db);title('幅度响应（单位:dB)');
axis([0 1 -60 10]);xlabel('频率单位:pi');