实验二离散卷积.docVIP

课程名称： 数字信号处理 实验项目： 离散信号的卷积和 专业班级：10电子本 姓名： 学号： 一、实验目的 1、掌握两个离散信号卷积和的计算方法和编程技术。 2、进一步熟悉用MATLAB描绘二维图像的方法。 二、实验原理与计算方法 两个离散序列x(n)与y(n)的卷积和f(n)定义为 　　 由于通常信号处理中所碰到的都是有始信号或有限时间信号，因此在实际计算卷积和时，求和是在有限范围内进行的。计算过程中上下限的选取和所得结果的分布区间取决于参与卷积的两个序列，下面将分别进行讨论： 1、两个从n = 0开始的序列和的卷积和 　　 (1) 上式右边因子u(n)表示卷积和的结果也是一个从n = 0开始的序列。 2、从n = n1开始的序列和从n = n2开始的序列的卷积和，其中n1和n2为任意整数。 (2) 上式右边因子u(n-n1-n2)表示卷积和是一个从n = n1+n2开始的序列。 3、从n = n1开始的长度为N1的加窗序列和从n = n2开始的长度为N2的加窗序列的卷积和，其中 　　　 则　　　 (3) 所得卷积和也是一个加窗序列，从n = n1+ n2开始，长度为N1+ N2-1。 MATLAB提供了一个内部函数conv(x,h)用来计算两个有限长度序列的卷积，该函数得到的卷积结果默认从n=0开始，因此当参与卷积的两个序列的起始位置不是n=0时，则由该函数得到的计算结果将出现错误，此时需要重新定义结果的位置向量。 由以上卷积运算的原理可知，两有限长序列卷积后仍为有限长序列，长度为两序列长度之和减1，结果的起始位置为两序列起始位置之和，截止位置为两序列截止位置之和。据此，可以得到卷积结果的位置向量。 三、实验内容 (1) 计算两个从n = 0开始的序列和的卷积和，其中A=40，a = 0.5。取50个样值点，作出序列、及卷积和f(n)的图像。 k1=0:49; f1=ones(1,50); k2=k1; f2=40*exp(-0.5*k2).*f1 %f1，f2两个信号的向量表示 [f,k]=dconv(f1,f2,k1,k2) %子函数调用 (2) 计算两个有限长序列和的卷积和，其中 　　 作出序列x(n)、y(n)及卷积和f(n)的图像。 k1=-10:9; f1=ones(1,20); k2=5:24; f2=(k2/8).*f1 %f1，f2两个信号的向量表示 [f,k]=dconv(f1,f2,k1,k2) %子函数调用 （3）已知两个序列信号如图所示，试用matlab求两个信号的卷积 ，画出波形 k1=0:3; f1=[2 4 -1 3]; k2=-2:2; f2=ones(1,5);%f1，f2两个信号的向量表示 [f,k]=dconv(f1,f2,k1,k2)%子函数调用 （4）已知两个连续时间信号如下所示，试用matlab求两个信号的卷积 ，画出波形 t1=0~10 t2=0~20 k1=0:3; f1=[2 4 -1 3]; k2=-2:2; f2=ones(1,5);%f1，f2两个信号的向量表示 [f,k]=dconv(f1,f2,k1,k2)%子函数调用 （5） 已知两个连续时间信号如下所示，试用matlab求两个信号的卷积 ，画出波形 f1(t)= e-2t t=0~10 f2(t)=1 t=0~5 p=0.01; k1=0:p:10; f1=exp(-2*k1); k2=0:p:5; f2=ones(1,length(k2)) [f,k]=sconv(f1,f2,k1,k2,p) 四、实验结论 1、掌握了两个离散信号卷积和的计算方法和编程技术。 2、进一步熟悉了用MATLAB描绘二维图像的方法。