《數字信号处理上机实验指导书》09级.docVIP

《数字信号处理》 上机实验指导书 光电信息与通信工程学院 通信工程专业 北京信息科技大学 2011年10月 目 录 实验一 离散傅里叶变换的性质及应用 1 实验二 IIR数字滤波器设计 4 实验三 FIR数字滤波器设计 6 实验一 离散傅里叶变换的性质及应用 一、实验目的 1．了解DFT的性质及其应用 2．熟悉MATLAB编程特点 二、实验仪器及材料 计算机，MATLAB软件 三、实验内容及要求 1．用三种不同的DFT程序计算的256点离散傅里叶变换，并比较三种程序计算机运行时间。 （1）编制用for loop语句的M函数文件dft1.m，用循环变量逐点计算； （2）编写用MATLAB矩阵运算的M函数文件dft2.m，完成下列矩阵运算： （3）调用fft库函数，直接计算； （4）分别调用上述三种不同方式编写的DFT程序计算序列的离散傅里叶变换，并画出相应的幅频和相频特性，再比较各个程序的计算机运行时间。 2．利用DFT实现两序列的卷积运算，并研究DFT点数与混叠的关系。 （1）已知两序列： ，用MATLAB生成随机输入信号x(n)，n的取值为0～2； （2）用直接法（即用线性卷积的定义计算，见下式）计算线性卷积y(n)=x(n)*h(n)的结果，并以图形方式表示结果； 其中：序列和序列 （3）用MATLAB编制利用DFT计算线性卷积y(n)=x(n)*h(n)的程序；分别令圆周卷积的点数为L=5，6，8，10，以图形方式表示结果。 （4）对比直接法和圆周卷积法所得的结果。 四、思考题 直接计算N=256点的DFT与用FFT计算256点，理论上速度差别应有多少？ 说明为什么在实验内容2中L=6，8，10时采用FFT法对两序列计算圆周卷积的结果与线性卷积y(n)=x(n)*h(n)的结果相同，而与L=5时计算结果不同？ 五、实验报告要求 1．列出本实验编写的所有文件及各项实验结果曲线，加注必要的说明； 2．对实验结果作理论计算，解释实验结果； 3．回答思考题； 4．总结实验体会及实验中存在的问题。 六、参考程序 1．dft1.m：用for循环计算DFT。 function[Am,pha] = dft1(x) N = length(x); w = exp(-j*2*pi/N); for k=1:N sum = 0; for n = 1:N sum = sum+x(n)*w^((k-1)*(n-1)); end Am(k) = abs(sum); pha(k) = angle(sum); end 2．dft2.m：用MATLAB矩阵计算DFT。 function[Am,pha] = dft2(x) N = length(x); n = [0:N-1];k = [0:N-1]; w = exp(-j*2*pi/N);nk = n’*k; wnk = w.^(nk); Xk = x*wnk; Am = abs(Xk); pha = angle(Xk); 3．dft3.m：调用FFT库函数计算DFT。 function[Am,pha] = dft3(x) Xk = fft(x); Am = abs(Xk); pha = angle(Xk); 4．lab_dft.m：调用M函数文件的M程序文件示例。 x = [ones(1,8),zeros(1,504)]; t = cputime; [Am1,pha1] = dft1(x); t1 = cputime-t; 5．本实验所用的部分MATLAB函数： x=rand(r,c) : 用MATLAB随机生成输入信号x（n），范围为0～1,r行,c列。 y=conv(x1,x2): 将序列x1(n)与序列x2(n)做线性卷积，结果赋给序列y(n)。 y=fft(x,L) : 将序列x(n)做L点快速傅立叶变换，结果赋给序列y(n)。 y=ifft(x,L) : 将序列x(n)做L点快速傅立叶反变换，结果赋给序列y(n)。 stem(x,y) : 以x为横坐标，y为纵坐标画离散曲线。 实验二 IIR数字滤波器设计 一、实验目的 1、熟悉用双线性变换法设计IIR数字滤波器的原理与方法。 2、掌握数字滤波器的计算机软件实现方法。 3、通过观察对实际心电图信号的滤波作用，学习数字滤波器在实际中的应用。 二、实验仪器及材料 计算机，MATLAB软件 三、实验内容及要求 1．设计巴特沃斯低通数字滤波器对人体心电信号进行滤波 (1) 人体心电图信号在测量过程中会受到工业高频干扰，所以必须经过低通滤波处理，才能作为判断心脏功能的有用信息。以下为一个实际心电图信号采样序列x(n)，其