冲击噪声下的任意阵列相干源测向方法.pdf

《冲击噪声下的任意阵列相干源测向方法.pdf》由会员分享，可在线阅读，更多相关《冲击噪声下的任意阵列相干源测向方法.pdf（15页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201910324483.3 (22)申请日 2019.04.22 (71)申请人 哈尔滨工程大学 地址 150001 黑龙江省哈尔滨市南岗区南 通大街145号哈尔滨工程大学科技处 知识产权办公室 (72)发明人 高洪元杜亚男孙贺麟刁鸣 谢婉婷马雨微王世豪 (51)Int.Cl. G01S 3/14(2006.01) (54)发明名称 一种冲击噪声下的任意阵列相干源测向方 法 (57)摘要 本发明公开一种冲击噪声下的任意阵列相 干源测向方法， 包括： 建立采样信号模型； 构造真 。

2、实阵列动态随机加权低阶协方差矩阵； 定义内插 变换矩阵T， 构建虚拟阵列协方差矩阵； 获得前后 向空间平滑修正后的数据协方差矩阵和噪声协 方差矩阵， 预白化处理得到动态随机加权低阶协 方差矩阵； 估计信源个数， 对动态随机加权协方 差矩阵进行特征分解， 确定信号子空间和噪声子 空间； 构建动态随机加权低阶协方差-空间平滑- MUSIC测向方法的谱估计公式， 进行谱峰搜索， 找 出极大值点对应的角度， 输出任意阵列相干源测 向结果。 本发明能够对任意阵列的信源来波方向 进行有效估计， 可在高斯噪声、 弱冲击噪声和强 冲击噪声下进行测向， 解相干性能优， 应用范围 广泛。 权利要求书2页 说明书7。

3、页 附图5页 CN 110007266 A 2019.07.12 CN 110007266 A 1.一种冲击噪声下的任意阵列相干源测向方法， 其特征在于， 包括以下步骤： 步骤一： 建立任意阵列相干信源的采样信号模型： 对于一个阵元数为N的任意阵列， 以原点为参考点， 假设第k个阵元的位置为(xk,yk)， 第 i个远场窄带信号从 i方向入射到N元任意阵列， 其中， k1,2,.,N； i1,2,.,M； i表示 第 i 个 远 场 窄 带 入 射 信 号 的 方 位 角 ， 则 在 t 时 刻 第 k 个 阵 元 接 收 到 的 信 号 为 其中，ki(xkcos i+yksin i)/c；。

4、 k1,2,.,N； i1,2,.,M； w02 c/ ， c为光速， 为波长； si(t)为t时刻的第i个信源的入射信号， 当入射信号为相干信 号时， 第个生成源生成的相干源组为其中，I0 为生成源总数，为第 个生成源生成个信号的权重系数矢量， 则有 nk(t)为t时刻 第k个阵元的噪声信号， nk(t)是满足S S稳定分布的复冲击噪声； k1,2,.,N， 则阵列接 收的第l次快拍相干信号模型可表示为式中， A( )a( 1) ,a( 2) ,a( M)NM为阵列流型矩阵， 其中, 为虚数单位； 1, 2,., M为来波方向方位角矢 量； s(l)为信号矢量；为满足S S稳定分布的复冲击阵。

5、列噪声矢 量； T表示转置； 步骤二： 利用阵列接收到的数据构造真实阵列的动态随机加权低阶协方差矩阵： 利 用 阵 列 接 收 的 数 据 构 造 真 实 阵 列 的 动 态 随 机 加 权 低 阶 协 方 差 矩 阵 式中，其中， k1,2,.,N， L为最大快拍数； 为动态变化系数； b为加权系数； l为0,1之间均匀分布 的随机数， H表示共轭转置； 步骤三： 定义内插变换矩阵T， 构建阵列信号内插变换后的虚拟阵列低阶协方差矩阵： 首先对虚拟变换的方位角度区域进行划分， 假设信号位于区域 内， 将区域均匀划分 为和分别为虚拟变换方位角度区域的左右边界， 为 虚 拟 变 换 方 位 角 度。

