如题所述
出现这样的问题,最关键的是把ode45()的位置设置错误,不能放在其调用的函数内。这样的结果就到底是谁调用谁,分不清自己。
ode45()正确的方法是:
y0=[0;0]; %从orbit()中移出来
[t,x]=ode45(@orbit,[0 12],y0) %从orbit()中移出来
其他,orbit()本身还有问题,写的比较乱,orbit()要求反馈变量dy含糊不清。
还有什么疑问,可以再交流或通过私信的方式交流也可以。交流时,最后提供完整的代码,这样好帮你分析问题所在。