%?function?[hs?As]?=?sunpoint（picName）
%%通過(guò)經(jīng)緯度，當(dāng)?shù)赜^測(cè)時(shí)間計(jì)算太陽(yáng)高度角和方位角；
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?分別是觀測(cè)時(shí)的年月日時(shí)值和分值
w?=?30+35/60+6/3600;??%w為當(dāng)?shù)鼐暥龋ǔ啥紈?=?30+35/60+6/3600）
D?=103;E=59+23/60;%D、?E?分別是觀測(cè)處經(jīng)度的度值和分度值（成都D?=103;E=59+23/60）
%%（1）求時(shí)差Eq
%%時(shí)差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語(yǔ)言中求出不
%大于?X?的最大整數(shù)的標(biāo)準(zhǔn)函數(shù)。
%%N?為積日?所謂積日即日期在年內(nèi)的順序號(hào)?例如?平年12月?31日的積日為?365閏年則為366。
%%當(dāng)給定年月日分別為?NF??Y?R?時(shí)?相應(yīng)的積日N?可通過(guò)以下MAT?LAB偽代碼來(lái)實(shí)現(xiàn):
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;