一种基于功率谱的线性调频信号参数估值方法

摘要

本发明公开了一种基于功率谱的线性调频信号参数估值方法，该方法包括：检测待估值的线性调频信号功率谱P(x)；提取所述线性调频信号功率谱P(x)的形状特征CO(P(x))；计算所述形状特征上升沿的斜率β，并根据调频斜率k和所述上升沿斜率β的关系计算所述线性调频信号的调频斜率k；根据所述形状特征上升沿的峰值频率f1与起始频率fs的频率差Δf和所述调频斜率k的关系计算所述频率差Δf；计算所述线性调频信号的起始频率fs和终止频率fe，并计算带宽B；根据脉宽τ和所述带宽B的关系，计算所述脉宽τ。本发明克服现有技术中低信噪比环境下线性调频信号估值结果受噪声影响大而且计算量大的缺陷，实现未知先验信息条件下的线性调频信号参数的估值。

主权项

1.一种基于功率谱的线性调频信号参数估值方法，其特征在于，该方法包括：步骤1、检测待估值的线性调频信号功率谱P(x)；步骤2、提取所述线性调频信号功率谱P(x)的形状特征CO(P(x))；其中，通过以下过程获得所述形状特征：a、选取x轴上的扁平结构元素g₁，长度为L₁，结构元素定义域$D\left[g_1\right]=[x-\frac{L_1}{2},x+\frac{L_1}{2}],$b、首先对P(x)进行形态学膨胀运算得到P₁(x)，表达式如下所示：$P_1\left(x\right)=(P\oplus g_1)\left(x\right)=\max \{P\left(z\right):z\in D[g_1\left]\right\}$其含义为以x点为中心，宽度等于扁平结构元素长度L₁的区域内P(x)幅值的最大值；再对P₁(x)进行形态学腐蚀运算得到P₂(x)，表达式如下所示：P₂(x)＝(P₁Θg₁)(x)＝min{P₁(z):z∈D[g₁]}其含义为以x点为中心，宽度等于扁平结构元素长度L₁的区域内P₁(x)幅值的最小值；最后对P₂(x)重复进行一次腐蚀运算和一次膨胀运算，表达式如下所示：$\mathrm{CO}\left(P\left(x\right)\right)=(\left(P_2\Theta g_1\right)\oplus g_1)\left(x\right);$步骤3、计算所述形状特征上升沿的斜率β，并根据调频斜率k和所述上升沿斜率β的关系计算所述线性调频信号的调频斜率k，其具体包括：a、设定线性调频信号调频斜率k变化区间为[50，2000]，步长为50；b、计算不同k值对应的所述线性调频信号功率谱形状特征，搜索所述线性调频信号功率谱形状特征上升沿峰值对应的频率f₁和半峰值对应的频率f_0.5，根据定义计算出每个所述线性调频信号功率谱形状特征对应的上升沿的斜率β，则获得一一对应的k-β数据点；c、选择指数函数k(β)＝a(β^b)+c对所述k-β数据点进行数据拟合；d、求得系数为：a＝1.4176，b＝-1.959,c＝15.4941；得到调频斜率k与上升沿斜率β的关系如下式：k(β)＝1.4176(β^-1.959)+15.4941根据调频斜率k和所述上升沿斜率β的关系计算所述线性调频信号的调频斜率k。步骤4、根据所述形状特征上升沿的峰值频率f₁与起始频率f_s的频率差Δf和所述调频斜率k的关系计算所述频率差Δf，其具体包括：a、设定线性调频信号调频斜率变化区间为[50，2000]，步长为50；b、计算不同k值对应的线性调频信号功率谱形状特征和每个所述线性调频信号功率谱形状特征对应的上升沿频率差Δf，则获得一一对应的Δf-k数据点；c、选择指数函数Δf(k)＝a(k^b)+c对所述Δf-k数据点进行数据拟合；d、求得系数为：a＝1.007,b＝0.4803，c＝-1.0466；得到调频斜率k与频率差Δf的关系如下式：Δf(k)＝1.007(k^0.4803)-1.0466根据所述形状特征上升沿的峰值频率f₁与起始频率f_s的频率差Δf和调频斜率k的关系计算所述频率差Δf；步骤5、计算所述线性调频信号的起始频率f_s和终止频率f_e，并计算带宽B，其具体包括：a、搜索所述形状特征上升沿和下降沿的峰值处频率f₁和f₁′；b、根据所述起始频率f_s和频率差Δf的关系Δf＝f₁-f_s计算起始频率：f_s＝f₁-Δfc、根据所述终止频率f_e和频率差Δf的关系Δf＝f_e-f₁′计算终止频率：f_e＝f₁′+Δfd、根据公式B＝|f_e-f_s|，则所述线性调频信号的带宽：B＝f_e-f_s；步骤6、根据脉宽τ和所述带宽B的关系，计算所述脉宽τ。