CN109212512B - 一种具有空间相关性的水声阵列海洋环境噪声仿真方法 - Google Patents

Abstract

本发明公开了一种具有空间相关性的海洋环境噪声仿真方法，包括如下步骤：（1）设定仿真参数；（2）根据水声阵列参数求解随频率变化的海洋环境噪声空间相关函数，将离散频率进行子频带划分；（3）计算各子频带内相关系数参考值；（4）生成不相关的宽带高斯白噪声；（5）依照步骤2划分的子频带，对步骤4生成的不相关宽带噪声进行滤波，生成不相关带限噪声；（6）依照步骤3中得到的空间相关系数参考值，生成相关带限噪声。（7）将相关带限噪声相加，生成相关宽带噪声；（8）计算海洋环境噪声AR滤波器系数；（9）将步骤7得到的相关宽带噪声作为输入激励步骤8得到的滤波器，生成具有空间相关性的水声阵列海洋环境噪声。

Description

一种具有空间相关性的水声阵列海洋环境噪声仿真方法

技术领域

本发明涉及水声阵列信号仿真技术领域，特别是涉及一种具有空间相关性的水声阵列海洋环境噪声仿真方法。

背景技术

声纳等水声系统利用水听器阵列接收水中声信号，实现水中目标检测、定位、识别等功能。随着现代计算机性能的不断提高，越来越多的水声技术研究人员利用计算机仿真方式生成水声阵列接收信号，对水声信号处理算法以及系统进行测试评估，以降低海上试验成本，提高测试评估效率。为了达到这一目的，水声阵列接收信号仿真系统需要尽可能地还原海上试验情景，生成的仿真信号需要具有与海上实测信号良好的吻合度。水声阵列接收信号是随机信号，因此仿真生成的水声阵列接收信号的统计特性需要尽可能地与海上实测信号的统计特性相吻合。其中，水声阵列各阵元接收的海洋环境噪声的时间和空间统计特性就是水声阵列信号仿真的重要部分。

为提高信号处理的增益，现有的声纳系统通常利用水听器阵列接收水下声信号，并采用波束形成方法进行目标方位测量和目标跟踪。在波束形成方法中，一般假设各阵元接收的背景噪声服从高斯分布，且是不相关的，这样波束形成可获得最大的空间增益[1]。因此，大多数现有水声阵列信号仿真系统忽略了各阵元背景噪声之间的相关性[1][2][3]。但是，作为水声阵列主要背景噪声源之一的海洋环境噪声具有一定的空间相关性[4]，这使得阵列信号处理算法的仿真测试和实际使用之间存在性能上的差异。

海洋环境噪声是水声阵列背景噪声的重要组成部分，是一种宽带噪声。一方面，由于海洋环境噪声自身的多变性以及水声阵列工作环境的复杂性，建立起一套易于信号仿真生成使用的宽带噪声理论模型和应用模型难度很大。另一方面，受到计算能力的限制，对于海洋环境噪声仿真，在满足符合时频特征要求的情况下，很难兼顾阵元间的海洋环境噪声空间相关性。因此，建立既满足符合时频特征要求，又符合空间相关特性要求的宽带海洋环境噪声仿真生成方法，对于提升水声阵列信号仿真的逼真性和有效性具有重要意义。

本发明构建了水声阵列宽带海洋环境噪声阵元间空间相关序列，根据空间相关序列的特性对离散频率进行了频带分解，并构造了子频带内恒定空间相关系数矩阵。利用分频带方法生成符合各子频带相关特性的带限白噪声，再合成空间相关宽带高斯白噪声；然后将相关宽带高斯白噪声作为AR(Auto Regressive，自回归)滤波器激励信号，采用AR建模方式生成宽带有色海洋环境噪声生成，利用AR滤波器的线性特性，将高斯白噪声激励间的相关性线性映射为不同阵元海洋环境噪声间的相关性，实现具有空间相关性的水声阵列海洋环境噪声时域序列仿真生成。

发明内容

