《数 字 信 号 处 理》 实 验 指 导 书.ppt

数字信号处理实验,电工电子实验中心,前言,本课程实验包含四个必修实验：实验一离散系统时域分析为验证性实验，使学生掌握计算离散卷积方法和MATLAB语言实现；实验二FIR数字滤波器的设计为综合性，使学生掌握用窗函数法和频率采样法设计FIR数字滤波器的原理和方法，了解不同窗函数对滤波器性能的影响；实验三IIR数字滤波器的设计为综合性，可使学生掌握脉冲响应不变法和双线性变换法设计IIR数字滤波器的原理和方法，了解双线性变换法和脉冲响应不变法的特点和区别。实验四FFT频谱分析及应用为验证性，可进一步加深对FFT的理解，熟悉应用FFT对典型信号进行频谱分析的方法。,目录,实验一离散系统时域分析实验二FIR数字滤波器的设计实验三IIR数字滤波器的设计实验四FFT频谱分析及应用附录：MATLAB基本操作及常用命令,实验一：离散系统时域分析,实验学时：2学时实验类型：验证实验要求：必修一、实验目的1、学习MATLAB语言的编程和调试技巧；2、掌握笔算离散卷积方法及其MATLAB语言实现。,二、实验内容,1、设某LTI的单位脉冲响应（1）判断此系统是否可实现；（2）当输入为矩形脉冲时，求此LTI的输出；（3）用MATLAB实现，并画出图形。2、计算卷积,三、实验原理与方法和手段,一个离散时间系统，输入信号为x(n)，输出信号为y(n)，运算关系用T表示，则输入与输出的关系可表示为y(n)Tx(n)。（1）线性时不变系统的输入输出关系可通过单位脉冲响应h(n)表示：y(n)=x(n)*h(n)=式中*表示卷积运算。,三、实验原理与方法和手段,（2）线性时不变系统的实现可物理实现的线性时不变系统是稳定的、因果的。这种系统的单位脉冲响应是因果的（单边）且绝对可和的，即：在MATLAB语言中采用conv实现卷积运算即：Y=conv(x,h),它默认从n=0开始。,四、实验组织运行要求,1、学生在进行实验前必须进行充分的预习，熟悉实验内容；2、学生根据实验要求，读懂并理解相应的程序；3、学生严格遵守实验室的各项规章制度，注意人身和设备安全，配合和服从实验室人员管理；4、教师在学生实验过程中予以必要的辅导，独立完成实验；,五、实验条件,1、具有WINDOWS98/2000/NT/XP操作系统的计算机一台；2.、MATLAB编程软件。,六、实验步骤,在“开始-程序”菜单中，找到MATLAB程序，运行启动；进入MATLAB后，首先熟悉界面；在CommandWindow中输入参考程序，并执行；记录运行结果图形，并与笔算结果对照。再根据实验内容，逐一修改程序，执行，记录结果并分析。,七、思考题,结合信号与系统课程所学，思考离散时间系统的线性卷积公式与连续时间系统的卷积公式的异同？,八、实验报告要求,1、报告中要给出实验的MATLAB程序，并对每个语句给出注释，说明语句作用；2、简述实验目的和原理；3、给出用笔算时卷积和conv计算线性卷积对照图；4、给出收获和体会。,九、参考程序,程序1x=ones(1,10);x1=ones(1,10),zeros(1,40);N1=length(x);n1=0:N1-1;N2=50;n2=0:N2-1;h=0.8.n2;y=conv(x,h);N=N1+N2-1;n=0:N-1;subplot(3,1,1);stem(n2,x1);subplot(312);stem(n2,h);subplot(313);stem(n,y);,九、参考程序,程序2如果、的起点不为0，则采用conv_m计算卷积；编写conv_m函数：functiony,ny=conv_m(x,nx,h,nh)%改进卷积程序nyb=nx(1)+nh(1);nye=nx(length(x)+nh(length(h);ny=nyb,nye;y=conv(x,h);在命令窗口输入：x=3,11,7,0,-1,4,2;nx=-3:3;h=2,3,0,-5,2,1;nh=-1:4;y,ny=conv_m(x,nx,y,ny)可得到结果：y(n)=6,31,47,6,-51,-5,41,18,-22,-3,8,2，-4”为运算提示符，表示MATLAB正处在准备状态。当在提示符后输入一段运算式并按Enter键后，MATLAB将给出计算结果，然后再次进入准备状态。（3）工作间窗口在主窗口的“View”菜单下选择“Workspace”命令，可以打开或关闭MATLAB的工作间窗口，工作间窗口是MATLAB的重要组成部分。,MATLAB的程序类型有三种1）脚本M文件在命令窗口中输入并执行，它所用的变量都要在工作空间中获取，不需要输入输出参数的调用，退出MATLAB后就释放了。2）程序M文件以.m格式进行存取，包含一连串的MATLAB指令和必要的注解。需要在工作空间中创建并获取变量，也就是说处理的数据为命令窗口中的数据，没有输入参数，也不会返回参数。程序运行时只需在命令窗口中键入其名称即可（不需要后缀名）。3）函数M文件（常称为M函数）与在命令窗口中输入命令一样，函数接受输入参数，然后执行并输出结果。用help命令可以显示它的注释说明。,M函数,具有标准的基本结构：（1）函数定义行（关键字function）functionout1,out2,.=filename(in1,in2,.)输入和输出（返回）的参数个数分别由nargin和nargout两个MATLAB保的变量来给出。（2）第一行帮助行以（%）开头，作为lookfor指令搜索的行（3）函数体说明及有关注解以（%）开头，用以说明函数的作用及有关内容（4）函数体语句函数体内使用的除返回和输入变量这些在function语句中直接引用的变量以外的所有变量都是局部变量，即在该函数返回之后，这些变量会自动在MATLAB的工作空间中清除掉。如果希望这些中间变量成为在整个程序中都起作用的变量，则可以将它们设置为全局变量。,例如，实现加减法运算的M函数调用方法如下：首先编写一个M函数：（打开一个程序编辑器“Untitled”）Functionoutl,out2=addec(inl，in2)outl=inl+in2；out2=inl-in2；保存的文件名为addec（默认后缀名为.m）然后，在MATLAB的命令窗口调用这个函数文件：a,b=addec(10,20)disp(a,b)%在命令窗口输出a、b的值显示结果为：3010,3、基本平面图形绘制命令plot,1）plot(X,Y)%当X,Y均为实数向量，且为同维向量，则以x的元素为横坐标，y的元素为纵坐标，绘制曲线。若X，Y均为同维同型实数矩阵，X=X(i)，Y=Y(i)，其中X(i),Y(i)为列向量，则plot(X,Y)依次画出plot(X(i),Y(i)，矩阵有几列就有几条线；2）plot(X1,Y1,X2,Y2,)，