下面是一種隨經(jīng)緯度變化的日出日落時(shí)間計(jì)算方法，我成功運(yùn)用在一智能路燈控制器中，希望對(duì)需要的朋友有幫助。

已知：日出日落時(shí)太陽的位置h＝-0.833°，要計(jì)算地的地理位置，經(jīng)度Long，緯度G1at，時(shí)區(qū)zone，UTo為上次計(jì)算的日出日落時(shí)間，第一次計(jì)算時(shí)UTo＝180°。

(1)先計(jì)算出從格林威治時(shí)間公元2000年1月1日到計(jì)算日天數(shù)days;

(2)計(jì)算從格林威治時(shí)間公元2000年1月1日到計(jì)算日的世紀(jì)數(shù)t，

    則        t＝(days+UTo／360)／36525;

(3)計(jì)算太陽的平黃徑   L=280.460+36000.770×t;

(4)計(jì)算太陽的平近點(diǎn)角

    G＝357.528+35999.050×t

(5)計(jì)算太陽的黃道經(jīng)度

    λ＝L+1.915×sinG+0.020xsin(2G);

(6)計(jì)算地球的傾角  ε＝23.4393-0.0130×t;

(7)計(jì)算太陽的偏差  δ＝arcsin(sinε×sinλ);

(8)計(jì)算格林威治時(shí)間的太陽時(shí)間角GHA：

    GHA=UTo-180-1.915×sinG-0.020×sin(2G) +2.466×sin(2λ)-0.053×sin(4λ)

(9)計(jì)算修正值e：

    e=arcos{[ sinh-sin(Glat)sin(δ)]/cos(Glat)cos(δ)}

(10)計(jì)算新的日出日落時(shí)間

UT＝UTo-(GHA+Long±e);

其中“+”表示計(jì)算日出時(shí)間，“-”表示計(jì)算日落時(shí)間;

(11)比較UTo和UT之差的絕對(duì)值，如果大于0.1°即0.007小時(shí)，把UT作為新的日出日落時(shí)間值，重新從第(2)步開始進(jìn)行迭代計(jì)算，如果UTo和UT之差的絕對(duì)值小于0.007小時(shí)，則UT即為所求的格林威治日出日落時(shí)間;

(12)上面的計(jì)算以度為單位，即180°=12小時(shí)，因此需要轉(zhuǎn)化為以小時(shí)表示的時(shí)間，再加上所在的時(shí)區(qū)數(shù)Zone，即要計(jì)算地的日出日落時(shí)間為

T=UT/15+Zone

上面的計(jì)算日出日落時(shí)間方法適用于小于北緯60°和南緯60°之間的區(qū)域，如果計(jì)算位置為西半球時(shí)，經(jīng)度Long為負(fù)數(shù)。