线阵单脉冲和差波束测角与误差分析仿真

资 源 简 介

本项目实现了一个基于8阵元均匀线阵的单脉冲和差波束测角MATLAB程序。该程序的核心功能是利用全阵列同时形成两个波束——和波束（Sum Beam）与差波束（Difference Beam），通过比较和差通道的输出来解算目标的方位角。程序详细模拟了相控阵雷达的测角过程，重点在于分析测角精度与扫描角度的关系。仿真结果展示了一个关键特性：随着测量角度偏离阵列法线方向（角度增大），测角误差会逐渐增大。程序对造成这一误差的原因进行了深入分析，主要归因于单脉冲测角算法在大角度下的近似误差，以及波束在进行宽角扫描时主瓣宽度变宽（即“波束变胖”效应），导致方向性增益下降和鉴角线性度变差。该代码适用于雷达原理教学及相控阵测角算法的性能评估。

详 情 说 明

线阵单脉冲和差波束测角仿真系统

项目简介

这是一个基于MATLAB开发的相控阵雷达仿真项目，专注于演示和分析8阵元均匀线阵（ULA）的单脉冲和差波束测角技术。该系统不仅模拟了静态（法线方向）的测角过程，还重点实现了宽角扫描条件下的性能评估。通过蒙特卡洛仿真，程序深入揭示了波束扫描角度对波束宽度（"波束变胖"效应）以及最终测角精度（RMSE）的影响，直观展示了相控阵雷达在偏离法线方向时性能下降的物理规律。

该代码非常适合用于雷达原理教学、阵列信号处理算法验证以及相控阵测角系统的性能预估。

功能特性

• 物理参数建模：仿真了工作在10GHz中心频率、半波长间距的8阵元均匀线阵。
• 和差波束形成：实现了基于静态权重的数字波束形成，包括均匀加权的和波束与反相加权的差波束。
• 相位扫描仿真：支持-60°到+60°范围内的波束电扫描，通过对基础权重施加相位偏转实现波束指向控制。
• 波束宽度分析：自动计算不同扫描角度下的3dB波束宽度，验证了波束宽度随扫描角增大而变宽的理论（$1/costheta$ 规律）。
• 鉴角曲线生成：绘制法线方向的单脉冲鉴角曲线（S曲线），并提取线性区的鉴角斜率。
• 测角精度评估：通过蒙特卡洛方法（500次实验/角度）计算不同扫描角下的均方根误差（RMSE），分析信噪比和扫描角对精度的综合影响。

系统要求

• MATLAB R2016a 或更高版本
• 无需额外工具箱（仅使用核心MATLAB函数）

使用方法

直接运行主程序脚本即可。 程序运行后将：

1. 在控制台输出仿真步骤、系统参数摘要以及关键性能指标对比（如法线与60度扫描时的波束宽度和RMSE）。
2. 弹出一个包含四个子图的综合结果窗口，展示波束图、鉴角曲线、波束宽度变化趋势及测角误差趋势。

代码实现细节与算法分析

1. 系统初始化与阵列建模

程序首先定义了光速、频率（10GHz）、波长及阵元数量（N=8）。阵元位置采用了中心化设置（即以阵列几何中心为原点），这种设置有助于保持信号处理时的相位对称性，简化和差波束的相位特性分析。

2. 基础波束权重设计

程序实现了两种基础波束形成器：

• 和波束（Sum Beam）：采用全阵列均匀加权（全1向量），并进行了归一化处理，用于检测目标存在并作为相位基准。
• 差波束（Difference Beam）：采用左右子阵反相加权（左半阵为-1，右半阵为+1），在波束指向方向形成零陷，用于精确定位目标角度。

3. 法线方向性能仿真

在扫描循环开始前，程序首先计算了阵列法线（0度）方向的基准性能：

• 方向图绘制：遍历-90°到+90°计算阵列因子，生成和、差波束的幅度方向图。
• S曲线计算：在主瓣范围内扫描微小角度差，计算单脉冲比（$text{Real}(Delta / Sigma)$）。
• 斜率提取：通过线性拟合提取零点附近的鉴角斜率 $k_0$，该值作为后续大角度测角的基准参考。

4. 宽角扫描与误差分析（核心逻辑）

这是程序最核心的部分，通过循环遍历 -60° 到 60° 的扫描范围，执行以下操作：

A. 相位扫描实现 利用移相原理，将基础的和/差权重向量与目标扫描角度的导向矢量进行点乘，生成具有特定指向的扫描权重 $w_{scan}$。

B. 波束变宽效应分析 对于每一个扫描角度，程序在波束指向附近进行局部精细搜索，寻找功率下降3dB的上下边界，从而数值计算出当前的半功率波束宽度（HPBW）。仿真结果不仅存储数据，还在结果图中与 $1/cos(theta)$ 的理论曲线进行对比，验证了随扫描角增大，有效孔径减小导致波束变胖的物理现象。

C. 蒙特卡洛测角仿真 为了评估测角精度，程序在每个扫描点执行500次蒙特卡洛实验：

• 目标生成：在波束中心附近随机生成一个微小偏离的真实目标。
• 噪声叠加：模拟加性高斯白噪声（AWGN），信噪比设定为20dB。
• 单脉冲测角：计算接收信号的和差比 $mr$。
• 角度解算与修正：利用公式 $theta_{est} = theta_{scan} + frac{mr}{k_{approx}}$ 进行角度估计。特别地，代码中采用了 $k_{approx} = k_0 cdot cos(theta_{scan})$ 对鉴角斜率进行理论修正，以适应大角度扫描时的增益变化。
• RMSE计算：统计偏差的均方根值，量化测角误差。

5. 结果可视化

程序最终生成四幅子图：

1. 和差波束方向图：展示法线方向的标准和、差波束形态。
2. S曲线：展示单脉冲比与角度偏差的线性关系及拟合斜率。
3. 波束宽度趋势：展示仿真计算的波束宽度随扫描角增大而指数级上升的趋势，并与理论值吻合。
4. RMSE趋势：展示随着扫描角度增加，测角均方根误差显著增大的现象，直观反映了相控阵雷达大角度测角性能恶化的特性。

6. 辅助功能

代码包含一个独立的子函数 get_steering_vector，基于窄带信号模型 $exp(j cdot k cdot d cdot sintheta)$ 计算任意角度的均匀线阵导向矢量，保证了主程序的简洁性。