6、 变 化 步 长 ， 则 在 区 域内 的 真 实 阵 列 流 行 矩 阵 为 为虚拟阵列的阵元个数， 则 在同一区域内， 假设虚拟阵列为阵元间距相等的均匀线阵， 虚拟阵列的阵列流型矩阵为 其中， 权利要求书 1/2 页 2 CN 110007266 A 2 根据真实阵列流型矩阵与虚拟阵列流型矩阵定义 变换矩阵则经过内插变换后的虚拟阵列低阶协方差矩阵为 其中， ()-1表示矩阵的求逆运算； H表示共轭转置； 步骤四： 对虚拟阵列的低阶协方差矩阵进行前后向空间平滑处理， 获得前后向空间平 滑修正后的数据协方差矩阵和噪声协方差矩阵， 通过预白化处理得到虚拟阵列的动态随机 加权低阶协方差矩阵： 将虚。

7、拟变换后阵元数为的虚拟均匀线阵分成相互交错的p个子阵， 每个子阵的阵元数 为即 有对 于 第个 子 阵 ，定 义 两 个维 数 据 矩 阵 ，即 和其中， 0为零矩阵；为维单位矩阵； 是维反对角线为1的置换矩阵；前后向空间平滑修正后的数据协方差 矩阵为前后向空间平滑修正后的噪声协方差矩阵为 其中， 2为噪声方差； *表示共轭； H表示共轭转 置， 对做预白化处理得到虚拟阵列的动态随机加权低阶协方差矩阵 步骤五： 根据最小描述长度准则估计信源个数， 对虚拟阵列的动态随机加权低阶协方 差矩阵进行特征分解， 确定信号子空间和噪声子空间： 根据最小描述长度准则估计信源个数， 对虚拟阵列的动态随机加权低。

8、阶协方差矩阵 进行特征分解:其中， US是由大特征值对应的特征矢量张成的信 号子空间； S是由大特征值构成的对角矩阵； UN是由小特征值对应的特征矢量张成的噪声 子空间； N是由小特征值构成的对角矩阵； 步骤六： 构建动态随机加权低阶协方差-空间平滑-MUSIC测向方法的谱估计公式， 在虚 拟变换的方位角度范围内进行谱峰搜索， 找出极大值点对应的角度， 输出任意阵列相干源 测向结果： 构建动态随机加权低阶协方差-空间平滑-MUSIC测向方法的谱估计公式为 其中，和分别为虚拟变 换方位角度区域的左右边界， 在虚拟变换方位角度范围内根据所构建的谱估计公式进行谱 峰搜索， 找出极大值点对应的角度就是。

9、信号的入射方向， 输出任意阵列相干信源测向结果。 权利要求书 2/2 页 3 CN 110007266 A 3 一种冲击噪声下的任意阵列相干源测向方法 技术领域 0001 本发明提出了一种冲击噪声下的任意阵列的相干源测向方法， 特别是一种冲击噪 声环境下基于动态随机加权低阶协方差的任意阵列相关源测向方法， 该方法涉及阵列信号 处理领域。 背景技术 0002 测向也称为波达方向角(Direction of Arrival,DOA)估计一直是阵列信号处理 领域的热点内容， 在通信、 雷达和声纳等系统中有着广泛应用。 目前大多数测向方法都是基 于等距均匀线阵进行的， 这主要是因为等距均匀线阵具有范德。

