数学建模-公式编辑和论文格式.doc

个人收集整理 勿做商业用途
1、
2、
3、
4、
5、

6、
 300

100
400

200

300

600

500
200
400
300
600

700
500

7、
EXCEL 制 表 题
    第1题［共一题]:
    请点击打开指定的文件,将sheet1中的数据复制到Sheet2中，对sheet2作如下操作(sheet1中的数据已经保护,不可以直接操作；并请在修改的过程中随时保存您所做的修改)：
1、将标题“职工工资表”设置为通栏、居中、红色、18号字.
2、第二行小标题设置为水平和垂直居中、绿色、16号字。
3、计算每位职工的实发工资及全体职工的基本工资和实发工资的合计值。
4、按实发工资从高到低排序.
5、为表格加边框线,外边框为粗线,内边框为细线.
6、按姓名和实发工资作“三维簇状条形图”，图形放在表格下方。

8、请将建筑1021班成绩按名次从高到低排名,总成绩=每门分数*学分和/总学分

2009高教社杯全国大学生数学建模竞赛C题
卫星和飞船的跟踪测控
卫星和飞船在国民经济和国防建设中有着重要的作用，对它们的发射和运行过程进行测控是航天系统的一个重要组成部分，理想的状况是对卫星和飞船(特别是载人飞船)进行全程跟踪测控。
测控设备只能观测到所在点切平面以上的空域,且在与地平面夹角3度的范围内测控效果不好，实际上每个测控站的测控范围只考虑与地平面夹角3度以上的空域。在一个卫星或飞船的发射与运行过程中，往往有多个测控站联合完成测控任务，如神州七号飞船发射和运行过程中测控站的分布如下图所示：

图片来源 http://www.gov。cn/jrzg/2008-09/24/content_1104882。htm
请利用模型分析卫星或飞船的测控情况,具体问题如下：
1. 在所有测控站都与卫星或飞船的运行轨道共面的情况下至少应该建立多少个测控站才能对其进行全程跟踪测控？
2。如果一个卫星或飞船的运行轨道与地球赤道平面有固定的夹角,且在离地面高度为H的球面S上运行。考虑到地球自转时该卫星或飞船在运行过程中相继两圈的经度有一些差异，问至少应该建立多少个测控站才能对该卫星或飞船可能飞行的区域全部覆盖以达到全程跟踪测控的目的？
3。 收集我国一个卫星或飞船的运行资料和发射时测控站点的分布信息，分析这些测控站点对该卫星所能测控的范围.

卫星的跟踪测控问题简析
摘要:本文主要利用计算几何与图论的有关知识，分析和解决了卫星的跟踪测控问题,并应用仿真模拟手段对地球自转、非自转、不同轨道、不同纬度的卫星跟踪测控进行了较为详细的讨论。针对问题一、问题二给出了求解过程和结果，并提出了优化改进模型，针对第三问通过神舟七号卫星参数和测控站点分布数据分析了测控站点对该卫星所能测控的范围。
(1）问题一考虑到所有测控站点与卫星的运行轨道共面，测控点个数只与卫星轨道的高度相关，将问题一转化为考虑卫星不同轨道高度、地球无自转情况下，测控站点最少测控模型求解问题。从而得到：低轨道卫星(小于500km）至少需要10个测控站，中轨道卫星（500~2000km）至少需要5个测控站，高轨道卫星(2000~20000km）至少需要3个测控站,太阳同步卫星（700~1000km）至少需要7个测控站，地球同步卫星轨道高度（35786km)远大于7651。5km，至少需要3个测控站。最后对模型1进行了优化，增加了卫星发射过程中卫星测控站点个数确定。
(2)问题二在地球自转情况下，本文首先对问题进行细化,基于不同纬度、不同轨道高度的情况分别考虑，然后利用蜂窝理论将问题转化为求解测控点有效测控面积覆盖问题,得到各种情况下所需测控站点个数。轨道倾角为0—30°时，低轨道23个，中轨道7个，高轨道3个；轨道倾角为30°-60°时，低轨道39个，中轨道12个，高轨道5个；轨道倾角为60°—90°时，低轨道45个,中轨道14个，高轨道6个。
最后针对模型2计算所得的测控站个数可能会出现测控盲区情形，利用有效等六边形无缝拼接方法，提出了进一步优化方案,得出了相对于模型2更为合理的站点个数。
（3）问题三本文通过获取神州七号卫星运行资料和测控点的分布信息，首先根问题二中得到的模型给出了全程测控神舟七号卫星的测控站点个数:28个，说明了神舟七号卫星16个站点并不能够完全测控;并由测控站点位置信息，利用蜂窝理论中等六边形有效面积，给出了测控站点有效测控范围平面图。

关键词：卫星全程跟踪测控 蜂窝理论 有效测控范围 计算几何 仿真模拟

一、问题重述
本题研究的是对卫星和飞船进行全程跟踪测控问题。其中每个测控站的测控范围只考虑与地平面夹角3度以上的空域，题目要求完成以下三个问题：
1、在只考虑测控站都与卫星或飞船的运行轨道共面的情况下，确定能进行全程跟踪卫星的测控站点最少数目。
2、考虑地球自转，同时卫星或飞船运行轨道与地球赤道平面有固定的夹角，确定能进行全程跟踪卫星的测控站点最少数目。
3、收集我国一个卫星或飞船的运行资料和测控站点的分布信息，利用问题1和问题2确定的模型分析这些站点对该卫星所能测控的范围。
二、问题1
1。问题假设
1) 假设每个测控站都正常工作；
2) 假设卫星或飞船的运行轨道为圆形;
3) 假设卫星或飞船不出故障；
4) 假设不考虑外太空各种空气阻力等因素；
5) 假设不考虑发射过程的监测；
6) 假设地球半径为6378km;
7) 假设卫星的发射是通过光电经纬仪光轴与平台棱镜精确对准；
8） 假设卫星发射时是东南方向.
2．符号说明
符号
符号说明

