四阶统计量和时差定位中的应用
- 格式:pdf
- 大小:125.34 KB
- 文档页数:5
雷达与对抗 2002年 第2期 27
四阶统计量在时差定位中的应用
刘东霞,林象平
(西安电子科技大学电子对抗研究所,陕西西安710071)
摘 要:实际时差定位中需要从有高斯噪声的环境中估计信号的时延,但基于互相 关或互功率谱的时延估计方法对噪声十分敏感,本文作了基于四阶统计量的时延 估计,并给出了几种求解延时的方法,可以很好地解决这一问题。
关键词:信号处理;时延估计;互相关
中图分类号:TN911.6文献标识码:A文章编号:1009—0401(2002)02—0027—05
Application of Fourth—Order Statistics in
the Time Delay Estimation
LIU Dong—xia,LIN Xiang—ping
(Research institute ofElectronic Counterm ̄-ures,XMian University,Xi’an 711XI71,Ckina)
Abstract:It is required tO estimate the time delay of signals in the Gaussian noise surroundings in the practical location.The time delay estimation(TDE)method based on cross—correlation and cross—power spectra is sensitive tO noise.The TDE method based on fourth—order statistics is dis—
cussed in this paper and several approaches tO find the time delay are given,which can efficiently solve the problem that traditional methods cannot do. Key words:signal processing;time delay estimation;cross—correlation
1 引 言
高阶统计量分析是近几年国内外信号处理领域内的一个前沿课题。高阶统计量广泛应用
于所有需要考虑非高斯性、非最小相位、有色噪声、非线性或循环平稳性的各类问题中。我们 知道,功率谱是极为常用的信号处理工具。但是,功率谱只反映了过程的幅度信息,即只能完
整地描述一个已知均值过程,有明显的局限性,这种局限性也限制了它的应用。与功率谱相 比,高阶谱揭示了更多有关过程的统计信息,有更多的独特性质,如能自动抑制高斯噪声,保留
了信号的时延特性等,可以解决许多功率谱等信号处理技术所无法解决的问题。所以,高阶统 计量特别是高阶累积量及其谱是一种更为有效的信号处理工具。将该技术用于电子对抗领
域,充分利用高阶谱的丰富信息和优异性质,可以解决电子对抗所面临的技术难题。下面介绍
收稿日期:2002 02.28 作者简介:刘东霞(1978一),女,河南漯河人,西安电子科技大学电子工程学院硕士在读,目前主要从事信 号处理与仿真研究。
维普资讯 http://www.cqvip.com 28 雷达与对抗 2002年 第2期
与本文应用有关的高阶统计量特别是四阶统计量的优异性质。
(1)零均值的实平稳高斯随机过程{z(n)},其高阶累积量及其谱为零,而功率谱不为零。
c, (r)=E{z(71)z(7z+r);:r(r) (1)
h(rl,r2,…, 一1)三0, 是≥3 (2)
S2 (W)≠0 (3)
Sh(Wl,W2,…,W 一1);0, 是≥0 (4)
从该性质可以看出,高阶谱对高斯过程不敏感,而功率谱对高斯过程是敏感的。因此可以 用高阶谱来自动抑制各种高斯噪声。这一性质在雷达信号处理及电子对抗中是极为有用的。
(2)如果随机变量{z }与随机变量{Y i=1,…, 独立,则有 cure(zl+Yl,…,z +Y )=cum(zl,…,z )+cum(Yl,…,Y ) (5)
这一性质给出了累积量的另一名称“半不变量”。由上述两个性质就可以得到一个重要的
结论:如果一非高斯信号是在与之独立的加性高斯色噪声中被观测的话,观测过程的高阶累积 量将与非高斯信号过程的高阶累积量恒等。因此,当使用高阶累积量作分析工具时,理论上可
完全抑制高斯色噪声的影响。 (3)累积量谱抑制线性相移,互累积量可保留线性相移。
(4)对于概率密度函数(pdf)是对称分布的零均值平稳随机序列{z(n)},它的三阶谱恒 为零,而其四阶谱不为零,即 S3 (Wl,W2);0 (6)
S4 (Wl,W2,W3)≠0 (7)
由此可见,四阶谱可以用来处理具有对称分布的pdf信号,而三阶谱却不能。在通信、雷
达及电子对抗系统中的信号可能会有pdf对称分布的特点,并且,对于某些信号而言,三阶谱 的值可能非常小,而四阶谱的值则要大得多,用四阶谱处理这类信号更为有效。所以,四阶谱 在电子对抗中占有重要的地位。
2 四阶统计量在时差定位中的应用
关于高阶统计量在时差定位中的应用,有些书中曾有过此方面的介绍,但它们是基于双谱 的时延估计,不能处理具有pdf对称分布的信号。而四阶累积量谱有处理pdf对称分布信号
的优点。所以,本文将四阶统计量应用于时差定位中,并给出了估算时差的几种方法。 时差定位是利用平面或空间中的多个侦察站,测量出同一个信号到达各侦察站的时间差, 由此确定辐射源在平面或空间中的位置。其定位精度依赖于合理的布站方式和测量到达时差
的精度。可见到达时差的测量是很关键的技术。
考虑一时差定位系统,它由A、B和C三个侦察站组成,可以是三个固定位置的地面侦察
站,也可以是相互位置确知的三个机载侦察站。被定位的雷达可以是和侦察站在同一平面的
雷达,也可以是低仰角雷达。侦察站均采用无方向性天线,覆盖宽的瞬时视野。三个站同时接 收一部雷达辐射的信号。假设在同一平面上,三站和辐射源E的位置如图1所示。
以A站为参考站,三个站接收到的信号分别为 维普资讯 http://www.cqvip.com 刘东霞等 四阶统计量在时差定位中的应用 29
A( )=S( )+硼A( )
.27B(P/)= ( —DI)+WB( )
c( ):S(n—D2)+叫c(n)
其中,{S( )}是零均值的平稳非高斯信号,D为待估计的
时差,WA( )、WB( )、Wc(,z)分别为三个侦察站接收到的 噪声,它们是零均值的高斯过程,可以是不相关、部分相关
或者是相关,并与 (n)独立,这就是四阶累积量的独到之 处。现在利用观测数值.27A(P/)、.27B(P/)、.27c(P/)来估计时延
D1、D2。 A
B 参考方向 C
图1 平面上的时差定位示意图
解决时延估计的基本思想是将 A( )、 B(n)、 c( )作相对移动,最佳匹配将发生在移
位量为D时。考虑到高斯噪声的四阶统计量为零,则有接收信号的四阶累积量:
c4.^ ^( 2 l,DI 2,m 3):C4,ssss( 21, 2,DI 3)+c4.u几眦删(;'711, 2,r ̄i 3)
C4,ssss( l, 2,m 3)
C4,ABAA( f1,DI 2,m 3):C4,ssss( 21一D1,17"1 2, 3)+c4.tc^ u( 21, 2, 3)
C4, ( 1一D1,DI 2, 3)
C4,AcAA( 2l, 2,m 3)=C4,ssss( 2l—D2,17.1 2, 3)+c4.u几眦删( 2 l,DI 2,” 3)
:c4.砧 ( 1一D2,DI 2, 3)
其对应的累积量谱为
同理 S4.q (WI,W2,ZU3)=
exp[一 (硼1 21+叫2 22+叫3 2 3)]
=S4㈣(叫1,W2,叫3)
 ̄4,ABAA(W ,Wz,硼,)=∑∑∑C4,ABAA( , :,m 3) l 一 2… 3… exp[一 (W1 1+W2 2+W3 3)]
exp[一J.(7_U1 21+叫2 2 2+叫3 3)]
exp[一 (硼1( +D1)+7.02 22+叫3 3)]
:S4. (叫1,叫2,w3)exp(一jwiD1) M A A f ∞ ,~∑一:∞ 一∑一。 研
优 娜 ∞ 一∑一,:∞ ~∑一。 研= 一~
∑维普资讯 http://www.cqvip.com 30 雷达与对抗 2002年 第2期
exp[一 (叫l l+叫2 2+叫3 £3)]
:S4. (ZU1,ZU2,w3)exp(一jwlD2)
基于四阶统计量求时延: 方法1 令
Gl(叫1)= —f f (27c) j— S4,A8AA(ZU1,ZU2,W3)
一 S4,^^^^(ZU1,ZU2,W3)
Gz c = 肛 S,AcM(7d)7d)w3S4 (ZU1,ZU2,ZO3) dw2dw3:exp[一 (叫lD1)]
dw2dw3:exp[一 (叫lD2)]
取Gl(叫1)、G2(叫1)的傅里叶逆变换,有
gl(n)=2*3(n—D1) .
g2(n)=2.8(n—D2) 于是Dl就是gl(n)取最大值时的n值,D2就是g2(n)取最大值时的n值。 方法2 只用一个互三谱同时估计出两个时延参数D。和D2。因为有
s ,A鼢(叫。,叫 ,叫,)=∑∑∑c .舳 ( 。, , ,) I… 2… 3… exp[一 (叫1 1+叫2优2+叫3m3)]
:S4. (ZU1,ZU2,w3)exp[一 (ZU1D2+w2D2)]
令
G(叫l,叫2)= S4,ABCA(7d)1,7d)2,ZO3) S4,^^^^(ZU1,ZU2,ZO3) =exp[一J(叫lD2+叫2D2)] (8)
则它的二维傅里叶逆变换为 r r g(nl,n2)=I I G(叫l,叫2)exp Ej(叫l nl+叫2 n2)]dwldw2 一一一 =(27c) ( l—Dl, 2一D2)
该函数在Dl、D2处取峰值。 方法3
一维傅里叶逆变换法
g。(叫,n):l Jl: G( ,叫)exp ( n) I:27c (n g (叫,n):l-r: G(叫, )exp ( n) l=27c (n D1)
D2)
其中G( ,叫)由式(8)定义。于是,gl(叫,n)和g2(叫,n)分别在Dl、D2处取峰值。
方法4 最大发生频度法 令Dl 是由gl(叫,n)在频率叫 ∈[一兀,兀]取峰值时的n值,假定对几个频率叫 ,i=1,2,
维普资讯 http://www.cqvip.com