针对现有水声阵列信号仿真技术中存在的问题，本发明提供一种具有空间相关性的水声阵列海洋环境噪声仿真方法，该方法先求解分频带空间相关系数矩阵，生成符合各子频带相关特性的带限白噪声，然后将各子频带带限白噪声合成得到相关宽带高斯白噪声，并激励AR滤波器，最终仿真生成具有空间相关特性的水声阵列宽带海洋环境噪声，为达此目的，本发明提供一种具有空间相关性的水声阵列海洋环境噪声仿真方法，包括如下步骤：

(1)设定水声阵列海洋环境噪声仿真参数；

(2)根据水声阵列参数求解M个阵元两两之间随频率变化的海洋环境噪声空间相关函数序列Q_uv(f_p)，1≤u<v≤M，M为阵元数，p＝0,1,…,N_f-1，N_f为离散频率总点数；并依照水声阵列中间距最小的两个阵元对应的空间相关函数序列对离散频率序列f_p进行频带分解，得到H个子频带B_j，j＝1,2…,H；

(2.1)设定的水声阵列阵元坐标，代入(1)～(3)式计算水声阵列阵元两两之间的间距序列s_uv，阵元连线俯仰角序列γ_uv以及阵元连线的方位角序列ζ_uv，再代入(4)式计算对应的海洋环境噪声空间相关函数序列Q_uv(f_p)，其中z_uv＝z_v-z_u，x_uv＝x_v-x_u；

s_uv＝||p_v-p_u||₂ (19)

式中表示方向概率密度函数，i为虚数单位，θ表示俯仰角，表示方位角；

(2.2)依据得到的海洋环境噪声空间相关函数序列Q_uv(f_p)，选取s_uv的最小值对应的空间相关函数序列作为参考序列进行频带划分；

首先，依照空间相关函数序列的单调性，将离散频率序列f_p划分为K个区间；前K-1个区间内，各离散频率对应的相关函数值单调变化，且最大值和最小值的差不小于ΔQ；第K个区间内，各离散频率对应的相关函数最大值和最小值的差小于ΔQ；第1个离散频率区间表示为第2个离散频率区间表示为第k个离散频率区间表示为第K个离散频率区间表示为且有s_K＝N_f-1，k＝1,2,…,K；

然后，在每个频率区间内按照(5)式所示的离散频率序列下标间隔ΔN_k对各离散频率区间再进行等间隔划分，每个频率区间划分为N_k个子频带，N_k取值如(6)式所示，总共划分为H个子频带：

(5)式和(6)式中round表示四舍五入取整，Q_kmax为第k个区间内相关函数序列的最大值，Q_kmin为第k个子区间内相关函数序列的最小值，H和K的满足(7)式的关系；

(3)根据空间相关函数序列Q_uv(f_p)，计算各子频带内相关系数的参考值

(4)生成M个长度为N，且互不相关的高斯白噪声序列X_l(n)，l＝1,2,…,M，n＝0,1,2,…,N-1；

(5)依照步骤2中的频带分解结果，对步骤4中M个互不相关的宽带高斯白噪声序列X_l(n)进行带通滤波，生成H组带限高斯白噪声序列，每组中的M个带限高斯白噪声序列互不相关；

(6)依照步骤3中得到的各子频带内空间相关系数参考值，生成H组相关带限高斯白噪声序列，第j组中的M个带限高斯白噪声序列分别对应于M个阵元在B_j频带内的带限相关高斯白噪声序列，且两两之间具有特定的相关系数；

(7)将与每个阵元对应的H个相关带限高斯白噪声序列相加，生成与M个阵元对应的相关宽带高斯白噪声序列Y_l(n)；

(8)根据设定的海洋环境噪声功率谱参数，求解海洋环境噪声AR滤波器系数；

(9)用步骤7中得到的与各阵元对应的相关宽带高斯白噪声序列Y_l(n)作为输入，激励步骤8中得到的AR滤波器，生成具有空间相关性的水声阵列海洋环境噪声。

作为本发明进一步改进，步骤1设定如下参数：

声速c，采样频率f_s，离散频率点数N_f，频率间隔Δf＝f_s/N_f，水声阵列阵元数量M，依照坐标系设定各阵元坐标(x_l，y_l，z_l)，阵元坐标用矩阵P表示，向量p_l表示第l个阵元的坐标，取第一个阵元位置为坐标原点，即矩阵P可用(8)式表示；设定空间相关函数划分精度ΔQ，0＜ΔQ≤1；

