连续时间信号卷积运算及MATLAB实现,连续时间信号卷积运算的MATLAB实现

连续时间信号卷积运算的MATLAB 实现

一、实验目的

(1)理解掌握卷积的概念及物理意义。

(2)理解单位冲激响应的概念及物理意义。

二、实验原理

根据前述知识，连续信号卷积元素按定义为

()()()()()1212f t f t f t f f t d τττ

∞

-∞

=*=

-?

卷积计算可以通过信号分段求和来实现，即

()()()()()()()1212120

lim

k f t f t f t f f t d f k f t k τττ∞

∞

?→=-∞

-∞

=*=

-=??-???

∑

?

如果只求当t n =?(n 为整数)时()f t 的值()

f n ?，则由上式可得

()()()()()1212k k f n f k f n k f k f n k ∞

∞

=-∞

=-∞

?=

???-?=??

??-??????

∑

∑

式中的()()

12k f k f n k ∞

=-∞???-?∑

实际上就是连续信号

()

1f t 和

()

2f t 经等时间间隔?均匀抽样的离散序列

()1f k ?和

()

2f k ?的卷积和。当?足够小时，

()

f n ?就是卷积积分的结

果——连续时间信号()

f t 的较好的数值近似。

三、实验内容

1、实验参考程序 以下是MA TLAB 实现连续信号卷积的通用函数sconv()。 function [f,k]=sconv(f1,f2,k1,k2,p) f=conv(f1,f2); f=f*p;

k0=k1(1)+k2(1);

k3=length(f1)+length(f2)-2; k=0:p:k3*p; subplot(2,2,1) plot(k1,f1) title('f1(t)') xlabel('t') ylabel('f1(t)') subplot(2,2,2) plot(k2,f2) title('f2(t)') xlabel('t')