10、蒙德(Vandermonde)矩阵形 式， 从而便于数学上的处理。 正是因为等距均匀线阵的特殊阵列结构使得大多数测向方法 只有在这种条件下才能实现， 如前后向空间平滑算法、 矩阵分解等解相干算法等。 0003 然而在许多实际应用环境中， 阵列的结构尺寸往往受到诸多限制， 不可能构造理 想的均匀阵列， 因而研究任意阵列的测向方法具有重要意义。 使用内插变换技术可以实现 对任意阵列的估计， 可以将任意阵列变换成一等距均匀线阵或者是一系列具有相同结构的 阵列。 0004 以MUSIC算法为代表的子空间类高分辨DOA估计算法对于非相干或相关程度较小 的空间信源具有良好的分辨性能， 且运算量较小， 但是。

11、它们优良的分辨性能会随着空间信 源间相关程度的增加而逐渐恶化， 甚至失效。 而在实际的测向环境中， 到达接收阵列的信号 源往往是部分相关或完全相关(即相干)的， 因此如何对空间相干信源去相关， 一直是国内 外学者研究的方向和目标。 0005 根据对已有的技术文献发现， Michael Rubsamen等在 IEEE Transactions on Signal Processing (2009,Vol.57,No.2,pp.588-599)上发表的 “Direction-of-Arrival Estimation for Nonuniform Sensor Arrays:From Manifo。

12、ld Separation to Fourier Domain MUSIC Methods” 提出了基于傅里叶变换思想的root-MUSIC算法， 避免了复杂的谱 峰搜索。 但是这种方法只适用于高斯噪声下任意阵列的独立信号情况， 不能实现相干信号 DOA估计， 且在冲击噪声背景下测向失效。 安春莲等在 哈尔滨工程大学学报 (2013 , Vol.34,No.4,pp.517-523)上发表的 “独立信号与相干信号并存的任意阵列测向方法” 提出 了一种信号子空间测量模型， 并将其用于独立信号与相干信号同时存在的情况， 实现了任 意阵列测向， 但是该方法无法在冲击噪声背景下有效测向。 0006 已。

13、有的文献表明， 使用内插技术可以实现任意阵列到均匀线阵或均匀圆阵的虚拟 变换， 但是至今无有效方法实现冲击噪声背景下的任意阵列相干信源测向， 因此需要设计 一种适用于冲击噪声背景的高性能任意阵列的相干信源测向方法。 本发明设计了一种冲击 噪声环境下基于动态随机加权低阶协方差的任意阵列测向方法， 该方法在复杂的冲击噪声 环境下， 利用内插技术和动态随机加权-空间平滑-MUSIC方法对任意阵列的信源进行有效 估计。 该方法不仅可以实现任意阵列测向， 而且在高斯噪声、 弱冲击噪声和强冲击噪声等其 他恶劣环境下均能获得较好的测向结果。 说明书 1/7 页 4 CN 110007266 A 4 发明内容。

14、 0007 针对上述现有技术， 本发明要解决的技术问题是提供一种可以实现任意阵列测 向， 而且在高斯噪声、 弱冲击噪声和强冲击噪声等其他恶劣环境下均能获得较好的测向结 果的冲击噪声下的任意阵列的相干源测向方法。 0008 为解决上述技术问题， 本发明提供一种冲击噪声下的任意阵列的相干源测向方 法， 包括以下步骤： 0009 步骤一： 建立任意阵列相干信源的采样信号模型： 0010 对于一个阵元数为N的任意阵列， 以原点为参考点， 假设第k个阵元的位置为(xk, yk)， 第i个远场窄带信号从 i方向入射到N元任意阵列， 其中， k1,2,.,N； i1,2,.,M； i表示第i个远场窄带入射信。

15、号的方位角， 则在t时刻第k个阵元接收到的信号为 其中，ki(xkcos i+yksin i)/c； k1,2,.,N； i1,2,.,M； w02 c/ ， c为光速， 为波长； si(t)为t时刻的第i个信源的入射信号， 当入射信号为相干信 号时， 第个生成源生成的相干源组为其中，I0 为生成源总数，为第个生成源生成个信号的权重系数矢量， 则有 nk(t)为t时 刻第k个阵元的噪声信号， nk(t)是满足S S稳定分布的复冲击噪声； k1,2,.,N， 则阵列 接 收 的 第 l 次 快 拍 相 干 信 号 模 型 可 表 示 为式 中 ， A( )a( 1) ,a( 2) ,a( M)N。

