Matlab 版 (精华区)

发信人: seesea (马虎学究), 信区: Matlab
标  题: [合集]Lorenz system 如何求得？
发信站: 哈工大紫丁香 (2003年01月22日21:33:10 星期三), 站内信件


────────────────────────────────────────
 xls (忙人马大哈)                     于 2002年12月29日19:38:03 星期天 说道:

Lorenz system
   x'=10*(y-x)
   y'=28*x-y-x*z
   z'=-8/3*z+x*y
如何用Matlab 求得解？
3X!
※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: 202.118.239.17]

────────────────────────────────────────
 zjliu (Robusting)                    于 2002年12月29日19:56:17 星期天 说道:

是不是得加上边值条件？
利用pde工具箱
找本书好好看看，慢慢写程序。
【 在 xls (忙人马大哈) 的大作中提到: 】
: Lorenz system
:    x'=10*(y-x)
:    y'=28*x-y-x*z
:    z'=-8/3*z+x*y
: 如何用Matlab 求得解？
: 3X!

────────────────────────────────────────
 xls (忙人马大哈)                     于 2002年12月29日20:36:36 星期天 说道:

谢谢热心的zjliu.
是不是新版的Matlab6.5中有chaos 工具箱了？
【 在 zjliu (Robusting) 的大作中提到: 】
: 是不是得加上边值条件？
: 利用pde工具箱
: 找本书好好看看，慢慢写程序。
: 【 在 xls (忙人马大哈) 的大作中提到: 】
: : Lorenz system
: :    x'=10*(y-x)
: :    y'=28*x-y-x*z
: :    z'=-8/3*z+x*y
: : 如何用Matlab 求得解？
: : 3X!

────────────────────────────────────────
 jq (冰城唐三藏)                      于 2002年12月29日20:46:23 星期天 说道:

这个用ode算法直接解不行吗？
【 在 xls (忙人马大哈) 的大作中提到: 】
: Lorenz system
:    x'=10*(y-x)
:    y'=28*x-y-x*z
:    z'=-8/3*z+x*y
: 如何用Matlab 求得解？
: 3X!

────────────────────────────────────────
 zjliu (Robusting)                    于 2002年12月29日20:54:52 星期天 说道:

我不清楚有没有chaos 工具箱，
你可以看看本版1717文章的说明
【 在 xls (忙人马大哈) 的大作中提到: 】
: 谢谢热心的zjliu.
: 是不是新版的Matlab6.5中有chaos 工具箱了？
: 【 在 zjliu (Robusting) 的大作中提到: 】
: : 是不是得加上边值条件？
: : 利用pde工具箱
: : 找本书好好看看，慢慢写程序。

────────────────────────────────────────
 skyfly (飞天)                        于 2002年12月30日12:09:59 星期一 说道:

我在一本书中看到有一段代码，是画它的图形不知对你有没有帮助。
clear
clc
A=[-8.1/3 0 0
    0 -10 10
    0 28.4 -1.1];
y=[0 1.0 0.7]';
h=0.01;
p=plot3(y(1),y(2),y(3),'.',...
    'EraseMode','none','MarkerSize',8);
axis([0 50 -25 25 -25 25]);
hold on
grid on
for n=1:1000
    A(1,3)=y(2);
    A(3,1)=-y(2);
    ydot=A*y;
    y=y+h*ydot;
    set(p,'Xdata',y(1),'YData',y(2),'ZData',y(3),...
    'color',[0.1 0 0.1]);
    drawnow;
    n=n+1;
end
【 在 xls (忙人马大哈) 的大作中提到: 】
: Lorenz system
:    x'=10*(y-x)
:    y'=28*x-y-x*z
:    z'=-8/3*z+x*y
: 如何用Matlab 求得解？
: 3X!

────────────────────────────────────────
 figure (美丽人生☆遥望☆务正业中....)  于 2002年12月30日14:22:44 星期一 说道:

好同志
！
【 在 skyfly (飞天) 的大作中提到: 】
: 我在一本书中看到有一段代码，是画它的图形不知对你有没有帮助。
: clear
: clc
: A=[-8.1/3 0 0
:     0 -10 10
:     0 28.4 -1.1];
: y=[0 1.0 0.7]';
: h=0.01;
: p=plot3(y(1),y(2),y(3),'.',...
:     'EraseMode','none','MarkerSize',8);
: axis([0 50 -25 25 -25 25]);

────────────────────────────────────────