数字信号处理project2.doc

数字信号处理Project21.实验题目Design a linear-phase FIR lowpass filter using the window design technique.The specifications are given as follows:stopband edge: 5 rad/sec passband edge: 3 rad/sec maximum passband attenuation: 0.3dBminimum stopbandband attenuation: 40dBsampling frequency: 15 rad/sec(a). Use each of the following windows for the design: Hamming, Hann, Blackman, and Kaiser. Show the ideal impulse response, the actual impulse response, and plot the gain response of the designed filters for each case. Give the main design steps. Comment on your results. Do not use the function fir1 of Matlab.(b). Repeat (a) except for the ideal impulse response, by using the function fir1 of Matlab. Compare results of (a) with those of (b).2. 实验设计(a) 题：1.确定通带截止频率由于题中给定的频率为实际频率，所以先对给定频率进行归一化。，通带的截止频率为：2.确定不同窗下的长度M，N归一化的过渡带宽为对于Hann、Hamming、Blackman窗用下面的函数计算M，N。c通过查表得,分别为3.11,3.32,5.56。对于Kaiser窗，M、N用下面的公式计算，其中为最小阻带衰减为：40dB。3.计算理想低通的脉冲响应4.使用窗函数产生窗系数winHamming:，Hann:，Blackman:，Kaiser:因为40 0;这样可以求得当n=0时原函数值。(也可以将n加上eps %-系统精度，使n最大程度接近0而不等0） ht=hd.*win; display(,num2str(ht); %打印显示冲激响应系数 %wvtool(ht); subplot(1,3,1); plot(n,hd,.-) %画理想冲激响应系数图 title(hd_labels(filter_kind); xlabel(n,FontSize,12); ylabel(hd,fontsize,12); grid on subplot(1,3,2); plot(n,ht,.-) %画冲激响应系数图 title(ht_labels(filter_kind); xlabel(n,FontSize,12); ylabel(ht,fontsize,12); grid on h,w=freqz(ht,1,512); W=w/pi; H=20*log10(abs(h); subplot(1,3,3); %画增益响应 hold on title(freq_labels(filter_kind); plot(W,H); xlabel(pi (omega/omega,s); ylabel(增益(dB); grid on %为了更直观的看到所设计FIR是否达到指标，可分别在图中标出相关点 %从而可以对不同窗函数进行对比。标示过程如下： dotp=round(mean(find(wwp-0.1&wws-0.1&wwp-0.1&wws-0.1&wws+0.1);%找到点ws的坐标 plot(W(dotp),H(dotp),.r,MarkerSize,25);%对应处打上显目的红点 plot(W(dots),H(dots),.r,MarkerSize,25); text(W(dotp),H(dotp),num2str(round(W(dotp)*1000)/1000),. ,num2str(round(H(dotp)*1000)/1000),dB,FontSize,15);%标示wp在频响图%中的位置，从中可以%看来相应增益 text(W(dots),H(dots),num2str(r