设定海洋环境噪声功率谱参数σ、f_m、f_m0和λ。

作为本发明进一步改进，步骤3具体包括如下步骤：

(3.1)根据步骤2中得到的频带划分结果，对每个Q_uv(f_p)进行频带划分，得到各子频带内的空间相关函数序列其中

(3.2)根据式(10)求解空间各子频带内的相关系数参考值

作为本发明进一步改进，步骤4具体包括如下步骤：

(4.1)设定高斯白噪声数学期望μ，方差σ_g ²；

(4.2)采用Box-Muller算法生成M个点数为N的宽带高斯白噪声序列X_l(n)，l＝1,2,…,M，n＝0,1,2,…,N-1，并组成矩阵X_N×M＝[X₁,X₂,…,X_M]，该矩阵的每一列为N点的宽带高斯白噪声序列构成的列向量，矩阵各列之间不相关。

作为本发明进一步改进，步骤5具体包括如下步骤：

(5.1)根据步骤2的频带划分结果，采用窗函数法设计FIR滤波器，生成H个FIR带通滤波器系数序列h_j，h_j的通带为子频带B_j；

(5.2)将步骤4中生成的矩阵X_N×M的各列分别与带通滤波器系数序列h_j进行卷积运算，如(11)式所示，得到H个带限高斯白噪声矩阵每个带限高斯白噪声矩阵的各列之间是不相关的；

(11)式中表示矩阵的第l列，X_l表示矩阵X_N×M的第l列，表示卷积运算。

作为本发明进一步改进，步骤6具体包括如下步骤：

(6.1)根据步骤3中得到的各子频带内相关系数参考值构成如(12)式所示的子频带内空间相关系数矩阵

(6.2)对空间相关系数矩阵进行Cholesky分解，将空间相关系数矩阵分解为上三角矩阵的转置与其自身的乘积，如(13)式所示；

(6.3)将不相关带限高斯白噪声矩阵与上三角矩阵相乘，如(14)式所示，得到相关带限高斯白噪声矩阵 的第i₁列和i₂列之间的相关系数等于空间相关系数矩阵的第i₁行、i₂列的元素

作为本发明进一步改进，步骤7包括如下步骤：

对步骤6得到的H个相关带限高斯白噪声矩阵依照(15)式进行相加，得到相关宽带高斯白噪声矩阵Y_N×M；

作为本发明进一步改进，步骤8包括如下步骤：

(8.1)采用三参数模型法模拟海洋环境噪声的功率谱序列E(f_p)，如(16)所示；

(8.2)建立q_n阶AR滤波器的Yule-Walker方程，如(17)所示；

其中a(l_n)和b₀为滤波器系数，δ(k_n)为冲击函数；r_x(k_n)为E(f_p)对应的自相关函数序列，l_n＝1,2,…,q_n，k_n＝0,1,…,N-1；

(8.3)采用Levison-Durbin算法求解(17)式，得到q_n阶AR滤波器系数a(l_n)和b₀。

作为本发明进一步改进，步骤9包括如下步骤：

利用步骤7得到的相关宽带高斯白噪声Y_N×M激励步骤8求得的滤波器系数a(l_n)和b₀，将Y_N×M中的每一列Y_l(n)分别代入(18)式，得到输出序列Z_l(n)；

Z_l(n)组成的矩阵Z_N×M＝[Z₁(n),Z₂(n),…,Z_M(n)]即为具有空间相关性的宽带海洋环境噪声，Z_N×M的每一列对应为水声阵列的一个阵元接收的海洋环境噪声序列，各列间具有指定的相关特性。

本发明公开的一种具有空间相关性的宽带海洋环境噪声仿真生成的方法，首先计算水声阵列宽带海洋环境噪声阵元间空间相关序列，根据空间相关序列的特性对离散频率进行了频带分解，生成子频带内恒定空间相关系数矩阵。利用分频带方法生成符合各子频带相关特性的带限白噪声，再合成相关宽带高斯白噪声。宽带有色海洋环境噪声生成采用AR建模方式，将相关宽带高斯白噪声作为AR滤波器激励信号，利用AR滤波器的线性特性，将高斯白噪声激励间的相关性线性映射为海洋环境噪声间的相关性，实现具有空间相关性的水声阵列海洋环境噪声时域序列生成。

