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

1、GPS卫星运动及定位matlab仿真作者姓名：专业班级：指导老师摘要全球定位系统是具有全球性、全能性、全天候优势的导航定位、定时和测速系统，现在在全球很多领域获得了应用。GPS卫星的定位是一个比较复杂的系统，其包含参数众多，如时间系统、空间坐标系统等。此次设计是针对卫星运动定位的matlab仿真实现，因要求不高，所以对卫星运动做了理想化处理,摄动力对卫星的影响忽略不计（所以为无摄运动），采用开普勒定律及最小二乘法计算其轨道参数,对其运动规律进行简略分析，并使用matlab编程仿真实现了卫星的运功轨道平面、运动动态、可见卫星的分布及利用可见卫星计算出用户位置。通过此次设计，对于GPS卫星有了初步

2、的认识，对于静态单点定位、伪距等相关概念有一定了解。关键字：GPS卫星无摄运动伪距matlab仿真ThemovementandlocationofGPSsatelliteonMATLABAbstract：Globalpositioningsystemisaglobal,versatility,all-weatheradvantageofnavigationandpositioning,timingandspeedsystem,nowtherehasmanaypplicationinmanyfields.GPSsatellitepositioningisacomplexsystem,whichi

3、ncludesmanyparameters,suchastimeandspacecoordinatessystem.Thisdesignisbasedonthematlabsimulationofsatellitemotionandlocation,becausedemandisnothigh,sotodotheidealizedsatellitemovement,andignorethedisturbedmotion(socallitnon-disturbedmotion).UsingtheKeplerandleast-squaremethodforcalculatingtheparamet

4、ersoforbitalmotion,forthecharacteristicsofmotiontomakeasimpleanalysis,andusethematlabsimulationtoprogramachievetheorbitalplaneofsatellite,thedynamicmotion,thedistributionofvisiblesatellitesandusingvisiblesatellitestocalculatetheusers'home.ThroughthedesignhaveprimaryunderstandingfortheGPSsatellit

5、e,andunderstandingthestaticsingle-point,pseudorangeandsoon.Keywords:GPSsatellitenon-disturbedmotionpseudorangematlabsimulation目录第一章前言11.1课题背景11.2本课题研究的意义和方法21.3GPS前景2第二章GPS测量原理42.1伪距测量的原理42.1.1 计算卫星位置52.1.2 用户位置的计算52.1.3最小二乘法介绍52.2载波相位测量原理6第三章GPS的坐标、时间系统103.1坐标系统103.1.1天球坐标系103.1.2地球坐标系123.2时间系统133.

6、2.1世界时系统143.2.2原子时系统153.2.3动力学时系统163.2.4协调世界时163.2.5GPS时间系统16第四章卫星运动基本定律及其求解184.1开普勒第一定律184.2开普勒第二定律194.3开普勒第三定律204.4卫星的无摄运动参数204.5真近点角的概念及其求解214.6卫星瞬时位置的求解22第五章GPS的MATLAB仿真255.1卫星可见性的估算255.2 GPS卫星运动的MATLAB仿真26结论38致谢40参考文献41附录42第一章前言1.1 课题背景GPS系统的前身为美军研制的一种子午仪卫星定位系统(Transit),1958年研制，64年正式投入使用。该系统用5到

7、6颗卫星组成的星网工作，每天最多绕过地球13次，并且无法给出高度信息，在定位精度方面也不尽如人意。然而，子午仪系统使得研发部门对卫星定位取得了初步的经验，并验证了由卫星系统进行定位的可行性，为GPS系统的研制埋下了铺垫。由于卫星定位显示出在导航方面的巨大优越性及子午仪系统存在对潜艇和舰船导航方面的巨大缺陷。美国海陆空三军及民用部门都感到迫切需要一种新的卫星导航系统13。1973年12月，美国国防部批准它的陆海空三军联合研制新的卫星导航系统:NAVSTAR/GPS。它是英文“NavigationSatelliteTimingandRanging/GlobalPositioningSystem”的

8、缩写词。其意为“卫星测时测距导航/全球定位系统”,简称GPS。这个系统向有适当接受设备的全球范围用户提供精确、连续的三维位置和速度信息，并且还广播一种形式的世界协调时(UTC)。通过遍布全球的(21+3)GPS导航卫星，向全球范围内的用户全天候提供高精度的导航、跟踪定位和授时服务。目前，GPS已在地形测量，交通管理，导航，野外勘探，空间宇宙学等诸多领域得到了广泛的应用11。目前全球共有4大GPS系统，分别是：美国GPS，由美国国防部于20世纪70年代初开始设计、研制，于1993年全部建成。1994年，美国宣布在10年内向全世界免费提供GPS使用权，但美国只向外国提供低精度的卫星信号。欧盟“伽利

