资源描述:
《日出日落时间计算文档》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、日出日落时间计算:是指假设地球为球形,并允许大气折射影响存在的条件下,圆盘形太阳的上沿刚好到达地平线的时刻,这相当于以观察者所处的地平面为参考平面太阳处在-0.833°的位置(计算时间按角度计算,即时间角,180°=12小时)。计算日出时,首先假设太阳的位置即日出时间为12h=180°,采用这个假设的位置计算日出位置,并用这个新位置采用迭代法重新计算太阳的位置,直到满足要求的计算精度为止。输入:该地的地理位置,如经度,纬度,时区,当日的日期,RiseSet(日出日落标志)输出:当日的日出时间(Rise
2、Set=1时)和日落时间(RiseSet=0时)。程序步骤:已知:日出日落时太阳的位置h=-0.833°,要计算地的地理位置,经度Long,纬度G1at,时区zone,UTo为上次计算的日出日落时间,第一次计算时UTo=180°。(1)先计算出从格林威治时间公元2000年1月1日到计算日天数Days;(2)计算从格林威治时间公元2000年1月1日到计算日的世纪数t,则 t=(Days+UTo/360)/36525;(3)计算太阳的平黄径 L=280.460+36000.770×t;(
3、4)计算太阳的平近点角G=357.528+35999.050×t(5)计算太阳的黄道经度λ=L+1.915×sinG+0.020xsin(2G);(6)计算地球的倾角 ε=23.4393-0.0130×t;(7)计算太阳的偏差 δ=arcsin(sinε×sinλ);(8)计算格林威治时间的太阳时间角GHA:GHA=UTo-180-1.915×sinG-0.020×sin(2G) +2.466×sin(2λ)-0.053×sin(4λ)(9)计算修正值e: e=arcos{[ sinh-si
4、n(Glat)sin(δ)]/cos(Glat)cos(δ)}(10)计算新的日出日落时间UT=UTo-(GHA+Long±e);其中“+”表示计算日出时间,“-”表示计算日落时间;(11)比较UTo和UT之差的绝对值,如果大于0.1°即0.007小时,把UT作为新的日出日落时间值,重新从第(2)步开始进行迭代计算,如果UTo和UT之差的绝对值小于0.007小时,则UT即为所求的格林威治日出日落时间;(12)上面的计算以度为单位,即180°=12小时,因此需要转化为以小时表示的时间,再加上所在的时区数
5、Zone,即要计算地的日出日落时间为T=UT/15+Zone(13)需要把以小时表示的时间转换为时间点表示的时间,即几点几分,把小时的小数部分乘以60再除以100,加上整数部分就行了。上面的计算日出日落时间方法适用于小于北纬60°和南纬60°之间的区域,如果计算位置为西半球时,经度Long为负数。验证:北京地区纬度:39.917经度:116.38东八区日期:2010.8.24日出时间:05:337696(05:33:04)日落时间:18:597195(19:00:26)日期:2010.6.30日出时间
6、:04:487101(04:47:20)日落时间:19:470659(19:48:34)日期:2009.2.24日出时间:06:562161(06:53:22)日落时间:17:599357(18:01:41)日期:2009.4.10日出时间:05:456831(05:42:38)日落时间:18:469744(18:48:39)2009-04-10[星期五] 日出:05:43:472009-04-10[星期五] 日落:18:47:09括号里边的数据为http://www.xdanger.com/vend
7、ors/sun-time/查得。括号前边的数据为计算得到的数据,只表示到分钟。北京经度:116.46 纬度:39.92(以下是最近七天的日出日落时刻)2010-12-17[星期五] 日出:07:29:412010-12-17[星期五] 日落:16:50:372010-12-18[星期六] 日出:07:30:182010-12-18[星期六] 日落:16:50:582010-12-19[星期日] 日出:07:30:542010-12-19[星期日] 日落:16:51:212010-12-20[星期一]
8、 日出: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平谷经度:117.10 纬度:40.13(以下是最近七