9.2 连续时间信号卷积运算的 MATLAB 实现
一、实验目的
(1)理解掌握卷积的概念及物理意义。
(2)理解单位冲击响应的概念及物理意义。
二、实验原理
根据前述知识,连续信号的卷积运算定义为
(9-
1)
卷积计算可以通过信号分段求和来实现,即
(9-2)
如
果只求当时的值,则由上式可得
(9-3)
式
(9-3) 中 的 实 际 上 就 是 连
续信号经等时间间隔均匀抽样的离散序列的卷积和。当足够小时,就是卷积积
分的结果——连续时间信号的较好的数值近似。
三、实验内容与方法
用 MATLAB 实现连续信号 f
1
(t)和 f
2
(t)卷积的过程如下:
(1)将连续信号 f
1
(t)和 f
2
(t)以时间间隔 Δ 进行取样,得到离散序列
f
1
(kΔ)、f
2
(kΔ);
(2)构造 f
1
(kΔ)、f
2
(kΔ)与相对应的时间向量 k1 和 k2;
(3)调用 conv()函数计算卷积积分 f(t)的近似向量 f(nΔ);
(4)构造 f(nΔ)对应的时间向量 k。
下面是利用 MATLAB 实现连续信号卷积运算的通用函数 sconv(),该程序
在计算出卷积积分的数值近似值的同时,还绘出 f(t)的波形图。需要注意的是,
程序中是如何构造 f(t)的对应时间向量 k 的?另外,程序在绘制 f(t)波形图采用
的是 plot 命令而不是 stem 命令。
function [f,k]=sconv(f1,f2,k1,k2,p)