9、略”，1999年欧洲提出计划，准备发射30颗卫星，组成“伽利略”卫星定位系统。俄罗斯“格洛纳斯”，尚未部署完毕。始于上世纪70年代，需要至少18颗卫星才能确保覆盖俄罗斯全境;如要提供全球定位服务，则需要24颗卫星。中国“北斗”2003年我国北斗一号建成并开通运行，不同于GPS，“北斗”的指挥机和终端之间可以双向交流。四川大地震发生后，北京武警指挥中心和四川武警部队运用“北斗”进行了上百次交流。北斗二号系列卫星今年起将进入组网高峰期，预计在2015年形成由三十几颗卫星组成的覆盖全球的系统。1.2 本课题研究的意义和方法GPS系统是一个很庞大的系统，包含了天文，地理，计算机，电磁学，通信学，信息学

10、等等。通过本文对GPS的学习研究，最重要的还是要学习其原理：卫星运动原理；卫星定位原理；卫星跟踪原理等等。通过基础原理的学习，一方面，可以使我们更进一步的理解卫星运动，定位的实现方法；通过仿真，进一步了解简单定位的方法及其在仿真平台上的实现途径；另一方面，也可以培养我们自学的能力，训练仿真模拟的技巧和方法。至今，基本上完成了课题的要求，通过不断的注入既定参数，可以更加详细，直观的理解基本的定位原理和实现方法！1.3 GPS前景GPS导航定位以其定位精度高、观测时间短、测站间无需通视、可提供三维坐标、操作简便、全天候作业、功能多、应用广泛等特点着称。用GPS信号可以进行海、空和陆地的导航、导弹的

11、制导、大地测量和工程测量的精密定位、时间的传递和速度的测量等。对于测绘领域，GPS卫星定位技术已经用于建立高精度的全国性的大地测量控制网，测定全球性的地球动态参数；用于建立陆地海洋大地测量基准，进行高精度的海岛陆地联测以及海洋测绘；用于检测地球板块运动状态和地壳形变；用于工程测量，成为建立城市与工程控制网的主要手段。用于测定航空航天摄影瞬间相机位置，实现仅有少量的地面控制或无地面控制的航测快速成图，导致地理信息系统、全球环境遥感监测的技术革命4。目前，GPS、GLONASS、INMARSAT等系统都具备了导航定位功能，形成了多元化的空间资源环境。这一多元化的空间资源环境，促使国际民间形成了一个

12、共同的策略，即一方面对现有系统充分利用，一方面积极筹建民间GNSS系统，待2011年左右，GNSS纯民间系统建成，全球将形成GPS/GLONASS/GNSS三足鼎立之势，才能从根本上摆脱对单一系统的依赖，形成国际共有、国际共享的安全资源环境。世界才可以将卫星导航作为单一导航手段的最高应用境界。国际民间的这一策略，反过来又影响和迫使美国对其GPS使用政策作出更开放的调整。多元化的空间资源环境的确立，给GPS的发发展应用创造了一个前所未有的良好的国际环境。第二章GPS测量原理GPS导航系统的基本原理是测量出已知位置的卫星到用户接收机之间的距离，然后综合多颗卫星的数据就可知道接收机的具体位置。要达到

