太阳位置计算

太阳位置计算

ID:37627282

大小:327.00 KB

页数:14页

时间:2019-05-26

太阳位置计算_第1页
太阳位置计算_第2页
太阳位置计算_第3页
太阳位置计算_第4页
太阳位置计算_第5页
资源描述:

《太阳位置计算》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、假如给出当前位置的经度、纬度和准确的天文时间,用什么方式计算出太阳当前的位置呢?例如算出类似于方位角和高度角的数据?实际中使用希腊字符表示经纬度,不过为了输入方便,我改用XYZ表示.X是地方纬度,Y是太阳赤纬,t是时角,z=90度-h高度h,方位角A,则sinh=sinXsinY+cosXcosYcostsinA=cosYsint/coshcosA=(sinXsinh-sinY)/(cosXcosh)具体的可以在网上搜示意图或这本书《天文学新概论》苏宜编著华中理工大学出版社2000年8月第一版P52的章节第24章太阳位

2、置计算[许剑伟于家里2008-3-30下午]一、低精度计算: 当计算精度要求为0.01度,计算太阳位置时可假设地球运动是一个纯椭圆,也就说忽略月球及行星摄动,计算表达如下。 设JD是儒略日数,可以用第7章表述的方法计算。T为J2000起算的儒略世纪数:   T=(JD-2451545.0)/36525 计算时要保留足够的小数位数,5位小数是不够的(除非所需的太阳黄经的精度要求不高),注意,T表达为儒略世纪数,所以T误差0.00001相当于0.37日。 接下来, 太阳几何平黄经:Lo=280°.46645+36000°.

3、76983*T+0°.0003032*T^2(Date平分点起算) 太阳平近点角:M=357°.52910+35999°.05030*T-0°.0001559*T^2-0°.00000048*T^3 地球轨道离心率:e=0.016708617-0.000042037*T-0.0000001236*T^2 太阳中间方程:C=+(1°.914600-0°.004817*T-0°.000014*T*T)*sin(M)+(0°.019993-0°.000101*T)*sin(2M)+0°.000290*sin(3M) 那么,太

4、阳的真黄经是:Θ=Lo+C 真近点角是:v=M+C 日地距离的单位是"天文单位",距离表达为:R=1.000001018(1-e^2)/(1+e*cos(v))……24.5式 式中的分子部分的值变化十分缓慢。它的值是:0.99971901800年0.99972041900年0.99972182000年0.99972322100年 太阳黄经Θ可由上述的方法算出,它是Date黄道分点坐标中的真几何黄经,需通过计算地心坐标星体位置也可算出。 要取得Date黄道坐标中太阳的视黄经λ,还应对Θ进行章动修正及光行差修正。如果精度要

5、求不高,可用下式修正:   Ω=125°.04-1934°.136*T   λ=Θ-0°.00569-0°.00478*sin(Ω) 某此时候,我们需要把太阳黄经转到J2000坐标中,在1900-2100年范围内可利用下式进行:Θ2000=Θ-0°.01397*(year-2000) 如果还想取得更高的转换精度(优于0.01度),那么你可以使用第25章的方法进行坐标旋转。 Date黄道坐标中的太阳黄纬不超过1".2,如果对精度要求不是很高,可以置0。因此,太阳的地心赤经α及赤纬δ可以用下式(24.6式,24.7式)计算

6、,式中ε是黄赤交角(由21章的21.2式计算)。 tanα=cosεsinΘ/cosΘ……24.6式 sinδ=sinεsinΘ……24.7式 如果要想得到太阳的视赤经及赤纬,以上二式中的Θ应换为λ,ε应加上修正量:   +0.00256*cos(Ω)[译者注]:实际上就是对Θ补上黄经章动及光行差,ε补上交角章动后再转到赤道坐标中。也可在赤道坐标中补章动及光行差,但公式不同。 公式24.6当然可以转为:tan(α)=cos(ε)*tan(Θ),接下来,我们要注意α与Θ应在同一象限。然而,如果你使用计算机语中有ATN2函

7、数(C语言是atan2),那最好保持24.6式不变,这样就可直接利用ATN2函数算出α,即:α=ATN2(cos(ε)*sin(Θ),cos(Θ))例24.a——计算1992-10-13,0点,即力学时TD=JDE2448908.5时刻的太阳位置。我们算得:T=-0.072183436Lo=-2318°.19281=201°.80719M=-2241°.00604=278°.99396e=0.016711651C=-1°.89732Θ=199°.90987=199°54'36"R=0.99766Ω=264°.65λ=1

8、99°.90897=199°54'32"εo=23°26'24".83=23°.44023(由21章的21.2式算得)ε=23°.43999α视=-161°.61918=+198°.38082=13h.225388=13h13m31s.4δ视=-7°.78507=-7°47'06"使用VSOP87行星理论计算出的的正确值是:(请与

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。