我正在写一个脚本
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(