设某理想带阻滤波器的频率响应为 试用窗函数设计长度为55的带阻滤波器 满足阻带衰减60dB。请帮忙

正确答案：程序清单如下： wlp=0．4*piwhp：0．6*pi； N=55：n=0：N一1 WC=[wlp／pi一6／Nwhp／pi+6／N； hnl=fir1(N—1WC′stop′boxcar(N))； hn2=fir1(N一1WC′stop′banning(N))； hn3=fir1(N一1WC′stop′blackman(N))； subplot(321)； stem(nhn′．′)； title(′h(n)′)； xlabel(′n′) hw1=fft(hn1512)； hw2=fft(hn2512)； hw3=fft(hn3512)； W=[0：511*2／512； subplot(322)； plot(w20*logl0(abs(hw1)))；xlabe1(′w／pi′) subplot(324)；plot(w20*logl0(abs(hw2)))；xlabe1(′w／pi′) subplot(326)； plot(w20*logl0(abs(hw3)))； xlabe1(′w／pi′) 程序运行结果如图7．8所示。 根据题意当N=55时阻带的衰减为55dB设计时选用了3个窗函数进行带阻滤波器的设计从程序的运行结果可以看出只有选用blackman窗时其阻带最小衰减可以达到＞60dB。
程序清单如下：wlp=0．4*pi,whp：0．6*pi；N=55：n=0：N一1WC=[wlp／pi一6／N,whp／pi+6／N；hnl=fir1(N—1,WC,′stop′,boxcar(N))；hn2=fir1(N一1,WC,′stop′,banning(N))；hn3=fir1(N一1,WC,′stop′,blackman(N))；subplot(3,2,1)；stem(n,hn,′．′)；title(′h(n)′)；xlabel(′n′)hw1=fft(hn1,512)；hw2=fft(hn2,512)；hw3=fft(hn3,512)；W=[0：511*2／512；subplot(3,2,2)；plot(w,20*logl0(abs(hw1)))；xlabe1(′w／pi′)subplot(3,2,4)；plot(w,20*logl0(abs(hw2)))；xlabe1(′w／pi′)subplot(3,2,6)；plot(w,20*logl0(abs(hw3)))；xlabe1(′w／pi′)程序运行结果如图7．8所示。根据题意,当N=55时,阻带的衰减为55dB,设计时选用了3个窗函数进行带阻滤波器的设计,从程序的运行结果可以看出,只有选用blackman窗时,其阻带最小衰减可以达到＞60dB。