数字信号处理课程设计IIR数字滤波器以实现Word下载.docx

1、%计算阶数N和3Db处的频率 N1 %在命令框中显示NB1,A1=butter（N1,wc1）;%求巴特沃斯型数字滤波器w=linspace（0,2*pi,1000）;%在0-2*pi内将w分成1000个点h1=freqz（B1,A1,w）;%调用freqz函数，求解幅频响应magh1=abs（h1）;%求出幅频特性pha1=angle（h1）;%求出相频特性subplot（321）;plot（w/pi,magh1,r）;%画图title（低通滤波器幅频特性曲线xlabel（模拟频率subplot（322）;plot（w/pi,pha1,低通滤波器相频特性曲线n=0:72;%选取序列长度w=p

2、i*n/73;%分出与n相对应的wxn=1+cos（pi*n/4）+cos（2*pi*n/3）;%原序列subplot（323）;stem（n,xn,.%画出原序列xn的时域曲线nylabel（xnX=abs（xn*exp（-1*j*n*w）;%原序列的序列傅里叶变换subplot（3,2,4）;stem（w,X,%画出序列傅里叶变换的曲线xn的频域曲线wXy1=filter（B1,A1,xn）;%利用filter函数恢复出滤波后的函数subplot（325）;stem（n,y1,.m低通滤波后的时域曲线y1Y1=abs（y1*exp（-1*j*n%恢复后的函数的序列傅里叶变换subplot（

3、3,2,6）;stem（w,Y1,低通滤波后的频域曲线Y1%（2）设计数字高通滤波器figure（2）;wp2=0.5;ws2=0.4;ap2=1;as2=30;N2,wc2=buttord（wp2,ws2,ap2,as2）; N2 B2,A2=butter（N2,wc2,high w=linspace（0,2*pi,1000）;h2=freqz（B2,A2,w）;magh2=abs（h2）;pha2=angle（h2）;plot（w/pi,magh2,高通滤波器幅频特性曲线plot（w/pi,pha2,高通滤波器相频特性曲线y2=filter（B2,A2,xn）;stem（n,y2,高通滤波

4、后的时域曲线y2Y2=abs（y2*exp（-1*j*nstem（w,Y2,高通滤波后的频域曲线Y2%（3）设计数字带通滤波器figure（3）;wp3=0.1,0.4;ws3=0.05,0.5;ap3=1;as3=30;N3,wc3=buttord（wp3,ws3,ap3,as3）; N3 B3,A3=butter（N3,wc3）;h3=freqz（B3,A3,w）;magh3=abs（h3）;pha3=angle（h3）;plot（w/pi,magh3,带通滤波器幅频特性曲线plot（w/pi,pha3,带通滤波器相频特性曲线y3=filter（B3,A3,xn）;stem（n,y3,带通

5、滤波后的时域曲线y3Y3=abs（y3*exp（-1*j*nstem（w,Y3,带通滤波后的频域曲线Y3%（4）设计数字带阻滤波器figure（4）;wp4=0.05,0.5;ws4=0.1,0.4;ap4=1;as4=30;N4,wc4=buttord（wp4,ws4,ap4,as4）; N4 B4,A4=butter（N4,wc4,stoph4=freqz（B4,A4,w）;magh4=abs（h4）;pha4=angle（h4）;plot（w/pi,magh4,带阻滤波器幅频特性曲线plot（w/pi,pha4,带阻滤波器相频特性曲线y4=filter（B4,A4,xn）;stem（n,

6、y4,带阻滤波后的时域曲线y4Y4=abs（y4*exp（-1*j*nstem（w,Y4,带阻滤波后的频域曲线Y45.测试数据及必要的理论分析与比较由figure1可以看出低通滤波器的幅频特性曲线和相频特性曲线，原信号的时域曲线和频域曲线以及滤波后的信号的时域曲线和频域曲线，通过对比，可以很清楚的看出该低通滤波器实现了题目所要求的滤波功能；由figure2可以看出高通滤波器的幅频特性曲线和相频特性曲线，原信号的时域曲线和频域曲线以及滤波后的信号的时域曲线和频域曲线，通过对比，可以很清楚的看出该高通滤波器实现了题目所要求的滤波功能；由figure3可以看出带通滤波器的幅频特性曲线和相频特性曲线，

