卫星轨道问题

1南京师范大学2015数学建模校内赛问题卫星轨道问题自1957年世界第一颗人造卫星发射升空以来，人造卫星在军事、气象预报、车辆导航、大地测量等方面有越来越多的应用。大部分卫星应用的前提是精确地确定卫星的轨道，根据牛顿定律，人造卫星绕地球沿椭圆轨道周而复始运行，但由于地球的非球形摄动力、日月引力、太阳辐射压、大气阻力等影响，卫星的实际轨道不会重复，确定卫星的精确位置并不容易。为描述卫星位置，需要建立适当的坐标系。本题惯性坐标系定义如下，取地球中心O为原点，地球自转轴取为z轴，指向北极为正向，x轴由O指向春分点，再按右手系确定y轴，建立直角坐标系。地心O在绕太阳椭圆轨道上运动，所以理论上该坐标系是非惯性系。但地球公转周期远大于卫星的观测弧段时长，故本题在短时间内认定该系为惯性坐标系，该坐标系不随地球旋转。由于受到岁差章动极移等影响，实际的春分点和地球自转轴都发生变化，国际大地测量学协会（IAG）和国际天文学联合会（IAU）决定，以2000年1月1日12时UTC为标准历元，定义协议天球坐标系，其z轴和x轴分别指向2000年1月1.日12时的瞬时平北天极和瞬时平春分点，这个坐标系称为J2000平赤道平春分点坐标系，将此坐标系记为ICRS，求解卫星运动方程时总是在该坐标系中进行。卫星的观测和应用是在地球上进行的，因此还需定义和地球固连的坐标系，和天球坐标系类似，国际上有一个协议地球坐标系，最新版本为ITRF2008，而美国国防制图局定义了一个地心地固坐标系WGS84坐标系，由于WGS84坐标系和ITRF相差很小，而GPS卫星使用WGS84坐标系，因此本题使用WGS84坐标系作为地球坐标系，WGS84直角坐标系和ICRS的转换公式见附录。以卫星观测站P为原点，P点的法线为z轴（指向天顶为正），以子午线方向为x轴（向北为正），y轴与x、z轴垂直（向东为正），可以建立站心地平直角坐标系以观测测站P为原点，至卫星s的距离r、卫星的方位角A、卫星的高度角h可以建立站心地平球坐标系。这些坐标系的转换关系见附录中的公式。请研究如下问题1.问题一已知卫星在时刻T0惯性坐标系ICRS中位置向量为（11000.0,39000.0,-4900.0）km,速度为（-1.5,1.0，-0.1）km/s计算卫星在该时刻的开普勒轨道根数，以及从T0时刻开始后一小时每次10分钟间隔的惯性系ICRS位置向量。2.问题二设一个观测站坐标北纬35.4251度，西经116.8893度，大地高为978.98,(单位km)测到某卫星的距离和角度2Datautc方位角A（度）俯仰角h(度)距离(km)2015/03/0208:00:00-165.837940.4421969.452015/03/0209：00:00-158.316170.665120514.34计算卫星在两次观测时刻的卫星在坐标系ICRS中的位置向量，用这两组位置向量计算卫星轨道根数（第一个观测时刻），分析岁差、章动、极移对计算卫星轨道的影响。3.问题3设有卫星轨道观测数据，数据格式如下日期时分秒wgs84xwgs84ywgs84z2015/3/205134273780.417-3405983.549-4581931.3010133649554.012-1759990.152-5867671.8第一列为日期，第二、三、四列分别为小时、分、秒（UTC时间），第五、六、七列分别为WGS84坐标系x,y,z坐标，数据文件为data1.xlsx。请估计卫星的实际轨道，给出卫星在1点到12点处整点处的坐标（ICRS坐标系和WGS84坐标系）4.问题4卫星数据同第三问。设有北纬32.05度，东经118.84度，大地高为200米的卫星观测站，确定卫星可由该观测站最长可观测时段（卫星高于当地地平线），如果仅仅只有时间连续一段可观测数据，如何外推卫星轨道？3附录卫星轨道根数轨道倾角i：inclinationangle赤道平面与卫星轨道平面间的夹角升交点卫星由南往北飞行轨迹在赤道平面上交点。升交点赤经Ω：therightascensionofascendingnode：春分点轴向东度量到升交点的角度轨道半长轴a：semi-majoraxis偏心率e：eccentricity(0≤e1)焦距与轨道半长轴之比。近地点幅角ω:argumentoftheperigee在轨道平面内升交点和近地点与地心连线间的夹角。平近点角M:meananomaly若卫星通过近地点的时刻为tp，卫星的平均角速度为n，则时刻t的平近点角M=n（t-tp）。坐标系天球的概念以地球质心M为球心，以任意长为半径的假想球体称为天球。天文学中常将天体沿天球半径方向投影到天球面上，再根据天球面上的参考点、线、面来确4定天体位置。（1）天轴与天极地球自转轴的延伸直线为天轴，天轴与天球面的交点称为天极，交点Pn为北天极，位于北极星附近，Ps为南天极。（2）天球赤道面与天球赤道通过地球质心M且垂直于天轴的平面称为天球赤道面，与地球赤道面重合。天球赤道面与天球面的交线称为天球赤道。（3）天球子午面与天球子午圈包含天轴的平面称为天球子午面，与地球子午面重合。天球子午面与天球面的交线为一大圆，称为天球子午圈。天球子午圈被天轴截成的两个半圆称为时圈。（4）黄道地球绕太阳公转的轨道面称为黄道面。黄道面与赤道面的夹角ε称为黄赤交角，约为23.5°。黄道面与天球面相交成的大圆叫黄道，也就是地球上的观测者见到的太阳在天球面上的运行轨道。由于地球自转，对于地面上的观测者来说，天球赤道面不动而黄道面每日绕天轴旋转一周。又由于地球绕太阳公转，直观上看，太阳在黄道上每日自西向东运行约1°，每年运行一周。（5）黄极通过天球中心且垂直于黄道面的直线与天球面的两个交点称为黄极，靠近北天极Pn的交点Πn称为北黄极，Πs称为南黄极。（6）春分点当太阳在黄道上从天球南半球向北半球运行时，黄道与天球赤道的交点称为春分点，也就是春分时刻太阳在天球上的位置，如图中的Γ。春分之前，春分点位于太阳以东。春分过后，春分点位于太阳以西。春分点与太阳之间的距离每日改变约1°。5天球坐标系常用的天球坐标系有天球空间直角坐标系和天球球面坐标系。天球空间直角坐标系的坐标原点位于地球质心。z轴指向北天极Pn，x轴指向春分点Γ，y轴垂直于xMz平面，与x轴和z轴构成右手坐标系。在天球空间直角坐标系中，任一天体的位置可用天体的三维坐标（x，y，z）表示。sinsincoscoscosrzyx天球球面坐标系的坐标原点也位于地球质心。天体所在天球子午面与春分点所在天球子午面之间的夹角称为天体的赤经，用α表示；天体到原点M的连线与天球赤道面之间的夹角称为赤纬，用δ表示；天体至原点的距离称为向径，用r表示。这样，天体的位置也可用三维坐标（α，δ，r）唯一地确定。天体的天球空间直角坐标系和球面坐标系是同一天体位置的不同表达方式。两种表达方式可通过下面的式子进行转换。22222yxzxyzyxrarctanarctan岁差与章动的影响地球绕自转轴旋转，在无外力矩作用时，其旋转轴指向应该不变。但由于日月对地球赤道隆起部分的引力作用，使得地球自转受到外力矩作用而发生旋转轴的进动现象，即从北天极上方观察时，北天极绕北黄极在圆形轨道上沿顺时针方向缓慢运动，致使春分点每年西移50.2″，25800年移动一周。这种现象叫岁差。在岁差影响下的北天极称为瞬时平北天极，相应的春分点称为瞬时平春分点。瞬时平北天极绕北黄极旋转的圆称为岁差圆。图5-3岁差与章动瞬时平北天极瞬时北天极岁差、章动叠加章动椭圆岁差圆Πn6协议天球坐标系由上可知，北天极和春分点是运动的，这样，在建立天球坐标系时，z轴和x轴的指向也会随之而运动，给天体位置的描述带来不便。为此，人们通常选择某一时刻作为标准历元，并将标准历元的瞬时北天极和真春分点作章动改正，得z轴和x轴的指向，这样建立的坐标系称为协议天球坐标系。国际大地测量学协会（IAG）和国际天文学联合会（IAU）决定，从1984年1月1日起，以2000年1月1.5日为标准历元。也就是说，目前使用的协议天球坐标系，其z轴和x轴分别指向2000年1月1.5日的瞬时平北天极和瞬时平春分点。将此坐标系记为ICRS.J2000.0为了便于区别，z轴和x轴分别指向某观测历元的瞬时平北天极和瞬时平春分点的天球坐标系称为平天球坐标系，z轴和x轴分别指向某观测历元的瞬时北天极和真春分点的天球坐标系称为瞬时天球坐标系。为了将协议天球坐标系的坐标转换为瞬时天球坐标系的坐标，须经过如下两个步骤的坐标转换。（1）将协议天球坐标系的坐标转换为瞬时平天球坐标系的坐标以TMtzyx,,和TCISzyx,,分别表示天体在瞬时平天球坐标系和协议天球坐标系中的坐标，因两坐标系原点同为地球质心，所以只要将协议天球坐标系的坐标轴旋转三次，便可转换为瞬时平天球坐标系的坐标，转换公式如下：zzyzMtICRSxxyRRRyzz式中ζ、θ、z为坐标系绕z轴和x轴旋转的角度，其值由观测历元与标准历元之间的时间差计算。=PzyzRzRR为岁差矩阵。（2）将瞬时平天球坐标系的坐标转换为瞬时天球坐标系的坐标以Ttzyx,,表示瞬时天球坐标系的坐标，则转换公式如下：MtxzxtzyxRRRzyx式中ε为观测历元的平黄赤交角；Δψ，Δε分别为黄经章动和交角章动。=NxzxRRR，称为章动矩阵协议地球坐标系地球坐标系确定卫星位置用天球坐标系比较方便，而确定地面点位则用地球坐标系比较方便。最常用的地球坐标系有两种，一种是地球空间直角坐标系，另一种是大地坐标系。地球空间直角坐标系的坐标原点位于地球质心（地心坐标系）或参考椭球中7心（参心坐标系），z轴指向地球北极，x轴指向起始子午面与地球赤道的交点，y轴垂直于xoz面并构成右手坐标系。大地坐标系是用大地经度L、大地纬度B和大地高H表示地面点位的。过地面点P的子午面与起始子午面间的夹角叫P点的大地经度。由起始子午面起算，向东为正，叫东经（0°~180°），向西为负，叫西经（0°~-180°）。过P点的椭球法线与赤道面的夹角叫P点的大地纬度。由赤道面起算，向北为正，叫北纬（0°~90°），向南为负，叫南纬（0°~-90°）。从地面点P沿椭球法线到椭球面的距离叫大地高。同一地面点在地球空间直角坐标系中的坐标和在大地坐标系中的坐标可用如下两组公式转换。BHeNzLBHNyLBHNxsinsincoscoscos2122221eNBzHyxBNezBxyLsinsinarctanarctan式中e——子午椭圆第一偏心率，可由长短半径按式2222abae/算得，BeaN221sin/算得。公式中的B必须用迭代的方法求解。Matlab中有函数作变换，将经纬度转为直角坐标[x,y,z]=geodetic2ecef(phi,lambda,h,ellipsoid)其中phi为纬度，lambda为经度，h为大地高，ellipsoid=[a,e]，a为长半轴，e为偏心率将直角坐标转为经纬度使用函数[phi,lambda,h]=ecef2geodetic(x,y,z,ellipsoid)注意matlab使用的是弧度地极移动与协议地球坐标系由于地球不是刚体，在地幔对流以及其它物质迁移的影响下，地球自转轴相对于地球体发生移动，这种现象叫地极移动，简称极移。在建立地球坐标系时，如果使z轴指向某一观测时刻的地球北极，这样的地球坐标系称为瞬时地球坐标系。显然，瞬时地球坐标系并未与地球固连，因而，地面点在瞬时地球坐标系中的位置也是变化的。为了比较简明地描述地极移动规律，国际纬度局根据1900.0至1905.0期