CN111008361B - 一种电离层参数重构方法 - Google Patents

一种电离层参数重构方法 Download PDF

Info

Publication number
CN111008361B
CN111008361B CN201911173417.7A CN201911173417A CN111008361B CN 111008361 B CN111008361 B CN 111008361B CN 201911173417 A CN201911173417 A CN 201911173417A CN 111008361 B CN111008361 B CN 111008361B
Authority
CN
China
Prior art keywords
reconstruction
ionospheric
phi
observation
ionosphere
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
Application number
CN201911173417.7A
Other languages
English (en)
Other versions
CN111008361A (zh
Inventor
王健
程华叶
黄翔东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201911173417.7A priority Critical patent/CN111008361B/zh
Publication of CN111008361A publication Critical patent/CN111008361A/zh
Application granted granted Critical
Publication of CN111008361B publication Critical patent/CN111008361B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明具体公开一种电离层参数重构方法，包括以下步骤：基于观测站点以及重构位置的地磁坐标，利用下式计算观测站点间及其与重构位置间的电离层距离；基于观测站点电离层距离矩阵以及重构位置与观测站点电离层距离向量，利用克里格重构方程组确定重构权重；基于观测站点所得的电离层参数值，利用线性加权方法基于所确定的重构权重，计算重构位置的电离层参数值。本发明建立的区域电离层参数重构方法，能够满足电离层研究的需求，计算复杂度较低，且能提高了电离层参数的重构精度。

Description

一种电离层参数重构方法
技术领域
本发明涉及电离层技术领域，特别是涉及一种区域电离层参数重构方法。
背景技术
电离层是近地空间环境重要组成部分，约是地面60公里以上到磁层顶之间的整个空间。它是人类生存的近地空间环境的一部分，充分了解电离层，是人类对自身生存环境认识和利用的重要基础。
电离层参数，如临界频率、传输因子等，存在全球空域变化。影响电离层的因素众多，既有来自太阳，行星际和磁层空间的环境变化，也包括热层的中性成分组成和浓度的变化和低层大气的影响。随着科技的进步，人类越来越依赖于现代的通讯和导航系统。有许多的通讯系统利用电离层反射电波信号进行远距离信息传播。更直接地，在中高层大气、电离层这些太空区域，正是人造卫星、航天飞机与空间站飞行的地方，而它们也恰是灾害性空间天气的直接发生地。因此，建立有效的区域化电离层参数模型，不仅对无线电通信有着重要的作用，而且在无线电定位和超视距雷达的应用中也是很重要的。
电离层的参数重构是对已知的部分探测站数据，实现电离层参数局部区域连续的实时分布，为相关科学的理论研究和工程应用的发展提供技术支撑。电离层参数的重构是目前国内外的研究热点，研究动力源于更高精度的需求，以保证对空间天气的感知、电波环境预报、警报和效应评估能力的准确评估。
发明内容
本发明的目的是针对现有技术中存在的技术缺陷，而提供一种区域电离层参数重构方法，是一种基于地磁坐标的区域电离层参数重构方法，用于实现对电离层参数的高精度重构。
为实现本发明的目的所采用的技术方案是：
一种区域电离层参数重构方法，包括以下步骤：
基于观测站点以及重构位置的地磁坐标，利用下式计算观测站点间及其与重构位置间的电离层距离；
d＝ae·acos{cosλAcosλBcos[SF·(φA-φB)]+sinφAsinφB}
式中，d表示电离层距离，ae为地球半径，取6371.2，(φA，λA)和(φB，λB)分别为A、B两点的地磁经纬度，SF为经验尺度因子，取1.2，
基于观测站点(φn，λn)电离层距离矩阵D以及重构位置(φ0，λ0)与观测站点(φn，λn)电离层距离向量D0，利用克里格重构方程组确定重构权重；
基于观测站点所得的电离层参数值，利用线性加权方法基于所确定的重构权重，计算重构位置的电离层参数值；
其中，所述的利用克里格重构方程组确定重构权重的步骤如下：
构建观测站点(φn，λn)电离层距离矩阵D，
计算重构位置(φ0，λ0)与观测站点(φn，λn)电离层距离向量D0,
利用克里格重构方程确定重构位置(φ0，λ0)与观测站点(φn，λn)间的重构权值向量W，
W满足/>
其中，所述的计算重构位置的电离层参数值的计算公式如下：
式中，p0表示重构位置(φ0，λ0)的电离层参数值，pn为观测站点(φn，λn)的电离层参数。
本发明建立的区域电离层参数重构方法，能够满足电离层研究的需求，计算复杂度较低，且能提高了电离层参数的重构精度。
附图说明
图1是本发明的区域电离层参数重构方法的流程图。
具体实施方式
以下结合附图和具体实施例对本发明作进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明是利用地磁坐标对电离层参数进行重构，给出区域重构精度更高方法，可为电子信息系统的可用频率预测、通信效能保障提供技术支撑。
如图1所示，本发明区域电离层参数重构方法，包括以下步骤：
步骤A：读取电离层观测站获得的电离层参数。
具体电离层参数的观测方式可包括垂直探测、斜向探测等手段，其中的电离层参数包括临界频率、峰值高度、传输因子等。
步骤B：计算观测站点及重构位置的地磁坐标。
其中，步骤B的计算方法为行业内通知内容，在此不再赘述。
步骤C：计算观测站点间及其与重构位置间的电离层距离，构建变异函数。
上述的电离层距离的计算方法选用加权大圆距离计算方法，具体可表示为：
d＝ae·acos{cosλAcosλBcos[SF·(φA-φB)]+sinφAsinφB} (1)
式中，ae为地球半径，取6371.2，(φA，λA)和(φB，λB)分别为A、B两点的地磁经纬度，SF为经验尺度因子，取1.2。
其中的A和B可为观测站点，亦可为重构位置点。
步骤D：利用克里格重构方程组确定重构权重：
其中，所述的步骤D利用克里格重构方程组确定重构权重的具体步骤如下：
步骤D1：构建观测站点电离层距离矩阵D。对于N个观测站点(φn，λn)，n＝1,2,…,N，其电离层距离可表示为
式中，任意两点的电离层距离利用公式(1)得到。
步骤D2：计算重构位置(φ0，λ0)与观测站点的电离层距离向量D0。具体可表示为
式中，任意两点的电离层距离利用公式(1)得到。
步骤D3：利用克里格重构方程确定重构位置与N个观测站点间的权值
W满足：
步骤E：计算重构位置的电离层参数值：
基于观测站点所得的电离层参数值，利用线性加权方法计算重构位置的电离层参数值p0，即
式中，pn为观测站点(φn，λn)的电离层参数。
以上所述仅是本发明的优选实施方式，应当指出的是，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (1)

1.一种区域电离层参数重构方法，其特征在于，包括以下步骤：
基于观测站点以及重构位置的地磁坐标，利用下式计算观测站点间及其与重构位置间的电离层距离；
d＝ae·acos{cosλAcosλBcos[SF·(φA-φB)]+sinφAsinφB}
式中，d表示电离层距离，ae为地球半径，取6371.2，(φA，λA)和(φB，λB)分别为A、B两点的地磁经纬度，SF为经验尺度因子，取1.2；
基于观测站点(φn，λn)电离层距离矩阵D以及重构位置(φ0，λ0)与观测站点(φn，λn)电离层距离向量D0，利用克里格重构方程组确定重构权重；
基于观测站点所得的电离层参数值，利用线性加权方法基于所确定的重构权重，计算重构位置的电离层参数值；
所述的利用克里格重构方程组确定重构权重的步骤如下：
构建N个观测站点(φn，λn)电离层距离矩阵D，n＝1,2,…,N，
计算重构位置(φ0，λ0)与观测站点(φn，λn)电离层距离向量D0,
利用克里格重构方程确定重构位置(φ0，λ0)与观测站点(φn，λn)间的重构权值向量W，
所述的计算重构位置的电离层参数值的计算公式如下：
式中，p0表示重构位置(φ0，λ0)的电离层参数值，pn为观测站点(φn，λn)的电离层参数。
CN201911173417.7A 2019-11-26 2019-11-26 一种电离层参数重构方法 Active CN111008361B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911173417.7A CN111008361B (zh) 2019-11-26 2019-11-26 一种电离层参数重构方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911173417.7A CN111008361B (zh) 2019-11-26 2019-11-26 一种电离层参数重构方法

Publications (2)

Publication Number Publication Date
CN111008361A CN111008361A (zh) 2020-04-14
CN111008361B true CN111008361B (zh) 2023-08-18

Family

ID=70112748

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911173417.7A Active CN111008361B (zh) 2019-11-26 2019-11-26 一种电离层参数重构方法

Country Status (1)

Country Link
CN (1) CN111008361B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111735686A (zh) * 2020-05-16 2020-10-02 厦门大学 一种用于钢骨混凝土结构性能分析的基于备选模型的Kriging预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108462545A (zh) * 2018-01-29 2018-08-28 武汉小石科技有限公司 一种基于单接收站的电离层foF2参数重构方法
CN109490641A (zh) * 2019-01-05 2019-03-19 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种中纬地区偶发e层短波场强的计算方法
CN109683196A (zh) * 2018-11-15 2019-04-26 天津大学青岛海洋技术研究院 一种电离层与地震前兆相关量时空特性分析方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108462545A (zh) * 2018-01-29 2018-08-28 武汉小石科技有限公司 一种基于单接收站的电离层foF2参数重构方法
CN109683196A (zh) * 2018-11-15 2019-04-26 天津大学青岛海洋技术研究院 一种电离层与地震前兆相关量时空特性分析方法
CN109490641A (zh) * 2019-01-05 2019-03-19 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种中纬地区偶发e层短波场强的计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
冯晓哲.短波单站定位中的准抛物电离层参数实时修正.《 电子信息对抗技术》.2008,全文. *

Also Published As

Publication number Publication date
CN111008361A (zh) 2020-04-14

Similar Documents

Publication Publication Date Title
Petrie et al. A review of higher order ionospheric refraction effects on dual frequency GPS
Gulyaeva et al. GIM-TEC adaptive ionospheric weather assessment and forecast system
Goudarzi et al. Assessing horizontal positional accuracy of Google Earth imagery in the city of Montreal, Canada
CN101295014B (zh) 基于gnss的远距离高精度实时/快速定位方法和系统
Nayak et al. Modulation of local wind-waves at Kalpakkam from remote forcing effects of Southern Ocean swells
Wang et al. Introduction to GPS geodetic infrastructure for land subsidence monitoring in Houston, Texas, USA
CN102865852B (zh) 一种对数字高程模型数据高程校准方法及装置
CN104729486A (zh) 一种基于似大地水准面精化的无验潮水深测量方法
CN104101888A (zh) 一种适合于中国区域的gnss电离层延迟精确建模方法
CN105844427A (zh) 一种风暴潮灾害风险精细化评估的计算方法
CN111123345B (zh) 一种基于gnss测量的经验电离层模型数据驱动方法
CN104392113A (zh) 一种近海海面冷空气大风风速的估算方法
CN111008361B (zh) 一种电离层参数重构方法
Gu et al. Surface floating objects moving from the Pearl River Estuary to Hainan Island: An observational and model study
Loftis et al. Integrated ocean, earth, and atmospheric observations for resilience planning in Hampton roads, Virginia
CN103322965A (zh) 一种惯性导航系统横卯酉面曲率半径测量方法
Di et al. Study on the regional ASF prediction method based on the ordinary kriging interpolation
Janowski et al. The Analyzes of PDOP factors for a Zigbee ground–based augmentation systems
Torres et al. Simulation of storm surge in northeast coast of the US; a closer look at the wind forcing
Zheng et al. Preliminary inquiry on the linear relationship between the height of the station and the ground height error retrieved by GNSS-IR with low-cost smart electronic equipment
Menegbo Determination of orthometric elevations using GNSS-derived height with the EGM2008 geoid height model
Dawod et al. Accuracy assessment of the PGM17 global geopotential model: a case study of Egypt and Northeast Africa
Ifadis et al. Space to earth geodetic observations: approaching the atmospheric effect
He et al. Analysis of Sea Surface Multipath and Impact on GNSS Precision Positioning
Khalid Coastal Inundation Risk Assessment Using Multi-Mission Satellite Altimeter and Digital Elevation Model, MSc

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