空时二维相关非高斯杂波的建模与仿真

王崇辉　邹鲲

摘 要： 在雷达信号处理中，脉冲之间存在时间相关性，不同距离单元的脉冲还存在空间相关性，因此在杂波建模中必须考虑其空时二维相关性。基于球不变随机矢量的杂波建模方法能够独立控制杂波的幅度概率密度函数和杂波序列相关性，但该方法不能同时描述杂波的空时二维相关性，通过对该方法的分析，给出获得空时二维相关的杂波建模方法，并通过计算机仿真，验证了该方法的有效性。

关键词： 杂波； 二维相关性； 球不变随机矢量； 雷达信号

中图分类号： TN958?34； TP391.4 文献标识码： A 文章编号： 1004?373X（2015）15?0025?04

0 引 言

杂波的建模在雷达信号处理领域具有重要的理论和实践意义[1]，这是因为雷达信号处理的目的是在杂波中检测有用信号，而只有对杂波的统计特性充分了解，并建立合适的杂波统计模型，才能设计合适的信号检测算法，并对算法进行性能评估。随着雷达工作频段的提高，在低掠射角（Grazing Angle）的条件下，尤其是在海面环境中，雷达杂波的幅度统计特性明显偏离了瑞利分布，幅度概率密度函数（Probability Density Function，PDF）的尾部增大，表明杂波幅度高的部分概率增加[2]。研究表明，针对瑞利分布的常规信号检测器在非高斯分布噪声下的检测性能会明显下降。因此针对非高斯分布背景下的信号检测算法是目前雷达信号处理领域研究的一个重点。

大量的实测数据和杂波产生的物理机理的研究表明，某个距离单元上的雷达杂波的基带信号可以表示为复合高斯过程[3]或球不变随机矢量（Spherically Invariant Random Vector，SIRV）。基于SIRV的杂波建模的优点在于能够独立控制非高斯杂波的幅度PDF和杂波序列的相关性。文献[4?5]给出了基于SIRV的非高斯相关杂波的建模方法，利用该方法可以产生Weibull，K分布等常见非高斯分布相关杂波。但在实际雷达工作环境中，杂波序列不仅仅表现为时间的相关性，在不同距离单元上还表现为空间相关性[6]，这是由于雷达波束照射范围内杂波大尺度的变化特性产生的。而文献[5]给出的方法不能够实现空间相关，这是因为杂波建模中的随机变量[s]和零均值高斯矢量都是独立产生的。

1 基于SIRV的非高斯杂波建模

雷达杂波建模的要求是杂波序列的边缘概率密度函数和功率谱满足指定要求。产生非高斯相关序列的方法主要有无记忆非线性变换[7]（Zero Memory Non Linear， ZMNL）和SIRV两种方法。ZMNL方法是对相关高斯序列进行非线性变换，这种非线性变换是无记忆的，可以产生某些非高斯幅度分布，但非线性变换后的杂波功率谱会展宽，因此该方法适合描述具有弱相关的杂波序列，从而限制了该方法的应用；而基于SIRV的杂波建模是目前最具潜力的杂波模型，该模型的最大优点是可以独立控制杂波序列的相关性和幅度PDF。SIRV可以表示为一个非负的实的随机变量与一个复高斯矢量的乘积：

[c=τ×g] （1）

式中：随机变量[τ]的PDF为[f（τ），]在SIRV中通常称之为结构（texture）分量；零均值[N]维矢量[g]满足联合高斯分布[N（0，M），]其中[M]为协方差矩阵。

基于SIRV的建模方法就是要产生满足PDF为[f（τ）]的随机变量[τ]和满足协方差矩阵为[M]的零均值高斯矢量。要获得协方差矩阵为[M]的零均值高斯矢量，主要方法是考虑到[M]是Hermitian的，因此其对角化为：

[M=UΛUH] （2）

式中：矩阵[U]为酉矩阵；矩阵[Λ]为对角矩阵；对角线上的元素就是协方差矩阵的特征值，而矩阵[U]的列向量由这些特征值对应的特征矢量构成；上标[H]表示矩阵的共轭转置。由此可以得到矩阵[A=UΛ12UH，]假定零均值白色高斯矢量为[w，]那么就可以得到矢量：

[g=Aw] （3）

