首页 > 专利技术

一种高分辨率海杂波建模与仿真方法

1.本发明属于雷达信号处理领域，具体涉及一种幅度分布满足广义pareto(gp)分布的海杂波仿真算法，适用于弹载、机载、星载、岸基对海探测雷达系统的回波建模与仿真。在杂波仿真过程中，利用球不变随机过程(spherically invariant random process，sirp)法独立控制海杂波功率谱以及幅度模型的产生，增强了算法的参数调整灵活性和模型适应性。

背景技术：

2.海杂波的幅度特性在海面目标检测以及雷达系统仿真等应用中起着不可或缺的作用。然而，获取实测海杂波数据的途径是十分有限的，一方面雷达海杂波数据集一般是保密的，另一方面实测海杂波数据的收集需要耗费大量的人力、物力和财力。因此，研究具备时空相关的高分辨率海杂波仿真方法，能够为海杂波特性分析及目标检测算法研究提供充足的样本支撑。
3.随着雷达技术的不断发展，现代体制的雷达分辨率不断提升，这导致雷达海杂波具有极强的重拖尾和非高斯特性，严重脱离了传统的瑞利分布模型。大量研究表明，gp分布不仅能够较好的拟合高分辨率海杂波的幅度特性，而且能够描述海杂波的时空相关性，因此开展时空相关的gp分布杂波仿真对于最优雷达系统的仿真与设计具有重大意义。

技术实现要素：

