声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1561|回复: 1

[综合讨论] 为何无法得到 lorenz 的分岔图

[复制链接]
发表于 2008-4-1 19:29 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 ChaChing 于 2010-8-31 00:56 编辑

我用论坛的程序
clear all; global r
t0=[0 150];%积分时间
y0=[1,0,0];
%bifurcation
for r=1:0.05:30 %r的变化精度
[t,y]=ode45('Lorenz',t0,y0);
[Xmax]=getmax(y);
plot(r,Xmax,'b','markersize',1); hold on
clear Xmax
end
xlabel('r'); ylabel('Xmax')

子程序1:
function dy = Lorenz(t,y)
global r
dy=zeros(3,1); dy(1)=-10*(y(1)-y(2));
dy(2)=-y(1)*y(3)+r*y(1)-y(2); dy(3)=y(1)*y(2)-8*y(3)/3;

子程序2
function [Xmax] = getmax(y)
a=length(y); j=1;
//for i=(a-1)/2:a
b=(y(i,1)-y(i-2,1))/2; c=(y(i,1)+y(i-2,1))/2-y(i-1,1);
if y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)&c==0
Xmax(j)=y(i-1,1); j=j+1;
elseif y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)
Xmax(j)=y(i-1,1)-b^2/(4*c); j=j+1;
end
end
怎么得不到分岔图,都是直线,请各位好心的高手具体写一下lorenz的分岔图.

[ 本帖最后由 eight 于 2008-4-1 20:14 编辑 ]
回复
分享到:

使用道具 举报

发表于 2010-8-30 17:22 | 显示全部楼层
回复 xiangwei27 的帖子

这个程序好像缺少一个end,具体得到什么图像我还没运行出来呢,比较慢
   
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 11:45 , Processed in 0.071046 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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