计算匿名函数的雅可比行列式 - MATLAB

Question

我正在尝试在 Matlab 中求解非线性颂歌系统，如下所示。

editparams %parameters
Tend = 50;
Nt = 100;


% Define RHS functions

RHS = @(t,x) ModelRHS(t,x,param); %anonymous function defining the set of equations

%Execution-----------------------------------------------------------------
x0 =[0.04;0.75;0.85]; %Initial condition

t = linspace(0,Tend,Nt); %TSPAN
[t x] = ode45(RHS, t, x0);

现在，我需要找到系统的稳定状态，我正在尝试为此创建一个函数。我想我应该使用雅可比行列式来计算稳态。我的方程位于一个匿名函数中，在下面的代码中定义为 f 。但是，我意识到 jacobian 不适用于匿名函数（或者可能有一种方法可以处理匿名函数）。所以我想我应该将匿名函数转换为符号函数并尝试一下。但 我仍然很难完成这项工作。所以非常感谢任何帮助！

function SS = SteadyState(param, E)

f = @(t,x)ModelRHS(t,x,param,E); %set of odes

SymbolicSystem = sym(f); %trying to make the anonymous function symbolic

SymbolicJacobian = jacobian(SymbolicSystem',x); %jacobian 
Jacob = matlabFunction(SymbolicJacobian,x);

end

另外，如果除了找到雅可比行列式之外还有其他方法，也请告诉我。 我尝试使用“fsolve”来计算稳态，如下所示：

f = @(t,x)ModelRHS(t,x,param);
x0 =[0.04;0.75;0.85]';

options = optimoptions('fsolve','Display','iter'); % Option to display output
SS = fsolve(f,x0,options); % Call solver

但它返回了一个错误

    Not enough input arguments.

Error in @(t,x)ModelRHS(t,x,param)

Error in fsolve (line 242)
            fuser = feval(funfcn{3},x,varargin{:});

Caused by:
    Failure in initial objective function evaluation. FSOLVE cannot continue.