連续时间信号卷积运算的实现(上海电力学院).docxVIP

时间连续信号卷积运算的MATLAB实现2020022班陈彬彬实验目的1、理解掌握卷积的概念及物理意义。2、理解单位冲击响应的概念及物理意义。实验原理根据前述知识,连续信号的卷积运算定义为 （9-1）卷积计算可以通过信号分段求和来实现，即(9-2)如果只求当时的值，则由上式可得 (9-3)式(9-3)中的实际上就是连续信号经等时间间隔均匀抽样的离散序列的卷积和。当足够小时，就是卷积积分的结果——连续时间信号的较好的数值近似。三、实验内容与方法1、用MATLAB实现连续信号f1(t)和f2(t)卷积的过程如下：将连续信号f1(t)和f2(t)以时间间隔Δ进行取样，得到离散序列 f1(kΔ)、f2(kΔ)；2、构造f1(kΔ)、f2(kΔ)与相对应的时间向量k1和k2；①调用conv()函数计算卷积积分f(t)的近似向量f(nΔ)；②构造f(nΔ)对应的时间向量k。下面是利用MATLAB实现连续信号卷积运算的通用函数sconv()，该程序在计算出卷积积分的数值近似值的同时，还绘出f(t)的波形图。需要注意的是，程序中是如何构造f(t)的对应时间向量k的？另外，程序在绘制f(t)波形图采用的是plot命令而不是stem命令。function [f,k]=sconv(f1,f2,k1,k2,p)%计算连续信号卷积积分f(t)=f1(t)*f2(t)%f:卷积积分f(t)对应的非零值向量%k:f(t)的对应时间向量%f1:f1(t)的非零样值向量%f2:f2(t)的非零样值向量%k1:f1(t)的对应时间向量%k2: 序列f2(t)的对应时间向量%p:取样时间间隔f=conv(f1,f2); %计算序列f1和f2的卷积和ff=f*p;k0=k1(1)+k2(1); %计算序列f非零样值的起点位置k3=length(f1)+length(f2)-2; %计算卷积和f的非零样值的宽度k=k0:p:k3*p; %确定卷积和f非零样值的时间向量subplot(2,2,1)plot(k1,f1) %在子图1绘出f1(t)时域波形图title(f1(t))xlabel(t)ylabel(f1(t))subplot(2,2,2)plot(k2,f2) %在子图2绘出f2(t)时域波形图title(f2(t))xlabel(t)ylabel(f2(t))subplot(2,2,3)plot(k,f); %画出卷积f(t)的波形h=get(gca,position);h(3)=2.5*h(3);set(gca,position,h) %将第三个子图的横坐标扩展为原来的2.5倍Title(f(t)=f1(t)*f2(t))xlabel(t)ylabel(f(t))四、实验过程1、对书中的例题进行仿真：在MATLAB命令行窗口键入＞＞edit指令或者通过“File”菜单中的“New”子菜单下的“M－File”命令或者单击工具栏上的新建按扭。将连续信号卷积运算的通用函数sconv()的程序输入在内，然后将文件保存，扩展名设置为．Ｍ。然后在命令行中键入以下内容：当p=0.5p=0.5;k1=0:p:2;f1=0.5*k1;k2=0:p:2;k2=k1[f,k]=sconv(f1,f2,k1,k2,p) 当p=0.01时p=0.01;k1=0:p:2;f1=0.5*k1;k2=0:p:2;k2=k1f2=f1; [f,k]=sconv(f1,f2,k1,k2,p) 2、试求下列图片的卷积波形列出编程步骤：p=0.01;k1=0:p:1;f1=ones(1,length(k1));k2=-1:p:1;f2= (k2+1).*(k20)+(-k2+1).*(k2=0); [f,k]=sconv(f1,f2,k1,k2,p) function?[f,k]=sconv(f1,f2,k1,k2,p)3、试求下列图片的卷积波形p=0.01;k1=-0.5:p:1;f1=ones(1,length(k1));k2=0:p:2;f2= 0.5*k2; [f,k]=sconv(f1,f2,k1,k2,p) 4、试求下列图片的卷积波形p=0.01;k1=-2:p:2;f1=(k1==-2)+(k1==2);k2=-3:p:3;f2=(k2+3).*(k2-2)+(-k2-1).*(k2=-2).*(k2=-1)+(k2-1).*(k2=1).*(k2=2)+(-k2+3).*(k22);[f,k]=sconv(f1,f2,k1,k