北交大DSP研究性学习数字滤波器设计2014要点.docxVIP

《数字信号处理》课程研究性学习报告 姓名 学号 同组成员 指导教师 时间 2014年5月18日星期日 数字滤波器设计专题研讨 【目的】 (1) 掌握IIR和FIR数字滤波器的设计方法及各自的特点。 (2) 掌握各种窗函数的时频特性及对滤波器设计的影响。 (3) 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。 【研讨题目】 基本题 分析矩形窗、汉纳窗、哈明窗、布莱克曼窗、凯泽窗的频域特性，并进行比较。 【题目分析】 本题分析不同的窗函数的频域特性，预计可以看出不同的窗有不同的过渡带大小和不同的旁瓣宽度，可以满足对不同的设计设计要求。 【仿真结果】 【结果分析】 各种窗有何特点? 计算过程中，不同的窗都采用了相同的长度（均为10），进行fft计算的长度均为512点，从以上结果可以看出，矩形窗的幅度最大，主瓣最窄，旁瓣幅度大；其余几个窗函数的主瓣幅度均小于矩形窗的主瓣幅度，主瓣宽度约为矩形窗主瓣宽度的2倍多；其中汉纳窗和哈明窗的旁瓣有较小幅度，而布莱克曼窗和凯泽窗旁瓣基本为0。 【阅读文献】 [1]数字信号处理.陈后金 【发现问题】 (专题研讨或相关知识点学习中发现的问题)： matlab提供的函数产生的窗为N行1列的列向量，不是1行N列的一行数据，计算时必须将其进行转置。 使用fft计算的结果是原序列周期化后求出的频谱的主值区间，必须使用fftshift将其整理，然后通过向左平移L/2才能求出原序列的频谱。 【问题探究】 在谱分析中如何选择窗函数，在滤波器设计中如何选择窗函数？ 在谱分析中窗函数的选择需要考虑泄露现象，当需要分辨的频率为?f时，若使用矩形窗，其宽度必须大于fsam?f。当要求分辨较小的频率而计算量较小可以选择其他窗函数来改善。 在FIR滤波器的设计中，使用矩形窗会由于其锐截止而产生Gibbs现象，会有9%的过冲，为了改善效果，可以在满足过渡带宽度和增益要求的前提下适当选择其他种类的窗函数。 【仿真程序】 %·????÷??′°oˉêyμ??μóòì?D? clear all; L=512;%L±?Dè′óóú??D??í?yμ?3¤?è￡?2μ??yêy′??Y N=10;%ò????÷??′°μ?3¤?è %??D?′° figure(1) w1=zeros(1,50); w2=ones(1,N); w3=zeros(1,50); wh1=[w1 w2 w3]; WH1=fftshift(fft(wh1,L)); w=(0:L-1)-L/2; plot(w,abs(WH1)); title(??D?′°μ?·ù?è?×); %oo?é′° figure(2) wh2=hann(N); WH2=fftshift(fft(wh2,L)); w=(0:L-1)-L/2; plot(w,abs(WH2)); title(oo?é′°μ?·ù?è?×); %1t?÷′° figure(3) wh3=hamming(N); WH3=fftshift(fft(wh3,L)); w=(0:L-1)-L/2; plot(w,abs(WH3)); title(1t?÷′°μ?·ù?è?×); %2?à3???ü′° figure(4) wh4=blackman(N); WH4=fftshift(fft(wh4,L)); w=(0:L-1)-L/2; plot(w,abs(WH4)); title(2?à3???ü′°μ?·ù?è?×); %?-?ó′° figure(5) beta=6; wh5=kaiser(N,beta); WH5=fftshift(fft(wh5,L)); w=(0:L-1)-L/2; plot(w,abs(WH5)); title(?-?ó′°μ?·ù?è?×); 【研讨题目】中等题 Dhexian.wav是对频率为293.66, 369.99, 440Hz的D大调和弦以8000Hz抽样所得的数字音乐信号，试设计数字滤波器从和弦中分离出369.99Hz的音符。要求： （1）设计IIR数字带通滤波器，通过实验研究不同、过渡带、对滤波器设计的影响，确定本题最合适的滤波器指标。 （2）