太阳位置的计算

太阳位置的计算：参考2007年中国天文年历（基本历表采用DE405）

一、时间计算

1.求出儒略日JD

在计算时，采用J2000.0（2000年1月1.5日TDB，即2000年1月1日12h，即儒略日JD2451545.0TDB）作为起算点，也就是基准历元。

下面的方法对正数年和负数年都有效，负的儒略日数除外。

设Y为给定年份，M为月份，D为该月日期（可以带小数）。

若M2，Y和M不变.

若 M1 或2，以Y1代Y，以M12代M。

换句话说，如果日期在1月或2月，则被看作是在前一年的13月或14月。

对格里高利历，有 AINT(

对儒略历，取B0。

要求的儒略日JD为：

JDINT(365.25(Y4716))INT(30.6001(M1))DB1524.5 YA)，B2AINT() 1004

2.从J2000.0 起算的儒略世纪数T

TJD2451545.0

36525

二、太阳平均轨道根数

1.太阳赤道偏心率

es0.01670862-0.00004204T-0.000000124T2

2.太阳平近点角

s'''''''''M03573144.76129596581.04T-0.562T2-0.012 T3

s''''''M0357.5291129596581.04T-0.562T2-0.012 T3

3.平黄赤交角0

'''''''''0232621.44846.815T0.00059T20.001813T3

''''''023.439291146.815T0.00059T20.001813T3

4.月亮升交点的平黄经n

''''n125.044555016962890''.5431T7.4722T20.007702T3

5.对于当天平春分点的太阳几何平黄经Ls

'''''''2Ls

02802759.21129602771.36T1.093T

''''2Ls

0280.4664472129602771.36T1.093T

6.轨道半长轴

a1.00000102

三、其他所需的天体轨道根数

1.日月平角距D

''''''D297.850195471602961601.2090T6.3706T20.006593T3

2.金星的平近点角Mv

Mv50.4160941.[1**********]d0.0013835T2

3.火星的平近点角Mm

Mm19.3730410.[1**********]d0.0001759T2

其中，d是相应的儒略日。

dJD2451545.0

四、太阳轨道根数的改正

1.太阳平近点角的摄动改正Ms

''''Ms6.4sin(251.420.2T)1.9sin(207.5150.3T)

2.黄赤交角章动改正

''9.2052331cosn

3.太阳黄经章动改正

''17.2sinn

4.中心差v

es351113103s4v(2e)sinMs(es2es4)sin2Mses3sin3Msesin4Ms 44241296s

5.太阳平黄经的主要摄动项

（1）长周期项

''''Ll6.4sin(251.420.2T)1.9sin(207.5150.3T)

（2）月亮摄动项

''Ln6.6sinD

（3）大行星主要摄动项

''''''Lp5.5cos(148.32Mv2Ms)4.8cos(299.1MvMs)2.0cos(343.92Mm2Ms)

6.太阳黄经光行差的改正

''20.4955Lg rs

1es2

其中太阳向径r1.[1**********] s1ecosfs

太阳真近点角fMsv

五、太阳轨道参数

1.太阳平近点角

sMsM0Ms

2.黄赤交角

0

3.太阳视黄经

sLs

0vLlLnLpLg

六、太阳赤道视坐标

由于太阳的黄纬很小，几乎为零，因此可以采用球面直角三角形公式直接将太阳视黄经转换为视赤经和视赤纬。

1.太阳赤经 RAarctan(costans)

2.太阳赤纬 Decarcsin(sinsins)

1es2

其中：当coss0时，太阳向径r1.[1**********] s1ecosfs

参考文献：

[1]张 英，现代化的天文定位——计算机代替天文定位中的查表，中国水运，第07卷第01 期

[2]戴冉，基于新历元的高精度太阳坐标模型，大连海事大学学报，第22卷第2期

[3]中国科学院紫金山天文台，2007年中国天文年历，北京：科学出版社，2007.1

[4] Ibrahim Reda *, Afshin Andreas，2004，Solar position algorithm for solar radiation applications，Solar Energy 76 (2004) 577–589

