数字信号处理实验报告37020new

《数字信号处理实验报告37020new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、《数字信号处理》实验报告专业电子信息工程学号091308304姓名丁宝宝实验一利用FFT实现快速卷积一、实验目的1．加深理解FFT在实现数字滤波（或快速卷积）中的重要作用，更好的利用FFT进行数字信号处理。2．掌握循环卷积和线性卷积两者之间的关系。二、实验原理用FFT来快速计算有限长度序列的线性卷积。这种方法就是先将输入信号x(n)通过FFT变换为它的频谱采样值，然后再和FIR滤波器的频响采样值H(k)相乘，H(k)可事先存放在存储器中，最后再将乘积H(k)X(k)通过快速傅里叶变换（简称IFF

2、T）还原为时域序列，即得到输出y(n)。现以FFT求有限长序列的卷积及求有限长度序列与较长序列间的卷积为例来讨论FFT的快速卷积方法。序列x(n)和h(n)的长差不多。设x(n)的长为N1，h(n)的长为N2，要求用FFT完成这一卷积的具体步骤如下：①为使两有限长序列的线性卷积可用其循环卷积代替而不发生混叠，必须选择循环卷积长度，若采用基2-FFT完成卷积运算，要求（为整数）。②用补零方法使x(n)和h(n)变成列长为N的序列。③用FFT计算x(n)和h(n)的N点离散傅里叶变换。④完成X(k)

3、和H(k)乘积，⑤用FFT计算的离散傅里叶反变换得三、主要实验仪器及材料微型计算机、Matlab6.5教学版。四、实验内容1．数字滤波器的脉冲响应为。输入序列x(n)可选下列几种情况①。②。③。2．实验前，预先编制一个应用FFT实现数字滤波器的通用程序。3．上机独立调试，并打印或记录实验结果。4．将实验结果与预先笔算的结果比较，验证其正确性。五、思考题1．N同时取8、16、32时，线性卷积和循环卷积的结果有何不同，为什么？2．分析直接计算线性卷积和利用FFT计算线性卷积的时间。附程序xn1=[1

4、,1/2,1/4,1/8,1/16,1/32,1/64,1/128];%建立x1(n)序列xn2=[1,1,1,1,1,1,1,1];%建立x2(n)序列N=length(xn1);n=0:N-1;k=0:N-1;Xk1=fft(xn1,N);%由x1(n)的FFT求X1(k)Xk2=fft(xn2,N);%由x2(n)的FFT求Yk=Xk1.*Xk2;%Y(k)=X1(k)X2(k)yn=ifft(Yk,N)%由Y(k)的IFFT求y(n)yn=abs(yn)subplot(2,3,1),st

5、em(n,xn1,'k');title('x1(n)');subplot(2,3,2),stem(n,xn2,'k');title('x2(n)');subplot(2,3,3),stem(n,yn,'k');title('y(n)');subplot(2,3,4),stem(k,abs(Xk1),'k');title('X1(k)');subplot(2,3,5),stem(k,abs(Xk2),'k');title('X2(k)');subplot(2,3,6),stem(k,abs(Yk)

6、,'k');title('Y(k)');实验二数字滤波器的结构一、实验目的（1）加深对数字滤波器分类与结构的了解；（2）明确数字滤波器的基本结构及其相互间的转换方法；（3）掌握用MATLAB进行数字滤波器各种结构相互间转换的子函数及程序编写方法。二、实验原理一个离散LSI系统可用系统函数来表示；也可用差分方程来表示：当至少有一个不为0时，则在有限z平面上存在极点，表示一个IIR数字滤波器；当全都为0时，系统不存在极点，表示一个FIR系统。IIR数字滤波器的基本结构分为直接Ⅰ型、直接Ⅱ型、级联型和

7、并联型。FIR数字滤波器的基本结构分为横截型、级联型、并联型、、线性相位型和频率抽样型。三、实验仪器微型计算机、MATLAB四、实验内容(1)已知一个IIR系统的系统函数为将其从直接型转换为级联型和并联型结构，并画出各种结构的流程图。(2)已知一个FIR系统的系统函数为将其从横截型转换为级联型结构，并画出各种结构的流程图。一、实验程序b=[8,-4,11,-2];%输入系统函数b参数a=[1,-1.25,0.75,-0.125];%输入系统函数a参数[sos,g]=tf2sos(b,a)%由直接

8、型转换为级联型[r,p,k]=residuez(b,a)%由直接型转换为并联型二、试验结果分析实验三用双线性变换法设计IIR数字滤波器一、实验目的（1）加深对双线性变换法设计FIR数字滤波器的基本原理的理解；（2）掌握用双线性变换法设计数字低通、高通、带通滤波器的设计；（3）了解MATLAB有关双线性变换法设计的常用子函数。二、实验原理1.设计思想s平面与z平面之间满足以下映射关系：s平面的虚轴单值地映射于z平面的单位圆上，s平面的左半平面完全映射到z平面的单位圆内。双线性变换不存在混叠问题。双