数字信号处理1

《数字信号处理》课程实践报告频率采样法设计FIR数字滤波器频率采样法设计FIR数字滤波器一、实验原理：设所希望得到的滤波器的理想频率响应为，那么由系统函数的内插公式为：（10-1）其中是频率采样值（10-2）为我们提供了一条从频域设计逼近的途径。即令（10-3）则对的逼近为（10-4）至少在采样点的频率上，两者可以具有相同的频响，即（10-5）当需要设计线性相位FIR滤波器时，还必须注意，采样值的幅度和相位一定要遵循线性相位滤波器幅度与相位的四种不同的约束关系。设计所得的频率响应逼近于理想频率响应的程度与理想特性有关，如果越平缓，则越逼近；反之，如果变化越剧烈，采样点之间的理想特性变化大，则内插值与理想值的误差就越大，因而在理想特性的每一个不连续点附近都会出现肩峰与起伏，不连续性越大，出现的肩峰和起伏也越大。为了解决这一问题，在理想特性不连续点的边缘加过渡的采样点，这样虽然加宽了过滤带，但缓和了边缘上两采样点之间的突变。因而将有效地减少起伏振荡，提高阻带的最小衰减。因此频率采样法的设计过程为：给定为了检验所设计的对的逼近性能，我们可以如下做检验：在做FT中，我们可以取点数大于N的DFT代替之。

二、.设计步骤（1）由给定的理想滤波器给出和。（2）由式QUOTE，求得。（3）根据求得和。三、实验内容1.采用频率采样设计法设计FIR数字低通滤波器，满足以下指标取N=20，过渡带没有样本。取N=40，过渡带有一个样本，T=0.39。取N=60，过渡带有两个样本，T1=0.5925，T2=0.1009。分别讨论采用上述方法设计的数字低通滤波器是否能满足给定的技术指标。实验代码与实验结果：（1）N=20过渡带没有样本matlab程序：N=20;alpha=(N-1)/2;l=0:N-1;wl=(2*pi/N)*l;Hrs=[1,1,1,zeros(1,15),1,1];*对理想幅度函数取样得到取样样本Hdr=[1,1,0,0];wdl=[0,0.25,0.25,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(k)h=ifft(H,N);*计算h(n)w=[0:500]*pi/500;H=freqz(h,1,w);*计算幅度响应[Hr,wr]=zerophase(h);*计算幅度函数subplot(221);plot(wdl,Hdr,wl(1:11)/pi,Hrs(1:11),'o');axis([0,1,-0.1,1.1]);xlabel('\omega(\pi)');ylabel('Hr(k)');subplot(222);stem(l,h,'filled');axis([0,N-1,-0.1,0.3]);xlabel('n');ylabel('h(n)');subplot(223);plot(wr/pi,Hr,wl(1:11)/pi,Hrs(1:11),'o');axis([0,1,-0.2,1.2]);xlabel('\omega(\pi)');ylabel('Hr(w)');subplot(224);plot(w/pi,20*log10((abs(H)/max(abs(H)))));axis([0,1,-50,5]);grid;xlabel('\omega(\pi)');ylabel('dB');（2）N=40过渡带有一个样本，T=0.39Matlab程序：N=40;alpha=(N-1)/2;l=0:N-1;wl=(2*pi/N)*l;Hrs=[1,1,1,1,1,0.39,zeros(1,29),0.39,1,1,1,1];*设置过渡带样本Hdr=[1,1,0.39,0,0];wdl=[0,0.2,0.25,0.3,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=ifft(H,N);w=[0:500]*pi/500;H=freqz(h,1,w);[Hr,wr]=zerophase(h);subplot(221);plot(wdl,Hdr,wl(1:21)/pi,Hrs(1:21),'o');axis([0,1,-0.1,1.1]);xlabel('\omega(\pi)');ylabel('Hr(k)');subplot(222);stem(l,h,'filled');axis([0,N-1,-0.1,0.3]);xlabel('n');ylabel('h(n)');subplot(223);plot(wr/pi,Hr,wl(1:21)/pi,Hrs(1:21),'o');axis([0,1,-0.2,1.2]);xlabel('\omega(\pi)');ylabel('Hr(w)');subplot(224);plot(w/pi,20*log10((abs(H)/max(abs(H)))))axis([0,1,-80,5]);grid;xlabel('\omega(\pi)');ylabel('dB');（3）N=60过渡带有两个样本T1=0.5925，T2=0.1009Matlab程序：N=60;alpha=(N-1)/2;l=0:N-1;wl=(2*pi/N)*l;Hrs=[1,1,1,1,1,1,1,0.5925,0.1099,zeros(1,43),0.1099,0.5925,1,1,1,1,1,1];*设置过渡带样本Hdr=[1,1,0.5925,0.1099,0,0];wdl=[0,0.2,0.2+1/30,0.3-1/30,0.3,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=ifft(H,N);w=[0:500]*pi/500;H=freqz(h,1,w);[Hr,wr]=zerophase(h);subplot(221);plot(wdl,Hdr,wl(1:31)/pi,Hrs(1:31),'o');axis([0,1,-0.1,1.1]);xlabel('\omega(\pi)');ylabel('Hr(k)');subplot(222);stem(l,h,'filled');axis([0,N-1,-0.1,0.3]);xlabel('n');ylabel('h(n)');subplot(223);plot(wr/pi,Hr,wl(1:31)/pi,Hrs(1:31),'o');axis([0,1,-0.2,1.2]);xlabel('\omega(\pi)');ylabel('Hr(w)');subplot(224);plot(w/pi,20*log10((abs(H)/max(abs(H)))))axis([0,1,-120,5]);grid;xlabel('\omega(\pi)');ylabel('dB');（4）由实验结果第四个图可知，当时，阻带增益都没有达到-50dB,阻带增益有所减低，所以设计结果不能满足最初的设计要求。2.采用频率采样技术设计下面的高通滤波器对于高通滤波器，N必须为奇数。选择N=33，过渡带有两个样本，过渡带样本最优值为T1=0.1095，T2=0.598。Matlab程序：N=33;alpha=(N-1)/2;l=0:N-1;wl=(2*pi/N)*l;Hrs=[zeros(1,11),0.1095,0.598,ones(1,7),0.598,0.1095,zeros(1,11)];Hdr=[0,0,0.1095,0.598,1,1];wdl=[0,20/33,22/33,24/33,26/33,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=ifft(H,N);w=[0:500]*pi/500;H=freqz(h,1,w);h=real(h);[Hr,wr]=zerophase(h);subplot(221);plot(wdl,Hdr,wl(1:17)/pi,Hrs(1:17),'o');axis([0,1,-0.1,1.1]);xlabel('\omega(\pi)');ylabel('Hr(k)');subplot(222);stem(l,h,'filled');axis([0,N-1,-0.1,0.3]);xlabel('n');ylabel('h(n)');subplot(223);plot(wr/pi,Hr,wl(1:17)/pi,Hrs(1:17),'o');axis([0,1,-0.2,1.2]);xlabel('\omega(\pi)');ylabel('Hr(w)');subplot(224);plot(w/pi,20*log10((abs(H)/max(abs(H)))));axis([0,1,-1