来自 Matlab 在 Maxima 上的等效 lsim 函数答案

【问题标题】：Equivalent lsim function from Matlab on Maxima来自 Matlab 在 Maxima 上的等效 lsim 函数
【发布时间】：2018-05-04 01:45:32
【问题描述】：

所以我试图在 Maxima 上找到 lsim 函数的等效项，并试图重现此示例中的结果：https://www.mathworks.com/help/control/ref/lsim.html。我正在使用 COMA 库，定义了一个方波函数，并将其乘以与上面链接的第一个示例相同的传递函数，使用 plot2d，但结果与预期完全不同。有人知道 lsim 在 Maxima 上使用的类似表格吗？

我的千里马代码：

load(coma);

h(t) := (1 + signum(t))/2;
unitBox(t) := h(t+1/2) - h(t - 1/2);
squareWave(t, t0, p, d) :=
       unitBox(mod((t-t0)/p,1)/2/d);
ft(t) := ilt((s-1)/(s^2+s+5), s, t);

resp(t) := ft(t)*squareWave(t,0,4,0.5);
plot2d(
  resp(t),
  [t,0,15]
);

【问题讨论】：

标签： matlab simulation maxima

【解决方案1】：

我假设ft 代表单位脉冲响应，而您希望resp 是对方波的响应，对吗？如果是这样，你不想将resp定义为单位脉冲响应和方波的卷积吗？

我对 COMA 不是很熟悉，所以我必须查找如何计算卷积。如果您同意这是需要的，那么我们可以解决这个问题。

顺便说一句，如果你能准确地展示实际和预期的结果，而不是假设他们会解决这个问题，其他人会更好地理解问题。我能够弄清楚，但总的来说，您可以通过明确说明预期和实际情况，让其他人更容易帮助您。

【讨论】：