复杂非线性系统的数值微分

Question

我有一个非线性稳定空间矩阵。我需要解决差分，如下图所示： 我在pic上解释了更多（这里dD/dx）

Answer 1

您可以使用MATLAB Symbolic Toolbox进行符号化差异，但您必须告诉它执行正确的偏导数。因变量（x）有三个分量，因此不能相对于x取D的偏导数，但可以采用x1，x2和x3的偏导数。

 syms t x1 x2 x3  x1_dot x2_dot x3_dot

 x1_dot = diff(x1, t)
 x2_dot = diff(x2, t)
 x3_dot = diff(x3, t)
 D(x1,x2,x3) = [... your function matrix ... ]
 dD_x1 = diff(D,x1)
 dD_x2 = diff(D,x2)
 dD_x3 = diff(D,x3)