process_mfoldns.m

计算动力学系统的分岔图

 
 function result=process_mfoldNS(x,jac)
%
%
global cds lpmds
[x1,p] =rearr(x);p1=n2c(p);
nphase=size(x1,1);
nphase = size(x1,1);
jac1=jac;
% calculate eigenvalues and eigenvectors
[V,D] = eig(jac1);
% find pair of complex eigenvalues
d=diag(D); 
idx1=0;idx2=0;
    for k=1:nphase
       for j=k+1:nphase
           if  (d(k)== conj(d(j))) &(imag(d(k))>0.0001)% & (abs(1-d(k)*d(j))<0.0000001)
               idx1=k;
               idx2=j;
           end
            end
        end
if idx1==0
       result='Neutral saddle';
       return;
end 
result=real(idx1);
return;
% ---------------------------------------------------------------
function [x,p] = rearr(x0)
% [x,p] = rearr(x0)
% Rearranges x0 into coordinates (x) and parameters (p)
global cds lpmds
p = lpmds.P0;
p(lpmds.ActiveParams) = x0(end-1:end);
x = x0(1:end-2);