24阵元水声相移波束形成仿真系统

资 源 简 介

本项目实现了一个针对水声信号处理应用的高精度相移波束形成仿真程序。该系统模拟由24个等间隔阵元组成的接收阵列，利用窄带相移原理对接收到的声学信号进行空间滤波。通过计算阵流矢量与权重矢量的乘积，实现在-90度到90度范围内生成181个扫描波束，每个波束的角分辨率为1度。程序核心实现了从时域信号到频域相移的转换过程，对各阵元接收到的微弱水声信号施加特定的相位延迟，使得目标方向的信号能在相干叠加后获得最大增益。该仿真不仅涵盖了基础的波束合成逻辑，还包含了阵列响应函数的推导以及空间指向性分布的计算，能够有效评估阵

详 情 说 明

相移波束形成仿真系统 (24阵元)

项目介绍

本仿真系统是一个基于MATLAB开发的24阵元等间隔线性阵列（ULA）相移波束形成程序。该项目模拟了水声环境中的信号接收与处理过程，通过窄带相移原理实现在特定空间方向上的增益极大化。系统能够精细地模拟声波到达阵列时的空间相位差，并通过权重补偿和信号叠加，实现从各阵元接收到的混杂噪声信号中提取特定方向目标信号的目的。该系统适用于水下声呐阵列指向性分析、目标方位估计（DOA）及空间增益评估等研究领域。

功能特性

• 高阶采样与窄带仿真：采用50kHz采样率对5kHz中心频率的声学信号进行处理，严格遵循半波长阵元间距（d = λ/2）准则。
• 动态波束扫描：在-90度到90度的空间范围内进行全覆盖扫描，生成181个扫描波束，具备1度的角分辨率。
• 信噪比受控模拟：支持自定义SNR（信噪比）参数，生成复数高斯白噪声以模拟真实的水平复杂声场环境。
• 自动性能评估：程序可自动计算并输出波束指向误差、主瓣-3dB带宽以及第一旁瓣级（SLL）等关键技术指标。
• 双模可视化展示：同时提供直角坐标系下的功率谱分布图和符合声呐显示习惯的极坐标空间指向性图。

系统要求

• 软件环境：MATLAB R2016b 或更高版本。
• 工具箱要求：Signal Processing Toolbox（用于findpeaks函数进行旁瓣分析）。
• 硬件要求：标准桌面或笔记本电脑，无需特殊计算加速器。

实现逻辑与功能说明

main.m 脚本按照标准的信号处理流程分为以下五个核心步骤：

1. 参数初始化 设置仿真基础环境，包括定义24个阵元、确定1500m/s的声速、计算5kHz主频对应的波长及其半波长间距。同时预设目标入射角（如15度）和信噪比（如15dB）。

2. 目标信号与噪声建模

• 复解析信号生成：产生复指数形式的窄带信号。
• 阵列流形构建：根据目标入射角度，计算各阵元相对于参考阵元的空间时延，进而导出相位延迟，形成阵列流形矢量（a_target）。
• 空间信号矩阵：将单信道信号扩展至24阵元，并叠加复高斯随机噪声，模拟空间接收矩阵。

3. 相移波束形成核心算法

• 权重矢量计算：针对每一个扫描角度（-90°:1°:90°），计算对应的相移补偿权重矢量。
• 空间相干叠加：执行权重矢量的共轭转置与接收信号矩阵的乘法，实现空间滤波。该操作使得目标方向的信号同相叠加，而其他方向信号相互抵消。
• 能量统计：对每个波束输出的时间序列计算均方根功率。

4. 技术指标自动分析

• 归一化处理：将功率转换为对数刻度（dB），并以最大值进行归一化。
• 主瓣分析：通过寻找峰值左右两侧降低3dB的交叉点，计算主瓣的宽度。
• 旁瓣检测：利用峰值搜索算法识别除了主波束之外的最大响应点，计算其与主瓣的电平差值。

5. 结果可视化输出

• 直角坐标图：展示-90度到90度范围内的归一化功率谱曲线。
• 极坐标图：将波束响应映射至极坐标空间，直观展现阵列的空间指向性能。

关键技术与算法细节

• 相移原理：系统利用窄带信号的特性，将空间几何距离引起的“时延”转化为频域中的“相移”，通过 $w = exp(-j cdot 2pi cdot f_0 cdot tau)$ 进行权重补偿。
• 阵列流形（Array Manifold）：代码通过 $dist cdot sin(theta) / c$ 建立空间几何模型，准确描述了波前到达线性阵列的物理过程。
• 信噪比控制：通过计算噪声功率对复高斯随机数进行缩放，确保仿真环境在不同干扰强度下的鲁棒性验证。
• 离散波束合成：采用循环迭代方式遍历空间角度，这种实现方式清晰地展示了常规束形成器（CBF）的工作机制，利于算法理解和参数调试。
• 辅助报告系统：通过 num_warm 等逻辑处理字符转换，并在控制台即时打印包含指向精度、带宽和旁瓣级的技术报告，实现了仿真与评估的一体化。