MATLAB：状态矩阵的雅可比随时间使用 jacobianest

编程问答 2022-04-28

问题描述

我正在使用 DERIVESTsuite 工具箱来查找与状态数组相关的雅可比。

我有一个带输入的函数

  1. x- 状态数组 (n x m)。其中 n 是状态数，m 是时间节点数。
  2. u- 控制向量 (1 x m)
  3. t- 时间向量 (1 x m)
  4. constants- 常量结构

我目前让我的代码循环遍历每个时间节点，为每个时间节点提供一个 nxn 矩阵：

for i = length(t):-1:1    % Jacobian of state vector at each node
    A(:,:,i) = jacobianest(@(x) F(x,u(i),t(i),constants),x(:,i));
end

有没有办法在不循环每个时间节点的情况下做到这一点？我知道下面的代码给了我一个矩阵 A，它是 (nm x nm)。

A = jacobianest(@(x) F(x,u,t,x);

我可以从中提取我正在寻找的内容还是有更好的方法来使用 jacobianest？

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

derivativematlabmatlab