我要投搞

标签云

收藏小站

爱尚经典语录、名言、句子、散文、日志、唯美图片

当前位置:欢乐棋牌下载手机版 > 非电离层 >

基于双频非差无电离层延迟组合观测量的导航卫星(GPS)精密星历

归档日期:06-10       文本归类:非电离层      文章编辑:爱尚语录

  基于双频非差无电离层延迟组合观测量的导航卫星(GPS)精密星历计算算法..

  2014年10 UrbanGeotechnical Investigation SurveyingOct.2014 No.5 文章编号:1672-8262(2014)05-11-04 中图分类号:P228 文献标识码:B 基于双频非差无电离层延迟组合观测量的 导航卫星( GPS) 精密星历计算算法与实践 倡收稿日期:2014—06—04作者简介:苏贝(1984—),男,工程师,主要从事卫星定位的技术工作。 (哈尔滨市勘察测绘研究院,黑龙江 哈尔滨150010) 摘要:翔实介绍了一种导航卫星精密定轨基本算法。 在开源代码基础上,实现了非差载波相位观测方程的解算和扩 展卡尔曼滤波器以及其他一些与定轨相关的精密改正。 应用区域跟踪站数据计算得到了导航卫星位置,与IGS 合精密星历相比较,其三维均方差略大于5cm,这个结果可以满足区域非差分用户的精密位置服务需求。关键词:精密定轨;扩展卡尔曼滤波;区域连续运行参考站 精密定轨理论的研究主要开始于20世纪50 年代 末,经过半个多世纪的发展已日趋完善,定轨精度随着 卫星跟踪技术的发展,测量精度的提高以及轨道动力 学模型的演变和精细化得到了明显的提升 在这之中,IGS扮演了重要的角色,它汇集了世界 知名的数据分析中心,包括JPL(Jet Propulsion Labora tory,美国喷气推进实验室)、CODE( Center OrbitDetermination Europe,欧洲定轨中心)、 GFZ GeoForschungs Zentrum,德国地学研究中心)、MIT( Massa chusetts Institute Technology,美国麻省理工学院)、SIO(Scripps Institution Oceanography,美国斯克里普斯海洋学研究所)、USNO(U.S.Naval Observatory,美 国海军天文台)、NOAA( National Oceanic Atmospheric Administration,美国国家海洋和大气管理局) 等,计算并发布精密星历,为地球科学研究、多学科应 用、商业开发及教育提供支撑。 在此基础上,各家数据 分析中心形成了业界知名的卫星导航定位数据分析软 GIPSY,CODE BERNESE,MIT 2基本算法从解微分方程的角度来看,精密定轨就是将一个 常微分方程初值问题转化为边值问题,由边值条件来 确定初值 2畅1观测模型导航卫星精密定轨的基本观测量是L1 和L2 载波 相位非差无电离层延迟线性组合观测量电离层延迟可 抵消掉。 基本观测方程如下: L1/L2载波波长( L1/L2载波频率(Hz);φ L1 L2为L1/L2 载波相位观测 tron为信号对流层路径延迟(m);N LC 为无电离层延迟线性组合的载波相位偏差(m); pcvr为卫星和接收机的天线相位中心变化(m); rel为相对论效应改正(m); phw为相位缠绕效应改正 LC为无电离层延迟线性组合测量噪声(m)。 式(1)里的站星几何距离ρ 通过解光时方程得到的, 该方程考虑了卫星运动和接收机钟差导致的误差,为: 为观测站的空间位置(m);U(t) 为惯性系到地固系的转换矩阵; apcr为卫星和接收机的天线位置改正(m); sdisp为测站位移(m)。 地固坐标通过转换矩阵可以转换成惯性系坐标, 考虑岁差章动及地球自转,转换矩阵为: 2014年10 UT1)+r(tUTC +(UT1-UTC)) 轴的旋转矩阵;GAST为格林尼治真恒星 时;GMST 为格林尼治平恒星时;x 为地极偏移;UT1-UTC 为地球自转角偏; 为章动经线分量和倾角;r 为世界时和恒星时之比。 卫星位置r 来表达,通常是采用GPS s,时间同步方程可简化为式(8),该方程仅考虑卫星速度和地球点质量效应。 LC相对于时刻t 卫星位置r ants-U(t-dt GME2r (11)式中,GM 时刻卫星s的速度(m/s);Φ rr 时刻到t时刻的卫星状态 转移矩阵。 应用映射函数NMF,对流层延迟V tron 表示如下: dry(El)ZHD+M wet (El)(1+G (ZTD-ZHD)(12) ZHD= 0.0022768P 1-0.00266cos2φ-0.00028H(13) =cotElsinAz(14) gN =cotElcosAz (15) 式中,ZTD 为天顶总对流层延迟(m);ZHD 为对流层大气水平梯度参数;Mdry 、Mwet 为映射函数的干、湿分量;Az、El 为平均海水面处的大气压(hPa);φ、H 为测站的纬度、大地高(m)。 卫星天线相位中线变化表示如下: satao(16) sateci为星固坐标系到惯性系的转换; satao为以星固坐标表示的天线偏移;θ 为相对于地面 站的天底角。 接收机天线相位中心变化表示如下: pcvr1(El)+C pcvr2(El) (19) localecef为测站坐标系到地固坐标系的转 ecc为测站标志到天线参考点的距离(m); apcr2为天线 载波相位中心的偏差 pcvr2为L1、L2 载波相位中心偏差(m)。 由地球潮汐引起的测站位移表示如下: solid为地球固体潮引起的位移(m); polar为极潮引起的位移(m)。 相对论效应和相位缠绕效应表示如下: cos-1 (24)式中,D、D′为卫星、接收机的有效偶极向量;x,y,z 为本地接收机单位向量(东、北、高度);x′,y′,z′为卫星 体坐标单位向量;k 为卫星至接收机的单位向量。 归纳了观测模型中精密测量改正的相关信息。精密测量改正项表1 项目 模型 卫星天线偏移和相位中心变化 Block II /IIA:(0 .279,0,1.023)mBlock IIR:(0,0,0)m 与天底角相关的相位中心变化 接收机天线偏移和相位中心变化 IGS_05 .PCV相位缠绕效应 Wu 模型 相对论效应 卫星时钟变化 测站位移 固体潮 海潮负荷 极移 参考国际地球自转服务决议 地球自转参数变化 参考国际地球自转服务产品 2畅2卫星轨道模型 一般而言,卫星轨道的高度决定其所受扰动的程 就精密定轨而言,任何大于10-10 量级的加速度都需要考虑进来。 因此,计算 GPS 卫星轨道时,模 型必须考虑到至少8 阶重力位影响及相关潮汐效应, 还有日月引力、太阳辐射压和相对论效应。 卫星轨道方程如下式所示,为卫星位置和速度的 常微分方程。 其中,状态转移矩阵采用方程数值积分 方式解算。 列出了计算时所采用的一些精密卫星轨道模型。 geop为重力位加速度改正(m/s srp为太阳辐射压加速度 rel为由相对论效应产生的加速度(m/s 精密卫星轨道模型表2项目 模型 重力位 EGM08 模型加上潮汐效应改正 潮汐效应 固体潮 海潮 参考国际地球自转服务决议三体引力 月亮和太阳视为点质量 行星星历 参考 JPL 采用的模型 太阳辐射压模型 参考 CODE 采用的模型 蚀模型 由地影和月影产生的半影或全影 相对论效应 参考国际地球自转服务决议 2畅3其他状态转移模型 其他的状态转移模型如下:卫星钟差采用一阶高 斯马尔科夫;接收机钟为白噪声;对流层参数为随机游 走;地球自转参数为随机游走。 在一个弧段内,载波相 位偏差认为是一个固定值。 如果在该弧段探测到了周 跳,则将其重新初始化。 2畅4参数估计 参数估计采用扩展卡尔曼滤波的方法。 下列方程 组分别表示观测更新和时序更新。 除了前向滤波外, 计算时还可采用后向滤波及其更平滑的方法。 Kk =Pk(-)H k(HkPk(-)H xk(-)+Kk{zk xk(-)]}Pk(+)=(I-Kk Hk )Pk(-) (27) (28)对于GPS 卫星定轨而言,卡尔曼滤波器所估计的 状态向量,包括卫星位置和速度、卫星时钟、接收机时 钟、对流层参数。 此外,如必要,卫星的太阳辐射压和 地球自转参数也能纳入进来。 在估计中,卫星和接收 机的时钟当作相对于参考时刻而变化的时钟来处理, 并固定参考时刻为0。 输入观测向量zk 后,tk 时刻的状态向量 xk通过卡 尔曼滤波估计过程逐步得到。 3畅1HRBCORS哈尔滨双星导航服务系统(GNSS)由哈尔滨市勘 察测绘研究院于2007 年建成,一期包含5 个连续运行 参考站点,分别是阿城、哈西、对青、方台、双城。 每站 配备一台双频72 通道接收机,型号为 NetR5,可同时 接收GPS 和GLONASS 卫星信号,是本次计算的数据 来源,如图1 所示。 图1哈尔滨市连续运行参考站分布图 3畅2轨道精度 在上面叙述的情况下,以开源的GAMIT 软件代码 为基础,添加了非差相位观测方程解算模块,修改了部 分参数估计算法,得到了 GPS 卫星的定轨结果,然后 与IGS 综合精密星历比较,得出结论。 后者由 IGS 织应用全球范围内的GPS站网,对所有分析中心给出 的结果加权平均得到,其中每家分析中心应用不同的 精密分析软件包进行数据分析,由此得到的轨道精度 一般能够达到 cm以内。 列出了卫星轨道的估计精度,以IGS综合精 密星历为标准,给出了所有卫星位置均方跟误差的平 均值。 比较了所有卫星的轨道估计误差。卫星轨道精度表3 三维 径向 2014年10 本文叙述了一种GPS卫星精密定轨算法,结合实 测数据计算了卫星星历,与IGS 综合精密星历相比较, 轨道外符合精度略大于 cm,即卫星位置的三维均方差约为 cm。作者下一步可继续进行精密钟差估计 的工作,为区域非差分用户提供精密定位和授时服务。 参考文献 [1]刘林.航天器轨道理论[M].北京:国防工业出版社,2000. [2]施闯,赵齐乐,李敏等.北斗卫星导航系统的精密定轨与 定位研究[J].中国科学:地球科学,2012,42:845 ~861. 精密定轨进展及其数学模型[J].大地测量与地球动力学,2007,27:112 ~118. [4]周建华,杨龙,徐波等.一种导航卫星中长期轨道预报方 法[J].测绘学报,2011,40(S):39 ~45. [5]李敏.多模GNSS 融合精密定轨理论及其应用研究[D]. 武汉:武汉大学,2011. [6]韩保民.基于星载GPS 的低轨卫星几何法定轨理论研究 [D].武汉:中科院测量与地球物理研究所,2003. [7]叶世榕.GPS 非差相位精密单点定位理论与实现[D]. 武汉:武汉大学,2002. [8]Oliver M.etc..Satellite Orbits: Models,Methods Applications[M].Springer-Verlag,Heidelberg,2000. Algorithms NavigationSatellites(GPS) Precise Ephemeris Based IndifferenceIonosphere-free Linear Combination L1/L2 Measurements Su Bei (Harbin Institute GeotechnicalInvestigation Surveying,Harbin150010,China) Abstract:It’s presented basickind navigationsatellite orbit determination. Based opensource codes,the software non-differentialcarrier phase,to estimate parameters using extended kalman filter otherprecise corrections associated orbitdetermination. And satelliteposition regionalcontinuously operating reference stations.Compared IGSfinal orbit,the estimated satellite position 3D RMS error slightlylarger than 5cm.This result can meet regionalnon-differential users preciselocation services. Key words:precise orbit determination;extended kalman filter;regional continuously operating reference stations (上接第10 PreciseRelative Positioning Single-frequencyGNSS Receivers Yan Wei,Wang Wei,Zhang Qi,Huang Yong (Tianjin Institute Mapping,Tianjin300381,China) Abstract:The rapid precisepositioning low-costsingle-frequency receivers over widearea has become one researchhot spot presentsatellite geodesy paper,amethod real-timepre cise relative positioning using single-frequency practicalsingle-frequency GPS observables from different length baselines bothIGS CORSnetwork.The result support positioningaccuracy ourmethod shortbaselines, while about5millimeter cmrespectively verticaldirections. Key words:single-frequency;GNSS;precise relative positioning;real-time

  基于双频非差无电离层延迟组合观测量的导航卫星(GPS)精密星历计算算法与实践

本文链接:http://halomarie.com/feidianliceng/3.html

上一篇:没有了

下一篇:电离层非相干散射软件雷达的研究