基于空间局部估计的海上短波频率数据重构方法
【专利摘要】本发明涉及一种基于空间局部估计的海上短波频率数据重构方法，属于无线通信频率预测【技术领域】。其特征是利用短波频谱管理设备采集的样本数据，计算各分离距离内样本点对的平均距离和实验变异函数值。利用实验变异函数点对数据对指数模型变异函数进行拟合，得到理论指数变异函数参数。根据短波频率数据重构区域样本点求解空间局部估计方程组，获得各样本点参与数据重构的权重系数。建立待重构点在空间局部的频率值重构公式，求出待重构点的短波通信频率数据。本发明充分考虑了短波频率数据的空间相关性，根据舰船实际使用的数据进行频率数据重构，数据重构精度较准确，能够有效缩短自适应短波通信设备的选频时间，提高海上短波中远程通信的效果。
【专利说明】基于空间局部估计的海上短波频率数据重构方法
【技术领域】
[0001] 本发明属于无线通信频率预测【技术领域】，涉及到一种基于空间局部估计的海上短 波频率数据重构方法。
【背景技术】
[0002] 短波通信是一种主要借助电离层反射天波实现远距离通信的通信方式，是海上舰 船实现中远程通信的重要手段之一。频率的选用制约着短波通信的发展，实现短波通信频 率的优选是发挥短波通信设备性能，提高短波通信效果的关键。
[0003] 目前，海上舰船在短波通信频率优选过程中，主要使用中、长期预测和实时探测两 种手段和方法。中、长期预测方法依据电离层特性参量的时空变化和太阳活动性指数的历 史数据，推断出最高可用频率等参数的中值。缺点是该方法预测的中值与实际通信过程中 的最高可用频率值有较大偏差。频率实时探测方法是对短波频段进行扫描探测，实时处理 收端信号的频率、误码率、信噪比等参数，按照信道质量给频率打分排队，获得可以实时利 用的频率数据。缺点是需要发射无线电探测信号，容易造成电磁环境的恶化，影响通信效 果；探测体制与通信体制不一致时，频率数据的参考价值不高。
【发明内容】
[0004] 本发明要解决的技术问题是运用空间局部重构方法，以海上舰船实际通信频率数 据为重构样本，实现利用已知海域频率数据重构未知海域数据的目的，为海上短波通信频 率选择提供一种可操作性强的新方法。
[0005] 本发明的技术方案是：
[0006] -种基于空间局部估计的海上短波频率数据重构方法，首先以频管设备采集的海 上短波通信频率数据为样本，计算各分离距离内样本点对的平均距离和实验变异函数值； 利用实验变异函数点对数据对指数模型变异函数进行拟合，得到理论指数变异函数参数； 根据短波频率数据重构区域样本点求解空间局部
[0007] 估计方程组，获得各样本点参与数据重构的权重系数。建立待重构点在空间局部 的频率值重构公式，求出待重构点的短波通信频率数据。
[0008] 根据分离距离区间将样本点对分类，计算不同分离距离区间内样本点对的平均距 离，计算落在每段分离距离内的样本点对的实验变异函数值，由各平均距离和实验变异函 数值组成实验变异函数点对。
[0009] 利用得到的实验变异函数点对，对指数模型变异函数进行最小二乘法拟合，删除 严重偏离拟合曲线的特异值点，利用新的实验变异函数点对进行再次拟合，拟合出理论指 数变异函数参数，得到理论指数变异函数。
[0010] 选取参与短波频率数据重构的样本点，计算样本点间及各样本点与待重构点之间 的电离层距离及变异函数值。
[0011] 列出并求解空间局部估计方程组，求出各样本点参与数据重构的权重系数。
[0012] 建立待重构点在空间局部的频率值重构公式，求出待重构点的短波通信频率数 据。
[0013] 本发明充分考虑了短波频率数据的空间相关性，根据舰船实际使用的数据进行频 率数据重构，数据重构精度较准确，能够有效缩短自适应短波通信设备的选频时间，提高海 上短波中远程通信的效果。
【具体实施方式】
[0014] 以下结合技术方案详细叙述本发明的具体实施例。
[0015] 步骤一：利用短波频谱管理设备采集重构区域的样本数据，检查样本数据质量，要 求样本数据要素完整，必须包含的数据要素有：通信收发两端的地理位置信息、通信时间信 息、频率数据信息。
[0016] 步骤二：运用a
【权利要求】
1. 一种基于空间局部估计的海上短波频率数据重构方法，其特征在于如下步骤： 步骤一：利用短波频谱管理设备采集重构区域的要素完整的样本数据，要素完整的样 本数据包含：通信收发两端的地理位置信息、通信时间信息、频率数据信息； 步骤二：运用
汁算步骤一获得的样本数据之间点 对点的电离层距离，个数为N个；其中，（x^y)为收信点A的经、纬度坐标，（Xj，yj)为收信 点B点的经、纬度坐标，SF为尺度因子，根据重构区域的不同纬度范围选取尺度因子值，纬 度在[0°?20° )范围内时，SF = 2. 5,纬度在[20°?40° )范围内时，SF = 2,纬度在 [40° ?90。）范围内时，SF= 1.5; 步骤三：确定步骤二中N个电离层距离的最大值dmax，将分离距离区间h确定为0? 2、2?4、4?6、…、2X ([dmax/2]_l)?2X [dmax/2];按照分离距离区间将N个样本点的 电离层距离分类； 步骤四：计算各分离距离区间h内样本点的电离层距离的平均距离h'，运用
计算样本数据的实验变异函数值，其中Y#(h'）为 实验变异函数，N(h)是分隔距离为h时样本点的个数，Z(Xi，yi)和Z(xi+1，y i+1)是样本点在 空间位置（Xi，yi)，（xi+1，yi+1)的短波通信频率值；由平均距离h'和实验变异函数值Y #(h'） 构成实验变异函数点对（h'，Y#(h'））； 步骤五：依据实验变异函数点对（h'，Y#(h'）），对指数模型变异函数
进行最小二乘法拟合，拟合过程中首先将曲线模型转换成线 性模型，再利用最小二乘法原理对模型参数进行估计，拟合出变异函数及曲线；删除严重偏 离拟合曲线的特异值点，利用新的实验变异函数点对进行再次拟合，得到理论指数变异函 数； 步骤六：选取参与短波频率数据重构的样本点n，计算样本数据点对点之间及各样本 点与待重构点之间的电离层距离，根据电离层距离计算样本数据点对点间及待重构点与各 样本点间的变异函数值； 步骤七：由步骤六中获得的变异函数值构成空间局部估计矩阵表达式：
'其中Y (Xi，χΡ为样本数据点 (Xi, χ」）（i, j = 1，2,. . .，η)之间的变异函数值，γ (Xi, xQ)为样本点xji = 1，2,. . .，η)与待 重构点X(i之间的变异函数值，λ = 1，2,...，η)为重构权重系数，μ为拉格朗日乘子； 步骤八：求解各样本点参与数据重构的权重系数\1和μ ; 步骤九：建立待重构点在空间局部的频率值重构公式：
，求出 待重构点（x〇，y。）处的短波通信频率数据
）为样本点在空间位置 (Xi，yi)的短波通信频率值。
【文档编号】H04B17/00GK104125023SQ201410321642
【公开日】2014年10月29日 申请日期:2014年7月8日 优先权日:2014年7月8日
【发明者】张海勇, 徐池, 韩东, 李梁, 任重, 李利 申请人:中国人民解放军海军大连舰艇学院