matlab航天工程工具箱相信对航天卫星仿真的人能有用.zip

matlab航天工程工具箱,相信对航天,卫星仿真的人能有用。。。

function oe=cart2orb（rrvvmu）
%CART2ORB  Cartesian coordinates to orbital elements.
%   OE = ORB2CART（RVMU） where R is the position vector （1x3）
%   and V is the velocity vector （1x3） at a given time point （epoch）.
%   MU is G*M for the central object
%   and OE is the orbital elements which has the form:
%
%   OE = [a e i W w M]
%
%   However only the first five elements will be used by this function.
%
%   Where:
%      a [m]   : semi-major axis
%      e []    : eccentricity
%      i [deg] : inclination
%      W [deg] : longitude of the ascending node
%      w [deg] : argument of periapsis
%      M [deg] : mean anomaly at epoch
%
%   M is the mean anomaly at epoch and is here the same as the mean anomaly
%   for the position R and velocity V given above. If R is the location of the
%   periapsis then the mean anomaly will be M = 0.
%
%   The xy-plane is the reference plane and the zero point of longitude
%   （or direction of vernal equinox） is in the direction of the x-axis.
%
%   See also ORB2CART TLE2ORB MA2TP MSATTRACK.

% Copyright （c） 2003-04-18 B. Rasmus Anthin.

r2d=180/pi;
rr=rr（1:）;
vv=vv（1:）;
hh=cross（rrvv）;
nn=cross（[0 0 1]hh）;
r=norm（rr）;
v=norm（vv）;
h=norm（hh）;
n=norm（nn）;
ee=（（v^2-mu/r）*rr-（rr*vv‘）*vv）/mu;
e=norm（ee）;

W=acos（nn（1）/n）;
if nn（2）<0 W=2*pi-W;end
i=atan（n/hh（3））;
a=mu*r/（2*mu-r*v^2）;
E=atan2（r*v/（e*sqrt（mu*a））（a-r）/（a*e））;
M=E-e*sin（E）;
w=acos（nn*ee‘/n/e）;
if ee（3）<0 w=2*pi-w;end

oe=[a e i*r2d W*r2d w*r2d M*r2d];

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2003-09-02 23:00  space\
     目录           0  2003-09-02 23:00  space\atmos\
     文件       69319  2003-09-02 23:00  space\atmos\atmdata.m
     文件       57680  2003-09-02 23:00  space\atmos\atmdata.mat
     文件       68440  2003-09-02 23:00  space\atmos\atmdata.txt
     文件        2832  2003-09-02 23:00  space\atmos\atmosphere.m
     文件        4872  2003-09-02 23:00  space\atmos\atmosphere.mat
     文件        3351  2003-09-02 23:00  space\atmos\stdatm.m
     文件        4913  2003-09-02 23:00  space\atmos\stdatm.txt
     文件        1516  2003-09-02 23:00  space\cart2orb.m
     文件          54  2003-09-02 23:00  space\cb_cancel.m
     文件         156  2003-09-02 23:00  space\cb_param_eng.m
     文件         244  2003-09-02 23:00  space\cb_param_ok.m
     文件         800  2003-09-02 23:00  space\checkinstall.m
     文件         596  2003-09-02 23:00  space\colors.cfg
     文件        2266  2003-09-02 23:00  space\contents.m
     文件         895  2003-09-02 23:00  space\daylight.m
     文件         387  2003-09-02 23:00  space\evalccm.m
     文件        3048  2003-09-02 23:00  space\evalcg.m
     文件        4734  2003-09-02 23:00  space\evalcp.m
     文件         372  2003-09-02 23:00  space\evaldcm.m
     文件         648  2003-09-02 23:00  space\greg2jul.m
     文件         365  2003-09-02 23:00  space\hohmann.m
     文件         392  2003-09-02 23:00  space\info.ins
     文件        1853  2003-09-02 23:00  space\install.m
     文件         264  2003-09-02 23:00  space\julepoch.m
     文件        1796  2003-09-02 23:00  space\kepler1.m
     文件         587  2003-09-02 23:00  space\loadcad.m
     文件        2048  2003-09-02 23:00  space\loadeng.m
     文件         974  2003-09-02 23:00  space\ma2tp.m
     目录           0  2003-09-02 23:00  space\maps\
............此处省略317个文件信息