1. CBF波束形成与DOA估计基础
第一次接触DOA估计时，我被各种算法弄得晕头转向，直到遇到CBF（常规波束形成）这个"老实人"。它不像MUSIC、ESPRIT那些高深算法，而是用最直观的物理原理解决问题——就像用耳朵判断声音方向一样自然。
想象你在操场上闭眼听同学喊你。如果声音从正前方来，两耳同时听到；如果从侧面来，靠近声源的那只耳朵会先听到。CBF就是这个原理的数学版：通过计算阵列天线各阵元接收信号的相位差，反推信号来源方向（DOA）。实际项目中，我用它做过声呐目标定位，效果稳定得像老式收音机调频——虽然分辨率不高，但绝不会突然罢工。
核心公式就藏在时延差里：Δτ = d·sinθ/c。其中d是阵元间距，θ是入射角，c是波速。这个时延转换成相位差就是φ=2πfΔτ。在MATLAB里，我们通过构造导向矢量（也叫权重向量）来补偿这个相位差：
a = exp(-1j*2*pi*f0*d*sin(theta)/c*(0:M-1)');这个看似简单的复指数运算，其实完成了三件事：1)计算理论时延 2)转为相位差 3)生成补偿系数。当补偿角度θ与实际DOA匹配时，各阵元信号会同相叠加，能量出现峰值——这就是CBF的"智能"之处。
2. MATLAB仿真环境搭建
建议直接用R2020b以上版本，信号处理工具箱会省去很多麻烦。我曾在旧版遇到过awgn函数不兼容的问题，调试到怀疑人生。新建脚本时务必先执行这三连：
clc; clear; close all;阵列参数设置是第一个坑。根据空间采样定理，阵元间距d≤λ/2（λ是波长），否则会出现栅瓣——就像数码相机拍条纹衬衫时产生的莫尔条纹。但间距也不是越小越好，我曾把d设为λ/4，结果波束宽度直接翻倍。经验公式：
f0 = 75e3; % 信号频率 c = 1500; % 声速(水下场景) d = c/f0/2; % 最优间距信号生成部分要注意采样率fs≥3f0（香农定理底线）。但实际中我习惯用5倍：
fs = 5*f0; T = 0.01; % 脉宽 N = T*fs; % 采样点数 tsig = (0:N-1)/fs; sig = exp(1j*2*pi*f0*tsig); % 复信号添加噪声时推荐awgn函数，它的信噪比参数是相对于信号功率的。有个冷知识：用'measured'选项可以自动计算输入信号功率，避免手动换算：
X_noisy = awgn(X, snr, 'measured');3. 阵列信号建模关键步骤
建模仿真时，新手常犯的错误是直接复制理论公式。比如接收信号建模，正确的矩阵运算应该是：
X = exp((0:M-1)'*-1j*2*pi*f0*d*sin(theta_doa/180*pi)/c) * sig;这里有几个细节：
(0:M-1)'生成列向量，实现阵元间相位递推- 角度θ必须转为弧度制，我见过有人漏掉
pi/180导致峰值出现在0.5度位置 - 矩阵乘法
*会自动扩展sig向量，等效于各阵元接收相同信号但有相位差
有一次我误用了.（点乘）代替，结果得到一堆杂乱频谱。后来才明白点乘会导致元素逐一相乘，破坏了阵列的相位关系。
导向矢量的生成也有讲究。扫描角度范围建议用linspace而非-90:1:90，前者能精确控制点数：
theta_scan = linspace(-90, 90, 181); % 1度间隔 a = exp(-1j*2*pi*f0*d*sind(theta_scan).*(0:M-1)'/c);注意这里用sind替代sin，避免频繁的弧度转换。.*实现矩阵外积，一次性生成所有扫描角度的导向矢量。
4. 波束形成与能量谱计算
核心循环其实可以向量化运算，提升10倍速度。这是我优化后的写法：
Y = a' * X; % 矩阵乘法一次性完成所有角度补偿 P = sum(abs(Y).^2, 2)/N; % 按行求和得功率曾经我在这里踩过坑：忘记对补偿后的信号取模平方，直接用了Y*Y'，导致出现负功率值。正确的能量计算应该是|Y|²。
波束图绘制时，归一化处理很重要：
P_dB = 10*log10(P/max(P)); % 转为dB尺度 plot(theta_scan, P_dB); xlabel('角度(°)'); ylabel('归一化功率(dB)');如果发现主瓣宽度异常，检查两个参数：1)阵元数M 2)间距d与波长比。我曾用16阵元仿真时，波束宽度比理论值宽了15%，最后发现是d设置成了0.6λ。
5. 参数影响与性能分析
通过修改M值对比效果最直观。下面是我实验的数据记录：
| 阵元数M | 理论波束宽度 | 实测3dB宽度 | 旁瓣电平 |
|---|---|---|---|
| 4 | 28.6° | 29.2° | -11.3dB |
| 8 | 14.3° | 14.8° | -13.2dB |
| 16 | 7.1° | 7.3° | -13.8dB |
可见阵元数加倍，分辨率近似线性提升。但硬件成本也倍增，工程中需要权衡。
信号频率的影响更微妙。保持d=λ/2时，提高频率会使波束变窄：
f0 = [50e3 75e3 100e3]; beamwidth = 50.7 ./ (M*d*f0/c); % 度数但频率过高会导致传播损耗增大，实际项目中选择75kHz是个折中方案。
6. 工程实践中的注意事项
硬件部署时，阵元位置误差会严重影响性能。有次实测发现峰值偏移5度，最后查出是第三个水听器安装偏差2mm。建议在仿真中加入位置扰动：
pos_err = 0.01*d*randn(M,1); % 1%间距的随机误差 real_pos = (0:M-1)'*d + pos_err;多目标场景下，CBF可能出现峰值合并现象。解决方法之一是先粗扫再精扫：
% 第一阶段：5度间隔粗扫 theta_coarse = -90:5:90; % 第二阶段：峰值附近1度间隔 theta_fine = max_angle-10:1:max_angle+10;最后提醒：所有角度扫描算法都受限于瑞利限，两个目标的角度差需大于波束宽度才能分辨。这是物理规律，算法再优化也无法突破。
