GPS卫星运动及定位matlab仿真_毕业设计论文

毕业设计论文 GPS卫星运动及定位matlab仿真 摘要 全球定位系统是具有全球性、全能性、全天候优势的导航定位、定时和测速系统，现在在全球很多领域获得了应用。 GPS卫星的定位是一个比较复杂的系统，其包含参数众多，如时间系统、空间坐标系统等。此次设计是针对卫星运动定位的matlab仿真实现，因要求不高，所以对卫星运动做了理想化处理，摄动力对卫星的影响忽略不计（所以为无摄运动），采用开普勒定律及最小二乘法计算其轨道参数，对其运动规律进行简略分析，并使用matlab编程仿真实现了卫星的运功轨道平面、运动动态、可见卫星的分布及利用可见卫星计算出用户位置。 通过此次设计，对于GPS卫星有了初步的认识，对于静态单点定位、伪距等相关概念有一定了解。 关键字：GPS卫星 无摄运动 伪距 matlab仿真 The movement and location of GPS satellite on MATLAB Abstract：Global positioning system is a global, versatility, all-weather advantage of navigation and positioning, timing and speed system, now there has many application in many fields. GPS satellite positioning is a complex system, which includes many parameters, such as time and space coordinates system. This design is based on the matlab simulation of satellite motion and location, because demand is not high, so to do the idealized satellite movement, and ignore the disturbed motion ( so call it non-disturbed motion ).Using the Kepler and least-square method for calculating the parameters of orbital motion, for the characteristics of motion to make a simple analysis, and use the matlab simulation to program achieve the orbital plane of satellite, the dynamic motion, the distribution of visible satellites and using visible satellites to calculate the users’ home. Through the design have primary understanding for the GPS satellite, and understanding the static single-point, pseudorange and so on. Key words：GPS satellite non-disturbed motion pseudorange matlab simulation 目录 1第一章 前言 11.1 课背景 21.2 本课题研究的意义和方法 21.3 GPS前景 4第二章 GPS测量原理 42.1 伪距测量的原理 52.1.1 计算卫星位置 52.1.2 用户位置的计算 52.1.3 最小二乘法介绍 62.2 载波相位测量原理 10第三章 GPS的坐标、时间系统 103.1 坐标系统 103.1.1 天球坐标系 123.1.2 地球坐标系 133.2 时间系统 143.2.1 世界时系统 153.2.2 原子时系统 163.2.3动力学时系统 163.2.4协调世界时 163.2.5 GPS时间系统 18第四章 卫星运动基本定律及其求解 184.1开普勒第一定律 194.2开普勒第二定律 204.3 开普勒第三定律 204.4 卫星的无摄运动参数 214.5 真近点角的概念及其求解 224.6 卫星瞬时位置的求解 25第五章 GPS的matlab仿真 255.1 卫星可见性的估算 265.2 GPS卫星运动的matlab仿真 38结 论 40致谢 41参考文献 42附录 第一章 前言 1.1 课题背景 GPS系统的前身为美军研制的一种子午仪卫星定位系统(Transit），1958年研制，64年正式投入使用。该系统用5到6颗卫星组成的星网工作，每天最多绕过地球13次，并且无法给出高度信息，在定位精度方面也不尽如人意。然而，子午仪系统使得研发部门对卫星定位取得了初步的经验，并验证了由卫星系统进行定位的可行性，为GPS系统的研制埋下了铺垫。由于卫星定位显示出在导航方面的巨大优越性及子午仪系统存在对潜艇和舰船导航方面的巨大缺陷。美国海陆空三军及民用部门都感到迫切需要一种新的卫星导航系统[13]。 1973年12月 ,美国国防部批准它的陆海空三军联合研制新的卫星导航系统: NAVSTAR/GPS。它是英文“Navigation Satellite Timing and Ranging/Global Positioning System” 的缩写词。其意为 “卫星测时测距导航/全球定位系统”,简称 GPS。这个系统向有适当接受设备的全球范围用户提供精确、 连续的三维位置和速度信息 ,并且还广播一种形式的世界协调时(U TC) 。通过遍布全球的（21+3）GPS导航卫星，向全球范围内的用户全天候提供高精度的导航、跟踪定位和授时服务。目前，GPS已在地形测量，交通管理，导航，野外勘探，空间宇宙学等诸多领域得到了广泛的应用[11]。 目前全球共有4大GPS系统，分别是： 美国 GPS ，由美国国防部于 20 世纪 70 年代初开始设计、 研制 ，于1993 年全部建成。1994 年 ，美国宣布在 10 年内向全世界免费提供 GPS使用权 ，但美国只向外国提供低精度的卫星信号。 欧盟 “伽利略”，1999 年 欧洲提出计划 ,准备发射 30 颗卫星 ，组成 “伽利略” 卫星定位系统。 俄罗斯 “格洛纳斯”，尚未部署完毕。始于上世纪 70年代，需要至少 18 颗卫星才能确保覆盖俄罗斯全境;如要提供全球定位服务 ，则需要 24 颗卫星。 中国“北斗”2003 年我国北斗一号建成并开通运行 ，不同于 GPS， “北斗” 的指挥机和终端之间可以双向交流。四川大地震发生后 ,北京武警指挥中心和四川武警部队运用 “北斗” 进行了上百次交流。北斗二号系列卫星今年起将进入组网高峰期 ，预计在 2015 年形成由三十几颗卫星组成的覆盖全球的系统。 1.2 本课题研究的意义和方法 GPS系统是一个很庞大的系统，包含了天文，地理，计算机，电磁学，通信学，信息学等等。通过本文对GPS的学习研究，最重要的还是要学习其原理：卫星运动原理；卫星定位原理；卫星跟踪原理等等。通过基础原理的学习，一方面，可以使我们更进一步的理解卫星运动，定位的实现方法；通过仿真，进一步了解简单定位的方法及其在仿真平台上的实现途径；另一方面，也可以培养我们自学的能力，训练仿真模拟的技巧和方法。 至今，基本上完成了课题的要求，通过不断的注入既定参数，可以更加详细，直观的理解基本的定位原理和实现方法！ 1.3 GPS前景 GPS导航定位以其定位精度高、观测时间短、测站间无需通视、可提供三维坐标、操作简便、全天候作业、功能多、应用广泛等特点著称。 用GPS信号可以进行海、空和陆地的导航、导弹的制导、大地测量和工程测量的精密定位、时间的传递和速度的测量等。对于测绘领域，GPS卫星定位技术已经用于建立高精度的全国性的大地测量控制网，测定全球性的地球动态参数；用于建立陆地海洋大地测量基准，进行高精度的海岛陆地联测以及海洋测绘；用于检测地球板块运动状态和地壳形变；用于工程测量，成为建立城市与工程控制网的主要手段。用于测定航空航天摄影瞬间相机位置，实现仅有少量的地面控制或无地面控制的航测快速成图，导致地理信息系统、全球环境遥感监测的技术革命[4]。 目前，GPS、GLONASS、INMARSAT等系统都具备了导航定位功能，形成了多元化的空间资源环境。这一多元化的空间资源环境，促使国际民间形成了一个共同的策略，即一方面对现有系统充分利用，一方面积极筹建民间GNSS系统，待2011年左右，GNSS纯民间系统建成，全球 将形成GPS/GLONASS/GNSS三足鼎立之势，才能从根本上摆脱对单一系统的依赖，形成国际共有、国际共享的安全资源环境。世界才可以将卫星导航作为单一导航手段的最高应用境界。国际民间的这一策略，反过来又影响和迫使美国对其GPS使用政策作出更开放的调整。多元化的空间资源环境的确立，给GPS的发发展应用创造了一个前所未有的良好的国际环境。 第二章 GPS测量原理 GPS导航系统的基本原理是测量出已知位置的卫星到用户接收机之间的距离 ,然后综合多颗卫星的数据就可知道接收机的具体位置。要达到这一目的 ,卫星的位置可以根据星载时钟所记录的时间在卫星星历中查出。而用户到卫星的距离则通过纪录卫星信号传播到用户所经历的时间 ,再将其乘以光速得到(由于大气层电离层的干扰 ,这一距离并不是用户与卫星之间的真实距离 ,而是伪距( PR) :当GPS卫星正常工作时 ,会不断地用 1 和 0 二进制码元组成的伪随机码(简称伪码)发射导航电文。GPS系统使用的伪码一共有两种 ,分别是民用的 C/ A 码和军用的 P( Y)码。C/ A 码频率 1. 023MHz ,重复周期一毫秒 ,码间距 1 微秒 ,相当于 300m; P 码频率10. 23MHz ,重复周期266. 4 天 ,码间距0. 1 微秒 ,相当于 30m。而 Y码是在 P码的基础上形成的 ,保密性能更佳。 GPS导航系统卫星部分的作用就是不断地发射导航电文。然而 ,由于用户接受机使用的时钟与卫星星载时钟不可能总是同步 ,所以除了用户的三维坐标 x、 y、 z外 ,还要引进一个Δt 即卫星与接收机之间的时间差作为未知数 ,然后用 4 个方程将这 4个未知数解出来。所以如果想知道接收机所处的位置 ,至少要能接收到 4 个卫星的信号。 2.1 伪距测量的原理 GPS定位采用的是被动式单程测距。它的信号发射书机由卫星钟确定，收到时刻是由接收机钟确定，这就在测定的卫星至接收机的距离中，不可避免地包含着两台钟不同步的误差和电离层、对流层延迟误差影响，它并不是卫星与接受机之间的实际距离，所以称之为伪距。 伪距定位法是利用全球卫星定位系统进行导航定位的最基本的方法，其基本原理是：在某一瞬间利用GPS接收机同时测定至少四颗卫星的伪距，根据已知的卫星位置和伪距观测值，采用距离交会法求出接收机的三维坐标和时钟改正数。它的优点是速度快、无多值性问题，利用增加观测时间可以提高定位精度；缺点是测量定位精度低，但足以满足部分用户的需要。 2.1.1 计算卫星位置​​​ 读入导航电文后 首先根据需要调用广播轨道 1 至广播轨道 5 上的数据 然后依次计算卫星的平均角速度 归化时间 平均近点角 需要注意的是进行真近点角计算时要同时计算正弦和余弦以得到正确象限内的角 计算经校正的升交点精度时需要用到地球旋转速率 在 WGS-84 中这一常数为[10]: Ωe=7.2921151467 10-5rdd/s。 (2-1) 2.1.2 用户位置的计算 首先利用近似的用户位置与伪距观测值计算出一个近似伪距 利用该近似伪距可以计算出部分值 然后计算出系数 并生成一个Nx4 的矩阵 为参与运算的卫星数 最后算出用户位移的坐标上述过程根据需要可以将计算出的用户坐标作为近似值 反复迭代直至符合精度要求[10]。 当在进行迭代的过程当中，如果所给定的用户的初始值越接近用户的实际值，则迭代的次数就越少。当我们可见的卫星多于四颗的时候，我们可以用以下介绍的最小二乘法原理带到上面的公式当中去计算。 2.1.3 最小二乘法介绍 当我们在一个地方同时可见的卫星如果多于四颗（GPS卫星的轨道设计和运动时间的安排使得用户在地球的任意位置（两极个别地点除外），都能够看到4——11颗的卫星）的时候，我们可以用最小二乘法去解算未知数，这样，充分的利用了已知的数据信息，使得结果的偏差最小化。 例如：对于下面的方程： (2-2) 如果令 a1 b1 c1 x H= a2 b2 ,C= c2 ,X= (2-3) a3 b3 c3 y 使用最小二乘法，用C,H表示X.令 = （a1x+b1y-c1）^2+ (a2x+b2y-c2)^2+(a3x+b3y-c3)^2 (2-4) = 2(a1x+b1y-c1)a1+2(a2x+b2y-c2)a2+2(a3x+b3y-c3)a3 (2-5) = 2(a1x+b1y-c1)b1+2(a2x+b2y-c2)b2+2(a3x+b3y-c3)b3 (2-6) 整理得到 (a1^2+a2^2+a3^2)x+(a1b1+a2b2+a3b3)y=a1c1+a2c2+a3c3 (2-7) (a1^2+a2^2+a3^2)y+(a1b1+a2b2+a3b3)x=b1c1+b2c2+b3c3 (2-8) 写成矩阵的形式就是： a1 a2 a3 a1 b1 x a1 a2 a3 c1 a2 b2 == c2 b1 b2 b3 a3 b3 y b1 b2 b3 c3 哪么就能够得到 x == (2-9) y 2.2 载波相位测量原理 载波相位观测方法是GPS接收机用接收到的卫星载波（L1:154f0,19.032cm;L2:120f0,24.42cm）与本地接收机产生的本振参考载波产生的相位差来计算的.(GPS所接收到得载波相位是不连续的，所以在进行相位测量的时候，先要进行解调工作，把调制在载波上面的测距码和导航电文去掉，通过码相关等方法重新获取载波)。以 表示k接收机在时刻tk所接收到的第j颗卫星接收到的载波相位的值；以 表示k接收机在时刻tk本地载波信号的相位值，则接收机在接收机钟面时刻tk时观测j卫星所取得的相位观测量可写为[7]： (2-10) 图2-1 载波相位测量原理图 如上图1所示，在初始时刻t0,测得小于一周的相位差为 ，其整周数为 ，此时包含整周数的相位观测值为： == + == - + (2-11) 接收机继续跟踪卫星信号，不断地测量小于一周的相位差 ，并利用整波计数器记录从t0到ti时间内的整周数变化量INT( ),只要卫星j从to到tj时间内信号没有中断，则整周模糊数 就为一个常数，任意时刻ti卫星到k接收机的相位差为： INT( ) (2-12) 这样，观测量就包含了相位差的小数部分和累计的整数部分的整周数。 载波相位的观测方程 假设在GPS系统时刻Ta(卫星a时刻)卫星Sj发射的载波信号相位为 ，经过传播的延迟 后，在GPS系统时刻Tb（接收机tb时刻到达接收机）。接收机产生的本地载波相位为 ，根据（2-10）得到：在Tb时刻，载波相位的观测值为 (2-13) ，考虑到卫星与接收机和系统时间的差值，Ta=ta+ ,Tb=tb+ , 则有： (2-14) 由于卫星和接收机的频率都比较稳定，所以在一个小的时间间隔里面，我们可以近似的理解在时间[t,t+ t]内，有 (2-15) （上式是考虑在频率前提下，所以没有在f前面乘以2 ） ，因为卫星到接收机有一个传播的延迟，即Tb=Ta+ ,所以有： (2-16) 由2-15和 2-16 代入到 2-14 得到： (2-17) ，如果同时考虑到传播延迟电离层和对流层的影响（ ），有 = ， 为卫星到接收机的距离。则有： (2-18) 将 2-12 代入上个式子，考虑到 ，得到以米为单位的测量的载波相位为： (2-19) 这样，上式即为接收机k对卫星k的载波相位的以米为单位的观测方程式。 第三章 GPS的坐标、时间系统 3.1 坐标系统 GPS定位测量当中，要用到两种坐标系，即天球坐标系和地球坐标系。天球坐标系是指坐标原点和各坐标轴的指向在空间是保持不变的，可以很方便的描述卫星的运动和状态。而地球坐标系则是与地球体相关联的坐标系，用于描述地面测量站的位置。下面就天球坐标系和地球坐标系做简要的说明。 3.1.1 天球坐标系 天球就是指的是以地球质心为中心，半径无穷大的理想球体。在这个系统当中，我们会涉及到几个参考点，线，面[9]。 1：天轴和天极：天轴是指地球自转的延伸直线，天轴和地球表面的交点叫做天极P,与地球北极相对应的是北天极Pn,与地球南极相对应的是南天极Ps.天极并不是固定的，有岁差和章动的影响，这个时候叫做真天极，而无岁差和章动影像的天极叫做平天极。 2：天球赤道：通过地球质心并与天轴垂直的平面与地球表面的交线叫做天球赤道。 3：天球子午面：包含天轴并通过天球面上任意一点的平面。 4：黄道：地球绕太阳公转的轨道平面和天球表面相交的大圆，黄道平面和天球赤道平面的夹角叫做黄赤交角，为23.5度。 5：黄极：过天球中心垂直于黄道平面的直线与天球表面相交的点，它分为黄北极和黄南极，分别用Kn和Ks表示。 6：春分点：指太阳由南向北运动的时候，所经过的天球黄道和天球赤道的交点。春分点和天球赤道面是建立天球坐标系的基准点和基准面。 图3-1 春分点 7：岁差和章动：岁差指的是平北天极以北黄极为中心，以黄赤交角为半径的一种顺时针圆周运动。 章动指的是真北天极绕平北天极做得椭圆型运动。 图3-2 岁差和章动 天球坐标系分为两种：真天球坐标系和平天球坐标系. 真天球坐标系的原点为地球的质心M，Z轴指向真北天极Pn,X轴指向春分点,Y轴垂直于XMZ平面。 平天球坐标系的原点为地球的质心M，Z轴指向平北天极Pn,X轴指向春分点,Y轴垂直于XMZ平面。 上述两种坐标系的差别在于他们选取了不同的北天极的位置，故要是由平天极坐标系到真天极坐标系的转换，就必须考虑岁差和章动旋转所影响的情况。换句话说就是要考虑到岁差旋转和章动旋转地影响。 3.1.2 地球坐标系 地球坐标系也可以分为两种：即平地球坐标系和真地球坐标系。 1 平地球坐标系：它的地极位置采用国际协议地极原点CIO(由1900到1905年测定的平均纬度所确定的平均地极位置)。 原点：地球质心M。 Z轴：指向CIO。 X轴：指向格林威治起始子午面与地球平赤道的交点。 Y轴：垂直于XMZ平面。 2 真地球坐标系 原点：地球质心M。 Z轴：指向地球瞬时极。 X轴：指向格林威治起始子午面与地球瞬时真赤道的交点。 Y轴：垂直于XMZ平面。 图3-3 世界地心坐标系 瞬时真天球坐标到瞬时真地球坐标的转换 这两种坐标的差异就在于X轴的指向是不同的。前者指向的是真春分点，而后者指向的是格林威治起始子午面与地球瞬时真赤道的交点。两者之间的夹角称为对应的平格林威治起始子午面的真春分点时角Ω。故仅仅需要绕Z轴旋转这个角度Ω，就能够做到二者的相互转换。相应的转换转动矩阵为： CosΩ -SinΩ 0 Rz(Ω )= SinΩ CosΩ 0 3-1 0 0 1 综合上面的，可以得到以下的结论：在GPS 定位系统所用的空间坐标系统当中，我们一般采用天球坐标去研究卫星的空间运动，而采用地球坐标去研究地面监控站点，他们之间的转换问题一般可以按照下面的步骤来分析： 3.2 时间系统 时间系统是卫星定位测量过程中的一个重要概念。现时的GPS测量的方法 是通过接收和处理GPS卫星发射的无线电信号，以确定用户接收机和观测卫星 间的距离，然后通过一定的数学方法以确定接收机所在的具体位置，为得到接 收机和卫星的准确距离，必须获得无线电信号从卫星传输至接收机这一过程中 的精确时间，因而利用卫星技术进行精密的定位和导航，必须要获得高精度的 时间信息，这需要一个精确的时间系统。现行的卫星定位测量中与之紧密相关 的时间系统有三种:世界时，原子时和动力学时。 3.2.1 世界时系统 以地球自转为基准的一种时间系统。根据不同的空间参考点，又可分为恒星 时，太阳时，世界时三种。 1.恒星时 选定春分点(地球赤道平面与其绕太阳公转轨道的一个交点)作为参考点，由 该点的周日视运动所确定的时间，即为恒星时(siderealTime，sT)。规定从春分 点连续两次经过本地子午圈的时间间隔为一恒星日，其1/24为一恒星时，由于 其定义涉及到地方子午圈，因而恒星时具有地方性，又称地方恒星时。当从格 林尼治子午线上观测时，所得的恒星时称为格林尼治恒星时。由于地球自转受 岁差、章动的影响，春分点的空间位置并不唯一，有真春分点和平春分点之分， 这导致恒星时可分为真恒星时和平恒星时，因而对格林尼治恒星时有格林尼治 真恒星时(GAST)和格林尼治平恒星时(GMsT)这两者之间的关系为: （3-2） 其中，△笋为黄经章动，‘为黄经交角。 2.太阳时 以真太阳周日视运动所确定的时间称为真太阳时。但据天体运动的开普勒定律，太阳视运动的速度不是均匀的，以真太阳作为观察地球自转的参考点，不符合时间系统的基本要求，因而假定了一个参考点，其在天球上的视运动速度，等于真太阳周年运动的平均速度，这个假定的参考点，在天文学上被称为平太阳。以平太阳连续两次经过本地子午圈的时间间隔，定义为一个平太阳日，其1/24为一平太阳时(MeansolarTim。，MT)。与恒星时一样，平太阳时也具有地方性，常称地方平太阳时。 3.世界时 以地球上格林尼治子午圈所对应的平太阳时且以平子夜起算时间系统，称为世界时(universalTime，uT)。世界时与平太阳时的尺度标准完全一致，仅仅是起算点有所不同。若有编表示平太阳相对格林尼治子午圈的时角，定义有世界时UTO可表示为 : (3-3) 由于地球自转的不均匀性，使地球自转轴产生了极移现象因而UTO并不均匀，为补偿这一缺陷，国际天文联合会在世界时中引入地轴极移修正△兄和地球自转变化的季节性改正参数双由此可得世界时UTI和uTZ: (3-4) (3-5) 其中观测瞬时地极相对国际协议地极原点(CIO)的极移修正△兄的表达式为 (3-6) 式中X’，厂为观测瞬间的极移分量;凡，汽分别为天文经度和纬度。地球自转速度的季节性变化改正△兀有如下的经验公式: (3-7) t为自本年起始日起算的年小数部分(即为计算时年积日与该年全年积日的比例)。上述修正并不能完全消除地球自转速度变化率和地球自转季节性变化的影响，故而UT:并不是严格均匀的时间系统。 3.2.2 原子时系统 原子时以物质内部原子跃迁时所辐射和吸收的电磁波频率来定义的，其秒长定义为:位于海平面上的艳原子侧”基态两个超精细能级，在零磁场中跃迁辐射振荡9，192，631，770周所持续的时间，为1原子秒，该原子时秒作为国际制秒(sI)的时间单位。原子时的起点是定在1958年1月1日0时0分0秒(UT2)，但与之又有微小误差，关系为: (3-8) 原子时具有很高的稳定性和复现性，是现时段最为理想的时间系统。许多国家都建立了各自的原子时系统，国际时间局为消除差异，对100座时钟作了对比分析，利用数据处理推算出了统一的原子时系统——国际原子时(hiternationalAiomicTime，TAD。在目前的导航定位系统中，均采用了原子时作为其高精度的时间基准。 3.2.3动力学时系统 动力学时(DynamicTime，DT)是天体力学中用以描述天体运动的时间单位。当以太阳系质心建立起天休运动方程时，所采用的时间参数称为质心力学时(BaryeeniricDynamicTime，TDB);当以地球质心建立起天体运动方程时，所采用的时间参数称为地球力学时仃℃仃estrialDynamicTime，TDT)。TDT所采用的基本单位为sl，与原子时一致。国际天文学联合会定义1977年1月1日TAI与TDT的严格关系为: TDT=TAI+32.184(s) (3-9) 3.2.4协调世界时 原子时尺度均匀稳定，但与人类日常生活紧密相关的是以地球自转为基础的世界时，在很多的科学研究中均采用的是世界时。世界时受地球速度长期性渐慢的影响，逐渐比原子时慢，为避免两都之间误差的扩大，自1972年起，国际上开始采用一种以原子时秒子为基础，在时刻上尽量接近于世界时的一种折衷的时间系统，称为协调世界时(eoordinateuniversalTime，uTe)。其引入了闰秒的概念，当协调时与世界时的时刻相差超过士0.9(s)时，便于协调时中引入闰秒士l(s)，闰秒一般于12月31日或6月30日加入。协调时与TAI的关系如下: TAI=UTC+n×1(s) (3-10) 其中，n为调整参数，其值由国际地球自转服务组织(lERs)发布。uTc是目前几乎所有国家发布时号的标准，相互之前的同步误差约为士0.2ms。 3.2.5 GPS时间系统 全球定位系统(GPS)为保证导航和定位精度，建立了专门的时间系统—GPS时间系统(GPST)。其隶属于原子时系统，秒长采用国际制秒sI，但不同于TAI，两都之前的关系为: TAI-GPST=19(s) (3-11) 据协调时与TAI的关系可得: GPST=UTC+n×1-19(s) (3-12) 第四章 卫星运动基本定律及其求解 卫星在空间绕地球运动的时候，除了受到地球重力场引力的作用外，还受到了太阳，月亮和其它的天体引力以及太阳光压，大气的阻力和地球潮汐力的影响。卫星的实际运动轨道非常的复杂，很难用非常精确的数学模型加以描述。在各种力作用对卫星影响的过程当中，以地球的引力场的作用最大，而其它力的影响则相对的小得多。通常把作用到卫星上的力按其影响的大小分成两部分：一类是中心力；一类是摄动力，也称为非中心力。假定地球为均匀球体的地球引力，称为在心力，它决定了卫星运动的基本规律和基本特征，由此决定地球的轨道，可以视为理想的轨道。非中心力包括地球非球形对称的作用力，日、月引力，大气阻力，光辐射压力以及地球的潮汐力等。摄动力的作用，使卫星偏离了既定的理想轨道。而在它影响下，卫星的运动称为卫星的受摄运动。而上述理想状态的卫星运动则称为无摄运动。卫星在地球的引力场当中所做的无摄运动，也称为开普勒运动，其规律可以由开普勒三大定律来描述。[7] 4.1开普勒第一定律 开普勒第一定律：卫星运动的轨道是个椭圆，而该椭圆的一个焦点和地球的质心重合。 这一个定律表明了，在中心引力的作用下，卫星绕地球轨道运行的轨道面，是一个通过地球质心的精致平面。轨道椭圆一般称期为开普勒椭圆，其形状和大小都不变。在轨道上，卫星离地球质心远的一点叫做远地点，近的一点就做近地点。轨道图形可以表示为如下图5： ms 远地点P’ 近地点P 卫星绕地球质心运动的轨道方程为： R= (4-1) 在该式当中，R是卫星的地心距离；as为开普勒椭圆的长半径；es为开普勒椭圆的偏心率；fs为真近点角，它描述了任意时刻，卫星在轨道上面相对于近地点的位置，是时间的函数，其定义见上图所示。 开普勒定义定律阐述了卫星运动轨道的基本形态及其与地心的关系。 4.2开普勒第二定律 开普勒第二定律：卫星的地心向径，即地球质心与卫星质心间的距离向量，在相同的时间内所扫过的空间面积是相等的。（如下图6所示） （图2） 与任何其它的运动物体一样，在轨道上面运动的卫星，也具有两种的能量：位能和动能。位能就是指仅仅受到地球重力场的影响，其大小和卫星的在轨高度有关。在近地点其位能最小，而在远地点其位能最大。卫星在任一个时刻t所具有的位能为 (G为万有引力常量，M为地球的质量,ms为卫星的质量)。动能则是由卫星的运动所引起的，其大小是卫星的运动速度的函数。如果取卫星的运动的速度为vs,则其动能为 。根据能量守恒定律，卫星的势能与动能的总量是不变的，即 =常量 (4-2) 因此，当卫星运行到近地点的时候，其动能最大；在远地点的时候，其动能最小，由此，开普勒第二定律所包涵的内容是：卫星在椭圆轨道上的运行速度是不断变化的，在近地点处的速度最大，而在远地点的速度最小。 4.3 开普勒第三定律 开普勒第三定律：卫星运动周期的平方与轨道椭圆长半径的立方之比为一个常数，而该常数等于地球引力常数和地球质量的乘积GM的倒数。 其数学表达式为： (4-3) 在这个式子当中， 为卫星的运行周期.如果我们假设卫星的平均角速度为N,则有: N= (4-4) 于是，开普勒第三定律 4-2就可以写成： (4-5) 或者表示为常用的形式：N= (4-6) 显然，当开普勒的长半径确定了过后，卫星运动的平均角速度就得到了确定，且保持不变。 4.4 卫星的无摄运动参数 卫星的无摄运动，一般的可以由下面的6个参数（图7）来描述： As----------------------------- 卫星轨道的长半径 Es ----------------------------卫星轨道的偏心率 Ω-------------------升交点的赤径 i--------------------卫星轨道面的倾角 ωs--------------近地点角距，即升交点与近地点的夹角 fs-----------------卫星的真近点角，在轨道平面上为卫星与进地点的地心角距。 图4-3开普勒轨道参数 当这6个参数一旦确定后，卫星在任意瞬时的相对于地球的空间位置及其速度，就被唯一的确定了！ 4.5 真近点角的概念及其求解 在描述卫星无摄运动的6个参数当中，只有fs是关于时间的函数，其他的都是一般的参数。所以，计算卫星瞬时的位置的关键，计算出参数fs，并由此确定卫星的空间位置及其和时间的关系。 为此，需要引进两个参数Es和Ms去计算真近点角。 Es:偏近点角，如果定义过卫星质心做平行与椭圆短半轴的直线，M’为该直线与近地点到椭圆中心连线的交点，则椭圆平面上近地点P到M’的圆弧所对应的圆心角就是Es。 Ms:平近点角。它是一个假设量，如果卫星在轨道运行的平速度为n，则平近点角定义为： Ms=n(t-t0) (4-7) t0为卫星过近地点的时刻，t为观察卫星的时刻。 由上面的式子知道，卫星的平近点角仅仅为卫星平均速度的时间的函数，对于一个确定的卫星来说，这个参数可以认为是常数。 其中Ms与Es有关系如下： Ms=Es-essinEs (4-8) 为了计算卫星的瞬时速度， 需要确定卫星运行的真近点角fs。由于有以下的关系成立： ascosEs=rcosfs+ases (4-9) 于是将上式带入到 (4-1)中就得到： Cosfs= (4-10) 或者得到以下常用的形式： Tan( )= (4-11) 4.6 卫星瞬时位置的求解 对于任意的观测时刻，根据卫星的平均运行速度n,根据4-9，4-10，4-11，便可以唯一确定真近点角fs。这样，卫星于任一观测历元t，相对于地球瞬间空间的位置便可以随之确定。 若以直角坐标的原点 与地心M重合， 轴指向近地点且垂直于轨道的平面， 轴在轨道平面上垂直 轴构成右手关系。于是，卫星任意时刻的轨道坐标可以表示成为： cosfs = sinfs (4-12) 0 而由上面的分析，可以得到： = cosE- = SinE = 0 (4-13) 而要把这个轨道坐标系坐标表示成为天球坐标的话，由于他们的坐标原点都是地球的质心，但是坐标轴的指向是不相同的，为了使他们的坐标轴相同，应该将坐标系（ , , ）依次做下面的变化旋转： （1） 绕 轴顺时针旋转角度ωs，使 轴的指向由近地点变为升交点。 （2） 绕 轴顺时针旋转角度i，使 轴与Z轴相同。 （3） 绕 轴顺时针旋转角度Ω，使 轴指向春分点。 实现上述三步的旋转矩阵分别为R1，R2，R3，即： Cosωs -Sinωs 0 R1= Sinωs Cosωs 0 (4-14) 0 0 1 1 0 0 R2= 0 Cosi -Sini (4-15) 0 Sini Cosi CosΩ -SinΩ 0 R3= SinΩ CosΩ 0 (4-16) 0 0 1 于是得到了在天球坐标系下面的卫星位置坐标可以表示为： X CosE-es Y =R3R2R1 =R3R2R1as SinE (4-17) Z 0 利用转换关系Rz(Θg)得到相应的地球坐标系的坐标了！ CosΘg SinΘg 0 而Rz(Θg) = -SinΘg CosΘg 0 (4-18) 0 0 1 第五章 GPS的matlab仿真 5.1 卫星可见性的估算 当初，卫星星座的设计要求在全球范围内任何时候，任何位置都必须保证至少四颗以上的卫星导航信号。换句话说，并不是所有的卫星都能够被一个用户所看见。一颗卫星信号能否被接收与下列因素是有关系的[7]： 1 > 地球是否影响了该GPS卫星信号的传播。下图8说明了地球对GPS信号的影响。如果卫星处于图中的阴影部分，则对图中的飞机是不能够接收到该卫星的信号的! 2> GPS接收机是否位于该GPS卫星发射天线的范围内。这种情况主要针对的是航天器上面的GPS接收机。GPS卫星信号的发射张角大约为21.3度，大于卫星到水平面的张角13.9度，这样就保证了一些飞行高度较高的航天器在高空可以更多的接收到GPS卫星的信号。但是对于那些超出发射角的飞行器就收不到信号了。 3>利用卫星—地心—用户之间的张角，可以估算出来可卫星信号能否被收到。如果这个角度小于90度的话，就可以收到，反之不能够收到。特别的，当这个角度刚刚为90度的时候，我们一般认为是收不到的。 5.2 GPS卫星运动的matlab仿真 程序主体见附录 在进行仿真之前，有几个子程序段需要说明一下： function plot3c(x,y,z,color) switch (color) case 0 plot3(x,y,z,'w--'); case 1 plot3(x,y,z,'r--'); case 2 plot3(x,y,z,'g--'); case 3 plot3(x,y,z,'c--'); case 4 plot3(x,y,z,'m--'); case 5 plot3(x,y,z,'y--'); case 6 plot3(x,y,z,'b--'); case 7 plot3(x,y,z,'k--'); otherwise end 它主要用来画出卫星轨道的曲线，用不同的颜色：w代表的白色,其他的依次代表的是红色,绿色,青绿色.品红色,黄色,蓝色,黑色. function boxplot3(x,y,z,Lx,Ly,Lz,color) x0=x-Lx/2;y0=y-Ly/2;z0=z-Lz/2; x=[x0 x0 x0 x0 x0+Lx x0+Lx x0+Lx x0+Lx]; y=[y0 y0 y0+Ly y0+Ly y0 y0 y0+Ly y0+Ly]; z=[z0 z0+Lz z0+Lz z0 z0 z0+Lz z0+Lz z0]; index=zeros(6,5); index(1,:)=[1 2 3 4 1]; index(2,:)=[5 6 7 8 5]; index(3,:)=[1 2 6 5 1]; index(4,:)=[4 3 7 8 4]; index(5,:)=[2 6 7 3 2]; index(6,:)=[1 5 8 4 1]; for k=1:6 plot3c(x(index(k,:)),y(index(k,:)),z(index(k,:)),color) hold on end 它主要是用来表示用户的空间位置的。 function drawearth(time) %time 是参数 %利用这个参数,可以绘制一个看起来是旋转的地球 r=6400; time=0; j1=[0:pi/10:2*pi]; w1=[-pi/2:pi/10:pi/2]; L1=length(w1); L2=length(j1); for n=1:L1 z=ones(L2,1); z=z*r*sin(w1(n)); temp=r*cos(w1(n)); x=temp*sin(j1); y=temp*cos(j1); plot3(x,y,z); hold on; grid; end %figure(3); unit=ones(1,1); z0=ones(1,1); x0=ones(1,1); y0=ones(1,1); for n=1:L2 %n=7; for m=1:L1 temp=w1(m); temp2=j1(n)+time*pi/12; z=r*sin(temp); x=r*cos(temp)*sin(temp2); y=r*cos(temp)*cos(temp2); z1=unit*z; x1=unit*x; y1=unit*y; z0=[z0 z1]; x0=[x0 x1]; y0=[y0 y1]; end z0(:,1)=[]; x0(:,1)=[]; y0(:,1)=[]; plot3(x0,y0,z0); axis equal; axis off; hold on; boxplot3(0,0,0,100,100,100,7); %标示出来地球球心的位置。方便观察 %以下是在同一个坐标系当中标识出空间直角坐标系的三个与天球相交的三个坐标点，以便观察。红色表示X轴，绿色Y轴，褐红色Z轴。 tempx=6400; tempy=0; tempz=0; cube=100; boxplot3(tempx,tempy,tempz,cube,cube,cube,1); tempx=0; tempy=6400; tempz=0; cube=100; boxplot3(tempx,tempy,tempz,cube,cube,cube,2); tempx=0; tempy=0; tempz=6400; cube=100; boxplot3(tempx,tempy,tempz,cube,cube,cube,4); end function drawsatellite(moveX,moveY,moveZ,color); Length=500; Width=500; Height=500; boxplot3(moveX,moveY,moveZ,Length,Width,Height,color); boxplot3(moveX,moveY,moveZ,Length,Width*5,Height/10,color); 它的作用就是画出一个在上述地球坐标图当中的卫星。 仿真程序一：绘制卫星的轨道平面 程序见附录： 图5-2 卫星轨道平面仿真图 仿真程序二：单颗卫星不同时刻的动态仿真 程序见附录： 图5-3 单颗卫星动态仿真图 由上图10可知，卫星的天空瞬时位置是随着时间的变动而发生变化的。对于我们在地球上面的一个用户来说，一天当中的不同时刻看到的卫星是不相同的。这个仿真程序的功能实际上就是仿真了在一个轨道上面的卫星在不同的时候（这里以一个小时为一个观测时元，进行动态的在屏幕上显示器位置）通过改变时间可以显示出不同时间内的卫星的瞬时位置），而且，通过matlab当中的三维旋转图标，我们可以很清晰地从不同角度看到卫星和坐标原点（这里用黑点表示的的相对位置的变化）。哪么，在其他轨道面上的卫星的运动也可以类似的模拟出来，这里不在重复！ 仿真程序三：卫星在某个时刻的全轨道平面的分布和可见卫星 程序见附录： 图5-4 全轨道平面的图形 其在轨坐标分别如下： SatellitePosition = 1.0e+004 * 1.7746 1.7572 0.7365 0.0001 -1.2161 0.9732 2.1091 0.0001 2.2883 0.2975 -1.4132 0.0001 -2.3882 -0.8836 0.7869 0.0001 -0.3681 -2.5255 0.7012 0.0001 -0.1323 2.7059 -0.0000 0.0001 -1.5424 0.4273 2.1704 0.0001 1.4582 0.4926 -2.1144 0.0001 -2.4090 0.6583 0.7365 0.0001 1.0759 1.1599 -2.1757 0.0001 1.3324 -1.8178 1.4392 0.0001 -1.3225 1.7688 -1.3835 0.0001 1.2127 -0.9774 2.1091 0.0001 -1.8766 0.3838 -1.9058 0.0001 -2.4302 -0.9295 -0.7516 0.0001 0.4443 1.8187 -1.8463 0.0001 -0.6718 2.4562 -0.7869 0.0001 1.4000 -1.3073 1.9058 0.0001 -1.4992 0.4223 -2.1091 0.0001 -0.8706 -2.0297 -1.3835 0.0001 -1.7722 -0.6021 -1.9217 0.0001 -0.2481 -1.5812 -2.1704 0.0001 -0.5421 1.7915 1.9217 0.0001 1.9376 -1.5756 -0.7365 0.0001 这个程序仿真了在某个时刻（在程序里面是在时刻timenow=0）的全部24颗卫星的轨道图形，以及对于用户来说在这个时刻可以看到的卫星。在程序当中，我们假定了用户的位置坐标是（6400，3352，5410）。通过改变用户的不同位置，可以在同一个时刻看到不同（4到11颗）的可见卫星。对于在这个程序当中的可见卫星的及时在轨坐标如下： 图5-5 用户可见的卫星分布 SatellitePosition = 1.0e+004 * 1.7746 1.7572 0.7365 0.0001 -1.2161 0.9732 2.1091 0.0001 2.2883 0.2975 -1.4132 0.0001 -2.3882 -0.8836 0.7869 0 -0.3681 -2.5255 0.7012 0 -0.1323 2.7059 -0.0000 0.0001 -1.5424 0.4273 2.1704 0.0001 1.4582 0.4926 -2.1144 0 -2.4090 0.6583 0.7365 0 1.0759 1.1599 -2.1757 0 1.3324 -1.8178 1.4392 0.0001 -1.3225 1.7688 -1.3835 0 1.2127 -0.9774 2.1091 0.0001 -1.8766 0.3838 -1.9058 0 -2.4302 -0.9295 -0.7516 0 0.4443 1.8187 -1.8463 0 -0.6718 2.4562 -0.7869 0 1.4000 -1.3073 1.9058 0.0001 -1.4992 0.4223 -2.1091 0 -0.8706 -2.0297 -1.3835 0 -1.7722 -0.6021 -1.9217 0 -0.2481 -1.5812 -2.1704 0 -0.5421 1.7915 1.9217 0.0001 1.9376 -1.5756 -0.7365 0.0001 仿真程序四：用可见卫星计算用户位置 程序见附录 多于四颗的卫星用最小二乘法逼近计算。对于前面的假设用户位置（6400，3352，5410）进行计算，结果如下： 图5-6 用户位置的计算 由图13可见，根据最小二乘法的原理，计算出来了用户在一个时刻的位置坐标，并把它表示在了这个天球坐标系当中。 calculaterecord = 1.0e+003 * 0 0 0 7.3399 3.8857 6.1388 6.3034 3.2976 5.3346 6.3914 3.3470 5.4034 6.4009 3.3525 5.4107 6.4001 3.3520 5.4101 先后经过了六次的迭代算法，吧用户的计算位置一步一步的逼近了用户的实际位置，根据部同的精度要求，我们运算的量的大小也有不同。这点，可以根据程序当中的参数Error的设定而有所出入。 图5-7 用户位置的计算 通过对用户仿真计算的进一步放大，可以很清楚地看到：在每一次迭代的过程当中，我们都实时的把每一步迭代的值也表示在了同一个坐标轴里面：目的就是很直观的反应出我们基于最小二乘法的伪距算法的主要思路。如下图所示： 图5-8 用户位置的计算 上图15当中的黑色箭头从上到下一次表示的是第一次到第六次计算出来的用户的位置，第七箭头表示的用户的实际位置（用白色的柱体表示）。 结 论 通过这次的设计，学习到了很多，因为GPS卫星运动及定位以前没有接触过，所以本次设计难度很大，但在最后都还是通过查找资料或者老师同学帮助也都一一克服。 GPS卫星是一个很复杂的系统，其内容很多，比如它的坐标系、时间系统等，其坐标系分为天球坐标和地球坐标。天球坐标描述卫星的运动位置和状态，地球系统则描述地面测站的位置。GPS卫星作为一个高空动态已知点，其位置是随时间不断变化的，因此，在给出卫星运行位置的同时，必须给出相应的瞬时时刻。所以时间系统是GPS卫星计算里面很重要的一环。要计算卫星的运动轨迹及位置，所需要的参数很多，卫星在空间绕地球运动时，除了受到地球重力场的引力作用外，还受到太阳、月亮和其他天体引力以及太阳光压、大气阻力和地球潮汐等因素影响。通常把作用于卫星上的各种力，按其影响的大小分为两类：中心力和摄动力。因为此次设计要求的精度不高，所以在设计中摄动力忽略不计。卫星在地球引力场中无摄运动，也称为开普勒运动。所以本次设计中主要用开普勒定律计算其位置轨迹，在运算中还加入了最小二乘法计算参数数据。此次设计中卫星的动态仿真，计算中主要遇到的问题是对于特定轨道的卫星在不同时刻的近地点幅角的计算。所以，根据卫星轨道的星座参数，找出在某时刻的近地点幅角（一般是零时刻），然后在此基础上以时间为基础进行叠加，确定出每个时间点的瞬时近地点幅角值。 基于面向对象和程序设计思想，利用matlab编程了GPS定位仿真程序，该程序从自定义的星座数据文件中建立起卫星星座。设计一个子函数，命名为calculateuserposition.m。其功能为其他程序得出一个全局变量satelliteposotion，把它和其他程序放入同一个文件夹中，以便其他程序使用。此次设计程序主要包含四个模块： 1> 绘制卫星的轨道平面，以便于直观了解其轨道； 2> 单颗卫星的动态仿真，了解其不同时刻的位置变化，通过matlab当中的三维旋转图标，可以很清晰的从不同角度看到卫星和坐标原点； 3> 卫星在某时刻的全轨道平面分布和可见卫星，该程序仿真了在某时刻的全部24颗卫星的轨道图形，以及对于用户来说在该时刻可见卫星位置； 4> 利用可见卫星计算用户位置，利用最小二乘法原理，计算用户某一时刻的位置坐标。 通过matlab仿真卫星运动及轨迹，对其有了更加深入的认识和了解，对GPS也有了更直观的认识，对于以后遇到GPS的应用时有了一个简单的基础。 致谢 在论文完成之际，我怀着激动和喜悦的心情向所有关心和帮助过我的老师和同学表示万分的感谢！ 首先要感谢的是帮助我完成本次设计的葛青老师。在本次设计中，我遇到了许多知识和技术方面的疑问，葛老师都一一帮我解答，并提出了许多具有创新意义的建议。她渊博的知识，严谨的工作态度，对工作尽心尽职的精神和平易近人的品格都深深感染了我，将会让我受益终身。在此，我向葛老师表示深深的感谢，并致以崇高的敬意！ 同时还要感谢我的班主任孟玲华老师、周伟老师。大学四年期间，在大一大二时孟老师作为我们班主任，她给以我无私的关怀，是她帮助我成长，让我从一个无知的高中生脱变到具有独立思考能力的大学生，并提供给我许多锻炼自身能力的机会。在大三大四的时候周伟老师出任我们班主任，他作为一个专业课老师，给了我在学习上很大的帮助。在后此外还要感谢所有教育和关怀过我的老师们，还要感谢成都理工大学核技术与自动化工程学院的领导们四年来给予的多方面的关心和帮助! 感谢所有的朋友和同学，让我度过了一个快乐而又充实的大学时光！ 最后要感谢我的父母和家人，他们总是尽最大能力提供给我最好的条件，让我无忧无虑的去学习。是他们的鼓励和支持，是他们深深的关爱，让我得以顺利的完成了学业。 参考文献 [1] 周忠谟,易杰军.GPS 卫星测量原理与应用 [M] .北京:测绘出版社,1997. [2] 李征航,黄劲松.GPS 测量与数据处理 [M] .武汉:武汉大学出版社,2005. [3] Elliott D.Kaplan. Understanding GPS: principles and applications [M] .邱致和,王万义,译.北京:电子工业出版社,2002. [4] 马中元.GPS卫星全球定位系统简介[J]. 科技经济市场 2002（8） [5] 郭秋英,胡振琪.GPS 卫星坐标的计算 [J] .全球定位系统,2006 (3) . [6] 夏振纯,金志华,田蔚风.基于MATLAB 的 GPS 选星分析及实验研究 [J] .计算机仿真,2005,22(1) :51-54. [7] 杨俊，武奇生编著.GPA基本原理及其Matlab仿真[M].西安：西安电子科技大学出版社 [8]张勤，李家权等编著.GPS测量原理及应用[M].北京：科学出版社 [9]汤长存.GPS卫星定位仿真分析与系统程序设计[D].同济大学工学硕士学位论文 [10]刘琼。岳昊，吴北平.MATLAB在GPS卫星几何分布计算中的应用[C].全国测绘科技信息网中南分网第二十四次学术信息交流会论文集 [11]毛敏.GPS卫星定位系统简介[Z].陕西国防工业职业技术学院学报第二十卷第一期。 [12]郭秋英，胡振琪.GPS卫星坐标的计算[J].2006 [13]朱明.全球卫星定位系统GPS及在现代交通中的应用[J]. [14]楼益栋，施闯，葛茂荣，赵齐乐.GPS卫星实时精密定轨及初步结果分析[J].2008.8 [15]刘永法.运用GPS卫星定位系统，提高车辆利用率[J].2008.第二十七期 [16]赵齐乐. GPS导航星座及低轨卫星的精密定轨理论和软件研究[D] .武汉: 武汉大学 ,2004 [17]Chen K. Real2Time Precise Point Positioning andIt s Potential Applications [ C] . ION GNSS 17th In2ternational Technical Meeting of the Satellite Divi2sion , Long Beach , CA ,2004 [18]Liu Jingnan , Ge Maorong. PANDA Sof tware andIt s Preliminary Result of Positioning and Orbit Det2ermination[J ] . Wuhan University Journal of NaturalSciences , 2003 ,8 (B2) : 603-609 [19] Kouba J , Heroux P. GPS Precise Point PositioningUsing IGS Orbit Projuct s[J ] . GPS Solution , 2001 ,5 (2) : 12-28 [20]www.cdut.edu.cn 附录 程序一（绘制卫星的轨道平面）代码 function drawsatelliteorbit a=26560;% 卫星轨道的长半轴. e=0.02;%e是椭圆的偏心率. E=[0:0.1:2*pi]; x=a*(cos(E)-e); y=a*sqrt((1-e^2))*sin(E); z=0*E; DtoR=2*pi/360; A1=[32.8 92.8 152.8 212.6 272.8 332.8];%卫星星座数据. for k=1:6 A=A1(k)*DtoR;%升交点的经度 B=55*DtoR;%轨道的倾角 C=pi/100;%近地点幅角 %总共有6个卫星轨道平面 R1=[cos(A) -sin(A) 0; sin(A) cos(A) 0; 0 0 1;]; R2=[1 0 0; 0 cos(B) -sin(B); 0 sin(B) cos(B);]; R3=[cos(C) -sin(C) 0; sin(C) cos(C) 0; 0 0 1;]; L1=length(E); R312=R1*R2*R3; Ans=R312*[x;y;z;]; x1=Ans(1,:); y1=Ans(2,:); z1=Ans(3,:); plot3c(x1,y1,z1,k); boxplot3(0,0,0,200,200,200,7); hold on; axis equal; axis off; end 程序二（单颗卫星不同时刻的动态仿真）代码 clear; clc;close all; DtoR=2*pi/360; jiaostep=360/24*DtoR; j=0; a=26560;% 卫星轨道的长半轴. e=0.02;%e是椭圆的偏心率. E=[0:0.1:2*pi]; x=a*(cos(E)-e); y=a*sqrt((1-e^2))*sin(E); z=0*E; drawearth(0); hold on; for time=1:12 axis on; A1=32.8 ;%卫星星座数据. A=A1*DtoR;%升交点的经度 B=55*DtoR;%轨道的倾角 C=pi/100;%近地点幅角 R1=[cos(A) -sin(A) 0; sin(A) cos(A) 0; 0 0 1;]; R2=[1 0 0; 0 cos(B) -sin(B); 0 sin(B) cos(B);]; R3=[cos(C) -sin(C) 0; sin(C) cos(C) 0; 0 0 1;]; L1=length(E); R312=R1*R2*R3; Ans=R312*[x;y;z;]; x1=Ans(1,:); y1=Ans(2,:); z1=Ans(3,:); plot3c(x1,y1,z1,2); hold on; axis equal; axis on; grid on; ctable=10; A=A1*DtoR; B=55*DtoR; C=ctable*DtoR+time*2*pi/24; %近地点幅角 x=a*(cos(C)-e); y=a*sqrt((1-e^2))*sin(C); z=0*C; R1=[cos(A) -sin(A) 0; sin(A) cos(A) 0; 0 0 1;]; R2=[1 0 0; 0 cos(B) -sin(B); 0 sin(B) cos(B);]; R3=[cos(C) -sin(C) 0; sin(C) cos(C) 0; 0 0 1;]; L1=length(E); R312=R1*R2*R3; Ans=R312*[x;y;z;]; x1=Ans(1,:); y1=Ans(2,:); z1=Ans(3,:); drawsatellite(x1,y1,z1,6); M(time)=getframe; M(time+1) = getframe; end axis on; 程序三（卫星在某个时刻的全轨道平面的分布和可见卫星）代码 clear;clc;close all; a=26560; % 卫星轨道的长半轴. e=0.02; temp=0;%便于把卫星的指标点顺序放到矩阵指标当中。 %e是椭圆的偏心率. E=[0:0.1:2*pi]; x=a*(cos(E)-e); y=a*sqrt((1-e^2))*sin(E); z=0*E; timenow=0;%单位是小时 global SatellitePosition global ttum SatellitePosition=ones(24,4); figure(1); drawearth(0); ttum=ones(24,4); hold on; DtoR=2*pi/360; A1=[32.8 92.8 152.8 212.6 272.8 332.8]; drawsatelliteorbit; ctable=[10 50 160 260; %平均近地角 80 180 220 320; 10 130 250 340; 50 150 170 300; 100 210 310 340; 120 140 240 350;]; simple=1; for k=1:6 A=A1(k)*DtoR; %升交点经度 B=55*DtoR; %轨道倾角 for m=1:4 C=ctable(k,m)*DtoR+timenow*2*pi/24; %近地点幅角 x=a*(cos(C)-e); y=a*sqrt((1-e^2))*sin(C); z=0*C; R1=[cos(A) -sin(A) 0; sin(A) cos(A) 0; 0 0 1;]; R2=[1 0 0; 0 cos(B) -sin(B); 0 sin(B) cos(B);]; R3=[cos(C) -sin(C) 0; sin(C) cos(C) 0; 0 0 1;]; L1=length(E); R312=R1*R2*R3; Ans=R312*[x;y;z;]; x1=Ans(1,:); y1=Ans(2,:); z1=Ans(3,:); drawsatellite(x1,y1,z1,k); temp=temp+1; SatellitePosition(temp,:)=[x1 y1 z1 1]; ttum(temp,:)=[x1 y1 z1 1]; hold on; end end %用余弦定理计算地心到用户和用户到卫星证件的夹角要是两者证件的夹角小于90°,就认为改卫星是不可见的. earthcenterpos=[0 0 0]; userposition=[6400 3352 5410]; %假设一个用户的位置 for k=1:24 temp=SatellitePosition(k,1:3)-userposition; Dist1=temp*temp'; temp=userposition-earthcenterpos; Dist2=temp*temp'; temp=SatellitePosition(k,1:3)-earthcenterpos; Dist3=temp*temp'; jiajiao=acos((Dist2+Dist3-Dist1)/2/sqrt(Dist3)/sqrt(Dist2)); if(jiajiao>=pi/2) SatellitePosition(k,4)=0; %ttum(k,:)=SatellitePosition(k,1:3); end end % ttum figure(2) drawearth(0); hold on; drawsatelliteorbit; j=1; for k=1:6 for m=1:4 if(SatellitePosition(j,4)==1) tempx=SatellitePosition(j,1); tempy=SatellitePosition(j,2); tempz=SatellitePosition(j,3); drawsatellite(tempx,tempy,tempz,k); end j=j+1; end end tempx=userposition(1); tempy=userposition(2); tempz=userposition(3); cube=100; boxplot3(tempx,tempy,tempz,cube,cube,cube,0); %标记出用户的位置 程序四（用可见卫星计算用户的位置）代码（mian3） %用余弦定理计算地心到用户和用户到卫星证件的夹角要是两者证件的夹角小于90°,就认为该卫星是不可见的. earthcenterpos=[0 0 0]; userposition=[6400 3352 5410]; %假设一个用户的位置 for k=1:24 SatellitePosition(k,:)=ttum(k,:); end for k=1:24 temp=SatellitePosition(k,1:3)-userposition; Dist1=temp*temp'; temp=userposition-earthcenterpos; Dist2=temp*temp'; temp=SatellitePosition(k,1:3)-earthcenterpos; Dist3=temp*temp'; jiajiao=acos((Dist3+Dist2-Dist1)/2/sqrt(Dist3)/sqrt(Dist2)); if(jiajiao>=pi/2) SatellitePosition(k,4)=0; end end figure(1) drawearth(0); hold on; drawsatelliteorbit; tempx=userposition(1); tempy=userposition(2); tempz=userposition(3); cube=55; boxplot3(tempx,tempy,tempz,cube,cube,cube,7); hold on; %先画出用户的位置，以便和后面的做比较 j=1; for k=1:6 for m=1:4 if(SatellitePosition(j,4)==1) tempx=SatellitePosition(j,1); tempy=SatellitePosition(j,2); tempz=SatellitePosition(j,3); drawsatellite(tempx,tempy,tempz,k); end j=j+1; end end tempx=userposition(1); tempy=userposition(2); tempz=userposition(3); hold on; SatellitePosition=[SatellitePosition;userposition 0]; %用变量caluserpos来记录迭代计算的结果. caluserpos=calculateuserposition(SatellitePosition); [m,n]=size(caluserpos); for(k=3:m); tempx=caluserpos(k,1); tempy=caluserpos(k,2); tempz=caluserpos(k,3); cube=50; boxplot3(tempx,tempy,tempz,cube,cube,cube,0); hold on; end 以下是一个和上面程序配套使用的功能函数，把它命名为calculateuserposition.m。并且要和上面的程序放在一个文件夹当中。 function caluserposition=calculateuserposition(SatellitePosition) global SatellitePosition %单独运行此程序时候，为调用前面程序而得的satelliteposotion结果，应该声明satelliteposotion为全局变量，即应该在satelliteposition前面加上global r1=6400; %地球半径 c=300000; deltat=0.00001; satelliteposnew=ones(1,3); vissatnum=0; calculateok=1; for k=1:24 if(SatellitePosition(k,4)==1) vissatnum=vissatnum+1; satelliteposnew=[satelliteposnew;SatellitePosition(k,1:3)]; end%if end%for satelliteposnew(1,:)=[]; if(vissatnum<4) calculateok=0; caluserposition=[0 0 0 0]; return end prange=ones(1,vissatnum); userpos=SatellitePosition(25,1:3); for n=1:vissatnum prange(1,n)=sqrt((satelliteposnew(n,:)-userpos)*(satelliteposnew(n,:)-userpos)')+c*deltat; end calculaterecord=[0 0 0]; xyz0=[0 0 0]; deltat0=0; wxyz=satelliteposnew; Error=1000; computertime=0; while((Error>0.00001)&(computertime<1000)) computertime=computertime+1; r=ones(1,vissatnum); for n=1:vissatnum r(1,n)=sqrt((wxyz(n,:)-xyz0)*(wxyz(n,:)-xyz0)')+deltat0*c; end deltap=r-prange; A=ones(vissatnum,3); for n=1:vissatnum A(n,:)=(wxyz(n,:)-xyz0)./r(1,n); end%for H=[A ones(vissatnum,1)]; deltax=inv(H'*H)*H'*deltap'; tempdeltax=deltax(1:3) Error=max(abs(tempdeltax)); xyz0=xyz0+deltax(1:3,:)'; if(computertime<1000) calculaterecord=[calculaterecord;xyz0]; end deltat0=deltax(4,1)/(-c); end%while calculaterecord; if(computertime==1000) caluserposition=[99999999999 99999999999 999999999999]; else caluserposition=[xyz0;calculaterecord]; end 毕业设计（论文）原创性声明和使用授权说明 原创性声明 本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得 及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。 作 者 签 名：　　 　 　　日　 期：　　 　 　　 ​​​​​​​​​​​​ 指导教师签名：　　 　 　　 日　　期：　　 　 　　 使用授权说明 本人完全了解 大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。 作者签名：　　 　 　　 日　 期：　　 　 　　 ​​​​​​​​​​​​ 学位论文原创性声明 本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。 作者签名： 日期： 年 月 日 学位论文版权使用授权书 本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权　　 　 　大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。 涉密论文按学校规定处理。 作者签名： 日期： 年 月 日 导师签名： 日期： 年 月 日 指导教师评阅书 指导教师评价： 一、撰写（设计）过程 1、学生在论文（设计）过程中的治学态度、工作精神 □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、学生掌握专业知识、技能的扎实程度 □ 优 □ 良 □ 中 □ 及格 □ 不及格 3、学生综合运用所学知识和专业技能分析和解决问题的能力 □ 优 □ 良 □ 中 □ 及格 □ 不及格 4、研究方法的科学性；技术线路的可行性；设计的合理性 □ 优 □ 良 □ 中 □ 及格 □ 不及格 5、完成毕业论文（设计）期间的出勤情况 □ 优 □ 良 □ 中 □ 及格 □ 不及格 二、论文（设计）质量 1、论文（设计）的整体结构是否符合撰写？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、是否完成指定的论文（设计）任务（包括装订及附件）？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 三、论文（设计）水平 1、论文（设计）的理论意义或对解决实际问题的指导意义 □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、论文的观念是否有新意？设计是否有创意？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 3、论文（设计说明书）所体现的整体水平 □ 优 □ 良 □ 中 □ 及格 □ 不及格 建议成绩：□ 优 □ 良 □ 中 □ 及格 □ 不及格 （在所选等级前的□内画“√”） 指导教师： （签名） 单位： （盖章） 年 月 日 评阅教师评阅书 评阅教师评价： 一、论文（设计）质量 1、论文（设计）的整体结构是否符合撰写规范？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、是否完成指定的论文（设计）任务（包括装订及附件）？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 二、论文（设计）水平 1、论文（设计）的理论意义或对解决实际问题的指导意义 □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、论文的观念是否有新意？设计是否有创意？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 3、论文（设计说明书）所体现的整体水平 □ 优 □ 良 □ 中 □ 及格 □ 不及格 建议成绩：□ 优 □ 良 □ 中 □ 及格 □ 不及格 （在所选等级前的□内画“√”） 评阅教师： （签名） 单位： （盖章） 年 月 日 教研室（或答辩小组）及教学系意见 教研室（或答辩小组）评价： 一、答辩过程 1、毕业论文（设计）的基本要点和见解的叙述情况 □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、对答辩问题的反应、理解、表达情况 □ 优 □ 良 □ 中 □ 及格 □ 不及格 3、学生答辩过程中的精神状态 □ 优 □ 良 □ 中 □ 及格 □ 不及格 二、论文（设计）质量 1、论文（设计）的整体结构是否符合撰写规范？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、是否完成指定的论文（设计）任务（包括装订及附件）？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 三、论文（设计）水平 1、论文（设计）的理论意义或对解决实际问题的指导意义 □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、论文的观念是否有新意？设计是否有创意？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 3、论文（设计说明书）所体现的整体水平 □ 优 □ 良 □ 中 □ 及格 □ 不及格 评定成绩：□ 优 □ 良 □ 中 □ 及格 □ 不及格 （在所选等级前的□内画“√”） 教研室主任（或答辩小组组长）： （签名） 年 月 日 教学系意见： 系主任： （签名） 年 月 日 学位论文原创性声明 本人郑重声明：所呈交的学位论文，是本人在导师的指导下进行的研究工作所取得的成果。尽我所知，除文中已经特别注明引用的内容和致谢的地方外，本论文不包含任何其他个人或集体已经发表或撰写过的研究成果。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式注明并表示感谢。本人完全意识到本声明的法律结果由本人承担。 学位论文作者（本人签名）： 年 月 日 学位论文出版授权书 本人及导师完全同意《中国博士学位论文全文数据库出版章程》、《中国优秀硕士学位论文全文数据库出版章程》(以下简称“章程”)，愿意将本人的学位论文提交“中国学术期刊（光盘版）电子杂志社”在《中国博士学位论文全文数据库》、《中国优秀硕士学位论文全文数据库》中全文发表和以电子、网络形式公开出版，并同意编入CNKI《中国知识资源总库》，在《中国博硕士学位论文评价数据库》中使用和在互联网上传播，同意按“章程”规定享受相关权益。 论文密级： □公开 □保密（___年__月至__年__月）(保密的学位论文在解密后应遵守此协议) 作者签名：_______ 导师签名：_______ _______年_____月_____日 _______年_____月_____日 独 创 声 明 本人郑重声明：所呈交的毕业设计(论文)，是本人在指导老师的指导下，独立进行研究工作所取得的成果，成果不存在知识产权争议。尽我所知，除文中已经注明引用的内容外，本设计（论文）不含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体均已在文中以明确方式标明。 本声明的法律后果由本人承担。   作者签名: 二〇一〇年九月二十日   毕业设计（论文）使用授权声明 本人完全了解滨州学院关于收集、保存、使用毕业设计（论文）的规定。 本人愿意按照学校要求提交学位论文的印刷本和电子版，同意学校保存学位论文的印刷本和电子版，或采用影印、数字化或其它复制手段保存设计（论文）；同意学校在不以营利为目的的前提下，建立目录检索与阅览服务系统，公布设计（论文）的部分或全部内容，允许他人依法合理使用。 （保密论文在解密后遵守此规定）   作者签名: 二〇一〇年九月二十日 致 谢 时间飞逝，大学的学习生活很快就要过去，在这四年的学习生活中，收获了很多，而这些成绩的取得是和一直关心帮助我的人分不开的。 首先非常感谢学校开设这个课题，为本人日后从事计算机方面的工作提供了经验，奠定了基础。本次毕业设计大概持续了半年，现在终于到结尾了。本次毕业设计是对我大学四年学习下来最好的检验。经过这次毕业设计，我的能力有了很大的提高，比如操作能力、分析问题的能力、合作精神、严谨的工作作风等方方面面都有很大的进步。这期间凝聚了很多人的心血，在此我表示由衷的感谢。没有他们的帮助，我将无法顺利完成这次设计。 首先，我要特别感谢我的知道郭谦功老师对我的悉心指导，在我的论文书写及设计过程中给了我大量的帮助和指导，为我理清了设计思路和操作方法，并对我所做的课题提出了有效的改进方案。郭谦功老师渊博的知识、严谨的作风和诲人不倦的态度给我留下了深刻的印象。从他身上，我学到了许多能受益终生的东西。再次对周巍老师表示衷心的感谢。 其次，我要感谢大学四年中所有的任课老师和辅导员在学习期间对我的严格要求，感谢他们对我学习上和生活上的帮助，使我了解了许多专业知识和为人的道理，能够在今后的生活道路上有继续奋斗的力量。 另外，我还要感谢大学四年和我一起走过的同学朋友对我的关心与支持，与他们一起学习、生活，让我在大学期间生活的很充实，给我留下了很多难忘的回忆。 最后，我要感谢我的父母对我的关系和理解，如果没有他们在我的学习生涯中的无私奉献和默默支持，我将无法顺利完成今天的学业。 四年的大学生活就快走入尾声，我们的校园生活就要划上句号，心中是无尽的难舍与眷恋。从这里走出，对我的人生来说，将是踏上一个新的征程，要把所学的知识应用到实际工作中去。 回首四年，取得了些许成绩，生活中有快乐也有艰辛。感谢老师四年来对我孜孜不倦的教诲，对我成长的关心和爱护。 学友情深，情同兄妹。四年的风风雨雨，我们一同走过，充满着关爱，给我留下了值得珍藏的最美好的记忆。 在我的十几年求学历程里，离不开父母的鼓励和支持，是他们辛勤的劳作，无私的付出，为我创造良好的学习条件，我才能顺利完成完成学业，感激他们一直以来对我的抚养与培育。 最后，我要特别感谢我的导师***老师、和研究生助教***老师。是他们在我毕业的最后关头给了我们巨大的帮助与鼓励，给了我很多解决问题的思路，在此表示衷心的感激。老师们认真负责的工作态度，严谨的治学精神和深厚的理论水平都使我收益匪浅。他无论在理论上还是在实践中，都给与我很大的帮助，使我得到不少的提高这对于我以后的工作和学习都有一种巨大的帮助，感谢他耐心的辅导。在论文的撰写过程中老师们给予我很大的帮助，帮助解决了不少的难点，使得论文能够及时完成，这里一并表示真诚的感谢。 致 谢 这次论文的完成，不止是我自己的努力，同时也有老师的指导，同学的帮助，以及那些无私奉献的前辈，正所谓你知道的越多的时候你才发现你知道的越少，通过这次论文，我想我成长了很多，不只是磨练了我的知识厚度，也使我更加确定了我今后的目标：为今后的计算机事业奋斗。在此我要感谢我的指导老师——***老师，感谢您的指导，才让我有了今天这篇论文，您不仅是我的论文导师，也是我人生的导师，谢谢您！我还要感谢我的同学，四年的相处，虽然我未必记得住每分每秒，但是我记得每一个有你们的精彩瞬间，我相信通过大学的历练，我们都已经长大，变成一个有担当，有能力的新时代青年，感谢你们的陪伴，感谢有你们，这篇论文也有你们的功劳，我想毕业不是我们的相处的结束，它是我们更好相处的开头，祝福你们！我也要感谢父母，这是他们给我的，所有的一切；感谢母校，尽管您不以我为荣，但我一直会以我是一名农大人为荣。 通过这次毕业设计，我学习了很多新知识，也对很多以前的东西有了更深的记忆与理解。漫漫求学路，过程很快乐。我要感谢信息与管理科学学院的老师，我从他们那里学到了许多珍贵的知识和做人处事的道理，以及科学严谨的学术态度，令我受益良多。同时还要感谢学院给了我一个可以认真学习，天天向上的学习环境和机会。 即将结束*大学习生活，我感谢****大学提供了一次在农大接受教育的机会，感谢院校老师的无私教导。感谢各位老师审阅我的论文。 毕业设计（论文）原创性声明和使用授权说明 原创性声明 本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得 及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。 作 者 签 名：　　 　 　　日　 期：　　 　 　　 ​​​​​​​​​​​​ 指导教师签名：　　 　 　　 日　　期：　　 　 　　 使用授权说明 本人完全了解 大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。 作者签名：　　 　 　　 日　 期：　　 　 　　 ​​​​​​​​​​​​ 学位论文原创性声明 本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。 作者签名： 日期： 年 月 日 学位论文版权使用授权书 本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权　　 　 　大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。 涉密论文按学校规定处理。 作者签名： 日期： 年 月 日 导师签名： 日期： 年 月 日 独 创 声 明 本人郑重声明：所呈交的毕业设计(论文)，是本人在指导老师的指导下，独立进行研究工作所取得的成果，成果不存在知识产权争议。尽我所知，除文中已经注明引用的内容外，本设计（论文）不含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体均已在文中以明确方式标明。 本声明的法律后果由本人承担。   作者签名: 年 月 日   毕业设计（论文）使用授权声明 本人完全了解**学院关于收集、保存、使用毕业设计（论文）的规定。 本人愿意按照学校要求提交学位论文的印刷本和电子版，同意学校保存学位论文的印刷本和电子版，或采用影印、数字化或其它复制手段保存设计（论文）；同意学校在不以营利为目的的前提下，建立目录检索与阅览服务系统，公布设计（论文）的部分或全部内容，允许他人依法合理使用。 （保密论文在解密后遵守此规定）   作者签名: 年 月 日 基本要求：写毕业论文主要目的是培养学生综合运用所学知识和技能，理论联系实际，独立分析，解决实际问题的能力，使学生得到从事本专业工作和进行相关的基本训练。毕业论文应反映出作者能够准确地掌握所学的专业基础知识，基本学会综合运用所学知识进行科学研究的方法，对所研究的题目有一定的，论文题目的范围不宜过宽，一般选择本学科某一重要问题的一个侧面。 毕业论文的基本教学要求是： 1、培养学生综合运用、巩固与扩展所学的基础理论和专业知识，培养学生独立分析、解决实际问题能力、培养学生处理数据和信息的能力。2、培养学生正确的理论联系实际的工作作风，严肃认真的科学态度。3、培养学生进行社会调查研究；文献资料收集、阅读和整理、使用；提出论点、综合论证、总结写作等基本技能。 毕业论文是毕业生总结性的独立作业，是学生运用在校学习的基本知识和基础理论，去分析、解决一两个实际问题的实践锻炼过程，也是学生在校学习期间学习成果的综合性总结，是整个教学活动中不可缺少的重要环节。撰写毕业论文对于培养学生初步的科学研究能力，提高其综合运用所学知识分析问题、解决问题能力有着重要意义。 毕业论文在进行编写的过程中，需要经过开题报告、论文编写、论文上交评定、论文答辩以及论文评分五个过程，其中开题报告是论文进行的最重要的一个过程，也是论文能否进行的一个重要指标。 撰写意义：1.撰写毕业论文是检验学生在校学习成果的重要措施，也是提高教学质量的重要环节。大学生在毕业前都必须完成毕业论文的撰写任务。申请学位必须提交相应的学位论文，经答辩通过后，方可取得学位。可以这么说，毕业论文是结束大学学习生活走向社会的一个中介和桥梁。毕业论文是大学生才华的第一次显露，是向祖国和人民所交的一份有份量的答卷，是投身社会主义现代化建设事业的报到书。一篇毕业论文虽然不能全面地反映出一个人的才华，也不一定能对社会直接带来巨大的效益，对专业产生开拓性的影响。但是，实践证明，撰写毕业论文是提高教学质量的重要环节，是保证出好人才的重要措施。 2.通过撰写毕业论文，提高写作水平是干部队伍“四化”建设的需要。党中央要求，为了适应现代化建设的需要，领导班子成员应当逐步实现“革命化、年轻化、知识化、专业化”。这个“四化”的要求，也包含了对干部写作能力和写作水平的要求。 3.提高大学生的写作水平是社会主义物质文明和精神文明建设的需要。在新的历史时期，无论是提高全族的科学文化水平，掌握现代科技知识和科学管理方法，还是培养社会主义新人，都要求我们的干部具有较高的写作能力。在经济建设中，作为领导人员和机关的办事人员，要写指示、、总结、调查报告等应用文；要写说明书、广告、解说词等说明文；还要写科学论文、经济评论等议论文。在当今信息社会中，信息对于加快经济发展速度，取得良好的经济效益发挥着愈来愈大的作用。写作是以语言文字为信号，是传达信息的方式。信息的来源、信息的收集、信息的储存、整理、传播等等都离不开写作。 论文种类：毕业论文是学术论文的一种形式，为了进一步探讨和掌握毕业论文的写作规律和特点，需要对毕业论文进行分类。由于毕业论文本身的内容和性质不同，研究领域、对象、方法、表现方式不同，因此，毕业论文就有不同的分类方法。 按内容性质和研究方法的不同可以把毕业论文分为理论性论文、实验性论文、描述性论文和设计性论文。后三种论文主要是理工科大学生可以选择的论文形式，这里不作介绍。文科大学生一般写的是理论性论文。理论性论文具体又可分成两种：一种是以纯粹的抽象理论为研究对象，研究方法是严密的理论推导和数学运算，有的也涉及实验与观测，用以验证论点的正确性。另一种是以对客观事物和现象的调查、考察所得观测资料以及有关文献资料数据为研究对象，研究方法是对有关资料进行分析、综合、概括、抽象，通过归纳、演绎、类比，提出某种新的理论和新的见解。 按议论的性质不同可以把毕业论文分为立论文和驳论文。立论性的毕业论文是指从正面阐述论证自己的观点和主张。一篇论文侧重于以立论为主，就属于立论性论文。立论文要求论点鲜明，论据充分，论证严密，以理和事实服人。驳论性毕业论文是指通过反驳别人的论点来树立自己的论点和主张。如果毕业论文侧重于以驳论为主，批驳某些错误的观点、见解、理论，就属于驳论性毕业论文。驳论文除按立论文对论点、论据、论证的要求以外，还要求针锋相对，据理力争。 按研究问题的大小不同可以把毕业论文分为宏观论文和微观论文。凡届国家全局性、带有普遍性并对局部工作有一定指导意义的论文，称为宏观论文。它研究的面比较宽广，具有较大范围的影响。反之，研究局部性、具体问题的论文，是微观论文。它对具体工作有指导意义，影响的面窄一些。 另外还有一种综合型的分类方法，即把毕业论文分为专题型、论辩型、综述型和综合型四大类： 1．专题型论文。这是分析前人研究成果的基础上，以直接论述的形式发表见解，从正面提出某学科中某一学术问题的一种论文。如本书第十二章例文中的《浅析领导者突出工作重点的方法与艺术》一文，从正面论述了突出重点的工作方法的意义、方法和原则，它表明了作者对突出工作重点方法的肯定和理解。2．论辩型论文。这是针对他人在某学科中某一学术问题的见解，凭借充分的论据，着重揭露其不足或错误之处，通过论辩形式来发表见解的一种论文。3．综述型论文。这是在归纳、总结前人或今人对某学科中某一学术问题已有研究成果的基础上，加以介绍或评论，从而发表自己见解的一种论文。4．综合型论文。这是一种将综述型和论辩型两种形式有机结合起来写成的一种论文。如《关于中国民族关系史上的几个问题》一文既介绍了研究民族关系史的现状，又提出了几个值得研究的问题。因此，它是一篇综合型的论文。 写作步骤：毕业论文是高等教育自学考试本科专业应考者完成本科阶段学业的最后一个环节，它是应考者的 总结 性独立作业，目的在于总结学习专业的成果，培养综合运用所学知识解决实际 问题 的能力。从文体而言，它也是对某一专业领域的现实问题或 理论 问题进行 科学 研究 探索的具有一定意义的论说文。完成毕业论文的撰写可以分两个步骤，即选择课题和研究课题。 首先是选择课题。选题是论文撰写成败的关键。因为，选题是毕业论文撰写的第一步，它实际上就是确定“写什么”的问题，亦即确定科学研究的方向。如果“写什么”不明确，“怎么写”就无从谈起。 教育部自学考试办公室有关对毕业论文选题的途径和要求是“为鼓励理论与工作实践结合，应考者可结合本单位或本人从事的工作提出论文题目，报主考学校审查同意后确立。也可由主考学校公布论文题目，由应考者选择。毕业论文的总体要求应与普通全日制高等学校相一致，做到通过论文写作和答辩考核，检验应考者综合运用专业知识的能力”。但不管考生是自己任意选择课题，还是在主考院校公布的指定课题中选择课题，都要坚持选择有科学价值和现实意义的、切实可行的课题。选好课题是毕业论文成功的一半。 第一、要坚持选择有科学价值和现实意义的课题。科学研究的目的是为了更好地认识世界、改造世界，以推动社会的不断进步和发展 。因此，毕业论文的选题，必须紧密结合社会主义物质文明和精神文明建设的需要，以促进科学事业发展和解决现实存在问题作为出发点和落脚点。选题要符合科学研究的正确方向，要具有新颖性，有创新、有理论价值和现实的指导意义或推动作用，一项毫无意义的研究，即使花很大的精力，表达再完善，也将没有丝毫价值。具体地说，考生可从以下三个方面来选题。首先，要从现实的弊端中选题，学习了专业知识，不能仅停留在书本上和理论上，还要下一番功夫，理论联系实际，用已掌握的专业知识，去寻找和解决工作实践中急待解决的问题。其次，要从寻找科学研究的空白处和边缘领域中选题，科学研究。还有许多没有被开垦的处女地，还有许多缺陷和空白，这些都需要填补。应考者应有独特的眼光和超前的意识去思索，去发现，去研究。最后，要从寻找前人研究的不足处和错误处选题，在前人已提出来的研究课题中，许多虽已有初步的研究成果，但随着社会的不断发展，还有待于丰富、完整和发展，这种补充性或纠正性的研究课题，也是有科学价值和现实指导意义的。 第二、要根据自己的能力选择切实可行的课题。毕业论文的写作是一种创造性劳动，不但要有考生个人的见解和主张，同时还需要具备一定的客观条件。由于考生个人的主观、客观条件都是各不相同的，因此在选题时，还应结合自己的特长、兴趣及所具备的客观条件来选题。具体地说，考生可从以下三个方面来综合考虑。首先，要有充足的资料来源。“巧妇难为无米之炊”，在缺少资料的情况下，是很难写出高质量的论文的。选择一个具有丰富资料来源的课题，对课题深入研究与开展很有帮助。其次，要有浓厚的研究兴趣，选择自己感兴趣的课题，可以激发自己研究的热情，调动自己的主动性和积极性，能够以专心、细心、恒心和耐心的积极心态去完成。最后，要能结合发挥自己的业务专长，每个考生无论能力水平高低，工作岗位如何，都有自己的业务专长，选择那些能结合自己工作、发挥自己业务专长的课题，对顺利完成课题的研究大有益处。 致 谢 这次论文的完成，不止是我自己的努力，同时也有老师的指导，同学的帮助，以及那些无私奉献的前辈，正所谓你知道的越多的时候你才发现你知道的越少，通过这次论文，我想我成长了很多，不只是磨练了我的知识厚度，也使我更加确定了我今后的目标：为今后的计算机事业奋斗。在此我要感谢我的指导老师——***老师，感谢您的指导，才让我有了今天这篇论文，您不仅是我的论文导师，也是我人生的导师，谢谢您！我还要感谢我的同学，四年的相处，虽然我未必记得住每分每秒，但是我记得每一个有你们的精彩瞬间，我相信通过大学的历练，我们都已经长大，变成一个有担当，有能力的新时代青年，感谢你们的陪伴，感谢有你们，这篇论文也有你们的功劳，我想毕业不是我们的相处的结束，它是我们更好相处的开头，祝福你们！我也要感谢父母，这是他们给我的，所有的一切；感谢母校，尽管您不以我为荣，但我一直会以我是一名农大人为荣。 通过这次毕业设计，我学习了很多新知识，也对很多以前的东西有了更深的记忆与理解。漫漫求学路，过程很快乐。我要感谢信息与管理科学学院的老师，我从他们那里学到了许多珍贵的知识和做人处事的道理，以及科学严谨的学术态度，令我受益良多。同时还要感谢学院给了我一个可以认真学习，天天向上的学习环境和机会。 即将结束*大学习生活，我感谢****大学提供了一次在**大接受教育的机会，感谢院校老师的无私教导。感谢各位老师审阅我的论文。 � EMBED Equation.3 ��� Sj(t0) � EMBED Equation.3 ��� N0 k N0 Int(φ) � EMBED Equation.3 ��� 通常的相位测量或相位差测量只是测出一周以内的相位值，实际测量中，如果对整周进行计数，则自某一初始取样时刻（t0）以后就可以取得连续的相位观测值。 Sj(ti) 平天球坐标 岁差，章动影响 真天球坐标 旋转春分点时角Ω 真地球坐标 极移旋转 平地球坐标 fs 图4-1 卫星的椭圆运行轨道 bs M as 图4-2 卫星地心向径在相同的时间间隔内扫描的面积 图5-1 可见卫星的测量原理 E 飞机 h Rn 地球 质心 用户位置 地心位置 PAGE 40 _1301236348.unknown _1301295950.unknown _1301831481.unknown _1301832108.unknown _1301943362.unknown _1305374614.unknown _1305290151.unknown _1301832245.unknown _1301831984.unknown _1301297078.unknown _1301313864.unknown _1301313937.unknown _1301831479.unknown _1301337102.unknown _1301313924.unknown _1301313846.unknown _1301313849.unknown _1301297253.unknown _1301296371.unknown _1301296759.unknown _1301296907.unknown _1301296056.unknown _1301296167.unknown _1301294722.unknown _1301295917.unknown _1301295925.unknown _1301295666.unknown _1301295756.unknown _1301294997.unknown _1301295660.unknown _1301294989.unknown _1301294264.unknown _1301294558.unknown _1301294628.unknown _1301294302.unknown _1301292403.unknown _1301292708.unknown _1301292709.unknown _1301292450.unknown _1301292460.unknown _1301292471.unknown _1301292428.unknown _1301292377.unknown _1301235867.unknown _1301236094.unknown _1301236101.unknown _1301235905.unknown _1301235406.unknown _1301235787.unknown _1301235827.unknown _1301235473.unknown _1301235543.unknown _1300812138.unknown _1300819049.unknown _1300819073.unknown _1300820722.unknown _1300812381.unknown _1300812121.unknown _1300812126.unknown _1300629931.unknown _1300629548.unknown