与现有技术相比，本发明公开的方法具有以下优点：本发明公开的方法将高斯白噪声激励间的相关性线性映射为海洋环境噪声间的相关性，使得海洋环境噪声空间相关模型清晰，易于仿真系统编程实现；采用频带分解与合成的方法，降低计算量，应用简单直接；本发明公开方法生成的海洋环境噪声通道间具有指定的相关性，更贴近真实的水声阵列海洋环境噪声。

附图说明

图1为本发明方法的实施流程图；

图2为本发明方法建立的空间坐标系；

图3为实施例中的海洋环境噪声空间相关函数序列；

图4为实施例中的空间相关系数参考值；

图5为实施例中的海洋环境噪声功率谱；

图6至图10分别为实施例中仿真得到的水声阵列阵元1至阵元5空间相关宽带海洋环境噪声功率谱；

图11为实施例中仿真得到的空间相关系数与理论值对比图；

图12为实施例中仿真得到的空间相关系数多次平均值与理论值对比图。

具体实施方式

下面结合附图与具体实施方式对本发明作进一步详细描述：

本发明提供一种具有空间相关性的水声阵列海洋环境噪声仿真方法，该方法先求解分频带空间相关系数矩阵，生成符合各子频带相关特性的带限白噪声，然后将各子频带带限白噪声合成得到相关宽带高斯白噪声，并激励AR滤波器，最终仿真生成具有空间相关特性的水声阵列宽带海洋环境噪声。

实施例1：

下面结合附图和具体实施方式，进一步阐明本发明。

一种具有空间相关性的水声阵列宽带海洋环境噪声仿真生成方法，如图1所示，包括如下步骤：

步骤1，设定仿真参数

设定声速c＝1500m/s，采样频率f_s＝15kHz，离散频率点数N_f＝150000，频率间隔Δf＝f_s/N_f＝0.1Hz，水声阵列阵元数M＝5，ΔQ取值为0.1。本实施例设定水声阵列为等间距线列阵，阵元坐标矩阵P表示如下：

设定海洋环境噪声功率谱参数f_m＝200，f₀＝500，λ＝0，σ＝0.01。

步骤2，根据水声阵列阵元坐标，求解M个阵元两两之间随频率变化的海洋环境噪声空间相关函数序列Q_uv(f_p)，并选取序列s_uv最小值对应的空间相关函数序列作为参考序列，进行频带划分，包括如下步骤：

(2.1)依据设定的水声阵列阵元坐标P，代入(2)式～(4)3计算水声阵列阵元间距序列s_uv，阵元连线俯仰角序列γ_uv以及阵元连线的方位角序列ζ_uv，结果如表1所示。

表1阵元间距、连线俯仰角、连线方位角求解结果

表1结果代入(5)式求得阵元间空间相关函数序列Q_uv(f_p)。本实施例中，概率密度函数如(19)式所示。

其中π/4≤θ₁≤π/2，u₁>0，α₁以及a₁满足(20)式，I₀(·)是零阶第一类贝塞尔函数。

本例中，取θ₁＝π/3,u₁＝100，α₁＝10。联合(2)～(4)、(19)～(20)式求解得到水声阵列空间相关函数序列Q_uv(f_p)，组如图3所示。

(2.2)依据得到的海洋环境噪声空间相关函数序列Q_uv(f_p)，选取s_uv的最小值对应的空间相关函数序列作为参考序列进行频带划分。首先，依照空间相关函数序列的单调性，将离散频率序列f_p划分为K个区间；前K-1个区间内，各离散频率对应的相关函数值单调变化，且最大值和最小值的差不小于ΔQ；第K个区间内，各离散频率对应的相关函数最大值和最小值的差小于ΔQ。第1个离散频率区间表示为第2个离散频率区间表示为第k个离散频率区间表示为第K个离散频率区间表示为且有s_K＝N_f-1，k＝1,2,…,K。本实施例中求得的K＝3，H＝15，频带划分结果如表2所示。

表2频带划分结果

步骤3，求解各阵元间各子频带内相关系数参考值，包括如下步骤：

