signal-processing - 用倍频程在时域中绘制卷积

Question

这个绘图但 conv 的结果是一个新长度的向量，因此t在绘图中包含plot(t, z1) %doesn't work!.

t = [-5:.1:10];
unit = @(t) 1.*(t>=0);

h1 = @(t) (3*t + 2).*exp(-3*t).*unit(t);
z1 = conv(unit(t), h1(t));

plot(z1);

我想要一个作为时间函数的卷积信号图。

Answer 1

您需要添加shape参数。这是规格：

— 函数文件：conv (a, b) — 函数文件：conv (a, b, shape) 对两个向量 a 和 b 进行卷积。

输出卷积是一个长度等于length (a) + length (b) - 1的向量。当a和b是两个多项式的系数向量时，卷积表示乘积多项式的系数向量。

可选的形状参数可能是

shape = "full" 返回全卷积。（默认） shape = "same" 返回卷积的中心部分，大小与 a 相同。

所以像这样卷积：

z1 = conv(unit(t), h1(t), "same");

您将获得与原版相同的时间单位。