地球赤道半径

卫星或飞船距地垂直高度
(3)
测控站个数

测控站距卫星或飞船的最远距离

卫星或飞船的速度

卫星或飞船从发射到进入轨道的时间

卫星或飞船经过一个测控站所花的时间

一个测控站所监测卫星或飞船的最大运行弧长

卫星或飞船轨道周长
3。问题分析
根据问题1中所有测控站都与卫星或飞船的运行轨道共面的信息，则问题一转化为考虑不同卫星飞行轨道高度、地球无自转下,测控点最少分布问题。
根据三角形内切圆的性质以及测控设备可以监控与地面夹角3度以上的区域（如图1所示）,由此通过计算可知:当卫星或飞船运行高度高于7651。5千米的情况下,最少利用三个测控站便可完全监控卫星或飞船的飞行情况.
当卫星或飞船运行低于7651。5千米时，三个测控站便出现了监测盲点（不能完全监测卫星或飞船的运行情况），需要增加观测站个数（如图1所示）.
本问题通过对不同轨道高度的卫星进行划分，逐个分析其运行测控点个数。

图1
4。模型建立与求解
因为对于不同轨道高度卫星来说，地球同步卫星高度较高（35786km）,且位于赤道正上空，在赤道上设立3个测控站便可全程跟踪监测。
其余不同轨道卫星卫星类比于地球同步卫星，测控站分布在卫星运行轨道面与地球赤道面相交的交线上。
根据问题1的分析，当测控点与卫星运行轨道共面的情况下，测控点个数只与卫星轨道的高度相关，因此本文利用卫星的最低轨道来进行最少卫星个数进行确定（如图2所示）.
首先确定在不同测控站个数情况下，卫星或飞船的最低高度模型1（如式2。1、2.2所示）。
2。1
则可得：
2.2

图 2

分析式2。2可以发现高度与测控点个数负相关，计算可得不同卫星轨道下测控站点个数（见表1）。
表1不同测控站点个数下卫星轨道的最低高度
测控站的个数
卫星或飞船的最低高度（h，单位km）
3
…
3
7651．5
4
3140．7
5
1817．7
6
1216．5
7
883．34
8
678．68
9
541．3
10
444．4
11
374．66
12
319．04
13
276．84
14
243．22
15
215．94
对表1分析可知，当卫星或飞船运行的高度越低，需要的观测点越多。根据不同轨道卫星的高度,可知需要的最少测控站个数。
按照不同高度，卫星可以划分为以下几种类型，如表2所示。
表2 不同类型卫星高度表
按照高度划分的卫星类型
高度（Km）
地球同步轨道卫星
35786km
太阳同步卫星
700～1000km
高轨道卫星
2000～20000km
中轨道卫星
500～2000km
低轨道卫星
小于500km
由上述模型以及数据可知:低轨道卫星至少需要10个测控站，中轨道卫星至少需要5个测控站,高轨道卫星至少需要3个测控站，太阳同步卫星至少需要7个测控站，地球同步卫星轨道高度远大于7651.5km，至少需要3个测控站;
5。模型的优化
对于问题1的上述模型，并没有考虑卫星发射过程，这在一定程度会影响测控点的个数,本小节对模型1进行一定的优化。
由图3卫星发射过程示意图可得，在卫星进入轨道之前，有一个斜向加速和惯性加速过程，在此过程中，卫星存在一个相对于发射点的距离偏差。若该距离大于位于发射点的测控点的最大测控距离，则需要进行测控站点的调整。

