用八度音阶在时域中绘制卷积

Question

此绘图但转化的结果是新长度的向量，因此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

您需要添加形状参数。这是spec：

  

- 函数文件：conv（a，b） - 函数文件：conv（a，b，shape）   卷积两个向量a和b。

     

输出卷积是长度等于长度（a）+的向量   长度（b） - 1.当a和b是2的系数向量时   多项式，卷积表示系数向量   乘积多项式。

     

可选的shape参数可能是

     

shape =“full”返回完整的卷积。 （默认）shape =“相同”   返回卷积的中心部分，其大小与a。相同。

如此卷入：

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

你将得到与原版相同的时间单位。