(3.1)根据表2所示的频带划分结果，对每个Q_uv(f_p)进行频带划分，得到子得到各子频带内的空间相关函数序列

(3.2)根据(10)式求解空间相关函数序列各子频带内的相关系数参考值本实施例求得的空间相关系数参考值如表3所示，作图如图4所示。

表3空间相关系数参考值

步骤4包括如下步骤：

(4.1)设定高斯白噪声数学期望μ＝0，方差σ_g ²＝1.

(4.2)采用Box-Muller算法，模拟生成5个点数为15000的宽带高斯白噪声序列X_l(n)，l＝1,2,…,5，n＝0,1,2,…,14999，组成矩阵X_15000×5＝[X₁,X₂,…,X_M]，该矩阵的每一列为15000点的宽带高斯白噪声序列构成的列向量，矩阵各列之间不相关。

步骤5包括如下步骤：

(5.1)根据表2的频带划分结果，采用函数法设计FIR滤波器，窗函数选用汉明窗，每个滤波器为512阶，生成15个FIR带通滤波器h_j。

(5.2)将步骤(4)中生成的矩阵X_15000×5的各列分别与带通滤波器系数序列h_j进行卷积运算，得到15个带限高斯白噪声矩阵每个带限高斯白噪声矩阵的各列之间是不相关的。

步骤6，生成多通道相关带限高斯白噪声激励信号，包括如下步骤：

(6.1)根据表3各子频带内相关系数参考值生成子频带内空间相关系数矩阵例如本实施例中第3个子频带内空间相关系数矩阵如(21)所示。

(6.2)对空间相关系数矩阵进行Cholesky分解，将空间相关系数矩阵分解为上三角矩阵的转置与其自身的乘积，例如本实施例中第3个子频带内上三角矩阵如(22)所示。

(6.3)将不相关带限高斯白噪声矩阵与上三角矩阵相乘，得到相关带限高斯白噪声矩阵

步骤7包括如下步骤：

对步骤6得到的15个相关带限高斯白噪声矩阵依照(15)式进行相加，得到相关宽带高斯白噪声矩阵Y_15000×5。

步骤8包括如下步骤：

(8.1)采用三参数模型法模拟海洋环境噪声的功率谱序列E(f_p)，功率谱E(f_p)如图5所示。

(8.2)建立q_n阶AR滤波器的Yule-Walker方程，r_x(k_n)为E(f_p)对应的自相关函数序列，可利用功率谱E(f_p)通过维纳-辛钦定理求得，l_n＝1,2,…,q_n，k_n＝0,1,…,15000-1；本实施例中，q_n＝8。

(8.3)采用Levison-Durbin算法求解Yule-Walker方程，得到q_n阶AR滤波器系数如表4所示。

表4.AR模型系数

步骤9包括如下步骤：

利用步骤7得到的相关宽带高斯白噪声激励Y_15000×5以及表4的系数a(l_n)，l_n＝1,2,…,8和b₀，将Y_15000×5中的每一列Y_l(n)分别代入(18)式，得到输出序列Z_l(n)。

本例仿真生成的具有空间相关特性的海洋环境噪声序列长度为1秒，对5个阵元的功率谱进行多次平均，结果如图6至图10所示。

为验证水声阵列宽带海洋环境噪声的相关特性，对海洋环境噪声依据表2频带划分进行分频带相关系数求解。仿真共进行10次，仿真结果如图11所示，曲线为理论的相关系数，散点为每次仿真的空间相关系数求解结果，可以看出每次仿真结果在理论值附近波动。对10次仿真的每个子频带内的相关系数进行平均，对比理论的相关系数与仿真的平均海洋环境噪声相关系数如图12所示，可以看出二者之间的相关性一致。

仿真结果表明，采用以上方法能够有效地实现具有相关性的水声阵列宽带海洋环境噪声仿真生成，当空间相关系数仿真精度取为0.1时，仿真生成的海洋环境噪声的阵元间相关特性以及频谱特性与理论值一致。

以上所述，仅是本发明的较佳实施例而已，并非是对本发明作任何其他形式的限制，而依据本发明的技术实质所作的任何修改或等同变化，仍属于本发明所要求保护的范围。