13、这一目的，卫星的位置可以根据星载时钟所记录的时间在卫星星历中查出。而用户到卫星的距离则通过纪录卫星信号传播到用户所经历的时间，再将其乘以光速得到(由于大气层电离层的干扰，这一距离并不是用户与卫星之间的真实距离，而是伪距(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导航

14、系统卫星部分的作用就是不断地发射导航电文。然而，由于用户接受机使用的时钟与卫星星载时钟不可能总是同步，所以除了用户的三维坐标X、y、z外，还要引进一个At即卫星与接收机之间的时间差作为未知数，然后用4个方程将这4个未知数解出来。所以如果想知道接收机所处的位置，至少要能接收到4个卫星的信号。2.1伪距测量的原理GPS定位采用的是被动式单程测距。它的信号发射书机由卫星钟确定，收到时刻是由接收机钟确定，这就在测定的卫星至接收机的距离中，不可避免地包含着两台钟不同步的误差和电离层、对流层延迟误差影响，它并不是卫星与接受机之间的实际距离，所以称之为伪距。伪距定位法是利用全球卫星定位系统进行导航定位的最基

15、本的方法，其基本原理是：在某一瞬间利用GPS接收机同时测定至少四颗卫星的伪距，根据已知的卫星位置和伪距观测值，采用距离交会法求出接收机的三维坐标和时钟改正数。它的优点是速度快、无多值性问题，利用增加观测时间可以提高定位精度；缺点是测量定位精度低，但足以满足部分用户的需要。2.1.1计算卫星位置读入导航电文后首先根据需要调用广播轨道1至广播轨道5上的数据然后依次计算卫星的平均角速度归化时间平均近点角需要注意的是进行真近点角计算时要同时计算正弦和余弦以得到正确象限内的角计算经校正的升交点精度时需要用到地球旋转速率在WGS-84中这一常数为a:。（2-1）2.1.2用户位置的计算首先利用近似的用户位

16、置与伪距观测值计算出一个近似伪距利用该近似伪距可以计算出部分值然后计算出系数并生成一个Nx4的矩阵为参与运算的卫星数最后算出用户位移的坐标上述过程根据需要可以将计算出的用户坐标作为近似值反复迭代直至符合精度要求a。当在进行迭代的过程当中，如果所给定的用户的初始值越接近用户的实际值，则迭代的次数就越少。当我们可见的卫星多于四颗的时候，我们可以用以下介绍的最小二乘法原理带到上面的公式当中去计算。2.1.3最小二乘法介绍当我们在一个地方同时可见的卫星如果多于四颗（GPS卫星的轨道设计和运动时间的安排使得用户在地球的任意位置（两极个别地点除外），都能够看到411颗的卫星）的时候，我们可以用最小二乘法去

17、解算未知数，这样，充分的利用了已知的数据信息，使得结果的偏差最小化。例如：对于下面的方程：alx+bly=cl<a2x+b2y=c2a3x+b3y=c3(2-2)如果令厂alH=a2a3bOb2b3I(2-3)使用最小二乘法，用C,H表示X.令alx+bly-cl)"2+(a2x+b2y-c2厂2+(a3x+b3y-c3厂2(2-4)dx=2(alx+bly-cl)al+2(a2x+b2y-c2)a2+2(a3x+b3y-c3)a3(2-5)2(alx+bly-cl)bl+2(a2x+b2y-c2)b2+2(a3x+b3y-c3)b3(2-6)整理得到(a2+a2“2+a3“2

18、)x+(albl+a2b2+a3b3)y二alcl+a2c2+a3c3(2-7)(a2+a2“2+a3“2)y+(albl+a2b2+a3b3)x二blcl+b2c2+b3c3(2-8)写成矩阵的形式就是:厂albl厂albl哪么就能够得到HtH-1Htc(2-9)2.2载波相位测量原理载波相位观测方法是GPS接收机用接收到的卫星载波量的时候，先要进行解调工作，把调制在载波上面的测距码和导航电文去掉，通过码相关等方法重新获取载波以®k"(Ll:154f0,19.032cm;L2:120f0,24.42cm)与本地接收机产生的本振参考载波产生的相位差来计算的.(GPS所接收到

19、得载波相位是不连续的，所以在进行相位测表示k接收机在时刻tk所接收到的第j颗卫星接收到的载波相位的值；以®、丿表示k接收机在时刻tk本地载波信号的相位值，则接收机在接收机钟面时刻tk时观测j卫星所取得的相位观测量可写为7j(t)=申j(t)-申(t)图2-1载波相位测量原理图N0j，此时包含整周数的相位观测值为：kA90+N0j=9kGt0)-9kGt0)+N0j如上图1所示，在初始时刻to,测得小于一周的相位差为A®0,其整周数为(2-11)接收机继续跟踪卫星信号，不断地测量小于一周的相位差A9(t)，并利用整波计数器记录从to到ti时间内的整周数变化量INT(申)，只要

20、卫星j从to到tj时间内信号没有中断，则整周模糊数NOj就为一个常数，任意时刻ti卫星到INT(申)(2-12)k接收机的相位差为：G)=屮jG)一屮c)+n0j+ikiki0这样，观测量就包含了相位差的小数部分和累计的整数部分的整周数。载波相位的观测方程假设在GPS系统时刻Ta(卫星a时刻)卫星Sj发射的载波信号相位为载波相位的观测值为，考虑到卫星与接收机和系统时间的差值，Ta=ta+收机)。接收机产生的本地载波相位为9(tb)，根据9(t)，经过传播的延迟A后，在GPS系统时刻Tb(接收机tb时刻到达接d二e(t)-ej(t)ba(2-13)(210)得到：在Tb时刻,5t,Tb=tb+t

21、,ab则有：=e(T-51)j(T-51)bbaa(2-14)由于卫星和接收机的频率都比较稳定，所以在一个小的时间间隔里面，我们可以近似的理解在时间t,t+At内，有e(t+At)(t)+fAt(上式是考虑在频率前提下，所以没有在f前面乘以2兀)(2-15)，因为卫星到接收机有一个传播的延迟，即Tb=Ta+at,所以有：e(Tb)(Ta)+fA(2-16)由2-15和2-16代入到2-14得到：2-17)，如果同时考虑到传播延迟电离层和对流层的影响(5P和5P)，有p-5-512Hppc，p为卫星到接收机的距离。则有：-5p-5p)+-2f5ta12(2-18)将2-12代入上个式子，考虑到九

