M函数文件：apolloeq.m
funtion xdot=apolloeq(t,x,flag,w);
xdot=[x(2);
sin(w*t)-132.509*x(1)+263.018*x(3)-0.0651*x(2)+1.198799*e(-6)*x(4);
x(4);
3.93999*e4*x(1)-3.94055*e4*x(3)+0.035915*(x(2)-x(4))];


>> w=0:0.01:40;
T=2*pi;
hold on
for j=1:length(w)
[t,x]=ode45('apolloeq',[0:T./100:70*T],(0.1;0;3e-6;0),[],w(j));
plot(d(j),x(500:100:1400,2),'linewidth',5);
title('分岔图二')
xlabel('参数 F0')
ylabel('位移')
end
大神帮忙看看,这个程序运行老是出错，“[t,x]=ode45('apolloeq',[0:T./100:70*T],(0.1;0;3e-6;0),[],w(j));
                                            |
错误: 圆括号或方括号不对称或异常。”个人觉得是x0的初值条件有问题 返回小木虫查看更多