%?function?[hs?As]?=?sunpoint（picName）
%%通過經緯度，當地觀測時間計算太陽高度角和方位角；
clear;
clc;
%?readDir?=?‘d:\我的文檔\桌面\1‘;
%?readPath?=?[readDir?‘\*.jpg‘];
%?readList?=?dir（readPath）;
%?picName?=?readList.name;
%?NF?=?str2num（picName（1:4））?
%?Y?=?str2num（picName（5:6））
%?R?=?str2num（picName（7:8））
%?S?=?str2num（picName（9:10））
%?F?=?str2num（picName（11:12））;
NF??=?2011;?Y?=3;R?=?25;S?=7;F=2;%%NFYRS?F?分別是觀測時的年月日時值和分值
w?=?30+35/60+6/3600;??%w為當地緯度（成都w?=?30+35/60+6/3600）
D?=103;E=59+23/60;%D、?E?分別是觀測處經度的度值和分度值（成都D?=103;E=59+23/60）
%%（1）求時差Eq
%%時差Eq?=?0.0028?-?1.9857sina+?9.9059sin2a-7.0924cosa?-?0.6882cos2a;
%%式中:?a為日角其求解公式為:a=?2*pi*?t/?365.2422;
%%t?由3部分組成?即?t?=?N?+N1?-?N?0;
N0?=?79.6764?+?0.2422*（?NF?-?1985）?-?floor（（?NF?-?1985）/4）;?%式中:?NF為年份?floor（?X）?為?MAT?LAB語言中求出不
%大于?X?的最大整數的標準函數。
%%N?為積日?所謂積日即日期在年內的順序號?例如?平年12月?31日的積日為?365閏年則為366。
%%當給定年月日分別為?NF??Y?R?時?相應的積日N?可通過以下MAT?LAB偽代碼來實現:
A?=?NF/4;%%NF為年份，Y為月，R為日
B=?A?-?floor（A）?;?
C=?32.8;
if?（Y<=2）
????C=?30.6;?
end
if?（?B==?0）?&?（Y>?2）?
????C=?31.8;