太阳方位角

　　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 （从正北，顺时针方向度数）

原文：http://www.cnblogs.com/doublemystery/p/4362097.html