我想使用 numjac 计算我的系统的数值雅可比行列式。
它适用于常规语法:
[Jac, FAC] = numjac(@(t,y)model(t,y,X,VAR),ttot,ytot',dydt,jac_tol,FAC);
但是我不需要完整的雅可比,而只需要关于某些 y 的导数。我可以计算完整的然后选择所需的行,但这是浪费计算时间。如何仅在需要的地方计算雅可比?
我想使用 numjac 计算我的系统的数值雅可比行列式。
它适用于常规语法:
[Jac, FAC] = numjac(@(t,y)model(t,y,X,VAR),ttot,ytot',dydt,jac_tol,FAC);
但是我不需要完整的雅可比,而只需要关于某些 y 的导数。我可以计算完整的然后选择所需的行,但这是浪费计算时间。如何仅在需要的地方计算雅可比?