22、=，得到以米为单位的测量的载波相位为:j-(p25kP"5P)+:(ft-51)+NjAabk(2-19)这样，上式即为接收机k对卫星k的载波相位的以米为单位的观测方程式。第三章GPS的坐标、时间系统3.1坐标系统GPS定位测量当中，要用到两种坐标系，即天球坐标系和地球坐标系。天球坐标系是指坐标原点和各坐标轴的指向在空间是保持不变的，可以很方便的描述卫星的运动和状态。而地球坐标系则是与地球体相关联的坐标系，用于描述地面测量站的位置。下面就天球坐标系和地球坐标系做简要的说明。3.1.1天球坐标系天球就是指的是以地球质心为中心，半径无穷大的理想球体。在这个系统当中，我们会涉及到几个参考点

23、，线，面9。1：天轴和天极：天轴是指地球自转的延伸直线，天轴和地球表面的交点叫做天极P,与地球北极相对应的是北天极Pn,与地球南极相对应的是南天极Ps.天极并不是固定的，有岁差和章动的影响，这个时候叫做真天极，而无岁差和章动影像的天极叫做平天极。2：天球赤道：通过地球质心并与天轴垂直的平面与地球表面的交线叫做天球赤道。3：天球子午面：包含天轴并通过天球面上任意一点的平面。4：黄道：地球绕太阳公转的轨道平面和天球表面相交的大圆，黄道平面和天球赤道平面的夹角叫做黄赤交角，为23.5度。5：黄极：过天球中心垂直于黄道平面的直线与天球表面相交的点，它分为黄北极和黄南极，分别用Kn和Ks表示。6：春分点

24、：指太阳由南向北运动的时候，所经过的天球黄道和天球赤道的交点。春分点和天球赤道面是建立天球坐标系的基准点和基准面。图3-1春分点7：岁差和章动：岁差指的是平北天极以北黄极为中心，以黄赤交角为半径的一种顺时针圆周运动。章动指的是真北天极绕平北天极做得椭圆型运动。图3-2岁差和章动天球坐标系分为两种：真天球坐标系和平天球坐标系.真天球坐标系的原点为地球的质心M，Z轴指向真北天极Pn,X轴指向春分点,Y轴垂直于XMZ平面。平天球坐标系的原点为地球的质心M,Z轴指向平北天极Pn,X轴指向春分点,Y轴垂直于XMZ平面。上述两种坐标系的差别在于他们选取了不同的北天极的位置，故要是由平天极坐标系到真天极坐标

25、系的转换，就必须考虑岁差和章动旋转所影响的情况。换句话说就是要考虑到岁差旋转和章动旋转地影响。3.1.2地球坐标系地球坐标系也可以分为两种：即平地球坐标系和真地球坐标系。1平地球坐标系：它的地极位置采用国际协议地极原点CIO(由1900到1905年测定的平均纬度所确定的平均地极位置)。原点：地球质心M。Z轴：指向CIO。X轴：指向格林威治起始子午面与地球平赤道的交点。Y轴：垂直于XMZ平面。2真地球坐标系原点：地球质心M。Z轴：指向地球瞬时极。X轴：指向格林威治起始子午面与地球瞬时真赤道的交点。Y轴：垂直于XMZ平面。图3-3世界地心坐标系瞬时真天球坐标到瞬时真地球坐标的转换相应的转换转动矩阵

26、为：cOsQ01这两种坐标的差异就在于X轴的指向是不同的。前者指向的是真春分点，而后者指向的是格林威治起始子午面与地球瞬时真赤道的交点。两者之间的夹角称为对应的平格林威治起始子午面的真春分点时角Q。故仅仅需要绕Z轴旋转这个角度Q,就能够做到二者的相互转换。CosQ-SinQRz(Q)=SinQ3-100综合上面的，可以得到以下的结论：在GPS定位系统所用的空间坐标系统当中，我们一般采用天球坐标去研究卫星的空间运动，而采用地球坐标去研究地面监控站点，他们之间的转换问题一般可以按照下面的步骤来分析：旋转春分点时角岁差，章动影响极移旋转3.2时间系统时间系统是卫星定位测量过程中的一个重要概念。现时的

27、GPS测量的方法是通过接收和处理GPS卫星发射的无线电信号，以确定用户接收机和观测卫星间的距离，然后通过一定的数学方法以确定接收机所在的具体位置，为得到接收机和卫星的准确距离，必须获得无线电信号从卫星传输至接收机这一过程中的精确时间，因而利用卫星技术进行精密的定位和导航，必须要获得高精度的时间信息，这需要一个精确的时间系统。现行的卫星定位测量中与之紧密相关的时间系统有三种:世界时，原子时和动力学时。3.2.1世界时系统以地球自转为基准的一种时间系统。根据不同的空间参考点，又可分为恒星时，太阳时，世界时三种。1.恒星时选定春分点(地球赤道平面与其绕太阳公转轨道的一个交点)作为参考点，由该点的周日

28、视运动所确定的时间，即为恒星时(siderealTime，sT)。规定从春分点连续两次经过本地子午圈的时间间隔为一恒星日，其1/24为一恒星时，由于其定义涉及到地方子午圈，因而恒星时具有地方性，又称地方恒星时。当从格林尼治子午线上观测时，所得的恒星时称为格林尼治恒星时。由于地球自转受岁差、章动的影响，春分点的空间位置并不唯一，有真春分点和平春分点之分这导致恒星时可分为真恒星时和平恒星时，因而对格林尼治恒星时有格林尼治真恒星时(GAST)和格林尼治平恒星时(GMsT)这两者之间的关系为：GAST-GMSTAcos(32)其中，笋为黄经章动，'为黄经交角。2太阳时以真太阳周日视运动所确定的

29、时间称为真太阳时。但据天体运动的开普勒定律，太阳视运动的速度不是均匀的，以真太阳作为观察地球自转的参考点，不符合时间系统的基本要求，因而假定了一个参考点，其在天球上的视运动速度，等于真太阳周年运动的平均速度，这个假定的参考点，在天文学上被称为平太阳。以平太阳连续两次经过本地子午圈的时间间隔，定义为一个平太阳日，其1/24为一平太阳时(MeansolarTim。，MT)。与恒星时一样，平太阳时也具有地方性，常称地方平太阳时。3.世界时以地球上格林尼治子午圈所对应的平太阳时且以平子夜起算时间系统，称为世界时(universalTime，uT)。世界时与平太阳时的尺度标准完全一致，仅仅是起算点有所不

30、同。若有编表示平太阳相对格林尼治子午圈的时角，定义有世界时UTO可表示为由于地球自(3-3)转的不均匀性，使地球自转轴产生了极移现象因而UTO并不均匀，为补偿这一缺陷，国际天文联合会在世界时中引入地轴极移修正兄和地球自转变化的季节性改正参数双由此可得世界时UTI和uTZ：-(3-4)(3-5)UT2其中观测瞬时地极相对国际协议地极原点(CIO)的极移修正兄的表达式为(3-6)式中X'，厂为观测瞬间的极移分量；凡，汽分别为天文经度和纬度。地球自转速虧、O.022xsin2/-0.012xcosO.OOtixsin<O.OOTxcos4度的季节性变化改正兀有如下的经验公式：(3-7)

31、t为自本年起始日起算的年小数部分(即为计算时年积日与该年全年积日的比例)。上述修正并不能完全消除地球自转速度变化率和地球自转季节性变化的影响，故而UT:并不是严格均匀的时间系统。3.2.2原子时系统原子时以物质内部原子跃迁时所辐射和吸收的电磁波频率来定义的，其秒长定义为:位于海平面上的艳原子侧”基态两个超精细能级，在零磁场中跃迁辐射振荡9,192,631，770周所持续的时间，为1原子秒，该原子时秒作为国际制秒(si)的时间单位。原子时的起点是定在1958年1月1日0时0分0秒(UT2)，但与之又有微小误差，关系为：zir-(72；-0.0039(5)(3-8)原子时具有很高的稳定性和复现性，

32、是现时段最为理想的时间系统。许多国家都建立了各自的原子时系统，国际时间局为消除差异，对100座时钟作了对比分析，利用数据处理推算出了统一的原子时系统一一国际原子时(hiternationalAiomicTime，TAD。在目前的导航定位系统中，均采用了原子时作为其高精度的时间基准。3.2.3动力学时系统动力学时(DynamicTime，DT)是天体力学中用以描述天体运动的时间单位。当以太阳系质心建立起天休运动方程时，所采用的时间参数称为质心力学时(BaryeeniricDynamicTime,TDB);当以地球质心建立起天体运动方程时，所采用的时间参数称为地球力学时仃°C仃estri

33、alDynamicTime,TDT)。TDT所采用的基本单位为sl,与原子时一致。国际天文学联合会定义1977年1月1日TAI与TDT的严格关系为:TDT=TAI+32.184(s)(3-9)3.2.4协调世界时原子时尺度均匀稳定，但与人类日常生活紧密相关的是以地球自转为基础的世界时，在很多的科学研究中均采用的是世界时。世界时受地球速度长期性渐慢的影响，逐渐比原子时慢，为避免两都之间误差的扩大，自1972年起，国际上开始采用一种以原子时秒子为基础，在时刻上尽量接近于世界时的一种折衷的时间系统，称为协调世界时(eoordinateuniversalTime,uTe)。其引入了闰秒的概念，当协调时

34、与世界时的时刻相差超过士0.9(s)时，便于协调时中引入闰秒士l(s)，闰秒一般于12月31日或6月30日加入。协调时与TAI的关系如下：TAI=UTC+nX1(s)(3-10)其中，n为调整参数，其值由国际地球自转服务组织(lERs)发布。uTc是目前几乎所有国家发布时号的标准，相互之前的同步误差约为士0.2ms。3.2.5GPS时间系统全球定位系统(GPS)为保证导航和定位精度，建立了专门的时间系统一GPS时间系统(GPST)。其隶属于原子时系统，秒长采用国际制秒si，但不同于TAI，两都之前的关系为：TAI-GPST=19(s)(3-11)据协调时与TAI的关系可得：GPST=UTC+n

35、X1-19(s)(3-12)第四章卫星运动基本定律及其求解卫星在空间绕地球运动的时候，除了受到地球重力场引力的作用外，还受到了太阳，月亮和其它的天体引力以及太阳光压，大气的阻力和地球潮汐力的影响。卫星的实际运动轨道非常的复杂，很难用非常精确的数学模型加以描述。在各种力作用对卫星影响的过程当中，以地球的引力场的作用最大，而其它力的影响则相对的小得多。通常把作用到卫星上的力按其影响的大小分成两部分：一类是中心力；一类是摄动力，也称为非中心力。假定地球为均匀球体的地球引力，称为在心力，它决定了卫星运动的基本规律和基本特征，由此决定地球的轨道，可以视为理想的轨道。非中心力包括地球非球形对称的作用力，日

36、、月引力，大气阻力，光辐射压力以及地球的潮汐力等。摄动力的作用，使卫星偏离了既定的理想轨道。而在它影响下，卫星的运动称为卫星的受摄运动。而上述理想状态的卫星运动则称为无摄运动。卫星在地球的引力场当中所做的无摄运动，也称为开普勒运动，其规律可以由开普勒三大定律来描述。4.1开普勒第一定律开普勒第一定律：卫星运动的轨道是个椭圆，而该椭圆的一个焦点和地球的质心重合。ms近地点PMse2)远地点P'卫星绕地球质心运动的轨道方bs这一个定律表明了，在中心引力的作用下，卫星绕地球轨道运行的轨道面，是一个通过地球质心的精致平面。轨道椭圆一般称期为开普勒椭圆，其形状和大小都不变。在轨道上，卫星离地球质

37、心远的一点叫做远地点，近的一点就做近地点。轨道图形可以表示为如下图5：(4-1)a(i-S1+'_cos图4-1卫星的椭圆运行轨道在该式当中，R是卫星的地心距离；as为开普勒椭圆的长半径；es为开普勒椭圆的偏心率；fs为真近点角，它描述了任意时刻，卫星在轨道上面相对于近地点的位置，是时间的函数，其定义见上图所示。开普勒定义定律阐述了卫星运动轨道的基本形态及其与地心的关系。4.2开普勒第二定律开普勒第二定律：卫星的地心向径，即地球质心与卫星质心间的距离向量，在相同的时间内所扫过的空间面积是相等的。（如下图6所示）（图与任何其它的运动物体一样，在轨2)I道上面运动的卫星,具有两种的能量:位

38、能和动能。位能就是指仅仅受到地球重力场的影响，其大小和卫星的在轨高度有关。在近地点其位能最小，而在远地点其位能最大。卫星在任一个时刻t所具有的位能为rG为万有引力常GMm量，M为地球的质量,ms为卫星的质量）。动能则是由卫星的运动所引起向径在其同的是卫星的运扫描的度的函数。如果取卫星的运动的速度为'则其动能为2mv。根据能量守恒定律，卫星的势能与动能的总量是不变的，即1 GMm2 S2s+,=常量(4-2)因此，当卫星运行到近地点的时候，其动能最大；在远地点的时候，其动能最小，由此，开普勒第二定律所包涵的内容是：卫星在椭圆轨道上的运行速度是不断变化的，在近地点处的速度最大，而在远地点的

39、速度最小。4.3开普勒第三定律开普勒第三定律：卫星运动周期的平方与轨道椭圆长半径的立方之比为一个常数，而该常数等于地球引力常数和地球质量的乘积GM的倒数。其数学表达式为：(4-T24.2s=As_g3）在这个式子当中，Ts为卫星的运行周期如果我们假设卫星的平均角速度为N,则有:2.N=(rad/s)(4-4)于是，开普勒第三定律4-2就可以写成：Ts2(NTs)2AssG(4-5)或者表示为常用的形式：N=:G1As3(4-6)显然，当开普勒的长半径确定了过后，卫星运动的平均角速度就得到了确定,且保持不变。4.4卫星的无摄运动参数卫星的无摄运动，一般的可以由下面的6个参数（图7）来描述：As卫

40、星轨道的长半径Es卫星轨道的偏心率q升交点的赤径i卫星轨道面的倾角3S近地点角距，即升交点与近地点的夹角fs卫星的真近点角，在轨道平面上为卫星与进地点的地心角距。图4-3开普勒轨道参数当这6个参数一旦确定后，卫星在任意瞬时的相对于地球的空间位置及其速度，就被唯一的确定了！4.5真近点角的概念及其求解在描述卫星无摄运动的6个参数当中，只有fs是关于时间的函数，其他的都是一般的参数。所以，计算卫星瞬时的位置的关键，计算出参数fs,并由此确定卫星的空间位置及其和时间的关系。为此，需要引进两个参数Es和Ms去计算真近点角。Es:偏近点角，如果定义过卫星质心做平行与椭圆短半轴的直线，M'为该直线

41、与近地点到椭圆中心连线的交点，则椭圆平面上近地点P到M'的圆弧所对应的圆心角就是Es。Ms:平近点角。它是一个假设量，如果卫星在轨道运行的平速度为n,则平近点角定义为：Ms=n(t-t0)(4-7)tO为卫星过近地点的时刻，t为观察卫星的时刻。由上面的式子知道，卫星的平近点角仅仅为卫星平均速度的时间的函数，对于一个确定的卫星来说，这个参数可以认为是常数。其中Ms与Es有关系如下：Ms二Es-essinEs(4-8)为了计算卫星的瞬时速度，需要确定卫星运行的真近点角fs。由于有以下的关系成立：ascosEs二rcosfs+ases(4-9)于是将上式带入到(4-1)中就得到：cosEs-

42、eCosfs=尸(4-10)1-ecosEss或者得到以下常用的形式：Tan(f)=*£tan(号)(4-11)4.6卫星瞬时位置的求解对于任意的观测时刻，根据卫星的平均运行速度n,根据4-9,4-10,4-11，便可以唯一确定真近点角fs。这样，卫星于任一观测历元t,相对于地球瞬间空间的位置便可以随之确定。若以直角坐标的原点与地心M重合，E轴指向近地点且垂直于轨道的s平面，耳轴在轨道平面上垂直E轴构成右手关系。于是，卫星任意时刻的轨道ss坐标可以表示成为：cosfs“=sinfss(4-12)COSE-eSinEfl-0(413(1)(2)(3)sss恳轴顺时针旋转角度3S,s恳轴

43、顺时针旋转角度i,彳s恳轴顺时针旋转角度Q,s轴的指向由近地点变为升交点。s轴与Z轴相同。s实现上述三步的旋转矩阵分别为R1,吏J轴指向春分点。sR2,R3,即：COS3SSin«sR1=Sin«sCos«s(414)_0R2=CosiSini而由上面的分析，可以得到匕=s而要把这个轨道坐标系坐标表示成为天球坐标的话，由于他们的坐标原点都是地球的质心，但是坐标轴的指向是不相同的，为了使他们的坐标轴相同，应该将坐标系（E，耳,匚）依次做下面的变化旋转：1(415)SiniCosi厂CosQSinQR3=SinQCosQ(4-16)于是得到了在天球坐标系下面的卫星位置

44、坐标可以表示为：_X_飞-sCosE-esY二R3R2R1ns=R3R2R1asSinEjl-es2(4-17)Z0s利用转换关系Rz(0g)得到相应的地球坐标系的坐标了!rCos®gSing0=Sin®gCosg(4-18)_001-第五章GPS的matlab仿真5.1卫星可见性的估算当初，卫星星座的设计要求在全球范围内任何时候，任何位置都必须保证至少四颗以上的卫星导航信号。换句话说，并不是所有的卫星都能够被一个用户所看见。一颗卫星信号能否被接收与下列因素是有关系的,1地球是否影响了该GPS卫星信号的传播。下图8说明了地球对GPS信号的影响。如果卫星处于图中的阴影部分，则

45、对图中的飞机是不能够接收到该卫星的信号的！2GPS接收机是否位于该GPSE发射天线的范围内。这种情况主要针对的是航天器上面的GPS接收机。GPS卫星信号的发射张角大约为卫星到水平面的张角13.921.3度，大于可以更多的接收信号了。，这样就保证了一些飞行高度较高的攵到GPS卫星的信号。但是对于那些超出发射角的J航天器在高空飞行器就收不到3利用卫星一地心一用户之间的张角，可以估算出来可卫星,'言号能否被收到。如果这个角度小于90度的话，就可以收到，反之不能够收到。特别的，当图5-1可见卫星的测量原理这个角度刚刚为90度的时候，我们一般认为是收不到的。5.2GPS卫星运动的matlab仿真

46、程序主体见附录在进行仿真之前，有几个子程序段需要说明一下functionplot3c(x,y,z,color)switch(color)case0plot3(x,y,z,'w-');case1plot3(x,y,z,'r-');case2plot3(x,y,z,'g_，)；case3plot3(x,y,z,'c-，);case4plot3(x,y,z,，m，）;case5plot3(x,y,z,'y-，);case6plot3(x,y,z,，b-，);case7plot3(x,y,z,，k-，);otherwiseend它主要用来画出卫星

47、轨道的曲线，用不同的颜色：w代表的白色，其他的依次代表的是红色，绿色，青绿色品红色，黄色，蓝色，黑色.functionboxplot3(x,y,z,Lx,Ly,Lz,color)x0=x-Lx/2;y0=y-Ly/2;z0=z-Lz/2;x=x0x0x0x0x0+Lxx0+Lxx0+Lxx0+Lx;y=y0y0y0+Lyy0+Lyy0y0y0+Lyy0+Ly;z=z0z0+Lzz0+Lzz0z0z0+Lzz0+Lzz0;index=zeros(6,5);index(1,:)=12341;index(2,:)=56785;index(3,:)=12651;index(4,:)=43784;ind

48、ex(5,:)=26732;index(6,:)=15841;fork=1:6plot3c(x(index(k,:),y(index(k,:),z(index(k,:),color)holdonend它主要是用来表示用户的空间位置的。functiondrawearth(time)%time是参数%利用这个参数，可以绘制一个看起来是旋转的地球r=6400;time=0;jl二0:pi/10:2*pi;wl二-pi/2:pi/10:pi/2;L1=length(w1);L2=length(j1);forn=1:L1z=ones(L2,1);z=z*r*sin(w1(n);temp=r*cos(w1

49、(n);x=temp*sin(jl);y=temp*cos(jl);plot3(x,y,z);holdon;grid;end%figure(3);unit=ones(1,1);z0=ones(1,1);x0=ones(1,1);y0=ones(1,1);forn=1:L2%n=7;form=1:L1temp二wl(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=z0z1;x0=x0x1;y0=y0y

50、1;endzO(:,l)=;xO(:,l)=;yO(:,l)=;plot3(x0,y0,z0);axisequal;axisoff;holdon;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=O;cube=100;bo

51、xplot3(tempx,tempy,tempz,cube,cube,cube,2);tempx=0;tempy=0;tempz=6400;cube=100;boxplot3(tempx,tempy,tempz,cube,cube,cube,4);endfunctiondrawsatellite(moveX,moveY,moveZ,color);Length=500;Width=500;Height=500;boxplot3(moveX,moveY,moveZ,Length,Width,Height,color);boxplot3(moveX,moveY,moveZ,Length,Width*

52、5,Height/10,color);它的作用就是画出一个在上述地球坐标图当中的卫星。仿真程序一：绘制卫星的轨道平面程序见附录：地球质图5-2卫星轨道平面仿真图仿真程序二：单颗卫星不同时刻的动态仿真程序见附录：图5-3单颗卫星动态仿真图由上图10可知，卫星的天空瞬时位置是随着时间的变动而发生变化的。对于我们在地球上面的一个用户来说，一天当中的不同时刻看到的卫星是不相同的。这个仿真程序的功能实际上就是仿真了在一个轨道上面的卫星在不同的时候（这里以一个小时为一个观测时元，进行动态的在屏幕上显示器位置）通过改变时间可以显示出不同时间内的卫星的瞬时位置），而且，通过matlab当中的三维旋转图标，我们

53、可以很清晰地从不同角度看到卫星和坐标原点（这里用黑点表示的的相对位置的变化）。哪么，在其他轨道面上的卫星的运动也可以类似的模拟出来，这里不在重复！仿真程序三：卫星在某个时刻的全轨道平面的分布和可见卫星程序见附录：图5-4全轨道平面的图形其在轨坐标分别如下：SatellitePosition=1.0e+004*1.77461.75720.73650.0001-1.21610.97322.10910.00012.28830.2975-1.41320.0001-2.3882-0.88360.78690.0001-0.3681-2.52550.70120.0001-0.13232.7059-0.000

54、00.0001-1.54240.42732.17040.00011.45820.4926-2.11440.0001-2.40900.65830.73650.00011.07591.1599-2.17570.00011.3324-1.81781.43920.0001-1.32251.7688-1.3835q0.00011.2127-0.97742.10910.0001-1.87660.3838-1.90580.0001-2.4302-0.9295-0.75160.00010.44431.8187-1.84630.0001-0.67182.4562-0.78690.00011.4000-1.30731.90580.0001-1.49920.4223-2.10910.0001-0.8706-2.0297-1.38350.0001-1.7722-0.6021-1.92170.0001-0.2481-1.5812-2.17040.0001-0.54211.79151.92170.00011.9376-1.5756-0.73650.0001这个程序仿真了在某个时刻（在程序里面是在时刻timenow=