[发明专利]一种基于粒子群算法的炮点纠偏方法

权利要求书

1.一种基于粒子群算法的炮点纠偏方法，包括以下步骤：

S10，估算近地表直达波的速度；

S20，根据近地表直达波的速度算出直达波的理论计算初至，通过比较理论计算初至与实际拾取初至的误差来建立炮点偏移校正的误差目标函数；

S30，利用粒子群算法反演出误差目标函数的最优解，从而获得炮点偏移校正量，进而确定炮点纠偏的位置。

2.根据权利要求1所述的炮点纠偏方法，其特征在于，步骤S10中，采用最小二乘法拟合炮间距与实际拾取初至，估算出近地表直达波的速度。

3.根据权利要求2所述的炮点纠偏方法，其特征在于，步骤S10中，通过下式的拟合函数拟合炮间距与实际拾取初至，估算出近地表直达波的速度：

h_i＝vt_i+c

其中，h_i为第i道的炮检距，t_i为第i道的实际拾取初至，v为近地表直达波的速度，i的取值范围为1到N，N为接收直达波的检波器的个数,c为常数。

4.根据权利要求1所述的炮点纠偏方法，其特征在于，步骤S20中，所述炮点偏移校正的误差目标函数F为下式：

F=Σi=1N(ti-τi)2]]>

其中，t_i为第i道的实际拾取初至，τ_i为第i道的理论计算初至，N为接收直达波的检波器的个数。

5.根据权利要求1所述的炮点纠偏方法，其特征在于，步骤S30中，所述误差目标函数的最优解是指满足迭代终止条件时的群体最佳位置。

6.根据权利要求5所述的炮点纠偏方法，其特征在于，所述步骤S30包括以下步骤：

S30.1初始化粒子群算法参数；

S30.2计算每个粒子的误差目标函数值；

S30.3对每个粒子的误差目标函数值进行判断，选取每个粒子在自身计算历史中函数值最小的位置为个体最佳位置，选取所有粒子中函数值最小的粒子位置为群体最佳位置；

S30.4根据获得的个体最佳位置和群体最佳位置，通过速度与位置的更新算式更新每个粒子的速度与位置；

S30.5判断迭代次数或所述群体最佳位置是否满足迭代终止条件；

若满足，结束迭代，所述群体最佳位置为炮点实际位置，根据炮点初始位置和炮点实际位置计算出炮点偏移校正量；

否则，返回步骤S30.2。

7.根据权利要求6所述的炮点纠偏方法，其特征在于，步骤S30.1中，所述粒子群算法参数包括：炮点初始位置、炮点搜索范围、粒子群种群数、粒子位置、粒子速度、认知速度权重和社会速度权重。

8.根据权利要求7所述的炮点纠偏方法，其特征在于，所述炮点搜索范围以炮点初始位置为中心；所述粒子速度由惯性速度、认知速度和社会速度向量生成。