zhangdexuan 发表于 2007-7-1 09:26

急求关于Duffing方程的庞加莱截面生成程序

急求关于Duffing方程的庞加莱截面图生成程序,如果谁能提供给我
将给您一定的Q币作为报酬,谢谢
我的QQ402559853

[ 本帖最后由 无水1324 于 2007-7-14 11:39 编辑 ]

zhangdexuan 发表于 2007-7-1 09:26

回复 #1 zhangdexuan 的帖子

急用哦,欢迎大家帮忙

无水1324 发表于 2007-7-1 09:49

借用网上一个网友的程序改了一下，但是不是你要的，你可以照着写


function df=dafen(t,x,flag,force)
df=;

％％％Poincare截面绘制
clear
ff=0.222;
options=odeset('RelTol',1e-7);
tt=3*pi/2;
=ode45(@dafen,,,options,[],ff);
plot(x(100:100:end,1),x(100:100:end,2),'.')
%pojialai
hold on
i=2000:100:3000
plot(x(i,1),x(i,2),'*')

dc1hawk 发表于 2007-7-1 10:16

思路很简单，在稳态状态下，每隔周期（自己选择）取点，作图就成了

zhangdexuan 发表于 2007-7-1 10:27

回复 #3 无水1324 的帖子

不是的哦,这个是生成相图的程序

zhangdexuan 发表于 2007-7-1 10:29

回复 #4 dc1hawk 的帖子

% Poincare_section[绘制庞加莱截面图]
betaa=0.25;
F=1.093;
v=2/3;
Poin=inline(['[x(2);',...
'-2*betaa*x(2)-x(2).^2.*sin(x(1))+F*cos(v*t);',...
'v]'],...
't','x','flag','betaa','F','v');

% Poincare_section[绘制庞加莱截面图]
=ode45(Poin,,,[],betaa,F,v);
x(:,2)=mod(x(:,2),2*pi)-pi;

phi0=pi*2/3; % 选择phi=2*pi/3这个截面
for k=1:round(max(x(:,3))/2/pi);
d=x(:,3)-(k-1)*2*pi-phi0;
=sort(abs(d));
x1l=x(K(1),1);
x1r=x(K(2),1);
x2l=x(K(1),2);
x2r=x(K(2),2);
x3l=x(K(1),3);
x3r=x(K(2),3);
if abs(P(1))+abs(P(2))<3e-16;
X1(k)=x1l;
X2(k)=x2l;
else
Q=polyfit(,,1);
X1(k)=polyval(Q,(k-1)*2*pi-phi0);
Q=polyfit(,,1);
X2(k)=polyval(Q,(k-1)*2*pi-phi0);
end
end
plot(X1,X2,'.');
xlabel('\theta','fontsize',14);
ylabel('d\theta/dt','fontsize',14);


问题：上面的程序中运用的方程是% 单摆方程[不显含时间t的自治系统]
% 方程如下:
% dθ/dt=ω,
% dω/dt=-2*β*-ω^2*sin(θ)+F*cos(vt)
% dψ/dt=v
betaa=0.25;
F=1.093;
v=2/3;
P2=-2*betaa*x(2)-x(2).^2.*sin(x(1))+F*cos(v*t);
dx=;


现在我想把这单摆脱方程换成达芬方程如下所示
function df=dafen(t,x,flag,force)
df=;
即

帮我换进去可以吗？程序都是现成的

无水1324 发表于 2007-7-1 10:32

原帖由 zhangdexuan 于 2007-7-1 10:27 发表 http://www.chinavib.com/forum/images/common/back.gif
不是的哦,这个是生成相图的程序


你看一下红色字体部分

mechanic05 发表于 2008-1-11 16:31

回复 3楼 的帖子

新年好！:handshake
请问这个程序中为什么tt=3*pi/2.一般计算周期的选择的原则是什么？

mechanic05 发表于 2008-1-11 16:33

回复 7楼 的帖子

取 tt=２＊pi/1.2 行吗？

xiaoqiu810818 发表于 2008-1-11 16:37

不错

octopussheng 发表于 2008-1-14 08:01

回复 9楼 的帖子

可以的！

咕噜噜 发表于 2008-1-14 09:24

这个截面是自己任意选取的

无水1324 发表于 2008-1-14 13:32

回复 10楼 的帖子

xiaoqiu810818可以参考这个程序了

xiaoqiu810818 发表于 2008-2-19 16:25

回复无水

谢谢无水了,自治系统和非自治系统我想截面的选取应该不一样吧,如Rossler等.

无水1324 发表于 2008-2-19 19:39

回复 14楼 的帖子

他们之间有差异
页: [1] 2
查看完整版本: 急求关于Duffing方程的庞加莱截面生成程序