一种适用于电离层暴期间短波链路的最高可用频率计算方法技术

本发明专利技术公开了一种适用于电离层暴期间短波链路的最高可用频率计算方法，包括如下步骤：步骤1、在全国范围内优选5条优质斜测链路作为背景链路；步骤2、利用国内14个电波观测站的垂测数据fof2、foE和M（3000）F2因子；步骤3、分别获取历史电离层暴期间5条背景链路中点的3个垂测数据和斜测链路最高可用频率MUF数据；步骤4、将上述步骤3中插值得到的电离层暴期间链路中点垂测数据作为ITU‑R P533建议短波最高可用频率计算公式中的输入。本发明专利技术所公开的适用于电离层暴期间短波链路的最高可用频率计算方法，利用此方法实现了电离层暴期间短波链路最高可用频率的获取，基于ITU‑R P.533建议短波链路最高可用频率计算公式，结合电离层暴期间的电离层状态信息，建立一种适用于电离层暴期间短波链路的最高可用频率计算方法。

【技术实现步骤摘要】
一种适用于电离层暴期间短波链路的最高可用频率计算方法
本专利技术涉及短波通信和雷达
，尤其涉及一种适用于电离层暴期间短波链路的最高可用频率计算方法。
技术介绍
从地面向上发射的高频无线电信号能被电离层反射，这是短波天波信号的主要传播方式。高频无线电信号在经由电离层反射时，存在一个与反射电子密度相关的最高可用频率(MUF)，如果其频率高于此值，无线电波将穿透电离层不再返回地面。通过对垂测电离图的判读，获取fof2、M(3000)F2因子等参数，即可获得3000km路径上观测站上方反射的MUF。同时，电离层斜向探测亦是通过电离层的反射实现两点之间的短波探测。可见，短波链路的MUF与电离层电子密度分布密切相关。近年来，电离层平静时MUF的统计特性及其预测技术已趋于成熟。其中，以国际电信联盟无线电通信组推荐的ITU-RP.533建议最具权威性。该建议以ITU-RP.1239建议给出的参考电离层为背景，通过射线路径分析和曲线拟合等方法，建立了短波链路参数预测方法。该建议适用于电离层平静期短波频率场强、可靠性和兼容性的预测，即可以预测F2层临界频率、3000km传播因子等参数的月中值。由此可见，ITU-RP.533建议仅适用于电离层平静期的短波链路参数预测。而电离层受化学、动力学和电动力学等多种机制的影响，表现出复杂的变化特性，当发生电离层暴时，由于电离层参数显著偏离其背景值，会导致短波链路MUF的相应变化，电离层暴期间短波链路MUF计算误差分布和累计误差分布如图1所示。可见，现有技术中对电离层暴期间电离层参数的变化特性和可预报性，以及其造成的短波链路参数变化，还缺乏规律性的认识。
技术实现思路
本专利技术所要解决的技术问题就是提供一种适用于电离层暴期间短波链路的最高可用频率计算方法。本专利技术采用如下技术方案：一种适用于电离层暴期间短波链路的最高可用频率计算方法，其改进之处在于，包括如下步骤：步骤1、在全国范围内优选5条优质斜测链路作为背景链路；步骤2、利用国内14个电波观测站的垂测数据fof2、foE和M(3000)F2因子，通过克令格插值算法获取5条背景链路中点位置的垂测数据fof2、foE和M(3000)F2因子，上述fof2指电离层F2层临界频率、foE指电离层E层临界频率；步骤3、分别获取历史电离层暴期间5条背景链路中点的3个垂测数据和斜测链路最高可用频率MUF数据，并对背景链路数据进行预处理，作为下述步骤4中的训练样本；步骤4、将上述步骤3中插值得到的电离层暴期间链路中点垂测数据作为ITU-RP533建议短波最高可用频率计算公式中的输入，将斜测数据作为输出，利用粒子群算法优化公式中的20个常数参数，设定其优化次数或优化误差，以获得最优常数参数；步骤5、将获得的20个最优常数参数替换原公式中的常数参数，以国内14个电波观测站的垂测数据插值获得指定链路中点的垂测数据作为输入，建立一种适用于电离层暴期间短波链路的最高可用频率计算方法。进一步的，所述步骤3具体包括如下步骤：步骤31、将fof2相对偏差df作为电离层暴的判定标准，其中，式中，fof2为观测值，fof2med为月中值，当df不小于0.3时即认为是正相电离层暴，当df不大于-0.2时即为负相电离层暴；步骤32、利用国内14个电波观测站的垂测数据和克令格插值算法，获取某一链路中点3个垂测数据，并根据以上电离层暴的判定标准，分别筛选某一时间段内电离层暴期间链路中点的fof2、foE和M(3000)F2因子3个垂测数据和该斜测链路的最高可用频率数据；步骤33、对于同一时刻收发点互换的两条链路MUF的差值不超过3MHz，即认为该组数据质量可靠，可作为样本数据。进一步的，所述步骤4具体包括如下步骤：步骤41、设置提取ITU-RP.533建议短波链路最高可用频率计算公式中的20个常数参数，对于大圆距离小于4000km的短波链路一跳MUF，ITU-RP.533建议提供的链路参数预测方法如下所示：其中：fH：控制点处的电子磁旋频率；Cd＝X1-X2·Z-X3·Z2-X4·Z3+X5·Z4+X6·Z5+X7·Z6；其中，Z＝1-2d/dmax；dmax＝X8+(X9+X10/x2-X11/x4+X12/x6)(1/B-X13)；其中：d：传播路径长度；C3000：3000km时的Cd值；x＝foF2/foE或X20，取较大者；其中，X1至X20即为预先通过粒子群算法获取的常数参数；步骤42、利用粒子群算法获取最优20个常数参数，利用步骤2中插值获取的电离层暴期间链路中点垂测和斜测数据训练样本，并通过离子群算法对20个常数参数进行优化，设置优化次数或误差阈值，最终获得20个最优常数参数如下：本专利技术的有益效果是：本专利技术所公开的适用于电离层暴期间短波链路的最高可用频率计算方法，利用此方法实现了电离层暴期间短波链路最高可用频率的获取，基于ITU-RP.533建议短波链路最高可用频率计算公式，结合电离层暴期间的电离层状态信息，将电离层垂测数据插值获得的链路中点垂测数据作为输入，利用粒子群算法优化ITU-RP.533建议最高可用频率计算公式中的20个常数参数以获取最优常数参数，建立一种适用于电离层暴期间短波链路的最高可用频率计算方法。附图说明图1是电离层暴期间短波链路MUF计算误差分布和累计误差分布；图2是本专利技术实施例1所公开的计算方法的流程示意图；图3是使用本专利技术实施例1所公开的计算方法于2017年2月20日至21日正相电离层暴期间短波链路MUF计算结果与实测结果对比；图4是使用本专利技术实施例1所公开的计算方法于2017年3月22日负相电离层暴期间短波链路MUF计算结果与实测结果对比。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图和实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。实施例1，如图2所示，本实施例公开了一种适用于电离层暴期间短波链路的最高可用频率计算方法，包括如下步骤：步骤1、在全国范围内优选5条优质斜测链路作为背景链路；步骤2、利用国内14个电波观测站的垂测数据fof2、foE和M(3000)F2因子，通过克令格插值算法获取5条背景链路中点位置的垂测数据fof2、foE和M(3000)F2因子，上述fof2指电离层F2层临界频率、foE指电离层E层临界频率；步骤3、分别获取历史电离层暴期间5条背景链路中点的3个垂测数据和斜测链路最高可用频率MUF数据，并对背景链路数据进行预处理，作为下述步骤4中的训练样本；所述步骤3具体包括如下步骤：步骤31、将fof2相对偏差df作为电离层暴的判定标准，其中，式中，fof2为观测值，fof2med为月中值，当df不小于0.3时即认为是正相电离层暴，当df不大于-0.2时即为负相电离层暴；步骤32、利用国内14个电波观测站的垂测数据和克令格插值算法，获取某一链路中点3个垂测数据，并根据以上电离层暴的判定标准，分别筛选某一时间段内电离层暴期间链路中点的fof2、foE和M(3000)F2因子3个垂测数据和该斜测链路的最高可用频率数据；步骤33、对于同一时刻收发点互换的两条链路MUF的差值不超过3MHz，即认为该组数据质量可靠，可作为样本数据。本文档来自技高网...

