MATLAB作稳态之后的庞加莱映射
微分方程为：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=;
>> tspan=;
>> =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');
但是出的图不太正确，请大家帮忙修改一个 建议去数学与非线性科学版块找找, 并看看这帖
建议提问的网友分清 编程问题 和 专业问题 http://forum.vibunion.com/thread-36746-1-1.html
页:
[1]