[发明专利]基于Radon变换与拉普拉斯算子的微多普勒抑制方法

权利要求书

1.基于Radon变换与拉普拉斯算子的微多普勒抑制方法，其特征在于包括以下步骤：

步骤一：对原始信号进行短时傅里叶变换；

步骤二：对短时傅里叶变换的时频图进行radon变换，并将变换域中的极值点标记为可能的刚体成分；

步骤三：使用拉普拉斯算子对图像的radon变换域进行处理，得到变换域的陡峭程度信息；

步骤四：确定radon变换域中产生可能的刚体成分的峰值点；

步骤五：根据步骤四中刚体成分产生的峰值点在radon变换域中的坐标位置确定刚体在短时傅里叶时频图中的位置参数；

步骤六：根据步骤五得到的刚体成分参数设计带通滤波器，对原信号进行滤波，滤除微多普勒信号成分，得到纯净的刚体信号成分；

所述步骤一中短时傅里叶变换公式如下：

其中，窗函数w(i)长度为Mw，s(i)是原信号，M是信号长度，w(i-m)是窗函数；

所述步骤二中radon变换的公式如下：

其中，f(x,y)是信号短时傅里叶变换后的时频谱图，x是时间，y是频率，δ是冲激函数，α是积分投影的角度，β是投影轴上的位置。

2.根据权利要求1所述的基于Radon变换与拉普拉斯算子的微多普勒抑制方法，其特征在于所述Mw取信号长度M的八分之一。

3.根据权利要求1所述的基于Radon变换与拉普拉斯算子的微多普勒抑制方法，其特征在于所述步骤三中拉普拉斯算子表达式如下

4.根据权利要求1所述的基于Radon变换与拉普拉斯算子的微多普勒抑制方法，其特征在于所述步骤四中确定radon变换域中产生可能的刚体成分的峰值点的方法为：

步骤四一：对变换域的值进行排序，设置经验百分比为1％，取变换域中的函数值为前1％的点；

步骤四二：点函数值大于等于其邻域的函数值即为radon变换域中产生刚体成分的峰值点。

5.根据权利要求1所述的基于Radon变换与拉普拉斯算子的微多普勒抑制方法，其特征在于所述步骤四中确定radon变换域中产生刚体成分的峰值点通过判断步骤二中可能的刚体峰值点位置经过步骤三拉普拉斯算子处理后在该点还是否为峰值确定，若处理后在该点还是峰值，则此位置为刚体成分产生的峰值点。

6.根据权利要求1所述的基于Radon变换与拉普拉斯算子的微多普勒抑制方法，其特征在于所述原始信号为雷达回波中所有含有微多普勒成分的距离单元信号。