[求助]网上找到的现成算法，求天亮时间的计算代码 [中华农历论坛]

中华农历论坛 → 知识讨论区 → 历法知识 → [求助]网上找到的现成算法，求天亮时间的计算代码

共有14438人关注过本帖树形打印主题：[求助]网上找到的现成算法，求天亮时间的计算代码

esnake

小大 1楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

等级：论坛游民帖子：40 积分：531 威望：0 精华：0 注册：2005/10/27 21:24:00

[求助]网上找到的现成算法，求天亮时间的计算代码 发帖心情

Post By：2013/3/18 10:48:00

这段代码是从网上找到的，里面只有日出的时间计算，要是想计算天亮的时间，该怎样写函数呢？

好象是许老师的代码，求教许老师，能否赐教？先谢过！

下载信息 [文件大小：下载次数： ]
点击浏览该文件:code.txt

支持(0) 中立(0) 反对(0) 单帖管理 | 引用 | 回复

浪-淘-沙

小大 2楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

等级：版主帖子：2068 积分：4263 威望：5 精华：4 注册：2008/11/13 21:03:00

Post By：2013/3/18 11:17:00

通常是假定天亮时刻为太阳的中心线在地平线之下6度。

支持(0) 中立(0) 反对(0) 单帖管理 | 引用 | 回复

esnake

小大 3楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

等级：论坛游民帖子：40 积分：531 威望：0 精华：0 注册：2005/10/27 21:24:00

Post By：2013/3/18 11:39:00

function sheng(jd,L,fa,TZ){ //jd儒略日平午，L地理经度，fa地理纬度，TZ时区，返回太阳升起时间 jd-=TZ; var T=jd/36525; var J = sunLat(T), sinJ=sin(J), cosJ=cos(J); //太阳黄经以及它的正余弦值 var gst = 2*Math.PI*(0.7790572732640 + 1.00273781191135448*jd) //恒星时（子午圈位置） + (0.014506 + 4612.15739966*T + 1.39667721*T*T)/rad; var E = (84381.4060 -46.836769*T)/rad; //黄赤交角 var A=Math.atan2( sinJ*cos(E), cosJ); //太阳赤经 var D=Math.asin ( sin(E)*sinJ ); //太阳赤纬 var cosH0 = (sin(-50*60/rad) - sin(fa)*sin(D) ) / ( cos(fa)*cos(D) ); //日出的时角计算 if(cosH0>=1||cosH0<=-1) return 0; var H0 = -Math.acos( cosH0 ); //升点时角（日出） var H = gst-L-A; //太阳时角 return jd - rad2rrad(H-H0)/6.28+TZ; } 这是计算日出的代码