卫星和飞船的跟踪测控-C题-朱

2009高教社杯全国大学生数学建模竞赛

承诺书

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A

我们的参赛报名号为（如果赛区设置报名号的话）：

所属学校（请填写完整的全名）：湖南城市学院

参赛队员 (打印并签名) ：1.

2.

3.

指导教师或指导教师组负责人 (打印并签名)：

日期：2009年 9 月日赛区评阅编号（由赛区组委会评阅前进行编号）：

2009高教社杯全国大学生数学建模竞赛编号专用页

赛区评阅编号（由赛区组委会评阅前进行编号）：

赛区评阅记录（可供赛区评阅时使用）：

全国统一编号（由赛区组委会送交全国前编号）：

全国评阅编号（由全国组委会评阅前进行编号）：

卫星和飞船的跟踪测控

摘要

对于问题一：我们根据已知问题一条件的假设和测控站的测控范围只考虑与地平面夹角3度以上的空域的条件，可以画出卫星的轨道图以及测控站的测控范围，同时测控站的测控范围可由其圆心角来表示，则测控站的数量可由360。与其圆心角的商得到。根据正弦定理，可得到轨道半径与测控范围的关系，继而得到测控站的数量个数，其个数至少为3个。

对于问题二：由问题二的假设，而且假设每个测控站的测控范围都相同

的条件下，我们想到把对卫星或飞船进行测控的测控站的个数问题看成在一个大的球壳（或部分球壳）上割去等大小的球盖问题，求的是割的最少次数。我们把球壳沿着一个平面方形切割成无数的圆环，再将圆环剪断成带状，按原来的顺序铺在一个平面上，中点在Y 轴上，即成一个椭圆形。再用问题一得出的弧长为半径的圆去覆盖。又因为在周长相同的几何图形中，圆的面积最大。然而如果我们用不相重叠的圆形来覆盖整个平面时，将会存在一些缝隙没被盖住。另一方面，在一个固定的圆面积中，使用规则六边形可以覆盖尽可能大的面积而不留缝隙。所以我们用等边六边形的网去覆盖椭圆即可得出所需的测控站最少个数，其最少个数应为012t

i i N N N ==+∑个。

对于问题三：由于问题三是求实际测控站的测控范围，由于各测控站的位置都不是很均匀，因而难以直接计算出来，于是我们把飞船的球带轨道近似成一个矩形，用经度和纬度建立坐标。把经度与纬度的乘积看做面积。那么在其上的测控站的测控区域为22πθ。再根据经纬度将各个站添加上去，求出每个站与其他站的重合测控区域。用总数乘以22πθ再减去重合区域即等到有效区域，再除以总面积可得测控率。

关键词： 测控站 跟踪测控 测控率

一、问题的提出

1.1背景

随着科学技术的提高，卫星和飞船的应用越来越广泛，同时卫星和飞船

在国民经济和国防建设中有着重要的作用。对它们的发射和运行过程进行测控是

航天系统的一个重要组成部分，理想的状况是对卫星和飞船（特别是载人飞船）进行全程跟踪测控，然而在现实生活中，我们以现在的技术还不能达到，所以在卫星或飞船的发射与运行过程中，往往采取以多个测控站联合来完成测控任务，因此，测控站的选取也决定了能否对卫星和飞船进行全程跟踪测控。

1.2问题重述

卫星和飞船在国民经济和国防建设中有着重要的作用，对它们的发射和运行过程进行测控是航天系统的一个重要组成部分，理想的状况是对卫星和飞船（特别是载人飞船）进行全程跟踪测控。

测控设备只能观测到所在点切平面以上的空域，且在与地平面夹角3度的范围内测控效果不好，实际上每个测控站的测控范围只考虑与地平面夹角3度以上的空域。在一个卫星或飞船的发射与运行过程中，往往有多个测控站联合完成测控任务，如神州七号飞船发射和运行过程中测控站的分布如下图所示：

