yanzi 发表于 2009-4-19 21:30

关于duffing方程相图程序的问题

以下是我摘自liliangbiao 空间求duffing方程的相图程序，帮我解释一下

r=length(t)-3000:length(t);这句是什么意思？
非常感谢！！

% Author: Thomas Lee
% E-mail: lixf1979@126.com
% Corresponding: School of Mathematics, Physics and Software Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China



function dx=duffing(t,x)；

mu=0.168;

%F=0.22;

%F=0.198;
%F=0.18;
F=0.15;
w=1;
dx=;

function xiangtu1

x0=;
tspan=;
=ode45('duffing',tspan,x0);

r=length(t)-3000:length(t);

plot(x(r,1),x(r,2),'r.','markersize',1)

无水1324 发表于 2009-4-20 08:23

r=length(t)-3000:length(t);
就是取积分后面的3000点画图

yanzi 发表于 2009-4-21 21:10

回复 沙发 无水1324 的帖子

那个点数是不是可以任意改动？ 本人愚钝

chenvy 发表于 2009-4-21 21:55

回复 板凳 yanzi 的帖子

应该可以，去掉前面一些不规则点就行
再者，可能是为了图像清晰一点吧

无水1324 发表于 2009-4-21 22:00

回复 地板 chenvy 的帖子

恩，楼上的解释是完全正确的

yanzi 发表于 2009-4-22 15:53

非常感谢大家的帮忙

yanzi 发表于 2009-5-5 21:56

回复 6楼 yanzi 的帖子

最近发现如果这个点数变化的话，相图就会有所变化，而且点数越大圈数越多，不明白到底取多少是正确的？？大家帮帮忙吧 非常感谢！！

linyxff 发表于 2009-5-6 12:51

回复 7楼 yanzi 的帖子

我也遇到这个问题，我的理解是此时参数影响下，相图不是周期或2周期等，已经呈现混沌或准周期，所以取点越多图越密集，只要取合适点保证相图清晰反映现象就可以了。可能解释不对，请高手确认一下~~O(∩_∩)O

20052354 发表于 2009-6-8 19:10

回复 楼主 yanzi 的帖子

zgegezenmesuna
页: [1]
查看完整版本: 关于duffing方程相图程序的问题