7、原信号的时域曲线和频域曲线以及滤波后的信号的时域曲线和频域曲线，通过对比，可以很清楚的看出该带通滤波器实现了题目所要求的滤波功能；由figure4可以看出带阻滤波器的幅频特性曲线和相频特性曲线，原信号的时域曲线和频域曲线以及滤波后的信号的时域曲线和频域曲线，通过对比，可以很清楚的看出该带阻滤波器实现了题目所要求的滤波功能；（六）第六题（1）用Hanning窗设计一线性相位带通数字滤波器，要求：N=15, 。观察它的实际3dB和20dB带宽。N=45，重复这一设计，观察幅频和相位特性的变化，注意长度N变化的影响；（2）分别改用矩形窗和Blackman窗，设计（1）中的带通滤波器，观察并记录窗函数

8、对滤波器幅频特性的影响，比较三种窗的特点；总结窗的不同长度和不同窗对滤波器的影响本题主要考察利用窗函数法设计线性相位的数字滤波器，与之前的设计方法并无太大区别，先是确定滤波器的技术指标（归一化），然后通过调用窗函数的语句进行设计，只不过这里需要根据窗函数的不同选择相应的过渡带宽及选择信号的长度和滤波器的阶数；同样，带通滤波器的截止频率需要用区间形式写出；对于不同的窗函数，过渡带宽会不同。关键部分应该有确定截止频率（注意要归一化），还有选择不同的窗函数对应的过渡带宽；另外就是阶数和信号长度的关系：信号长度=阶数+1，这是从资料书上看来的，所以每个滤波器的阶数都是按照这个等式来确定的。%第六题%（

9、1）用汉宁窗设计线性相位带通滤波器%当长度N1为15时N1=15;w1=0.3*pi;w2=0.5*pi;wc=0.3,0.5;%归一化后的通带截频B1=8*pi/N1;%过渡带宽M1=N1-1;%单位脉冲响应的h（n）的长度N=M+1,此处的M1是阶数；n1=0:M1;hn1=fir1（M1,wc,hanning（N1）;%调用firl函数求解单位脉冲响应h1,w1=freqz（hn1,1,512）;%求其频率响应%求其幅频特性%求其相频特性subplot（221）;plot（w1/pi,10*log（magh1）;%画出其幅频特性grid;N=15时,汉宁窗的幅频特性maghsubplot

10、（222）;plot（w1/pi,pha1）;N=15时,汉宁窗的相频特性pha%当长度N2为45时N2=45;B2=8*pi/N2;M2=N2-1;n2=0:M2;hn2=fir1（M2,wc,hanning（N2）;h2,w2=freqz（hn2,1,512）;subplot（223）;plot（w2/pi,10*log（magh2）,gN=45时,汉宁窗的幅频特性subplot（224）;plot（w2/pi,pha2,N=45时,汉宁窗的相频特性%（2）用矩形窗设计线性相位带通滤波器%当长度N3为15时N3=15;B3=4*pi/N3;M3=N3-1;n3=0:M3;hn3=fir1（

11、M3,wc,boxcar（N3）;h3,w3=freqz（hn3,1,512）;plot（w3/pi,10*log（magh3）;N=15时,矩形窗的幅频特性plot（w3/pi,pha3）;N=15时,矩形窗的相频特性N4=45;B4=4*pi/N4;M4=N4-1;n4=0:M4;hn4=fir1（M4,wc,boxcar（N4）;h4,w4=freqz（hn4,1,512）;plot（w4/pi,10*log（magh4）,N=45时,矩形窗的幅频特性plot（w4/pi,pha4,N=45时,矩形窗的相频特性%（3）用Blackman窗设计线性相位带通滤波器N5=15;B5=12*pi

12、/N5;M5=N5-1;n5=0:M5;hn5=fir1（M5,wc,blackman（N5）;h5,w5=freqz（hn5,1,512）;magh5=abs（h5）;pha5=angle（h5）;plot（w5/pi,10*log（magh5）;N=15时,Blackman窗的幅频特性plot（w5/pi,pha5）;N=15时,Blackman窗的相频特性N6=45;B6=12*pi/N6;M6=N6-1;n6=0:M6;hn6=fir1（M6,wc,blackman（N6）;h6,w6=freqz（hn6,1,512）;magh6=abs（h6）;pha6=angle（h6）;plot

