Matlab code for poltting the Poicanre sections of autonomous chaotic systems

下面的这个程序是修改于某个博客上的程序，若侵犯了版权，请及时告知，以便我及时修改！谢谢！振动资讯rvFYdF
Y
% M File for systems, which can be modified by any autonomous chaotic systems, should be saved in the Work file of your Matlab!!!!!!!振动资讯 s~s,@o
function dx=Chu(t,x);
{Evx5P"Cui%o0% Chu system [A new autonomous chaotic system, which has many similar properties and dissimilar properties with Lorenz families]振动资讯%wa&BxP/DiOB
% dx=Chu(t,[x;y;z;a;b;c])
U+xM}U*};f4{[X"N0% t-time，x,y,z- Parameters，a,b,c-coefficients
l g'i1I%bcP4H,|H0% eg: dx=Chu(0,[0.11;0.11;0.11;5;16;1])%Of course, you can verified the other values of parameter b if you like,
}adHv/dM0Wj0%some typical are given in the listed paper below.振动资讯^bXfQ0l
%% Please refer one of our published papers as follows:
Rsv%M[0%褚衍东， 李险峰， 张建刚，常迎香.一个新自治混沌系统的计算机仿真与电路模拟.
t-W F |$Sz$q%Ik0%《四川大学学报：自然科学版》-2007年44卷3期 -596-602.振动资讯P wb7p:]C
~
lH
dx(1,1)=x(4)*(x(2)-x(1));
i1})Yzj!^V%^a0dx(2,1)=x(1)*x(3)-x(2);振动资讯T-W!K5V*Ni0J3h
dx(3,1)=x(5)-x(1)*x(2)-x(6)*x(3);振动资讯,JV3LA;k"C
dx(4,1)=0;
4CF?(f*pm$Lg+s6R0dx(5,1)=0;
%?jU&~'^+a J0dx(6,1)=0;振动资讯4rM5EN_5pCU!xX

8DD DT a4x0oyh0

H\)``|0% Another M File, which can be pasted in the Command windows or saved in the Work file of your Matlab.
.co`_%` ]0振动资讯X$Z:|P r6AP
Z=[];
4til!xz`4m4k0[T,Y]=ode45('Chu',[0,5000],[0.1;0.1;0.1;5;16;1]);
&dw]4f7`c @%r4n0for k=1:length(Y)
:w J;Yp'HS0B$Z/n p9_0if abs(Y(k,3)-0)<1e-2; % Be careful of 0, that is z=0, namely, on the plane of x-y, the Poincare sections we chose in this paper,振动资讯U1AN7xfK6v
%which can be selected any other Poincare section if you like, e.g. x,y,z=0 or others!振动资讯Hl$]_]Tqj&nn`e
Z=[Z Y(k,1)+i*Y(k,2)];
9uS w6N0H;Pr9G0end
~3LK'x4O0end
$]4jJBHl0plot(Z,'.','markersize',2)
}O$v.[_\0title('Chu system's Poincare map on the plane of z=0')振动资讯)I4c h5Zn+Q(|
xlabel('x'),ylabel('y')'振动资讯2q7_As6EALG I
振动资讯_hO7ol7Dw:U
untitled.png