- 1、本文档共14页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
营销研究营销研究太阳位置计算
假如给出当前位置的经度、纬度和准确的天文时间,用什么方式计算出太阳当前的位置呢?例如算出类似于方位角和高度角的数据?
实际中使用希腊字符表示经纬度,不过为了输入方便,我改用XYZ表示.X是地方纬度,Y是太阳赤纬,t是时角,z=90度-h
高度h,方位角A,则
sinh=sinXsinY+cosXcosYcost
sinA=cosYsint/cosh
cosA=(sinXsinh-sinY)/(cosXcosh)
具体的可以在网上搜示意图
或这本书《天文学新概论》 苏宜 编著 华中理工大学出版社 2000年8月第一版 P52的章节第24章 太阳位置计算
[许剑伟 于家里 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°.76983*T + 0°.0003032*T^2 (Date平分点起算)
? 太阳平近点角: M = 357°.52910 + 35999°.05030*T - 0°.0001559*T^2 -0°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)
? 那么,太阳的真黄经是:Θ = Lo + C
? 真近点角是: v = M + C
? 日地距离的单位是天文单位,距离表达为:R = 1.000001018 (1-e^2) / (1+e*cos(v)) ……24.5式
? 式中的分子部分的值变化十分缓慢。它的值是:0.9997190 1800年0.9997204 1900年0.9997218 2000年0.9997232 2100年
? 太阳黄经Θ可由上述的方法算出,它是Date黄道分点坐标中的真几何黄经,需通过计算地心坐标星体位置也可算出。
? 要取得Date黄道坐标中太阳的视黄经λ,还应对Θ进行章动修正及光行差修正。如果精度要求不高,可用下式修正:
??? Ω = 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式)计算,式中ε是黄赤交角(由21章的21.2式计算)。
? tanα = cosεsinΘ / cosΘ ……24.6式? sinδ = sinεsinΘ ……24.7式
? 如果要想得到太阳的视赤经及赤纬,以上二式中的Θ应换为λ,ε应加上修正量:??? +0.00256*cos(Ω)
[译者注]:实际上就是对Θ补上黄经章动及光行差,ε补上交角章动后再转到赤道坐标中。也可在赤道坐标中补章动及光行差,但公式不同。
? 公式24.6当然可以转为:tan(α) = cos(ε)*tan(Θ),接下来,我们要注意α与Θ应在同一象限。然而,如果你使用计算机语中有ATN2函数(C语言是atan2),那最好保持24.6式不变,这样就可直接利用ATN2函数算出α,即:α = ATN2( cos(ε)*sin(Θ),cos(Θ) )例24.a——计算1992-10-13,0点,即力学时TD=JDE 2448908.5时刻的太阳位置。
我们算得:
T = -0.072183436Lo= -2318°.19281 = 201°.80719M = -2241°.00604 = 278°.99396e = 0.016711651C = -1°.89732Θ= 199°.90987 = 199°54 36R = 0.99
您可能关注的文档
- 天然药物化学(中药资源专业整理).pdf
- 天然药物化学习题集.pdf
- 天然气集输第一章.pptx
- 天王封魔版最新完整暴率.doc
- 天然气集输 第二章.ppt
- 天籁公爵.ppt
- 天线增益的计算.doc
- 天线安装.docx
- 天线考试总结.doc
- 天线基础.ppt
- [吉林]2023年吉林省直中小学“强师”招聘34人笔试历年高频考点(难、易错点)附带答案详解.docx
- [台州]浙江台州市黄岩区市场监督管理局招聘编外工作人员笔试历年高频考点(难、易错点)附带答案详解.docx
- [台州]浙江台州市建设工程质量检测中心招聘编制外用工4人笔试历年高频考点(难、易错点)附带答案详解.docx
- [台州]浙江台州椒江区大陈镇综合执法队招聘笔试历年高频考点(难、易错点)附带答案详解.docx
- [合肥]2023年安徽合肥市公安局第二批招聘警务辅助人员543人笔试历年高频考点(难、易错点)附带答案详解.docx
- [吉安]2023年江西吉安市大学生乡村医生专项招聘笔试历年高频考点(难、易错点)附带答案详解.docx
- [吉安]2024年江西吉安市中心人民医院招聘卫生专业技术人才73人笔试历年高频考点(难、易错点)附带答案详解.docx
- [合肥]2023年安徽医科大学管理岗和专业技术辅助岗招聘46人笔试历年高频考点(难、易错点)附带答案详解.docx
- [吉安]江西吉安市青原区人民医院人才引进笔试历年高频考点(难、易错点)附带答案详解.docx
- [台州]浙江台州温岭市交通运输局招聘编外人员笔试历年高频考点(难、易错点)附带答案详解.docx
文档评论(0)