SIRP法相干相关K分布雷达杂波的建模与仿真.doc

SIRP法相干相关K分布雷达杂波的建模与仿真 gjj_所谓杂波仿真，实际上就是要生成一系列在幅度上服从特定的概率密度分布（pdf）的相关随机序列，常见的杂波仿真方法有两种：零记忆非线性变换法（ZMNL）和 球不变随机过程法（SIRP）。ZMNL方法的基本思想是：首先产生相关的高斯随机过程，然后经过某种非线性变换得到所求的相关随机序列。这种方法的缺点就是输入序列与输出序列间有复杂的非线性关系,因此必须寻找输入序列与输出序列的相关函数间的非线性对应关系。SIRP方法的基本思想是：产生一个相关的高斯随机过程，然后用具有所要求的单点概率密度函数的随机序列进行调制。这种方法的缺点则是受所求的序列的阶数及自相关函数的限制,同时这种方法的计算量非常大,不易形成快速算法。ISAR是一种相干雷达，其海杂波必然是相干且时空相关的。对于相干相关杂波，以往的方法都是将非相干的ZMNL方法加以推广得到相干的ZMNL模型。这种方法得以应用的一个前提是已知非线性变换前后杂波相关系数的非线性关系，然而对于相干相关K分布杂波却很难找到这样一种非线性变换，于是我们采取SIRP方法来仿真ISAR的海杂波。K分布适用于描述高分辨雷达的非均匀杂波，多用于对海杂波的模拟。K分布可以由一个均值是慢变化的瑞利分布来表示，其中这个慢变化的均值服从分布。K分布的概率密度函数为：其中，是形状参数，是尺度函数，是伽马函数，是第二类修正贝赛尔函数。杂波平均功率，和之间的关系可表示为：对于大多数杂波来说，形状参数的取值范围是，对于较小的的取值，如时，杂波有较长的托尾，时的分布接近于瑞利分布。图1给出了K分布杂波序列的实现结构。图1 相干相关K分布杂波SIRP方法图中，为一复高斯白噪声，线性滤波器由的相关函数设计决定，为一与相互独立的实高斯噪声，线性滤波器必须使得输出的高斯序列具有高度的相关性(相关函数接近于1)，ZMNL变换使得输出的的概率密度函数（pdf）为杂波的特征pdf。对于K分布来说，服从广义分布，该分布的定义如下：要用图1所示的模型产生K分布杂波，需要产生符合广义K分布的并设计线型滤波器1和线型滤波器2。滤波器1的设计比较简单，它使输出具有所要产生杂波的功率谱，设计方法同ZMNL法的滤波器设计。由于我们对的相关函数不感兴趣，因此，可将滤波器2设计为一带宽很窄的低通滤波器，使得非线性变换随机序列的功率谱足够窄。下面以一例MATLAB仿真说明上述产生K分布杂波的过程 1。具体程序见“Matlab程序”文件夹K_distribution.m和nonline_eq_sirp.m。 例：产生杂波的幅度概率密度函数的参数为，功率谱密度为高斯谱，其3dB带宽为40Hz的K分布杂波，滤波器1的设计采用傅立叶级数展开法，模拟的杂波的功率谱密度采用Burg法估计得到。一、 复高斯白噪声的产生二、 滤波器的设计傅立叶级数展开法2这种方法是通过将所希望的网络的频率特性展成傅立叶级数的方法求滤波器加权系数的，故称这种方法为傅立叶级数展开法。众所周知，非递归滤波器可由以下差分方程来描述：式中：表示滤波器的第个输入；表示滤波器的第个输出；为滤波器加权系数。滤波器的传递函数可通过变换求出：频率响应为：（6）式为数字滤波器的频率响应，令，将（6）转化为模拟滤波器的频率响应：其中为将模拟滤波器转化成数字滤波器时的抽样间隔（这里的的单位也是频率单位，因为是在频域抽样），为其抽样频率，为模拟滤波器的频域周期。又已知，杂波归一化的高斯谱密度为：希望在输入白噪声时，有：显然，所设计滤波器应有高斯响应：将其展成傅立叶级数 这里要知道是以Fs为周期的频谱函数。：又由于为偶函数，所以：其中：对式(7)取绝对值 对比式(11).，根据谱的偶函数特性知，式(12)中的便等于式(7)中的，即非递归滤波器频率响应的傅立叶级数展开式的系数，就是该滤波器的加权系数。由于频率响应是给定的，于是使问题简单了。 为了求系数，改变变量，将的傅立叶变换写成：将式(10)代入，得：当有限时，傅立叶级数的系数 这里应用了周期信号的傅立叶变换与周期信号的傅立叶级数展开系数之间的关系3。式中为抽样间隔，至少应为杂波频谱宽度的倒数。这样，在高斯谱已知的情况下，非递归滤波器的加权系数就由完全确定了。该滤波器的主要特点：(1) 首先，它具备非递归滤波器的优点，结构简单，运算速度快。(2) 便于硬件实现，特别适用于雷达模拟器。(3) 要得到一个较好的响应，N值应大于8。(4) 这种方法对于输出序列的长度没有限制，取决于输入序列的长度。这对雷达系统的性能测试具有重要意义，例如对虚警概率进行测试时，应给出足够长的序列，如虚警概率时，其长度应大于。(5) 我们知道描述数字滤波器的差分方程时稳态情况下的差分方程，在输入序列小于它的阶数时，输出序列仍处于暂态期，它们不满足给定的统计特性。因此，在将其用于雷达模拟器时必须控制暂态输出。下面设计一个非递归高斯滤波器，程序见“Matlab程序”文件夹fseries.m。由例子知，设采样频率（杂波模拟时的脉冲重复频率）。假定在傅立叶展开式中，经计算，系数的数值给在表1中。实际上，在表中也给出了的数值，可以看出，对频率响应的贡献已经是很小了。表1 用傅立叶级数展开法获得的非递归滤波器的加权系数C00.138 473C70.014 470C10.260 755C80.005 862C20.217 645C90.002 105C30.161 043C100.000 670C40.105 635C110.000 189C50.061 426C120.000 047C60.031 665C130.000 010说明(1) (2) 采样频率根据此式所计算的功率谱密度曲线与理论值的差异如图2所示。从图2可以看出，当N取9时，说得到的功率谱曲线1与理论高斯谱模型曲线2重合在一起了；当N取3时，功率谱密度曲线3要比理论高斯曲线2宽，并且在高端有小的起伏振荡。计算表明，对于当N8时，再增加谐波次数，效果并不明显，由表1看到，这时因为系数的贡献太小的原因。图2 滤波器输出的归一化的功率谱密度曲线三、 滤波器的设计以及广义分布变量的产生4对于滤波器2的设计，由于通常非线性变换会使随机序列的功率谱展宽，所以，应使滤波器2的带宽足够窄，不致使非线性变换后输出的随机序列的功率谱展宽，以致序列的相关函数不再接近于1，因而对的要求主要是针对带宽的。由于我们并不对控制的相关函数感兴趣，因此，可将滤波器2设计为一带宽很窄的低通滤波器，使得非线性变换随机序列的功率谱仍足够窄。在本文例子的设计中，对于滤波器2，我们采用5阶butterworth低通滤波器。应用SIRP的K分布产生模型，必须产生广义分布随机变量，由于的平方即为伽马分布，所以可以产生伽马随机变量，再对它求平方根得到。文献5 感谢大树文献服务网：友情提供此篇文献的下载给出了伽马分布随机序列产生的方法，然而比较复杂。对于图1所示的结构，文献6给出了非线性变换的表达式如下：式中：为不完全欧拉函数(Incomplete Eulerian Function) 这个函数的具体表达形式我还没有找到，谁找到的话请共享一下.。文献7给出的表达式如下：式中为不完全伽马函数。为滤波器的输出，是K分布形状参数，为K分布尺度参数。由式(2)及可知，式(18)可变化为：为标准正态随机变量的尾部面积，即有：将式(20)代入式(19)，并应用概率密度在全区间积分为1，得：式中：是误差函数，定义为：因此，产生变量的问题转化为求式(21)的问题，这是一个非线性方程，可以用二分法求解。为了提到仿真时的速度，也可以提前对此非线性变换进行拟合，实际仿真时调用拟合得到的曲线关系即可。 后记 关于空域相关的海杂波仿真正在进一步研究之中。回忆自己当初接触海杂波到现在写这个总结，中间经历了很多，有迷惑，也有失落，但在研学论坛的帮助下，在与我认识的同学(比如hahnc和一位电子科技大学的同学（抱歉，不知道你研学上的id）)的讨论中，我逐渐对海杂波有了初步的认识，只是最初步的认识。为了感谢研学论坛，也为了总结一下以前的认识，特写了这篇总结(主要是摘抄，再组合，J)，聊以共享。希望各位研友再接再厉，在自己的研究道路上有新的成就。参考文献1. 罗军辉，et al. Matlab 7.0在数字信号处理中的应用. 北京：机械工业出版社，2005:180-184.2. 杨万海. 雷达系统建模与仿真. 西安：西安电子科技大学出版社，2007：111-113.3. 郑君里，et al. 信号与系统（第二版）.北京：高等教育出版社，2000：145-149.4. 吕雁，史林. SIRP法相干相关K分布雷达杂波的建模与仿真. 现代雷达, 2002，24(2): 13-16. 5. Oliver, C. J.; Tough, R. J. A. On The Simulation of Correlated K-Distributed Random Clutter. Opt. Acta, 1986,33(3):223-250.6. E. Conte, M. Longo, and M. Lops. Modelling