16、M为阵列流型矩阵， 其中, j为虚数单位； 1, 2,., M为来波方向方位角矢 量； s(l)为信号矢量；为满足S S稳定分布的复冲击阵列噪声矢 量； T表示转置； 0011 步骤二： 利用阵列接收到的数据构造真实阵列的动态随机加权低阶协方差矩阵： 0012 利 用阵 列接收 的 数 据 构 造 真 实阵 列的 动态随 机 加权 低阶 协方 差 矩阵 式 中 ，其 中， k1,2,.,N， L为最大快拍数； 为动态变化系数； b为加权系数； l为0,1之间均匀分 布的随机数， H表示共轭转置； 0013 步骤三： 定义内插变换矩阵T， 构建阵列信号内插变换后的虚拟阵列低阶协方差矩 阵： 00。

17、14首先对虚拟变换的方位角度区域进行划分， 假设信号位于区域内， 将区域均匀 划分为和分别为虚拟变换方位角度区域的左右 说明书 2/7 页 5 CN 110007266 A 5 边界 ，为虚拟变换方位角度变化步长 ， 则在区域内的真实阵列流行矩阵为 为虚拟阵列的阵元个数， 则在 同一区域内， 假设虚拟阵列为阵元间距相等的均匀线阵， 虚拟阵列的阵列流型矩阵为 其中， 根据真实阵列流型矩阵与虚拟阵列流型矩阵定义 变换矩阵则经过内插变换后的虚拟阵列低阶协方差矩阵为 其中， ()-1表示矩阵的求逆运算； H表示共轭转置； 0015 步骤四： 对虚拟阵列的低阶协方差矩阵进行前后向空间平滑处理， 获得前后。

18、向空 间平滑修正后的数据协方差矩阵和噪声协方差矩阵， 通过预白化处理得到虚拟阵列的动态 随机加权低阶协方差矩阵： 0016将虚拟变换后阵元数为的虚拟均匀线阵分成相互交错的p个子阵， 每个子阵的 阵 元 数 为即 有对 于 第个 子阵 ， 定 义 两 个维 数 据 矩阵 ， 即 和其中， 0为零矩阵；为维单位矩阵； 是维反对角线为1的置换矩阵；前后向空间平滑修正后的数据协方差 矩阵为前后向空间平滑修正后的噪声协方差矩阵为 其中， 2为噪声方差； *表示共轭； H表示共轭转 置， 对做预白化处理得到虚拟阵列的动态随机加权低阶协方差矩阵 0017 步骤五： 根据最小描述长度准则估计信源个数， 对虚拟。

19、阵列的动态随机加权低阶 协方差矩阵进行特征分解， 确定信号子空间和噪声子空间： 0018 根据最小描述长度准则估计信源个数， 对虚拟阵列的动态随机加权低阶协方差矩 阵进行特征分解:其中， US是由大特征值对应的特征矢量张成 的信号子空间； S是由大特征值构成的对角矩阵； UN是由小特征值对应的特征矢量张成的 噪声子空间； N是由小特征值构成的对角矩阵； 0019 步骤六： 构建动态随机加权低阶协方差-空间平滑-MUSIC测向方法的谱估计公式， 在虚拟变换的方位角度范围内进行谱峰搜索， 找出极大值点对应的角度， 输出任意阵列相 干源测向结果： 0020 构建动态随机加权低阶协方差-空间平滑-MU。