【技术保护点】
一种适用于电离层暴期间短波链路的最高可用频率计算方法，其特征在于，包括如下步骤：步骤1、在全国范围内优选5条优质斜测链路作为背景链路；步骤2、利用国内14个电波观测站的垂测数据fof2、foE和M(3000)F2因子，通过克令格插值算法获取5条背景链路中点位置的垂测数据fof2、foE和M(3000)F2因子，上述fof2指电离层F2层临界频率、foE指电离层E层临界频率；步骤3、分别获取历史电离层暴期间5条背景链路中点的3个垂测数据和斜测链路最高可用频率MUF数据，并对背景链路数据进行预处理，作为下述步骤4中的训练样本；步骤4、将上述步骤3中插值得到的电离层暴期间链路中点垂测数据作为ITU‑R P533建议短波最高可用频率计算公式中的输入，将斜测数据作为输出，利用粒子群算法优化公式中的20个常数参数，设定其优化次数或优化误差，以获得最优常数参数；步骤5、将获得的20个最优常数参数替换原公式中的常数参数，以国内14个电波观测站的垂测数据插值获得指定链路中点的垂测数据作为输入，建立一种适用于电离层暴期间短波链路的最高可用频率计算方法。

【技术特征摘要】
1.一种适用于电离层暴期间短波链路的最高可用频率计算方法，其特征在于，包括如下步骤：步骤1、在全国范围内优选5条优质斜测链路作为背景链路；步骤2、利用国内14个电波观测站的垂测数据fof2、foE和M(3000)F2因子，通过克令格插值算法获取5条背景链路中点位置的垂测数据fof2、foE和M(3000)F2因子，上述fof2指电离层F2层临界频率、foE指电离层E层临界频率；步骤3、分别获取历史电离层暴期间5条背景链路中点的3个垂测数据和斜测链路最高可用频率MUF数据，并对背景链路数据进行预处理，作为下述步骤4中的训练样本；步骤4、将上述步骤3中插值得到的电离层暴期间链路中点垂测数据作为ITU-RP533建议短波最高可用频率计算公式中的输入，将斜测数据作为输出，利用粒子群算法优化公式中的20个常数参数，设定其优化次数或优化误差，以获得最优常数参数；步骤5、将获得的20个最优常数参数替换原公式中的常数参数，以国内14个电波观测站的垂测数据插值获得指定链路中点的垂测数据作为输入，建立一种适用于电离层暴期间短波链路的最高可用频率计算方法。2.根据权利要求1所述的适用于电离层暴期间短波链路的最高可用频率计算方法，其特征在于，所述步骤3具体包括如下步骤：步骤31、将fof2相对偏差df作为电离层暴的判定标准，其中，式中，fof2为观测值，fof2med为月中值，当df不小于0.3时即认为是正相电离层暴，当df不大于-0.2时即为负相电离层暴；步骤32、利用国内14个电波观测站的垂测数据和克令格插值算法，获取某一链路中点3个垂测数据，并根据以上电离层暴的判定标准，分别筛选某一时间段内电离层暴期间链路中点的fof2、foE和M(3000)F2因子3个垂测数据和该斜测链路的最高可用频率数据；步骤...

【专利技术属性】
技术研发人员：王飞飞，盛冬生，孙树计，刘玉梅，
申请(专利权)人：中国电波传播研究所中国电子科技集团公司第二十二研究所，
类型：发明
国别省市：山东,37