中华农历论坛知识讨论区历法知识 → [求助]就许建伟先生的《天文算法》一书中的疑点请教各位大侠


  共有11197人关注过本帖树形打印

主题:[求助]就许建伟先生的《天文算法》一书中的疑点请教各位大侠

帅哥哟,离线,有人找我吗?
看不懂算法
  1楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:新手上路 帖子:5 积分:216 威望:0 精华:0 注册:2014/11/2 16:10:00
[求助]就许建伟先生的《天文算法》一书中的疑点请教各位大侠  发帖心情 Post By:2014/11/2 17:01:00

    许建伟先生以及其同好们翻译的《天文算法》一书很好,解决了大部分业余爱好者的难题。

    看了这本书,就想自己学着编写程序(自学的编写程序,不是专业出身),因为许建伟先生提供的是javascript版本的源文件,而我看不懂js语言,于是就想自己编代码来实现农历的计算功能。

    在参考《天文算法》一书的过程中,问题来了。里面的知识很多,我【无法确定哪些章节】是【农历的天文计算】需要掌握的知识点。

    我只想学习【农历的每年开始日期计算,农历的每月的开始日期计算,每个节气的时间计算(精确到秒%)】

    若有大侠肯把那些章节指出来,不胜感激,小生在此三拜了。

    下面,我就贴出自己编写的部分章节的VBA代码(可以在Excel中运行)

  'VBA制作万年历的程序思路
'1、先算出儒略日JD,在【第七章】
‘A = INT(Y/100)  B = 2 - A + INT(A/4)
‘JD = INT(365.25(Y+4716))+INT(30.6001(M+1))+D+B-1524.5 (7.1)
Function rlr(nian, yue, ri, shi, fen, miao) '用年月日时分秒计算儒略日JD
a = Int(nian / 100)
b = 2 - a + Int(a / 4)
If yue = 1 Or yue = 2 Then
nian = nian - 1
yue = yue + 12
End If
fen = fen + miao / 60
shi = shi + fen / 60
ri = ri + shi / 24
rlr = Int(365.25 * (nian + 4716)) + Int(30.6001 * (yue + 1)) + ri + b - 1524.5
End Function



'''''''''''''''''''''''''以下不是必须的代码
'闰年的判断
Function rn(nian) '闰年平年的判断K
If nian Mod 400 = 0 Then
rn = 1
Else
If nian Mod 100 = 0 Then
rn = 2
Else
If nian Mod 4 = 0 Then
rn = 1
Else
rn = 2
End If
End If
End If
End Function
'星期几的判断
Function xqj(rlr) '星期几的判断
xqj = (rlr + 1.5) Mod 7
End Function
'年内序数日的判断
Function nnxsr(nian, yue, ri) '年内序数日的计算,一年内的第几日
k = rn(nian)
nnxsr = Int(275 * yue / 9) - k * Int((yue + 9) / 12) + ri - 30
End Function


支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
帅哥哟,离线,有人找我吗?
看不懂算法
  2楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:新手上路 帖子:5 积分:216 威望:0 精华:0 注册:2014/11/2 16:10:00
  发帖心情 Post By:2014/11/2 17:06:00

‘2、再求出力学时TD和世界时UT的差异ΔT,'UT = TD –ΔT
Function shicha1(rlr) '计算力学时和世界时的秒差ΔT,ΔT = -15 + (JD-2382148)2/41048480
a = rlr - 2382148
b = a * a
c = 41048480
shicha1 = -15 + b / c
End Function

Function shicha2(nian) '计算力学时和世界时的秒差的另一种方法ΔT
a = nian - 1810
b = a * a
shicha2 = -15 + 0.00325 * b
End Function

Function shicha3(rlr) 'ΔT = 102.3 + 123.5T+32.5T2
a = sjs(rlr)
a2 = a * a
shicha3 = 102.3 + 123.5 * a + 32.5 * a2
End Function

‘世界时=力学时减去时差
Function ut(td, shicha1) 'UT = TD –ΔT
ut = ta - shicha1
End Function

'计算儒略世纪数T
Function sjs(rlr) '计算儒略世纪数T
sjs = (rlr - 2451545#) / 36525
End Function


支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
帅哥哟,离线,有人找我吗?
看不懂算法
  3楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:新手上路 帖子:5 积分:216 威望:0 精华:0 注册:2014/11/2 16:10:00
  发帖心情 Post By:2014/11/2 17:06:00

'3、恒星时与格林尼治时间的关系 '本文所述的时间指世界时UT,世界时(0h UT)表示世界时0点,(3h UT)表示世界时3点。如果没有特另申明,均指格林尼治子午圈时间。 '格林尼治子午圈的恒星时可按如下计算: 'θo = 280.46061837 + 360.98564736629*(JD-2451545.0) + 0.000387933*T^2 - T^3/38710000 Function hxs(rlr) '计算任意时刻的格林尼治平恒星时θo sjs1 = sjs(rlr) sjs2 = sjs1 * sjs1 sjs3 = sjs2 * sjs1 hxs1 = 280.46061837 + 360.98564736629 * (rlr - 2451545#) + 0.000387933 * sjs2 - sjs3 / 38710000 hxs = jdzh(hxs1) '此时的数据是很大的度数,把±360之外的度数转换为360之内 hxs = hxs / 15 '每15度是一个小时,360/24=15 h = Int(hxs) If h < 10 Then h = "0" & h m = Int((hxs - h) * 60) If m < 10 Then m = "0" & m s = (hxs - h - (m / 60)) * 3600 If s < 1 Then s = "0" & s hxs = h & "h" & m & "m" & s End Function

支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
帅哥哟,离线,有人找我吗?
看不懂算法
  4楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:新手上路 帖子:5 积分:216 威望:0 精华:0 注册:2014/11/2 16:10:00
  发帖心情 Post By:2014/11/2 17:08:00


这格式怎么这么乱呢?谁教教我怎么编排格式?
复制进去都是排好的格式,发表了就乱了。

支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
帅哥哟,离线,有人找我吗?
浪-淘-沙
  5楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:版主 帖子:2068 积分:4263 威望:5 精华:4 注册:2008/11/13 21:03:00
  发帖心情 Post By:2014/11/3 21:56:00

以下是引用看不懂算法在2014-11-2 17:08:00的发言:
这格式怎么这么乱呢?谁教教我怎么编排格式? 复制进去都是排好的格式,发表了就乱了。

需要把回车符自己加进去。估计与论坛的版本有关。
像我这样,每换一个段落,插入一个回车符。
回车符就是/b br,都要用尖括号刮起来。

支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
帅哥哟,离线,有人找我吗?
浪-淘-沙
  6楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:版主 帖子:2068 积分:4263 威望:5 精华:4 注册:2008/11/13 21:03:00
  发帖心情 Post By:2014/11/3 21:58:00

如果你计算的结果与《寿星天文历》的数据差很多,你要查看一下程序的代码,是哪里有错漏?

支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
帅哥哟,离线,有人找我吗?
zhanglei1953
  7楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:论坛游民 帖子:34 积分:350 威望:0 精华:0 注册:2014/4/12 20:35:00
  发帖心情 Post By:2014/11/6 10:23:00

我已在论坛里发了EXCEL表计算节气时刻的程序,可以参考一下。
http://bbs.nongli.net/dispbbs_2_64668.html

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

返回版面帖子列表

[求助]就许建伟先生的《天文算法》一书中的疑点请教各位大侠








签名