声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1367|回复: 1

[编程技巧] MATLAB作稳态之后的庞加莱映射

[复制链接]
发表于 2012-5-29 09:24 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
微分方程为:
dy=[x(3)
       x(4)
       0.1*sin(10*t)-2.5*x(3)+1.5*x(4)-x(1)+0.5*x(2)
       0.15*sin(10*t)-2.5*x(4)+1.5*x(3)-x(2)+0.5*x(1)];
这是我做的程序
function poincare(tstop)
>> y0=[0;0;0;0];
>> tspan=[0:1:50];
>> [t,y]=ode45('appollo',tspan,y0);
>> count=find(t>15);
>> y=y(count,:);
>> n=length(y(:,1));
>> for i=2:n-1
   if y(i-1,1)+eps<y(i,1) && y(i,1)>y(i+1,1)+eps
   plot(y(i,1),y(i,3),'.');
   hold on
   end
   end
>> xlabel('y1');ylabel('y3');
>> n=length(y(:,2));
>> for i=2:n-1
   if y(i-1,2)+eps<y(i,2) && y(i,2)>y(i+1,2)+eps
   plot(y(i,2),y(i,4),'.');
   hold on
   end
   end
>> xlabel('y2');ylabel('y4');
但是出的图不太正确,请大家帮忙修改一个
回复
分享到:

使用道具 举报

发表于 2012-5-30 00:09 | 显示全部楼层
建议去数学与非线性科学版块找找, 并看看这帖
建议提问的网友分清 编程问题 和 专业问题 http://forum.vibunion.com/thread-36746-1-1.html
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 00:33 , Processed in 0.054311 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表