如何在 Octave 上绘制函数微分图？

Question

我从一开始就想学习三角学。 所以我可以通过将这些命令作为命令来绘制曲线。

x1 = linspace(-2*pi, 2*pi)
y1 = cos(x1)
plot(x1, y1)

但如果我尝试这样做，

x1 = linspace(-2*pi, 2*pi)
y2 = diff(cos(x1))
plot(x1, y1) 

我会收到一条错误消息，指出

error: __plt2vv__: vector lengths must match error: called from     __plt__>__plt2vv__ at line 487 column 5     __plt__>__plt2__ at line 247 column 14     __plt__ at line 112 column 18     plot at line 229 column 10

“但是，我确实获得了 y2 = diff(cos(x1)) 的值，由于此错误，我无法绘制任何内容。”

注意：我知道 sin(x) 的微分是 cos(x)，我可以只使用 cos(x)，但是不能'我这样做吗？

尝试了 diff(cos(90)) 我得到了

ans = \[\](0x0)

Answer 1

diff 比其输入短一个元素，因为它是差异 [v(2)-v(1), v(3)-v(2), ...]< /代码>。

您需要执行类似 y2 = [NaN, diff(cos(x1))]; 的操作来将数组填充到与 x1 相同的大小。

请注意，这当然会导致偏移，您需要一个 中心差分 方法而不是前向差分方法（如 diff）来保持事物的良好对齐。