chenvy 发表于 2009-2-14 19:35

求助：改进一下这个程序

这是我弄得一个程序，老师说我这个poincare图求解有问题，希望那位高手能指教，不甚感激
function dx=sbvp(t,x);
global a
global c
global A1
dx=;


functionpoincare
global a
global c
global A1;
c=0.1;A1=0.75;a=0.7;
%pi=3.14159265;
=ode45(@sbvp,,);
subplot(2,2,1),plot(x(160000:200000,1),x(160000:200000,2)),title('x(1) x(2)')
subplot(2,2,2),plot(t(190000:200000),x(190000:200000,1)),title('t x(1)')
subplot(2,2,3),plot(t(190000:200000),x(190000:200000,2)),title('t x(2)')
%m=zeros(2000,1);
%n=zeros(2000,1);
=ode45(@sbvp,,);
subplot(2,2,4),plot(y1(:,2),y1(:,1),'b.','markersize',3),title('poincare')

无水1324 发表于 2009-2-15 13:37

回复 楼主 chenvy 的帖子

对这个应该是有问题的，你求解的时候步长太大了。

T=2*pi;
dt=T/200;

=ode45(@sbvp,,);
subplot(2,2,4),plot(y1(40000:200:end,1),y1(40000:200:end,2),'b.','markersize',3),title('poincare')

chenvy 发表于 2009-2-15 21:43

谢谢，能告诉我为什么吗？麻烦了

jww100 发表于 2009-2-17 13:32

太~~~~~长了
页: [1]
查看完整版本: 求助：改进一下这个程序