求助：改进一下这个程序

chenvy

这是我弄得一个程序，老师说我这个poincare图求解有问题，希望那位高手能指教，不甚感激
function dx=sbvp(t,x);
global a
global c
global A1
dx=[x(1)-1/3*x(1)^3-x(2)+A1*cos(t);c*(x(1)+a-0.8*x(2))];


function  poincare
global a
global c
global A1;
c=0.1;A1=0.75;a=0.7;
%pi=3.14159265;
[t,x]=ode45(@sbvp,[0,200000],[0.01,0.1]);
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);
[t,y1]=ode45(@sbvp,[0:2*pi:200000],[0.01,0.1]);
subplot(2,2,4),plot(y1(:,2),y1(:,1),'b.','markersize',3),title('poincare')

无水1324

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

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

[t,y1]=ode45(@sbvp,[0:dt:2100*T],[0.01,0.1]);
subplot(2,2,4),plot(y1(40000:200:end,1),y1(40000:200:end,2),'b.','markersize',3),title('poincare')

chenvy

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

jww100

[0,200000]太~~~~~长了