tle 根据卫星的位置和速度

根据卫星的位置和速度，求取实时卫星轨道的两行根数，即6个开普勒根数。计算结果精度很高。

clc;  close all;  clear all;
u=39860047e7;
e0=1e-6;
%r=[3047351.25796702-155873.158056562-6741613.595401350];
%v=[-6214.447953-2776.927754-2744.855816];
r=[-1600590.778830-5398039.858220-4340077.951502];

%r=r+30
v=[-233.7507886134455093.214922408990-5681.570113212520];
h=cross（rv）;
a=1/（2/norm（r）-sum（v.*v）/u）;
i=atan2（sqrt（h（1）^2+h（2）^2）h（3））;
e=sqrt（（sum（r.*v）/sqrt（u*a））^2+（norm（r）*sum（v.*v）/u-1）^2）;
if i    i=0;
    omiga=0;
    if e        M=atan2（r（2）r（1））;
        e=0;
        w=0;
    end
    if e>=e0
        w=atan2（r（2）r（1））-atan2（sqrt（1-e^2）*sum（r.*v）/sqrt（u*a）（（norm（r）*sum（v.*v）/u-1）-e^2））

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1300  2012-03-12 14:08  tle.m

----------- ---------  ---------- -----  ----

                 1300                    1