【MATLAB源码-第186期】matlab基于MLE算法的8天线阵列DOA估计仿真，对比粗估计、精确估计输出RMSE对比图。

操作环境：

MATLAB 2022a

1、算法描述

第一部分：基本概念与系统设置

方向到达估计（Direction of Arrival, DOA）是信号处理中一项重要的技术，主要用于确定信号的到达方向。这种技术在雷达、无线通信和声纳等领域中有广泛的应用。DOA估计的核心目的是从接收到的信号中提取出信号源的位置信息。

在基于网格搜索的最大似然估计（Maximum Likelihood Estimation, MLE）方法中，首先需要构建一个天线阵列系统。这个系统通常由多个均匀排列的天线组成，每个天线都能接收来自同一信号源的信号。由于信号在不同天线之间传播的距离存在差异，因此会产生相位延迟，这些相位延迟包含了关于信号源方向的信息。

在实际操作中，首先需要定义一个搜索范围，这个范围内包含了可能的到达角度。然后，系统将这个范围划分为多个小的网格。每一个网格点都代表一个潜在的到达角度，系统将依次评估这些角度。

第二部分：网格搜索与成本函数

在基于网格搜索的最大似然估计（MLE）方法中，系统的核心是一个成本函数，这个函数能够衡量某个特定角度假设下的似然性有多高。对于每一个角度网格点，系统会计算一个所谓的“方向矢量”，这个矢量是一个复数数组，描述了在该角度假设下，信号在每个天线上的预期相位。

粗估计

在初步的DOA估计过程中，首先进行的是粗略估计。在这一阶段，搜索网格的尺寸较大，例如每个网格可能代表一度或更多。系统将逐一评估这些网格点，计算对应的成本函数。这个成本函数通常是基于预测信号和实际接收信号之间的误差来定义的，误差越小，成本函数值越小，表示该角度的似然性越高。选择成本函数值最小的角度作为粗略的DOA估计结果。这一步虽然不够精确，但可以快速缩小潜在的目标角度范围，为后续的精细估计奠定基础。

精细估计

粗略估计之后，系统进入精细估计阶段。在这一阶段，搜索的网格尺寸将显著减小，例如每个网格点可能只代表0.1度或更少。系统在粗估计得到的角度附近进行更密集的搜索，再次计算每个网格点的成本函数。由于网格更细，这一步能够提供更高精度的角度估计。精细估计虽然计算量更大，但通过减小网格尺寸，能够更准确地定位信号源的真实方向。

Cramer-Rao 下边界 (CRLB)

除了实际的估计过程，理论分析中常用的一个重要工具是Cramer-Rao下边界（CRLB）。CRLB为估计器的性能提供了一个理论下限，即在给定的信号条件下，任何无偏估计器的方差不可能低于这个界限。CRLB的计算通常依赖于信号模型的细节，包括信号和噪声的统计特性。在DOA估计中，CRLB可以帮助我们评估估计方法的效率，并对比不同方法在相同条件下的性能。