日出日落時(shí)間計(jì)算文檔

1、日出日落時(shí)間計(jì)算：是指假設(shè)地球?yàn)榍蛐?，并允許大氣折射影響存在的條件下，圓盤形太陽的上沿剛好到達(dá)地平線的時(shí)刻，這相當(dāng)于以觀察者所處的地平面為參考平面太陽處在-0.833°的位置（計(jì)算時(shí)間按角度計(jì)算，即時(shí)間角，180°=12小時(shí)）。計(jì)算日出時(shí)，首先假設(shè)太陽的位置即日出時(shí)間為12h=180°，采用這個(gè)假設(shè)的位置計(jì)算日出位置，并用這個(gè)新位置采用迭代法重新計(jì)算太陽的位置，直到滿足要求的計(jì)算精度為止。輸入：該地的地理位置，如經(jīng)度，緯度，時(shí)區(qū)，當(dāng)日的日期，RiseSet（日出日落標(biāo)志）輸出：當(dāng)日的日出時(shí)間（RiseSet=1時(shí)）和日落時(shí)間（RiseSet=0時(shí)）。程序步驟：已知

2、：日出日落時(shí)太陽的位置h-0.833°,要計(jì)算地的地理位置,經(jīng)度Long,緯度G1at,時(shí)區(qū)zone,UTo為上次計(jì)算的日出日落時(shí)間,第一次計(jì)算時(shí)UTo180°。 (1)先計(jì)算出從格林威治時(shí)間公元2000年1月1日到計(jì)算日天數(shù)Days; (2)計(jì)算從格林威治時(shí)間公元2000年1月1日到計(jì)算日的世紀(jì)數(shù)t, 則        t(Days+UTo360)36525; (3)計(jì)算太陽的平黃徑   L=280.460+36000.770×t; (4)計(jì)算太陽的平近

3、點(diǎn)角 G357.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：   &#

4、160; e=arcos sinh-sin(Glat)sin()/cos(Glat)cos() (10)計(jì)算新的日出日落時(shí)間 UTUTo-(GHA+Long±e); 其中“+”表示計(jì)算日出時(shí)間,“-”表示計(jì)算日落時(shí)間; (11)比較UTo和UT之差的絕對值,如果大于0.1°即0.007小時(shí),把UT作為新的日出日落時(shí)間值,重新從第(2)步開始進(jìn)行迭代計(jì)算,如果UTo和UT之差的絕對值小于0.007小時(shí),則UT即為所求的格林威治日出日落時(shí)間; (12)上面的計(jì)算以度為單位,即180°=12小時(shí),因此需要轉(zhuǎn)化為以小時(shí)表示的時(shí)間,再加上所在的時(shí)區(qū)數(shù)Zo

5、ne,即要計(jì)算地的日出日落時(shí)間為 T=UT/15+Zone (13)需要把以小時(shí)表示的時(shí)間轉(zhuǎn)換為時(shí)間點(diǎn)表示的時(shí)間，即幾點(diǎn)幾分，把小時(shí)的小數(shù)部分乘以60再除以100，加上整數(shù)部分就行了。上面的計(jì)算日出日落時(shí)間方法適用于小于北緯60°和南緯60°之間的區(qū)域,如果計(jì)算位置為西半球時(shí),經(jīng)度Long為負(fù)數(shù)。驗(yàn)證：北京地區(qū)緯度：39.917 經(jīng)度：116.38 東八區(qū)日期：2010.8.24日出時(shí)間：05:337696（05:33:04）日落時(shí)間：18:597195（19:00:26）日期：2010.6.30日出時(shí)間：04:487101（04:47:20）日落時(shí)間：19:470659（

6、19:48:34）日期：2009.2.24日出時(shí)間：06:562161（06:53:22）日落時(shí)間：17:599357（18:01:41）日期：2009.4.10日出時(shí)間：05:456831（05:42:38）日落時(shí)間：18:469744（18:48:39）2009-04-10 星期五日出：05:43:472009-04-10 星期五日落：18:47:09括號里邊的數(shù)據(jù)為括號前邊的數(shù)據(jù)為計(jì)算得到的數(shù)據(jù)，只表示到分鐘。北京經(jīng)度：116.46緯度：39.92（以下是最近七天的日出日落時(shí)刻）2010-12-17 星期五日出：07:29:412010-12-17 星期五日落：16:50:372010-

7、12-18 星期六日出：07:30:182010-12-18 星期六日落：16:50:582010-12-19 星期日日出：07:30:542010-12-19 星期日日落：16:51:212010-12-20 星期一日出：07:31:282010-12-20 星期一日落：16:51:472010-12-21 星期二日出：07:32:002010-12-21 星期二日落：16:52:142010-12-22 星期三日出：07:32:302010-12-22 星期三日落：16:52:442010-12-23 星期四日出：07:32:582010-12-23 星期四日落：16:53:15平谷經(jīng)度：

8、117.10緯度：40.13（以下是最近七天的日出日落時(shí)刻）2010-12-17 星期五日出：07:27:452010-12-17 星期五日落：16:47:242010-12-18 星期六日出：07:28:232010-12-18 星期六日落：16:47:462010-12-19 星期日日出：07:28:592010-12-19 星期日日落：16:48:092010-12-20 星期一日出：07:29:332010-12-20 星期一日落：16:48:342010-12-21 星期二日出：07:30:052010-12-21 星期二日落：16:49:022010-12-22 星期三日出：07:

9、30:352010-12-22 星期三日落：16:49:312010-12-23 星期四日出：07:31:032010-12-23 星期四日落：16:50:03經(jīng)度：116.46緯度：39.92（以下是最近七天的日出日落時(shí)刻）2012-02-24 星期五日出：06:55:342012-02-24 星期五日落：17:59:212012-02-25 星期六日出：06:54:082012-02-25 星期六日落：18:00:292012-02-26 星期日日出：06:52:422012-02-26 星期日日落：18:01:362012-02-27 星期一日出：06:51:152012-02-27 星期一日落：18:02:432012-02-28 星期二日出：06:49