图3 卫星发射过程示意图
根据图2中的测控站点与测控卫星的几何信息，可得到测控站点距卫星最远距离x与卫星高度、地球半径的关系式（如下式2。3）：
2.3
解得测控站点距卫星最远距离x为:

通过几何运算,可得到一个测控站所监测卫星或飞船的最大运行弧长:

经上计算得卫星或飞船经过一个测控站所花的时间T：

通过卫星或飞船经过一个测控站所花的时间T与卫星或飞船从发射到进入轨道的时间t的比较
当Tt时，卫星不会进入盲区；
当T接近或小于t时,卫星会进入盲区，需要增加测控站
本文这里采用神舟七号卫星数据,对提出的优化模型进行验证。
,
神七的入轨飞行速度7820。185米/秒，神七上升段飞行时间583.828秒，
可知：

=秒　　
如果按照上面的分析则神七在发射需到进入轨道需要再增加测控站才能进行全程测控。
而实际上运载火箭发射航天飞行器的飞行轨道方案是一次主动段就直接入轨(消耗能量较多）。
例如：神州七号飞船先被送入近地点200公里到远地点347公里的椭圆轨道，飞船在椭圆轨道飞行第1至5圈，由于大气阻力的影响,每圈轨道降低近1公里,飞船远地点高度从347公里降为343公里，第5圈,飞船远地点点火变轨,抬升近地点，轨道由椭圆轨道变成高度343公里圆轨道。
三、问题2
1.问题假设
1) 假设每个测控站都正常工作;
2) 假设卫星或飞船的运行轨道为圆形；
3) 假设卫星或飞船不出故障；
4) 假设卫星在空中的时间足够长;
5) 假设卫星在运行过程之中没有太阳风暴所引起的电磁干扰的影响；
6) 假设卫星卫星在发射时气象条件（温度、湿度、平均风速和高空风速)合适；
2.符号设定
符号
符号说明

地球赤道半径

卫星或飞船距地高度
(3）
测控站个数

测控站距卫星或飞船的最远距离

卫星或飞船在空中时间

卫星或飞船的轨道与赤道的夹角（轨道倾角)

一个测控站在空中形成的圆形测控区域半径

一个测控站在空中形成的有效测控区域面积

球缺曲面面积

所需测控站覆盖的面积
3。问题分析
根据对问题二的分析，由于地球的自转，当卫星在空中时间足够长的情况下，卫星在地球上的投影将会形成一个区域（见图4），而测控的目标是空中的卫星，从而考虑将投影区域向空中移动（移动距离为卫星高度h）形成一个“卫星”区域,为了达到测控的目的，测控的范围将要完全覆盖这一“卫星”区域。

