声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1752|回复: 0

[稳定性与分岔] 我运行sixchaos.m文件不出图像是什么原因

[复制链接]
发表于 2016-6-2 15:51 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
以下是函数sixchaos的sixchaos.m文件
function dxdt=sixchaos(t,x)
%定义一个六维混沌系统
A=-15;
B=20;
C=-20;
D=20;
E=-5;
F=-7;
G=12;
H=-25;
I=-20;
J=10;
K=5;
dxdt=zeros(6,1);
dxdt(1)=A*x(1)+B*x(6)+x(2).*x(3).*x(4).*x(5).*x(6);
dxdt(2)=C*(x(1)+x(2))-x(1).*x(3).*x(4).*x(5).*x(6);
dxdt(3)=D*(x(2)-x(5))+E*x(3)+x(1).*x(2).*x(4).*x(5).*x(6);
dxdt(4)=F*x(4)+G*x(5)+H*x(6)-x(1).*x(2).*x(3).*x(5).*x(6);
dxdt(5)=-x(4)+I*x(5)+x(1).*x(2).*x(3).*x(4);
dxdt(6)=G*x(5)+K*x(6)-x(1).*x(2).*x(3).*x(4);
以下是sixchaosdriver.m文件。
options=odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-4 1e-4 1e-4 1e-4]);
[t,x]=ode45(@sixchaos,[0 200],[0.02,0.01,0.03,0.04,0.05,0.06],options);

figure(1)
subplot(231);
plot(x(:,1),x(:,1));
title('x(1)-x(2)平面相图')

subplot(232);
plot(x(:,1),x(:,3));
title('x(1)-x(3)平面相图')

subplot(233);
plot(x(2,:),x(3,:));
title('x(2)-x(3)平面相图')

subplot(234);
plot(t,x(:,1));
title('x(1)时域波形')

subplot(235);
plot(t,x(:,2));
title('x(2)时域波形')

subplot(236);
plot(t,x(:,3));
title('x(3)时域波形')

figure(2)
plot3(x(:,1),x(:,2),x(:,3));
title('x(1)-x(2)-x()3相图')


不知道是什么原因在运行sixchaos.m文件时不出图像?

回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-9-20 15:04 , Processed in 0.054510 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表