20、SIC测向方法的谱估计公式为 其中，和分别为虚拟 说明书 3/7 页 6 CN 110007266 A 6 变换方位角度区域的左右边界， 在虚拟变换方位角度范围内根据所构建的谱估计公式进行 谱峰搜索， 找出极大值点对应的角度就是信号的入射方向， 输出任意阵列相干信源测向结 果。 0021 本发明的有益效果： 与现有技术相比， 本发明具有以下优点： 0022 (1)本发明在高斯噪声、 弱冲击噪声和强冲击噪声等复杂环境下对任意阵列接收 的快拍数据进行处理， 在任意阵列结构下均能对来波方向进行有效估计。 0023 (2)本发明对某个虚拟变换的方位角度范围进行划分， 求出虚拟变换角度范围内 真实阵列的。

21、动态随机加权协方差矩阵和虚拟变换后的虚拟动态随机加权协方差矩阵， 实现 任意阵列到均匀线阵的虚拟变换， 再结合空间平滑MUSIC测向方法， 能在冲击噪声等复杂环 境下， 对任意阵列的相干源、 独立源、 相关源与独立源同时存在的混合源进行有效估计， 具 有较好的解相干能力， 测向性能优， 应用范围广泛。 附图说明 0024 图1为基于动态随机加权协方差的空间平滑MUSIC任意阵列相干源测向方法基本 框图； 0025 图2为特征指数为0.8时， DSW-SS-MUSIC和SS-MUSIC波达方向估计仿真对比曲线； 0026 图3为特征指数为1.3时， DSW-SS-MUSIC和SS-MUSIC波达。

22、方向估计仿真对比曲线； 0027 图4为特征指数为1.8时， DSW-SS-MUSIC和SS-MUSIC波达方向估计仿真对比曲线； 0028 图5为特征指数为2.0即高斯噪声下， DSW-SS-MUSIC和SS-MUSIC波达方向估计仿真 对比曲线。 具体实施方式 0029 下面结合附图对本发明具体实施方式进行进一步说明。 0030 本发明通过使用内插预处理技术和构建动态随机加权低阶协方差矩阵， 设计了一 种适用于冲击噪声等复杂环境的任意阵列测向方法。 实现步骤如下： 建立任意阵列相干信 源的采样信号模型； 利用阵列接收到的数据构造真实阵列的动态随机加权低阶协方差矩 阵； 定义内插变换矩阵T，。

23、 构建阵列信号内插变换后的虚拟阵列协方差矩阵； 对虚拟阵列的 协方差矩阵进行前后向空间平滑处理， 获得前后向空间平滑修正后的数据协方差矩阵和噪 声协方差矩阵， 通过预白化处理得到虚拟阵列的动态随机加权低阶协方差矩阵； 根据最小 描述长度准则估计信源个数， 对虚拟阵列的动态随机加权协方差矩阵进行特征分解， 确定 信号子空间和噪声子空间； 构建动态随机加权低阶协方差-空间平滑-MUSIC测向方法的谱 估计公式， 在虚拟变换的方位角度范围内进行谱峰搜索， 找出极大值点对应的角度， 输出任 意阵列相干源测向结果。 本发明能够对任意阵列的信源来波方向进行有效估计， 可在高斯 噪声、 弱冲击噪声和强冲击噪。

24、声下进行测向， 解相干性能优， 应用范围广泛。 0031 如图1所示， 步骤一： 建立任意阵列相干信源的采样信号模型。 0032 对于一个阵元数为N的任意阵列， 以原点为参考点， 假设第k个阵元的位置为(xk, yk)， 第i个远场窄带信号从 i方向入射到N元任意阵列， 其中， k1,2,.,N； i1,2,.,M； i表示第i个远场窄带入射信号的方位角。 则在t时刻第k个阵元接收到的信号为 说明书 4/7 页 7 CN 110007266 A 7 其中，ki(xkcos i+yksin i)/c； k1,2,.,N； i1,2,.,M； w02 c/ ， c为光速， 为波长； si(t)为t。

