http://www.pveducation.org/pvcdrom/properties-of-sunlight/suns-position
北京,纬度39,经度-116.
计算 Time Correction Factor (TC)的时候,
经度和lstm 取的都是度数, 经度(东经为负,西为正)
最后公式中的alpha为 elevation.
fai为 维度(北纬为正,南为负)
以北京2015/3/24为例,d=83
B=360/365.0*(83-81)/180*pi
lstm=15*(-8)
eot=9.87*sin(2*B)-7.53*cos(2*B)-1.5*sin(B)
tc=4*(-116-lstm)+eot
lst=9.75(当前时间 09:45分)+tc
hra=15/180.0*pi*(lst-12)
posai=23.45/180.0*pi*sin(360/365.0*(83-81))
elev=asin(sin(39.0/180.0*pi)*sin(posai)+cos(39.0/180.0*pi)*cos(posai)*cos(hra))
alt=39.0/180.0*pi
azi=acos((sin(uk)*cos(alt)-cos(uk)*sin(alt)*cos(hra))/cos(elev))
degree=azi/pi*180. 结果为113.71 (从正北,顺时针方向度数)