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