飞码网-免费源码博客分享网站

点击这里给我发消息

从位置/速度矢量看轨道元素|-Matlab源码下载

文件名称-飞码网: 从位置/速度矢量看轨道元素
文件类型-飞码网: 压缩包
文件大小-飞码网: 4.2 KB
作 者-飞码网: 你的眼睛会笑弯成一条桥
下载点数-飞码网: 0
下载权限-飞码网: 普通会员
整理时间-飞码网: 2021-04-06 23:30:30
文件说明-飞码网:

vec2orbElem(rs,vs,mus)将
具有引力参数(mus)的物体的位置(rs)和速度(vs)转换为Keplerian轨道元素。

输入:
rs:3n x 1堆叠的初始位置向量:
[r1(1); r1(2); r1(3); r2(1); r2(2)r2(3); ...; rn(1) ; rn(2); rn(3)]
或位置向量的3 xn矩阵。
vs:3n x 1个堆叠的初始速度矢量或3 xn个矩阵
mus:引力参数(G * m_i)其中G是
引力常数,m_i是第i个物体的质量。
如果所有向量都代表同一物体,则mus可能是
标量。
输出:
a:半长轴
e:偏心距
E:偏心距异常
I:倾角
ω:
中心点论点Ω:上升节的经度
P:轨道周期
tau:
根尖周越时间A,B:方向矩阵(参见Vinti,1998) )

所有单位必须是互补的,即,如果位置以AU为单位,并且
时间以天为单位,则dx0必须以AU /天为单位,mus必须以
AU ^ 3 / day ^ 2(这些是solarSystemData.mat中的单位)。

solarSystemData.mat中的数据是从JPL的系统Web
界面(http://ssd.jpl.nasa.gov/?horizo​​ns)下载的。它包括
行星,太阳和冥王星的位置(因为我
在2006年之前上过小学)。带卫星的行星的位置是
重心。

示例:
%太阳系统Oribtal元素
ssdat = load('solarSystemData.mat');
rs = ssdat.p0(1:end-3)-repmat(ssdat.p0(end-2:end),9,1);
vs = ssdat.v0(1:end-3)-repmat(ssdat.v0(end-2:end),9,1);
mus = ssdat.mus(1:9)+ ssdat.mus(10);
[a,e,E,I,omega,Omega,P,tau,A,B] = vec2orbElem(rs,vs,mus);
%转换回来:
r = A * diag(cos(E)-e)+ B * diag(sin(E));
rdot =(-A * diag(sin(E))+ B * diag(cos(E)))* ...
diag(sqrt(mus(:)。'。* a。^-3)./(1 -e。* cos(E)));

下载地址-飞码网: [ 下载地址1 ]
提 取 码-飞码网:
下载说明-飞码网: ⊙本网站所有资料都是由飞码网或其它作者上传,若资料造成侵权,请联系我们删除。
⊙如果文件出现下载问题,或链接失效,请在评论区反馈或者直接联系我们!
⊙如果下载网页缓冲慢,可能是网页浏览量比较大,请稍后再试!
⊙如果遇到积分扣除错误或其它问题,请直接点击qq图标联系我们!
⊙本站所提供的资料均由飞码网在网络上收集作为学习用途,请勿商用!
赞 ()

相关推荐

内容页底部广告位3
留言与评论(共有 0 条评论)
   
验证码: