线性调频(linear frequency modulation,LFM)信号作为一种雷达,声呐等系统中常用的信号形式,对它进行检测和方位估计具有重要的意义。然而,奈奎斯特采样定理指出,在采用阵列对信号进行方位估计时,要求相邻阵元间距小于信号波长一半,高中心频率使LFM信号很难满足这个条件,而空间欠采样将会带来阵元间相位差的模糊值并因此难以估计出信号的方位。本文主要讨论解决多个LFM信号方位估计问题的方法。
大量文献对LFM信号方位估计问题进行了研究。其中,文献[1]提出了基于分数阶Fourier变换的单个LFM信号DOA估计方法,但没有考虑到多个LFM信号方位估计问题;文献[2]研究了基于欠采样下的单个超宽带LFM信号调制斜率和初始频率的估计方法,但是没有涉及到欠采样下多个LFM信号方位估计问题。文献[3, 4]从数学角度深入研究了欠采样问题。文献[3]提出了重构单个实数的广义稳健中国剩余定理,并利用该定理解决了单源窄带信号的方位估计问题。文献[4]提出了同时估计多个实数的广义稳健中国剩余定理,并利用它解决了时间欠采样下多源窄带信号频率估计的问题。
受到文献[3, 4]的启发,本文提出了一种空间欠采样下多个LFM信号方位估计的新方法,该方法以同时估计多个实数的广义稳健中国剩余定理(generalized robust Chinese remainder theorem,GRCRT)为基础,当被估计的多个LFM信号参数满足一定条件时,可以在较低信噪比下稳健的估计出信号的方位角。 1 空间欠采样下多个LFM信号方位估计的问题分析图 1所示的是由L个阵元组成的均匀线列阵,P个来自远场的宽带信号的方位角的集合为θ={θ1,θ2,…,θP},1≤p≤P,θ1<θ2<…<θP。d为相邻阵元间距,则第l个阵元接收到的信号表示为
式中,Ap为信号幅度,为了讨论方便假设信号的幅度相同且幅值均为A,fp,kp分别为第p个信号的初始频率和调制斜率(已知),τl(θp)(1≤l≤L,1≤p≤P)为第l个阵元接收到第p个信号的延迟时间,τ1=0,w(t)是与信号无关的均值为0,方差为1的复高斯白噪声。
用采样频率fs对信号xl(t)进行采样,可以得到
1.1 LFM信号方位的估计受到文献[2]中多个线性调频信号调制斜率估计方法的启发,本文希望把所有阵元上的信号转化为单频信号后再进行方位估计。方法如下:把每个阵元上采样后的信号xl(n)延时τ后,取τ=0.4T[2],与xl(n)共轭相乘,可以得到:
式中,N=τ*fs。从上面公式可以看出,右边第1项是多个频率为kpτ的正弦信号的叠加,需要估计它的信号的相位谱中包含的延迟信息nl进而估计出信号的方位;第2项是各个信号的交叉项;第3项和第4项是信号与噪声的交叉项;第5项是噪声的自相关项。假设信号与噪声是不相关的,而且高斯白噪声本身是相互独立的,因此经过把延迟后的信号与原信号共轭相乘基本达到了把非平稳的线性调频信号转化为平稳的单频信号的目的,但目前(4)式中右边第1项相位的估计还会受到信号交叉项的影响。
1.2节将分析忽略交叉项影响的条件在交叉项影响可以忽略的情况下,得到信号r1(n)和rl(n)可以近似认为是多个频率为kpτ的正弦信号的叠加,对它们做FFT得到每个阵元上信号频谱,再对各自频谱进行谱峰搜索,找到最大的P个值并求出这些P个值对应的相位,从(4)式中可以得到第1个阵元上的信号r1(n)的P个谱峰对应的相位分别是,第l个阵元上的信号rl(n)的P个谱峰对应的相位分别是
在得到相位估计值后,为了估计出P个信号的入射方位角,需要得到阵元间转化后的同频率信号的相位差,这就要考虑到信号频率和相位的对应问题。相同频率谱峰相位之差即为P个信号的相位差估计值,观察(4)式可以看出,在非空间欠采样的情况下,只要相同频率谱峰相位相减即可得到第l个阵元上来自第p个方位的信号与第1个阵元上来自第p个方位的信号的相位差准确值,又因为kp和N均是已知的,就可以得到nl,继而因为nl=dlsinθpfs/c,dl为阵元间距,可以直接求出第p个方位角,但在空间欠采样时得到的是相位差的模糊值,所以需要利用广义稳健剩余定理对它们进行解模糊进而得到信号的真实方位角。本文第3节将具体分析空间欠采样下LFM信号方位估计问题。
1.3 信号交叉项对方位估计的影响分析通过前面的分析可以看出,估计多个LFM信号的方位受到信号交叉项的影响,如果信号交叉项可以忽略,则可以得到每个阵元上的相位之后求得阵元间的相位差,进而估计出多个LFM信号的方位。现在讨论如何忽略交叉项的影响。
对于(4)式,包含每个阵元上信号相位信息的复正弦矩形脉冲信号的幅度谱的最大值可以表示为:
交叉项中任何一个线性调频信号的幅度谱的最大值可以表示为:
从而可以得到。因此只要保证即可忽略信号交叉项的影响。
2 同时估计多个实数的广义稳健中国剩余定理描述设有k个正实数集合Φ={N1,N2,…,NK},N1<N2<…<NK和L个模数M1,M2,…,ML并且Ml满足:
式中,M=gcd(M1,M2,…,ML)是L个模数M1,M2,…,ML的最大公约数,Γi和Γj互质,Γ1<Γ2<…<ΓL,1≤i≠j≤L。在实际应用中也可以选择互质的一组数的整数倍作为模数。
对于1≤l≤L,定义实数集合Φ模Ml的余数集合为:
但在许多应用场合,例如数字信号处理和雷达系统中,获得的rkl是有误差的,即
当已知有误差的余数和模数Ml时,则需要估计的实数Nk可以表示为:
式中,nkl表示模Ml的商,0≤≤Ml-1,|-rkl|=|εkl|≤εM,εM表示的是余数的误差上限。
对得到的有误差的余数从小到大进行排序,即使
式中,
定义,则可以得到。
定义,则可以得到,即可以得到一个方程组:
根据上述定义,同时估计多个实数的广义稳健中国剩余定理可以表述为:
定理1设,则Nk模Ml,1≤k≤K,1≤l≤L后的余数可以由余数集合确定。确定的原则如下:
定理2 假设定理1的条件成立,模数满足(7)式条件,Γi和Γj互质,1≤i≠j≤L。如果
则求解(12)式中的方程组可以得到唯一一组解nkl。
获得了和(10)式中的nkl后,整数Nk的估计值为:
估计的误差上限是:
定理1和定理2的证明参考文献[3]。
3 空间欠采样下基于同时估计多个实数GRCRT的多个LFM信号方位估计定义dij表示的是第i和第j个阵元的间距,即dij=|i-j|d,d为相邻阵元间距,当相邻阵元间距大于所有入射信号的半波长时(对于线性调频信号,信号波长定义为,fp表示线性调频信号的中心频率),入射信号的方位角存在模糊区间,对于第p个信号,其模糊区间Ψp可以表示为
当入射信号的方位角θp处在模糊区间内时,间距为dij的2个阵元接收到第p个信号的相位差π,因而在把线性调频信号转化为频率为kpτ的单频信号后不同阵元上同频率单频信号之间的相位差Φ′p,dij也是模糊的,并且是以2π为模糊的观测值,即
式中,为转化后的单频信号的波长。
通过FFT运算得到相位差的估计值的范围是为未知的模糊数,为估计的误差。(17)式转化为目标方位角形式为
式中,d0为一个固定值。
将GRCRT应用到空间欠采样下的信号方位估计,其中的参数对应关系为:,分别对应GRCRT中的。dij和d0的选取应满足GRCRT定理中对模数Ml的要求。
通过以上分析,P个线性调频信号方位角估计的步骤如下:
(1)从均匀线列阵中选择S组阵元,每组包含2个阵元,每组阵元的选择需满足GRCRT定理的要求;
(2)对其中1组阵元上的信号进行1次采样,采样后信号延迟与原信号共轭相乘把信号近似转化为多个频率为kpτ的单频信号的叠加,采用FFT算法和谱峰搜索,得到每个阵元上信号的频谱,相同频率的谱峰对应的相位之差即为对应的P个信号的1组相位差估计值;
(3)对其他组阵元上信号重复第2步直至获得P个信号的S个相位差,并且建立转化后的多个单频信号频率和相位差模糊估计值的对应关系;
(4)利用GRCRT定理分别估计出P个信号的无模糊方位角。
4 仿真实验及结果设有2个幅度相同且都来自远场直线传播的LFM信号,其初始频率,调制斜率和入射方位角(fp,kp,θp) (1≤p≤2)分别为:550MHz,5.5×107 MHz/s,60°),(600 MHz,6×107 MHz/s,65°),信号持续时间为T=10 μs,信号波长用中心频率近似计算可得:λ1=0.36 m,λ2=0.33 m(传播速度为c=3×108 m/s)并且满足大于1的条件;相邻阵元间距d=0.5 m,由于d>λk/2(k=1,2),根据(16)式,如果信号1和信号2的方位角分别在区间内时,阵元间对应的相位差是模糊的,可见仿真选取的信号方位角分别位于各自的方位角模糊区间内。
根据上述信号参数及GRCRT,在空间欠采样下分析两组阵元上的信号,第1组是1号阵元和6号阵元,第2组是1号阵元和5号阵元,即d16=5d,d15=4d。设定常数d0=80d,则有d0/d16=16,d0/d16=20(对应GRCRT定理中M为4,Γ1和Γ2为4和5)。
为了验证在空间欠采样下所提方法进行多个宽带LFM信号方位估计的性能,仿真实验给出了信噪比变化时方位角的估计性能曲线。图 2a为估计的信号方位角与信噪比的关系曲线;图 2b为方位角估计的均方根误差与信噪比之间的关系曲线。图 2a和图 2b均是10 000次Monte-Carlo实验的结果。图中结果表明,在信噪比大于-5 dB并且参数的选择满足定理要求时,能利用同时估计多个实数的GRCRT定理稳健地重构多个宽带LFM信号的方位角。
5 结 论本文提出了一种空间欠采样下多个LFM信号方位无模糊估计的方法。当被估计的LFM信号参数满足一定条件时,可以忽略信号交叉项对信号方位估计的影响,这时利用同时估计多个实数的GRCRT可以估计出无模糊的方位。仿真实验验证了所提出方法的有效性和稳健性,新提出的方法适用于空间欠采样下多个LFM信号的方位估计。
[1] | 张艳红,齐林,穆晓敏,陶然. 基于分数阶傅里叶变换的WLFM信号DOA估计[J]. 信号处理,2005,21(4A):57-60 Zhang Yanhong, Qi Lin, Mu Xiaomin, Tao Ran. DOA Estimation of WLFM Signal Based on Fractional Fourier Transform[J]. Signal Processing, 2005, 21(4A):57-60 (in Chinese) |
Cited By in Cnki (6) | |
[2] | 沈显祥,叶瑞青,唐斌. 基于欠采样的宽带LFM信号参数估计[J]. 电波科学学报,2007,22(1):43-46 Shen Xianxiang, Ye Ruiqing, Tang Bin. An Algorithm for Estimation of Wideband LFM Signal Parameters Based on Sub-Sampling[J]. Chinese Journal of Radio Science, 2007, 22(1):43-46 (in Chinese) |
Cited By in Cnki (2) | Click to display the text | |
[3] | Li Xiaowei, Liang Hong, Xia Xianggen. A Robust Chinese Remainder Theorem with Its Applications in Frequency Estimation from Undersampled Waveforms[J]. IEEE Trans on Signal Processing, 2009, 57(11):4314-4322 |
Click to display the text | |
[4] | Liang Hong, Zhang Heng, Jia Ning. A Generalized Robust Chinese Remainder Theorem for Multiple Numbers and Its Application in Multiple Frequency Estimation with Low Sampling Rates[C]//Proc Int Conf Signal Processing, Communication and Computing,2011 |