卷积的物理意义

卷积的物理意义

广东梅州市梅县区-观测场

1

二．连续时间信号卷积

1．连续时间信号卷积)()()(

21

tftftf；

连续时间信号卷积

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);%计算序列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,7)

plot(k,f);%画卷积f(t)的时域波形

h=get(gca,'position');

h(3)=2.5*h(3);

set(gca,'position',h)%将第三个子图的横坐标范围

2

扩为原来的2.5倍

title('f(t)=f1(t)*f2(t)')

xlabel('t')

ylabel('f(t)')