13、（w6/pi,10*log（magh6）,N=45时,Blackman窗的幅频特性plot（w6/pi,pha6,N=45时,Blackman窗的相频特性 由figure1可以看出，N=15时的汉宁窗设计的滤波器的通带比N=45时的通带宽，旁瓣的衰减也要小一些，所以，阶数越大，所设计出来的滤波器效果越好，通带稳定，衰减快。同理，根据fugure2和figure3的矩形窗和Blackman窗也可以得出相同的结论。 通过对三个图不同窗函数设计的滤波器进行比较，可以看出：矩形窗设计的滤波器通带最稳定，衰减最快；汉宁窗次之；Blackman窗滤波效果相对最差。所以，矩形窗更能较好的实现滤波效果。（七）

14、选做题7音乐信号处理： 1）获取一段音乐或语音信号，设计单回声滤波器，实现信号的单回声产生。给出加入单回声前后的信号频谱。 2）设计多重回声滤波器，实现多重回声效果。给出加入多重回声后的信号频谱。 3）设计均衡器，使得得不同频率的混合音频信号，通过一个均衡器后，增强或削减某些频率区域。设计思想为：先从网上下载一段格式为wav的音乐信号，选择一个较长的长度，设计一个滤波器，然后用该滤波器对音乐信号进行滤波，实质上就是将信号拉长后，与原信号就行一次叠加，便可设计出单回声滤波器；给滤波器加入一个阶次，让滤波器进行多次滤波然后再与原信号叠加，便可得到多重回声滤波器。设计单回声滤波器的关键在于该滤波器的

15、系统函数的分子部分，要通过补零的方式将滤波器进行拉长，最好还要加入一个固定常数项；设计多重回声滤波器时，系统函数分子部分补的零根据阶数的不同会相应增加，最后的常数项也要变为相应的阶次幂，而且分母部分也应该进行补零。%第七题%（1）单回声滤波器%addpath（D:Progame Files/MATLAB/R2011a/binx,fs,bits=wavread（Prettyboy.wav,216）;%读出音频文件%x,fs,bits=wavread（yue.wavwavplay（x,fs）;%播放音频文件fsbitspause（1）;a1=0.8;R1=5000;B1=1,zeros（1,R1-

16、1）,0.8;%系统函数分子A1=1;y=filter（B1,A1,x）;%滤波器函数wavplay（y,fs）;x1=fft（x）;plot（abs（x1）;单回声前信号幅频响应plot（angle（x1）;单回声前信号相频响应y1=fft（y）;%快速傅里叶变换plot（abs（y1）,m单回声后信号幅频响应subplot（324）;plot（angle（y1）,单回声后信号相频响应h1,w1=freqz（B1,A1）;%求滤波器的幅频响应函数plot（abs（h1）,单回声滤波器幅频响应y2=ifft（h1）;subplot（326）;stem（abs（y2）,.g单回声滤波器冲激响应%

17、多重回声滤波器pause（2）;pause（0）;N=5;a2=0.8;R2=8000;B2=1,zeros（1,N*R2-1）,-0.8N;A2=1,zeros（1,R2-1）,-0.8;z=filter（B2,A2,x）;wavplay（z,fs）;z1=fft（z）;多重回声前信号幅频响应多重回声前信号相频响应plot（abs（z1）,多重回声后幅频响应plot（angle（z1）,多重回声后相频响应h2,w2=freqz（B2,A2）;plot（abs（h2）,多重回声滤波器幅频响应z2=ifft（h2）;stem（abs（z2）,多重回声滤波器冲激响应由figure1可以看出，声音信号通过单回声滤波器后整个频率域上的幅值都发生了一些变化，变得比以前稍大一些，这是因为滤波后的信号是原信号与回声信号叠加起来的；同理，由figure2可以看出，声音信号通过多重回声滤波器后整个频率域上的幅值也都发生了一些变化，变得比以前更大。 通过对比figure1和figure2可以看出，单回声滤波器的幅频响应幅度比相应的多重回声滤波器的幅度小，因为多重回声滤波器是经过N次的单回声滤波器滤波后叠加的效果；而且单回声滤波器的冲激响应只有一个值，多