中华农历论坛知识讨论区历法知识 → 我自己制作的万年历


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

主题:我自己制作的万年历

帅哥哟,离线,有人找我吗?
春光
  51楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:版主 帖子:898 积分:7341 威望:10 精华:1 注册:2005/3/17 11:37:00
  发帖心情 Post By:2008/4/4 9:00:00

再次谢谢xjw01兄真心无私的劳动与奉献。

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


加好友 发短信
等级:新手上路 帖子:14 积分:575 威望:0 精华:0 注册:2007/3/28 10:57:00
  发帖心情 Post By:2008/4/5 18:21:00

楼上的xjw01先生真厉害!希望你能用你的算法做一个公能全面的万年历程序,而不仅仅只是计算节气和朔的时刻。最好还能公开源程序。

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


加好友 发短信
等级:蝙蝠侠 帖子:866 积分:3824 威望:3 精华:20 注册:2008/3/20 22:14:00
  发帖心情 Post By:2008/4/5 22:36:00

第20章 岁差
  [许剑伟 清明节下午 在家翻译 翻译工具google.cn语言]
  地球的自转轴的方向并不是真的固定在空间的某一方向上。随着时间的推移,它缓慢的漂移,就像陀螺一样,这种漂移称为岁差。这个效果源于月球和太阳对地球赤道突起部分的引力。
  由于岁差,北极(目前位于αUrsae Minoris星或Polaris星附近)缓慢的围绕黄极转动,周期是26000年。因这一原因,春风点(黄道与赤道的交点)每年退行50"。
  另外,行星的黄道在空间中也不是固定的。由于行星对地球的摄动,它回绕一个交点缓慢的移动,每世纪47"。
  黄道面、赤道面及春风点,是黄道坐标及赤道坐标系统中的基本平面和起点。黄道坐标及赤道坐标也是天球中的两个重要的坐标系统:黄道坐标(黄经λ及黄纬β)及赤道坐标(赤经α及赤纬δ)。由于岁差,固定的恒星的坐标却是不断改变的。因此,星库中,赤经赤纬表中总要给定某一个历元时间,如1900.0或1950.0或2000.0。
  在这一章中,我们还考虑赤经α及赤纬δ转换问题,即从一个历元坐标中转到另一个历元坐标中。这里只考虑平位置及岁差,而星体的视位置将在第22章中考虑。
  如果精度要求不高,两个历元相差不远,并且如果星体没有太靠近天极,下面的公式可以用来计算赤经及赤纬的每年度岁差:
    Δα = m + n*sin(α)*tan(δ)      Δδ = n*cos(α)   ……20.1式
  式中m和n是两个随时间不同速度缓慢变化的量。他们表达为:
    m = 3s.07496 + 0s.00186*T
    n = 1s.33621 - 0s.00057*T
    n = 20".0431 - 0".0085*T
  T是J2000.0起算的儒略世纪数。以下是某些历元下的m和n的值:
    历元    m(单位s)  n(单位s)  n(角秒)
   1700.0    3.069     1.338     20.07
   1800.0    3.071     1.337     20.06
   1900.0    3.073     1.337     20.05
   2000.0    3.075     1.336     20.04
   2100.0    3.077     1.336     20.03
   2200.0    3.079     1.335     20.03
  计算Δα时,n值的单位是s。记住:1s=15"。
  如果是恒星的情况,计算星体位置时,应把恒星自行效果加到20.1式。

  例20.a ——Reulus(αLeonis)在历元J2000.0时的坐标是:
      αo = 10h 08m 22s.3    δo = +11°58'02"
    每年自行速度是:
      -0s.0169    (赤经)
      +0".006     (赤纬)
    请把该坐标转为历元1978.0的分点坐标。
  解:依题意有:
    αo = 152°.093    m = 3s.075
    δo = +11°.967    n = 1s.336 = 20".04
  用20.1式可算出:
    Δα = +3s.208,   Δδ = -17".71
  我们还应加上每年的自行速度,得到每所变化速度为:
    Δα = +3s.191,   Δδ = -17".70
  在2000.0到1978.0这-22年里,总变化量是:
   α的: +3s.191 * (-22) = -70s.2 = -1m 10s.2
   δ的: -17".70 * (-22) = +389"  = +6'29"
  最的得到:
    α = αo - 1m10s.2 = +10h 07m 12s.1
    δ = δo + 6'29"   = +12°.04'31"
