数字信号处理上机指导.docVIP

《数字信号处理 上机实验指导书 电子信息工程专业 北京信息 2011年1月目 录实验一 离散傅里叶变换的性质及应用 3 实验二 IIR数字滤波器设计 5 实验一 离散傅里叶变换的性质及应用 一、实验目的 1．了解DFT的性质及其应用 2．熟悉MATLAB编程特点 二、实验仪器及材料 计算机，MATLAB软件 三、实验内容及要求 1． DFT性质。 （1）已知序列DFT的幅度、相位、实部及虚部，并用图形表示相应的 。 （）和 ，验证 *（）序列．利用DFT实现两序列的卷积运算，并研究DFT点数与混叠的关系。 （1）已知两序列 （2）用直接法（即用线性卷积的定义计算，见下式）计算线性卷积y(n)=x(n)*h(n)的结果，并以图形方式表示结果； （3）利用DFT计算圆周卷积；分别令圆周卷积的点数为L=，，以图形方式表示结果。 （4）对比直接法和圆周卷积法所得的结果。．利用DFT。 请问信号的最高频率是多少？ （2）假定抽样频率，编程实现记录数据点数为N=16,32,64 三种情况下的谱分析。分别画出其幅度谱，讨论三种情况幅度谱的区别？幅度谱与频率分辨率 四、思考题 说明为什么L=8，时采用FFT法对两序列计算圆周卷积的结果不同？与线性卷积y(n)=x(n)*h(n)的结果说明为什么、实验报告要求 1．列出本实验编写的所有文件及各项实验结果曲线，加注必要的说明； 2．对实验结果作理论计算，解释实验结果； 3．总结实验体会及实验中存在的问题、本实验所用的MATLAB函数 (x) : 对复数取模（对实数取绝对值）。(x) : 对复数取相角。 (x) : 对复数取实部。 (x) : 对复数取虚部。 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．设计巴特沃斯低通数字滤波器人体心电图信号在测量过程中会受到工业高频干扰，所以必须经过低通滤波处理，才能作为判断心脏功能的有用信息。以下为一个实际心电图信号采样序列x(n)，其中存在高频干扰，采样周期T=1秒。在实验中，以x(n)作为输入序列，滤除其中干扰成分。 x(n)=[-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0] 对序列x(n)用FFT做频谱分析，生成x(n)的频谱图。 用设计一个巴特沃斯低通IIR数字滤波器H(z)。 设计指标参数为：在通带内频率低于0.2π时，最大衰减小于dB； 在阻带内 [0.3π, π]频率区间上，最小衰减大于dB。 数字滤波器H(z)滤波器的幅频响应曲线。 用所设计的滤波器对实际心电图信号采样序列x(n)进行滤波处理，编写程序滤波后的序列y(n)并分别出滤波前后的心电图信号波形图和频谱图。cheby1.m，了解调用格式。 编程设计教材习题6-2，求模拟滤波器H(s)的表达式，画出该模拟滤波器频响。bilinear实现教材习题6-5，求数字滤波器H(z)。、思考题 对比滤波前后心电图信号x(n), y(n),说明滤波器的作用。、实验报告要求1．列出本实验编写的所有文件及各项实验结果曲线，加注必要的说明； 2．总结实验体会及实验中存在的问题、本实验所用的MATLAB函数 length(x) ：求序列x长度。 y=fft(x,L)将序列x(n)做L点快速傅立叶变换，结果赋给序列y(n)。 [n,Wn] = buttord(Wp,Ws,Rp,Rs,s) ：计算模拟Butterworth滤波器的最小阶次n和截止频率为Wn。Wp,Ws为模拟角频率，单位为rad/sRp,Rs的单位为 [b,a] = butter(n,Wn,s)：设计模拟截止频率为Wn（rad/s）的n阶 Butterworth低通滤波器，返回值为模拟滤波器的系数。 [n,Wn] = buttord(Wp,Ws,Rp,Rs)