GPS天空视图

matlab编写的一个GPS天空视图，可以根据输出的年月日UTC时和历书文件显示GPS天空视图

function [ElevationAzimuth]=Angle_calculation（enu）

X_unit=[100];
Y_unit=[010];
Z_unit=[001];

enu_unit=enu/norm（enu）;
phi=acos（dot（enu_unitZ_unit））*180/pi;

Elevation=90-phi;
a=[enu_unit（1）enu_unit（2）0];
a=a/norm（a）;
theta=acos（dot（aX_unit））*180/pi;
% if（enu_unit（1）>=0&&enu_unit（2）>=0）
%    Azimuth=90-theta; 
% elseif（enu_unit（1）<0）
%     Azimuth=（90-theta）+360;
% elseif（enu_unit（2）<0）   
%     Azimuth=90+theta;
% end   
if（enu_unit（1）>=0&&enu_unit（2）>=0）
   Azimuth=90-theta; 
elseif （enu_unit（1）>=0&&enu_unit（2）<0）
    Azimuth=90+theta;
    elseif （enu_unit（1）<0&&enu_unit（2）<0）
    Azimuth=90+theta;
        elseif （enu_unit（1）<0&&enu_unit（2）>=0）
    Azimuth=450-theta;
end
    

end

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

     文件      10374  2013-08-13 21:12  历书解算及天顶视图\2013_08_13_1753_405504.alm

     文件      10370  2013-08-02 20:31  历书解算及天顶视图\Almanac.alm

     文件      10374  2013-08-02 20:31  历书解算及天顶视图\Almanac_backup.alm

     文件        733  2013-08-02 20:31  历书解算及天顶视图\Angle_calculation.m

     文件       2214  2013-08-02 20:31  历书解算及天顶视图\CompSatPos.m

     文件        234  2013-08-02 20:31  历书解算及天顶视图\CompSatPosFromAlm.m

     文件        760  2013-08-03 16:33  历书解算及天顶视图\inputParameter.m

     文件      10370  2013-08-02 20:31  历书解算及天顶视图\new\Almanac.alm

     文件       1206  2013-08-02 20:31  历书解算及天顶视图\new\llh2xyz.m

     文件        451  2013-08-02 20:31  历书解算及天顶视图\new\ReadAlm.m

     文件       2661  2013-08-02 20:31  历书解算及天顶视图\new\skyplot.m

     文件       3043  2013-08-02 20:31  历书解算及天顶视图\new\svposalm.m

     文件       1609  2013-08-02 20:31  历书解算及天顶视图\new\test.m

     文件       1565  2013-08-02 20:31  历书解算及天顶视图\new\xyz2enu.m

     文件       1615  2013-08-02 20:31  历书解算及天顶视图\new\xyz2llh.m

     文件       5828  2013-08-03 13:26  历书解算及天顶视图\ReadAlm.m

     文件       1639  2015-06-02 10:58  历书解算及天顶视图\ReadAlm_GPS.m

     文件      18731  2013-10-22 09:39  历书解算及天顶视图\read_BD_eph.m

     文件       1581  2013-10-21 22:23  历书解算及天顶视图\Satellite_Location_Calculate.m

     文件       5625  2013-08-03 16:35  历书解算及天顶视图\skyplot.m

     文件       1920  2013-08-02 20:31  历书解算及天顶视图\skyplotun.m

     文件       2753  2015-06-02 11:11  历书解算及天顶视图\test.m

     文件       1658  2013-08-03 15:05  历书解算及天顶视图\utc2gps.m

     文件       1582  2013-08-02 20:31  历书解算及天顶视图\wgslla2xyz.m

     文件       1690  2013-08-02 20:31  历书解算及天顶视图\wgsxyz2enu.m

     文件       1195  2013-08-02 20:31  历书解算及天顶视图\wgs_84_parameters.m

     文件      27127  2013-08-02 20:31  历书解算及天顶视图\历书解算及天顶视图.docx

     目录          0  2013-08-04 09:52  历书解算及天顶视图\new

     目录          0  2015-06-02 11:37  历书解算及天顶视图

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

............此处省略2个文件信息