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];