4.要解决的技术问题
5.由于sirp算法能够独立控制海杂波序列功率谱以及幅度分布模型的仿真过程，更易于推广。因此，本发明采用sirp法来仿真具有高斯型功率谱的一维海杂波序列，利用递推算法在一维杂波序列的基础上仿真具有空时相关特性的二维海杂波序列。
6.技术方案
7.一种高分辨率海杂波建模与仿真方法，其特征在于步骤如下：
8.步骤1：构建海杂波幅度分布模型为：
[0009][0010]
其中：c表示符合gp分布的海杂波时间序列；τ表示慢变的纹理分量，符合逆gamma分布，概率密度函数pdf表示式为u表示快变的散斑分量，符合瑞利分布，条件pdf表示式为因此，gp分布的pdf被定义为：
[0011][0012]
其中：v表示gp分布的形状参数，决定其拖尾特性；b表示gp分布的尺度参数，决定其能量特性；
[0013]
步骤2：构建高斯型功率谱的傅里叶级数展开序列
[0014]
根据高斯型功率谱的函数表达式：
[0015][0016]
其中：s0表示任意常数；f表示频率；fd表示多普勒频偏；σc表示功率谱谱宽；利用傅里叶级数展开，构建高斯型功率谱的傅里叶级数系数b；
[0017]
步骤3：构建具有高斯型功率谱的幅度时间序列
[0018]
根据步骤2中所求高斯型功率谱的傅里叶级数系数b，具有高斯型功率谱的幅度时间序列u被定义为：
[0019]
x1＝b*u
[0020]
其中：x1表示具有高斯型功率谱的幅度时间序列，u表示均值为0，方差为1的复白高斯噪声序列；
[0021]
步骤4：构建具有空间相关性递减的高斯型功率谱幅度时间序列组
[0022]
根据步骤3构建相互独立的高斯型功率谱的幅度时间序列组gm(k),m＝1,2
……
，利用递推公式：获得m+1个距离单元空间衰减因子λ＝0.8的高斯型功率谱序列组x
m+1
(k),m＝1,2
……
；
[0023]
其中gm(k),m＝1,2
……
满足如下准则：
[0024]
(1)gm(k)是均值为0且能量与x1相同的复高斯随机向量：
[0025]
e(gm(k))＝0
[0026][0027]
其中：e表示对序列取均值；
[0028]
(2)gm(k)与xm(k)是相互独立的随机向量(rv)：
[0029]
e(xm(k)g
m+n
(k+l))＝0
[0030]
(3)gm(k)与xm(k)的自相关函数值相同：
[0031]rxx
(l)＝e(xm(k)xm(k+l))
[0032]
＝r
gg
(l)＝e(gm(k)gm(k+l))
[0033]
其中：k表示任意常数，r
xx
(l)表示序列xm(k)的自相关函数，r
gg
(l)表示序列gm(k)的自相关函数；
[0034]
(4)gm(k)在空间上不存在相关性：
[0035]
e(gm(k)g
m+l
(k))＝0,l≠0
[0036]
其中：l表示任意常数；
[0037]
步骤5：构建具有逆gamma分布幅度特性的时间序列
[0038]
高斯分布随机变量nv与逆gamma分布随机变量igv之间的非线性关系可以被定义为：
[0039][0040]
其中：表示不完全gamma函数；表示误差函数；u1(k)表示均值为0方差为1的实nv，z(k)表示形状参数为v的igv，即具有逆gamma分布幅度特性的时间序列；
[0041]
步骤6：构建时空相关的gp分布海杂波
[0042]
根据步骤1中gp分布海杂波的建模过程以及步骤4和5所分别构建的x
m+1
(k),m＝1,2
……
和形状参数为v的igv，具有时空相关的gp分布海杂波被定义为：
[0043][0044]
其中：u
m+1
(k)表示具有时空相关性的gp分布海杂波数据集；m+1表示海杂波数据集的第m+1距离单元；.*表示m+1个空间相关性递减的高斯谱序列x
m+1
(k),m＝1,2
……
分别与igvz(k)相乘；
[0045]
步骤7：仿真海杂波数据的空间相关性分析
[0046]
所述仿真海杂波数据的空间相关性为：
[0047][0048]
其中：r
xx
(0)表示海杂波序列的能量；xm(k)和x
m+l
(k)分别表示具有高斯型功率谱序列组的第m个和第m+l个距离单元的数据，根据计算结果可归纳任意两个距离单元仿真海杂波数据之间的空间相关性以λ＝0.8为衰减因子，呈指数倍衰减；
[0049]
步骤8：仿真海杂波数据的时间相关性分析
[0050]
所述仿真海杂波数据的时间相关性为：
[0051][0052]
其中：r
m+1,m+1
(l)表示第m+1个距离单元海杂波数据的自相关函数，r
xx
(l)表示初始距离单元的海杂波数据自相关函数；根据计算结果，任意距离单元与初始距离单元海杂波数据的时间相关性相同；
[0053]
步骤9：仿真海杂波序列的幅度特性分析
[0054]
根据步骤6中仿真的时空相关的海杂波数据，gp分布的幅度特性可以利用参数估计进行验证；gp分布的二/四阶矩估计被定义为：
[0055][0056]
其中，和分别表示仿真海杂波数据的2阶原点矩和4阶原点矩；和分别表示形状参数和尺度参数的估计值。
[0057]
一种计算机系统，其特征在于包括：一个或多个处理器，计算机可读存储介质，用于存储一个或多个程序，其中，当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现上述的方法。
[0058]
一种计算机可读存储介质，其特征在于存储有计算机可执行指令，所述指令在被执行时用于实现上述的方法。
[0059]
有益效果
[0060]
本发明提出的一种高分辨率海杂波建模与仿真方法，采用sirp算法以及空间相关性递推算法仿真了具有时空相关特性的二维海杂波时间序列，克服了实测海杂波数据收集困难以及保密的壁垒，有助于海上目标检测器以及雷达系统的设计与仿真。
[0061]
本发明基于sirp法仿真了时空相关的gp分布海杂波，相较于零记忆非线性变换法(zero memory non-linearity，zmnl)能够独立控制功率谱和幅度模型的产生过程，更易于推广。
[0062]
实验结果表明：仿真海杂波数据具有较好的幅度与时空相关特性，有助于海上目标探测和雷达系统的仿真。
附图说明
[0063]
附图仅用于示出具体实施例的目的，而并不认为是对本发明的限制，在整个附图中，相同的参考符号表示相同的部件。
[0064]
图1为sirp法仿真gp分布海杂波数据的流程。
[0065]
图2为一维gp分布海杂波幅度仿真序列及其拟合效果：(a)一维gp分布海杂波幅度序列；(b)一维仿真海杂波数据的拟合效果。
[0066]
图3为一维gp分布海杂波数据的功率谱。
[0067]
图4为二维时空相关gp分布海杂波数据集仿真：(a)一维gp分布海杂波幅度序列；(b)时空相关海杂波序列的三维图。
[0068]
图5为二维时空相关gp分布海杂波数据集的功率谱。
[0069]
图6为二维时空相关gp分布海杂波数据的空间相关衰减特性曲线图。
具体实施方式
[0070]
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图和实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。此外，下面描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
[0071]
本发明提出了一种时空相关gp分布的海杂波数据仿真算法，主要包括以下步骤：
[0072]
步骤1：构建海杂波的建模过程。根据海杂波的形成机理，海杂波幅度分布模型可被定义为：
[0073][0074]
其中：c表示符合gp分布的海杂波时间序列；τ表示慢变的纹理分量，符合逆gamma分布，概率密度函数(pdf)表示式为u表示快变的散斑分量，符合瑞利分布，条件pdf表示式为因此，gp分布的pdf被定义为：
[0075][0076]
其中：v表示gp分布的形状参数，决定其拖尾特性；b表示gp分布的尺度参数，决定其能量特性。在本发明中，仿真的海杂波数据的形状参数v＝3，尺度参数b＝1。
[0077]
步骤2：构建高斯型功率谱的傅里叶级数展开序列。根据高斯型功率谱的函数表达式：
[0078][0079]
其中：s0表示任意常数；f表示频率；fd＝100表示多普勒频偏；σc＝40表示功率谱谱
宽，利用傅里叶级数展开，构建高斯型功率谱的傅里叶级数系数b。
[0080]
步骤3：构建具有高斯型功率谱的幅度时间序列。根据步骤2中所求高斯型功率谱的傅里叶级数系数b，具有高斯型功率谱的幅度时间序列u被定义为：
[0081]
x1＝b*u
[0082]
其中：x1表示具有高斯型功率谱的幅度时间序列，u表示均值为0，方差为1的复白高斯噪声序列。
[0083]
步骤4：构建具有空间相关性递减的高斯型功率谱幅度时间序列组。根据步骤3构建相互独立的高斯型功率谱的幅度时间序列组gm(k),m＝1,2
……
，利用递推公式：可以获得m+1个距离单元空间衰减因子λ＝0.8的高斯型功率谱序列组x
m+1
(k),m＝1,2
……
。
[0084]
其中gm(k),m＝1,2
……
满足如下准则：
[0085]
(1)gm(k)是均值为0且能量与x1相同的复高斯随机向量：
[0086]
e(gm(k))＝0
[0087][0088]
其中：e表示对序列取均值。
[0089]
(2)gm(k)与xm(k)是相互独立的随机向量(rv)：
[0090]
e(xm(k)g
m+n
(k+l))＝0
[0091]
(5)gm(k)与xm(k)的自相关函数值相同：
[0092]rxx
(l)＝e(xm(k)xm(k+l))
[0093]
＝r
gg
(l)＝e(gm(k)gm(k+l))
[0094]
其中：k表示任意常数，r
xx
(l)表示序列xm(k)的自相关函数，r
gg
(l)表示序列gm(k)的自相关函数。
[0095]
(6)gm(k)在空间上不存在相关性：
[0096]
e(gm(k)g
m+l
(k))＝0,l≠0
[0097]
其中：l表示任意常数。
[0098]
步骤5：构建具有逆gamma分布幅度特性的时间序列。根据文献“modelling and simulation of non-rayleigh radar clutter”所揭示的高斯分布随机变量(nv)与gamma分布随机变量(gv)的非线性关系以及逆gamma分布随机变量(igv)与gv互为倒数关系，nv与igv之间的非线性关系可以被定义为：
[0099][0100]
其中：表示不完全gamma函数；表示误差函数；u1(k)表示均值为0方差为1的实nv，z(k)表示形状参数为v的igv。即分布幅度特性的时间序列。
[0101]
步骤6：构建时空相关的gp分布海杂波。根据步骤1中gp分布海杂波的建模过程以及步骤4和5所分别构建的x
m+1
(k),m＝1,2
……
和形状参数为v的igv，具有时空相关的gp分
布海杂波被定义为：
[0102][0103]
其中：u
m+1
(k)表示具有时空相关性的gp分布海杂波数据集；m+1表示海杂波数据集的第m+1距离单元；.*表示m+1个空间相关性递减的高斯谱序列x
m+1
(k),m＝1,2
……
分别与igvz(k)相乘。
[0104]
步骤7：仿真海杂波数据的空间相关性分析。根据步骤6中所揭示的海杂波仿真过程(空间相关的高斯谱序列分别调制同一参数的igvz(k))，仿真海杂波数据的空间相关性体现在高斯型功率谱序列组x
m+1
(k),m＝1,2
……
的空间相关性。因此，仿真海杂波数据的空间相关性为：
[0105][0106]
其中：r
xx
(0)表示海杂波序列的能量；xm(k)和x
m+l
(k)分别表示具有高斯型功率谱序列组的第m个和第m+l个距离单元的数据，根据计算结果可归纳任意两个距离单元仿真海杂波数据之间的空间相关性以λ＝0.8为衰减因子，呈指数倍衰减。
[0107]
步骤8：仿真海杂波数据的时间相关性分析。根据步骤6和步骤7，仿真海杂波序列的时间相关性同样体现在高斯型功率谱序列组x
m+1
(k),m＝1,2
……
的时间相关性。因此，仿真海杂波数据的时间相关性为：
[0108][0109]
其中：r
m+1,m+1
(l)表示第m+1个距离单元海杂波数据的自相关函数，r
xx
(l)表示初始距离单元的海杂波数据自相关函数。根据计算结果，任意距离单元与初始距离单元海杂波数据的时间相关性相同。
[0110]
步骤9：仿真海杂波序列的幅度特性分析。根据步骤6中仿真的时空相关的海杂波数据，gp分布的幅度特性可以利用参数估计进行验证。gp分布的二/四阶矩估计被定义为：
[0111][0112]
其中，和分别表示仿真海杂波数据的2阶原点矩和4阶原点矩；和分别表示形状参数和尺度参数的估计值。对仿真的每一距离单元海杂波数据进行参数估计，参数估计结果如表一所示。
[0113]
表一、参数估计结果
[0114][0115]
根据对形状参数和尺度参数估计结果，可以发现：所仿真的海杂波数据具有较好的幅度统计特性。
[0116]
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明公开的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。

