jiemianz.m

此程序用于画混沌系统的相图及其庞加莱截面

[T,Y]=ode45(@chen,[0:0.01:50],[0.1 0.1 0.21 0.31]);
figure
plot(Y(:,1),Y(:,3));
z=-40;
T=1;%厚度
k=1;
z1=z+T;
q=Y(:,3);
m=length(q);
flag=1;
for i=1:m
    if(q(i)>z&q(i)<z1)
        p(k,1)=Y(i,1);
        p(k,2)=Y(i,2);
        k=k+1;
        flag=0;
    end
end
if(flag)
   disp('截面上没有切点，请重新输入x，y，z轴截距：')
   break
end
a=p(:,1);
b=p(:,2);
figure
plot(a,b,'.')
%plot(p(:,1),p(:,2),'.')