computebase1.asv

计算动力学系统的分岔图

function  [Q0,evl,NSub] =  computeBase(A0, unstable_flag)

if unstable_flag
    
[V1,D1]=eig(A0);
%[i,j]=find(abs(D1)>1);
  NU=size(find(abs(D1)>1),1); 
  EIGS(A,B,NU)
%  [QU,RU]=schur(A0);
%  Q0=QU;
%  NSub=NU;
 %Q0=QU(:,1:Nsub);
    
else    
% Compute all eigenvalues and eigenvectors and check number of unstable  ones
[V1,D1]=eigs(A0);
 NU=size(find(abs(D1)>1),1); 
 [QS,RS]=schur(A0);
 Q0=QS;
 %Q0=QS(:,1:Nsub);
 NSub=NU
 
%    
% [Q1,R1]=schur(A1);
% [Q2,R2]=schur(Anm);
% Q1C=Q1(:,lu1+1:nphase);
% QnmC=Q2(:,lunm+1:nphase);
% x=X0(:,1,1);
% for i=1:n
% x1=feval(fpmds.func,0,x,p{:});
% end
% xx(:,1)=x-X0(:,:,1);  %f^i(x1)-x1=0
% 
% x=X(:,:,nmesh);
% for i=1:n
% x=feval(fpmds.func,0,x,p{:});
% end
% xx(:,nphase)=x-X0(:,:,nphase); %f^i(xn)-xn=0
% 
% for i=2:nmesh-1
%     x=X0(:,:,i);
%     for i=1:n
%      x=feval(fpmds.func,0,x,p{:});
%     end
%     xx(:,i)=x-X0(:,:,i);%f^i(xi-1)-xi=0
% end
% 
% if unstable_flag