Matlab 由轨道根数计算地心惯性坐标系下的状态向量：位置矢量r和速度矢量v

MAtlab 利用六个轨道根数计算飞行器在地心惯性坐标系的位置矢量r和速度矢量v：h比角动量 e偏心率 Omega升交点赤经 inc轨道倾角 omega近地点幅角 theta真近点角

function [rv] = coe2rv（heiomegaOmegathetavalid）
%coe2rv 由轨道根数计算地心惯性坐标系下的位置矢量r和速度矢量v
%输入：h比角动量  e偏心率  Omega升交点赤经  inc轨道倾角  omega近地点幅角 theta真近点角
% valid为显示过结果的有效数字
% 距离单位：千米，角度单位：度，时间单位：秒
miu = 398600;            % 地球引力常数 km^3/s^2
%近焦点坐标系下的位置矢量rx和速度矢量vx
rx = h^2/miu/（1 + e*cosd（theta）） * [cosd（theta）;sind（theta）;0];
vx = miu/h*[-sind（theta）;e + cosd（theta）;0];