Matlab 语法 (1:end-1) - 解释

这些是数组的索引。

如果您使用该diff函数将导数近似为有限差分，最好记住这些公式。一阶导数的近似值，二阶精确，是

diff 函数返回向量的相邻元素之间的差异。因此，diff 函数适用于估计相邻元素之间中点的导数。

您没有提供有关您的代码/数据的足够信息，因此我将尝试用一个数字示例来说明它。

给定多项式

f = @(x) 0.2 + 25*x − 200*x.^2 + 675*x.^3 − 900*x.^4 + 400*x.^5;

一系列等距的自变量和因变量值通过以下方式生成：

x = 0:0.1:0.8;
y = f(x);

导数的除法近似仅仅是 y 差和 x 差的向量除法：

dydx = diff(y)./diff(x)

为了绘制结果图，由于在每个区间的中点估计导数，我们必须为中点的 x 值生成一个向量：

n = length(x);
xm = (x(1:n−1)+x(2:n))./2;

为了与解析解进行比较，

xa = 0:.01:.8;
ya = 25 − 400*xa + 3*675*xa.^2 − 4*900*xa.^3 + 5*400*xa.^4;
plot(xm,dydx,'o',xa,ya)

这就是您的代码正在做的事情。jerk 是加速度的导数，但它是在时间向量中每个区间的中点估计的。您的代码正在计算这些中点。