25、时刻的i第个信源的入射信号， 当入射信号为相干信 号时， 第个生成源生成的相干源组为其中，I0 为生成源总数，为第个生成源生成个信号的权重系数矢量， 则有 nk(t)为t时 刻第k个阵元的噪声信号， nk(t)是满足S S稳定分布的复冲击噪声； k1,2,.,N。 则阵列 接 收 的 第 l 次 快 拍 相 干 信 号 模 型 可 表 示 为式 中 ， A( )a( 1) ,a( 2) ,a( M)NM为阵列流型矩阵， 其中, j为虚数单位； 1, 2,., M为来波方向方位角矢量； s (l)为信号矢量；为满足S S稳定分布的复冲击阵列噪声矢量； T 表示转置。 0033 步骤二： 利用阵列。

26、接收到的数据构造真实阵列的动态随机加权低阶协方差矩阵。 0034 利 用阵 列接收 的 数 据 构 造 真 实阵 列的 动态随 机 加权 低阶 协方 差 矩阵 式中，其中， k 1,2,.,N， L为最大快拍数； 为动态变化系数； b为加权系数； l为0,1之间均匀分布的随 机数， H表示共轭转置。 0035 步骤三： 定义内插变换矩阵T， 构建阵列信号内插变换后的虚拟阵列低阶协方差矩 阵。 0036首先对虚拟变换的方位角度区域进行划分， 假设信号位于区域内， 将区域均匀 划分为和分别为虚拟变换方位角度区域的左右 边界 ，为虚拟变换方位角度变化步长 ， 则在区域内的真实阵列流行矩阵为 为虚拟阵。

27、列的阵元个数。 则在 同一区域内， 假设虚拟阵列为阵元间距相等的均匀线阵， 虚拟阵列的阵列流型矩阵为 其中， 根据真实阵列流型矩阵与虚拟阵列流型矩阵定义 变换矩阵则经过内插变换后的虚拟阵列低阶协方差矩阵为 说明书 5/7 页 8 CN 110007266 A 8 其中， ()-1表示矩阵的求逆运算； H表示共轭转置。 0037 步骤四： 对虚拟阵列的低阶协方差矩阵进行前后向空间平滑处理， 获得前后向空 间平滑修正后的数据协方差矩阵和噪声协方差矩阵， 通过预白化处理得到虚拟阵列的动态 随机加权低阶协方差矩阵。 0038将虚拟变换后阵元数为的虚拟均匀线阵分成相互交错的p个子阵， 每个子阵的 阵 元。

28、 数 为即 有对 于 第个 子阵 ， 定 义 两 个维 数 据 矩阵 ， 即 和其中， 0为零矩阵；为维单位矩阵； 是维反对角线为1的置换矩阵；前后向空间平滑修正后的数据协方差 矩阵为前后向空间平滑修正后的噪声协方差矩阵为 其中， 2为噪声方差； *表示共轭； H表示共轭转 置。 对做预白化处理得到虚拟阵列的动态随机加权低阶协方差矩阵 0039 步骤五： 根据最小描述长度准则估计信源个数， 对虚拟阵列的动态随机加权低阶 协方差矩阵进行特征分解， 确定信号子空间和噪声子空间。 0040 根据最小描述长度准则估计信源个数， 对虚拟阵列的动态随机加权低阶协方差矩 阵进行特征分解:其中， US是由大特。

29、征值对应的特征矢量张成 的信号子空间； S是由大特征值构成的对角矩阵； UN是由小特征值对应的特征矢量张成的 噪声子空间； N是由小特征值构成的对角矩阵。 0041 步骤六： 构建动态随机加权低阶协方差-空间平滑-MUSIC测向方法的谱估计公式， 在虚拟变换的方位角度范围内进行谱峰搜索， 找出极大值点对应的角度， 输出任意阵列相 干源测向结果。 0042 构建动态随机加权低阶协方差-空间平滑-MUSIC测向方法的谱估计公式为 其中，和分别为虚拟 变换方位角度区域的左右边界。 在虚拟变换方位角度范围内根据所构建的谱估计公式进行 谱峰搜索， 找出极大值点对应的角度就是信号的入射方向。 输出任意阵列。

