一种四阶多项式相位信号参数估计方法技术

本发明专利技术提出了一种四阶多项式相位信号参数估计方法，步骤如下：根据采集的四阶多项式相位信号计算瞬时自相关函数，对瞬时自相关函数傅里叶变换得到二维函数；根据二维函数的峰值的位置坐标分别估计三阶项系数和四阶项系数，构造信号的高阶相位补偿函数，利用高阶相位补偿函数对四阶多项式相位信号进行相位补偿得到高阶相位补偿信号；对高阶相位补偿信号做SoWVD变换，估计二阶多项式相位系数，构造二次项补偿函数实现二阶相位补偿得到二次项补偿信号；对二次项补偿信号做二维傅里叶变换，采用峰值检测算法估计一阶多项式系数。本发明专利技术通过瞬时自相关函数实现参数的降阶，可以采用FFT快速实现，在低计算复杂度条件下，有效提高输出信噪比。输出信噪比。输出信噪比。

【技术实现步骤摘要】
一种四阶多项式相位信号参数估计方法


[0001]本专利技术涉及相位信号参数估计的
，尤其涉及一种四阶多项式相位信号参数估计方法，可以实现高机动雷达目标检测和参数估计。

技术介绍

[0002]雷达技术中，为了实现更高的探测距离和速度分辨率，通常要延长观测时间，由于目标的机动性，回波相位表现为非均匀变化的特点，需要用多项式相位信号模型描述。电子侦察技术中，信号多呈现出时变的非平稳特点，也可以用多项式相位信号模型描述。无线通信技术中，收发双方的相对运动会导致时变的多普勒频率，接收信号呈现出变化的瞬时频率，也可以采用多项式相位信号模型建模。综上，多项式相位信号模型广泛应用于雷达、电子侦察、无线通信等技术中。
[0003]多项式相位信号的各阶参数具有明确的物理含义，例如雷达回波的一阶多项式系数对应目标的径向速度、二阶多项式系数对应目标的加速度、三阶多项式系数对应目标的加加速度等；因此，参数估计是多项式相位信号处理的核心需求之一。多项式的阶数越高，对信号的表达能力越强，能更精确地建模信号，但是阶数的提高会给参数估计带来更高的计算复杂度。正弦信号可以看做一阶多项式相位信号，可以用快速傅里叶变换算法(FFT)予以实现。非平稳信号用一阶多项式相位信号模型建模会出现模型的失配，造成明显的性能损失。雷达、声呐、超声成像中常用的线性调频信号可以看做二阶多项式相位信号，可以采用分数阶傅里叶变换、非均匀快速傅里叶变换(NUFFT)、吕氏分布(LVD)等算法实现参数估计。非平稳性更强的信号用二阶多项式相位信号建模依然会有模型失配现象，造成性能下降。三阶、四阶多项式模型可以建模非平稳性更强的信号。三阶多项式相位信号参数估计方法得到了大量研究，主要有离散调频率傅里叶变换(Discrete Chirp Fourier Transform,DCFT)、多项式Winger
‑
Ville分布(Polynomial Winger
‑
Ville Distribution,PWVD)等方法，但是表达能力更强的四阶多项式相位信号模型研究较少，主要原因是其参数较多，参数空间搜索计算复杂度很高。

技术实现思路

[0004]针对现有高阶多项式相位信号参数估计方法的参数较多，计算复杂度较高的技术问题，本专利技术提出一种四阶多项式相位信号参数估计方法，通过几种瞬时自相关函数实现参数的降阶，使大多数步骤都可以采用FFT实现，实现了计算复杂度和参数估计性能之间的平衡。
[0005]为了达到上述目的，本专利技术的技术方案是这样实现的：一种四阶多项式相位信号参数估计方法，其步骤如下：
[0006]步骤一：根据采集的四阶多项式相位信号计算瞬时自相关函数，对瞬时自相关函数做变尺度傅里叶变换算法和二维傅里叶变换得到二维函数；
[0007]步骤二：根据二维函数的峰值的位置坐标分别估计三阶项系数和四阶项系数，通
过三阶项系数和四阶项系数构造信号的高阶相位补偿函数，利用高阶相位补偿函数对四阶多项式相位信号进行相位补偿，得到高阶相位补偿信号；
[0008]步骤三：对高阶相位补偿信号做SoWVD变换，估计出二阶多项式相位系数，利用二阶多项式相位系数构造二次项补偿函数实现二阶相位补偿，得到二次项补偿信号；
[0009]步骤四：对二次项补偿信号做二维傅里叶变换，采用峰值检测算法估计一阶多项式系数。
[0010]优选地，所述多阶多项式相位信号为四阶多项式相位信号，离散化采样后的四阶多项式相位信号的数学模型为：
[0011][0012]式中，ρ为复幅度，a1为信号频率，a2为调频率，a3和a4分别为二次和三次调频率，m表示慢时间采样序号，且m＝
‑
M/2:M/2
‑
1，M表示脉冲个数，T
r
表示慢时间采样间隔，对应脉冲重复间隔，与脉冲重复频率PRF的关系为T
r
＝1/PRF。
[0013]优选地，计算瞬时自相关函数的方法为：
[0014]引入延时变量m
τ
，通过瞬时对称自相关函数R1(m,m
τ
)降低慢时间维四阶多项式相位信号的阶数；
[0015]利用共轭对称自相关函数消去瞬时对称自相关函数R1(m,m
τ
)的第一个指数项，得到共轭对称自相关函数其中，表示延时常量对应的索引号；
[0016]将瞬时对称自相关函数R1(m,m
τ
)与共轭对称自相关函数相乘消去第一个指数项得到函数x(m,m
τ
)；函数x(m,m
τ
)沿慢时间维计算得到瞬时自相关函数R3(m,m
τ'
,m
τ
)，其中，m
τ'
表示延时变量对应的索引号。
[0017]优选地，所述瞬时对称自相关函数
[0018][0019]其中，ρ1表示瞬时对称自相关函数R1(m,m
τ
)的复幅度，上式简化为：
[0020][0021]所述共轭对称自相关函数
为延时常量，()
*
表示取共轭；
[0022]化简得到共轭对称自相关函数：
[0023][0024]由于共轭对称自相关函数中引入了延时常量第二个指数项中的二维耦合项得到保留，则函数
[0025][0026]式(5)中与m
τ
T
r
存在耦合，为了便于采用FFT算法实现多项式相位系数估计，对函数x(m,m
τ
)沿慢时间维计算，得到瞬时自相关函数：
[0027][0028]优选地，所述变尺度傅里叶变换算法实现变量的尺度变换和变量替换；所述变尺度傅里叶变换算法实现变量的尺度变换的方法为：瞬时自相关函数：
[0029][0030]其中，延时变量ξ表示尺度伸缩因子常量，表示尺度变换之后m
τ'
对应的延时变量；
[0031]沿延时变量m
τ
做累加，消去延时变量m
τ
得函数：
[0032][0033]其中，表示沿延时变量m
τ
做累加；
[0034]变尺度傅里叶变换算法实现变量替换的方法为：函数
[0035][0036]其中，慢时间采样序号η为尺度变换因子常量，表示尺度变换之后m对应的延时变量；
[0037]对函数做二维傅里叶变换，得到二维函数：
[0038][0039]式中，ρ2表示能量峰幅度，分别表示示对应的频率变量采样点序号，N
f
表示频率采样点数，F0表示频率变量采样间隔且f
s
表示采样频率，δ()表示狄拉克函数。
[0040]优选地，估计三阶项系数和四阶项系数的方法为：瞬时自相关函数二维傅里叶变换后得到的二维函数能量峰在平面上的位置为则三阶、四阶多项式相位信号的系数估计值分别为
[0041]根据三阶、四阶多项式相位信号的系数估计值构造信号的高阶相本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种四阶多项式相位信号参数估计方法，其特征在于，其步骤如下：步骤一：根据采集的四阶多项式相位信号计算瞬时自相关函数，对瞬时自相关函数做变尺度傅里叶变换算法和二维傅里叶变换得到二维函数；步骤二：根据二维函数的峰值的位置坐标分别估计三阶项系数和四阶项系数，通过三阶项系数和四阶项系数构造信号的高阶相位补偿函数，利用高阶相位补偿函数对四阶多项式相位信号进行相位补偿，得到高阶相位补偿信号；步骤三：对高阶相位补偿信号做SoWVD变换，估计出二阶多项式相位系数，利用二阶多项式相位系数构造二次项补偿函数实现二阶相位补偿，得到二次项补偿信号；步骤四：对二次项补偿信号做二维傅里叶变换，采用峰值检测算法估计一阶多项式系数。2.根据权利要求1所述的四阶多项式相位信号参数估计方法，其特征在于，所述多阶多项式相位信号为四阶多项式相位信号，离散化采样后的四阶多项式相位信号的数学模型为：式中，ρ为复幅度，a1为信号频率，a2为调频率，a3和a4分别为二次和三次调频率，m表示慢时间采样序号，且m＝
‑
M/2:M/2
‑
1，M表示脉冲个数，T
r
表示慢时间采样间隔，对应脉冲重复间隔，与脉冲重复频率PRF的关系为T
r
＝1/PRF。3.根据权利要求1或2所述的四阶多项式相位信号参数估计方法，其特征在于，计算瞬时自相关函数的方法为：引入延时变量m
τ
，通过瞬时对称自相关函数R1(m,m
τ
)降低慢时间维四阶多项式相位信号的阶数；利用共轭对称自相关函数消去瞬时对称自相关函数R1(m,m
τ
)的第一个指数项，得到共轭对称自相关函数其中，表示延时常量对应的索引号；将瞬时对称自相关函数R1(m,m
τ
)与共轭对称自相关函数相乘消去第一个指数项得到函数x(m,m
τ
)；函数x(m,m
τ
)沿慢时间维计算得到瞬时自相关函数R3(m,m
τ'
,m
τ
)，其中，m
τ'
表示延时变量对应的索引号。4.根据权利要求3所述的四阶多项式相位信号参数估计方法，其特征在于，所述瞬时对称自相关函数其中，ρ1表示瞬时对称自相关函数R1(m,m
τ
)的复幅度，上式简化为：
所述共轭对称自相关函数所述共轭对称自相关函数为延时常量，()
*
表示取共轭；化简得到共轭对称自相关函数：由于共轭对称自相关函数中引入了延时常量第二个指数项中的二维耦合项得到保留，则函数式(5)中与m
τ
T
r
存在耦合，为了便于采用FFT算法实现多项式相位系数估计，对函数x(m,m
τ
)沿慢时间维计算，得到瞬时自相关函数：5.根据权利要求4所述的四阶多项式相位信号参数估计方法，其特征在于，所述变尺度傅里叶变换算法实现变量的尺度变换和变量替换；所述变尺度傅里叶变换算法实现变量的尺度变换的方法为：瞬时自相关函数：其中，延时变量ξ表示尺度伸缩因子常量，表示尺度变换之后m
τ'
对应的延时变量；沿延时变量m
τ
做累加，消去延时变量m
τ
得函数：
其中，表示沿延时变量m
τ
做累加；变尺度傅里叶变换算法实现变量替换的方法为：函数其中，慢时间采样序号η为尺度变换因子常量，表示尺度变换之后m对应的延时变量；对函数做二维傅里叶变换，得到二维函数：式中，ρ2表示能量峰幅度，分别表示分别表示对应的频率变量采样点序号，N
f
表示频率采样点数，F0表示频率变...

【专利技术属性】
技术研发人员：牛志永，庾骏，何红杰，
申请(专利权)人：郑州轻工业大学，
类型：发明
国别省市：