主权项 |
一种双线性相参的立方相位信号参数估计方法,其特征在于,所述方法包括如下步骤:步骤1,采集一段含有噪声的信号,用离散多分量立方相位信号z(n)表示,并计算所述离散多分量立方相位信号z(n)的总能量;n表示离散采样时间变量,<img file="FDA0001053400710000011.GIF" wi="286" he="103" />N为离散信号长度;所述离散多分量立方相位信号z(n)包含P个信号分量;步骤2,计算所述离散多分量立方相位信号z(n)的参数化双线性瞬时自相关函数R<sub>z</sub>(n,m);m表示离散延时变量;步骤3,构造参数化双线性瞬时自相关函数R<sub>z</sub>(n,m)的时间‑调频率分布函数G<sub>z</sub>(n,f<sub>m</sub>),并对所述时间‑调频率分布函数G<sub>z</sub>(n,f<sub>m</sub>)进行复数求模操作,并得到模值的整数次幂函数Κ<sub>z</sub>(n,f<sub>m</sub>);f<sub>m</sub>表示离散频域变量;步骤4,对所述模值的整数次幂函数Κ<sub>z</sub>(n,f<sub>m</sub>)沿f<sub>m</sub>轴做逆傅里叶变换,得到逆傅里叶变换后的函数Q<sub>z</sub>(n,b),再对逆傅里叶变换后的函数Q<sub>z</sub>(n,b)采用变尺度傅里叶变换完成沿离散采样时间轴的能量相参积累,得到能量相参积累函数Γ<sub>z</sub>(f<sub>[nb]</sub>,b);b表示离散频域变量f<sub>m</sub>对应的离散时域变量,f<sub>[nb]</sub>表示离散采样时间变量n对应的变尺度频域变量;步骤5,沿能量相参积累函数Γ<sub>z</sub>(f<sub>[nb]</sub>,b)的时间轴b做傅里叶变换,得到参数化调频率‑二次调频率分布函数Ψ<sub>z</sub>(f<sub>[nb]</sub>,f<sub>b</sub>),f<sub>b</sub>表示离散时域变量b对应的频域变量;步骤6,确定所述参数化调频率‑二次调频率分布函数Ψ<sub>z</sub>(f<sub>[nb]</sub>,f<sub>b</sub>)的最大峰值,得到最大峰值对应的第j个信号分量的调频率估计值和二次调频率估计值,j=1,…,P;步骤7,根据第j个信号分量的调频率估计值和二次调频率估计值对离散多分量立方相位信号z(n)进行多普勒扩散补偿,并确定最大峰值对应的第j个信号分量的中心频率估计值和幅度估计值;步骤8,根据第j个信号分量的调频率估计值、二次调频率估计值、中心频率估计值以及幅度估计值重构第j个信号分量;步骤9,从离散多分量立方相位信号z(n)中滤除第j个信号分量,得到剩余离散信号,并计算剩余离散信号的能量;步骤10,若剩余离散信号的能量与离散多分量立方相位信号的总能量的比值大于或者等于预设门限,则将离散多分量立方相位信号替换为剩余离散信号,并依次重复执行步骤2至9,直到剩余离散信号的能量与离散多分量立方相位信号的总能量的比值小于预设门限时,所述立方相位信号的参数估计结束。 |