显然[EggH=EAwwHAH=M，]由此可以根据指定的协方差矩阵构造出相应的高斯分布矢量。获得指定PDF的随机变量[τ]的方法通常采用选舍法[8]（Rejection Method），该方法可以产生任意指定PDF的随机变量。其基本方法是同时产生一对均匀分布的随机变量[（x，y），]其中随机变量[x]在[f（τ）]的支撑域内均匀分布，而随机变量[y]在[f（τ）]的值域内均匀分布。如果满足[y

以上方法可以用于某个距离单元内的杂波序列在相干处理间隔（Coherent Processing Interval，CPI）的脉冲回波杂波建模，其中序列长度[N]表示脉冲累积的数目。但如果要描述不同距离单元之间的脉冲相关性，上述方法就失效了，这是因为在SIRV杂波建模过程中，随机变量[τ]和高斯矢量[g]是统计独立的，两个随机变量[τ]之间，以及两个高斯矢量[g]对应元素之间也是统计独立的。

假定利用上述方法产生[K]个[N]维SIRV矢量，那么就构成了一个[N×K]的矩阵，表示为[C=[c1，c2，…，cK]，]其中[cK]为[N]维SIRV列矢量。假定[cK]满足协方差矩阵[R]为[N]维方阵，该协方差矩阵描述了杂波一个CPI内的时间相关性。同样，[cK]之间的协方差矩阵[S]为[K]维方阵，该协方差矩阵描述了杂波不同距离单元之间的空间相关性。而由上述SIRV产生方法，只能指定协方差矩阵[R，]而仿真得到的协方差矩阵[S]始终是单位矩阵，即上述杂波建模方法不能描述杂波的空间相关性。空间相关性可以表示为：

[Eck1ck2=Eτk1gk1nτk2gk2n=Eτ1τ2Egngn ≡0， 1≤n≤N；1≤k1，k2≤K] （4）

这说明，要使得杂波序列满足时间相关性的同时，还要求满足空间相关性，必须保证公式（4）右边两个因子不能恒等于零，也就是说产生SIRV过程中的随机变量[τ]要满足一定的相关性，而且高斯序列[gk]要满足空时二维相关性。

2 空时二维相关高斯序列的产生

建立空时二维相关的高斯序列，本质上就是建立具有二维相关的高斯分布矩阵。这种矩阵的获得方法很多，而本文给出一种较为直接的算法。假定如下的矢量化表示：

[vecGN×K=vecg1g2…gK=vecgH1gH2?gHN=g1g2?gK=gNK×1] （5）

该函数表示将[N×K]维矩阵[G]的列矢量首尾相连，构成一个[NK]维列矢量，如果假定[gk～N（0，R），][gn～N（0，S），]那么容易得到[g～N（0，T），]其中矩阵[T]为[R]和[S]的Kronecker乘积，是[NK]维方阵：

[TNK×NK=R?S=RSij] （6）

由于矩阵[T]是两个协方差矩阵的Kronecker乘积，因此也是Hermitian的，可以利用类似式（2），式（3）的方法，得到满足协方差矩阵[T]的高斯序列，然后利用式（5）的逆变换得到[N×K]维矩阵。

将变换前后的二维随机序列的相关系数进行绘制，如图1所示。在变换前，相关系数是一个delta脉冲形状，说明序列在两个方向上是不相关的。而经过上述处理之后，两个方向的相关性得到了增强。

3 相关[τ]分量的产生

产生指定PDF，具有一定相关性的[τ]变量的方法则要复杂得多，这是因为当要考虑建立具有相关性的[τ]随机序列时，必须指定多变量PDF。通常而言，[f（τ）]是比较复杂的，采用ZMNL方法是不合适的，有时甚至不可能实现。这是因为每一类型的[f（τ）]惟一对应一种ZMNL方法。本文给出了一种较为简单的方法，该方法的基本思想是考虑到满足某个分布的随机变量的排序不会改变其幅度分布，但会改变相关性，因此可以先获得无相关性的随机序列，经过适当的排序，得到相关的序列，具体方法如图2所示。

图1 相关系数对比

图2 相关[τ]分量的产生

首先是利用选舍法获得幅度分布满足[f（τ）]的随机变量。由选舍法的基本原理可知，这些随机变量之间是统计独立的，在时序上没有任何相关性。同时，还考虑到随机变量在时序上调换位置，而不改变变量的大小并不影响幅度分布的规律。因此可以对其重新排序，使得在时序上满足某种相关性。由此可见要获得这种相关性，本质上就是要获得时序排列的方法。

本文假定独立同分布的高斯分布随机序列是可以获得的，考虑到高斯分布的线性变换仍然是高斯分布。值得指出的是[f（τ）]一般都是非高斯的，因此利用选舍法得到的随机序列经过线性滤波器，输出的幅度分布可能不再满足[f（τ），]因此可以将独立同分布的高斯分布序列通过线性滤波器，得到相关的高斯分布随机序列。输出高斯序列的相关性与滤波器的频率响应存在对应关系。对该相关高斯序列进行排序，就可以得到相关序列与顺序序列之间的位置映射关系[P：]

[P：1，2，…，K?k1，k2，…，kK] （7）

随机序列的排序过程，就是位置的映射过程，这种映射是一一映射，也是满射，因此其逆映射是存在的，定义为[P-1。]将这个逆映射作用于排序后的[τ]序列，就可以得到满足某种相关性的[τ]序列，其相关性与高斯序列的相关性存在某种关系，这种关系的确定较为复杂，但通过计算机仿真验证发现，高斯序列的相关性越强，获得的[τ]序列的相关性也越强。

图3给出了满足某种PDF的随机序列在变换前后的波形图。从波形图可以看出，改变随机序列的位置，可以得到满足某种相关性条件，同时不改变幅度分布特性。

图3 位置映射前后的序列波形

为了考察[τ]序列经过位置逆映射之后的相关性与相关高斯序列的相关性之间的关系，假定高斯序列满足的相关系数为对称指数分布：[ρ（m）=exp（-bm）]，其中，参数[b]控制了序列的相关性，[b]越大序列的相关性越小，如图4中虚线所示。

图4 高斯序列和[τ]的相关系数比较

图4给出了[m=1，2，3]时的相关系数随着[b]的变化规律。利用前述方法得到相关[τ]序列，其相关系数随[b]的变化规律如图4中实线所示。可以发现，[τ]序列经过位置变换之后，具有与高斯序列相接近的相关性，即随着[b]的增大，[τ]序列的相关性也减弱。但与高斯序列相关性比较可以发现，其相关系数比高斯序列的相关性弱，这主要是由于在位置变换这种逆映射过程中，等价于某种非线性变换过程，从而导致了功率谱的展宽，使得相关性减弱。

4 结 论

杂波序列通常表现为空时二维相关性，在时间方向上反映了CPI内某个距离单元回波之间的相关性，而不同距离单元的回波相关性则反映了雷达天线波束照射区域内的场景结构，而这种结构也导致了杂波序列具有空间相关性。已有的基于SIRV的杂波建模方法，无法描述杂波序列的空间相关性，这是由其产生方法决定的。

本文通过设计空时二维相关的高斯序列，以及采用位置映射的方法获得空间相关的[τ]序列，可以产生基于SIRV模型的空时二维相关的杂波序列。计算机仿真分别验证了高斯序列和[τ]序列的相关性。但本文并没有给出[τ]序列的相关函数与高斯序列相关函数之间的显性关系，这部分内容将是下一步研究的重点。

参考文献

[1] WATTS S. Radar sea clutter modelling and simulation?recent progress and future challenges [C]// 2008 IET Seminar on Radar Clutter Modelling. London： IET， 2008，1?7.

[2] ANATASSOPOULOS V， LAMPROPOULOS G A， DROSOPOULOS A， et al. High resolution radar clutter statistics [J]. IEEE Transactions on Aerospace and Electronic System， 1999， 35（1）： 43?60.

[3] CONTE E， LONGO M. Characterisation of radar clutter as a spherically invariant random process [J]. IEE Proceedings F： Communications， Radar and Signal Processing， 1987， 134（2）： 191?197.

[4] RANGSWAMY M， WEINER D D， OZTUK A. Non?Gaussian random vector identification using spherically invariant random processes [J]. IEEE Transactions on Aerospace and Electronic Systems， 1993， 29（1）： 111?123.

[5] RANGSWAMY M， WEINER D D. Computer generation of correlated non?Gaussian radar clutter [J]. IEEE Transactions on Aerospace and Electronic Systems， 1995， 31（1）： 106?116.

[6] RAGHAVAN R S. A model for spatially correlated radar clutter [J]. IEEE Transactions on Aerospace and Electronic Systems， 1991， 27（2）： 268?275.

[7] MARIER L J. Correlated K?distributed clutter generation for radar detection and track [J]. IEEE Transactions on Aerospace and Electronic Systems， 1995， 31（2）：568?580.

[8] CASELLA G， BERGER R. Statistical inference [M]. 2nd ed. New York： Duxbury Press， 2001.