CN105738919B - 一种基于折半搜索算法的电离层穿刺点坐标计算方法 - Google Patents
一种基于折半搜索算法的电离层穿刺点坐标计算方法 Download PDFInfo
- Publication number
- CN105738919B CN105738919B CN201610088229.4A CN201610088229A CN105738919B CN 105738919 B CN105738919 B CN 105738919B CN 201610088229 A CN201610088229 A CN 201610088229A CN 105738919 B CN105738919 B CN 105738919B
- Authority
- CN
- China
- Prior art keywords
- mtd
- mrow
- mtable
- msubsup
- mtr
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于折半搜索算法的电离层穿刺点坐标计算方法，步骤S1是初始化折半搜索算法的参数；步骤S2是探测位置的计算、判定探测位置是否满足要求；步骤S3是更新边界条件，然后回到步骤S2；步骤S4是结束。本发明通过折半搜索算法计算电离层穿刺点的位置，克服了目前传统的三角公式计算电离层穿刺点只能得到近似坐标的弱点，有效的保证了电离层穿刺点始终在接收机与导航卫星之间的线段上，也保证了中心电离层模型是一个相对于地球椭球的固定高度椭球面，有利于高精度的电离层相关研究。并且，本发明提出的折半搜索算法计算电离层穿刺点坐标的精度优于现有技术。
Description
技术领域
本发明涉及电离层穿刺点坐标的计算，特别是涉及一种基于折半搜索算法的电离层穿刺点坐标计算方法。
背景技术
地球上的电离层是位于地球上高度约为60km到1000km被太阳和宇宙射线电离的高层大气，主要有D层、E层和F层组成，其中F层又分为F1层和F2层。由于电离层会对电磁波的传播产生一定影响，因此在利用导航卫星导航定位时，电离层误差是影响导航定位精度的主要误差因素之一。电离层误差对导航卫星信号产生的影响采用在信号传播路径上单位面积内的一个柱体所含电子总数表示，一般称为信号传播路径上的总电子含量(TEC)。利用这种方式求解信号传播路径上的总电子含量相对比较复杂，通常将整个电离层压缩为一个单层用以简化计算。这个单层被称为中心电离层，并且假设认为整个电离层的电子含量都集中在这个单层上。导航卫星信号传播路径与中心电离层的交点通常称为电离层穿刺点。
基于导航卫星对于电离层的相关研究都会用到电离层穿刺点的计算，如电离层单层建模，导航卫星和接收机的码间偏差计算等。因此获得电离层穿刺点的精度会影响电离层的相关研究。目前精确计算穿刺点坐标的方法主要是采用三角公式进行推导。这个方法是将中心电离层看作一个地球半径加上中心电离层高度的圆球面，所得的电离层穿刺点坐标是用地球椭球的大地坐标表示，所以这种方法得到的电离层穿刺点坐标只能是一种近似坐标，而且这种近似坐标并不能保证电离层穿刺点在在接收机与导航卫星之间的线段上。电离层穿刺点坐标采用的是地球椭球的大地坐标，相对应的中心电离层应该是一个相对于地球椭球的椭球面。同时随着对电离层研究的不断深入，目前计算电离层穿刺点的方法并不能满足当前研究需求，因此获取电离层穿刺点的高精度坐标将对今后电离层相关研究具有重要的意义。
发明内容
发明目的：本发明的目的是提供一种高精度的基于折半搜索算法的电离层穿刺点坐标计算方法。
技术方案：为达到此目的，本发明采用以下技术方案：
本发明所述的基于折半搜索算法的电离层穿刺点坐标计算方法，包括以下的步骤：
S1：初始化折半搜索算法的参数：初始化上边界、下边界和中心电离层的椭球高；
S2：探测位置的计算、判定探测位置是否满足要求：根据上边界和下边界计算探测位置的大地空间直角坐标，并将探测位置的大地空间直角坐标转换成大地坐标；将探测位置的椭球高与中心电离层的椭球高进行比较，判定探测位置是否满足要求：如果探测位置的椭球高与中心电离层的椭球高相等，则探测位置满足要求，进行步骤S4；否则，进行步骤S3；
S3：更新边界条件，然后回到步骤S2；
S4：结束。
进一步，所述步骤S1包括以下的步骤：
S1.1：初始化上边界
将导航卫星的大地空间直角坐标赋值给上边界的大地空间直角坐标如式(1)所示：
S1.2：初始化下边界
将接收机的大地空间直角坐标赋值给下边界的大地空间直角坐标如式(2)所示：
S1.3：初始化中心电离层的椭球高
对中心电离层的椭球高进行初始化。
进一步，所述步骤S2中探测位置的计算如式(3)所示：
其中，为上边界的大地空间直角坐标，为下边界的大地空间直角坐标，为探测位置的大地空间直角坐标。
进一步，所述步骤S2中的探测位置的大地空间直角坐标[X Y Z]T转换成大地坐标[B L H]T的方法包括以下的步骤：
S2.1：计算大地纬度B，确定迭代方程，如式(4)所示：
其中，a为地球椭球的长半轴，e为地球椭球的第一偏心率；
设定迭代初始值，如式(5)所示：
将式(5)带入式(4)中，进行多次迭代，直到最后两次计算得到的B值之差小于允许的误差才停止迭代；
S2.2：计算大地经度L，如式(6)所示：
S2.3：计算椭球高H，如式(7)所示：
其中，
进一步，所述步骤S3中边界条件更新的方法包括以下的步骤：
S3.1：如果探测位置的椭球高小于中心电离层的椭球高，则上边界的大地空间直角坐标不变，将探测位置的大地空间直角坐标赋值给下边界的大地空间直角坐标
S3.2：如果探测位置的椭球高大于中心电离层的椭球高，则下边界的大地空间直角坐标不变，将探测位置的大地空间直角坐标赋值给上边界的大地空间直角坐标
有益效果：本发明通过折半搜索算法计算电离层穿刺点的位置，克服了目前传统的三角公式计算电离层穿刺点只能得到近似坐标的弱点，有效的保证了电离层穿刺点始终在接收机与导航卫星之间的线段上，也保证了中心电离层模型是一个相对于地球椭球的固定高度椭球面，有利于高精度的电离层相关研究。并且，本发明提出的折半搜索算法计算电离层穿刺点坐标的精度优于现有技术。
附图说明
图1为本发明具体实施方式的电离层穿刺点的示意图；
图2为本发明具体实施方式的欧洲IGS观测站MEDI观测到的卫星轨迹图；
图3为本发明具体实施方式的方法流程图；
图4为本发明具体实施方式的方法A、方法B计算的电离层穿刺点到接收机的射线和到导航卫星的射线之间的夹角图。
具体实施方式
下面结合具体实施方式对本发明的技术方案作进一步的介绍。
本发明的方法采用的2012年5月29日欧洲的IGS观测站MEDI的观测数据和欧洲定轨中心CODE提供的GPS卫星的精密轨道数据(SP3文件)进行试验，具体的试验过程如下面所述：
为了方便，我们将采用三角公式推导计算电离层穿刺点的方法称为方法A，将基于折半搜索算法的电离层穿刺点坐标计算方法称为方法B。电离层穿刺点如图1所示。MEDI观测数据的采样间隔为30秒，设置卫星的高度截止角为10°，其观测到的导航卫星轨迹如图2所示。根据对MEDI观测数据的统计，统计结果显示PRN24的卫星没有观测值。试验中电离层的高度选取和CODE一致的506.7km，地球的半径为6378137m。
下表为具体实施方式的MEDI的坐标。
表1 MEDI的坐标
如图3所示，本发明方法的步骤如下：
S1：初始化折半搜索算法的参数
首先根据SP3文件计算MEDI每个观测时刻的GPS卫星的大地空间直角坐标再对边界条件初始化：
将导航卫星的大地空间直角坐标赋值给上边界的大地空间直角坐标如式(1)所示：
将接收机的大地空间直角坐标赋值给下边界的大地空间直角坐标如式(2)所示：
中心电离层的椭球高H初始化为506.7km。
S2：探测位置的计算、判定探测位置是否满足要求
根据上边界的大地空间直角坐标和下边界的大地空间直角坐标计算探测位置的大地空间直角坐标则探测位置的大地空间直角坐标为：
将探测位置的大地空间直角坐标转换成大地坐标利用探测位置的椭球高P.H和中心电离层的椭球高H相比较，判定探测位置是否满足要求：
1)若P.H＝H，则探测位置就是所求的电离层穿刺点位置，进行步骤S4；
2)若P.H≠H，则进行步骤S3。
S3：更新边界条件，然后回到步骤S2；
由于探测位置的椭球高和中心电离层的椭球高不相等P.H≠H，探测位置不满足要求，因此需要更新边界条件：
1)如果P.H<H，则上边界的大地空间直角坐标不变，将探测位置的大地空间直角坐标赋值给下边界的大地空间直角坐标
2)如果P.H>H，则下边界的大地空间直角坐标不变，将探测位置的大地空间直角坐标赋值给上边界的大地空间直角坐标
S4：结束。
为了验证两种方法所求的电离层穿刺点的正确性，只需要保证电离层穿刺点位于在接收机与导航卫星之间的线段上。因此计算接收机到穿刺点的射线和接收机到导航卫星的射线之间的夹角，若两条射线之间的夹角不等于零，说明计算的电离层穿刺点的位置不在接收机与GPS卫星之间的线段上；若两条射线之间的夹角等于零，说明计算的电离层穿刺点的位置在接收机与GPS卫星之间的线段上。这里采用余弦公式计算两条射线之间的夹角。计算得到两条射线夹角如图4所示。方法A计算的两条射线的夹角不等于零，夹角的范围在1°到7°之间，说明方法A所得到电离层穿刺点的位置不在接收机与GPS卫星之间的线段上。方法B计算的两条射线的夹角等于零，说明方法B所得到电离层穿刺点的位置在接收机与GPS卫星之间的线段上。
因此基于折半搜索算法的电离层穿刺点坐标计算方法获得的电离层穿刺点位于接收机与GPS卫星之间的线段上，具有很高的精度。由于折半搜索算法中是根据探测位置的椭球高和中心电离层的椭球高进行搜索的判定，保证了电离层穿刺点位于中心电离层上和接收机与GPS卫星之间的线段上，甚至可以认为就是电离层穿刺点的真实位置。
下表为具体实施方式的方法A和方法B计算的部分电离层穿刺点大地坐标以及两种方法之间的偏差。
表2 方法A和方法B计算的部分电离层穿刺点大地坐标以及两种方法之间的偏差
下表为具体实施方式的方法A和方法B计算的部分电离层穿刺点大地空间直角坐标以及两种方法之间的偏差。
表3 方法A和方法B计算的部分电离层穿刺点大地空间直角坐标以及两种方法之间的偏差
Claims (5)
1.一种基于折半搜索算法的电离层穿刺点坐标计算方法，其特征在于：包括以下的步骤：
S1：初始化折半搜索算法的参数：初始化上边界、下边界和中心电离层的椭球高；
S2：探测位置的计算、判定探测位置是否满足要求：根据上边界和下边界计算探测位置的大地空间直角坐标，并将探测位置的大地空间直角坐标转换成大地坐标；将探测位置的椭球高与中心电离层的椭球高进行比较，判定探测位置是否满足要求：如果探测位置的椭球高与中心电离层的椭球高相等，则探测位置满足要求，进行步骤S4；否则，进行步骤S3；
S3：更新边界条件，然后回到步骤S2；
S4：结束。
2.根据权利要求1所述的基于折半搜索算法的电离层穿刺点坐标计算方法，其特征在于：所述步骤S1包括以下的步骤：
S1.1：初始化上边界
将导航卫星的大地空间直角坐标
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>SAT</mi>
<mi>T</mi>
</msubsup>
赋值给上边界的大地空间直角坐标
<mrow>
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>U</mi>
<mi>T</mi>
</msubsup>
<mo>,</mo>
</mrow>
如式(1)所示：
<mrow>
<msubsup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>U</mi>
<mi>T</mi>
</msubsup>
<mo>:</mo>
<mo>=</mo>
<msubsup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mi>S</mi>
<mi>A</mi>
<mi>T</mi>
</mrow>
<mi>T</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
S1.2：初始化下边界
将接收机的大地空间直角坐标
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>SAT</mi>
<mi>T</mi>
</msubsup>
赋值给下边界的大地空间直角坐标
<mrow>
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>L</mi>
<mi>T</mi>
</msubsup>
<mo>,</mo>
</mrow>
如式(2)所示：
<mrow>
<msubsup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>L</mi>
<mi>T</mi>
</msubsup>
<mo>:</mo>
<mo>=</mo>
<msubsup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mi>S</mi>
<mi>A</mi>
<mi>T</mi>
</mrow>
<mi>T</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
S1.3：初始化中心电离层的椭球高
对中心电离层的椭球高进行初始化。
3.根据权利要求1所述的基于折半搜索算法的电离层穿刺点坐标计算方法，其特征在于：所述步骤S2中探测位置的计算如式(3)所示：
<mrow>
<msubsup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>P</mi>
<mi>T</mi>
</msubsup>
<mo>:</mo>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>U</mi>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>L</mi>
<mi>T</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>U</mi>
<mi>T</mi>
</msubsup>
为上边界的大地空间直角坐标，
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>L</mi>
<mi>T</mi>
</msubsup>
为下边界的大地空间直
角坐标，
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>P</mi>
<mi>T</mi>
</msubsup>
为探测位置的大地空间直角坐标。
4.根据权利要求1所述的基于折半搜索算法的电离层穿刺点坐标计算方法，其特征在于：所述步骤S2中的探测位置的大地空间直角坐标[X Y Z]T转换成大地坐标[B L H]T的方法包括以下的步骤：
S2.1：计算大地纬度B，确定迭代方程，如式(4)所示：
<mrow>
<mi>tan</mi>
<mi>B</mi>
<mo>:</mo>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msqrt>
<mrow>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
<mrow>
<mo>(</mo>
<mi>Z</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mi>a</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
<mi>tan</mi>
<mi>B</mi>
</mrow>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>tan</mi>
<mn>2</mn>
</msup>
<mi>B</mi>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
<msup>
<mi>tan</mi>
<mn>2</mn>
</msup>
<mi>B</mi>
</mrow>
</msqrt>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，a为地球椭球的长半轴，e为地球椭球的第一偏心率；
设定迭代初始值，如式(5)所示：
<mrow>
<msub>
<mi>tanB</mi>
<mn>1</mn>
</msub>
<mo>:</mo>
<mo>=</mo>
<mfrac>
<mi>Z</mi>
<msqrt>
<mrow>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
将式(5)带入式(4)中，进行多次迭代，直到最后两次计算得到的B值之差小于允许的误差才停止迭代；
S2.2：计算大地经度L，如式(6)所示：
<mrow>
<mi>L</mi>
<mo>:</mo>
<mo>=</mo>
<mi>arctan</mi>
<mfrac>
<mi>Y</mi>
<mi>X</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
S2.3：计算椭球高H，如式(7)所示：
<mrow>
<mi>H</mi>
<mo>:</mo>
<mo>=</mo>
<mfrac>
<msqrt>
<mrow>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mrow>
<mi>cos</mi>
<mi>B</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>N</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，
<mrow>
<mi>N</mi>
<mo>:</mo>
<mo>=</mo>
<mfrac>
<mi>a</mi>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>B</mi>
</mrow>
</msqrt>
</mfrac>
<mo>.</mo>
</mrow>
5.根据权利要求1所述的基于折半搜索算法的电离层穿刺点坐标计算方法，其特征在于：所述步骤S3中边界条件更新的方法包括以下的步骤：
S3.1：如果探测位置的椭球高小于中心电离层的椭球高，则上边界的大地空间直角坐
标
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>U</mi>
<mi>T</mi>
</msubsup>
不变，将探测位置的大地空间直角坐标
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>P</mi>
<mi>T</mi>
</msubsup>
赋值给下边界的大地空间
直角坐标
<mrow>
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>U</mi>
<mi>T</mi>
</msubsup>
<mo>;</mo>
</mrow>
S3.2：如果探测位置的椭球高大于中心电离层的椭球高，则下边界的大地空间直角坐
标
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>L</mi>
<mi>T</mi>
</msubsup>
不变，将探测位置的大地空间直角坐标
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>P</mi>
<mi>T</mi>
</msubsup>
赋值给上边界的大地空间
直角坐标
<mrow>
<msubsup>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>U</mi>
<mi>T</mi>
</msubsup>
<mo>.</mo>
</mrow>
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610088229.4A CN105738919B (zh) | 2016-02-17 | 2016-02-17 | 一种基于折半搜索算法的电离层穿刺点坐标计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610088229.4A CN105738919B (zh) | 2016-02-17 | 2016-02-17 | 一种基于折半搜索算法的电离层穿刺点坐标计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105738919A CN105738919A (zh) | 2016-07-06 |
| CN105738919B true CN105738919B (zh) | 2018-05-15 |
Family
ID=56245336
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610088229.4A Active CN105738919B (zh) | 2016-02-17 | 2016-02-17 | 一种基于折半搜索算法的电离层穿刺点坐标计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105738919B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108491616B (zh) * | 2018-03-19 | 2022-03-08 | 东南大学 | 一种基于椭球谐函数理论的电离层垂直总电子含量建模方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102928850A (zh) * | 2012-11-22 | 2013-02-13 | 桂林电子科技大学 | 一种广域电离层误差改正新方法 |
| CN103454695A (zh) * | 2013-08-20 | 2013-12-18 | 河海大学 | 一种gps电离层tec层析方法 |
| CN103760573A (zh) * | 2014-01-21 | 2014-04-30 | 北京北斗星通导航技术股份有限公司 | 电离层延迟获取方法和装置 |
| CN104536019A (zh) * | 2014-12-12 | 2015-04-22 | 中国电子科技集团公司第二十二研究所 | 一种利用电离层空间相关性实现的gnss电离层延迟修正方法 |
-
2016
- 2016-02-17 CN CN201610088229.4A patent/CN105738919B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102928850A (zh) * | 2012-11-22 | 2013-02-13 | 桂林电子科技大学 | 一种广域电离层误差改正新方法 |
| CN103454695A (zh) * | 2013-08-20 | 2013-12-18 | 河海大学 | 一种gps电离层tec层析方法 |
| CN103760573A (zh) * | 2014-01-21 | 2014-04-30 | 北京北斗星通导航技术股份有限公司 | 电离层延迟获取方法和装置 |
| CN104536019A (zh) * | 2014-12-12 | 2015-04-22 | 中国电子科技集团公司第二十二研究所 | 一种利用电离层空间相关性实现的gnss电离层延迟修正方法 |
Non-Patent Citations (1)
| Title |
|---|
| GPS航迹线的自适应折半查找化简算法;王静 等;《西安电子科技大学学报（自然科学版）》;20141031;第41卷(第5期);第155-160页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105738919A (zh) | 2016-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104061932B (zh) | 一种利用引力矢量和梯度张量进行导航定位的方法 | |
| CN102809376B (zh) | 一种基于等值线的辅助导航定位方法 | |
| CN103279642B (zh) | 无地面控制点的目标定位精度分析方法 | |
| CN108061889A (zh) | Ais与雷达角度系统偏差的关联方法 | |
| CN106970398A (zh) | 顾及卫星遮挡条件的卫星可见性分析及星历预报方法 | |
| CN105842679A (zh) | 一种国产卫星激光高度计在轨几何标定方法及系统 | |
| CN102565814A (zh) | 卫星导航系统的信号精度和定位服务可用性的评估方法 | |
| CN105301617A (zh) | 一种卫星导航系统中的整周模糊度有效性检验方法 | |
| CN104076381B (zh) | 实时精密单点定位方法 | |
| CN105044738A (zh) | 一种接收机自主完好性监视的预测方法及预测系统 | |
| CN109116392A (zh) | 复杂环境条件下多星座gnss融合高精度动态定位方法 | |
| CN106772483A (zh) | 一种基于cors系统的数据后处理方法及装置 | |
| CN106934113A (zh) | 适用于区域电离层垂直总电子含量建模的改进多项式模型的建模方法 | |
| CN107782317B (zh) | 一种倒vlbi深空单向导航定位方法 | |
| CN103823224A (zh) | 一种基于北斗卫星导航系统的分轨选星方法 | |
| CN105738919B (zh) | 一种基于折半搜索算法的电离层穿刺点坐标计算方法 | |
| CN105182367A (zh) | 一种提取bds卫星电离层穿刺点电子浓度的新方法 | |
| CN104933316A (zh) | 一种基于双参数混合正则化的电离层掩星反演方法 | |
| CN101581778B (zh) | 利用陀螺全站仪求取隐蔽点itrf框架坐标的方法 | |
| CN105842710A (zh) | 一种基于vrs差分原理的低轨双星时差频差精密修正方法 | |
| CN103776444B (zh) | 天空模式图对仿生偏振导航精度影响的云计算控制方法 | |
| CN103196451A (zh) | 一种基于Fisher信息矩阵的导航脉冲星选择方法 | |
| CN115793002A (zh) | 基于测向误差权重的双星联合无源定位方法 | |
| CN103837130B (zh) | 用于机载激光扫描系统的数据处理方法及装置 | |
| CN108008367A (zh) | 星载单航过InSAR系统电离层误差校正方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |