CBF & DOA估计

llnheu

于 2024-03-20 23:44:25 发布

阅读量925

点赞数 15

分类专栏： DOA学习文章标签：算法 matlab

本文链接：https://blog.csdn.net/weixin_63880614/article/details/136890402

版权

DOA学习专栏收录该内容

1 篇文章 0 订阅

订阅专栏

DOA 估计算法学习（一）

前言

这学期要完成一下水下目标探测与跟踪的科创项目，项目涉及到的一个重要方法就是DOA估计，想在学习的同时通过博客进行记录，方便总结回顾，首先就是利用CBF即波束形成进行DOA估计【20240320】。

一、CBF简介

CBF（Conventional Beamforming）即所谓常规波束形成，其原理是由于多个阵元接收到信号具有相位差，通过在接收端进行相位补偿求和，来确认波达方向的，波束形成器结构如下：

式中的W指的是权重，根据加权方式的不同W的值也不同，但常用的都是均匀加权。

对于波束形成器的理解如下图：

这里只选出了两个阵元来说明波束形成的原理，可以看到在接收信号的时候，由于角度的关系，1号阵元和2号阵元接收到的信号是有时延的，这个时延等于 $d*sin(\theta )/c$ ，对应到频域就是相位差。

如果直接将接收到的信号相加，由于相位差的原因信号一定是有衰减的，为了消除信号之间的相位差，需要对接收到的信号进行相位补偿。

但是根据时延可以看出，要想准确地补偿相位需要知道波达方向θ，而在实际应用中是不知道θ的（如果知道就不需要进行DOA估计了），所以进行枚举，从-90°试到90°，都去补偿信号，补偿完把每个阵元得到的信号相加，求总能量。

如果是正确的角度，信号能量就会达到最大，因为全部信号都没有相位差，相加没有能量损失；如果是其他角度，最后得到的信号能量值就有限。因此，假设以角度-90到90作为x轴，以信号能量作为y轴，就能得到一张清晰展现波达方向的图像，也就是波束图。

二、Matlab代码

1. 代码

%% 20240320 CBF进行DOA估计
% 以入射角为30°的信号为例，进行CBF的仿真实验
clc;
close all;
%% 构筑回波
f0 = 75e3;      % 信号频率
fs = 3*f0;      % 采样频率
T  = 0.01;      % 脉宽
N  = 1000;      % 快拍数
c  = 1500;      % 水中声速
tsig = 0:1/fs:(N-1)/fs;
sig = exp(j*2*pi*f0*tsig);
theta = 30;     % DOA角度

%% 接收阵列参数
d = c/f0/2;     % 这是阵元间距的基本要求，目的是为了避免栅瓣
M = 18;          % 阵元数，均匀线列阵
X = exp((0:M-1)'*-j*2*pi*f0*d*sin(theta/180*pi)/c)*sig; % 阵列接收到的信号，是M*N的矩阵，每个阵元接收到的信号相差固定相位
snr = 10;       % SNR信噪比设为10
X = awgn(X,snr); 
%% CBF
THETA = linspace(-90,90,181);
for i = 1:length(THETA)
    a = exp(-j*2*pi*f0*d*sin(THETA(i)/180*pi)/c*(0:M-1)'); %相位补偿向量
    Y = a'* X;      %对接收到的信号做相位补偿并求和，Y是1*N的向量
    b_CBF(i) = Y*Y'/N; %求平均能量
end
% [b I] = max(abs(b_CBF));
B_CBF = 10*log10(abs(b_CBF)/max(abs(b_CBF)));

plot(THETA,B_CBF);
title('波束图');
xlabel('角度（°）');
ylabel('总能量');