资源描述：
《信号与系统课件实验 2.连续时间信号卷积和离散时间信号卷积.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、二．连续时间信号卷积和离散时间信号卷积1．连续时间信号卷积；2．离散时间信号卷积连续时间信号卷积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:取样时间间隔f1=0.5*(0:0.01:2);f2=0.5*(0:0.01:2);k1=0:0.01:2;k2=0:0.01:2;p=0.01;f=conv(f1,f2
2、);%计算序列1与序列2的卷积和f=f*p;k0=k1(1)+k2(1);%计算序列f非零样值的起点位置k3=length(f1)+length(f2)-2;%计算卷积和f非零样值得宽度k=k0:p:k0+k3*p;%确定卷积和f非零样值的时间向量subplot(3,3,1)plot(k1,f1)%在子图1绘制f1(t)时域波形图title('f1(t)')xlabel('t')ylabel('f1(t)')subplot(3,3,4)plot(k2,f2)%在子图2绘制f2(t)时域波形图title('f2(t)')xlabel('t')ylabel('f2(t)')subplot(3,3
3、,7)plot(k,f);%画卷积f(t)的时域波形h=get(gca,'position');3h(3)=2.5*h(3);set(gca,'position',h)%将第三个子图的横坐标范围扩为原来的2.5倍title('f(t)=f1(t)*f2(t)')xlabel('t')ylabel('f(t)')离散时间信号卷积function[f,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=[1,2,1];f2=ones(1,5);k1=[-101];k2=-2:2;f=conv(f1,f2)%计算序列f1与f2的卷积和fk0=k1(1)+k2(1);%计算序列f非零样值的起点位置k3=length(f1)+length(f2)-2;%计算卷积和f的非零样值的宽度k=k0:k0+k3%确定卷积和f非零样值得序号向量subplot(3,3,1)stem(k1,f1)%在子图1绘制序列f1(k)时域波形图title('f1(n)')xlabel('n')ylabel('f1(n)')subplot(3,3,4)s
5、tem(k2,f2)%在子图2绘制序列f2(k)时域波形图3title('f2(n)')xlabel('n')ylabel('f2(n)')subplot(3,3,7)stem(k,f)%在子图3绘制序列f(k)时域波形图title('f1(n)与f2(n)的卷积和f(n)')xlabel('n')ylabel('f(n)')h=get(gca,'position');h(3)=2.5*h(3);set(gca,'position',h)%将第三个子图的横坐标范围扩为原来的2.5倍3
