圆阵自适应干扰抑制仿真系统

资 源 简 介

本系统专门用于研究和演示均匀圆阵（UCA）在复杂电磁环境下的自适应干抑制性能。系统首先构建圆阵的数学物理模型，计算阵列流形矢量，并结合空间信号源的方位角和俯仰角信息，模拟期望信号与多路干扰信号的接收过程。核心功能采用自适应波束形成算法（如最小方差无失真响应MVDR或LCMV准则）来计算最优权矢量，在保证期望方向增益的同时，在干扰方向自动形成深零陷。系统提供了极具深度的可视化功能，能够生成自适应方向图的精确数据，并从多个角度剖析阵列性能，包括在特定俯仰角或方位角下的二维平面投影图。此外，系统支持在两种坐标系下导出三维图形：一种是基于直角坐标系的增益分布图，用于直观观察旁瓣级和零陷深度；另一种是基于球坐标系的三维极坐标图，用于真实还原阵列在空间的物理辐射特性。该项目适用于雷达抗干扰、车载通信天线设计及声呐阵列处理等领域，能够有效提升复杂环境下的信号干扰噪声比（SINR）。

详 情 说 明

圆阵列自适应干扰抑制仿真系统

项目介绍

本系统是一个基于 MATLAB 开发的均匀圆阵（UCA）自适应干扰抑制仿真平台。其核心目标是模拟在复杂电磁环境下，通过构建圆阵的数学物理模型，利用自适应波束形成算法实现对期望信号的增强以及对多路强干扰信号的有效抑制。系统能够直观展示阵列流形对空间增益的影响，并通多维度的可视化图形分析系统性能，包括零陷深度、旁瓣电平以及物理空间的增益分布。

---

功能特性

1. 精确的阵列流形建模：针对均匀圆阵（UCA）建立物理模型，根据阵元数量和半波长布阵原则计算阵列半径，并实现方位角与俯仰角耦合的空间流形矢量计算。
2. 多路干扰抑制能力：支持设置多个干扰源的方位、俯仰及干扰噪声比（INR），模拟真实场景下的复杂干扰环境。
3. 自适应波束形成算法：采用最小方差无失真响应（MVDR）准则计算最优权矢量，在保证期望方向归一化增益的同时，在干扰方向自动形成深零陷。
4. 全空间增益扫描：系统对全方位（-180°至180°）和俯仰（0°至90°）空间进行网格化扫描，生成精细的增益分布矩阵。
5. 多维度性能可视化：

* 二维剖面图：提取期望方向的方位向和俯仰向切片，观察波束宽度与干扰处的下陷。 * 直角坐标系三维方向图：以方位和俯仰为底面，增益（dB）为高度，观察旁瓣分布。 * 球坐标系物理图：将增益映射为极径，还原阵列在真实三维物理空间的辐射特性，并标注信号与干扰的射线方向。

---

系统要求

• 开发环境：MATLAB R2016b 或更高版本。
• 必备工具箱：无（仅基于标准矩阵运算与绘图函数实现）。

---

实现逻辑与算法细节

#### 1. 阵列几何与流形构建 系统的基础是 $N=16$ 阵元的均匀圆阵。阵元等间隔分布在 X-Y 平面的圆周上。阵列流形矢量通过以下闭包函数实现： $A(theta, phi) = exp(j cdot frac{2pi R}{lambda} cdot cosphi cdot cos(theta - alpha_n))$ 其中，$alpha_n$ 是第 $n$ 个阵元的角位置，$theta$ 和 $phi$ 分别表示方位角和俯仰角。

#### 2. 信号环境模拟 系统通过统计特性直接构造接收信号的协方差矩阵 $R_{xx}$。总协方差矩阵由期望信号（$R_s$）、干扰信号（$R_i$）和加性高斯白噪声（$R_n$）组成： $R_{xx} = P_s a_s a_s^H + sum P_{im} a_{im} a_{im}^H + P_n I$ 这种建模方式避免了时域仿真的随机扰动，能够直接反映算法在理论上的最优性能。

#### 3. 最优权矢量计算 系统基于 MVDR（最小方差无失真响应）准则求解权矢量 $w$。该准则的优化目标是在约束期望方向增益为 1 的条件下最小化输出总功率，其闭合解为： $w = frac{R_{xx}^{-1} a_s}{a_s^H R_{xx}^{-1} a_s}$ 通过矩阵求逆和内积运算，得到对各阵元接收信号的复加权系数。

#### 4. 可视化处理细节

• 归一化增益计算：全空间扫描得到的功率数据经过归一化处理，并转化为分贝（dB）单位。为了视觉清晰度，系统将动态范围限制在 -60dB 到 0dB。
• 物理映射：在三维极坐标图中，系统利用球坐标系转换公式，将处理后的增益值映射为坐标点到原点的距离 $R_{text{disp}}$，从而生动展现出“波束”和“零陷”的空间形态。

---

使用方法

1. 参数配置：修改脚本开头的配置区，可自定义阵元数 $N$、期望信号位置 (theta_s, phi_s)、干扰位置 (theta_i, phi_i) 以及信噪比 SNR 和干噪比 INR。
2. 执行计算：运行脚本后，系统将自动完成协方差矩阵构造、权矢量求解及全空间 Gain 矩阵的计算。
3. 结果解读：

* 查看命令行输出的最优权矢量复数值。 * 分析图形窗口中的四幅图表：通过二维剖面图确认干扰方向的零陷深度；通过三维图确认主瓣是否指向期望方向，以及旁瓣级是否满足设计要求。