matlab连续时间信号卷积和离散时间信号卷积程序

1、 二连续时间信号卷积和离散时间信号卷积连续时间信号卷积f(t)=f(t)f(t)；12离散时间信号卷积f(n)=f(n)f(n)12连续时间信号卷积functionf,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:取样时间间隔f1=0.5*(0:0.01:2);f2=0.5*(0:0.01:2);k1=0:0.01:2;k2=0:0.0

2、1:2;p=0.01; f=conv(f1,f2);积和f=f*p;k0=k1(1)+k2(1);起点位置k3=length(f1)+length(f2)-2;度k=k0:p:k0+k3*p;的时间向量subplot(3,3,1)plot(k1,f1)形图title(f1(t)xlabel(t)ylabel(f1(t)subplot(3,3,4)plot(k2,f2)形图title(f2(t)xlabel(t)ylabel(f2(t)subplot(3,3,7)plot(k,f);h=get(gca,position);%计算序列1与序列2的卷%计算序列f非零样值的%计算卷积和f非零样值得宽%

3、确定卷积和f非零样值%在子图1绘制f1(t)时域波%在子图2绘制f2(t)时域波%画卷积f(t)的时域波形 # %将第三个子图的横坐标范围h(3)=2.5*h(3);set(gca,position,h)扩为原来的2.5倍title(f(t)=fl(t)*f2(t)xlabel(t)ylabel(f(t)吐) # #离散时间信号卷积functionf,k=dconv(f1,f2,k1,k2)%Thefunctionofcomputef=f1*f2%f:卷积和序列f(k)对应的非零样值向量%k:序列f(k)的对应序号向量%f1:序列f1(k)非零样值向量%f2:序列f2(k)非零样值向量%k1:

4、序列f1(k)的对应序号向量%k2:序列f2(k)的对应序号向量%计算序列f1与f2的卷积和f%计算序列f非零样值的起点位置%计算卷积和f的非零样值的宽度%确定卷积和f非零样值得序号向量f1=1,2,1;f2=ones(1,5);k1=-101;k2=-2:2;f=conv(f1,f2)k0=k1(1)+k2(1);k3=length(f1)+length(f2)-2;k=k0:k0+k3%在子图1绘制序列f1(k)时域波形图subplot(3,3,1)stem(k1,f1)title(f1(n)xlabel(n)ylabel(f1(n)subplot(3,3,4) stem(k2,f2)title(f2(n)xlabel(n)ylabel(f2(n)subplot(3,3,7)stem(k,f)title(f1(n)与f2(n)的卷积和f(n)xlabel(n)ylabel(f(n)h=get(gca,position);h(3)=2.5*h(3);set(gca,posi