sun_azimuth_data.m

根据经度纬度时间

%SUN_AZIMUTH_DATA Generate an ASCII File (Report) containing Sun 
%                 Azimuth Data.
%
%   Using this script, you can generate a set of data regarding Sun 
%   Azimuth Positions at various instants of time for any location 
%   on Earth. This data can be used in many applications, like, in 
%   determination of orientation of any plane's azimuth from True 
%   North.
%
%   The syntax is 'sun_azimuth_data' (no arguments). It prompts for 
%   Observer's geodetic coordinates, observer's time zone and starting 
%   time [according to local standard time]. Then, it asks the Time 
%   Period and intervals for which Sun Azimuth data has to be generated.
%   It creates a MS Word file containing Sun Azimuth data.
%
%   Acknowledgements for the script 'sun_position.m' by Vincent Roy, 
%   which returns the Sun Position data for a given observer location 
%   and time.
%   
%   History
%    25/05/2004 - Original Creation by Khalil Sultan ([e-mail])
%    15/06/2004 - Code Modified.
%    08/07/2004 - Code Re-modified and Uploaded on MATLAB Central - File Exchange.

disp(' ')
disp('WARNING:')
disp('Skipping any field, with a hyphen before it, will result ')
disp('in termination of execution of this script. Fields with ')
disp('some value, in square brackets after them, represent their')
disp('default values.')
disp(' ')

% Input the Location of the Observer
locstr = input('- Enter the Location of the Observer: ','s');
if  isempty(locstr) == 1
    return
end

% Input the Country of the Observer
countrystr = input('  Enter the Country of the Observer: ','s');
if  isempty(countrystr) == 1
    country_set = 0;
else
    country_set = 1;
end

% Input the Name of the GPS Receiver Set used for determining
% the Observer's Geodetic Coordinates
gpsstr = input('  Enter the GPS Receiver Set Title: ','s');
if  isempty(gpsstr) == 1
    gps_set = 0;
else
    gps_set = 1;
end

% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% The Geodetic Coordinates of the Observer.
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
disp(' ')

% Input the Longitude in dd.dddddd