带有周期激励的非自治动力系统的Poincare截面程序
查看( 138 ) /
评论( 24 )
%具体的讨论可以见我的个人空间！
%这是duffing方程的Poincare截面的程序，经过试验和调试是没有错误的！
%这是定义的duffing方程的式子！
%要调用poincare.m!
function dy=duffing(t,y);
global det
global w
global r
dy=[y(2);y(1)-y(1)^3-det*y(2)-r*cos(w*t)];
%Duffing 方程的Poincare截面程序！
%经过试验和调试是没有错误的！
%要调用duffing.m!
function poincare
global det
global r
global w
det=0.25;w=1;
r=0.3;
x0=[0;0];
tspan=[0:0.1*2*pi:50000];
[t,y]=ode45('duffing',tspan,x0);
m=zeros(7000,1);
n=zeros(7000,1);
for i=1:7000
m(i,1)=y(10*i,1);
n(i,1)=y(10*i,2);
end
plot(m,n,'k.','markersize',1);
[本帖最后由 liliangbiao 于 2008-4-24 15:46 编辑]
untitled.PNG
%这是duffing方程的Poincare截面的程序，经过试验和调试是没有错误的！
%这是定义的duffing方程的式子！
%要调用poincare.m!
function dy=duffing(t,y);
global det
global w
global r
dy=[y(2);y(1)-y(1)^3-det*y(2)-r*cos(w*t)];
%Duffing 方程的Poincare截面程序！
%经过试验和调试是没有错误的！
%要调用duffing.m!
function poincare
global det
global r
global w
det=0.25;w=1;
r=0.3;
x0=[0;0];
tspan=[0:0.1*2*pi:50000];
[t,y]=ode45('duffing',tspan,x0);
m=zeros(7000,1);
n=zeros(7000,1);
for i=1:7000
m(i,1)=y(10*i,1);
n(i,1)=y(10*i,2);
end
plot(m,n,'k.','markersize',1);
[本帖最后由 liliangbiao 于 2008-4-24 15:46 编辑]
untitled.PNG
TAG:
-
liliangbiao 发布于2008-04-24 15:49:52
-
其他图像
以下这些图像都可以由这个程序生成！
pioncare1.PNG
poincare2.PNG
poincare3.PNG
poincare4.PNG
-
octopussheng 发布于2008-04-24 15:52:44
-
如果是多频激励的话，截面图应该如何截取呢？
-
-
修改后的程序
修改后的程序还可以生成比较精确的分岔图程序！
1.PNG
2.PNG
23.PNG
-
无水1324 发布于2008-04-24 22:19:34
-
回复 3楼 的帖子
这个问题原来讨论过哈，不过还是没有解决，你现在做这个嘛？
-
-
有些东西需要用这个算一下，无水有建议没？
-
-
回复 6楼 的帖子
直接选取一个截面做
-
-
直接选取一个截面做，能否举个例子讲一下！
-
-
回复 8楼 的帖子
哦晕哈，我说的选取截面就是在相空间选取一平面去截相轨迹，求交点，做Poincare这个你不会不知道吧？
-
-
呵呵，原来你说的是这个意思，我之前没有理解，呵呵！
-
-
这一现象是典型的Chaotic bubbles不断裂化，向Periodic bubbles转化的过程！
-
ivy_1031发布于2008-05-21 21:24:12
-
matlab我只用来画图，对于它的编程我不是很熟！嘿嘿，所以想问问楼主，你的截面选取的原则是什么？
-
-
简单啊！！这个用的就是频闪采样法！
-
pmsl 发布于2008-05-26 15:26:50
-
为什么采样周期为10——“y(10*i,2);”
-
notemple发布于2008-07-03 17:39:52
-
请问2楼的几个图怎么改的？
再请教一个基本问题，Poincare截面图和相图（Phase portraits）有什么区别？
新人请教，谢谢
-
shipo507发布于2008-07-13 19:25:05
-
不懂
主动磁悬浮轴承--转子系统的poincare截面图是不是和这一样或差不多，请问各位大哥有源程序吗？小弟急需ing，谢谢
-
咕噜噜 发布于2008-07-13 19:57:18
-
回复 16楼 的帖子
不同的系统截面图是不一样的，不过程序类似
可以参考此帖和论坛其他帖子
-
zzglyp发布于2008-08-04 17:49:52
-
地震荷载
如果外部激励是地震荷载会有什么现象？
-
-
试用时不能通过
在matlab里面试用你的程序，提示有错误，是什么原因啊
??? Input argument "y" is undefined.
Error in ==> duffing at 5
dy=[y(2);y(1)-y(1)^3-det*y(2)-r*cos(w*t)];
-
-
这个没有错，按照提示分别保存在work里面！
-
-
回复 2楼 的帖子
如何形成这些图形，楼主可否明示，谢谢