简介本资源是一套面向电子工程专业学生与天线设计初学者的MATLAB阵列天线仿真实践资料聚焦线性、平面及圆形三类典型天线阵列的方向图建模与性能分析。资源以精简实用为特点共3个文件2个MATLAB脚本polar_dB.m用于极坐标方向图绘制ARRAYS.m实现多种阵列因子计算与可视化、1份中文说明文档READ ME.doc总大小仅19KB便于快速下载与本地运行。已有298人学习下载适合课程设计、毕业设计或自学巩固天线理论知识的用户。通过该包读者可直接运行代码观察阵元间距、相位激励与阵列构型对波束指向、旁瓣电平及主瓣宽度的影响并借助文档理解关键参数设置逻辑与结果解读方法有效打通理论公式与MATLAB实操之间的关键环节。1. 用 MATLAB 快速建模线性阵列与圆形阵列天线从方向图仿真到波束扫描实战你手头有一份名为Chapter06.rar的压缩包解压后发现是 MATLAB 天线阵列仿真实验材料——这很常见高校《天线原理》《电磁场与微波技术》课程第六章常以线性阵列Linear Array和圆形阵列Circular Array为典型载体训练学生理解阵因子、方向图乘积定理、波束指向控制等核心概念。但真正卡住多数人的不是公式推导而是如何在 MATLAB 中不依赖 Antenna Toolbox 就写出可复现、可调参、可验证的阵列辐射模型尤其当你的 MATLAB 版本是 R2018a 或更早未预装天线工具箱或需将代码嵌入已有项目、避免额外依赖时纯脚本实现阵列建模就成为刚需。本文聚焦“零工具箱依赖”的底层建模路径用向量运算构建阵元位置、用相位加权实现波束扫描、用极坐标绘图呈现方向图、用mesh和surf可视化三维空间辐射特性。适合通信工程、雷达信号处理、射频系统设计方向的工程师与高年级本科生——只要你会写for循环和复数运算就能跑通全部流程。2. 线性阵列建模从阵元坐标到方向图计算的完整推导链2.1 理解线性阵列物理结构与关键参数定义线性阵列ULA, Uniform Linear Array是最基础的天线阵列形式由 N 个相同天线单元沿一条直线等距排列构成。其性能由三个核心参数决定阵元数量 N、相邻阵元间距 d单位波长 λ、以及每个阵元的激励幅度与相位。间距 d 直接影响栅瓣出现条件d ≥ λ/2 时易产生栅瓣而激励相位差 Δφ 则控制主波束指向角 θ₀满足 sinθ₀ Δφ·λ/(2πd)。MATLAB 中不调用phased.ULA类而是用基础矩阵运算显式构造阵元位置向量和阵因子响应这能让你彻底掌握方向图形成的数学本质——即远场辐射场是各阵元辐射场的矢量叠加其归一化方向图函数为$$ AF(\theta) \sum_{n0}^{N-1} w_n e^{j k d n \cos\theta} $$其中 $k 2\pi/\lambda$ 是波数$w_n$ 是第 n 个阵元的复数激励权重含幅度与相位$\theta$ 是观察角度从阵列法线方向起算。这个公式就是所有后续代码的根基。2.2 手动构建阵元位置与激励权重矩阵我们以 8 元线性阵列为例设定阵元间距 d 0.5λ避免栅瓣工作频率 f 3 GHz对应 λ ≈ 0.1 m并施加线性相位渐变实现 30° 波束扫描。以下代码完全基于 MATLAB 基础语法无需任何工具箱% 参数定义 N 8; % 阵元数量 d_lambda 0.5; % 阵元间距单位波长 theta0_deg 30; % 期望波束指向角度 f0 3e9; % 工作频率Hz c 3e8; % 光速m/s lambda c / f0; % 波长m % 构建阵元位置向量沿 x 轴单位米 x_pos (0:N-1) * d_lambda * lambda; % 列向量尺寸 N×1 % 计算扫描所需相位差单位弧度 k 2*pi / lambda; delta_phi k * d_lambda * lambda * sind(theta0_deg); % 注意sind() 输入为度 % 构建激励权重向量等幅线性相位 w exp(1j * (0:N-1) * delta_phi); % 复数权重尺寸 N×1提示sind()函数比sin()更安全因它直接接受角度制输入避免sin(theta0_deg*pi/180)的手动换算错误。权重向量w是列向量与位置向量x_pos维度一致为后续矩阵运算铺平道路。2.3 计算并绘制二维方向图方位面方向图计算的核心是对每个观测角度 θ计算所有阵元到该方向的路径差引起的相位偏移再与激励权重相乘求和。使用向量化运算可避免低效循环% 定义观测角度网格-90° 到 90°步长 0.5° theta_deg -90:0.5:90; theta_rad theta_deg * pi/180; % 向量化计算阵因子AF(theta) sum(w .* exp(j*k*x_pos*cos(theta))) % 利用广播机制x_pos (N×1) 与 cos(theta_rad) (1×M) 相乘得 N×M 矩阵 AF w. * exp(1j * k * x_pos * cos(theta_rad)); % 结果为 1×M 行向量 % 归一化并转为 dB 刻度 AF_dB 20*log10(abs(AF) / max(abs(AF))); % 绘制方向图 figure(Name, Linear Array Pattern (Theta Scan)); plot(theta_deg, AF_dB, b-, LineWidth, 1.5); xlabel(Angle \theta (degrees)); ylabel(Normalized Pattern (dB)); title(sprintf(ULA: N%d, d%.1f\lambda, Beam at %.0f^\circ, N, d_lambda, theta0_deg)); grid on; ylim([-40, 0]);关键参数说明x_pos * cos(theta_rad)利用 MATLAB 广播规则自动计算每个阵元在每个角度下的路径差单位米再乘以波数k得到相位差。w. * ...行向量w.1×N左乘 N×M 相位矩阵结果为 1×M 向量即每个角度的阵因子值。20*log10(...)将线性幅度转为对数刻度便于观察旁瓣与零点。此段代码输出的方向图清晰显示主瓣峰值位于 30°旁瓣电平约 -13 dB符合理论预期8 元均匀激励 ULA 的第一旁瓣理论值为 -13.2 dB。3. 圆形阵列建模坐标变换与三维方向图可视化3.1 圆形阵列几何建模与激励策略差异圆形阵列Circular Array将 N 个阵元均匀分布在半径为 R 的圆周上其优势在于全向扫描能力与对称性。与线性阵列不同圆形阵列的阵元位置需用极坐标转换为直角坐标且波束扫描需同时控制方位角 φ 和俯仰角 θ激励权重设计更复杂。关键参数包括阵元数 N、圆半径 R单位波长、以及目标扫描方向φ₀, θ₀。当 R 过大如 R 0.75λ时阵元间耦合增强需考虑互耦效应——但本节聚焦理想无耦合模型用纯几何方法建模。3.2 构建圆形阵列位置矩阵与多维方向图计算% 圆形阵列参数 N_circ 12; % 阵元数量 R_lambda 0.8; % 圆半径单位波长 phi0_deg 45; % 方位角扫描目标度 theta0_deg_circ 60; % 俯仰角扫描目标度 % 构建圆周上阵元位置x, y, z 坐标单位米 angle_circ linspace(0, 2*pi, N_circ1); % 生成 N_circ 个角度去掉最后一个重复点 angle_circ angle_circ(1:end-1); x_circ R_lambda * lambda * cos(angle_circ); y_circ R_lambda * lambda * sin(angle_circ); z_circ zeros(N_circ, 1); % 全部位于 xy 平面 % 计算扫描所需激励相位针对目标方向 phi0, theta0 % 目标方向单位矢量[cos(phi0)*sin(theta0), sin(phi0)*sin(theta0), cos(theta0)] ux cosd(phi0_deg) * sind(theta0_deg); uy sind(phi0_deg) * sind(theta0_deg); uz cosd(theta0_deg); % 每个阵元到目标方向的路径差 [x,y,z] · [ux,uy,uz] path_diff x_circ.*ux y_circ.*uy z_circ.*uz; delta_phi_circ k * path_diff; % 相位差向量 % 激励权重等幅相位补偿 w_circ exp(1j * delta_phi_circ); % 定义三维观测网格方位角 φ: 0~360°, 俯仰角 θ: 0~180° phi_deg 0:5:360; theta_deg_3D 0:5:180; [Phi, Theta] meshgrid(phi_deg, theta_deg_3D); Phi_rad Phi * pi/180; Theta_rad Theta * pi/180; % 计算三维方向图对每个 (phi, theta)计算单位方向矢量与各阵元点积 ux_grid cos(Phi_rad) .* sin(Theta_rad); uy_grid sin(Phi_rad) .* sin(Theta_rad); uz_grid cos(Theta_rad); % 向量化计算path_diff_grid x_circ*ux y_circ*uy z_circ*uz结果为 N×M×P % 使用 bsxfun 或隐式扩展R2016breshape 位置向量为 N×1×1广播至 M×P x3 reshape(x_circ, [], 1, 1); y3 reshape(y_circ, [], 1, 1); z3 reshape(z_circ, [], 1, 1); path_diff_grid x3.*ux_grid y3.*uy_grid z3.*uz_grid; % 阵因子sum over n of w_n * exp(j*k*path_diff_n) AF_3D squeeze(sum(w_circ. .* exp(1j * k * path_diff_grid), 1)); % 归一化为 dB AF_3D_dB 20*log10(abs(AF_3D) / max(abs(AF_3D(:))));3.3 用surf和polarplot呈现多维度辐射特性圆形阵列方向图需多视角展示极坐标图看方位面切片三维曲面图看全空间分布。% 绘制方位面切片固定俯仰角 θ90°即 xy 平面 figure(Name, Circular Array Azimuth Pattern); theta_slice_idx find(theta_deg_3D 90); AF_azimuth AF_3D_dB(theta_slice_idx, :); polarplot(Phi_rad(theta_slice_idx,:), AF_azimuth, r-, LineWidth, 1.5); title(sprintf(Circular Array: N%d, R%.1f\\lambda, Beam at \\phi%.0f^\\circ, ... N_circ, R_lambda, phi0_deg)); % 绘制三维方向图使用 surf view figure(Name, 3D Radiation Pattern); surf(Phi, Theta, AF_3D_dB, EdgeColor, none); shading interp; colormap(jet); colorbar; xlabel(Azimuth \phi (degrees)); ylabel(Elevation \theta (degrees)); zlabel(Pattern (dB)); title(3D Radiation Pattern of Circular Array); view([-37.5, 30]); % 调整视角便于观察主瓣注意squeeze(sum(..., 1))将 N×M×P 三维数组沿第 1 维阵元维度求和得到 M×P 二维方向图矩阵。surf绘图时X/Y 轴为角度网格Z 轴为 dB 值直观反映能量在球面上的分布密度。4. 阵列性能对比与关键参数调优指南4.1 线性 vs 圆形阵列方向图特性与适用场景对照表特性线性阵列ULA圆形阵列Circular Array扫描自由度单维仅方位角或仅俯仰角双维方位角 俯仰角主瓣宽度3dB≈ 0.886λ/(N·d)rad取决于半径 R 与 NR↑ 则主瓣↓但栅瓣风险↑旁瓣电平均匀激励-13.2 dBN8-10 ~ -12 dB典型 N12零点数量N-1 个在特定角度更多零点方向图更复杂硬件实现难度低单一直线布线中需精密圆周定位与馈电网络MATLAB 实现复杂度低一维向量运算中需三维坐标与广播计算该对比并非优劣判断而是选型依据若系统只需水平扫描如车载雷达ULA 足够且高效若需全空域覆盖如无人机测向圆形阵列不可替代。4.2 三个必调参数及其对方向图的定量影响调整以下参数时务必同步重绘方向图并观察变化这是理解阵列物理本质的最有效方式参数默认值调整效果推荐调试范围验证方法阵元间距 d0.5λd 0.5λ主瓣展宽分辨率↓d 0.5λ出现栅瓣虚假主瓣0.4λ ~ 0.6λULA观察方向图是否出现 -20 dB 的额外峰阵元数量 N8N↑主瓣变窄分辨率↑旁瓣↓但计算量↑N↓主瓣宽旁瓣高4 ~ 32依精度需求测量 3dB 主瓣宽度deg激励幅度 taper均匀全1加窗如 Hamming旁瓣↓-40 dB但主瓣↑分辨率↓Hamming / Taylor / Binomial对比加窗前后旁瓣电平与主瓣宽度例如对 8 元 ULA 施加 Hamming 窗w_tapered hamming(N) .* w; % w 为原相位权重 % 后续用 w_tapered 替代 w 计算 AF执行后方向图旁瓣降至 -42 dB但主瓣宽度从 14° 增至 18°体现分辨率与旁瓣的固有折衷。5. 故障排查与高频问题现场解决5.1 “方向图主瓣不在目标角度” 的三步定位法这是最常遇到的问题根源几乎总在相位计算环节。按顺序检查确认角度制/弧度制混用MATLAB 三角函数sin()/cos()输入为弧度sind()/cosd()输入为角度。扫描角theta0_deg若误用于sin(theta0_deg)会导致相位差错一个数量级。验证波数 k 计算k 2*pi/lambda而lambda c/f0。若f0单位错为 MHz如3e6而非3e9λ 被放大 1000 倍k 错误导致相位差失效。检查阵元序号与相位符号ULA 中第 n 个阵元n0,1,...,N-1的相位应为n * delta_phi。若写成(n1) * delta_phi或-n * delta_phi主瓣将反向或偏移。快速验证将theta0_deg设为 0°此时所有阵元应同相方向图主瓣必须严格位于 0°。若不在此处上述三步必有一处出错。5.2 “三维方向图出现异常尖峰或 NaN” 的内存与精度修复当N_circ较大24且角度步长过细2°时path_diff_grid数组可能超出内存或触发浮点溢出。解决方案降维计算不生成完整三维网格改用arrayfun分块计算% 对每个 (phi, theta) 点单独计算避免大数组 AF_3D_dB zeros(numel(theta_deg_3D), numel(phi_deg)); for i 1:numel(theta_deg_3D) for j 1:numel(phi_deg) ux cosd(phi_deg(j)) * sind(theta_deg_3D(i)); uy sind(phi_deg(j)) * sind(theta_deg_3D(i)); uz cosd(theta_deg_3D(i)); path_diff x_circ.*ux y_circ.*uy z_circ.*uz; AF_val sum(w_circ .* exp(1j * k * path_diff)); AF_3D_dB(i,j) 20*log10(abs(AF_val)/max_abs_AF); end end精度加固max(abs(AF_3D(:)))若为 0罕见但可能log10(0)返回-Inf。添加保护max_abs_AF max(abs(AF_3D(:))); if max_abs_AF 0, max_abs_AF eps; end % eps 为最小正浮点数 AF_3D_dB 20*log10(abs(AF_3D) / max_abs_AF);5.3 将.rar实验包转化为可复用的函数模块Chapter06.rar中的脚本通常是面向教学的单次运行代码。要工程化复用需封装为函数function [AF_dB, theta_deg] calc_ula_pattern(N, d_lambda, theta0_deg, f0) % CALC_ULA_PATTERN 计算线性阵列方向图 % 输入N-阵元数, d_lambda-间距(波长), theta0_deg-扫描角(度), f0-频率(Hz) % 输出AF_dB-归一化方向图(dB), theta_deg-角度向量(度) c 3e8; lambda c/f0; k 2*pi/lambda; x_pos (0:N-1) * d_lambda * lambda; delta_phi k * d_lambda * lambda * sind(theta0_deg); w exp(1j * (0:N-1) * delta_phi); theta_deg -90:0.5:90; theta_rad theta_deg * pi/180; AF w. * exp(1j * k * x_pos * cos(theta_rad)); AF_dB 20*log10(abs(AF) / max(abs(AF)eps)); end调用示例[pat, ang] calc_ula_pattern(16, 0.45, 22.5, 5e9); plot(ang, pat);此举使Chapter06.rar不再是孤立实验而成为你天线仿真工具箱中的一个可靠组件。本文还有配套的精品资源点击获取 MATLAB零工具箱实现天线阵列建模与波束扫描
简介本资源是一套面向电子工程专业学生与天线设计初学者的MATLAB阵列天线仿真实践资料聚焦线性、平面及圆形三类典型天线阵列的方向图建模与性能分析。资源以精简实用为特点共3个文件2个MATLAB脚本polar_dB.m用于极坐标方向图绘制ARRAYS.m实现多种阵列因子计算与可视化、1份中文说明文档READ ME.doc总大小仅19KB便于快速下载与本地运行。已有298人学习下载适合课程设计、毕业设计或自学巩固天线理论知识的用户。通过该包读者可直接运行代码观察阵元间距、相位激励与阵列构型对波束指向、旁瓣电平及主瓣宽度的影响并借助文档理解关键参数设置逻辑与结果解读方法有效打通理论公式与MATLAB实操之间的关键环节。1. 用 MATLAB 快速建模线性阵列与圆形阵列天线从方向图仿真到波束扫描实战你手头有一份名为Chapter06.rar的压缩包解压后发现是 MATLAB 天线阵列仿真实验材料——这很常见高校《天线原理》《电磁场与微波技术》课程第六章常以线性阵列Linear Array和圆形阵列Circular Array为典型载体训练学生理解阵因子、方向图乘积定理、波束指向控制等核心概念。但真正卡住多数人的不是公式推导而是如何在 MATLAB 中不依赖 Antenna Toolbox 就写出可复现、可调参、可验证的阵列辐射模型尤其当你的 MATLAB 版本是 R2018a 或更早未预装天线工具箱或需将代码嵌入已有项目、避免额外依赖时纯脚本实现阵列建模就成为刚需。本文聚焦“零工具箱依赖”的底层建模路径用向量运算构建阵元位置、用相位加权实现波束扫描、用极坐标绘图呈现方向图、用mesh和surf可视化三维空间辐射特性。适合通信工程、雷达信号处理、射频系统设计方向的工程师与高年级本科生——只要你会写for循环和复数运算就能跑通全部流程。2. 线性阵列建模从阵元坐标到方向图计算的完整推导链2.1 理解线性阵列物理结构与关键参数定义线性阵列ULA, Uniform Linear Array是最基础的天线阵列形式由 N 个相同天线单元沿一条直线等距排列构成。其性能由三个核心参数决定阵元数量 N、相邻阵元间距 d单位波长 λ、以及每个阵元的激励幅度与相位。间距 d 直接影响栅瓣出现条件d ≥ λ/2 时易产生栅瓣而激励相位差 Δφ 则控制主波束指向角 θ₀满足 sinθ₀ Δφ·λ/(2πd)。MATLAB 中不调用phased.ULA类而是用基础矩阵运算显式构造阵元位置向量和阵因子响应这能让你彻底掌握方向图形成的数学本质——即远场辐射场是各阵元辐射场的矢量叠加其归一化方向图函数为$$ AF(\theta) \sum_{n0}^{N-1} w_n e^{j k d n \cos\theta} $$其中 $k 2\pi/\lambda$ 是波数$w_n$ 是第 n 个阵元的复数激励权重含幅度与相位$\theta$ 是观察角度从阵列法线方向起算。这个公式就是所有后续代码的根基。2.2 手动构建阵元位置与激励权重矩阵我们以 8 元线性阵列为例设定阵元间距 d 0.5λ避免栅瓣工作频率 f 3 GHz对应 λ ≈ 0.1 m并施加线性相位渐变实现 30° 波束扫描。以下代码完全基于 MATLAB 基础语法无需任何工具箱% 参数定义 N 8; % 阵元数量 d_lambda 0.5; % 阵元间距单位波长 theta0_deg 30; % 期望波束指向角度 f0 3e9; % 工作频率Hz c 3e8; % 光速m/s lambda c / f0; % 波长m % 构建阵元位置向量沿 x 轴单位米 x_pos (0:N-1) * d_lambda * lambda; % 列向量尺寸 N×1 % 计算扫描所需相位差单位弧度 k 2*pi / lambda; delta_phi k * d_lambda * lambda * sind(theta0_deg); % 注意sind() 输入为度 % 构建激励权重向量等幅线性相位 w exp(1j * (0:N-1) * delta_phi); % 复数权重尺寸 N×1提示sind()函数比sin()更安全因它直接接受角度制输入避免sin(theta0_deg*pi/180)的手动换算错误。权重向量w是列向量与位置向量x_pos维度一致为后续矩阵运算铺平道路。2.3 计算并绘制二维方向图方位面方向图计算的核心是对每个观测角度 θ计算所有阵元到该方向的路径差引起的相位偏移再与激励权重相乘求和。使用向量化运算可避免低效循环% 定义观测角度网格-90° 到 90°步长 0.5° theta_deg -90:0.5:90; theta_rad theta_deg * pi/180; % 向量化计算阵因子AF(theta) sum(w .* exp(j*k*x_pos*cos(theta))) % 利用广播机制x_pos (N×1) 与 cos(theta_rad) (1×M) 相乘得 N×M 矩阵 AF w. * exp(1j * k * x_pos * cos(theta_rad)); % 结果为 1×M 行向量 % 归一化并转为 dB 刻度 AF_dB 20*log10(abs(AF) / max(abs(AF))); % 绘制方向图 figure(Name, Linear Array Pattern (Theta Scan)); plot(theta_deg, AF_dB, b-, LineWidth, 1.5); xlabel(Angle \theta (degrees)); ylabel(Normalized Pattern (dB)); title(sprintf(ULA: N%d, d%.1f\lambda, Beam at %.0f^\circ, N, d_lambda, theta0_deg)); grid on; ylim([-40, 0]);关键参数说明x_pos * cos(theta_rad)利用 MATLAB 广播规则自动计算每个阵元在每个角度下的路径差单位米再乘以波数k得到相位差。w. * ...行向量w.1×N左乘 N×M 相位矩阵结果为 1×M 向量即每个角度的阵因子值。20*log10(...)将线性幅度转为对数刻度便于观察旁瓣与零点。此段代码输出的方向图清晰显示主瓣峰值位于 30°旁瓣电平约 -13 dB符合理论预期8 元均匀激励 ULA 的第一旁瓣理论值为 -13.2 dB。3. 圆形阵列建模坐标变换与三维方向图可视化3.1 圆形阵列几何建模与激励策略差异圆形阵列Circular Array将 N 个阵元均匀分布在半径为 R 的圆周上其优势在于全向扫描能力与对称性。与线性阵列不同圆形阵列的阵元位置需用极坐标转换为直角坐标且波束扫描需同时控制方位角 φ 和俯仰角 θ激励权重设计更复杂。关键参数包括阵元数 N、圆半径 R单位波长、以及目标扫描方向φ₀, θ₀。当 R 过大如 R 0.75λ时阵元间耦合增强需考虑互耦效应——但本节聚焦理想无耦合模型用纯几何方法建模。3.2 构建圆形阵列位置矩阵与多维方向图计算% 圆形阵列参数 N_circ 12; % 阵元数量 R_lambda 0.8; % 圆半径单位波长 phi0_deg 45; % 方位角扫描目标度 theta0_deg_circ 60; % 俯仰角扫描目标度 % 构建圆周上阵元位置x, y, z 坐标单位米 angle_circ linspace(0, 2*pi, N_circ1); % 生成 N_circ 个角度去掉最后一个重复点 angle_circ angle_circ(1:end-1); x_circ R_lambda * lambda * cos(angle_circ); y_circ R_lambda * lambda * sin(angle_circ); z_circ zeros(N_circ, 1); % 全部位于 xy 平面 % 计算扫描所需激励相位针对目标方向 phi0, theta0 % 目标方向单位矢量[cos(phi0)*sin(theta0), sin(phi0)*sin(theta0), cos(theta0)] ux cosd(phi0_deg) * sind(theta0_deg); uy sind(phi0_deg) * sind(theta0_deg); uz cosd(theta0_deg); % 每个阵元到目标方向的路径差 [x,y,z] · [ux,uy,uz] path_diff x_circ.*ux y_circ.*uy z_circ.*uz; delta_phi_circ k * path_diff; % 相位差向量 % 激励权重等幅相位补偿 w_circ exp(1j * delta_phi_circ); % 定义三维观测网格方位角 φ: 0~360°, 俯仰角 θ: 0~180° phi_deg 0:5:360; theta_deg_3D 0:5:180; [Phi, Theta] meshgrid(phi_deg, theta_deg_3D); Phi_rad Phi * pi/180; Theta_rad Theta * pi/180; % 计算三维方向图对每个 (phi, theta)计算单位方向矢量与各阵元点积 ux_grid cos(Phi_rad) .* sin(Theta_rad); uy_grid sin(Phi_rad) .* sin(Theta_rad); uz_grid cos(Theta_rad); % 向量化计算path_diff_grid x_circ*ux y_circ*uy z_circ*uz结果为 N×M×P % 使用 bsxfun 或隐式扩展R2016breshape 位置向量为 N×1×1广播至 M×P x3 reshape(x_circ, [], 1, 1); y3 reshape(y_circ, [], 1, 1); z3 reshape(z_circ, [], 1, 1); path_diff_grid x3.*ux_grid y3.*uy_grid z3.*uz_grid; % 阵因子sum over n of w_n * exp(j*k*path_diff_n) AF_3D squeeze(sum(w_circ. .* exp(1j * k * path_diff_grid), 1)); % 归一化为 dB AF_3D_dB 20*log10(abs(AF_3D) / max(abs(AF_3D(:))));3.3 用surf和polarplot呈现多维度辐射特性圆形阵列方向图需多视角展示极坐标图看方位面切片三维曲面图看全空间分布。% 绘制方位面切片固定俯仰角 θ90°即 xy 平面 figure(Name, Circular Array Azimuth Pattern); theta_slice_idx find(theta_deg_3D 90); AF_azimuth AF_3D_dB(theta_slice_idx, :); polarplot(Phi_rad(theta_slice_idx,:), AF_azimuth, r-, LineWidth, 1.5); title(sprintf(Circular Array: N%d, R%.1f\\lambda, Beam at \\phi%.0f^\\circ, ... N_circ, R_lambda, phi0_deg)); % 绘制三维方向图使用 surf view figure(Name, 3D Radiation Pattern); surf(Phi, Theta, AF_3D_dB, EdgeColor, none); shading interp; colormap(jet); colorbar; xlabel(Azimuth \phi (degrees)); ylabel(Elevation \theta (degrees)); zlabel(Pattern (dB)); title(3D Radiation Pattern of Circular Array); view([-37.5, 30]); % 调整视角便于观察主瓣注意squeeze(sum(..., 1))将 N×M×P 三维数组沿第 1 维阵元维度求和得到 M×P 二维方向图矩阵。surf绘图时X/Y 轴为角度网格Z 轴为 dB 值直观反映能量在球面上的分布密度。4. 阵列性能对比与关键参数调优指南4.1 线性 vs 圆形阵列方向图特性与适用场景对照表特性线性阵列ULA圆形阵列Circular Array扫描自由度单维仅方位角或仅俯仰角双维方位角 俯仰角主瓣宽度3dB≈ 0.886λ/(N·d)rad取决于半径 R 与 NR↑ 则主瓣↓但栅瓣风险↑旁瓣电平均匀激励-13.2 dBN8-10 ~ -12 dB典型 N12零点数量N-1 个在特定角度更多零点方向图更复杂硬件实现难度低单一直线布线中需精密圆周定位与馈电网络MATLAB 实现复杂度低一维向量运算中需三维坐标与广播计算该对比并非优劣判断而是选型依据若系统只需水平扫描如车载雷达ULA 足够且高效若需全空域覆盖如无人机测向圆形阵列不可替代。4.2 三个必调参数及其对方向图的定量影响调整以下参数时务必同步重绘方向图并观察变化这是理解阵列物理本质的最有效方式参数默认值调整效果推荐调试范围验证方法阵元间距 d0.5λd 0.5λ主瓣展宽分辨率↓d 0.5λ出现栅瓣虚假主瓣0.4λ ~ 0.6λULA观察方向图是否出现 -20 dB 的额外峰阵元数量 N8N↑主瓣变窄分辨率↑旁瓣↓但计算量↑N↓主瓣宽旁瓣高4 ~ 32依精度需求测量 3dB 主瓣宽度deg激励幅度 taper均匀全1加窗如 Hamming旁瓣↓-40 dB但主瓣↑分辨率↓Hamming / Taylor / Binomial对比加窗前后旁瓣电平与主瓣宽度例如对 8 元 ULA 施加 Hamming 窗w_tapered hamming(N) .* w; % w 为原相位权重 % 后续用 w_tapered 替代 w 计算 AF执行后方向图旁瓣降至 -42 dB但主瓣宽度从 14° 增至 18°体现分辨率与旁瓣的固有折衷。5. 故障排查与高频问题现场解决5.1 “方向图主瓣不在目标角度” 的三步定位法这是最常遇到的问题根源几乎总在相位计算环节。按顺序检查确认角度制/弧度制混用MATLAB 三角函数sin()/cos()输入为弧度sind()/cosd()输入为角度。扫描角theta0_deg若误用于sin(theta0_deg)会导致相位差错一个数量级。验证波数 k 计算k 2*pi/lambda而lambda c/f0。若f0单位错为 MHz如3e6而非3e9λ 被放大 1000 倍k 错误导致相位差失效。检查阵元序号与相位符号ULA 中第 n 个阵元n0,1,...,N-1的相位应为n * delta_phi。若写成(n1) * delta_phi或-n * delta_phi主瓣将反向或偏移。快速验证将theta0_deg设为 0°此时所有阵元应同相方向图主瓣必须严格位于 0°。若不在此处上述三步必有一处出错。5.2 “三维方向图出现异常尖峰或 NaN” 的内存与精度修复当N_circ较大24且角度步长过细2°时path_diff_grid数组可能超出内存或触发浮点溢出。解决方案降维计算不生成完整三维网格改用arrayfun分块计算% 对每个 (phi, theta) 点单独计算避免大数组 AF_3D_dB zeros(numel(theta_deg_3D), numel(phi_deg)); for i 1:numel(theta_deg_3D) for j 1:numel(phi_deg) ux cosd(phi_deg(j)) * sind(theta_deg_3D(i)); uy sind(phi_deg(j)) * sind(theta_deg_3D(i)); uz cosd(theta_deg_3D(i)); path_diff x_circ.*ux y_circ.*uy z_circ.*uz; AF_val sum(w_circ .* exp(1j * k * path_diff)); AF_3D_dB(i,j) 20*log10(abs(AF_val)/max_abs_AF); end end精度加固max(abs(AF_3D(:)))若为 0罕见但可能log10(0)返回-Inf。添加保护max_abs_AF max(abs(AF_3D(:))); if max_abs_AF 0, max_abs_AF eps; end % eps 为最小正浮点数 AF_3D_dB 20*log10(abs(AF_3D) / max_abs_AF);5.3 将.rar实验包转化为可复用的函数模块Chapter06.rar中的脚本通常是面向教学的单次运行代码。要工程化复用需封装为函数function [AF_dB, theta_deg] calc_ula_pattern(N, d_lambda, theta0_deg, f0) % CALC_ULA_PATTERN 计算线性阵列方向图 % 输入N-阵元数, d_lambda-间距(波长), theta0_deg-扫描角(度), f0-频率(Hz) % 输出AF_dB-归一化方向图(dB), theta_deg-角度向量(度) c 3e8; lambda c/f0; k 2*pi/lambda; x_pos (0:N-1) * d_lambda * lambda; delta_phi k * d_lambda * lambda * sind(theta0_deg); w exp(1j * (0:N-1) * delta_phi); theta_deg -90:0.5:90; theta_rad theta_deg * pi/180; AF w. * exp(1j * k * x_pos * cos(theta_rad)); AF_dB 20*log10(abs(AF) / max(abs(AF)eps)); end调用示例[pat, ang] calc_ula_pattern(16, 0.45, 22.5, 5e9); plot(ang, pat);此举使Chapter06.rar不再是孤立实验而成为你天线仿真工具箱中的一个可靠组件。本文还有配套的精品资源点击获取