ここでは太陽黄経を求める方法を説明します。

太陽黄経を求める

2000年1月1日12時(世界時)からのユリウス世紀
t = (観測時のユリウス通日 - 2451545) / 36525

太陽の視黄経
ψ = Σ_i=1¹⁵P_i * cos(Q_i * t + R_i) + Σ_i=16¹⁸P_i * t * cos(Q_i * t + R_i) + 36000.7695° * t + 280.4602°

係数表

i	P (°)	Q (°)	R (°)
1	1.9147	35999.05	267.52
2	0.0200	71998.10	265.10
3	0.0020	32964.00	158.00
4	0.0018	19.00	159.00
5	0.0018	445267.00	208.00
6	0.0015	45038.00	254.00
7	0.0013	22519.00	352.00
8	0.0007	65929.00	45.00
9	0.0007	3035.00	110.00
10	0.0007	9038.00	64.00
11	0.0006	33718.00	316.00
12	0.0005	155.00	118.00
13	0.0005	2281.00	221.00
14	0.0004	29930.00	48.00
15	0.0004	31557.00	161.00
16	-0.0048	35999.00	268.00
17	0.0048	1934.00	145.00
18	-0.0004	72002.00	111.00

太陽の幾何学的黄経
λ = ψ + 0.0057°

参考文献
太陽視黄経の計算式
太陽視赤緯，均時差計算に関する一考察 (秋田県立大学) p.8 (16スライド目)
http://news-sv.aij.or.jp/kankyo/s13/OLDHP/matsu0512ppt.pdf

太陽黄経を求めるプログラム

観測日0時(世界時)の太陽について、太陽の幾何学的黄経の位置変化を一ヶ月分表示します。

  年  月