CN113569433B - 一种太阳闪烁下深空电磁波信道传播建模方法 - Google Patents
一种太阳闪烁下深空电磁波信道传播建模方法 Download PDFInfo
- Publication number
- CN113569433B CN113569433B CN202111117746.7A CN202111117746A CN113569433B CN 113569433 B CN113569433 B CN 113569433B CN 202111117746 A CN202111117746 A CN 202111117746A CN 113569433 B CN113569433 B CN 113569433B
- Authority
- CN
- China
- Prior art keywords
- electromagnetic wave
- sun
- solar
- deep space
- space electromagnetic
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Plasma Technology (AREA)
Abstract
本发明公开了一种太阳闪烁下深空电磁波信道传播建模方法,该方法根据深空电磁波的传播路径,考虑了太阳、行星、电磁波传播的几何关系的影响,建立了受不同程度太阳风影响下的深空电磁波信道传播模型,引入太阳活动周期及空间天气对深空电磁波传播的影响。相比传统方法,本方法能够应用于复杂太阳活动情况下描述太阳闪烁对深空电磁波在太阳系行星间传播的影响,为深空探测信息有效传输提供良好的技术支撑。
Description
技术领域
本发明属于深空探测领域,具体涉及一种太阳闪烁下深空电磁波信道传播建模方法。
背景技术
深空电磁波传播是深空通信和深空探测的主要支撑技术手段,深空电磁波在宇宙空间传播时会受到太阳风中等离子体的影响,具体表现为太阳闪烁。电磁波在太阳闪烁的影响下产生幅度衰减、相位抖动、法拉第旋转等效应,在严重情况下甚至导致深空通信链路的中断。通常,太阳闪烁受深空电磁波传播几何路径、太阳风速度、太阳风中等离子体电子密度等因素的影响,其变化较复杂,并且受限于人类在深空探测过程所获得的电磁波信号传播特征的有限性,目前对太阳闪烁与深空电磁波相互作用机理的理解尚处于探索阶段。此外,目前描述衰减信道的一些模型如莱斯模型、瑞利模型等并不能有效刻画电磁波在太阳风等离子体中的传播特征,而Nakagami-m模型也仅适用于较弱的闪烁环境,其普适性和广域有效性尚未得以全面验证。因此有必要全面认知深空电磁波传播的太阳闪烁效应,构建准确可信的信道传播模型,为深空探测,特别是太阳系内行星探测器的遥测遥控通信提供良好的理论基础和技术支撑。
目前研究已较为广泛的电离层闪烁强度,其中一个衡量标准即电离层总电子含量变化指数,Rate of TEC Index(ROTI)。ROTI指数最早由美国NASA JPL实验室的研究人员PiXiaoqing提出,参见文献1,其于1997年在地球科学领域顶级期刊《地球物理快报》所发表的文章目前已被广泛引用(Web of Science 引用469次),此后相关领域学者均以ROTI指数作为电离层闪烁强弱的一个重要标准,即通过ROTI指数与设定门限比较,判断电离层闪烁强弱。
文献1:Pi, X., Mannucci, A. J., Lindqwister, U. J., & Ho, C. M.Monitoring of global ionospheric irregularities using the worldwide GPSnetwork. Geophysical Research Letters, 1997,24(18):2283-2286。
发明内容
本发明要解决的技术问题是:克服现有技术的不足,提供一种太阳闪烁下深空电磁波信道传播建模方法,该方法根据深空电磁波的传播路径,考虑了太阳、行星、电磁波传播的几何关系的影响,建立了受不同程度太阳风影响下的深空电磁波信道传播模型,引入太阳活动周期及空间天气对深空电磁波传播的影响。相比传统方法,本方法能够应用于复杂太阳活动情况下描述太阳闪烁对深空电磁波在太阳系行星间传播的影响,为深空探测信息有效传输提供良好的技术支撑。
本发明的技术方案如下:一种太阳闪烁下深空电磁波信道传播建模方法,具体实现步骤为:
步骤(1)确定深空电磁波传播信道的几何特征,即太阳、地球、行星之间的相对位置,判断行星与太阳的位置关系是上合、下合还是其它位置关系;其中,太阳、行星探测器、地球的上合几何关系为行星探测器和地球分别处于太阳两侧,太阳、行星探测器、地球的下合几何关系为行星探测器和地球位于太阳同侧;上合时太阳闪烁影响大,下合时太阳闪烁影响小。进一步考虑太阳、地球、行星探测器的夹角α SEP ,其为太阳地球连线与地球行星探测器连线之间的夹角,且α SEP 越小,太阳闪烁影响越大。
步骤(2)确定电磁波传播时的空间环境特征,太阳活动所处的周期状态,是否存在太阳耀斑和日冕抛射物活动;空间环境特征主要考虑太阳活动状态,用F10.7指数、太阳黑子数、太阳耀斑和太阳日冕抛射物的状态来描述;太阳活动所处的周期状态主要指太阳活动在太阳活动周期中属于高年、低年、最大年或者平静期。
步骤(3)将太阳活动状态输入WSA-ENLIL太阳风模型,通过WSA-ENLIL太阳风模型计算太阳风速度、太阳风等离子体的密度和温度。
步骤(4)利用步骤(3)计算得到的太阳风等离子体密度计算深空电磁波传播路径上的总电子含量;具体计算方法为:
步骤(5)根据所计算得到的总电子含量进一步计算总电子含量的变化指数,并根据总电子含量的变化指数判断太阳闪烁强度;首先计算相邻时间间隔的总电子含量变化:,其中TEC i+1 、TEC i 分别表示第i+1和第i时刻的总电子含量,∆t i 表示第i+1和第i时刻的时间间隔,ROT表示总电子含量变化。进而计算总电子含量变化指数,具体方法为:
步骤(6)根据太阳闪烁强度分别进行深空电磁波信道传播建模,若属于弱太阳闪烁,则采用α-μ分布构建深空电磁波信道传播模型,若属于强太阳闪烁,则采用κ-μ分布构建深空电磁波信道传播模型;其中两种模型的具体描述为:
其中参数α和μ可以通过下式计算:
步骤(7)利用经验太阳风等离子体密度模型计算深空电磁波传播路径上的总电子含量,并进一步计算总电子含量的变化指数;其中经验太阳风等离子体电子密度模型为:
其中R 0为太阳半径,R为相对于太阳中心的径向距离,a、b、c、d为待定系数。N e 为太阳风等离子体密度,则总电子含量为:
步骤(8)利用步骤(7)中计算得到的总电子含量,重复步骤(5)和步骤(6),并设定确信因子θ 1和θ 2,且θ 1+θ 2=1。
步骤(9)利用实际观测数据拟合计算得到确信因子θ 1和θ 2的具体赋值,并确定最终的深空电磁波的太阳闪烁信道传播模型,可表示为:
其中,θ 1和θ 2为在步骤8和9中的确信因子,f WSA-ENLIL (x)为由步骤(3)中WSA-ENLIL太阳风模型计算得到的深空电磁波的太阳闪烁信道传播模型,f emperical (x)为由步骤(7)中经验模型得到的深空电磁波的太阳闪烁信道传播模型。
本发明与现有技术相比的优点在于:
(1)相比传统的计算方法,本发明方法(如图1)充分考虑了太阳活动对太阳闪烁的影响,综合了太阳风WSA-ENLIL理论模型和经验的太阳风等离子体电子密度计算方法,获得了更准确可信的深空电磁场太阳闪烁信道模型。
附图说明
图1为本发明一种太阳闪烁下深空电磁波信道传播建模方法实现流程图。
具体实施方式
下面将结合附图及具体实施方式对本发明加以详细说明,需要指出的是,所描述的实施例仅旨在便于对本发明的理解,而不起任何限定作用。
本发明提供了一种太阳闪烁下深空电磁波信道传播建模方法,该方法根据深空电磁波的传播路径,考虑了太阳、行星、电磁波传播的几何关系的影响,建立了受不同程度太阳风影响下的深空电磁波信道传播模型,引入太阳活动周期及空间天气对深空电磁波传播的影响。相比传统方法,本方法能够应用于复杂太阳活动情况下描述太阳闪烁对深空电磁波在太阳系行星间传播的影响,为深空探测信息有效传输提供良好的技术支撑。
如图1所示,本发明具体实现步骤如下:
1、确定深空电磁波传播信道的几何特征,即太阳、地球、行星之间的相对位置,判断行星与太阳的位置关系是上合、下合还是其它位置关系;其中,太阳、行星探测器、地球的上合几何关系为行星探测器和地球分别处于太阳两侧,太阳、行星探测器、地球的下合几何关系为行星探测器和地球位于太阳同侧;上合时太阳闪烁影响大,下合时太阳闪烁影响小。所述步骤A中,进一步考虑太阳、地球、行星探测器的夹角α SEP ,其为太阳地球连线与地球行星探测器连线之间的夹角,且α SEP 越小,太阳闪烁影响越大。
2、确定电磁波传播时的空间环境特征,太阳活动所处的周期状态,是否存在太阳耀斑和日冕抛射物活动;空间环境特征主要考虑太阳活动状态,用F10.7指数、太阳黑子数、太阳耀斑和太阳日冕抛射物的状态来描述;太阳活动所处的周期状态主要指太阳活动在太阳活动周期中属于高年、低年、最大年或者平静期。
3、将太阳活动状态输入WSA-ENLIL太阳风模型,通过WSA-ENLIL太阳风模型计算太阳风速度、太阳风等离子体的密度和温度。
4、利用上面一步计算得到的太阳风等离子体密度计算深空电磁波传播路径上的总电子含量;具体计算方法为:
5、根据所计算得到的总电子含量进一步计算总电子含量的变化指数,并根据总电子含量的变化指数判断太阳闪烁强度;首先计算相邻时间间隔的总电子含量变化:,其中TEC i+1 、TEC i 分别表示第i+1和第i时刻的总电子含量,Δt i 表示第i+1和第i时刻的时间间隔,ROT表示总电子含量变化。进而计算总电子含量变化指数,具体方法为:
6、根据太阳闪烁强度分别进行深空电磁波信道传播建模,若属于弱太阳闪烁,则采用α-μ分布构建深空电磁波信道传播模型,若属于强太阳闪烁,则采用κ-μ分布构建深空电磁波信道传播模型;其中两种模型的具体描述为:
α-μ模型:,其中参数α和μ可以通过下式计算:,其中。其中r是电磁波信号振幅向量的单位向量,是中间变量,可用来计算,Γ(·)表示Gamma函数,E(·)表示数学期望,R表示电磁波信号振幅向量,β为待定参数,通常可取值为β=3和β=4。
7、利用经验太阳风等离子体密度模型计算深空电磁波传播路径上的总电子含量,并进一步计算总电子含量的变化指数;其中经验太阳风等离子体电子密度模型为:
其中R 0为太阳半径,R为相对于太阳中心的径向距离,a、b、c、d为待定系数,此处赋值为a=2.39×1014,b=-6,c=1.67×1012,d=-2.5。N e 为太阳风等离子体密度,则总电子含量为:
8、利用上一步中计算得到的总电子含量,重复步骤5和步骤6,并设定确信因子θ 1和θ 2,且θ 1+θ 2=1。此处设定θ 1=0.7,θ 2=0.3。
9、利用实际观测数据拟合计算得到确信因子θ 1和θ 2的具体赋值,并确定最终的深空电磁波的太阳闪烁信道传播模型,可表示为:
其中,θ 1和θ 2为在步骤8和9中的确信因子,f WSA-ENLIL (x)为由步骤3中WSA-ENLIL太阳风模型计算得到的深空电磁波的太阳闪烁信道传播模型,f emperical (x)为由步骤7中经验模型得到的深空电磁波的太阳闪烁信道传播模型。
综述,本发明提供了一种太阳闪烁下深空电磁波信道传播建模方法,该方法根据深空电磁波的传播路径,考虑了太阳、行星、电磁波传播的几何关系的影响,建立了受不同程度太阳风影响下的深空电磁波信道传播模型,引入太阳活动周期及空间天气对深空电磁波传播的影响。相比传统方法,本方法能够应用于复杂太阳活动情况下描述太阳闪烁对深空电磁波在太阳系行星间传播的影响,为深空探测信息有效传输提供良好的技术支撑。
以上所述仅为本发明的具体实施例,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种太阳闪烁下深空电磁波信道传播建模方法,其特征在于,包括以下步骤:
步骤A、确定深空电磁波传播信道的几何特征,即太阳、地球、行星探测器之间的相对位置,判断行星探测器与太阳的位置关系是上合、下合还是其它位置关系;
步骤B、确定电磁波传播时的空间环境特征,包括:太阳活动所处的周期状态,是否存在太阳耀斑和日冕抛射物活动;
步骤C、根据步骤B中所确定的空间环境特征,将太阳活动状态输入WSA-ENLIL太阳风模型,进而计算太阳风速度、太阳风等离子体密度和温度;
步骤D、利用步骤C中计算得到的太阳风等离子体密度计算深空电磁波传播路径上的总电子含量;
步骤E、根据所计算得到的总电子含量进一步计算总电子含量的变化指数,并根据总电子含量的变化指数判断太阳闪烁强度;
步骤F、根据太阳闪烁强度分别进行深空电磁波信道传播建模,若属于弱太阳闪烁,则采用α-μ分布构建深空电磁波信道传播模型,若属于强太阳闪烁,则采用κ-μ分布构建深空电磁波信道传播模型;
步骤G、利用经验太阳风等离子体密度模型计算深空电磁波传播路径上的总电子含量,并进一步计算总电子含量的变化指数;
所述步骤G中,利用经验太阳风等离子体密度模型计算深空电磁波传播路径上的总电子含量的方法为:
其中R 0为太阳半径,R为相对于太阳中心的径向距离,a、b、c、d为待定系数;N e 为步骤C中计算得到的太阳风等离子体密度,则总电子含量为:
步骤H、利用步骤G中计算得到的总电子含量变化指数为输入,重复步骤F,并设定确信因子θ 1和θ 2,且θ 1+θ 2=1;
步骤I、利用实际观测数据拟合计算得到确信因子θ 1和θ 2的具体赋值,并确定最终的深空电磁波的太阳闪烁信道传播模型。
2.根据权利要求1所述的一种太阳闪烁下深空电磁波信道传播建模方法,其特征在于:所述步骤A中太阳、行星探测器、地球的上合几何关系为行星探测器和地球分别处于太阳两侧,太阳、行星探测器、地球的下合几何关系为行星探测器和地球位于太阳同侧;上合时太阳闪烁影响大,下合时太阳闪烁影响小;所述步骤A中,进一步考虑太阳、地球、行星探测器的夹角α SEP,其为太阳地球连线与地球行星探测器连线之间的夹角,且α SEP越小,太阳闪烁影响越大。
3.根据权利要求1所述的一种太阳闪烁下深空电磁波信道传播建模方法,其特征在于:所述步骤B中,空间环境特征包括考虑太阳活动状态,用F10.7指数、太阳黑子数、太阳耀斑和太阳日冕抛射物的状态来描述;所述步骤C中,将太阳活动状态输入WSA-ENLIL太阳风模型,通过WSA-ENLIL太阳风模型计算太阳风速度、太阳风等离子体的密度和温度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111117746.7A CN113569433B (zh) | 2021-09-24 | 2021-09-24 | 一种太阳闪烁下深空电磁波信道传播建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111117746.7A CN113569433B (zh) | 2021-09-24 | 2021-09-24 | 一种太阳闪烁下深空电磁波信道传播建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113569433A CN113569433A (zh) | 2021-10-29 |
CN113569433B true CN113569433B (zh) | 2022-03-25 |
Family
ID=78174105
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111117746.7A Active CN113569433B (zh) | 2021-09-24 | 2021-09-24 | 一种太阳闪烁下深空电磁波信道传播建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113569433B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105116469A (zh) * | 2015-08-05 | 2015-12-02 | 中国电子科技集团公司第二十二研究所 | 一种中国低纬地区电离层闪烁发生概率预报方法 |
CN113093225A (zh) * | 2021-02-18 | 2021-07-09 | 北京航空航天大学 | 一种广域、局域融合的高精度电离层闪烁模型建立方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8181511B2 (en) * | 2008-07-18 | 2012-05-22 | Meier Robert R | Method and system of imaging electrons in the near earth space environment |
CN101672922A (zh) * | 2009-09-23 | 2010-03-17 | 贵州师范大学 | 一种弱太阳闪烁下的深空探测信道预测方法 |
CN112512924B (zh) * | 2018-04-30 | 2022-08-16 | J·弗朗西斯三世·达金 | 从近地轨道清除轨道太空碎片 |
CN113359164B (zh) * | 2021-06-03 | 2021-12-03 | 中国矿业大学 | 一种低频率gnss电离层闪烁因子的有效性验证方法 |
-
2021
- 2021-09-24 CN CN202111117746.7A patent/CN113569433B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105116469A (zh) * | 2015-08-05 | 2015-12-02 | 中国电子科技集团公司第二十二研究所 | 一种中国低纬地区电离层闪烁发生概率预报方法 |
CN113093225A (zh) * | 2021-02-18 | 2021-07-09 | 北京航空航天大学 | 一种广域、局域融合的高精度电离层闪烁模型建立方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113569433A (zh) | 2021-10-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tsugawa et al. | A statistical study of large‐scale traveling ionospheric disturbances using the GPS network in Japan | |
Cane et al. | Interplanetary coronal mass ejections in the near‐Earth solar wind during 1996–2002 | |
Nguyen et al. | Evaluation of tropical cyclone center identification methods in numerical models | |
Assink et al. | Evaluation of wind and temperature profiles from ECMWF analysis on two hemispheres using volcanic infrasound | |
Wang et al. | The similarity between optimal precursor and optimally growing initial error in prediction of Kuroshio large meander and its application to targeted observation | |
Delgado‐Fernandez et al. | Event‐scale dynamics of a parabolic dune and its relevance for mesoscale evolution | |
Viljanen et al. | Climatology of rapid geomagnetic variations at high latitudes over two solar cycles | |
Menk et al. | Monitoring the plasmapause using geomagnetic field line resonances | |
Şentürk | Investigation of global ionospheric response of the severe geomagnetic storm on June 22-23, 2015 by GNSS-based TEC observations | |
Zhao et al. | Influence of Walker circulations on East African rainfall | |
Huang et al. | Simultaneous upward and downward propagating inertia‐gravity waves in the MLT observed at Andes Lidar Observatory | |
Zou et al. | Electrodynamics of the high‐latitude trough: Its relationship with convection flows and field‐aligned currents | |
CN107870336B (zh) | 一种星载激光雷达的可穿透层层底迭代检测方法 | |
Unnikrishnan et al. | A study on chaotic behaviour of equatorial/low latitude ionosphere over Indian subcontinent, using GPS-TEC time series | |
CN113569433B (zh) | 一种太阳闪烁下深空电磁波信道传播建模方法 | |
Wu et al. | Diurnal variations in tropical cyclone intensification | |
Shreedevi et al. | Geomagnetic storm‐induced plasma density enhancements in the southern polar ionospheric region: A comparative study using St. Patrick's Day storms of 2013 and 2015 | |
Carlson et al. | ESR mapping of polar‐cap patches in the dark cusp | |
Wang et al. | Human detection through RSSI processing with packet dropout in wireless sensor network | |
Uwamahoro et al. | Storm time total electron content modeling over African low‐latitude and midlatitude regions | |
Garniwa et al. | Intraday forecast of global horizontal irradiance using optical flow method and long short-term memory model | |
Baker et al. | Diurnal variability in aeolian sediment transport at Gale crater, Mars | |
Fišer et al. | Medium-scale traveling ionospheric disturbances over Taiwan observed with HF Doppler sounding | |
Yiğit et al. | Characterization of the thermospheric mean winds and circulation during solstice using ICON/MIGHTI observations | |
Hu et al. | Investigation of gust‐seismic relationships and applications to gust detection |
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 |