MATLAB答案

我的代码卡在“忙”上；我第一次使用ode45和函数。

5个视图（过去30天）

阿曼达·威利斯 2021年2月3日

评论：星际漫游者 2021年2月3日

接受的答复：星际漫游者

我正在尝试解一个方程组，这是我第一次使用函数或ode45。我无法成功运行我的代码而不让它一直处于“忙碌”状态。我试着停下来，一步一步地看台词，但我猜它是被困在颂歌45里了。我非常感谢任何关于如何更好地利用这一点的见解或建议。

                    %变量
                   
                    G=6.67*10^-11；%公里^3/（千克*s^2）
                   
                    M=8；%千克
                   
                    T=.0121；%N
                   
                    mdot=8.0*10^-7；%千克/秒
                   
                    %初始条件
                   
                    r0=19178；rdot0=0；
                   
                    θ0=0；thetadot0=sqrt（（G*M）/r0^3）；
                   
                    IC=[r0，rdot0，θ0，θ0]；
                   
                    %调用ode45解算器
                   
                    [t，state_values]=ode45（@f[02600000]，IC）；
                   
                    r=状态_值（：，1）；
                   
                    图（1）
                   
                    绘图（t，r）
                   
                    xlabel(“时间”);
                   
                    ylabel (“高度”);
                   
                    作用sdot=f（t，状态值）
                   
                    %状态变量
                   
                    r=状态_值（1）；%存储前一个时间步骤的半径值
                   
                    rdot=状态_值（2）；
                   
                    θ= state_values (3);
                   
                    thetadot=状态_值（4）；
                   
                    %参数
                   
                    G=6.67*10^-11；%公里^3/（千克*s^2）
                   
                    M=8；%千克
                   
                    T=.0121；%N
                   
                    mdot=8.0*10^-7；%千克
                   
                    d2rdt2=r*thetadot^2-（G*M）/r^2；
                   
                    d2thetad2=（T/（M/（mdot*T））-2*rdot*thetadot；
                   
                    sdot=[rdot d2rdt2 thetadot d2thetad2]；%存储为列向量
                   
                    终止