30、相干信源测向结 果。 0043 为便于叙述， 将基于动态随机加权协方差的空间平滑MUSIC任意阵列相干源测向 方法简记为DSW-SS-MUSIC， 将用于比较的空间平滑MUSIC任意阵列相干源测向方法记为SS- MUSIC。 0044 冲击噪声下基于动态速随机加权协方差的任意阵列相干源测向方法的参数设置 情况如下： 阵元数M5， 方位角变化步长最大快拍数L2000， 动态变化系数 0.99， 加权系数b0.02， 虚拟阵列的阵元数子阵数p7， 每个子阵的阵元数 0045 SS-MUSIC任意阵列测向方法的参数设置情况参考刁鸣等在 弹舰与制导学报 说明书 6/7 页 9 CN 110007266。

31、 A 9 (Vol.25,No.4,pp:416-418)发表的 “基于虚拟阵列变换和修正MUSIC算法的相干源测向” 。 0046 图2中， 冲击噪声特征指数为0.8， 两个相干信源从 6 ,24 方向入射， 虚拟变 换的角度区域为0 30 ， 阵元位置坐标为0.5 (0,0),0.5 (4,3),0.5 (5,5),0.5 (7,7), 0.5 (9,10)， 广义信噪比为30dB。 从仿真图2中可以看出在强冲击噪声下， 本发明所设计的 DSW-SS-MUSIC测向方法能够估计出信源个数和信源的来波方向， 而SS-MUSIC测向方法无法 估计信源数目和信源来波方向。 0047 图3中， 冲。

32、击噪声特征指数为1.3， 两个相干信源从 6 ,24 方向入射， 虚拟变 换的角度区域为0 30 ， 阵元位置坐标为0.5 (0,0),0.5 (4,3),0.5 (5,5),0.5 (7,7), 0.5 (9,10)， 广义信噪比为30dB。 从图3中可以看出， 在冲击噪声指数为1.3的弱冲击噪声 下， DSW-SS-MUSIC测向方法能够准确估计出信源个数和信源的来波方向， 而SS-MUSIC测向 方法对信源数目估计不准确， 无法有效估计信源的来波方向。 0048 图4中， 冲击噪声特征指数为1.8， 两个相干信源从 6 ,24 方向入射， 虚拟变 换的角度区域为0 30 ， 阵元位置坐标。

33、为0.5 (0,0),0.5 (4,3),0.5 (5,5),0.5 (7,7), 0.5 (9,10)， 广义信噪比为20dB。 从图4中可以看出， 在冲击噪声指数为1.8的弱冲击噪声 下， DSW-SS-MUSIC测向方法能够准确估计出信源个数和信源的来波方向， 而SS-MUSIC测向 方法虽然能够估计出信源数目但是无法准确估计出信源的来波方向， 但是测向精度不高。 0049 图5中， 在高斯白噪声下， 两个相干信源从 6 ,24 方向入射， 虚拟变换的角度 区域为0 30 ， 阵元位置坐标为0.5 (0,0),0.5 (4,3),0.5 (5,5),0.5 (7,7),0.5 (9, 1。

34、0)， 广义信噪比为20dB。 从图5中可以看出， 在高斯噪声下DSW-SS-MUSIC测向方法和SS- MUSIC测向方法均能估计出信源数目和信源来波方向， 但是DSW-SS-MUSIC测向方法对信源 来波方向的估计更加准确。 说明书 7/7 页 10 CN 110007266 A 10 图1 说明书附图 1/5 页 11 CN 110007266 A 11 图2 说明书附图 2/5 页 12 CN 110007266 A 12 图3 说明书附图 3/5 页 13 CN 110007266 A 13 图4 说明书附图 4/5 页 14 CN 110007266 A 14 图5 说明书附图 5/5 页 15 CN 110007266 A 15 。