[发明专利]基于普通克里金插值的电磁频谱地图构建方法

权利要求书

1.一种基于普通克里金插值的电磁频谱地图构建方法，其特征在于：所述基于普通克里金插值的电磁频谱地图构建方法包括以下步骤：

1)由传感器获得的观测值计算实验变差函数得到实验变差函数在离散的h上的值；所述实验变差函数的具体计算方式是：

其中：

所述i表示传感器的序号；

所述z(s_i)表示位置s_i处的观测值；

所述z(s_i+h)表示与位置s_i间距离为h处的观测值；

所述N(h)表示相距h的传感器对的数量；

2)选择理论变差函数模型，对实验变差函数在离散的h上的值进行拟合，得到理论变差函数γ(h)；所述拟合是将函数定义域拓展到任意h＞0，进而计算任意h＞0对应的理论变差函数γ(h)的值；

3)根据步骤2)拟合得到的理论变差函数γ(h)，计算位置s_z处的普通克里金权重和拉格朗日乘子L(s_z)，z取1-M的自然数，M为位置点的总数量；

所述普通克里金权重和拉格朗日乘子L(s_z)的计算方式是：

其中N为传感器数量；

4)根据步骤3)计算得到的位置s_z处的普通克里金权重和拉格朗日乘子L(s_z)计算未布设传感器位置s_z处的观测值，所述未布设传感器位置s_z处的观测值的计算方式是：

所述未布设传感器位置s_z处的观测值在计算过程中的估计方差的计算方式是：

5)取z＝z+1，重复步骤1)至步骤4)；直至估计出区域内任意未布设传感器位置处的观测值，得到观测值在整个区域内的分布状况，完成电磁频谱地图构建。

2.根据权利要求1所述的基于普通克里金插值的电磁频谱地图构建方法，其特征在于：所述步骤2)中理论变差函数模型是球状模型、指数模型或高斯模型。

3.根据权利要求2所述的基于普通克里金插值的电磁频谱地图构建方法，其特征在于：所述步骤2)中理论变差函数模型是球状模型时，所述理论变差函数γ(h)的具体表达式是：

其中：

C₀、C和a为模型中的待估计值。

4.根据权利要求2所述的基于普通克里金插值的电磁频谱地图构建方法，其特征在于：所述步骤2)中理论变差函数模型是指数模型时，所述理论变差函数γ(h)的具体表达式是：

其中：

C₀、C和a为模型中的待估计值。

5.根据权利要求2所述的基于普通克里金插值的电磁频谱地图构建方法，其特征在于：所述步骤2)中理论变差函数模型是高斯模型时，所述理论变差函数γ(h)的具体表达式是：

其中：

C₀、C和a为模型中的待估计值。