时叫赤道轨道，卫星轨道在赤道赤道平面内，地球同步轨道卫星的轨道属于此类.
叫顺行轨道，多数卫星采用这种轨道，因为它可以利用地球自转速度，从而节省发射需要的能量，而且可以覆盖任务规定的区域,轨道倾角大则覆盖轨道区域大.
叫极轨道，在极轨道上的卫星可以观测整个地球，因此地球资源卫星、全球侦察卫星采用这种轨道。
轨道倾角的范围在内，其中地球同步卫星的倾角为，极地卫星为。根据倾角的范围将区域分别为(见表3）：
表3 根据轨道的倾角度数划分为不同的测控区域表
轨道倾角度数
测控区域

地球同步卫星区域

低纬度区域

中纬度区域

高纬度区域

极地卫星区域
每个测控站的测控范围为一个圆形区域，故我们采用圆形的蜂窝状分布（见示意图2，三个圆的圆心连线为正三角形，将每个圆形测控范围的有效范围为以R为边长的正六边形）,以达到测控站最少的目的。

图 6

4.模型建立与求解
1．地球同步卫星(轨道倾角）
地球同步卫星的倾角为，在赤道的上空,不受地球自转的影响，根据问题一的分析可知,地球同步卫星的最少测控站个数为3个。
2．根据表2、表3中对卫星运行状态的划分，本文首先按照纬度进行划分，再分别考虑不同的运行轨道参数。
经以上分析（如图7所示），可计算得测控站距卫星或飞船的最远距离x：
通过几何运算，可得：

所以：；

3。1
由此,计算可得每个测控站测控的正六边形有效范围，既:。

图 7

图 8
当地球自转时，轨道倾角卫星或飞船在地球上的投影汇聚成一个去掉两端的球缺（如图9所示），球缺的曲面面积为）;测控的目标是空中的卫星，从而考虑将投影区域向空中移动（移动距离为卫星高度h）形成一个“卫星区域"，此时球缺曲面面积为.则所需测控站覆盖的面积(s），是以r+h为半径的圆表面积减去去两个球缺的面积,既：.
测控站的个数为：
3。2
由问题1可知，.

示意图 9
将卫星或飞船的轨道倾角和高度代入上式3。1、3.2便可得到不同角度、不同高度的卫星或飞船所需测控站的个数。
表3是根据轨道倾角为、、,卫星高度分别为200km、500km、2000km、20000km时，计算的测控站个数。
表3 轨道倾角为、、各轨道高度所需测控站个数
轨道倾角
卫星高度(h/km）
测控站个数（Dn）
0～30°
500
23
2000
7
20000
3
30°~60°
500
39
2000
12
20000
5
60°～90°
500
45
2000
14
20000
6
5.模型优化
卫星轨道倾角在和时，由上面的模型计算所得的测控站个数可能会出现测控盲区。模型没有考虑卫星区域有可能不被正六边形所形成的无缝区域完全覆盖而是留有一定的盲区（图10)。

图10
当时，用无缝拼接的方法(图11）以最优分配测控站的个数。

图11
当时，用同样方法（图12)最优分布测控站的个数。

图12
此时对无缝拼接进行优化，沿水平垂直两方向，每个正六边形在覆盖区的有效面积乘以个数等于对应“卫星区域"在平面上区域的长和宽。
水平方向测控站的个数,，行数为奇数时取整，行数为偶数时取整后加上1.
垂直方向测控站的个数，

通过分析计算，可计算出轨道倾角30°、60°时，不同高度卫星或飞船所需测控站个数（见表4）。
表4 轨道倾角为30°、60°时各轨道高度所需测控站个数表
轨道倾角
卫星高度(h/km）
测控站个数（Dn）
0~30°
500
30
2000
13
20000
6
30°~60°
500
52
2000
16
20000
16
60°～90°
500
62
2000
19
20000
9
四、问题3
实验卫星采用神舟七号卫星。
1。资料收集
采用国家航天局提供的神舟七号卫星运行轨道参数以及测控点分布信息，具体数据见表5、6。

表5 神舟七号卫星运行轨道参数
运行周期
90min/圈
飞船高度
343km
轨道倾角
42。4°
运行圈数
45圈
表6 测控点位置信息（测控点坐标）
序号
测控站点
纬度
经度
1
主场站（北京站)
N
E
2
喀什站
N
E
3
和田站
N
E
4
东风站
N
E
5
青岛站
N
E
6
渭南站
N
E
7
厦门站
N
E
8
纳米比亚站
S
E
9
马林迪站
S
E
10
卡拉奇站
N
E
11
圣地亚哥
N
W
12
远一
远洋测量船随时间移动，位置信息不确定，这里给出的是近似平面位置信息。
13
远二
14
远三
15
远四
16
远五
表6中16个测控点数据分布平面图如图14所示。

图14 测控点数据分布平面图
2.模型建立与求解
2.1根据问题2的模型讨论神七的测控范围
由于地球每转一周，卫星在地球上的投影呈现如正弦曲线的分布形状，而当地球转动第二周时，卫星在地球上的投影与第一周形成的曲线,在同纬度有一个向右的偏移量（如图15所示），偏移量可以通过地球自转速度与卫星运行一周时间求的,则：

神七飞船在空中一共运行了45圈，飞船向右总的偏移量为为:

通过计算总偏移量与地球赤道周长的比值，
当时,则表明神州七号卫星在地球上的投影轨迹并没有形成一个完整周期；
当时，则表明卫星在地球上的投影轨迹并形成一个或者多个完整周期.
代入数据可得.说明飞船在地球上的投影轨迹有近3个完整的周期，基本上覆盖了倾角42。4°范围内的区域。
根据问题二中的模型,当卫星投影完全覆盖这一区域时,所需要的测控站个数为：45.而神七的测控站个数为16个,所以说并不能全程对神七进行跟踪测控。

图15 卫星星下点轨道平面图
2。2由测控站点求有效测控范围
根据表6测控点位置信息（测控点坐标）数据，本文采用蜂窝理论中的等六边形表示测控站点的有效测控范围，神七测控站点有效范围示意图如下图16所示.

图16 测控站点有效范围示意图
五、参考文献
［1］ 金永德，导弹与航天技术概论，哈尔滨，哈尔滨工业大学出版社,2002，P163—165。
[2] 褚桂柏，航天技术概论,北京，中国宇航出版社，2002,P23-27。
[3］ http://www。cnsa。gov。cn/n615708/index。html国家航天局。
［4] http://www.aerochina。
[5] http://www.holland—aerospace。com/荷兰国家航空航天协会网站.
［6］ 胡震宇，万慧琳，赵毅,2006年全国研究生数学建模竞赛A题优秀论文，南京师范大学.
［7］ 王庚，王敏生,现代数学建模方法,北京，科学出版社，2008。