1、1卫星和飞船的跟踪测控摘 要卫星和飞船在国民经济和国防建设中有着重要的作用，本文通过对卫星或飞船运行过程中测控站需要的数目进行求解，从而实现能够对卫星或飞船进行全程跟踪测控的目标。对于问题一，由于测控站都与卫星运行轨道共面，且测控站的测控范围只考虑与地平面夹角 3 度以上的空域，所以，我们首先考虑将卫星或者飞船的运行轨道理想化成圆形，建立其与地球共心的圆形轨道模型，此时，运用几何知识和正弦定理计算出至少应建立 12 个测控站。但是，在现实中卫星或飞船的轨道为椭圆形状，接着我们又给出了质点运行轨道为椭圆时的数学模型计算得出需要建立测控站数目的区间为 12至 16 个。问题二，我们利用每个测控站测
2、控的锥形区域与卫星或飞船轨道曲面相交的圆的内接多边形来覆盖整个卫星轨道曲面，就可以将需要这样内接多边形的个数近似的看作需要建立测控站的最少个数，这里我们只给出内接正四边形和正六边形两种数学模型，此时，计算出需要测控站的最少数目分别为 60 和 67 个。问题三，通过网络查询得到神舟七号的观测站位置和数目，以及飞船运行的倾角和高度等相关数据。通过线性拟合我们发现测控站的位置近似符合正弦曲线。最后，我们给出了模型优缺点的分析和评价，并提出了模型的改进的方向。关键字：卫星或飞船的跟踪测控；圆形轨道模型；圆锥测控模型；测控站点的数目21、问题重述1.1 背景资料现代航天工业中卫星和飞船的测控设备只能观
3、测到所在点切平面以上的空域，且在与地平面夹角 3 度的范围内测控效果不好，实际上每个测控站的测控范围只考虑与地平面夹角 3 度以上的空域。在一个卫星或飞船的发射与运行过程中，往往有多个测控站联合完成测控任务，因此需要分析卫星或飞船的测控情况。 1.2 需要解决的问题问题一：在所有测控站都与卫星或飞船的运行轨道共面的情况下至少应该建立多少个测控站才能对其进行全程跟踪测控。问题二：如果一个卫星或飞船的运行轨道与地球赤道平面有固定的夹角，且在离地面高度为 H 的球面 S 上运行。考虑到地球自转时该卫星或飞船在运行过程中相继两圈的经度有一些差异，问至少应该建立多少个测控站才能对该卫星或飞船可能飞行的区
4、域全部覆盖以达到全程跟踪测控的目的。问题三：收集我国一个卫星或飞船的运行资料和发射时测控站点的分布信息，分析这些测控站点对该卫星所能测控的范围。2、问题分析问题一主要是求在测控站与卫星或飞船的运行轨道共面的情况下至少应该建立多少个测控站才能对其进行全程跟踪测控的问题，我们首先应该考虑求出每个测控站能够测控的有效弧度，然后利用卫星或飞船运行的轨道周长除以每个测控站测控的有效弧度就可得到需要的测控站数目。图 1 卫星或飞船运行轨道示意图3问题二，卫星或飞船的运行轨道与地球赤道平面有固定的夹角，在离地面高度为H 的球面 S 上运行且地球自转时卫星或飞船在运行过程中相继两圈的经度有一些差异，这是一个求
5、球面积的问题，建立模型一将问题简化为求卫星或飞船的运行轨道面积，假设其运行的轨道是一个球面，计算出每个测控站在这个球面上的有效投影，然后求出需要的最少测控站个数。模型二在地球自转的影响下，卫星运行过程中星下点轨迹在地球表面形成一些“8”字型的轨迹，对称地分布在赤道两边。图 2 卫星或飞船对于地球和太阳的位置考虑到卫星星下点轨迹的密集程度问题，对于星下点轨迹圈数较少的卫星，可根据其相应的轨迹进行分析求解;由于测控站在卫星轨道面的测控区域是个圆形，因此对于星下点轨迹圈数较多较密的卫星，利用测控区域的圆内接正方形来覆盖整个卫星轨道面，从而计算出全程监控所需的最少的测控站的个数问题三，在网上搜集我国一
6、个卫星或飞船的运行资料和发射时测控站点的分布信息，从而进行分析这些测控站对该卫星所能测控的范围。3、问题假设1、假设卫星或飞船相对于地球极小可以看做质点2、假设地球是个规则球体，质量集中于地心3、假设外界引力对该系统可忽略不计4、忽略影响测控站布置的地理因素5、不考虑测控站周围地理环境和天气环境对卫星测控的影响44、符号的约定R 地球半径H 卫星或飞船距地面高度近地点高度1远地点高度2测控站测控范围与卫星运行轨道曲面相交的半径r第i行正方形覆盖的轨道面圆周长iC,3.=圆内接正方形每条边在卫星轨道面上所对应的圆弧长l轨道与赤道平面夹角q球帽面积1S卫星运行曲面的面积2测控站测控的范围与曲面交线
7、圆的面积3圆内接正六边形的面积4圆内接正六边形与圆的面积之比e5、模型的建立及求解5.1 问题一模型的建立及求解5.1.1 模型一：假设卫星或飞船运动轨道为圆在不考虑地球自转的条件下，地球自转时该卫星或飞船在运行过程中相继两圈的经度的差异可不予考虑。卫星或飞船从起飞时加速升空后经一系列加速变轨，最终的运行轨迹是一圆周。即最终卫星或飞船绕地球做匀速圆周运动。用卫星或飞船的运动轨迹所在的平面去切地球会得到一圆面。如图 l 所示：图 3 观测站对圆形卫星轨道覆盖范围示意图5RHCOAB12地球卫星轨道我们只需在如图 C 点建立一测控站即可测控 A 至 B 之间的劣弧区域，最小测控站数目即为需要覆盖卫
8、星轨道的这样的 C 点的个数，利用正弦定理解三角形 13sin(90)sinRHRO81B362nC按照此模型以神州七号飞船为例：地球半径为 6400 公里，飞船进入预定轨道运行稳定后距地球表面高度为 343 公里，相关数据代入，运用 MATLAB 计算得出，n=12，71.408,15.92OBCOB 即此时需要最少测控站的数目为 12 个。5.1.2 模型二：考虑到实际，按卫星或飞船运动轨道为椭圆由于在实际情况中飞船的运行轨道为椭圆形，如图 2 或下图，取椭圆近地点旁边的焦点为地球的圆心，椭圆轨道定位很麻烦，因此先估算，然后再精算61、以地心为圆心，地球半径与近地点之和为半径作圆，如图 4
9、、由于圆包含在椭圆区域之内，若能监控到圆周及以外空域，则定能监控到椭圆及以外空域，因此，在地球上均匀建站监控整个圆周。图 4 观测站对椭圆形卫星轨道覆盖范围示意图 1Rq2f1地心具体算法为： sin93iqRH12f7193sin93()arc80Rf H 1nf其中 是如图所示的圆心角， 角如图所示，以神舟七号为例，近地点高度 =20012,f q 1H公里，所以 ，用 MATLAB 软件解得 n=16 。1H2、以地心为圆心，地球半径与远地点之和为半径作圆，如图 5 所示。图 5 观测站对椭圆形卫星轨道覆盖范围示意图 2RH2f1焦点 地心BA由于大圆包含了椭圆区域，因此只要监控到大圆周
10、及以外空域，则未必能监控整个椭圆周。在地球上均匀建站监控整个圆周，其算法和 1）中相同： ，以神舟2HS七号为例，其远地点 公里，令 ，解得 n=12。2347H2H综上，椭圆轨道上的监控站应该在 12 至 16 个之间。下面我们进行椭圆轨道的监控站数精确计算。1、椭圆轨道的测控算法思想用逼近方法和迭代算法来实现。由于对椭圆监控不能像对待圆一样均匀分布站点，而 ，因此分布测控站是一件很棘手的事情。 可用前面模型中的算法， 可借12f 1f 2f用 的算法来算出 的近似值。根据物理学、力学知识及开普勒三大定律可得到 右2f侧的向径 ，a 为椭圆的长半轴，e 为离心率，f 为向径与椭圆水平长轴的(
11、)cosrf夹角，r 为球心到椭圆上的向径。所以椭圆轨道的监控算法为：8，12HRace1S0rR1 193sin93()ac8f r21(cosi kerf用逼近的方法可近似计算出 。14n5.2 对于问题二模型的建立与求解5.2.1、模型一：由于一个卫星或飞船的轨道与地球赤道平面有固定夹角，而求地球自转时该卫星或飞船在运行过程中相继两圈的经度有一些差异，因此，对卫星在地球表面的星下点轨迹进行分析。1、当卫星运行角度与地球自转角速度相同时，卫星沿运行轨道运行一周后星下点轨迹又回到起点，星下点轨迹见图 5，但其相继两圈的经度无变化，不合题意。2、当卫星运行角速度是地球自转 速度的两倍是卫星沿运
12、行古道运行两圈后星下点轨迹回到起点。3、当卫星运行角速度是地球自转速度的三倍时，卫星沿运行轨道运行三圈星下点轨迹回到起点。图 6 卫星或飞船地面轨迹卫 卫 卫 卫 卫 卫 卫 卫 卫0 60 120 180 240 300 360-90-60-300306090卫星运行角速度越大，卫星在地球表面上留下的星下点轨迹越密，从而形成了卫星或飞船地面轨迹，如图 6 所示。9由以上图形观察分析可以得出：星下点轨迹均匀地分布在赤道的两边，即北纬 a与南纬 a 之间，因此，在卫星或飞船的星下点轨迹较为简单的情况下可沿着星下点轨迹设立测控站，对于一般情况，特别是卫星或飞船的星下点较密时，只需测控北纬 a与南纬
13、 a 之间的区域，就可以实现对该卫星的全程监控。在实际测控范围与卫星轨道面的相交圆 D 中取圆内接正方形并用内接正方形覆盖所要测控的区域。首先根据轨道面与赤道面的夹角 a 以及每个内接正方形边长在地心所对的圆心角 ，可求出测控区域中正方形的行数 j，如图所示；其次根据每一行正2方形覆盖的轨道面的圆周长即可确定每一行所需的正方形的个数 i，从而计算出全程监控所需要的监控站的最少个数，如图 7 所示。图 7 圆内接正方形覆盖法Dr OCD模型求解过程：1、计算正方形的边长在地心所对的圆心角为 ，有2ACrtan2DO()sirRHcohin93arctnt(87ar)2Rh由于卫星围绕地球运动时，
14、轨道与赤道平面有一个夹角 a，测控范围由正方形覆盖，所以正方形的行数： 2j圆内接正方形每条边在卫星轨道面上所对应的圆弧长为： ()()lRhh102、当 j 为偶数时，各行正方形在南北纬之间对称分布，每一行正方形的总长度是正方形所在轨道面的圆的周长，北纬第一行正方形覆盖的轨道面圆周长，12(rcosCR）北纬第一行排列的正方形个数 ，第二行正方形覆盖的轨道面圆周长：nl=，第二排排列的正方形个数 ，依次类推：北纬第 k 行排2()cos(3CRrpb=+ 2Cnl=列的正方形个数 ，为实现全程测控需要的测控站总数kCnl= 21jkn=当 j 为奇数时，先用一行正方形覆盖赤道，使正方形的中心
15、在赤道上，然后将各行正方形在南北半球对称分布，所以赤道上排列的正方形的个数为 ，0()RHlp+北纬第一行正方形覆盖的轨道面圆周长为 ，北纬第一行排列的正12()cos2CRrpb=+方形个数 ，同理可以求出第 k 行上分布的测控站的个数 ，实现全程测1Cnl= kCnl=控的测控站的总数120jkn-=+根据神舟七号的运行数据，将轨道与赤道面夹角 42.2 度， 公里代入a34H得 j=4， n=60，即如果要实现对神舟七号的全程监控至少需要 60 个测控站对其进行测控。5.2.2、问题二中，为使计算方便，采用内接正方形覆盖所要测控的区域.，而圆内接正六边形的面积占圆面积的 82.74%，因
16、此可以考虑采用圆内接正六边形覆盖，以提高有效覆盖率，我们因此建立了模型二：设卫星或飞船绕地球飞行的倾角为 ，离地面的高度为 ,地球半径为 R，卫星qH或飞船飞行的轨迹为一环绕地球半径为 的球面去掉两端球帽的曲面，通过公式RH+计算出球帽的面积得到卫星或者飞船飞行的曲面面积，每个测控站测控的范围与曲面的交线为一个半径为 r 的圆，只有圆与圆之间有重叠才能够保证测控站的测控范围覆盖所有曲面，这样我们就取每个圆的内接正六边形作为每个测控站对卫星或飞船运行曲面的平均有效测控范围，通过圆内接六边形与圆的面积比率计算出圆内接正六边形的面积，用曲面的面积去除以正六边形的面积，就可以得到需要的最少测控站数目，具体计算方法如下：球帽面积公式： 12()1sin)SRHpq=+-球面面积计算公式： 24卫星运行曲面的面积： 21S-
