Octave符号包计算雅可比矩阵报错：下标索引必须为整数或布尔值

错误根源 #

你调用Octave symbolic包的jacobian函数时，未指定求偏导的目标变量。原代码最后一行jacobian([xdot;ydot])缺少变量参数，导致函数无法正确解析求导对象，进而触发下标错误。

修正后的代码 #

pkg load symbolic

syms x y

xdot = -x + 2*x^3 + y;
ydot = -x - y;

[xeq,yeq] = solve(xdot==0,ydot==0);

xeq = double(xeq);
yeq = double(yeq);

// 补充求导变量参数，明确对x和y求偏导
jacobian_matrix = jacobian([xdot;ydot], [x,y]);

后续扩展：计算平衡点处的雅可比数值 #

如果需要在每个平衡点处得到雅可比矩阵的具体数值，可通过subs函数代入平衡点坐标：

// 遍历所有平衡点，计算对应雅可比矩阵数值
for i = 1:length(xeq)
    jacobian_num = subs(jacobian_matrix, {x, y}, {xeq(i), yeq(i)});
    jacobian_num = double(jacobian_num);
    disp(["平衡点(", num2str(xeq(i)), ",", num2str(yeq(i)), ")处的雅可比矩阵："]);
    disp(jacobian_num);
end

内容的提问来源于stack exchange，提问作者yahskA