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

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

我正在使用 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？