adams.rar

adams.rar

function adams（x1y1z1hn）
a=35;b=3;c=28;
[xyz]=crk（x1y1z1h3）;
for i=4:n
fxi_0=a*（y（i）-x（i））;
fxi_1=a*（y（i-1）-x（i-1））;
fxi_2=a*（y（i-2）-x（i-2））;
fxi_3=a*（y（i-3）-x（i-3））;

fyi_0=（c-a）*x（i）-x（i）*z（i）+c*y（i）;
fyi_1=（c-a）*x（i-1）-x（i-1）*z（i-1）+c*y（i-1）;
fyi_2=（c-a）*x（i-2）-x（i-2）*z（i-2）+c*y（i-2）;
fyi_3=（c-a）*x（i-3）-x（i-3）*z（i-3）+c*y（i-3）;

fzi_0=x（i）*y（i）-b*z（i）;
fzi_1=x（i-1）*y（i-1）-b*z（i-1）;
fzi_2=x（i-2）*y（i-2）-b*z（i-2）;
fzi_3=x（i-3）*y（i-3）-b*z（i-3）;

x（i+1）=x（i）+h/24*（55*fxi_0-59*fxi_1+37*fxi_2-9*fxi_3）;
y（i+1）=y（i）+h/24*（55*fyi_0-59*fyi_1+37*fyi_2-9*fyi_3）;
z（i+1）=z（i）+h/24*（55*fzi_0-59*fzi_1+37*fzi_2-9*fzi_3）;

fxi=a*（y（i+1）-x（i+1））;
fyi=（c-a）*x（i+1）-x（i+1）*z（i+1）+c*y（i+1）;
fzi=x（i+1）*y（i+1）-b*z（i+1）;

x（i+1）=x（i）+h/24*（9*fxi+19*fxi_0-5*fxi_1+fxi_2）;
y（i+1）=y（i）+h/24*（9*fyi+19*fyi_0-5*fyi_1+fyi_2）;
z（i+1）=z（i）+h/24*（9*fzi+19*fzi_0-5*fzi_1+fzi_2）;
end

figure（1）;
plot3（xyz）;
title（‘Chen等式的三维相空间吸引子‘）;
%plot（x）;
grid;
figure（2）;
plot（xy）;
title（‘Chen等式的xy平面投影图‘）;
%plot（y）;
grid;
figure（3）;
plot（xz）;
title（‘Chen等式的xz平面投影图‘）;
%plot（z）;
grid;
figure（4）;
plot（yz）;
title（‘Chen等式的yz平面投影图‘）;
grid;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1199  2012-07-11 20:14  adams\adams.m

     文件        793  2012-07-11 20:15  adams\crk.m

     目录          0  2012-07-11 20:14  adams

----------- ---------  ---------- -----  ----

                 1992                    3