matlab中ode45如何设置,如何使用Matlab中的ode45修正赋值错误(ode45函数的第488行)

我正在写一个脚本

ode45

为了整合

卫星在火星附近双曲线轨道上的运动方程

.

我需要整合地球上的整个通道:从SOI半径开始(

576000km

)向行星前进,然后穿过大气层直到卫星到达

opposite

“大气边界(设置在

250km

从表面)。

当它接收到输入A时

tspan

比大约高

200000

秒(我需要大约

400000

秒),Matlab给出了以下信息:

'无法执行分配,因为左侧的大小为

4乘2,右边的尺寸是4乘5。

error happens in line 488 of Ode45

.

我搜索了一些类似的案例,但找不到任何东西,也不知道如何使用

条件断点

在一个复杂的函数上计算出一些东西,如ode45。我也尝试了不同的选择

odeset

“但没什么变化。我不知道错误可能在哪里。

这是脚本,其中我使用两个附加函数来获取一些参数:

Vinf=[2.7 4 6 8];

mu_m=42828;

R_msoi=.576e6;

[afas,efas,pfas]=parasint(Vinf);

an_vera0=zeros(1,length(Vinf));

Vr0=zeros(1,length(Vinf));

Vt0=zeros(1,length(Vinf));

for j=1:length(Vinf)

c_tstar0=(1/efas(j))*((pfas(