CN114024632B - 地-各向异性电离层波导vlf波传播特性的获取方法 - Google Patents

地-各向异性电离层波导vlf波传播特性的获取方法 Download PDF

Info

Publication number
CN114024632B
CN114024632B CN202111288035.6A CN202111288035A CN114024632B CN 114024632 B CN114024632 B CN 114024632B CN 202111288035 A CN202111288035 A CN 202111288035A CN 114024632 B CN114024632 B CN 114024632B
Authority
CN
China
Prior art keywords
point
coordinate system
earth
rectangular coordinate
vlf
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
CN202111288035.6A
Other languages
English (en)
Other versions
CN114024632A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202111288035.6A priority Critical patent/CN114024632B/zh
Publication of CN114024632A publication Critical patent/CN114024632A/zh
Application granted granted Critical
Publication of CN114024632B publication Critical patent/CN114024632B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/309—Measuring or estimating channel quality parameters
    • 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • 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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • H04B17/3912—Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region
    • 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)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Electromagnetism (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Signal Processing (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种地‑各向异性电离层波导VLF波传播特性的获取方法，属于电子通信技术领域。本发明方法利用坐标变换和地球的椭球率，采用二分法求解过球面上两定点的大圆所满足的方程，确定电磁波的传播路径；还提出了一种新的电子密度模型，通过国际地磁参考场和坐标变换以求得传播路径上各向异性电离层的介电系数，基于有限元方法计算地‑各向异性电离层波导VLF波的传播特性，与实测结果对比，对该电子密度模型进行修正，最终得到地‑各向异性电离层波导VLF波的传播特性。通过计算不同频率VLF波在接收点的幅度和相位的传播特性以模拟原始信号，通过模拟原始信号及与实测结果的比对，实现电离层监测、地震监测、达到通信、导航目的。

Description

地-各向异性电离层波导VLF波传播特性的获取方法
技术领域
本发明属于电子通信技术领域，具体涉及一种模拟地-各向异性电离层波导VLF波传播特性的方法。
背景技术
甚低频(Very Low Frequency，VLF)电磁波(3～30kHz)在通信和导航、授时、电离层监测、地质勘探和地震预测等方面有着广泛的应用前景。VLF电磁波在地球-电离层波导中传播时损耗小，幅度与相位稳定，并能穿透一定深度的海水和土壤。除了实验手段，常见的计算地-电离层波导VLF波传播的方法包括波导模理论方法和时域有限差分方法。与数值FDTD方法相比，频域有限元(Finite element method，FEM)方法可以采用非均匀三角形网格，边界适应性强，具有求解速度快等优点，迄今为止尚未见基于有限元法地-电离层波导的波传播特性研究。此外，在现有的解析解和数值仿真中，传播路径的确定多基于发射点和接收点满足的球面三角(或精度为3°的多个地理方位角)，将地球等效为一均匀圆球，未引入椭球率；尽管有文献基于波导模理论研究了地磁场引起电离层的各向异性，但均采用简化模型，未考虑地磁场沿传播路径的变化，这些因素一定程度上会影响地-电离层波导结构模型建立的准确性。
发明内容
本发明是针对背景技术存在的缺陷，提供一种地-各向异性电离层波导VLF波传播特性的获取方法。
本发明采用的技术方案如下：
一种地-各向异性电离层波导VLF波传播特性的获取方法，包括以下步骤：
S1已知VLF波的发射点和接收点的经纬度，以及VLF波实际传播路径上任意一个待求点的经度L3，计算在地电离层波导结构直角坐标系下VLF波实际传播路径所在的大圆路径。
S2计算VLF波实际传播路径上任意一点在任意时刻的电离层D区的电子密度Ne(h)。
S3选取VLF波实际传播路径上的若干预测点，计算所有预测点的地磁要素。
S4计算VLF波实际传播路径上所有预测点的各向异性的相对介电系数。
S5基于有限元方法进行数值建模仿真，获得接收点VLF波的幅度和相位的传播特性。
本发明方法利用坐标变换和地球的椭球率，采用二分法求解过球面上两定点的大圆所满足的方程，确定电磁波的传播路径；还提出了一种新的电子密度模型，通过国际地磁参考场和坐标变换以求得传播路径上各向异性电离层的相对介电系数，基于有限元方法计算地-各向异性电离层波导VLF波的传播特性，与实测结果对比，对该电子密度模型进行修正，最终得到地-各向异性电离层波导VLF波的传播特性。通过计算不同频率VLF波在接收点的幅度和相位的传播特性以模拟原始信号，通过模拟原始信号及与实测结果的比对，实现电离层监测、地震监测、达到通信、导航目的。
附图说明
图1为地球直角坐标系-地电离层波导结构直角坐标系的变换关系示意图；
图2为二维直角坐标系下地-电离层波导模型；
图3为地电离层波导结构直角坐标系和地磁场所在球坐标系的转换示意图。
具体实施方式
以下结合附图以及具体实施例对本发明进行详细说明。
本实施例通过计算获得的地-各向异性电离层波导VLF波传播特性，实现电离层监测、通信和导航的目的。
S1已知VLF波的发射点和接收点的经纬度，以及VLF波实际传播路径上任意一个待求点的经度L3，计算在地电离层波导结构直角坐标系下VLF波实际传播路径所在的大圆路径。
S1.1地球椭球坐标系-地球直角坐标系的转换。
在地球椭球坐标系下，已知发射点的经度L1和纬度u1，接收点的经度L2和纬度u2，利用公式(1)的地球椭球坐标系和地球直角坐标系的转换关系，得到地球直角坐标系下发射点的坐标(x1，y1，z1)和接收点的坐标(x2，y2，z2)。
xn＝acosun cosLn
yn＝acosun sinLn (1)
zn＝bsinun
其中，xn、yn和zn分别为空间中任意一点位置在地球直角坐标系下的三个分量，下标n＝1、2、3...，a和b分别为地球的赤道半径和极半径。
S1.2已知VLF波实际传播路径上任意一个待求点的经度L3，计算该待求点在地球直角坐标系下的坐标(x3，y3，z3)。
通过接收点、发射点和待求点的空间平面方程表示为：
f＝a1xn+b1yn+c1zn＝0 (2)
其中f为空间平面方程，系数a1、b1和c1满足：
a1＝y1z2-y2z1,b1＝x2z1-x1z2,c1＝x1y2-x2y1
已知实际传播路径上任意一个待求点的经度L3，采用二分法求解式(1)和式(2)，得到待求点的纬度u3，然后代入式(1)得到待求点在地球直角坐标系下的坐标(x3，y3，z3)。
具体的，将经度L3代入式(1)和式(2)，通过计算发射点和接收点纬度区间范围内中点对应的函数值f，并检验其正负号，将有根的纬度区间范围作为新的纬度区间范围；迭代计算新的纬度区间范围内中点对应的函数值f，直到纬度区间范围的宽度达到精度要求，将纬度区间的中点作为待求点的纬度u3。本实施例采用的二分法具有计算简单可靠，能够得到收敛的精确解的优势。
S1.3计算待求点在地电离层波导结构直角坐标系下的坐标(X″3,Y″3)，从而得到在地电离层波导结构直角坐标系下VLF波实际传播路径所在的大圆路径。
利用有限元方法进行建模时，需要将VLF波的实际传播路径投影到地电离层波导结构直角坐标系，因此需要进行相应的坐标变换。如图1所示，将地球直角坐标系oxyz先绕y轴旋转β角，再绕x轴旋转-α角，即可得到地电离层波导结构直角坐标系O″X″Y″Z″。具体的，地球直角坐标系-地电离层波导结构直角坐标系的变换关系为：
Figure GDA0003890503740000051
其中R为转换矩阵，X″n、Y″n、Z″n分别为传播路径上任意一点在地电离层波导结构直角坐标系下的三个分量，下标n＝1、2、3...，α和β为坐标变换的旋转角度。
令Z″n＝0，分别将(x1，y1，z1)和(x2，y2，z2)代入式(3)，即可求出α和β。
再代入地球直角坐标系下发射点的坐标(x1，y1，z1)、接收点的坐标(x2，y2，z2)、待求点的坐标(x3，y3，z3)，得到地电离层波导结构直角坐标系下的发射点的坐标(X″1,Y″1)、接收点的坐标(X″2,Y″2)、待求点的坐标(X″3,Y″3)，从而得到在地电离层波导结构直角坐标系下VLF波实际传播路径所在的大圆路径。
S2计算VLF波实际传播路径上任意一点在任意时刻的电离层D区的电子密度Ne(h)
在地电离层波导结构直角坐标系下，VLF波实际传播路径上任意一点的电离层D区的电子密度Ne(h)可以采用现有的电子密度指数模型或国际电离层参考模型(International Reference Ionosphere，IRI)计算得到。
其中，现有的电离层电子密度指数模型可以采用国际无线电咨询委员会(International Radio Consultative Committee，CCIR)推荐的模型：
Ne(h)＝1.43×107e-0.15He(β1-0.15)(h-H) (4)
其中Ne(h)表示离地高度h的电子密度，β1表示电子浓度随高度变化的梯度参数，单位km-1；H为低电离层距地面的参考高度，单位为km。
本实施例还提出了一种新的适用于低电离层D区的电子密度模型：
Figure GDA0003890503740000061
Figure GDA0003890503740000062
其中x为天顶角，系数R0、A、y0、A1、A2、x0、d随太阳天顶角、离地高度、地理位置经纬度和季节的变化而变化，这些系数的初值根据电子密度IRI模型提供的某一天内随高度和传播路径变化的电子密度数据，结合传播路径上该天24小时内太阳天顶角的变化，采用数据拟合的方式确定，从而得到初始的电子密度模型。
S3选取VLF波实际传播路径上的若干预测点，计算所有预测点的地磁要素
IGRF(International Geomagnetic Reference Field，国际地磁参考场)作为国际通用的全球地磁标准模型，描述了地球的主磁场及其长期变化。
在地电离层波导结构直角坐标系下，在电离层D区选取至少30层不同的离地高度，同一离地高度下选取VLF波实际传播路径的至少30个点作为预测点；然后根据IGRF全球地磁标准模型计算所有预测点的磁场三个分量、磁倾角、磁偏角这些地磁要素。
S4计算VLF波实际传播路径上所有预测点的相对介电系数
从地电离层波导结构直角坐标系到地磁场所在球坐标系的转换关系为：
Figure GDA0003890503740000071
其中Q为转换矩阵，
Figure GDA0003890503740000072
Aθ、Ar为地磁场所在球坐标系的坐标或电磁场矢量的三个分量，AX、AY、AZ为地-电离层波导结构直角坐标系的坐标或电磁场矢量的三个分量，θi为磁倾角，θd为磁偏角。
在地磁场所在球坐标系中，预测点地磁场H0的方向与球坐标系中坐标或电磁场Ar方向重合，电离层D区中各向异性的相对介电系数张量表达式为：
Figure GDA0003890503740000073
其中
Figure GDA0003890503740000074
Figure GDA0003890503740000075
Figure GDA0003890503740000076
ε′,ε″,q为计算过程的中间变量，没有具体物理意义；ω为工作角频率，ω0为等离子体频率，ωT为磁旋频率，ν为等离子体碰撞频率，ε0为真空介电常数，i为虚数单位。
在电离层D区中，离地高度h的等离子体频率满足：
Figure GDA0003890503740000081
Ne(h)为离地高度h的电子密度，m为电子质量。
磁旋频率满足：
ωT＝μ0eH0/m (11)
e为电子的电量，H0为预测点的地磁场，μ0为真空磁导率。
得到在地电离层波导结构直角坐标系O″X″Y″Z″下，电离层D区各预测点的各向异性的相对介电系数：
Figure GDA0003890503740000082
将本实施例提出的适用于电离层D区的电子密度模型、或现有电离层电子密度指数模型或IRI模型Ne(h)代入式(10)，国际地磁参考场IGRF得到的磁场强度代入式(11)，磁倾角和磁偏角代入式(11)，然后代入式(12)，即可得到电离层D区各预测点的各向异性的相对介电系数。
S5基于有限元方法进行数值建模仿真，获得接收点VLF波的幅度和相位。
基于有限元方法，根据待求点的经度和纬度，获取大圆路径上的横纵坐标进行地层、空气层和电离层的建模，设置工作频率、电偶极子源、散射边界条件、发射点和接收点的位置，引入随横纵坐标变换的电子密度，设置各向异性的相对介电系数，然后进行网格划分，对矩阵方程进行物理求解，进而获得接收点VLF波的幅度和相位。
S6优化适用于低电离层D区的电子密度模型
将基于有限元方法得到的接收点VLF波幅度和相位与实测结果进行对比，对本实施例步骤S2提出的适用于低电离层D区的初始电子密度模型进行参数修正，得到新的随离地高度、太阳天顶角、地理位置和季节变化的电子密度模型，迭代优化直到VLF波传播的仿真结果和实测结果在一天24小时内的吻合精度满足要求，确定这些系数的最终值，得到最终的适用于低电离层D区的电子密度模型。
通过计算不同频率VLF电磁波幅度和相位的传播特性以模拟原始信号，通过模拟原始信号及与实测结果的比对，实现电离层监测、地震监测、达到通信、导航等目的。

Claims (1)

1.一种地-各向异性电离层波导VLF波传播特性的获取方法，包括以下步骤：
S1已知VLF波的发射点和接收点的经纬度，以及VLF波实际传播路径上任意一个待求点的经度L3，计算在地电离层波导结构直角坐标系下VLF波实际传播路径所在的大圆路径；
S2计算VLF波实际传播路径上任意一点在任意时刻的电离层D区的电子密度Ne(h)；
S3选取VLF波实际传播路径上的若干预测点，计算所有预测点的地磁要素；
S4计算VLF波实际传播路径上所有预测点的各向异性的相对介电系数；
S5基于有限元方法进行数值建模仿真，获得接收点VLF波的幅度和相位的传播特性；
具体地，步骤S1中，计算在地电离层波导结构直角坐标系下VLF波实际传播路径所在的大圆路径，包括以下步骤：
S1.1地球椭球坐标系-地球直角坐标系的转换；
在地球椭球坐标系下，已知发射点的经度L1和纬度u1，接收点的经度L2和纬度u2，利用公式(1)的地球椭球坐标系和地球直角坐标系的转换关系，得到地球直角坐标系下发射点的坐标(x1，y1，z1)和接收点的坐标(x2，y2，z2)：
Figure FDA0003890503730000011
其中，xn、yn和zn分别为空间中任意一点位置在地球直角坐标系下的三个分量，下标n＝1、2、3...，a和b分别为地球的赤道半径和极半径；
S1.2已知VLF波实际传播路径上任意一个待求点的经度L3，计算该待求点在地球直角坐标系下的坐标(x3，y3，z3)；
通过接收点、发射点和待求点的空间平面方程表示为：
f＝a1xn+b1yn+c1zn＝0 (2)
其中f为空间平面方程，系数a1、b1和c1满足：
a1＝y1z2-y2z1,b1＝x2z1-x1z2,c1＝x1y2-x2y1
已知实际传播路径上任意一个待求点的经度L3，采用二分法求解式(1)和式(2)，得到待求点的纬度u3，然后代入式(1)得到待求点在地球直角坐标系下的坐标(x3，y3，z3)；
S1.3计算待求点在地电离层波导结构直角坐标系下的坐标(X″3,Y″3)，从而得到在地电离层波导结构直角坐标系下VLF波实际传播路径所在的大圆路径；
将地球直角坐标系oxyz先绕y轴旋转β角，再绕x轴旋转-α角，得到地电离层波导结构直角坐标系O″X″Y″Z″；地球直角坐标系-地电离层波导结构直角坐标系的变换关系为：
Figure FDA0003890503730000021
其中R为转换矩阵，X″n、Y″n、Z″n分别为传播路径上任意一点在地电离层波导结构直角坐标系下的三个分量，下标n＝1、2、3...，α和β为坐标变换的旋转角度；
令Z″n＝0，分别将(x1，y1，z1)和(x2，y2，z2)代入式(3)，求出α和β；
再代入地球直角坐标系下发射点的坐标(x1，y1，z1)、接收点的坐标(x2，y2，z2)、待求点的坐标(x3，y3，z3)，得到地电离层波导结构直角坐标系下的发射点的坐标(X″1,Y″1)、接收点的坐标(X″2,Y″2)、待求点的坐标(X″3,Y″3)，从而得到在地电离层波导结构直角坐标系下VLF波实际传播路径所在的大圆路径；
步骤S4中，计算VLF波实际传播路径上所有预测点的相对介电系数，包括以下步骤：
从地电离层波导结构直角坐标系到地磁场所在球坐标系的转换关系为：
Figure FDA0003890503730000031
其中Q为转换矩阵，
Figure FDA0003890503730000032
Aθ、Ar为地磁场所在球坐标系的坐标或电磁场矢量的三个分量，AX、AY、AZ为地-电离层波导结构直角坐标系的坐标或电磁场矢量的三个分量，θi为磁倾角，θd为磁偏角；
在地磁场所在球坐标系中，预测点地磁场H0的方向与球坐标系中坐标或电磁场Ar方向重合，电离层D区中各向异性的相对介电系数张量表达式为：
Figure FDA0003890503730000033
其中
Figure FDA0003890503730000034
Figure FDA0003890503730000035
Figure FDA0003890503730000036
ε′、ε″、q为计算过程的中间变量，没有具体物理意义；ω为工作角频率，ω0为等离子体频率，ωT为磁旋频率，ν为等离子体碰撞频率，ε0为真空介电常数，i为虚数单位；
在电离层D区中，离地高度h的等离子体频率满足：
Figure FDA0003890503730000037
Ne(h)为离地高度h的电子密度，m为电子质量；
磁旋频率满足：
ωT＝μ0eH0/m (11)
e为电子的电量，H0为预测点的地磁场，μ0为真空磁导率；
得到在地电离层波导结构直角坐标系O″X″Y″Z″下，电离层D区各预测点的各向异性的相对介电系数：
Figure FDA0003890503730000041
步骤S5中，获得接收点VLF波的幅度和相位的传播特性的方式为：基于有限元方法，根据待求点的经度和纬度，获取大圆路径上的横纵坐标进行地层、空气层和电离层的建模，设置工作频率、电偶极子源、散射边界条件、发射点和接收点的位置，引入随横纵坐标变换的电子密度，设置各向异性的相对介电系数，然后进行网格划分，对矩阵方程进行物理求解，进而获得接收点VLF波的幅度和相位。
CN202111288035.6A 2021-11-02 2021-11-02 地-各向异性电离层波导vlf波传播特性的获取方法 Active CN114024632B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111288035.6A CN114024632B (zh) 2021-11-02 2021-11-02 地-各向异性电离层波导vlf波传播特性的获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111288035.6A CN114024632B (zh) 2021-11-02 2021-11-02 地-各向异性电离层波导vlf波传播特性的获取方法

Publications (2)

Publication Number Publication Date
CN114024632A CN114024632A (zh) 2022-02-08
CN114024632B true CN114024632B (zh) 2023-03-28

Family

ID=80059706

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111288035.6A Active CN114024632B (zh) 2021-11-02 2021-11-02 地-各向异性电离层波导vlf波传播特性的获取方法

Country Status (1)

Country Link
CN (1) CN114024632B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114741839B (zh) * 2022-03-02 2024-04-30 西北工业大学 一种分析甚低频电磁波在地-电离层中传播的fdtd方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109858102A (zh) * 2019-01-04 2019-06-07 西安理工大学 一种结合iri模型的甚低频电波传播时变特性预测方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5759676B2 (ja) * 2010-03-17 2015-08-05 株式会社東芝 伝搬経路推定システム及び伝搬経路推定方法
US10852333B2 (en) * 2017-12-11 2020-12-01 X Development Llc Earth-ionosphere waveguide power transfer
US10739450B2 (en) * 2017-12-11 2020-08-11 X Development Llc Using phased arrays of waveguide couplers in waveguide cavities to detect loading
CN111984913A (zh) * 2020-07-10 2020-11-24 西安理工大学 实际地层和电离层情况下的vlf模方程根的求解方法
CN112036011B (zh) * 2020-08-05 2022-09-20 中国人民解放军海军工程大学 一种用于水下航行器的甚低频波通信传输分析方法及系统

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109858102A (zh) * 2019-01-04 2019-06-07 西安理工大学 一种结合iri模型的甚低频电波传播时变特性预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"地-电离层"波导结构中媒质变化对VLF电波场强的影响;蒲玉蓉等;《西安理工大学学报》;20191230(第04期);全文 *

Also Published As

Publication number Publication date
CN114024632A (zh) 2022-02-08

Similar Documents

Publication Publication Date Title
Rietveld et al. New capabilities of the upgraded EISCAT high‐power HF facility
Yu et al. A 3-D global earth-ionosphere FDTD model including an anisotropic magnetized plasma ionosphere
CN109858102B (zh) 一种结合iri模型的甚低频电波传播时变特性预测方法
Liu et al. Discrimination between dry and water ices by full polarimetric radar: Implications for China’s first martian exploration
Vertogradov et al. Chirp ionosonde–radiodirection finder as a new tool for studying the ionosphere and radio-wave propagation
CN114024632B (zh) 地-各向异性电离层波导vlf波传播特性的获取方法
CN101719802A (zh) 一种预测短波通信最大可用频率的装置及计算方法
CN114167365A (zh) 导轨式外场目标、环境雷达特征获取方法及系统
Zhang et al. Simulation of Martian near-surface structure and imaging of future GPR data from Mars
Gasdia et al. A new longwave mode propagator for the Earth–ionosphere waveguide
CN110414182B (zh) 引入天线方向图的探地雷达frtm算法
CN107491611A (zh) 一种基于微波暗室天线布局的抗干扰测试场景设计方法
CN112036011B (zh) 一种用于水下航行器的甚低频波通信传输分析方法及系统
Carrano et al. Wave-optics analysis of HF propagation through traveling ionospheric disturbances and developing plasma bubbles
Niknam et al. A review of grid-based, time-domain modeling of electromagnetic wave propagation involving the ionosphere
Zhang et al. The spatial distribution features of three Alpha transmitter signals at the topside ionosphere
CN114167505B (zh) 基于LoranC甚低频信号的低电离层D层探测系统及方法
Zhang et al. A detection performance analysis of Sanya incoherent scatter radar tristatic system
Danilogorskaya et al. On the determination of the effect of horizontal ionospheric gradients on ranging errors in GNSS positioning
Zhou et al. Comparisons of ionospheric electron density distributions reconstructed by GPS computerized tomography, backscatter ionograms, and vertical ionograms
Stathacopoulos et al. Geometric considerations in the design of communications circuits using field-aligned ionospheric scatter
Barabashov et al. Mathematical space-time model of a sky wave radio field
Ma et al. Simulation research on the influence of plasma bubbles on radio wave propagation
CN109490967A (zh) 一种用于由卫星磁异常反演居里等温面埋深的方法
Li et al. Intelligent Frequency Selection of the Sky-Wave Radar Based on Numerical Ray Tracing

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