因此我们需要建立模型去分析卫星或飞船的测控情况，需要解决的问题是：问题一：在所有测控站都与卫星或飞船的运行轨道共面的情况下至少应该建立多少个测控站才能对其进行全程跟踪测控？

问题二：如果一个卫星或飞船的运行轨道与地球赤道平面有固定的夹角，且在离地面高度为H的球面S上运行。考虑到地球自转时该卫星或飞船在运行过程中相继两圈的经度有一些差异，问至少应该建立多少个测控站才能对该卫星或飞船可能飞行的区域全部覆盖以达到全程跟踪测控的目的？

问题三：收集我国一个卫星或飞船的运行资料和发射时测控站点的分布信息，分析这些测控站点对该卫星所能测控的范围。

1.3问题分析

对于问题一：我们根据已知条件的假设，和测控站的测控范围只考虑与地平面夹角3度以上的空域的条件，可以画出卫星的轨道图以及测控站的测控范围，

同时测控站的测控范围可由其圆心角来表示，则测控站的数量可由360。与其圆心角的商得到。根据正弦定理，可得到轨道半径与测控范围的关系，继而得到测控站的数量个数。

对于问题二：由问题二的假设，而且假设每个测控站的测控范围都相同的条件下，我们想到把对卫星或飞船进行测控的测控站的个数问题看成在一个大的球壳（或部分球壳）上割去等大小的球盖问题，求的是割的最少次数。我们把球壳沿着一个平面方形切割成无数的圆环，再将圆环剪断成带状，按原来的顺序铺在一个平面上，中点在Y轴上，即成一个椭圆形。再用问题一得出的弧长为半径的圆去覆盖。又因为在周长相同的几何图形中，圆的面积最大。然而如果我们用不相重叠的圆形来覆盖整个平面时，将会存在一些缝隙没被盖住。另一方面，在一个固定的圆面积中，使用规则六边形可以覆盖尽可能大的面积而不留缝隙。所以我们用等边六边形的网去覆盖椭圆即可得出所需的测控站最少个数。

对于问题三：由于问题三是求实际测控站的测控范围，由于各测控站的位置都不是很均匀，因而难以直接计算出来，于是我们把飞船的球带轨道近似成一个矩形，用经度和纬度建立坐标。把经度与纬度的乘积看做面积。那么在其上的测

控站的测控区域为2

2πθ。再根据经纬度将各个站添加上去，求出每个站与其他站的重合测控区域。用总数乘以2

2πθ再减去重合区域即等到有效区域，再除以总面积可得测控率。

二、模型假设

（1）假设地球是一个理想的球形，不考虑其他影响建站的因素。

（2）假设每个测控站的测控范围都相同。

三、符号说明

符号说明：

H……………………………………为卫星到地球的表面高度。θ……………………………………为测控站到地点的圆心角的一半。

2θ……………………………………为测量圆心角。N………………………………………为测控站的数量。R………………………………………为地球的半径（且取平均半径6371km）。h………………………………………为卫星的轨道半径。

r………………………………………为正六边形的边长

i

l………………………………为测控站能测到s球面的圆弦长

s………………………………为第i个纬度所在平面分s球面所截平面的圆周长i

N……………………………为第i个纬度上的测控站的站点数

i

T………………………………为测控站的总数?………………………………为卫星或飞船与赤道平面的夹角n…………………………………相差多少个纬度

四、问题分析

五、模型的建立

模型一：根据题意可知，卫星的轨道如图所示：

所以我们根据正弦定理可得：

sin 93sin(87)

