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

yanzi

以下是我摘自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=[x(2);F*sin(w*t)-mu*x(2)+x(1)/2-x(1)^3/2];

function xiangtu1

x0=[0.9;0.01];
tspan=[0:0.01*2:100*pi];
[t,x]=ode45('duffing',tspan,x0);

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

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

无水1324

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

yanzi

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

chenvy

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

无水1324

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

yanzi

非常感谢大家的帮忙

yanzi

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

linyxff

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

20052354

zgegezenmesuna