ode_lin.m

来自「关于混沌系统的李氏指数计算等混沌系统中重要参数计算的代码」· M 代码 · 共 16 行

function [f] = ode_lin( t, X, P, n, neq, n_exp )
  f=zeros(neq,1);
 f(1,1)=-P(3)*(X(1)-X(2));
 f(2,1)=-X(1)*X(3)+P(1)*X(1)-X(2);
 f(3,1)=X(1)*X(2)-P(2)*X(3);
 Jcb=zeros(3);
 Jcb(1,1)=-P(3);
 Jcb(1,2)=P(3);
 Jcb(1,3)=0;
 Jcb(2,1)=-X(3)+P(1);
 Jcb(2,2)=-1;
 Jcb(2,3)=-X(1);
 Jcb(3,1)=X(2);
 Jcb(3,2)=X(1);
 Jcb(3,3)=-P(2);
 for k=1:n_exp,  f(k*n+1:(k+1)*n,1)=Jcb*X(k*n+1:(k+1)*n,1); end;