s43.m

数字信号处理的MATLAB仿真加实验指导

%a26.m
M=21;alpha=(M-1)/2;l=0:M-1;wl=(2*pi/M)*l;
Hrs=[1,1,1,1,zeros(1,13),1,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+0.5),alpha*(2*pi)/M*(M-k2-0.5)];
H=Hrs.*exp(j*angH);
h1=ifft(H,M);
n=0:1:M-1;
h=real(h1.*exp(j*pi*n/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+1/M,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+1/M,Hrs(1:11),'.');
axis([0 1 -0.2 1.3]);title('振幅响应H(w)');xlabel('频率（单位:pi）');
subplot(2,2,4);plot(w/pi,db);title('幅度响应（单位:dB）');
axis([0 1 -60 10]);xlabel('频率（单位:pi）');