|
楼主 |
发表于 2010-9-10 21:45
|
显示全部楼层
回复 咕噜噜 的帖子
我看了一下论坛里相关的帖子,内容如下以Duffing系统为例,根据前面给出的庞加莱映射定义,计算其映射(非原创)
function dx=duffing(t,X)
global F wd;
r=0.25;
x=X(1);
y=X(2);
psi=X(3);
dx=zeros(3,1);
dx(1)=y;
dx(2)=-r*y+1/2*x*(1-x^2)+F*sin(psi);
dx(3)=wd;
function duffing_poincare
global F wd;
wd=1;
F=0.28;
x0=[0;0.0001;0];
tspan=[0:0.1*2*pi:100000*pi];
[t,y]=ode45('duffing',tspan,x0);
m=[];
n=[];
for i=1000:50000
m(i,1)=y(10*i,1);
n(i,1)=y(10*i,2);
end
plot(m,n,'k.','markersize',2);
,我看就是把参数取某一具体值时所得的解取出来画图,并没有取什么截面,另外他的横纵坐标各代表什么呢,不胜感激
|
|