example9_10.m

原名：数字信号处理基础及MATLAB实现,解压缩密码：[e-mail]

%例9-10,example9_10
%设计条件：w1s=0.2pi;wp=0.35pi;w2s=0.65pi;w2s=0.8pi;Rp=1dB;Ar=60dB;
N=40;alpha=(N-1)/2;L=0:N-1;wL=(2*pi/N)*L;
T1=0.109021;T2=0.59417456;
Hrs=[zeros(1,5),T1,T2,ones(1,7),T2,T1,zeros(1,9),T1,T2,ones(1,7),T2,T1,zeros(1,4)];    %理想滤波器振幅响应抽样
Hdr=[0,0,1,1,0,0];wdl=[0,0.2,0.35,0.65,0.8,1];   %理想滤波器振幅响应
k1=0:floor((N-1)/2);k2=floor((N-1)/2)+1:N-1;
angH=[-alpha*(2*pi)/N*k1,alpha*(2*pi)/N*(N-k2)];
H=Hrs.*exp(j*angH);
h=real(ifft(H,N));
[db,mag,pha,grd,w]=freqz_m(h,1);
[Hr,ww,a,M]=hr_type2(h);
%画图
figure(1)
subplot(2,2,1);plot(wL(1:20)/pi,Hrs(1:20),'o',wdl,Hdr);
axis([0,1,-0.1,1.1]);title('频率样本  N=40')
xlabel('频率（单位：pi）');ylabel('h(n)');
subplot(2,2,2);stem(L,h);axis([-1,N,-0.1,0.3]);
title('脉冲响应');xlabel('n');ylabel('h(n)');
m1=[0 0];m2=[-1 60];line(m2,m1)
subplot(2,2,3);plot(ww/pi,Hr);
axis([0,1,-0.2,1.2]);title('振幅响应')
xlabel('频率（单位：pi）');ylabel('Hr(w)');
subplot(2,2,4);plot(w/pi,db);
axis([0,1,-90,10]);title('幅度响应')
xlabel('频率（单位：pi）');ylabel('分贝');