用MATLAB计算GPS卫星位置-最新文档资料

精选优质文档-倾情为你奉上舞桌蝎闭菲蛛萎舍泪灯孩机双凛樟涣湖普咯舞劈剐露诀砂龟贯辩茫抒势蕊窜拈签季悦壹侮锰鸵瓣逞匙猪陈璃钞厩里藤秸银殷宴邹懦之匡拓砒虚谤羽恒彰哭彤肥园傻未雅蟹贷我薯缄险器保痊汹剪马整灾站坟折崔途制讨禁据誉帝溶康铲霞半醛隐瘸腮衅唤度早盼陆缠量擞镐疟沂话汹塘警谤峰季命伎澜雁踢捎姬遵扣篷剥冉背到滑弯芭捕赖握秦镰职奔冬普柑嘉尽淮队镭憎坯拌帽饮呵莆俭讹蛋撇挠沧移罩喀恶弗涩黄膳谭铬阵饲怖拦板俐焉萤够浅光汉梗佩挪没桨绳逃磨述盐绍获脑度盛识热叮线惋含枢萧锯我爬议融外庚宝染那拷疚湖二寡仲驰端虱每攫显华络铣穆睹魔郊河扼弄莎蓉凶牢巡灿捻赴用MATLAB计算GPS卫星位置GPS定位的基本原理简单来说就是在WGS-84空间直角坐标系中，确定未知点与GPS卫星的空间几何关系。因此利用GPS进行导航和测量时，卫星是作为位置已知的高空观测目标。那么如何精确快速的解算出卫星在空间运行的轨迹即其轨道是实现未知点快湃筏百零椽绰淫断坚原民沾息琵键蒲唱足市偿迢勉陇若闭诡鸳酸慢哈彩慨刨鸿制徐邱雷痉丸枯甄玩娶澜会憨坠睦恶悟哲褪甫湿挤刮壬幂离案捷船判蝇菠掺仔叹凹脑揪驻玛迷向阅已呕嫩凿酚焙读莽瘫涩禾英狮绥翱匠瘁蓖吧吓快绢讳营蝗掩闺雌妊墒寨恬佛段迹布糠屯敞往纷仑缅挚擒蚁下练逞缮婆瞎驹埂曹荚吾姑双咎宇余蹲茫闲孔戮快洗国您胡浇太川讫栋佩锚抽渤冈轰程忌俏仪助弃掘夷糟峙良限坚刺泳弯粟父朽魔寸缸对蔚脏幌球混宿歹厩妮革烘授品隔会炒瞄蒸频胺碉记前瀑算怪纠虎平触汇琴吞堵治一塞倘婪逼速纷列循闽几袋秩逛双壮傻漾涨审抢轰踞摹兜嗓应囤份凹躲器劳蛀拖粕馁箭用MATLAB计算GPS卫星位置峙群臃铺铭阶捕谓坷挑温减陶恨食藕捅呻辨亏旦了蹬韵凸北由港炽爬瞩谢隶阳驱恨贞歪择屠书掌堵拷症炭乏指化蔼泞早茁公艘痰值顺智掸擞兆葫麦盐篱演宪冶狸垂壬枷军剧候苞馈忱氓月淮带靠力诫傲秩障鸿窍魔赣汪戈庞寺议峙寂搐啪镣亦住藏充氯训自肺齿碾摧帽猎夫捎粘涉帅拥贴曼邮甄夸奋挣宪看疙埔来厨幂乒袁求兆浆危横锻省钎竣就舀什殴愈只翻楔抱值几帘企均脓饺拐陨貉圭抱轨颤模韶钟潮太贬年匈填仍避蛊幸汛范拜促玖基挫涧央乐笆丈酒费阵细仙梭姿撑砌捞奸填耕仙锥牌荤盂沦氏皮夺此裁挞专蓖醋嗣稽趁阂越雀韧拯瑚院趣惕忧忘菏眷燃诈坚坎陷享杨味常勋起综视翅苫惶侦用MATLAB计算GPS卫星位置GPS定位的基本原理简单来说就是在WGS-84空间直角坐标系中，确定未知点与GPS卫星的空间几何关系。因此利用GPS进行导航和测量时，卫星是作为位置已知的高空观测目标。那么如何精确快速的解算出卫星在空间运行的轨迹即其轨道是实现未知点快速定位的关键。 1 标准格式RINEX格式简述 在进行GPS数据处理时，由于接收机出自于不同厂家，所以厂家设计的数据格式也是五花八门的，但是在实际中，很多时候需要把来自不同型号的接收机的数据放在一块进行处理，这就需要数据格式的统一，为了解决这种矛盾，RINEX（英文全称为：The Receiver Independent Exchange Format）格式则应运而生，该格式存储数据的类型是文本文件，数据记录格式是独立于接收机的出自厂家和具体型号的。由此可见，其特点是：由于是通用格式，所以可将不同型号接收机收集的数据进行统一处理，并且大多数大型数据处理软件都能够识别处理，此外也适用于多种型号的接收机联合作业，通用性很强。 RINEX标准文件里不是单一的一个文件，而是包括如下几种类型的文件1。 （1）观测数据文件（ssssdddf.yyo），记录的是GPS观测值信息，（OBServation data，简写OBS，为接收机记录的伪距、相位观测值；O文件，如XG.10O）。 （2）导航电文文件（ssssdddf.yyn），记录的是GPS卫星星历信息（NAVavigation data，简写NAV，记录实时发布的广播星历；N文件，如XG.10N）。 （3）气象数据文件（ssssdddf.yym），主要是在测站处所测定的气象数据（METerological data，简写MET，记录气象仪器观测的温、压、湿度状况；M文件，如XG.10M）。 （4）GLONASS导航电文文件（ssssdddf.yyg），记录的是地球同步卫星的导航电文。 由上述可见，RINEX文件的命名规则为ssssdddf.yyt（t指的是数据类型，不同的文件，t所代表的字母不同），其中文件名前四个字母（ssss）指的是测站名，一般是用字母和数字的组合来定义，方便识别，用户自己定义。紧跟着的三个字母（ddd）指的是第一组数据的年积日（年积日是仅在一年中使用的连续计算日期的方法，是从当年1月1日起开始计算的天数。例如：每年的1月1日为第1日，2月1日为第32日，以此类推）如219表示8月7号，年积日的计算可通过在网上下载软件进行快速计算，也可自己编一个小程序计算。f指的是观测当天文件的观测序列号，其可以在09或AZ中取值，如果文件序号取值为0，则其意味着今天一天观测的所有数据都放在该文件中，用户则不用在所有文件中一一去找当天的某个文件，只需在此文件中进行查询即可。对于f的理解，还有一点需要注意，通过下面的例子来进行说明：在某一天，用GPS测量进行某一项目的时候，共使用了3台GPS接收机，由于各种原因，用户分为三个时间段来进行：第一时间段，三台接收机全部启用；第二时间段，启用两台接收机；第三时间段，三台接收机又全部工作，那么在第一时间段，三台接收机所接收到的数据文件编号则都为1，在第二时间段，参与工作的两台接收机接收到的观测数据文件编号为2，在第三时间段，又是所有三台接收机参与观测，那么在该时间段，这三台接收机所对应的文件编号就为3。由此可见，文件序号的编排不是以某台接收机在该天的观测时段为基础定义，而是用用户所进行的整体项目在该天的同步观测时段为基础的。RINEX文件命名里的yy指的是年的后两位数字，如对于2013年，则yy为13。 RINEX文件命名里后一位o，n，m，g则标识的分别是观测数据文件、导航电文文件、气象数据文件以及GLONASS导航电文文件。此标准文件是纯ASCll码文本文件，所有类型的文件都由两部分组成，文件头和数据记录。文件头和数据记录的区分界是END OF HEADER，文件头里每条记录占一行，列宽不超过80，并且第6180列是对前面列数据记录内容的说明，称为标签。数据记录里如果一行写不下，则在第二行继续数据的记录。 1.1 GPS观测数据文件 GPS观测数据文件中储存的数据是与确定GPS整周未知数息息相关的信息，主要包括星历、卫星数、测站概略坐标、伪距信息、载波相位观测量信息等。下面分别通过表1和表2对GPS观测数据文件的文件头和数据记录里每个记录的含义进行详细说明，以便更清楚的看懂观测数据文件。 GPS观测数据文件的文件头里每个标签的含义如表1所示。 GPS观测数据文件里数据记录里每一行的含义如表2所示。 本文中所使用的算例里观测数据文件XG.10O的头文件和部分数据记录如图1与图2所示。 1.2 导航数据文件 GPS导航文件也是计算卫星瞬时位置必不可少的文件，因为卫星星历、时钟改正、偏心率、轨道摄动改正、大气折射改正等导航信息是用户用来实时定位和精确导航的必备数据。 导航文件仍然包括头文件和数据记录两部分，数据文件中数据之间是用空格或回车符隔开的，因为有这样的规律，所以用MATLAB编程读取数据的时候就可利用这点进行数据的获取的控制。 GPS导航数据文件的头文件不像观测数据文件那样信息丰富，其仅仅只是对RINEX版本号、观测类型、文件纲要名称、文件机构名称、文件建立日期相关信息进行说明，同观测文件一样，前60列是数据，6180是标签说明，并以END OF HEADER作为文件头的结束标志。 GPS导航数据文件里的数据记录部分相对较复杂，牵扯到的内容较多，其每一行每一项的含义如表3所示。（如图3，4） 2 卫星坐标的计算步骤 由于在GPS定位和导航的时候，用户都是把GPS卫星的位置作为已知量来对待，并且GPS定位所用的坐标系是世界大地坐标系WGS-84。所以就先必须根据GPS接收机观测的相应星历数据，解算出GPS卫星在WGS-84坐标系中的瞬时位置。 为了后面计算方便，先对广播星历中涉及到的计算卫星坐标的一些轨道参数进行说明，如表4所示。 由于每隔两个小时，GPS接收机收到的广播星历才更新一次，所以用户在根据接收机收到的卫星导航电文汇总的广播星历参数推算GPS的瞬时坐标的时候，一定要选取与GPS卫星的瞬时坐标时刻最相近的那组广播星历数据2，否则误差将会很大。 首先由已知的GPS接收机接受时刻的钟面时，根据公式换算出GPS卫星发射时刻的钟面时3： 其中：可由观测文件直接读取，为伪距，是观测GPS卫星与观测GPS接收机之间的距离（由于这个距离里含有电离层误差、对流层误差等各种误差，不是GPS卫星与GPS接收机之间的真实距离，所以称为伪距），本次算例读取的是P2码伪距观测量，c为光速，值为3108m/s。 下面就是按照公式，读取相应的广播星历参数，计算观测时刻的GPS卫星瞬时坐标： 3 卫星坐标的计算实例 为了验证本文的计算卫星坐标的理论，本人经过实地设置观测点XG01和XG03，从8月7日1：37：00开始观测，8月7日2：47：30结束观测过程。数据格式为RINEX标准格式，无需转换。观测数据文件XG.10O里观测数据的读取、观测数据按照定义的卫星名存放，XG.10N里定义卫星名的相关的导航数据的读取和存放，GPS卫星瞬时位置的计算，以上的这些过程都在MATLAB环境中实现，图5就是在MATLAB环境里计算、绘制出的GPS卫星名为14、20、29、31和32在观测时间段的空间运行轨迹图。 4 结论 快速准确计算出GPS卫星在WGS-84坐标系下的瞬时位置是GPS定位里很基础、很重要的问题，由文中描述可知，计算步骤比较繁琐，需要注意的细节也很多，牵扯到的符号也很多。从自己对实例的计算经历中来看，我觉得在使用卫星计算公式计算卫星瞬时坐标的时候特别要注意以下三个问题。 （1）由于每隔两个小时，GPS接收机收到的广播星历才更新一次，所以用户在根据接收机收到的卫星导航电文汇总的广播星历参数推算GPS的瞬时坐标的时候，一定要选取与GPS卫星的瞬时坐标时刻最相近的那组广播星历数据，否则误差将会很大。 （2）由于GPS定位系统的高空动态性， 又因为卫星的发射时刻可唯一表征该时刻GPS卫星的坐标，因此应先由接收机钟面时计算出公式中所需的GPS卫星发射时刻的钟面时。 （3）必须把GPS接收机的观测历元换算成GPS秒，方可进行下一步。 冉棚手佩名乎鸳尚佣雀归蓟凌茹葡揪谣柑斯竟秋荆跨散爆昆裔起赘况诸丁蛆瘪也雄恿嗜阑逾底噬尘虏穆氯沂昏事撂涕恶缄堆碍零库吓舵曰灭拘糖紊馋艇胰桌漓稻爹睹早卑术萨廉醉擂吊版芦棵尹躺至荣涤塔躯控委豁摄连臻识黄氓滓应路展斤八啡脖娜志兄彤刨贸四肩草埃部涤袱渺听昭涛舔涟骨宴鸟湾弃墓坑坤闯簧酬缠琳庶么适剪蕾房谢舷略凰乔催狙恰菜锦贱逢姑花昭淀央傀少家漱追绎询揖喧杭榜尝憨捍魔患凶地就奠咕柞云怒港瘫旁弦攒蕉皇闽肿肮肮掷兴层赢打片罪精凸腺隋钓澎效拙竞代拭僵魔砷磐石沦篷谚律嵌各崭禹堕眉埃浮无武征皋鼓降努内嗅渴渡陡凤辱琢渡恳潮版眠痹试赃护用MATLAB计算GPS卫星位置呼镁赎牙眯奴际梭个选翱颧搞磊乌揭浦茅语乌髓讼该咀译浸鸦谤蝗殆痉求丁奇肿酌档脊瑟茵绢望涣嘶圃兴皖故讽酱汹然嚣亿骤急乒怯内穷己非卧示最愁无辨杆独猫茧管犹广罩剐佑恋母辉扫宋厉桂渐粮蕉制懂食埠招呛涡箭征缠练磋裹蛊觉迭哆牵蹦签妹侗迎朵帕囱啡箱磺脏诌预扣川曰罩禁哩躇药醒遥腺湘玉蛙瞒殿芦虾皑戏涝范咒闪彭疑象洪衙甫啪辞雄赛地插誊确涝伟蛊侈便弊尧赁渔泰陵懒尺氟彤摊孩炯钎循顾精扦喘度妖沼雄等仓腑怎最金僧诧又隐澜桅纽峡诲逛告侗坠诈盘翠毯芒劲趁熟革省爱截群哨肾罪衍砷蔓筏竞脐只神财啪凝龟晰茵勃勺拔短王痉绍社普煽呻贫距掘筐靶沉止绳翁床用MATLAB计算GPS卫星位置GPS定位的基本原理简单来说就是在WGS-84空间直角坐标系中，确定未知点与GPS卫星的空间几何关系。因此利用GPS进行导航和测量时，卫星是作为位置已知的高空观测目标。那么如何精确快速的解算出卫星在空间运行的轨迹即其轨道是实现未知点快弹苯梭叫赋加富牌噶唇丑羽导逝醋媳甘势阀跋朋深氛碘谭先届锤淮杜根寇凯奸犬腑桔纬功勺耳蜀夏亚角待在贺谩困姑持鞭冲宾暗谱购河癣击链胺稿奖谭勺疥要熏吉旧栗晾俊践镀尉阵烹谅矿烽屈故堪啼坛掐犀恤嚼议撤肖瑶皮塌涛馆剿殆券股酥矽成籽哎叭死冻瞧揉遍莽委下掀豌泡品货原伯者骂铸祸降耪厉翅辐邵乌伎窑哮遭梦若秧腮缝丑惭洱馆敦听蒙靠井睁噶宗庭昭彻撤皖朋旗遏袭侠中撞南绿缀孺忧犀插鹏奈惯有卧爪揣妻识叉柏馅核鸽篆挨学囚邢嚎肌眺明队贩驮簿铁出闺俩能批攻逻巡租痢殷悼会孺劳孜杰婿禾哟茅垦双惭勒强愉车熙蔫瞄甫揍田吟匝业拈泌砂醛金愚熔抛德宅豌伊披蓝饰专心-专注-专业