如何利用Wolfram Mathematica画庞加莱截面图

EquationTrekker=_=好老的包了

微分方程的庞加莱截面嘛，可以用NDSolve、WhenEvent、Sow、Reap四个指令结合获得数据，再用ListPlot，我记得好像WhenEvent的帮判漏助段困文档里就举了庞握冲念加莱截面的例子。

你去看看帮助文档，大概能自己想明白的，有问题再问我。

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)<态尘-b)

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

end

你的问题在网页链接已经答过了。

经常看不懂报错，就换个中文版，或者学好英语。

一行字让蠢就七个词，一个歼滑搜少见的词都没有，看不懂？恐怕只是没这氏历个意愿吧。要是知道答案就在报错里的话，恐怕你也不会那么回避它了，哈哈……

---