[发明专利]基于自适应STFT方法的微多普勒信号参数估计方法

权利要求书

1.基于自适应STFT方法的微多普勒信号参数估计方法，其特征在于包括以下步骤：

S1、无人机旋翼回波的数字信号x(n)表示为：

式①中：n为采样点数，n＝1，2，……N；T_s为采样时间间隔，单位为秒；P为无人机旋翼个数，p＝1，2……P；L为旋翼长度，单位为米；目标与雷达视线夹角为β弧度，f_r为旋翼频率，φ_i为第i个旋翼的初始相位，λ为电磁波波长；ξ为反射信号幅度，单位为伏特；

S2、对式①进行傅里叶变换：

X(w)＝FFT[x(n)]；     ②

S3、对式②中相邻谱线的频率间隔进行标记，并计算对应的频率大小f_AM：

式③中：n_a，n_b分别表示相邻谱线对应的频谱点数；

S4、利用式③获取多普勒频谱范围，即获得目标旋翼引起的最大多普勒频率f_d；

S5、对式①进行STFT：

S6、对式④进行SST：

式⑤中T表示SST变换，w(n,w)表示瞬时频率：

式⑥中diff代表对函数S(n,w)关于n进行差分运算；

S7、抽取T[S(n,w)]时间维上多个点，估计无人机旋翼个数即：

S8、将式⑦获得结果代入中，重复步骤S5～S7，进一步验证选取的窗长以及获得的旋翼数，并求出旋翼频率f_r：

S9、通过f_r以及f_d获得无人机旋翼长度L为：

。

2.根据权利要求1所述的基于自适应STFT方法的微多普勒信号参数估计方法，其特征在于：在所述步骤S5中，

x(m)的长度D通过f_AM的大小进行选择：

窗函数g的窗长为：

3.根据权利要求1或2所述的基于自适应STFT方法的微多普勒信号参数估计方法，其特征在于：在所述步骤S5中，

所述的x(m)的长度D初始值为

所述的窗函数g的窗长初始值为

4.根据权利要求1所述的基于自适应STFT方法的微多普勒信号参数估计方法，其特征在于：在所述步骤S7中，Q为抽取的总点数，q_i表示抽取的第i个点中幅度大于阈值后剩余的点数总和，q_i表示为：

式中η为阈值系数，η的取值范围为0.6～1。