h R θ=-。。 模型二：此问题就相当于，在一个大的球壳（或球带）上割去等大小的球盖问题，求的是割的最少次数。于是我们把球壳沿着一个平面方形切割成无数的圆环，再将圆环剪断成带状，按原来的顺序铺在一个平面上，中点在Y 轴上，即成一个椭圆形。再用问题一得出的弧长为半径的圆去覆盖又因为在周长相同的几何图形中，圆的面积最大。然而如果我们用不相重叠的圆形来覆盖整个平面时，将会存在一些缝隙没被盖住。另一方面，在一个固定的圆面积中，使用规则六边形可以覆盖尽可能大的面积而不留缝隙。所以我们用等边六边形的网去覆盖椭圆，其图如下：

椭圆方程为

22222x h y ππ=-

由问题一可知:2,sin ,2l h r h s h θθπ===

∴ 我们由此可得出，在赤道上的测控站的站点数： 00023

s N l ≥? 于是，我们可得：相差θ角度即在θ纬度上时建立测控站

则 12(cos )s h πθ= 22(cos2)s h πθ=

32(cos3)s h πθ=

以此类推可得 2(cos )i s h i πθ=

所以即在θ纬度上建立测控站的站点数：11123

s N l ≥

? N ∈正整数 那么在i θ纬度上建立的测控站的站数为：23

i i i s N l ≥? N ∈正整数 离赤道最远的纬度应为卫星或飞船的轨道与赤道平面的夹角?，于是可得 ：

t θ

?= t ∈正整数 21t +即为分布测控站纬度的条数

则总的测控站最少个数为 012t

i i N

N N ==+∑ 模型三：

把飞船的球带轨道近似成一个矩形，用经度和纬度建立坐标。把经度与纬度的乘积看做面积。那么在其上的测控站的测控区域为22πθ。再根据经纬度将各个站添加上去，求出每个站与其他站的重合测控区域。用总数乘以22πθ再减去重合

区域即等到有效区域，再除以总面积可得测控率。

五、求解模型

模型一：

sin 93sin(87)

sin 93)sin 9387arcsin()h R R h

R h θθθ=--=∴-= 。。。

。。。由模型的建立，即sin(87

即θ=87。-sin 93arcsin()R h

。……………………………① 且h H R =+

3602o N θ≥…………………………………………………………②

且 N N +∈

∴我们就方程①进行讨论：

当h R 时，则θ取最大值为87。，所以我们由公式②，可得

2.07N ≥，又N N +∈，所以可得测控站的数量至少为3个。而且当测控站的数量为3个时，则θ为60o ，所以由方程①可得，h=2.2R , H=1.2R ，所以当H 1.2R ≥时，则测控站的数量随H 的增大，而增加 。

代入神七的数据，飞行高度为343km ，可计算得15.6o θ=

36011.542o

N θ

≥= 所以12N = 即应该设立12个测控站。

模型二：根据上述模型，我们带入神七的参数可得。

3t =

08N = 17N = 26N = 35N =

01244t

i i N N N ==+=∑

模型三：

我们通过网上收集了有关“神舟七号”载人飞船的运行资料和发射时测控站点的分布信息如下：

“神舟七号”载人飞船的测控站，除了五艘远望号测量船外，在国内，还有主场站、喀什站、和田站、东风站、青岛站、渭南站、厦门站等7个地面测控站；在国外，也有卡拉奇站（巴基斯坦）、马林迪站、圣地亚哥站（智利）、纳米比亚站等4个地面测控站（其分布图如下），这16个测控站组成了一个监控天网，即刻掌握着神舟七号的相关信息，并给神舟七号发送指令，全程保护神舟七号的安全。

下表数据为神舟七号飞船10个观测站所在经纬度:

该模型还未求解。

六、模型的评价及推广

模型二

还存在许多的缺陷，例如在我们将球面展开后进行覆盖，但是在复原成球形时，对接的两边已经被我们忽略掉了。还有是不是正六边形就能让测控站达到最少，正方形等我们未做考虑。

模型三

我们将球带看作矩形就存在误差，再以经纬度做为长度的单位，简化了计算的同时，也带来了较大的误差。

七、参考文献

八、附录