CN108491616A - 一种基于椭球谐函数理论的电离层垂直总电子含量建模方法 - Google Patents
一种基于椭球谐函数理论的电离层垂直总电子含量建模方法 Download PDFInfo
- Publication number
- CN108491616A CN108491616A CN201810223959.XA CN201810223959A CN108491616A CN 108491616 A CN108491616 A CN 108491616A CN 201810223959 A CN201810223959 A CN 201810223959A CN 108491616 A CN108491616 A CN 108491616A
- Authority
- CN
- China
- Prior art keywords
- ellipsoid
- model
- electron content
- total electron
- vertical total
- 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.)
- Granted
Links
- 239000005433 ionosphere Substances 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000000295 complement effect Effects 0.000 claims description 12
- 238000010606 normalization Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000013316 zoning Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 7
- 238000012937 correction Methods 0.000 description 6
- 239000010410 layer Substances 0.000 description 6
- 238000011160 research Methods 0.000 description 3
- 238000009499 grossing Methods 0.000 description 2
- NUHSROFQTUXZQQ-UHFFFAOYSA-N isopentenyl diphosphate Chemical compound CC(=C)CCO[P@](O)(=O)OP(O)(O)=O NUHSROFQTUXZQQ-UHFFFAOYSA-N 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000003574 free electron Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011541 reaction mixture Substances 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- 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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于椭球谐函数理论的电离层垂直总电子含量建模方法，包括以下步骤：S1：获取区域观测点的穿刺点地理经度、穿刺点地理纬度和穿刺点垂直方向的电子含量；S2：对地心纬度进行归化计算；S3：选定椭球谐函数的阶数，采用椭球谐函数建立区域电离层垂直总电子含量模型；S4：通过最小二乘法计算区域电离层垂直总电子含量模型的参数，对区域电离层垂直总电子含量模型进行修正。本发明与现有技术相比，更加准确的描述了电离层延迟信息随时间的非线性振动特性，较好的模拟了区域电离层的时空变化特征，有效提高了模型精度。
Description
技术领域
本发明涉及一种基于椭球谐函数理论的电离层垂直总电子含量建模方法。
背景技术
在GNSS电离层建模的相关研究和应用中，目前较为常用的电离层修正模型总体上分为两大类：其中一类是经验模型，包括GPS中使用的Klobuchar模型、IRI模型、Bent模型等，由于经验模型的建立是以长期的收集资料为基础的，因此部分模型结构较为复杂且适用性较差，通常情况下在电离层模型建模以及GNSS定位精度要求较高时不能满足要求；另一类是基于高精度GNSS双频观测数据构建的拟合电离层模型，一般常用的电离层模型包括下面三种：多项式模型(POLY)、球谐函数模型(SHF)、三角级数(TSF)模型。拟合模型相比于经验模型精度有一定的提高，但是，拟合模型的绝对精度仍然不是很高，球谐函数模型(SHF)是目前应用和研究较为广泛的一种拟合模型，也能够一定程度的对电离层进行描述，但众所周知，地球表面更接近于一个旋转椭球面，而不是一个球面，使用球谐函数建立电离层延迟模型会产生一定的偏差，地球扁率会对于高精度定位会产生不利影响，球谐函数模型因忽略扁率级误差己不能满足当前相对高精度定位和导航应用的要求，因此需要建立一个更加符合地球形状且可以广泛应用的模型，从而提供更高精度的服务，满足更高精度测量的需求。
发明内容
发明目的：本发明的目的是提供一种能够提高精度的基于椭球谐函数理论的电离层垂直总电子含量建模方法。
技术方案：为达到此目的，本发明采用以下技术方案：
本发明所述的基于椭球谐函数理论的电离层垂直总电子含量建模方法，包括以下步骤：
S1：获取区域观测点的穿刺点地理经度、穿刺点地理纬度和穿刺点垂直方向的电子含量；
S2：对地心纬度进行归化计算；
S3：选定椭球谐函数的阶数，采用椭球谐函数建立区域电离层垂直总电子含量模型；
S4：通过最小二乘法计算区域电离层垂直总电子含量模型的参数，对区域电离层垂直总电子含量模型进行修正。
进一步，所述步骤S2中，通过以下方法对地心纬度进行归化计算：
首先根据式(1)计算出归化后地心纬度的余角θ：
其中，E为线性偏心率，u为椭球坐标系中椭球的短轴，α为地心纬度的余角；
然后根据式(2)计算出归化后的地心纬度θ'：
进一步，所述步骤S3中，根据式(3)建立区域电离层垂直总电子含量模型：
其中，VTEC为区域电离层垂直总电子含量，θ为归化后地心纬度的余角，λ为地心经度；anm、bnm为未知的椭球谐函数系数，即待求的全球或区域性VTEC参数，n为椭球谐函数的度数，m为椭球谐函数的阶数；Pmn(cosθ)为经典拉格朗日函数。
进一步，所述anm、bnm通过式(4)计算得到：
其中，0≤m≤2，0≤n≤2，A如式(5)所示，l是VTEC矩阵；
其中，为n度m阶的归化拉格朗日函数，MC(n,m)为归化函数，即δ0m为Kronecker型的δ函数，为经典拉格朗日函数；N'为所选取的建立模型的数据数量，1≤i1≤N′，为第i1个归化后地心纬度的余角，为第i1个地心经度。
有益效果：本发明公开了一种基于椭球谐函数理论的电离层垂直总电子含量建模方法，基于椭球谐函数理论对地心纬度进行转化，与现有技术相比，更加准确的描述了电离层延迟信息随时间的非线性振动特性，较好的模拟了区域电离层的时空变化特征，有效提高了模型精度。
附图说明
图1为本发明具体实施方式中球谐函数模型在0点～2点的残差绝对值分布散点图，其中X轴为点位序号；Y轴为模型拟合结果残差绝对值，单位为TECU；
图2为本发明具体实施方式中椭球谐VTEC模型在0点～2点的残差绝对值分布散点图，其中X轴为点位序号；Y轴为模型拟合结果残差绝对值，单位为TECU；
图3为本发明具体实施方式中球谐函数模型在6点～8点的残差绝对值分布散点图，其中X轴为点位序号；Y轴为模型拟合结果残差绝对值，单位为TECU；
图4为本发明具体实施方式中椭球谐VTEC模型在6点～8点的残差绝对值分布散点图，其中X轴为点位序号；Y轴为模型拟合结果残差绝对值，单位为TECU；
图5为本发明具体实施方式中球谐函数模型在12点～14点的残差绝对值分布散点图，其中X轴为点位序号；Y轴为模型拟合结果残差绝对值，单位为TECU；
图6为本发明具体实施方式中椭球谐VTEC模型在12点～14点的残差绝对值分布散点图，其中X轴为点位序号；Y轴为模型拟合结果残差绝对值，单位为TECU；
图7为本发明具体实施方式中球谐函数模型在18点～20点的残差绝对值分布散点图，其中X轴为点位序号；Y轴为模型拟合结果残差绝对值，单位为TECU；
图8为本发明具体实施方式中椭球谐VTEC模型在18点～20点的残差绝对值分布散点图，其中X轴为点位序号；Y轴为模型拟合结果残差绝对值，单位为TECU。
具体实施方式
下面结合具体实施方式对本发明的技术方案作进一步的介绍。
本具体实施方式为了更好的描述区域电离层的时空变化特征，采用的数据来源于利用rtklib所下载的igs站发布欧洲地区的数据。包括func、zimm、sass、mobj等在内的40余个参考站，数据时间为2014年11月20日，区域为0°～30°E，30°～60°N。为了更好地描述修正效果，一天内共分为0点～2点、6点～8点、12点～14点、18点～20点共4个时间段进行建模工作。
本具体实施方式公开了一种基于椭球谐函数理论的电离层垂直总电子含量建模方法，包括以下步骤：
S1：获取区域观测点的穿刺点地理经度、穿刺点地理纬度和穿刺点垂直方向的电子含量；
S2：对地心纬度进行归化计算；
S3：选定椭球谐函数的阶数，采用椭球谐函数建立区域电离层垂直总电子含量模型；
S4：通过最小二乘法计算区域电离层垂直总电子含量模型的参数，对区域电离层垂直总电子含量模型进行修正。
步骤S1中，区域观测点的穿刺点地理经度和穿刺点地理纬度分别通过式(1)、式(3)得到：
式(1)中，λ为区域观测点的穿刺点地理经度，λ0为测站接收机的地理经度，Ψpp为地心张角，如式(2)所示，A为卫星方位角，为区域观测点的穿刺点地理纬度；
式(2)中，E为卫星高度角，R为地球半径，H为电离层单层模型的等效高度；
式(3)中，为测站接收机的地理纬度。
步骤S1中，区域观测点的穿刺点垂直方向的电子含量由下面方法所得：
假设地球平均半径为R，电离层薄层高度为H，在接收机处卫星的天顶距为z，在穿刺点处的天顶距为z′(天顶距与高度角互为余数)，则有如下关系式：
时钟的基本频率f0为10.23MHz，GPS上L波段的两种载波L1和L2的频率分别为：
fL1＝154×f0＝1575.42MHz (5)
fL2＝120×f0＝1227.60MHz (6)
采用伪距观测量时，根据每个历元两个频率上的伪距之差解算出天顶方向自由电子含量VTEC值，其观测方程如下：
式(7)中，i,j＝1,2，i≠j；P为伪距观测值，Pj为第j个信号的伪距观测值，Pi为第i个信号的伪距观测值，Pj-Pi为电离层延迟观测量，BS为伪距观测量的卫星硬件延迟偏差，为第j个信号的伪距观测量的卫星硬件延迟偏差，为第i个信号的伪距观测量的卫星硬件延迟偏差，BR为伪距观测量的接收机硬件延迟偏差，为第j个信号的伪距观测量的接收机硬件延迟偏差，为第i个信号的伪距观测量的接收机硬件延迟偏差，fi为载波L1的频率，fj为载波L2的频率，由此可解出电离层天顶方向总电子含量VTEC的表达式如下：
式(8)中，ΔPij是采用伪距观测时的伪距之差，为伪距观测值的卫星相对硬件延迟偏差，为伪距观测值的接收机相对硬件延迟偏差；
采用相位观测量时，其观测方程如下：
式(9)中，λ、N分别为波长和整周模糊度，λj为第j个信号的波长，λi为第i个信号的波长，Nj为第j个信号的整周模糊度，Ni为第i个信号的整周模糊度，bS为相位观测值的卫星硬件延迟偏差，为第j个信号的相位观测值的卫星硬件延迟偏差，为第i个信号的相位观测值的卫星硬件延迟偏差，bR为相位观测值的接收机硬件延迟偏差，为第j个信号的相位观测值的接收机硬件延迟偏差，为第i个信号的相位观测值的接收机硬件延迟偏差；
解方程可得：
式(10)中，LΔij为相位观测量组合，为载波相位观测值的卫星相对硬件偏差，为载波相位观测值的接收机硬件延迟偏差，Ambij为模糊度组合常数。
步骤S2中，通过以下方法对地心纬度进行归化计算：
首先根据式(11)计算出归化后地心纬度的余角θ：
其中，E为线性偏心率，u为椭球坐标系中椭球的短轴，α为地心纬度的余角；
然后根据式(12)计算出归化后的地心纬度θ'：
步骤S3中，根据式(13)建立区域电离层垂直总电子含量模型：
其中，VTEC为区域电离层垂直总电子含量，(u,θ',λ)为椭球坐标系里的坐标，u为椭球的短轴，θ为归化后地心纬度的余角，λ为地心经度；anm、bnm为未知的椭球谐函数系数，即待求的全球或区域性VTEC参数，n为椭球谐函数的度数，m为椭球谐函数的阶数；Pmn(cosθ)为经典拉格朗日函数。
anm、bnm通过式(14)计算得到：
其中，0≤m≤2，0≤n≤2，A如式(15)所示，l是VTEC矩阵；
其中，为n度m阶的归化拉格朗日函数，MC(n,m)为归化函数，即δ0m为Kronecker型的δ函数，为经典拉格朗日函数；N'为所选取的建立模型的数据数量，1≤i1≤N′，为第i1个归化后地心纬度的余角，为第i1个地心经度。
采用的数据来源于利用rtklib所下载的igs站发布的欧洲地区数据。数据时间为2014年11月20日(年积日为324)，区域为0°～30°E，30°～60°N。根据测站位置信息，利用相位平滑伪距方法提取相关电离层数据，即包括穿刺点(IPP)的经纬度、观测时间、VTEC值等信息，具体见情况见表1(a)、1(b)、1(c)、1(d)。
表1(a)部分0点～2点电离层建模数据
表1(b)部分6～8点电离层建模数据
| 测站 | 时 | 分 | 秒 | PRN | 地理纬度/° | 地理经度/° | VTEC/TECU |
| ebre | 6 | 0 | 0 | G29 | 43.312 | 6.519 | 7.261 |
| yebe | 6 | 0 | 0 | G29 | 43.35 | 3.462 | 6.951 |
| hers | 6 | 0 | 0 | G21 | 48.86 | 1.83 | 8.401 |
| hers | 6 | 0 | 0 | G29 | 52.237 | 6.949 | 5.887 |
| wsrt | 7 | 0 | 0 | G21 | 52.529 | 8.485 | 11.874 |
| wsrt | 7 | 0 | 0 | G29 | 52.410 | 20.385 | 14.987 |
| opmt | 7 | 0 | 0 | G21 | 48.872 | 4.467 | 11.486 |
| brux | 7 | 0 | 0 | G29 | 50.744 | 18.473 | 16.624 |
| brux | 8 | 0 | 0 | G21 | 51.992 | 8.855 | 20.975 |
| wsrt | 8 | 0 | 0 | G18 | 51.305 | 11.935 | 22.651 |
| wsrt | 8 | 0 | 0 | G21 | 53.813 | 11.028 | 19.761 |
| wsrt | 8 | 0 | 0 | G22 | 46.881 | 9.789 | 26.204 |
表1(c)部分12～14点电离层建模数据
| 测站 | 时 | 分 | 秒 | PRN | 地理纬度/° | 地理经度/° | VTEC/TECU |
| sass | 12 | 0 | 0 | G14 | 55.476 | 20.259 | 30.275 |
| sass | 12 | 0 | 0 | G32 | 53.367 | 12.34 | 32.794 |
| joze | 12 | 0 | 0 | G04 | 51.804 | 21.067 | 26.405 |
| joze | 12 | 0 | 0 | G14 | 52.915 | 26.453 | 25.709 |
| joze | 13 | 0 | 0 | G32 | 52.278 | 21.678 | 22.308 |
| joze | 13 | 0 | 0 | G11 | 49.808 | 20.806 | 25.084 |
| joze | 13 | 0 | 0 | G20 | 51.590 | 15.277 | 26.321 |
| joze | 13 | 0 | 0 | G01 | 51.833 | 20.126 | 25.013 |
| joze | 14 | 0 | 0 | G04 | 46.578 | 25.491 | 23.256 |
| svtl | 14 | 0 | 0 | G01 | 57.391 | 29.138 | 21.430 |
| svtl | 14 | 0 | 0 | G23 | 55.278 | 22.744 | 22.232 |
| crao | 14 | 0 | 0 | G23 | 41.739 | 29.019 | 31.821 |
表1(d)部分18～20点电离层建模数据
| 测站 | 时 | 分 | 秒 | PRN | 地理纬度/° | 地理经度/° | VTEC/TECU |
| joze | 18 | 0 | 0 | G10 | 51.387 | 19.471 | 6.955 |
| joze | 18 | 0 | 0 | G07 | 51.006 | 21.377 | 4.738 |
| joze | 18 | 0 | 0 | G05 | 55.643 | 11.692 | 6.403 |
| joze | 18 | 0 | 0 | G30 | 48.458 | 18.396 | 6.113 |
| mobj | 19 | 0 | 0 | G05 | 56.638 | 29.443 | 6.454 |
| svtl | 19 | 0 | 0 | G10 | 56.399 | 27.227 | 4.857 |
| svtl | 19 | 0 | 0 | G30 | 58.669 | 27.368 | 6.698 |
| crao | 19 | 0 | 0 | G05 | 46.986 | 27.965 | 9.431 |
| crao | 20 | 0 | 0 | G26 | 47.715 | 26.563 | 10.273 |
| crao | 20 | 0 | 0 | G05 | 44.305 | 28.681 | 13.439 |
| joze | 20 | 0 | 0 | G28 | 48.277 | 22.953 | 8.423 |
| joze | 20 | 0 | 0 | G07 | 53.526 | 27.205 | 2.710 |
本发明的研究采用的数据来源于利用RTKLIB所下载的IGS站发布欧洲地区的数据。包括func、zimm、sass、mobj等在内的40余个参考站，数据时间为2014年11月20日(年积日为324)，区域为0°～30°E，30°～60°N。由于时间所限，本文未能将1天中全部时段的数据都进行研究，因此本文抽取了等时间间距(间距为6小时)的部分数据进行研究，一天内共分为0点～2点、6点～8点、12点～14点、18点～20点共4个时间段进行建模工作。
根据建模区域大小，本发明决定建立的模型为2阶2次模型。
本研究根据测站位置信息，利用利用相位平滑伪距方法提取相关电离层数据，即包括穿刺点(IPP)的经纬度、观测时间、VTEC值等信息。在这段时间内共解算出1262条数据，具体为0点～2点共264条、6点～8点共231条、12点～14点共377条、18点～20点共390条。
为了方便对比两种VTEC模型的拟合效果，在建模结束并进行计算后，我们将采用中误差(root mean square error,RMSE)以及残差绝对值平均值(the mean of absoluteerror,Eabs)作为评价指标，具体公式如下：
上述式子中，N表示进行误差分析的数据个数；VTECpred表示VTEC的模型得到的结果；VTECtrue表示根据IGS站下载的数据利用相位平滑伪距方法提取的VTEC真值。
表1、2、3、4给出了4个时段内椭球谐VTEC模型同球谐函数模型的精度指标对比，表5两种VTEC模型前面数据中误差的对比。图1-图8给出了0～2点、6～8点、12～14点、18～20点四个时段数据残差绝对值散点分布图。表1、2、3、4如下所示：
表1 0点～2点两种模型修正效果对比
表2 6点～8点两种模型修正效果对比
表3 12点～14点两种模型修正效果对比
表4 18点～20点两种模型修正效果对比
通过表1、2、3、4我们可以得到以下分析结果：
①在0点～2点时段，除了在1点时刻球谐函数模型的残差绝对值平均值小于椭球谐VTEC模型，且只有0.006TECU，而在其他时段的所有误差指标中，椭球谐VTEC模型的误差都是要小于球谐函数模型的。残差绝对值平均值和中误差改善量都在0.05TECU～0.1TECU以内；
②在6点～8点时段，在6点和7点时段，相比于球谐函数模型，椭球谐函数模型的结果在残差绝对值平均值和中误差两项指标中都略有减小，残差绝对值平均值和中误差改善量都在0.1TECU左右。其中6点时段残差绝对值平均值减和中误差小了约0.14TECU和0.09TECU；7点时段残差绝对值平均值和中误差减小了约0.09TECU和0.03TECU。8点时段椭球谐VTEC模型比球谐函数模型中误差减小了0.5TECU左右，残差绝对值平均值减小了0.96TECU，与前两个时段相比，在8点时段，椭球谐函数模型的提升效果相对明显一些。
③在12点时段椭球谐函数VTEC模型比球谐函数VTEC模型在两项误差指标上都减小了超过0.5TECU；在13点时段则是0.4TECU左右；在14点时段两种模型在上述指标中差距不大，残差绝对值平均值球谐函数模型比椭球谐函数VTEC模型小0.012TECU，中误差则是大了0.004TECU。在12点和13点时段椭球谐函数VTEC模型相比于球谐函数模型拟合精度要高一些，而在14点时段二者的精度相当；
④在18点时段和19点时段椭球谐VTEC模型的改善效果相当有限，残差绝对值平均值和中误差的改善量都在0.05TECU左右，在20点时段，椭球谐函数VTEC模型的拟合结果残差绝对值平均值比球谐函数减小了0.312TECU，中误差减小了0.3TECU左右，改善效果相对明显。
两种VTEC模型前面数据中误差的对比见表5：
表5两种VTEC模型拟合数据中误差对比表
总体说来，利用椭球谐函数思想建立的电离层VTEC拟合模型在绝大部分的时段各项指标均优于球谐函数模型拟合的效果，只有在极个别时段(1点时段)的残差绝对值平均值一项椭球谐函数VTEC模型稍差于球谐函数VTEC模型，而在中误差这一项上基本上每个时段都要优于球谐函数VTEC模型(14点时段可以看作二者相当)，最终相比于球谐函数模型，椭球谐VTEC模型在拟合精度上提升约9.3％。
Claims (4)
1.一种基于椭球谐函数理论的电离层垂直总电子含量建模方法，其特征在于：包括以下步骤：
S1：获取区域观测点的穿刺点地理经度、穿刺点地理纬度和穿刺点垂直方向的电子含量；
S2：对地心纬度进行归化计算；
S3：选定椭球谐函数的阶数，采用椭球谐函数建立区域电离层垂直总电子含量模型；
S4：通过最小二乘法计算区域电离层垂直总电子含量模型的参数，对区域电离层垂直总电子含量模型进行修正。
2.根据权利要求1所述的基于椭球谐函数理论的电离层垂直总电子含量建模方法，其特征在于：所述步骤S2中，通过以下方法对地心纬度进行归化计算：
首先根据式(1)计算出归化后地心纬度的余角θ：
其中，E为线性偏心率，u为椭球坐标系中椭球的短轴，α为地心纬度的余角；
然后根据式(2)计算出归化后的地心纬度θ'：
3.根据权利要求1所述的基于椭球谐函数理论的电离层垂直总电子含量建模方法，其特征在于：所述步骤S3中，根据式(3)建立区域电离层垂直总电子含量模型：
其中，VTEC为区域电离层垂直总电子含量，θ为归化后地心纬度的余角，λ为地心经度；anm、bnm为未知的椭球谐函数系数，即待求的全球或区域性VTEC参数，n为椭球谐函数的度数，m为椭球谐函数的阶数；Pmn(cosθ)为经典拉格朗日函数。
4.根据权利要求3所述的基于椭球谐函数理论的电离层垂直总电子含量建模方法，其特征在于：所述anm、bnm通过式(4)计算得到：
其中，0≤m≤2，0≤n≤2，A如式(5)所示，l是VTEC矩阵；
其中，为n度m阶的归化拉格朗日函数，MC(n,m)为归化函数，即δ0m为Kronecker型的δ函数，为经典拉格朗日函数；N'为所选取的建立模型的数据数量，1≤i1≤N′，为第i1个归化后地心纬度的余角，为第i1个地心经度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810223959.XA CN108491616B (zh) | 2018-03-19 | 2018-03-19 | 一种基于椭球谐函数理论的电离层垂直总电子含量建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810223959.XA CN108491616B (zh) | 2018-03-19 | 2018-03-19 | 一种基于椭球谐函数理论的电离层垂直总电子含量建模方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108491616A true CN108491616A (zh) | 2018-09-04 |
| CN108491616B CN108491616B (zh) | 2022-03-08 |
Family
ID=63339932
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810223959.XA Active CN108491616B (zh) | 2018-03-19 | 2018-03-19 | 一种基于椭球谐函数理论的电离层垂直总电子含量建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108491616B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110851976A (zh) * | 2019-11-07 | 2020-02-28 | 东南大学 | 一种聚焦离子束辅助沉积结构轮廓确定方法 |
| CN111273335A (zh) * | 2019-12-20 | 2020-06-12 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | 一种基于垂测数据约束的电离层层析成像方法 |
| CN111796309A (zh) * | 2020-06-24 | 2020-10-20 | 中国科学院精密测量科学与技术创新研究院 | 导航卫星单频数据同步确定大气水汽和总电子含量的方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130326453A1 (en) * | 2012-06-01 | 2013-12-05 | Wistron Corp. | Circuit layout method for printed circuit board, eletronic device and computer readable recording media |
| US20140282286A1 (en) * | 2013-03-15 | 2014-09-18 | Globalfoundries Singapore Pte. Ltd. | Etch failure prediction based on wafer resist top loss |
| CN105301601A (zh) * | 2015-10-09 | 2016-02-03 | 中国科学院光电研究院 | 一种适用于全球区域的gnss电离层延迟三维建模方法 |
| CN105738919A (zh) * | 2016-02-17 | 2016-07-06 | 东南大学 | 一种基于折半搜索算法的电离层穿刺点坐标计算方法 |
| CN106405589A (zh) * | 2016-06-24 | 2017-02-15 | 西安科技大学 | 一种全球电离层格网模型的确定方法及装置 |
| CN106405576A (zh) * | 2016-08-25 | 2017-02-15 | 中国科学院国家授时中心 | 一种基于gnss载波相位观测值的电离层监测与建模方法 |
| CN107622304A (zh) * | 2017-08-21 | 2018-01-23 | 东南大学 | 一种基于bp神经网络技术的电离层球谐函数改进方法 |
-
2018
- 2018-03-19 CN CN201810223959.XA patent/CN108491616B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130326453A1 (en) * | 2012-06-01 | 2013-12-05 | Wistron Corp. | Circuit layout method for printed circuit board, eletronic device and computer readable recording media |
| US20140282286A1 (en) * | 2013-03-15 | 2014-09-18 | Globalfoundries Singapore Pte. Ltd. | Etch failure prediction based on wafer resist top loss |
| CN105301601A (zh) * | 2015-10-09 | 2016-02-03 | 中国科学院光电研究院 | 一种适用于全球区域的gnss电离层延迟三维建模方法 |
| CN105738919A (zh) * | 2016-02-17 | 2016-07-06 | 东南大学 | 一种基于折半搜索算法的电离层穿刺点坐标计算方法 |
| CN106405589A (zh) * | 2016-06-24 | 2017-02-15 | 西安科技大学 | 一种全球电离层格网模型的确定方法及装置 |
| CN106405576A (zh) * | 2016-08-25 | 2017-02-15 | 中国科学院国家授时中心 | 一种基于gnss载波相位观测值的电离层监测与建模方法 |
| CN107622304A (zh) * | 2017-08-21 | 2018-01-23 | 东南大学 | 一种基于bp神经网络技术的电离层球谐函数改进方法 |
Non-Patent Citations (1)
| Title |
|---|
| 刘聪: "基于改进ＡＲＭＡ模型的电离层球谐模型系数预报", 《测绘地理信息》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110851976A (zh) * | 2019-11-07 | 2020-02-28 | 东南大学 | 一种聚焦离子束辅助沉积结构轮廓确定方法 |
| CN110851976B (zh) * | 2019-11-07 | 2023-08-25 | 东南大学 | 一种聚焦离子束辅助沉积结构轮廓确定方法 |
| CN111273335A (zh) * | 2019-12-20 | 2020-06-12 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | 一种基于垂测数据约束的电离层层析成像方法 |
| CN111273335B (zh) * | 2019-12-20 | 2021-09-17 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | 一种基于垂测数据约束的电离层层析成像方法 |
| CN111796309A (zh) * | 2020-06-24 | 2020-10-20 | 中国科学院精密测量科学与技术创新研究院 | 导航卫星单频数据同步确定大气水汽和总电子含量的方法 |
| CN111796309B (zh) * | 2020-06-24 | 2023-04-18 | 中国科学院精密测量科学与技术创新研究院 | 导航卫星单频数据同步确定大气水汽和总电子含量的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108491616B (zh) | 2022-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113093242B (zh) | 一种基于球谐展开的gnss单点定位方法 | |
| EP2746811B1 (en) | Methods for generating accuracy information on an ionosphere model for satellite navigation applications | |
| Farrell et al. | Differential GPS reference station algorithm-design and analysis | |
| Dettmering et al. | Combination of different space-geodetic observations for regional ionosphere modeling | |
| Orús et al. | Performance of different TEC models to provide GPS ionospheric corrections | |
| CN104102822B (zh) | 一种多频gnss观测值随机特性建模方法 | |
| CN110568459B (zh) | 基于igs和cors站的区域电离层tec实时监测方法 | |
| CN106405576A (zh) | 一种基于gnss载波相位观测值的电离层监测与建模方法 | |
| CN108491616B (zh) | 一种基于椭球谐函数理论的电离层垂直总电子含量建模方法 | |
| CN107622304A (zh) | 一种基于bp神经网络技术的电离层球谐函数改进方法 | |
| CN116594046B (zh) | 基于低轨卫星信号多普勒误差补偿的运动目标定位方法 | |
| CN111123345B (zh) | 一种基于gnss测量的经验电离层模型数据驱动方法 | |
| US20230126547A1 (en) | Machine learning in gnss receivers for improved velocity outputs | |
| Brack et al. | Operational multi-GNSS global ionosphere maps at GFZ derived from uncombined code and phase observations | |
| Raquet et al. | Use of a Covariance Analysis Technique for Predicting Performance of Regional‐Area Differential Code and Carrier‐Phase Networks | |
| Gaglione et al. | A stochastic sigma model for GLONASS satellite pseudorange | |
| CN102830406A (zh) | 一种gps天线相位中心变化对绝对定位精度的修正方法 | |
| Karimi et al. | Combination of swarm, Jason-3, and GNSS observations to construct a new modeling of global ionospheric maps | |
| CN114594435B (zh) | 一种国产民用sar卫星几何检校与定位精度改进方法 | |
| Chen et al. | Accuracy evaluation of XUST’s global ionospheric products | |
| Sharifi et al. | Local ionospheric modeling using the localized global ionospheric map and terrestrial GPS | |
| CN114527500A (zh) | 一种室内外一体化定位方法、设备、介质及产品 | |
| Zhilinskiy et al. | Study on the influence of ephemeris and clock information on GLONASS standard positioning service | |
| Madhu Krishna et al. | Improving the position accuracy of rover receiver using differential positioning in indian regional navigation satellite system | |
| Elshambaky | A novel approach to enhancing the Klobuchar algorithm to mitigate the effect of ionospheric delay errors on static single-frequency receivers |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |