计算卫星角度程序，包括高度角、方位角

function angles=satang4（strSys1latAlongAlongsatstrSys2latBlongB）
%function angles=satang4（strSys1latAlongAlongsatstrSys2latBlongB）
%calculation of satellite angles: elevation azimuth and also relative azimuth and baseline orientation if input data
%belongs to two sites
%Sources:
%Digital satellite communications 2nd Ed. Author Tri Ha for satellite azimuth and elevation calculation
%http://williams.best.vwh.net/avform.htm#Rhumb for distance and relative azimuth between two coordinates
%tested using various satellites and earth stations and using http://www.sfwx.com/softworks/geoBearing.html and
%http://www.geocities.com/senol_gulgonul/linkbudget/azel.html
%differences are in part due to the selection of the earth radius and orbit heigth values and also to possible mistakes in the
%code...please point them out
%INPUTS
%strSys1: string: ‘E‘ for longitude in degrees east ‘W‘ for longitude in degrees West. no negative values
%strSys2: string: ‘GC‘ for GREAT CIRCLE     ‘RH‘ for RHUMB （CONSTANT AZIMUTH: recommended for site diversity purposes）  
%latXX: latitude of site A or B. + north - south
%longXX: longitude of site A and B: following strSys1.
%longsat: longitude of subsatellite point E or W according to srtSys
%OUTPUTS
%angles: a structure containing:
%       .elevation              elevation angle to satelite
%       .azimuth                azimuth to satelite
%       .poltilt & poltilt_c    feed rotation to compensate geographic position
                                %NOTE standing behind the antenna
                                %rotation is CLOCKWISE for a station
                                %located east of the satellite meridian
                                %and ANTICLOCKWISE for a station west of
                                %the satellite meridian
                                %poltilt_c is more reliable than poltilt.
%       .Az_AB                  true course azimuth in the direction of site B （rhumb line constant azimuth）
%       .distance               distance between the two sites
%       .baselineorientation    for site diversity BOA acording to ITU-R P618-8 definition.
%EXAMPLE OF USE
%angles=satang4（‘E‘-8.8313.250）          |   angles=satang（‘E‘-33.9318.470）
%angles =                                   |
%    elevation: 71.2774                     |        elevation: 45.7078
%     azimuth: 303.1005                     |        azimuth: 329.1038
%Meteosat satellite cities of cape town and luanda.
%
%angles=satang4（‘E‘6.4360-76.4304.5）
%angles = 
%        elevation: 64.3472
%        azimuth: 106.2730      Intelsat 8A city of medellin
%geographical constants

% earthradius=6378.155;   %according to TRI Ha book page 41
earthradius=6371;   %http://www.movable-type.co.uk/scripts/latlong.html
satheigth=41264.2;

 I=find（longA<0）;
 if isempty（I）==0
     longA（I）=360+longA（I）;
 end
 clear I
[Az El poltilt poltilt_c]=pointangle（strSys1