matlab绘制庞加莱截面_如何用mathematica或Matlab画弹性摆的庞加莱截面

满意答案

出门不带女人

2017.10.27

采纳率：42%    等级：13

已帮助：10867人

function ex2

options=odeset('RelTol',1e-7);

tt=2*pi/1.2;

Pm=0.1;

Pa=0.2;

w=1;

phi=0;

b=1;

e=0.1;

[t,x]=ode45(@dafen,[0:tt/100:80*tt],[0,0],options,[],e,b,Pa,Pm,phi,w);

plot(x(2000:end,1),x(2000:end,2),'-') %X=x-xxx %检验options %pojialai

hold on;

i=2000:100:3000;

plot(x(i,1),x(i,2),'*');

end

function df=dafen(t,x,flag,e,b,Pa,Pm,phi0,w)

fx1 = (x(1)-b)*(x(1)>b) + (x(1)+b)*(x(1)

df=[x(2); Pa*cos(w*t + phi0) + Pm - fx1 - 2*e*x(2) ];

end

00分享举报