太阳位置的计算：参考2007年中国天文年历（基本历表采用DE405）

一、时间计算

1.求出儒略日JD

在计算时，采用J2000.0（2000年1月1.5日TDB，即2000年1月1日12h，即儒略日JD2451545.0TDB）作为起算点，也就是基准历元。

下面的方法对正数年和负数年都有效，负的儒略日数除外。

设Y为给定年份，M为月份，D为该月日期（可以带小数）。

若M2，Y和M不变.

若 M1 或2，以Y1代Y，以M12代M。

换句话说，如果日期在1月或2月，则被看作是在前一年的13月或14月。

对格里高利历，有 AINT(

对儒略历，取B0。

要求的儒略日JD为：

JDINT(365.25(Y4716))INT(30.6001(M1))DB1524.5 YA)，B2AINT() 1004

2.从J2000.0 起算的儒略世纪数T

TJD2451545.0

36525

二、太阳平均轨道根数

1.太阳赤道偏心率

es0.01670862-0.00004204T-0.000000124T2

2.太阳平近点角

s'''''''''M03573144.76129596581.04T-0.562T2-0.012 T3

s''''''M0357.5291129596581.04T-0.562T2-0.012 T3

3.平黄赤交角0

'''''''''0232621.44846.815T0.00059T20.001813T3

''''''023.439291146.815T0.00059T20.001813T3

4.月亮升交点的平黄经n

''''n125.044555016962890''.5431T7.4722T20.007702T3

5.对于当天平春分点的太阳几何平黄经Ls

'''''''2Ls

02802759.21129602771.36T1.093T

''''2Ls

0280.4664472129602771.36T1.093T

6.轨道半长轴

a1.00000102

三、其他所需的天体轨道根数

1.日月平角距D

''''''D297.850195471602961601.2090T6.3706T20.006593T3

2.金星的平近点角Mv

Mv50.4160941.[1**********]d0.0013835T2

3.火星的平近点角Mm

Mm19.3730410.[1**********]d0.0001759T2

其中，d是相应的儒略日。

dJD2451545.0

四、太阳轨道根数的改正

1.太阳平近点角的摄动改正Ms

''''Ms6.4sin(251.420.2T)1.9sin(207.5150.3T)

2.黄赤交角章动改正

''9.2052331cosn

3.太阳黄经章动改正

''17.2sinn

4.中心差v

es351113103s4v(2e)sinMs(es2es4)sin2Mses3sin3Msesin4Ms 44241296s

5.太阳平黄经的主要摄动项

（1）长周期项

''''Ll6.4sin(251.420.2T)1.9sin(207.5150.3T)

（2）月亮摄动项

''Ln6.6sinD

（3）大行星主要摄动项

''''''Lp5.5cos(148.32Mv2Ms)4.8cos(299.1MvMs)2.0cos(343.92Mm2Ms)

6.太阳黄经光行差的改正

''20.4955Lg rs

1es2

其中太阳向径r1.[1**********] s1ecosfs

太阳真近点角fMsv

五、太阳轨道参数

1.太阳平近点角

sMsM0Ms

2.黄赤交角

0

3.太阳视黄经

sLs

0vLlLnLpLg

六、太阳赤道视坐标

由于太阳的黄纬很小，几乎为零，因此可以采用球面直角三角形公式直接将太阳视黄经转换为视赤经和视赤纬。

1.太阳赤经 RAarctan(costans)

2.太阳赤纬 Decarcsin(sinsins)

1es2

其中：当coss0时，太阳向径r1.[1**********] s1ecosfs

参考文献：

[1]张 英，现代化的天文定位——计算机代替天文定位中的查表，中国水运，第07卷第01 期

[2]戴冉，基于新历元的高精度太阳坐标模型，大连海事大学学报，第22卷第2期

[3]中国科学院紫金山天文台，2007年中国天文年历，北京：科学出版社，2007.1

[4] Ibrahim Reda *, Afshin Andreas，2004，Solar position algorithm for solar radiation applications，Solar Energy 76 (2004) 577–589