数字信号处理实验讲义-长师.docVIP

前言 数字信号处理的理论和技术一出现就受到人们的极大关注，发展非常迅速。20世纪60年代以来，随着信息科学和计算技术的迅速发展，数字信号处理的理论与应用飞跃发展，它的重要性日益在通信、信息技术、图像处理、遥感、声呐、雷达、生物医学、地震、语音处理等领域的应用中表现出来，并发展成为一门极其重要的学科。数字信号处理由于信号处理涉及大量的运算，可以说离开了计算机及相应的软件，就不可能解决任何稍微复杂的实际应用问题。Matlab是1984年美国MathWorks公司的产品，MATLAB语言具备高效、可视化及推理能力强等特点，它的推出得到了各个领域专家学者的广泛关注，其强大的扩展功能为各个领域的应用提供了基础，是目前工程界流行最广的科学计算语言。早在20世纪90年代中期，MATLAB就己成为国际公认的信号处理的标准软件和开发平台。从1996年后，美国新出版的信号处理教材就没有一本是不用MATLAB的。用科学计算语言MATLAB实现数字信号处理的方法和实践，通过实验用所学理论来分析解释程序的运行结果，进一步验证、理解和巩固学到的理论知识，从而达到掌握数字信号处理的基本原理和方法的目的。实验一 常见离散信号的MATLAB产生和图形显示一、实验目的１学会用MATLAB在时域中产生一些基本的离散时间信号。２了解信号的各种运算。二、实验原理（一）序列的产生由于MATLAB数值计算的特点，用它来分析离散时间信号与系统是很方便的。离散信号是数字信号处理的最基础的内容，由于内存有限，MATLAB无法表示无限序列。在MATLAB中，可以用一个列向量来表示一个有限长度的序列，但是这种表示方法没有包含采样位置的信息，要完全表示x(n)，要用x和n两个向量，例如x(n)={2,1,0,2,3,-1,2,3}，在MATLAB中表示为n=[-4,-3,-2,-1,0,1,2,3];x=[2,1,0,2,3,-1,2,3];当序列从n=0开始，则不需要采样位置信息，这时可以只用x来表示1．正弦序列 离散正弦序列的MATLAB表示与连续信号类似，只不过是用stem函数而不是用plot函数来画出序列的波形。下面就是正弦序列的MATLAB源程序。程序运行结果如图1.19所示。 %正弦序列实现程序 k=0:39; fk=sin(pi/6*k); stem(k,fk) 图1.19，可用MATLAB中的数组幂运算（即点幂运算）c*来实现。下面为用MATLAB编写绘制离散时间实指数序列波形的函数。 functiondszsu(c,a,k1,k2) %c：指数序列的幅度 %a：指数序列的底数 %k1：绘制序列的起始序号 %k2：绘制序列的终止序号 k=k1:k2; x=c*(a.^k); stem(k,x,filled) holdon plot([k1,k2],[0,0]) holdoff 利用上述函数，实现实指数波形MATLAB程序如下（其中值分别为）。 %离散时间实指数序列实现程序 subplot221; dszsu(1,5/4,0,20); xlabel(k); title(f1[k]); subplot222 dszsu(1,3/4,0,20); xlabel(k); title(f2[k]); subplot223; dszsu(1,-5/4,0,20); xlabel(k); title(f3[k]); subplot224; dszsu(1,-3/4,0,20); xlabel(k); title(f4[k]); 图1.20，当的绝对值大于1时，序列为随时间发散的序列，当的绝对值小于1时，序列为随时间收敛的序列。同时可见，当的值小于零时，其波形在增长或衰减的同时，还交替地改变序列值的符号。 对于离散时间虚指数序列，可用通过调用下列绘制虚指数序列时域波形的MATLAB函数。 function[]=dxzsu(n1,n2,w) %n1:绘制波形的虚指数序列的起始时间序号 %n2：绘制波形的虚指数序列的终止时间序号 %w：虚指数序列的角频率 k=n1:n2; f=exp(i*w*k); Xr=real(f) Xi=imag(f) Xa=abs(f) Xn=angle(f) subplot(2,2,1),stem(k,Xr,filled),title(实部); subplot(2,2,3),stem(k,Xi,filled),title(虚部); subplot(2,2,2),stem(k,Xa,filled),title(模); subplot(2,2,4),stem(k,Xn,filled),title(相角); 利用上述函数，实现