Example 9. 月の位置計算問題




Example 9. V9 p.239.
Computation of moon's latitude and longitude from the tables for a given date.
Example:(Cf. VII3 p.336).
      Date: Nabonassar 466, Thoth[I] 7/8, 2 equinoctial hours after midnight.

問題は「Nabonassar 466年1月7/8日(7日の夜という意味)夜半から2時間後における月の位置を求めよ」です。


月の位置の計算手順は以下です。(「Ptolemy's Almagest」p.239)
 1. 月のパラメータのそれぞれの1日の平均速度及び初期値(Epoch)は以下です。
     パラメータ    平均速度        初期値(Epoch)
      λo      (13;10,34,58,33,30,30°)    41;22°
      αo      (13; 3,53,56,17,51,59°)   268;49°
      ωo      (13;13,45,39,48,56,37°)   354;15°
      ηo      (12;11,26,41,20,17,59°)    70;37°

 2. 位置を計算する日までの日数を求める。(day)
 3. 期間毎の月平均運動表を読み取り、初期値も含め足し算する。
      Table of the moon's Longitude (λ)
      Table of the moon's Anomaly (α)
      Table of the moon's Latitudee (ω)
      Table of the moon's Elongation (η)

 4. 項目3で求めた値Elongation (η)を2倍し以下の[Table V8. Table of the complete Lunar
   Anomaly]の"Equaqtion for Apogee(3)"からパラメータ"c3"、"Sixtieths(6)"から
   パラメータ"c6"を読む。
   (2η<=180° ⇒c3は負。2η>180° ⇒c3は正)
 5. 項目4で求めた値c3を Anomaly(α)に加え真アノマリ(α=α+c3)を求める
 6. 項目5で求めた真アノマリ(αで以下の[Table V8. Table of the complete Lunar Anomaly]の
   "Epicyclic Equation(4)"からパラメータ"c4"、"Increment in Epicyclic(5)"から
   パラメータ"c5"を読む。
   (α<=180° ⇒c4,c5は負。α>180° ⇒c4,c5は正)
 7. 計算したc4,c5,c6よりc=c4+c5×(c6/60)を求める。
 8. 項目3で求めたλにcを加えたものが月黄経(λλ + c)となる。

 9. 項目3で求めたωにcを加えω(=ω+c)を求める。
 10. 項目9で求めたωで以下の[Table V8. Table of the complete Lunar Anomaly]の
   "Latitude(7)"から月の黄緯(β)を求める。
    (90°≦α≦270° ⇒βは負。α<90°及びα>270° ⇒βは正)
 
 
 例題の解法
 "Nabonassar 466年1月7/8日(7日の夜という意味)夜半後2時"より
 暦元からの日数はそれぞれ1づつ引いて、450年+15年+0月+6日+(14)時間
 それぞれの平均運動の表から秒以上を読み取ると、
 
 
項目
λ α ω η
初期値 41;22 268;49 354;15 70;37
450年 260;46,44 323;26,5 320;54,6 10;11,3
15年 140;41,33 250;46,52 70;41,48 144;20,22
0月 0 0 0 0
6日 79;3,30 78;23,24 79;22,34 73;8,40
14時 7;41,10 7;37,16 7;43,2 7;6,41
合計 169;35°
(169;34,57)
209;3°
(209;2,37)
112;56°
(112:56,30)
305;24°
(305;23,46)

(プログラムで計算する場合は1日の平均速度×(期間の合計)で差は出ない。)
 
 従って2η=250;48°
 Col.3: c3(2η)=-(13;6 -(13;6-13;3)/3×(250;48-249) )=-(13;6 -0;2 )=-13;4°
 Col.6: c6(2η)=37;50 -(37;50-36;14)/3×(250;48-249) =37;50-0;58=36;52°

 α=α+c3=209;3-13;4=195;59°
 Col.4: c4(α)=1;25 +(1;41-1;25)/3×(195;59-195) =1;25+0;5=1;30°
 Col.5: c5(α)=0;52 +(1;2-0;52)/3×(195;59-195) =0;52+0;3=0;55°
 従って
 c==c4+c5×c6=1;30+0;55×(36;52)/60=2;04°
 月黄経(λλ + c)=169;35 +2;04=171;39°【「Almagest」本文では171;30°】
 
 ω=ω+c=112;56+2;04=115;0°
 Col.7: 月の黄緯(β)= c7(ω)=-(2;2 +(2;16-2;2)/3×(115-114))= -(2;2 +0;5)=-2;7°
                                        【「Almagest」本文では-2 1/6°】


Table V8. Table of the complete Lunar Anomaly
 
No.
(1)

 No.
(2)
Equaqtion
for Apogee
(3)
Epicyclic
Equation
(4)
Increment
in Epicyclic
(5)
 
Sixtieths
(6)
 
Latitude
(7)
03600;00;00;00;05;0
63540;530;290;140;124;58
123481;460;570;280;244;54
183422;391;250;421;204;45
243363;311;530;562;164;34
303304;232;191;103;244;20
363245;152;441;234;324;3
423186;73;81;356;253;43
483126;583;311;458;183;20
543067;483;511;5410;222;56
603008;364;82;312;262;30
662949;224;242;1115;52;2
7228810;64;382;1817;441;33
7828210;484;492;2520;341;3
8427611;274;562;3123;240;32
9027012;04;592;3526;360;0
9326712;155;02;3728;120;16
9626412;285;12;3829;490;32
9926112;395;02;3931;250;48
10225812;484;592;3933;11;3
10525512;564;572;3934;371;17
10825213;34;532;3836;141;33
11124913;64;492;3837;501;48
11424613;94;442;3739;262;2
11724313;74;382;3541;22;16
12024013;44;312;3242;382;30
12323712;594;242;2844;32;43
12623412;504;162;2445;282;56
12923112;364;72;2046;533;8
13222812;163;572;1648;183;20
13522511;543;462;1149;323;32
13822211;293;352;550;453;43
14121911;23;231;5851;593;53
14421610;333;101;5153;124;3
14721310;02;571;4354;34;11
1502109;222;431;3554;544;20
1532078;382;281;2755;454;27
1562047;482;131;1956;364;34
1592016;561;571;1157;154;40
1621986;31;411;257;554;45
1651955;81;250;5258;354;50
1681924;111;90;4259;44;54
1711893;120;520;3159;264;56
1741862;110;350;2159;374;58
1771831;70;180;1059;494;59
1801800;00;00;060;05;0

(「Ptolemy's Almagest」p.238) 2016/03/19作成


2016/03/19 Up
2016/03/24 修正

Copyright(C) 2016 Shinobu Takesako
All rights reserved