Matlab code for poltting the Poicanre sections of autonomous chaotic systems

下面的这个程序是修改于某个博客上的程序，若侵犯了版权，请及时告知，以便我及时修改！谢谢！振动资讯(B-q-v1]N'i v H"\
% M File for systems, which can be modified by any autonomous chaotic systems, should be saved in the Work file of your Matlab!!!!!!!振动资讯!w OqG%kzE
function dx=Chu(t,x);
"r)| nK b"G}}0% Chu system [A new autonomous chaotic system, which has many similar properties and dissimilar properties with Lorenz families]
eT0B @X8X7I0% dx=Chu(t,[x;y;z;a;b;c])振动资讯#V#A+g|Hk_g
% t-time，x,y,z- Parameters，a,b,c-coefficients

L0S%VJ `v`0% 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,
ba!L[8{ivOf6{6W0%some typical are given in the listed paper below.
+n3G:q"h j]7NWF/l_0%% Please refer one of our published papers as follows:振动资讯h#M,rJ9X0V2QmW j
%褚衍东， 李险峰， 张建刚，常迎香.一个新自治混沌系统的计算机仿真与电路模拟.
DtC4m {1o's0%《四川大学学报：自然科学版》-2007年44卷3期 -596-602.
Z.teV/P|L,|l0dx(1,1)=x(4)*(x(2)-x(1));
^8@Fi;G k*X$I*T0dx(2,1)=x(1)*x(3)-x(2);振动资讯 uO7U)Hn4a
dx(3,1)=x(5)-x(1)*x(2)-x(6)*x(3);振动资讯:C!t\;c5h?
dx(4,1)=0;振动资讯3iJA'H)f,aK8|*h|
dx(5,1)=0;
y\i7t-}0dx(6,1)=0;振动资讯)f,}7Ho cr"c['~8g
振动资讯B.ofO?S)R ]9k

0q#m~$R6qvYiGv0% Another M File, which can be pasted in the Command windows or saved in the Work file of your Matlab.
8g]7HR8t0振动资讯BLE1jR-Gbt
Z=[];
-itE3~7Zv;a*Y0[T,Y]=ode45('Chu',[0,5000],[0.1;0.1;0.1;5;16;1]);
k6W3\q Xy;d0for k=1:length(Y)振动资讯h9Rm1@"SF
if 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,
l'oS/p8~
}0%which can be selected any other Poincare section if you like, e.g. x,y,z=0 or others!
3n`?}wqwE0Z=[Z Y(k,1)+i*Y(k,2)];
@$]N%FS D"r~P@${0end
Fr Yele(q:T0end振动资讯x[ U4ruY
plot(Z,'.','markersize',2)振动资讯"m"J,l9]9]1J-Y N
title('Chu system's Poincare map on the plane of z=0')
0n'S4?(fDrq s0xlabel('x'),ylabel('y')'振动资讯7p^9^-g {#~ i
振动资讯 L1~F8W [0WI\
untitled.png