太阳位置的计算:参考2007年中国天文年历(基本历表采用DE405)一、时间计算1.求出儒略日JD在计算时,采用J2000.0(2000年1月1.5日TDB,即2000年1月1日12h,即儒略日2451545.0JDTDB)作为起算点,也就是基准历元。下面的方法对正数年和负数年都有效,负的儒略日数除外。设Y为给定年份,M为月份,D为该月日期(可以带小数)。若2M,Y和M不变.若1M或2,以1Y代Y,以12M代M。换句话说,如果日期在1月或2月,则被看作是在前一年的13月或14月。对格里高利历,有()100YAINT,2()4ABAINT对儒略历,取0B。要求的儒略日JD为:(365.25(4716))(30.6001(1))1524.5JDINTYINTMDB2.从J2000.0起算的儒略世纪数T2451545.036525JDT二、太阳平均轨道根数1.太阳赤道偏心率20.01670862-0.00004204-0.000000124seTT2.太阳平近点角'''''''2''303573144.76129596581.04-0.562-0.012sMTTT''''2''30357.5291129596581.04-0.562-0.012sMTTT3.平黄赤交角0'''''''2''30232621.44846.8150.000590.001813TTT''''2''3023.439291146.8150.000590.001813TTT4.月亮升交点的平黄经n''2''3125.044555016962890''.54317.47220.007702nTTT5.对于当天平春分点的太阳几何平黄经0sL'''''''202802759.21129602771.361.093sLTT''''20280.4664472129602771.361.093sLTT6.轨道半长轴1.00000102a三、其他所需的天体轨道根数1.日月平角距D''''2''3297.850195471602961601.20906.37060.006593DTTT2.金星的平近点角vM250.4160941.602130343640.0013835vMdT3.火星的平近点角mM219.3730410.524020682190.0001759mMdT其中,d是相应的儒略日。2451545.0dJD四、太阳轨道根数的改正1.太阳平近点角的摄动改正sM''''6.4sin(251.420.2)1.9sin(207.5150.3)sMTT2.黄赤交角章动改正''9.2052331cosn3.太阳黄经章动改正''17.2sinn4.中心差v3243451113103(2)sin()sin2sin3sin444241296sssssssssseveMeeMeMeM5.太阳平黄经的主要摄动项(1)长周期项''''6.4sin(251.420.2)1.9sin(207.5150.3)lLTT(2)月亮摄动项''6.6sinnLD(3)大行星主要摄动项''''''5.5cos(148.322)4.8cos(299.1)2.0cos(343.922)vsvsmspLMMMMMM6.太阳黄经光行差的改正''20.4955gsLr其中太阳向径211.000001017781cosssseref太阳真近点角sfMv五、太阳轨道参数1.太阳平近点角0sssMMM2.黄赤交角03.太阳视黄经0sslnpgLvLLLL六、太阳赤道视坐标由于太阳的黄纬很小,几乎为零,因此可以采用球面直角三角形公式直接将太阳视黄经转换为视赤经和视赤纬。1.太阳赤经arctan(costan)sRA2.太阳赤纬arcsin(sinsin)sDec其中:当cos0s时,太阳向径211.000001017781cosssseref参考文献:[1]张英,现代化的天文定位——计算机代替天文定位中的查表,中国水运,第07卷第01期[2]戴冉,基于新历元的高精度太阳坐标模型,大连海事大学学报,第22卷第2期[3]中国科学院紫金山天文台,2007年中国天文年历,北京:科学出版社,2007.1[4]IbrahimReda*,AfshinAndreas,2004,Solarpositionalgorithmforsolarradiationapplications,SolarEnergy76(2004)577–589