Claims (9)

1.一种具有空间相关性的水声阵列海洋环境噪声仿真方法，包括如下步骤，其特征在于：

(1)设定水声阵列海洋环境噪声仿真参数；

(2)根据水声阵列参数求解M个阵元两两之间随频率变化的海洋环境噪声空间相关函数序列Q_uv(f_p)，1≤u<v≤M，M为阵元数，p＝0,1,…,N_f-1，N_f为离散频率总点数；并依照水声阵列中间距最小的两个阵元对应的空间相关函数序列对离散频率序列f_p进行频带分解，得到H个子频带B_j，j＝1,2…,H；

(2.1)设定的水声阵列阵元坐标，代入(1)～(3)式计算水声阵列阵元两两之间的间距序列s_uv，阵元连线俯仰角序列γ_uv以及阵元连线的方位角序列ζ_uv，再代入(4)式计算对应的海洋环境噪声空间相关函数序列Q_uv(f_p)，其中z_uv＝z_v-z_u，x_uv＝x_v-x_u；

s_uv＝||p_v-p_u||₂ (1)

式中表示方向概率密度函数，i为虚数单位，θ表示俯仰角，表示方位角；

(2.2)依据得到的海洋环境噪声空间相关函数序列Q_uv(f_p)，选取s_uv的最小值对应的空间相关函数序列作为参考序列进行频带划分；

首先，依照空间相关函数序列的单调性，将离散频率序列f_p划分为K个区间；前K-1个区间内，各离散频率对应的相关函数值单调变化，且最大值和最小值的差不小于ΔQ；第K个区间内，各离散频率对应的相关函数最大值和最小值的差小于ΔQ；第1个离散频率区间表示为第2个离散频率区间表示为第k个离散频率区间表示为第K个离散频率区间表示为且有s_K＝N_f-1，k＝1,2,…,K；

然后，在每个频率区间内按照(5)式所示的离散频率序列下标间隔ΔN_k对各离散频率区间再进行等间隔划分，每个频率区间划分为N_k个子频带，N_k取值如(6)式所示，总共划分为H个子频带：

(5)式和(6)式中round表示四舍五入取整，Q_kmax为第k个区间内相关函数序列的最大值，Q_kmin为第k个子区间内相关函数序列的最小值，H和K的满足(7)式的关系；

(3)根据空间相关函数序列Q_uv(f_p)，计算各子频带内相关系数的参考值

(4)生成M个长度为N，且互不相关的高斯白噪声序列X_l(n)，l＝1,2,…,M，n＝0,1,2,…,N-1；

(5)依照步骤2中的频带分解结果，对步骤4中M个互不相关的宽带高斯白噪声序列X_l(n)进行带通滤波，生成H组带限高斯白噪声序列，每组中的M个带限高斯白噪声序列互不相关；

(6)依照步骤3中得到的各子频带内空间相关系数参考值，生成H组相关带限高斯白噪声序列，第j组中的M个带限高斯白噪声序列分别对应于M个阵元在B_j频带内的带限相关高斯白噪声序列，且两两之间具有特定的相关系数；

(7)将与每个阵元对应的H个相关带限高斯白噪声序列相加，生成与M个阵元对应的相关宽带高斯白噪声序列Y_l(n)；

(8)根据设定的海洋环境噪声功率谱参数，求解海洋环境噪声AR滤波器系数；

(9)用步骤7中得到的与各阵元对应的相关宽带高斯白噪声序列Y_l(n)作为输入，激励步骤8中得到的AR滤波器，生成具有空间相关性的水声阵列海洋环境噪声。

2.根据权利要求1所述的一种具有空间相关性的水声阵列海洋环境噪声仿真方法，其特征在于：步骤1设定如下参数：

声速c，采样频率f_s，离散频率点数N_f，频率间隔Δf＝f_s/N_f，水声阵列阵元数量M，依照坐标系设定各阵元坐标(x_l，y_l，z_l)，阵元坐标用矩阵P表示，向量p_l表示第l个阵元的坐标，取第一个阵元位置为坐标原点，即矩阵P可用(8)式表示；设定空间相关函数划分精度ΔQ，0＜ΔQ≤1；

