[发明专利]无线电环境图场强参数估计算法

权利要求书

1.无线电环境图场强参数估计算法,其特征在于，包括如下步骤

步骤1、获取样本点坐标和样本点处采样值，以及待估点坐标，再由各点坐标求出两点间的距离并筛选出滞后距；

步骤2、根据滞后距、滞后距处坐标点对数量和对应坐标的区域化变量值，定义无线电环境下的变异函数，再将采样值代入变异函数获得实际变异函数值；

所述的步骤2，与传输损耗量纲匹配的无线电环境下的变异函数γ(h)为

n(h)为在滞后距h处坐标点对数量，z(x_i)为区域化变量在坐标x_i处的取值，z(x_i+h)是与坐标x_i距离为h的区域化变量的取值；

步骤3、根据Longley-Rice模型，引入不同距离条件下的对应损耗系数，由损耗系数表征出理论变异函数、考虑自由空间传输损耗的总理论变异函数；所述的步骤3，包括如下步骤：

步骤3.1、根据Longley-Rice模型的传播损耗参考中值，表示出传播损耗函数和考虑自由空间传输损耗的传播损耗函数；

步骤3.1.1、传播损耗函数L为L＝A_el+k₁d+k₂log_e(d/d_LS)，A_el为自由空间下视距，k₁、k₂为传播损耗系数，d_LS为临界距离，d为传播距离；

步骤3.1.2、考虑自由空间传输损耗的传播损耗函数L_f为L_f＝A_ref+32.45+201lgd+201lgf，32.45+20lgd+20lgf为自由空间传播损耗，f为发射频率，d为传播距离，A_ref为传播损耗参考中值；

步骤3.2、分别对应步骤3.1.1和步骤3.1.2的两个传播损耗函数写出理论变异函数和考虑自由空间传输损耗的总理论变异函数h为距离，a₁、a₂、a₃、a₄为待定系数，ε为极小常数，取θ＝14*10^∧3；

步骤4、根据样本点坐标的分布特征和无线电在样本点的传播特征，结合实际变异函数值，通过迭代算法获得理论变异函数和总理论变异函数的函数式，作为理论变异函数模型；

步骤5、由普通克里金方程组，结合理论变异函数模型，求解出待估点坐标的估值。

2.根据权利要求1所述的无线电环境图场强参数估计算法,其特征在于，所述的步骤4，包括如下步骤，

步骤4.1、根据样本点坐标的分布特征和无线电在样本点的传播特征，定义具有权重系数的适应度函数；

步骤4.2、由适应度函数结合实际变异函数值通过粒子群优化算法将理论变异函数和总理论变异函数拟合，获得理论变异函数和总理论变异函数内的待定系数，将理论变异函数式和总理论变异函数式作为理论变异函数模型。

3.根据权利要求2所述的无线电环境图场强参数估计算法,其特征在于，所述的步骤4.1，包括如下步骤

步骤4.1.1、根据样本点的特征，定义第一权重系数λ₁＝N/N_i，N_i代表对应滞后距处的样本点对数、N代表总的样本点对数，

定义第二权重系数代表实际变异函数值的平均值、γ(h_i)代表变异函数在h_i位置的值，

并定义第三权重系数为滞后距的平均值、h_i为第i个滞后距；

步骤4.1.2、由步骤4.1.1得出最终权重系数λ＝λ₁·λ₂·λ₃；

步骤4.1.3、将适应度函数与最终权重系数λ结合，得到

其中，F(j)表示第j个变量的适应度函数值，h_i,j表示第j个变量的第i个滞后距，γ^*(h_i,j)表示h_i,j位置的变异函数估计值，γ(h_i,j)表示在h_i,j位置的实际变异函数值。