《GPS伪距单点定位ppt课件.ppt》由会员分享，可在线阅读，更多相关《GPS伪距单点定位ppt课件.ppt（34页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用 GPS伪距单点定位经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用目录CONTENTS整体综述Research Background01数据读取Main Function02坐标转换Designing Scheme04卡尔曼滤波 Chart05时间转换Performance Analysis03经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金
2、额为消费者购买商品的价款或接受服务的费用整体综述 RINEX文件读取（导航和观测文件）计算卫星位置和进行校正进行卡尔曼滤波处理后输出结果。坐标转换和时间转换3经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用卫星数据RINEX(Receiver Independent Exchange fromat，与接收机无关的交换格式)是一种GPS测量应用中普遍采用的数据格式。一般来说，GPS接收机传输到电脑上的数据是以二进制文件的形式进行存储，此方法有存储效率高，各类信息全的优点。但是，不同厂商或不同型号接收机数据格式不同，
3、无法相互识别。RINEX格式文件就是这样的背景下，为方便进行各类GPS二进制数据间的转换和传输而产生的。RINEX从开始的1.0，到2.0，最新的版本是RINEX 3.01。下面介绍的内容主要介绍RINEX第二版RINEX由来4经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用卫星数据在RINEX 格式的第2 版中定义了6 种不同类型的数据文件,它们分别是：1、观测值文件(O文件，存放GNSS观测值)；2、导航电文文件(N文件，存放GNSS卫星导航电文)；3、气象数据文件(M文件，存放在测站处所测定的气象数据)；4
4、、GLONASS导航电文文件(G文件，存放GLONASS卫星导航电文)；5、GEO 导航电文文件(H文件，存放在增强系统中搭载有类GPS信182号发生器的地球同步卫星(GEO)的导航电文)；6、卫星和接收机钟文件(C文件，存放卫星和接收机时钟信息)；比较常见的是O、N、M三种文件，O、N两种文件在进行数据处理分析时通常是必须的。文件类型5经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用卫星数据RINEX文件有严格的文件结构和格式。分为两个部分，文件头（HEADER）和数据记录。每行不得超过80个字符，前60个字符
5、为实际内容，后20个为内容标签。1.GPS1.GPS导航电文文件导航电文文件（NN文件）文件）接收机锁定卫星并解出C/A码后，就能取得广播星历，即卫星坐标计算参数，在实时GPS应用中，它是必不可少的，大部分的工程网观测数据的后处理也采用广播星历。RINEX 2格式的广播星历文件如表1-1所示。2.2.观测文件观测文件观测文件观测文件（OO文件）文件）文件）文件）接收机测得的相位和伪距观测值均记录在观测文件中，如表1-3所示。文件结构6经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用卫星数据图1-1 RINEX 格
6、式观测值文件的结构7经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用卫星数据表表1-1 RINEX 2格式的广播星历文件格式的广播星历文件RINEX版本号版本号指本文件类型指本文件类型为广播星历为广播星历第二行是生成该文件的单位、第二行是生成该文件的单位、执行人、及文件形成日期。执行人、及文件形成日期。表示结束表示结束电离层改电离层改电离层改电离层改正参数正参数正参数正参数给出了给出了GPST与与UTC之间的跳秒数之间的跳秒数每每8行表示一颗卫星的广播星历，第一行表示一颗卫星的广播星历，第一行的第行的第1个数字是
7、卫星的个数字是卫星的PRN号码号码(伪伪随机编号随机编号)。对应量的含义如表。对应量的含义如表1-2表表1-2 RINEX 2格式广播星历数据含义格式广播星历数据含义8经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用导航数据 上表中的前几行为表头，表头中的第6080个字符是相应行的说明，如第一行的“2”是RINEX版本号，“NAVIGATION DATA”是指本文件类型为广播星历；第二行是生成该文件的单位、执行人、及文件形成日期。表头以“END OF HEADER”表示结束，表头结束有时也用空行表示。表头中注解为
8、“ION ALPHA”和“ION BETA”的两行指相应的参数是电离层改正参数(采用差分模型相对定位时，该参数无用)，注解为“LEAP SECONDS”的行给出了GPST与UTC之间的跳秒数，在表头结束前可插入无限多的注解行，注解行的说明为“COMMENT”。表头结束后，每8行表示一颗卫星的广播星历，第一行的第1个数字是卫星的PRN号码(伪随机编号)。9经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用10经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的
9、价款或接受服务的费用导航数据11经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用导航数据采用的是上海跟踪站（shao）2011年7月20号7：59的数据，其中有172个广播星历。12经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用观测数据表表1-3 RINEX 2格式的观测数据文件格式的观测数据文件接收机号码、类型、版本接收机号码、类型、版本天线号码、类型天线号码、类型天线天线L1相位中心与点位的关系相位中心与点位的关系(向上、向东
10、、向北的偏离量向上、向东、向北的偏离量)观测值种类数观测值种类数及观测值类型及观测值类型注解行注解行是观测各历元的观测数据，每历元数据由一个历元时间行和多个观测数据行组成，是观测各历元的观测数据，每历元数据由一个历元时间行和多个观测数据行组成，每个历元时间行的内容为：每个历元时间行的内容为：年，月，日，时，分，秒，质量标记，卫星数，卫星号码，卫星号码，年，月，日，时，分，秒，质量标记，卫星数，卫星号码，卫星号码，.，钟差，钟差测站近似坐标测站近似坐标13经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用观测数据 与
11、广播星历表一样，观测文件也有一个表头，表头也以与广播星历表一样，观测文件也有一个表头，表头也以“END OF HEADER”或空行表示结束，表头中每行的或空行表示结束，表头中每行的6080字符为给出本行内容的字符为给出本行内容的说明，如表说明，如表1-3所示，观测文件的表头中列出了所示，观测文件的表头中列出了RINEX版本号、形成文件版本号、形成文件的单位、人员、点号、点名、观测者与单位、接收机号码类型版本、天线的单位、人员、点号、点名、观测者与单位、接收机号码类型版本、天线号码类型、台站近似坐标、天线号码类型、台站近似坐标、天线L1相位中心与点位的关系相位中心与点位的关系(向上、向东、向上、
12、向东、向北的偏离量向北的偏离量)、采样间隔、波长因子、观测值种类数及观测值类型、第、采样间隔、波长因子、观测值种类数及观测值类型、第一个观测历元时刻、最后一个观测历元时刻、测得卫星的号码，说明为一个观测历元时刻、最后一个观测历元时刻、测得卫星的号码，说明为“COMMENT”的注解行等。的注解行等。表头结束后，是观测各历元的观测数据，表头结束后，是观测各历元的观测数据，每历元数据由一个历元时间每历元数据由一个历元时间行和多个观测数据行组成行和多个观测数据行组成，每个历元时间行的内容为：，每个历元时间行的内容为：年，月，日，时，分，秒，质量标记，卫星数，卫星号码，卫星号码，年，月，日，时，分，秒，
13、质量标记，卫星数，卫星号码，卫星号码，.，钟差，钟差14经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用观测数据15经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用观测数据采用的是上海跟踪站（shao）2011年7月20号7：59的数据，有2880个历元的观测数据。16经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用站心坐标系 站心坐标系是以测站为原点
14、的坐标系，站心坐标系指用户所在位置点P为坐标原点，三个坐标轴分别指向相互垂直的东向、北向和天向，因而站心坐标系又称东北天（ENU）坐标系。站心坐标系的一个重要应用在于计算卫星在用户处的观测矢量和仰角。如图，用户位置点P在地心地固坐标系中的坐标为（x,y,z），卫星坐标为。17经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用站心坐标系等效表达在以P点为原点的站心坐标系中的向量为：有了在用户处的卫星观测矢量，就可以接着计算出卫星相对于用户的方位角和仰角：18经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加
15、赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用站心坐标系19经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用地心坐标系（大地坐标系）空间大地坐标系是采用大地经度（L）、大地纬度(B)和大地高(H)来描述空间位置的。1.纬度是空间的点与参考椭球面的法线与赤道面的夹角.2.经度是空间中的点与参考椭球的自转轴所在的面与参考椭球的起始子午面的夹角，3.大地高是空间点沿参考椭球的法线方向到参考椭球面的距离20经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额
16、为消费者购买商品的价款或接受服务的费用地心坐标系空间直角坐标系向空间大地坐标系的转换公式为：WGS84坐标系的参数：已知量：a,f 可以求解出b与e 长半轴：a=6378137 WGS84椭球扁率：f=1/298.257223563 椭球扁率f=（a-b）/a 椭球第一曲率半径（或叫偏心率）：e*e=（a*a-b*b）/a*a第二曲率半径：e2*e2=（a*a-b*b）/b*b 21经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用地心坐标系22经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到
17、的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用GPS时与民用时间转换GPSTGPST与与UTCUTC的关系：的关系：相同点相同点:都属于原子时都属于原子时.不同点不同点:协调时在年末可能通过跳秒来保持与世界时接近协调时在年末可能通过跳秒来保持与世界时接近,而为了而为了保持导航的连续性保持导航的连续性,GPST,GPST不能跳秒不能跳秒.GPSTGPST的表述形式：的表述形式：GPSGPS周周+GPS+GPS周内秒周内秒 GPSGPS周周:从从19801980年年1 1月月6 6日日0 0时起起算的星期数时起起算的星期数.GPSGPS周内秒周内秒:为周内累计的秒数为周内累计的秒数,
18、取值在取值在0 0到到604800604800之间之间.23经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用JD与MJD儒略日（儒略日（JDJD）定义）定义:从公元前从公元前47124712年年1 1月月1 1日日1212时开始计算的天数时开始计算的天数.示例示例:1982:1982年年1 1月月1 1日日0 0时的儒略日为时的儒略日为2444970.52444970.520072007年年1010月月2626日日9 9时时3030分的儒略日为分的儒略日为2454399.8962454399.896MJDMJD（
19、简化儒略日）是采用（简化儒略日）是采用19581958年年1111月月1717日平子夜日平子夜(JD=2400000.5)(JD=2400000.5)作为计时起点的一种连续计时法作为计时起点的一种连续计时法.MJD=JD-2400000.5MJD=JD-2400000.524经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用UTC到GPST的转换运行运行:时间零点定义的为：1980年1月5日夜晚与1980年1月6日凌晨之间0点。25经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔
20、偿的金额为消费者购买商品的价款或接受服务的费用卡尔曼滤波卡尔曼滤波的基本思想在海图作业中，航海长通常以前一时刻的船位为基准，根据航向、船速和海流等一系列因素推算下一个船位，但是他并不轻易认为船位就一定在推算船位上，还要选择适当的方法，通过仪器得到另一个推算船位。观测和推算这两个船位一般不重合，航海长需要通过分析和判断选择一个可靠的船位，作为船舰当前的位置。卡卡尔尔曼曼滤滤波波思思想想26经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用卡尔曼滤波例子：假设我们要研究一个房间的温度，以一分钟为时间单位。根据我们的经验
21、判断，这个房间的温度是恒定温度是恒定的，但是对我们的经验不是完全相信，可能存在上下几度的偏差，我们把该偏差看做是高斯白噪声。另外，我们在房间里放一个温度计，温度计也不准确，测量值会与实际值存在偏差，我们也把这偏差看做是高斯白噪声。现在，我们要根据我们的经验温度经验温度和温度计的测量值测量值及它们各自的噪声各自的噪声来估算出房间的实际温度。27经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用卡尔曼滤波假如我们要估算 k 时刻的实际温度值。首先你要根据根据 k-1 k-1 时刻时刻的温度值，来预测的温度值，来预测 k
22、 k 时刻的温度（时刻的温度（K K时刻的经验温度）时刻的经验温度）。因为你相信温度是恒定的，所以你会得到 k 时刻的温度预测值是跟 k-1 时刻一样的，假设是 23 23 度（度（*公式一）公式一），同时该值（预测（预测值）的高斯噪声的偏差是值）的高斯噪声的偏差是 5 5 度度（5 是这样得到的：如果 k-1 时刻估算出的最优温度值的偏差是 3，你对自己预测的不确定度是 4 度，他们平方相加再开方，就是 5（*公式二）公式二）。然后，你从温度计那里得到了 k时刻的温度值，假设是25 25 度度，同时该值的偏差是 4 4 度度。28经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔
23、偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用卡尔曼滤波现在，我们用于估算K时刻房间的实际温度有两个温度值：估计值23度和测量值25度。究竟实际温度是多少呢？是相信自己还是相信温度计？究竟相信谁多一点？究竟相信谁多一点？我们需要用他们的均方误差来判断。因为，（*公式三）公式三），所以我们可以估算出估算出K K时刻的最优时刻的最优温度值温度值为：度（*公式四）公式四）。得到了K时刻的最优温度，下一步就是对K+1时刻的温度值进行最优估算，需要得到K时刻的最优温度（24.56）的偏差，算法如下：（*公式五）公式五）就这样，卡尔曼滤波器就不断的把均方误差递归，从而估算出最优的温度
24、值，运行速度快，且只保留上一时刻的只保留上一时刻的协方差协方差。29经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用卡尔曼滤波（1）状态的一步预测方程：（2）均方误差的一步预测：（3）滤波增益方程（权重）：（4）滤波估计方程（K时刻的最优值）：（5）均方误差更新矩阵（K时刻的最优均方误差）：30经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用卡尔曼滤波31经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用卡尔曼滤波32经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用卡尔曼滤波 利用卡尔曼滤波算法对三差观测值进行粗差及周跳的修复。实验结果证明卡尔曼滤波可以对含噪信号进行有效的降噪，经降噪后的信号具有更好的分布。33经营者提供商品或者服务有欺诈行为的，应当按照消费者的要求增加赔偿其受到的损失，增加赔偿的金额为消费者购买商品的价款或接受服务的费用THANKS恳请老师点评