-----------------------
贝塞尔年及儒略年(两种常用历元)

  国际天文学联合会决定,从1984年起,天文历表应使用以下系统。
  新的标准历元是2000年1月1日12时(TD时),对应JDE 2451545.0。该历元写作J2000.0。为了星体位置计算方便,实历每年的年首与标准历元J2000.0中的各年年首是不同的,J2000中每年的单位是整倍数的儒略年(一个儒略年是365.25日)。例如,J1986.0是J2000.0之前14*365.25天,对应JDE是2451545.0 - 14*365.25 = 2446431.50。
  J2000.0或J1986.0中的字母J,指明了一个标准日历的时间单位是儒略年。之前,星体位置使用的标准历元始于贝塞尔年。贝塞尔年是太阳年,年首是:从春风点起算黄经280度时(修正了光行差:-20".5),这一时刻总是在格里高公历年附近。贝塞尔年长度是太阳年,即365.2421988日(公元1900)。
  在新的系统中,为了区别基于贝塞尔年的旧历元,用字母B标识旧历元。例如:
    B1900.0 = JDE 2415020.3135 = 1900年1月0.8135
    B1950.0 = JDE 2433282.4235 = 1950年1月0.9235
  而
    J2000.0 = JDE 2451545.00 (正好是整数)
    J2050.0 = JDE 2469807.50 (正好是0.5)
    ……
  J1986.0或J2000.0等,它们后面的那个".0"表示年首。
------------------------
  严格的岁差计算方法:
  设T是J2000.0起算的儒略世纪数,设t是某一起始历元到终止历元之间的时间差,单位也是儒略世纪数。
  换句话说,如果(JD)o和(JD)是儒略日,分别对应初始和结束历元,就有:
    T = ( (JD)o - 2451545.0 )/36525
    t = ( (JD)  - (JD)o     )/36525
  那么我们有以下关于ζ、z、θ的数值表达式可用于精确坐标位置转换:
    ζ = ( 2306".2181 + 1".39656*T - 0".000139*T^2)*t
          + (0".30188 - 0".000344*T)*t^2 + 0".017998*t^3
     z = ( 2306".2181 + 1".39656*T - 0".000139*T^2)*t      ……20.2式
          + (1".09468 + 0".000066*T)*t^2 + 0".018203*t^3
    θ = ( 2004".3109 - 0".85330*T - 0".000217*T^2)*t
          - (0".42665 + 0".000217*T)*t^2 - 0".041833*t^3
  如果起始历元已经是J2000.0,那么T=0,则20.2式变为:
    ζ =  2306".2181*t + 0".30188*t^2 + 0".017998*t^3
     z =  2306".2181*t + 1".09468*t^2 + 0".018203*t^3      ……20.3式
    θ =  2004".3109*t - 0".42665*t^2 - 0".041833*t^3
  然后,这个严格的公式就可用于某一起始历元赤道坐标(αo和δo)与终止历元赤道坐标(α和δ)的转换:
    A = cos(δo)*sin(αo+ζ)
    B = cos(θ)*cos(δo)*cos(αo+ζ) - sin(θ)*sin(δo)
    C = sin(θ)*cos(δo)*cos(αo+ζ) + cos(θ)*sin(δo)    ……20.4式
    tan(α-z) = A/B      sin(δ) = C
  使和带两个参数的ATN2函数计算α-z,以保证得到正确的象限。也可以使用其它方法,见第1章。
  如果星体接近天极,使用cos(δ) = sqrt(A*A+B*B)代替sin(δ)=C。
  在把αo和δo转为α和δ之前,应先计算行星自行。
-------------------
  例20.b:——恒星θPersei在J2000.0的历元平赤道坐标为:
       αo = 2h 44m 11s.986     δo = +49°13'42".48
    该坐标中,它每年的自行速度是:
       +0s.03425   (赤经)
       -0".0895    (赤纬)
    把该坐标转到2028年11月19.19日TD历元平坐标中。
  解:初始历元J2000.0(JD 2451545.0),终止历元是 JD 2462088.69。因此,t = +0.288670500儒略世纪数,或28.8670500儒略年。
  我们先计算行星自行。在28.86705年中,变化:
    +0s.03425*28.86705 = +0s.989  (赤经)
    -0".0895 *28.86705 = -2".58   (赤纬)
  因此,J2000.0坐标中,历元2028年11月13.19日,该星的位置是:
    αo = 2h 44m 11s.986 + 0s.989 = 2h 44m 12s.975 = +41°.054063
    δo = +49°13'42".48 - 2".58  = +49°13'39".90 = +49°.227750
  从初始坐标为J2000.0起算,所以可能使用20.3式。t的值是+0.288670500,我们得到:
    ζ = +665".7627 = +0°.1849341
     z = +665".8288 = +0°.1849524
    θ = +578".5489 = +0°.1607080

     A = +0.43049405
     B = +0.48894849
     C = +0.75868586

  α-z = +41°.362262
    α = +41°.547214 = 2h 46m 11s.331
    δ = +49°.348483 = +49°20'54".54
--------------------
  练习:
  αUrsae Minoris在J2000.0平赤道坐标是:
    α = 2h 31m 48s.704    δ = +89°15'50".72
  该坐标中,星体的每年自行速度为:
    +0s.19877  (赤经)
    -0".0152   (赤纬)
  计算该星体的B1900.0,J2050.0及J2100.0平赤道坐标
  答案:
    B1900.0  α = 1h 22m 33s.90   δ = +88°46'26".18
    J2050.0       3h 48m 16s.43        +89°27'15".38
    J2100.0       5h 53m 29s.17        +89°32'22".18

  应当注意的是,式20.2及20.3适用在有限在的时间范围内,比如,如果我们计算32700年,我们发现该历元的αUMi将在赤纬-87°,这完全错误。
------------------------
  使用黄道坐标
  如果我们使用黄道坐标(黄经、黄纬)而不是赤道坐标(赤经、赤纬),那么可以使用以下严格的方法。
  T和t的含义与上述的相同,接下来计算:
    η = (47".0029 - 0".06603*T + 0".000598*T^2)*t
        +(-0".03302 + 0".000598*T)*t^2 +0".000060*t^3
    П = 174°.876384 + 3289".4789*T + 0".60622*T^2     ……20.5式
        -(869".8089 + 0".50491*T)*t + 0".03536*t^2
     P = (5029".0966 + 2".22226*T - 0".000042*T^2)*t
        +(1".11113 - 0".000042*T)*t^2 -0".000006*t^3
  参数η是起始历元黄道面到终元黄道面的夹角。
  如果起始历元是J2000.0,就有T=0,所以以上表达式化简为:
    η = 47".0029 *t  -0".03302*t^2  +0".000060*t^3
    П = 174°.876384 -869".8089*t   +0".03536 *t^2     ……20.6式
     p = 5029".0966*t +1".11113 *t^2 -0".000006*t^3
  那么,以下严格的公式,要用于起始历元黄道坐标转到终止历元黄道坐标:
    A' = cos(η)*cos(βo)*sin(П-λo) - sin(η)*sin(βo)
    B' = cos(βo)*cos(П-λo)
    C' = cos(η)*sin(βo) + sin(η)*cos(βo)*sin(П-λo)

    tan(p +П-λ) = A'/B'     sin(β) = C'
--------------------------
  旧的岁差要素
  正如以前所说的,为了星库及计算方便,现在使用J2000.0历元,单位是儒略年(365.25天)或儒略世纪(36525)。以前,贝塞尔年首由参考系时间决定,单位是太阳年世太阳年世纪数。
  然而,这不是旧系统(FK4)与新系统(FK5)的唯一区别。
  首先,FK4的赤经零点有个小错误(分点改正)。
  其次,在第22章还要讲到,地球公转造成恒星光行差,引起黄经偏移Δλ黄围编移Δβ:
    Δλ = -k*cos(Θ-λ)/cos(β) + e*k*cos(π-λ)/cos(β)
    Δβ = -k*sin(Θ-λ)*sin(β) + e*k*sin(π-λ)*sin(β)
  式中Θ是太阳经度。π是地球轨道近点经度,e是轨道离心率,K是光行差常数。
  现在,对于恒星,右边的第二项几乎是常数,因为e、π-λ、β随时间变化十分缓慢。正因如此,天文学上曾经把式中的第二项(也称E项修正)保留为恒星的平位置。
  目前,依赖地球椭圆轨道的E项修正不再属于恒星的平位置。相反,它们被折算到了视位置中(详见第22章)。
  基于Newcom岁差表达式的B1950.0平赤道分点坐标(FK4系统)转换到新的J2000.0的IAU系统(FK5系统),恒星平位置及自行的的转换步骤可以找到,例如在1984的天文年历中。
  岁差公式20.2及20.3适用于FK5系统。如果只有FK4位置及自行可供使用,那么可按下面方法着手计算恒星在FK5系统中的位置:
  1.使用Newcomb的岁差公式(见下面) 
  2.在平视位置转换时,E项修正应忽略。
  3.最后的恒星视赤经应加下如下修正:
      Δα = 0s.0775 + 0s.0850*T
    式中T是J2000.0起算的儒略世纪数
  Newcomb的岁差表达式是:
  设(JD)o和(JD)分别是起始历元和终止历元,那么:
    T = ( (JD)o - 2415020.3135 ) / 36524.2199
    t = ( (JD)  - (JD)o        ) / 36524.2199

    ζ = (2304".250 + 1".396*T)*t + 0".302*t^2 + 0".018*t^3
     z = ζ + 0"791*t^2 + 0".001*t^3
    θ = (2004".682 - 0".853*T)*t - 0".426*t^2 - 0".042*t^3
  如果起始历元是B1950.0,那么T=0.5,以上表达式变为:

    ζ = 2304".948*t + 0".302*t^2 + 0".018*t^3
     z = 2304".948*t + 1".093*t^2 + 0".019*t^3
    θ = 2004".255*t - 0".426*t^2 - 0".042*t^3


 


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


加好友 发短信
等级:超级版主 帖子:547 积分:4342 威望:0 精华:0 注册:2004/12/4 22:56:00
  发帖心情 Post By:2008/4/10 19:58:00

楼上的兄弟都是高人啊,图片点击可在新窗口打开查看我得好好学习一下!图片点击可在新窗口打开查看

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


加好友 发短信
等级:蝙蝠侠 帖子:866 积分:3824 威望:3 精华:20 注册:2008/3/20 22:14:00
  发帖心情 Post By:2008/4/10 21:18:00

第11章 恒星时与格林尼治时间
------------------------
  译者注:
  地球不断的自转着,天球子午圈时刻不断的变化着,我们必须找到适当的方法来标定子午圈在各时刻的位置。恒星时是天文学和大地测量学标定的天球子午圈位置的值,由于借用了时间的计量单位,所以常被误解为是一种时间单位。恒星时是根据地球自转来计算的,它的基础是恒星日。由于地球环绕太阳的公转运动,恒星日比平太阳日(也就是日常生活中所使用的日)短约1/365(相应约四分钟或一度)。
  本地恒星时的定义是一个地方的子午圈与天球的春分点之间的时角,各地方的经度不同,所以子午圈不同,因此地球上每个地方的恒星时都与它的经度有关。恒星时的参考点是春分点,所以春分点的变化也将对恒星时产生影响。由于地球的章动春分点在天球上并不固定,而是以18.6年的周期围绕着平均春分点摆动。因此恒星时又分真恒星时和平恒星时。真恒星时是通过直接测量子午线与实际的春分点之间的时角获得的,平恒星时则忽略了地球的章动。真恒星时与平恒星时之间的差异最大可达约0.4秒。
  一个地方的当地恒星时与格林尼治天文台的恒星时之间的差就是这个地方的经度。因此通过观测恒星时可以确定当地的经度(假如格林尼治天文台的恒星时已知的话)或者可以确定时间(假如当地的经度已知的话)。
  通过确定恒星时可以简化天文学的计算,比如通过恒星时和当地的纬度可以很方便地计算出哪些星正好在地平线以上。
  总之:
  以春分点作为基本参考点,由春分点周日视运动确定的时间,称为恒星时(简称ST)。某一地点的地方恒星时,在数值上等于春分点相对于这一地方子午圈的时角。
  以平太阳作为基本参考点,由平太阳周日视运动确定的时间,称为平太阳时(简称MT)。平太阳是美国天文学家纽康(S.Newcomb,1835 – 1909年)在十九世纪末做的一个假想参考点。

-----------------------------
  译者说明:一个时刻,通常由日期部分及时间部分组成。本文所述的时间指世界时UT,世(0h UT)表示世界时0点,(3h UT)表示世界时3点。如果没有特另申明,均指格林尼治子午圈时间。
  一个给定的世界时日期的(0h UT),格林尼治子午圈的恒星时可按如下计算。
  先计算当天(0h UT)对应的JD,它是以.5结束的数字。那么就有:
    T = (JD - 2451545.0) / 36525                 ……(11.1)
  那么(0h UT)的格林尼治平恒星时使用以下表达式计算,该式采用IAU1982的表达式:
    θo = 6h 41m 50s.54841 + 8640184s.812866*T   ……(11.2)
          + 0s.093104*T^2 - 0s.0000062*T^3
  上式表达为度单位,则公式改写为:
    θo = 100.46061837 + 36000.770053608*T       ……(11.3)
          + 0.000387933*T^2 - T^3/38710000
  重要的是,公式(11.2)及(11.3)仅在T对应(0h UT)时有效。
  为了找出任意UT时刻的恒星时,先计算UT*1.00273790935,再把计算结果加入(0h UT)时刻对应的平恒星时,即得到UT时刻的平恒星时。
  当平恒星时表达为度单位,也可使用下式直接得到任意UT刻的恒星时,式中的JD是任意时刻的儒略日数(不在要求JD对应在0点),T由公式(11.1)式计算。
    θo = 280.46061837 + 360.98564736629*(JD-2451545.0)
          + 0.000387933*T^2 - T^3/38710000
  如果要求高精度,那么要求计算机处理足够的有效数字位数。
  由(11.2)、(11.3)、(11.4)得到的恒星时是平恒星时,也就是平春分点起算的格林尼治时角(黄道与Date平赤道的交点)。
  要取得视恒星时(即真春点起算的格林尼治时角),应加上修正值Δψ*cos(ε),式中Δψ是黄经章动,ε是真黄赤交角(见21章)。这个修正值也就作赤经章动(也称作分点方程)。因为Δψ是个小量,所以ε误差10"也没关系。
  如果Δψ的单位是角秒(1度是3600秒),则以时秒(1小时是3600秒)为单位的修正值是:
    Δψ*cos(ε)/15
------------------------------
例11.a ——求1987年4月10日(0h UT)的“平”及“视”恒星时。
  解:
  1、求平恒星时
  该日期对应的儒略日是 JD 2446895.5,由公式(11.1)得:
    T = -0.127296372348
  那么我们利用(11.2)式可得:
    θo = 6h 41m 50s.54841 - 1099864.18158秒
  或者,可再简化:上式可以再加上86400秒(1天的秒数)的适当倍数,把结果转到0到86400秒以内:
  θo =  6h 41m 50s.54841 + 23335s.81842
      =  6h 41m 50s.54841 + 6h 28m 55s.81842
      = 13h 10m 46s.3668
  这就是所要求解的平恒星时。

  2、求视恒星时:
  由例21.a,我们得到此时的Δψ = -3".788及ε = 23°26'36".85,事实上这两个值是(0h TD)时刻的,而不是(0h UT),但我们忽略在ΔT = TD - UT期间内Δψ的微小差别。
  因此,赤经章动是 -3.788/15*cos(23°26'36".85) = -0s.2317,则所要计算的视恒星时为:
    13h 10m 46s.3668 -0s.2317 = 13h 10m 46s.1351
----------------------------------------------
例11.b ——1987年4月10日 19h 21m 00s UT时刻的平恒星时。

  首先,我们计算0h UT的平恒星时。我们得到 13h 10m 46s.3668(详见上一个例子)
  然后:
      1.00273790935 * (19h 21m 00s)
    = 1.00273790935 * 69660秒
    = 69850.7228秒
    = 19h 24m 10s.7228
  那么,要计算的平恒时是:
      13h 10m 46s.3668 + 19h 24m 10s.7228
    = 32h 34m 57s.0896
    =  8h 34m 57s.0896
  --------
  另外,我们也可以使用11.4式,直接计算。
  1987年4月10日 19h 21m 00s UT时,对应的儒略日数是:
        JD = 2446896.30625
  由(11.1)式,相应T的值为 -0.12727430。则,由公式(11.4)得:
        θo = -1677831°.2621266
  转到0到360度得:
        θo = 128°.7378734
  这就是要求计算的平恒星时,单位是度。除以15后,可变换为时角单位(因为1小时是15度):
        θo = 8h.58252489 = 8h 34m 57s.0896
  结果与上面的相同。
---------------------------------------------
第12章 坐标变换

  我们将使用到以下符号:
  α = 赤经。这个量一般表达为时间单位,也就是用时、分、秒表示。因此,在公式中使用到它时,应先转为“度”单位,必要时还要转为弧度单位。相反,如果使用公式计算出α,通常表达为“度”或“弧度”单位,应除转为小时单位(如果是度单位,除15即可),然后,还可以转为时分秒格式。
  δ = 赤纬。天赤道以北为正,以南为负。
  α1950 = 涉及B1950.0标准分点的赤经。
  δ1950 = 涉及B1950.0标准分点的赤纬。
  α2000 = 涉及J2000.0标准分点的赤经。
  δ2000 = 涉及J2000.0标准分点的赤纬。
  λ = 黄经。从春风点,沿黄道测量的经度。
  β = 黄纬。黄道以北为正,以南为负。
   l = 银经。
   b = 银纬。
   h = 地平纬度,地平线以上为正,以下为负。
   A = 地平经度(方位角)。由南向西测量。值得注意的是,航海家、气象学家的指南针方向(或地平经度),北方向为0度,东90度,南180度,西270度。但天文学家(注1)不同意,他们从南开始测量,因为时角也是从南开始测量的。因此,一个天体正好在子午圈的南方向,就有 A = H = 0度。
   ε= 黄赤交角。黄道与天赤道的夹角。平黄赤交角可用(21.2)式计算。然而,如果使用视赤经及视赤纬(受光行差及章动影响),那么计算时就要用到真黄赤交角ε+Δε(详见第21章)。如果α、δ是涉及J2000.0标准分点坐标的,那么该历元的ε就时ε2000=23°26'21".448 = 23°.439291。对于历元B1950.0标准分点坐标,ε1950=23°.4457889
  φ= 观测者(站)纬度,北半球为正,南半球为负。
  H = 本地时角,从南向本测量。

  如果θ是本地恒星时,θo是格林尼治恒星时,L是观者站经度(从格林尼治向西为正,东为负),那么本地时角计算如下:
    H = θ - α  或  H =θo - L - α
  如果α含章动效果,那么H也含章动(见11章)。

  要从赤道坐标转到黄道坐标,可以使用以下公式:
    tan(λ) = ( sin(α)*cos(ε) + tan(δ)*sin(ε) ) / cos(α) ……12.1式
    sin(β) = sin(δ)*cos(ε) - cos(δ)*sin(ε)*sin(α)       ……12.2式
----------------------------------------
  注1:William Chauvenet在它的《球面几何和实用天文学》(第5版,1981),卷I,第20页说到:地平经度的原点选取得很随意,所以他们计算出的方向也很随意。但天文学家通常选取地平的南点为原点,... 然而,航海家通常根据他们所在位置是北纬还是南纬来选择原点在北或在南。
  S.Newcom在它的《球面天文学概论》第95页中写道:“在实践中,可以从北点或南点测量,并且方向可以是东或西...”——所以说,伟大的美国天文学家没有特别选择。
----------------------------------------
  地理经度
  在这里,地理经度是从子午圈向西测量的,而不是向东。这个约定被多数天文学家认可长达1个世纪——见实例参考1—6。例如:华盛顿经度,D.C.,+77°04';奥地利维也纳经度是:-16°23'。
  我们不能理解,为什么IAU(国际天文联合会)最初决所有的行星地理经度从它们自转轴相反的方向测量,于是,1982年为地球修改了这个系统。我们将不跟从IAU的决定,我们将考虑西经为正。其它行星也遵照这个系统,向西测量为正,这正是为什么它们的中心子午圈经度和地球一样,随时间不断增加。
-----------------------------------------
  赤道坐标转到黄道坐标:
    tan(α) = ( sin(λ)*cos(ε) - tan(β)*sin(ε) ) / cos(λ) ……12.3式
    sin(δ) = sin(β)*cos(ε) + cos(β)*sin(ε)*sin(λ)       ……12.4式
  计算本地地平坐标:
    tan(A) = sin(H)/( cos(H)*sin(φ) - tan(δ)*cos(φ) )      ……12.5式
    sin(h) = sin(φ)*sin(δ) + cos(φ)*cos(δ)*cos(H)         ……12.6式
  如果希望地平经度A是从北点开始计算的,而不是上面的南点起算,则只须对(12.5)式算出的A加上180°即可。
  地平坐标转到赤道坐标:
    tan(H) = sin(A)/( cos(A)*sin(φ) +tan(h)*cos(φ) )
    sin(δ)= sin(φ)*sin(h) - cos(φ)*cos(h)*cos(A)
----------------------------------
  当前的银河系统坐标已在1959年IAU中定义了。在B1950.0标准赤道系统,银河(银河系)北极的坐标是:
    α1950 = 12h 49m = 192°.25,  δ1950 = +27°.4
  银经的原点在银道上,该点距银道与B1950.0赤道升点33度。
  这些值都是固定的惯例值,因此我们还须严格考虑它们与B1950.0赤道坐标的关系。
·从B1950.0标准分点赤道坐标转到银道坐标:
    tan(x) = sin(192.25-α) /
           ( cos(192.25-α)*sin(27.4) - tan(δ)*cos(27.4) )    ……12.7式
         l = 303° - x
    sin(b) = sin(δ)*sin(27.4) + cos(δ)*cos(27.4)*cos(192.25-α) …12.8式
   式中角度单位是度。
·从银道坐标转到B1950.0标准分点坐标:
    tan(y) = sin(l-123) /
           ( cos(l-123)*sin(27.4) - tab(b)*cos(27.4) )
        α = y + 125°
   sin(δ) = sin(b)*sin(27.4) + cos(b)*cos(27.4)*cos(l-123)
  如果给定恒星2000.0的平位置而不是1950.0平位置,那么在使用公式(12.7)及(12.8)之前,应将α2000转为α1950,将δ2000转为δ1950,详见第20章。
--------------------------------------
  公式(12.1)、(12.3)等,分别tan(λ)、tan(α)等的值,进而使用反正切函数算出λ、α等。然而,直接使用反正切函数将造成无法确定这些角度所在的象限,会产生不确定的180度角问题。可直接使用二参数函数ATN2,分子及分母作为入口参数传入。参见第1章的“正确象限”
--------------------------------------
  例12.a:——计算恒星Pollux(β Gem)的黄道坐标,它的赤道坐标是:
    α2000 = 7h 45m 18s.946,    δ2000 = +28°.026183
  解:
  使用这些值α=116°.328942,δ=+28°.026183,和ε=23°.4392911,由公式(12.1)及(12.2)得:
  tan(λ) = (+1.03403986)/(-0.44352398),因此λ=113°.215630
      β  = +6°.684170
  因为α、δ涉及J2000.0标准分点坐标,所以λ、β也涉及同样的分点。
--------------------------------------
  练习:
  利用上面计算出的λ和β值,利用12.3式及12.4式反算出α和δ
--------------------------------------
  例12.b:——计算在1987年4月10日(19h 21m 00s UT)时刻,在华盛顿U.S Naval天文台(经度 = +77°03'56" = 5h 08m 15s.7,纬度 = +38°55'17")金星的地平经度及纬度。此时金星的赤道视坐标是:α= 23h 09m 16s.641,δ= -6°43'11".61
  解:
  这是行星的视赤经及视赤纬。我们需要此刻的视恒星时。
  我们先计算格林尼治1987年4月10日(19h 21m 00s UT)的平恒星时,得到的值是:
    8h 34m 57s.0896(详见例11.b)
  利用第21章描述的方法,我们得到该时刻的以下值:
    黄经章动:  Δψ = -3".868
    真黄赤交角:  ε = 23°26'36".87
  格林尼治视恒星时是:
    θo = 8h 34m 57s.0896 + (-3.868/15)*cos(ε)秒 = 8h 34m 56s.853
  金星在华盛顿的时角:
    H = θo - L - α
      = ( 8h 34m 56s.853) - (5h 08m 15s.7) - (23h 09m 16s.641)
      = -19h 42m 35s.488 = -19h.7098578 = -295°.647867
      = +64°.352133
  由公式(12.5)和(12.6)得:
    tan(A) = (+0.9014712)/(+0.3636015)  因此 A = +68°.0337
        h  = +15°.1249
  所以,行星在地平线上面15度,在西南方向到正西方向之间。
  应注意到,公式(12.6)没有考虑大气折射的影响,也没有考虑行星视差。大气折射问题参见第15章。视差改正将在第39章中研究。
  作为一个练习,请计算Nova Serpentis 1978的银道坐标。它的赤道坐标是:
   α1950 = 17h 48m 59s.74,  δ1950 = -14°43'08".2
  答案: l = 12°.9593,  b = +6°.0463
--------------------------------------------
  黄道和地平线

  如果ε是黄赤交角,φ是观测站的纬度,θ是本地恒星时,那么在地平线(地面与天球交的大圆)与黄道相交的两点的黄经是(注意:只需计算1点,另一点与它相距180度):
    tan(λ) = (-cos(θ)) / (sin(ε)tan(φ) + cos(ε)*sin(θ))   ……(12.9)式
  黄道面与地平面的夹是:
    cos(I)  = cos(ε)*sin(φ) - sin(ε)*cos(φ)*sin(θ)         ……(12.10)式
  在一个恒星日期间,角度I在两个极限值之间。例如,纬度48°00'N,ε=23°26',那么I的两个极限值是:
    θ = 90°时, I = 90°- φ + ε = 65°26'
    θ =270°时, I = 90°- φ - ε = 18°34'
  要注意的是,I不是太阳周日运动的平面与地平面的夹角。译者注:而是周年运动平面(黄道面)与地平的夹角,当然上式还未涉及视差问题。
------------------------------------------
  例12.c:——ε=23°.44,φ=+51°,θ = (5h 00m) =75°,
  利用公式(12.9),得到:
    tan(λ) = -0.1879,因此 λ =169°21'和 λ = 349°21'
  利用(12.10)式得:
    I = 62°
------------------------------------------
  参考资料
1、《1835年航海年历和天文历书》,第508页(伦敦,1833)。
2、《1857年美国星历及航海年历》,第491页(华盛顿,1854)
3、《1960年天文历书》,第434页等(伦敦,1958)
4、W.Chauvenet,《球面学和实用天文学手册》,卷I,第317页等(费城,1891)
5、A.Danjon,《Astronomie Generale》,第46页(巴黎,1959)
6、S.Newcomb,《球面天文学概论》,第119页(纽约,1906)

 

 


 


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


加好友 发短信
等级:蝙蝠侠 帖子:866 积分:3824 威望:3 精华:20 注册:2008/3/20 22:14:00
  发帖心情 Post By:2008/4/10 21:34:00

第22章和46章

 

 下载信息  [文件大小:   下载次数: ]
图片点击可在新窗口打开查看点击浏览该文件:


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


加好友 发短信
等级:版主 帖子:898 积分:7341 威望:10 精华:1 注册:2005/3/17 11:37:00
  发帖心情 Post By:2008/4/16 9:08:00

学习了。

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


加好友 发短信
等级:蝙蝠侠 帖子:866 积分:3824 威望:3 精华:20 注册:2008/3/20 22:14:00
  发帖心情 Post By:2008/4/17 21:28:00

花了一些时间讨论一下天文坐标的变换

 

 下载信息  [文件大小:   下载次数: ]
图片点击可在新窗口打开查看点击浏览该文件:


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


加好友 发短信
等级:论坛游民 帖子:54 积分:1018 威望:0 精华:0 注册:2007/6/16 13:30:00
  发帖心情 Post By:2008/4/23 19:57:00

剑伟兄是越来越厉害了,呵

 

建议春光老大重开一贴,将剑伟兄的成果整理一下统一发出来


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


加好友 发短信
等级:论坛游民 帖子:54 积分:1018 威望:0 精华:0 注册:2007/6/16 13:30:00
  发帖心情 Post By:2008/4/23 20:28:00

以下是引用春光在2008-3-26 23:22:00的发言:
      其实,如果能找到英文版的,也总比没有强啊,毕意不是无米之炊,不是闭门造车,其实我国教育部门不应主张全民学外语,而是应重点重视一下培养一些职业翻译,就象美国一样,术业有专攻嘛,再重点翻译国外的一些好书,把国外最好的科学知识介绍到中国来。

        中国广大民众是不需要现代天文计算法吗?我想不是,因为我国使用的是农历,是一个天文年历性质的历法,计算朔望两弦,及节气等就是一个迫切而又巨大的需求。解决历算中最关键的现代历书天文学中的星历表计算理论书,是少之又少,我至今没有发现,而且我发现国内人写的技术专业方面的书,从来不细细地教您作,都是空洞的理论,在有些方面上绕来绕去的,实际使用中用到又很少。就是不讲实际如何怎样做,具体步骤如何,读过和不读差不多(我真实的体会,可能是以偏概全了),尤其是一些专业技术的书,更是如此,我想可能是国内编书的专家时间太忙吧,抽不出时间把书编好,不是我对国内的书有偏见,您想如果国内的书好用,谁还去费力看国外的书,还得需要外语基础;我读过国外的有关技术方面的书,特别是欧美国家作者写的书真是很好,读过以后,就能按书上的步骤就能做了,实际水平也就提高了一大截,而且作者还在书中告诉您一些细节需要注意的问题,注意读这方面国外高科技的书时尽量要读外文原版的,因为这是最准确的。总之,我认为如果是技术书最好是读国外的书。

严重同意春光的意见!

不要说别的,单单就现行农历的详细规则我就在网上找了好久,除了春光的那个100个问答里面有介绍外,根本找不到完整的中文资料,所有的中文资料全都是简单的介绍,没有任何实用价值,倒是一个法国的网站上面有英文的中国农历算法详细介绍,实在让人感到心寒啊

这是那个网站:

http://www.chinesefortunecalendar.com/clc/Default.htm


支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
总数 106 上一页 1 2 3 4 5 6 7 8 9 10 下一页 ..11

返回版面帖子列表

我自己制作的万年历








签名