技术特征：

1.一种高分辨率海杂波建模与仿真方法，其特征在于步骤如下：步骤1：构建海杂波幅度分布模型为：其中：c表示符合gp分布的海杂波时间序列；τ表示慢变的纹理分量，符合逆gamma分布，概率密度函数pdf表示式为u表示快变的散斑分量，符合瑞利分布，条件pdf表示式为因此，gp分布的pdf被定义为：其中：v表示gp分布的形状参数，决定其拖尾特性；b表示gp分布的尺度参数，决定其能量特性；步骤2：构建高斯型功率谱的傅里叶级数展开序列根据高斯型功率谱的函数表达式：其中：s0表示任意常数；f表示频率；f
d
表示多普勒频偏；σ
c
表示功率谱谱宽；利用傅里叶级数展开，构建高斯型功率谱的傅里叶级数系数b；步骤3：构建具有高斯型功率谱的幅度时间序列根据步骤2中所求高斯型功率谱的傅里叶级数系数b，具有高斯型功率谱的幅度时间序列u被定义为：x1＝b*u其中：x1表示具有高斯型功率谱的幅度时间序列，u表示均值为0，方差为1的复白高斯噪声序列；步骤4：构建具有空间相关性递减的高斯型功率谱幅度时间序列组根据步骤3构建相互独立的高斯型功率谱的幅度时间序列组g
m
(k),m＝1,2
……
，利用递推公式：获得m+1个距离单元空间衰减因子λ＝0.8的高斯型功率谱序列组x
m+1
(k),m＝1,2
……
；其中g
m
(k),m＝1,2
……
满足如下准则：(1)g
m
(k)是均值为0且能量与x1相同的复高斯随机向量：e(g
m
(k))＝0其中：e表示对序列取均值；(2)g
m
(k)与x
m
(k)是相互独立的随机向量(rv)：
e(x
m
(k)g
m+n
(k+l))＝0(3)g
m
(k)与x
m
(k)的自相关函数值相同：r
xx
(l)＝e(x
m
(k)x
m
(k+l))＝r
gg
(l)＝e(g
m
(k)g
m
(k+l))其中：k表示任意常数，r
xx
(l)表示序列x
m
(k)的自相关函数，r
gg
(l)表示序列g
m
(k)的自相关函数；(4)g
m
(k)在空间上不存在相关性：e(g
m
(k)g
m+l
(k))＝0,l≠0其中：l表示任意常数；步骤5：构建具有逆gamma分布幅度特性的时间序列高斯分布随机变量nv与逆gamma分布随机变量igv之间的非线性关系可以被定义为：其中：表示不完全gamma函数；表示误差函数；u1(k)表示均值为0方差为1的实nv，z(k)表示形状参数为v的igv，即具有逆gamma分布幅度特性的时间序列；步骤6：构建时空相关的gp分布海杂波根据步骤1中gp分布海杂波的建模过程以及步骤4和5所分别构建的x
m+1
(k),m＝1,2
……
和形状参数为v的igv，具有时空相关的gp分布海杂波被定义为：其中：u
m+1
(k)表示具有时空相关性的gp分布海杂波数据集；m+1表示海杂波数据集的第m+1距离单元；.*表示m+1个空间相关性递减的高斯谱序列x
m+1
(k),m＝1,2
……
分别与igvz(k)相乘；步骤7：仿真海杂波数据的空间相关性分析所述仿真海杂波数据的空间相关性为：其中：r
xx
(0)表示海杂波序列的能量；x
m
(k)和x
m+l
(k)分别表示具有高斯型功率谱序列组的第m个和第m+l个距离单元的数据，根据计算结果可归纳任意两个距离单元仿真海杂波数据之间的空间相关性以λ＝0.8为衰减因子，呈指数倍衰减；步骤8：仿真海杂波数据的时间相关性分析所述仿真海杂波数据的时间相关性为：
其中：r
m+1,m+1
(l)表示第m+1个距离单元海杂波数据的自相关函数，r
xx
(l)表示初始距离单元的海杂波数据自相关函数；根据计算结果，任意距离单元与初始距离单元海杂波数据的时间相关性相同；步骤9：仿真海杂波序列的幅度特性分析根据步骤6中仿真的时空相关的海杂波数据，gp分布的幅度特性可以利用参数估计进行验证；gp分布的二/四阶矩估计被定义为：其中，和分别表示仿真海杂波数据的2阶原点矩和4阶原点矩；和分别表示形状参数和尺度参数的估计值。2.一种计算机系统，其特征在于包括：一个或多个处理器，计算机可读存储介质，用于存储一个或多个程序，其中，当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现权利要求1所述的方法。3.一种计算机可读存储介质，其特征在于存储有计算机可执行指令，所述指令在被执行时用于实现权利要求1所述的方法。

技术总结

本发明涉及一种高分辨率海杂波建模与仿真方法，属于雷达信号处理领域。包括以下步骤：构建海杂波幅度分布模型；构建高斯型功率谱的傅里叶级数系数；构建具有高斯型功率谱的幅度时间序列；构建具有空间相关性递减的高斯型功率谱幅度时间序列组；构建具有逆gamma分布幅度特性的时间序列；构建时空相关的GP分布海杂波；分别仿真海杂波数据的空间相关性、时间相关性；最后对幅度特性进行验证。本发明采用SIRP算法以及空间相关性递推算法仿真了具有时空相关特性的二维海杂波时间序列，克服了实测海杂波数据收集困难以及保密的壁垒，有助于海上目标检测器以及雷达系统的设计与仿真。海上目标检测器以及雷达系统的设计与仿真。海上目标检测器以及雷达系统的设计与仿真。