设定海洋环境噪声功率谱参数σ、f_m、f_m0和λ。

3.根据权利要求1所述的一种具有空间相关性的水声阵列海洋环境噪声仿真方法，其特征在于：步骤3具体包括如下步骤：

(3.1)根据步骤2中得到的频带划分结果，对每个Q_uv(f_p)进行频带划分，得到各子频带内的空间相关函数序列其中

(3.2)根据式(10)求解空间各子频带内的相关系数参考值

4.根据权利要求1所述的一种具有空间相关性的水声阵列海洋环境噪声仿真方法，其特征在于：步骤4具体包括如下步骤：

(4.1)设定高斯白噪声数学期望μ，方差σ_g ²；

(4.2)采用Box-Muller算法生成M个点数为N的宽带高斯白噪声序列X_l(n)，l＝1,2,…,M，n＝0,1,2,…,N-1，并组成矩阵X_N×M＝[X₁,X₂,…,X_M]，该矩阵的每一列为N点的宽带高斯白噪声序列构成的列向量，矩阵各列之间不相关。

5.根据权利要求1所述的一种具有空间相关性的水声阵列海洋环境噪声仿真方法，其特征在于：步骤5具体包括如下步骤：

(5.1)根据步骤2的频带划分结果，采用窗函数法设计FIR滤波器，生成H个FIR带通滤波器系数序列h_j，h_j的通带为子频带B_j；

(5.2)将步骤4中生成的矩阵X_N×M的各列分别与带通滤波器系数序列h_j进行卷积运算，如(11)式所示，得到H个带限高斯白噪声矩阵每个带限高斯白噪声矩阵的各列之间是不相关的；

(11)式中表示矩阵的第l列，X_l表示矩阵X_N×M的第l列，表示卷积运算。

6.根据权利要求1所述的一种具有空间相关性的水声阵列海洋环境噪声仿真方法，其特征在于：

步骤6具体包括如下步骤：

(6.1)根据步骤3中得到的各子频带内相关系数参考值构成如(12)式所示的子频带内空间相关系数矩阵

(6.2)对空间相关系数矩阵进行Cholesky分解，将空间相关系数矩阵分解为上三角矩阵的转置与其自身的乘积，如(13)式所示；

(6.3)将不相关带限高斯白噪声矩阵与上三角矩阵相乘，如(14)式所示，得到相关带限高斯白噪声矩阵 的第i₁列和i₂列之间的相关系数等于空间相关系数矩阵的第i₁行、i₂列的元素

7.根据权利要求1所述的一种具有空间相关性的水声阵列海洋环境噪声仿真方法，其特征在于：步骤7包括如下步骤：

对步骤6得到的H个相关带限高斯白噪声矩阵依照(15)式进行相加，得到相关宽带高斯白噪声矩阵Y_N×M；

8.根据权利要求1所述的一种具有空间相关性的水声阵列海洋环境噪声仿真方法，其特征在于：步骤8包括如下步骤：

(8.1)采用三参数模型法模拟海洋环境噪声的功率谱序列E(f_p)，如(16)所示；

(8.2)建立q_n阶AR滤波器的Yule-Walker方程，如(17)所示；

其中a(l_n)和b₀为滤波器系数，δ(k_n)为冲击函数；r_x(k_n)为E(f_p)对应的自相关函数序列，l_n＝1,2,…,q_n，k_n＝0,1,…,N-1；

(8.3)采用Levison-Durbin算法求解(17)式，得到q_n阶AR滤波器系数a(l_n)和b₀。

9.根据权利要求1所述的一种具有空间相关性的水声阵列海洋环境噪声仿真方法，其特征在于：步骤9包括如下步骤：

利用步骤7得到的相关宽带高斯白噪声Y_N×M激励步骤8求得的滤波器系数a(l_n)和b₀，将Y_N×M中的每一列Y_l(n)分别代入(18)式，得到输出序列Z_l(n)；

Z_l(n)组成的矩阵Z_N×M＝[Z₁(n),Z₂(n),…,Z_M(n)]即为具有空间相关性的宽带海洋环境噪声，Z_N×M的每一列对应为水声阵列的一个阵元接收的海洋环境噪声序列，各列间具有指定的相关特性。