简介围绕Ti AWR2243级联板发射端波束形成相控阵的实践数据与Matlab代码包面向无线通信、雷达感知方向的学生和工程师用于学习相控阵波束控制、相位计算及硬件级联配置。压缩包共38个文件、64.36MB其中24个bin为采集/运行数据11个m脚本覆盖相位控制器计算、波束形成算法与结果绘图1个mat保存天线阵列中间变量另含vsdx示意图和txt说明便于按代码、数据、文档分类查阅数据覆盖静止/动目标及单/双目标场景可与脚本配合对照分析。已有2217人学习下载。通过运行脚本可复现从阵列建模、相位偏置计算到波束指向验证的全流程结合实测数据能观察不同方向下的波束成形效果、接收功率与信噪比变化并掌握AWR2243多通道独立控制的调试思路适合作为从仿真到硬件实践的系统参考。1. 把 AWR2243 级联板的发射波束转到指定角度实践前先想清楚这几件事很多人第一次拿到 AWR2243 级联板想跑发射端波束形成相控阵时下意识会以为只要在配置里把多个发射通道全打开雷达就能自动把能量集中到某个方向。真这样做过一次就会发现方向图要么乱跳要么主瓣宽度和预期差出好几倍。真正能落地的做法是先把每个发射通道的相位差算准、写进寄存器再用原始回波数据画出方向图来验证。这篇笔记带着你把这块板从原理、配置、数据读取做到方向图对比配套代码可以直接用于你自己的数据。尤其适合做车载雷达、交通监控和工业感知的工程师如果你手里只有一块单芯片雷达板这套思路同样能照搬只是通道数变少、波束变宽。2. 发射端波束形成的原理与硬件选型为什么相控阵不是简单开多路天线2.1 等距线阵里的相位差公式是整个实践的地基在开始配置之前先理解发射端波束形成到底在做什么。把阵列上的每个发射天线看成一个点源它们以同样的频率、同样的波形发射信号唯一不同的是各路信号之间存在一个初始相位差。在远场空间中任意一个观察方向上来自不同天线的信号会叠加。当这些信号的相位在某个方向恰好对齐叠加幅度达到最大在其他方向被部分抵消幅度变小。因为控制的是相位差所以把这种阵列叫作相控阵。假设阵列是等距线阵阵元编号为 0 到 N-1阵元间距为 d。对于从法线方向偏转角度 θ 的目标或观察点第 k 个阵元到该方向远场点的距离与第 0 个阵元相差 k·d·sin(θ)折算成相位就是 k·2π·d·sin(θ)/λ。要让波束指向 θ0就在发射端给第 k 个阵元补一个相反的相位增量 k·Δφ其中 Δφ 2π·d·sin(θ0)/λ。这样每个阵元在 θ0 方向上的空间相位差刚好被发射相位差抵消所有信号同相叠加合成波束就指向了 θ0。这里有几个容易混淆的点。第一d 必须按真实的天线口径来取不能随手填半个波长。很多板卡为了控制栅瓣把天线间距设计成略小于半波长比如 0.9mm 到 1.1mm 之间。计算相位时如果把 d 设成波长的二分之一在 77GHz 下就是 1.9mm 左右30° 指向的相位差会差出几十度方向图在实测里完全对不上。第二相位旋转器能提供的相位码位数有限连续算出来的相位值要先量化到硬件支持的相位码。量化误差会让旁瓣稍微抬高但不会让主瓣消失。这段为什么要写这么细因为后面所有配置文件、数据处理都依赖这个公式。我看到过不少人拿着板卡先跑默认配置再把波束指向参数随手填一个数最后测出来主瓣乱跳就开始怀疑硬件有问题。其实大多数情况是相位差公式里的 d 没对齐或者相位旋转器的量化方式没有换算对。把地基打好后面才好排查。2.2 AWR2243 级联板把几十路射频链路做进了同一块板AWR2243 本身是一颗工作在 76~81GHz 的毫米波射频收发芯片单颗芯片集成有 3 路发射通道和 4 路接收通道。所谓级联板就是把多颗这样的芯片通过同步参考时钟、本振分配网络和复位控制电路绑在一起对外暴露为一个更大的雷达前端。常见的 AWR2243 级联板为 4 片级联合计 12 路发射、16 路接收。这个数量的发射通道意味着相控阵的阵元数可以从单芯片的 3 跳到 12主瓣宽度收窄为原来的四分之一等效发射功率也明显提升。级联板最核心的价值不是“通道数量多”而是“相位相干”。发射波束形成要求每一路发射信号在同一个时刻维持确定的相位关系否则阵列合成不到一起。单片 AWR2243 内部的三个发射通道天然共享同一本振相位关系容易保证。跨芯片就要麻烦得多主芯片和从芯片必须同步锁定才能让本振相位差固定在一个已知常数上。级联板在硬件上已经把主从之间的同步信号和时钟树做好了用户不需要自己调锁相环只需要检查状态寄存器里的同步标志。对一线工程师来说选 AWR2243 级联板还有一个实际原因发射通道内置了相位旋转器不需要外接移相器。工作频率到了 77GHz在外面的射频链路里加移相器板级走线和寄生参数都很难控制成本也高。片上相位旋转器虽然分辨率有限但胜在稳定、可重复。配置一个角度本质上就是在写一组成品寄存器值不需要动任何模拟器件。这让自动测试和批量校准都变得可行。不过也应该说清楚级联板的优势同时带来调试复杂度。12 路发射通道并不是所有都排在一根直线上。有的板卡把 12 路分成两排或者夹在接收天线中间计算阵列方向图时必须知道每个天线的实际位置坐标。拿不到天线封装坐标时方向图拟合很难做。这也是我建议拿到板卡后先读设计文档而不是直接拿通用公式到处套的原因。真实坐标可以从板卡丝印、供应商设计包或者高分辨率照片里量出来。2.3 发射波束形成和 MIMO 是两回事别把资源用错不少做雷达的人一提到多路发射第一反应就是 MIMO即用时分或编码方式让不同发射通道发不同信号配合多路接收合成虚拟大阵列。MIMO 的发射通道实际上是“分时复用”的每个通道单独发接收端通过正交编码把各通道分离开。发射波束形成完全不同它要求所有发射通道在同一时刻、以相同波形、不同相位向同一个方向辐射目的是把能量集中起来提高这个方向上的信噪比。这两种用法会互相冲突。如果你想同时获得 MIMO 的角度分辨率和发射波束形成的增益在 AWR2243 级联板上通常做不到完全兼容因为发射通道同时发射就没有编码自由度无法再分离每条通道。工程上常见的折中是分成两个阶段先做 MIMO 模式探测目标方位再根据目标方位设置发射波束指向该角度用波束形成模式精测距离和速度。这样做的好处是既不用增加硬件又能在关键目标上多拿到几个 dB 的增益。选择级联板而不选单芯片的第三个理由就在这12 路发射的波束形成增益理想情况下接近 12 路功率合成即使扣除插损和相位误差实测也能拿到 9dB 以上的有效增益。对于远距离检测或弱目标回波这 9dB 可能就是能不能检测到的分界线。而单芯片 3 路发射增益只有 4.8dB 左右明显不够用。所以这个标题里的“相控阵”不是概念演示而是实打实要能让目标回波幅度上一个台阶。3. 让发射波束转起来从角度到相位表的最小闭环3.1 用一段 Python 把波束指向角换算成各通道相位基线先写一个最基础的工具函数输入目标角度和天线间距输出每个发射通道要设置的相位。以 77GHz 为例光速按 3e8 算波长约 3.9mm。下面的代码里我按等距线阵计算只做相位递增不处理天线位置排列的特殊情况如果你的板卡天线不是等距或不在一条直线上需要换成实际坐标再算。import numpy as np def calc_tx_phase(steer_deg, fc77.0e9, spacingNone, tx_num12): 计算等距线阵各发射通道的理论相位度 steer_deg: 波束目标指向角法线方向为 0 度 fc: 工作频率 spacing: 天线间距默认取半波长 tx_num: 发射通道数 c 3.0e8 wavelength c / fc if spacing is None: spacing 0.5 * wavelength # 相邻阵元相位差单位弧度 delta_rad 2 * np.pi * spacing * np.sin(np.deg2rad(steer_deg)) / wavelength delta_deg np.rad2deg(delta_rad) # 各通道相位按 0, delta, 2*delta ... 递推并归一到 [-180, 180) phase_deg [(i * delta_deg 180) % 360 - 180 for i in range(tx_num)] return np.array(phase_deg) # 指向 30 度通道数 12天线间距半波长 phase_table calc_tx_phase(30.0) for idx, phase in enumerate(phase_table): print(fTX{idx1:02d}: {phase:8.3f} deg)这段代码的关键参数有三个目标角度steer_deg、工作频率fc、天线间距spacing。如果你只是验证算法用默认的半波长就够了如果你想要实测方向图贴合理论曲线spacing必须换成板卡设计文件里的真实值。代码注释里也标明了这里算出来的是“理论相位”还不能直接写寄存器。后面要做量化。为什么要在这一步先用 Python 而不是直接写配置文件因为你可以快速批量生成不同角度下的相位表比如 0°、15°、30°、45°用同一段脚本一次性算完。手工在图形界面里填 12 个通道的相位填错一个字符都难察觉脚本生成再人工核对至少不会犯低级算术错误。另外注意脚本假设的是所有发射通道相位按顺序递增。如果板卡上天线编号和物理位置不是顺序排列你必须先把相位重新映射到物理天线坐标上。否则即使配置写进去了阵列合成出来的方向图仍是乱的。这块是最容易“翻车”的地方务心用天线坐标做一遍映射。3.2 相位量化与寄存器映射别拿着连续角度直接写寄存器AWR2243 发射通道里的相位旋转器是数字控制的它能支持的相位步进不是 1 度也不是 0.1 度而是由硬件和固件共同决定的一个离散间隔。常见实现里相位码从 0 到某个上限每个码对应一个相位但码和角度并不一定线性。你拿 Python 算出来 83.5 度硬件附近可能只有 82.6875 度和 84.375 度两档要选哪个得看固件手册里的相位码映射表。我一般会在拿到板卡后先做一个“相位码扫描”实验把同一个目标放在固定位置单独开某一个发射通道只改这一个通道的相位码记录目标回波的相位。这样就能测出当前固件下相位码和实际相位的对应关系。这个实验听起来繁琐但一次做完能省掉后面大量排查时间。很多团队跳过了这一步直接用角度值写寄存器结果方向图主瓣位置偏出好几度还在怀疑算法。下面的表格是一张典型的“目标角度到理论相位”的中间表供你核对自己的 Python 输出。里面假设 12 路发射、天线间距为半波长、工作频率 77GHz、波束指向 30 度发射通道理论相位度TX010.000TX0290.000TX03180.000TX04-90.000TX050.000TX0690.000TX07180.000TX08-90.000TX090.000TX1090.000TX11180.000TX12-90.000这张表看着很有规律但它只在半波长等距线阵下成立。如果你的天线间距变成 0.9mm30 度对应的相位差大约是 84 度马上就会打乱这个漂亮的循环。所以表格只适合用来目测脚本有没有算错不适合直接抄。真正写进寄存器的相位码必须用上一节说的扫描结果做量化。量化这一步常规做法是写一个映射函数def quantize_phase(phase_deg, phase_values_deg): phase_deg: 理论连续相位 phase_values_deg: 固件支持的所有相位离散值单位度 返回最接近的离散相位和对应索引 arr np.asarray(phase_values_deg) idx np.argmin(np.abs(arr - phase_deg)) return arr[idx], idx这个函数本身很简单但它暴露了一个实际问题如果硬件相位分辨率粗糙比如只有 22.5 度一档那么算出来的 90 度没问题84 度就会被量化到 90 度造成 6 度指向误差。所以不要等方向图测出来再骂硬件在量化之前就要评估这个误差是否能接受。3.3 初始化级联板并写入相位表关键参数和配置模板相位表计算好之后接下来要完成 AWR2243 级联板的初始化和配置。常见做法是通过 mmWave Studio 或配套的雷达驱动库按“复位 - 加载固件 - 配置 profile - 配置 chirp - 配置 TX - 配置 frame”的顺序操作。顺序不能乱尤其是 TX 相位配置要放在 chirp 启动之前。采集参数的设定直接影响后面数据处理常见的一套方向图验证参数如下参数值说明起始频率77 GHz避开频带边缘频率斜率40 MHz/us决定差拍频率和距离分辨率ADC 采样率20 Msps保证差拍频率不混叠采样点数256距离分辨率和数据量折中Chirp 数64做相干平均提高信噪比Frame 数16多次采集取统计值TX 通道数12全部打开同时发射RX 通道数16接收全部打开如果你用雷达驱动库直接写配置伪代码大致是这样。实际函数名以你拿到的 SDK 为准但逻辑不变// 伪代码初始化并配置发射波束 // 1. 复位并等待同步 reset_device(); wait_for_sync(); // 2. 配置射频前端的 profile struct profile_cfg prof; prof.start_freq 77.0e9; prof.slope 40.0e12; prof.samples 256; set_profile(prof); // 3. 配置发射通道全部开启 set_tx_channel(0xFFF); // 4. 写入 12 个发射通道的相位码 // 数组中第 k 个元素对应第 k1 路 TX float phase_deg[12] {0, 90, 180, -90, 0, 90, 180, -90, 0, 90, 180, -90}; for (int i 0; i 12; i) { set_tx_phase(i, phase_deg[i]); } // 5. 启动 frame start_frame();这里要重点解释的是set_tx_phase这一步。有些驱动封装里set_tx_phase接受的是相位码而不是角度你需要在调用之前把角度转成码。另一个容易踩的情况是部分固件会要求所有 TX 通道的相位在同一个 chirp 配置中一次性写入不能逐通道独立启动否则从片通道没有同步参考写入后不生效。你调试时如果发现某个通道怎么配都没反应先查是不是跨芯片通道的使能时序有问题。配置完成后先不急着测方向图可以做一个快速检查把发射波束指向 0 度全通道相位设为 0用手持角反射器在法线方向扫描看目标回波幅度然后把波束切到 30 度把角反射器移到 30 度位置如果回波明显增强说明发射波束形成的基本链路已经通了。这个检查只需要十几分钟但能过滤掉大部分配置错误。4. 数据和代码配套从 ADC 原始数据到发射方向图的验证4.1 采集配置怎么定原始 bin 文件怎么读方向图验证需要一组“波束指向不同角度”的采集数据。我不建议把所有角度的数据写进同一个文件那样后面解析时非常容易把索引搞混。常见做法是每个波束指向角单独存一个 bin 文件文件名直接包含角度比如data_angle_0.bin、data_angle_15.bin。采集时把角反射器固定在远场某个位置雷达转台或支架固定不动每次只改发射波束的相位配置重新采集。读取原始数据的代码要格外注意字节序。mmWave Studio 导出的 ADC 原始数据通常是 int16 交错存储即实部、虚部交替。直接用np.complex64读取偶尔会在对齐上出问题更稳妥的做法是先按 int16 读进来再用 view 转成复数。import numpy as np def load_adc_bin(file_path, frame_cnt, chirp_cnt, num_sample, num_rx): 读取 mmWave Studio 导出的 ADC 原始 bin 文件 返回形状: (frame_cnt, chirp_cnt, num_sample, num_rx) raw_iq np.fromfile(file_path, dtypenp.int16) # 相邻两个 int16 拼成一个复数小端模式 adc raw_iq.view(np.complex64).reshape(-1) expected frame_cnt * chirp_cnt * num_sample * num_rx if adc.size ! expected: raise ValueError(f数据长度不匹配: {adc.size} vs {expected}) data adc.reshape(frame_cnt, chirp_cnt, num_sample, num_rx) return data这里num_rx是接收通道数。如果你用的是 4 片级联板的 16 路接收num_rx就是 16如果你的板卡只接了 4 路接收那要改成 4。需要特别注意的是发射波束形成时 12 路 TX 是同时发射的接收数据不能区分出发射通道所以 bin 文件里没有 TX 维度。千万不要在 reshape 里强行加一个tx_num那样会解读出一堆乱数据。如果你的采集软件导出的不是int16交错而是float32或int16分离的 I/Q 两列读取方式要相应调整。判断方法很简单把raw_iq.size和文件大小对比一下文件大小等于实际样本数 * 2 * 2字节时就是 int16 交错。这个细节卡住过不少人宁可多花两分钟确认也不要闷头处理。4.2 距离维 FFT 找峰值用静态角反射器做目标方向图验证需要一个强而稳定的回波数学上最理想的目标是金属角反射器。它放在远场某个固定角度时回波可以看成点目标的散射。处理方式对每一帧、每一个 chirp 做距离维 FFT然后取目标所在距离 bin 的幅度再做平均值得到该波束指向下的回波强度。def range_fft(samples, nfft1024, fs20e6, slope40e12): samples: 一个 chirp 的 ADC 采样序列 nfft: FFT 点数 fs: ADC 采样率 slope: chirp 频率斜率 返回距离谱和距离向量 win np.hanning(samples.shape[-1]) spectrum np.fft.fft(samples * win, nfft) # 距离 bin 换算 freq_axis np.fft.fftfreq(nfft, d1/fs) dist freq_axis * 3e8 / (2 * slope) return spectrum, dist def extract_target_amplitude(data, target_bin, frame_cnt, chirp_cnt): 从多帧多 chirp 数据中提取目标距离 bin 的幅度均值 amp_sum 0.0 count 0 for f_idx in range(frame_cnt): for c_idx in range(chirp_cnt): spec, dist range_fft(data[f_idx, c_idx, :, 0]) # 先只取第 0 接收通道 amp_sum np.abs(spec[target_bin]) count 1 return amp_sum / count这段代码里我刻意先只用第 0 接收通道原因是实际验证中可以先把接收通道固定看某个信道避免 16 路接收相位不一致带来的额外影响。如果你想综合所有接收通道可以对 16 路幅度求平均但要注意接收通道本身的幅相不一致也会折算进方向图结果无法区分是发射问题还是接收问题。所以第一轮验证尽量看单一接收通道。target_bin怎么选根据距离公式先让角反射器放在已知距离比如 8 米用某个波束指向的采集数据做一次完整距离谱查找最大值所在 bin再把那个 bin 固定下来。后面的角度数据都取同一个 bin不要每帧都重新找最大值因为噪声会在不同帧之间跳动导致幅度曲线毛刺。4.3 发射方向图怎么画峰值幅度随指向角的变化曲线验证方向图时目标是固定不动的发射波束指向角在改变。当波束正好对准目标时回波幅度出现峰值波束偏离目标回波幅度下降。于是一组“幅度 vs 波束指向角”的曲线就近似等效为发射方向图。实际操作中目标角度和波束指向角通常设置成相等比如目标放在 30 度那么你就测一组波束指向从 -60 到 60 度的数据看峰值是否出现在 30 度附近。下面的代码把实测数据和理论方向图叠在一起画。理论方向图假设 12 路发射、等幅、天线间距半波长。import numpy as np import matplotlib.pyplot as plt def theoretical_pattern(scan_deg, steer_deg, tx_num12, d_lambda0.5): 等距线阵理想发射方向图 scan_deg: 扫描角 steer_deg: 波束指向角 d_lambda: 天线间距除以波长 scan np.deg2rad(scan_deg) steer np.deg2rad(steer_deg) scan_phase 2 * np.pi * d_lambda * np.sin(scan) steer_phase 2 * np.pi * d_lambda * np.sin(steer) # 每个发射通道的相位 扫描相位 - 指向相位 arr np.exp(1j * (np.arange(tx_num)[:, None] * scan_phase - np.arange(tx_num)[:, None] * steer_phase)) return np.abs(np.sum(arr, axis0)) # 实测数据angle_list 是波束指向角amp_list 是对应回波幅度 angle_list np.array([-60, -45, -30, -15, 0, 15, 30, 45, 60]) amp_list np.array([100, 210, 480, 780, 1150, 980, 2100, 730, 260]) amp_db 20 * np.log10(amp_list / amp_list.max()) # 归一化到 0dB # 理论方向图指向目标角度 30 度 theory_angle np.linspace(-90, 90, 181) theory theoretical_pattern(theory_angle, steer_deg30.0) theory_db 20 * np.log10(theory / theory.max() 1e-6) plt.figure(figsize(8, 5)) plt.plot(theory_angle, theory_db, -, labeltheory) plt.plot(angle_list, amp_db, o-, labelmeasured) plt.xlabel(beam pointing angle (deg)) plt.ylabel(normalized amplitude (dB)) plt.grid(True) plt.legend() plt.show()代码里我把实测和理论都归一化到了 0dB方便比较主瓣位置和宽度。首次对照时不要指望两条曲线重合得很完美因为天线间距、接收通道不一致、环境多径都会影响旁瓣。你要关注的核心指标有三个峰值位置是否在目标角度附近主瓣宽度是否接近理论值旁瓣电平是否明显偏高。如果峰值偏了回去查相位量和天线坐标如果主瓣宽了大概率是目标不满足远场条件或接收通道饱和如果旁瓣高查通道幅度一致性。这张图也是后续所有校准的基准。我习惯在每次调整硬件或驱动版本后重新画一条实测曲线存起来。方向图这种东西过一段时间再测可能因为环境变化、板卡老化发生缓慢漂移没有历史曲线就很难发现。把基线数据留好后面做排障会省很多事。5. 发射波束形成避坑指南从相位校准到近场判据的五个血泪点这章记录的都是我在调试过程中真实撞过的坑每一个都让项目至少多花了两三天。按“现象 - 原因 - 解决”来写你遇到类似问题时可以直接对照。5.1 波束指向偏移且随频率漂移现象目标放在 30 度波束也配置成 30 度但实测峰值出现在 22 度附近把起始频率从 77GHz 改到 78GHz偏移量又变了。原因最常见有两个一是没有做发射通道相位校准相位旋转器本身的初始相位不一致二是 Python 计算相位表时用了半波长的默认间距但板卡实际天线间距不是半波长。解决先用 0 度指向扫描一遍看峰值是否真的在 0 度。如果不在把偏差记为初始相位加到每个通道上再核对板卡天线坐标文档按实际 d 重算相位表。做完这两步峰值基本能回到目标角度附近。5.2 方向图旁瓣不对称甚至一侧抬高 5dB现象主瓣出现在目标角度附近但左右旁瓣高度差很明显看起来不像等距线阵理论方向图。原因各发射通道的输出幅度不一致通常靠近供电入口的通道输出偏高射频走线长的通道输出偏低低频通道之间的耦合也会让某些通道等效辐射效率变差。解决单独开每个通道在远场同一个位置测单通道回波幅度得到 12 个幅度权重。然后在相位配置里对每个通道乘上这个权重再重新测量方向图。如果做完后旁瓣仍然不对还要检查是不是某个发射通道根本没有工作用热成像或近场探针能快速定位。5.3 跨芯片通道相位跳变方向图主瓣碎裂现象方向图不是单峰而是在多个角度出现差不多高度的峰像一把梳子。原因级联板的主片和从片之间同步没有锁定或者启动时序里 frame 在同步建立之前就开始了。跨芯片的发射通道如果本振相位不固定相位差会随机跳变阵列自然合不出单一波束。解决初始化后先读主从芯片的状态寄存器确认同步锁定位为 1每次复位后加一个固定的等待时间让 PLL 和同步链路完全稳定再启动 chirp。如果还是不定检查参考时钟走线附近有没有强干扰源。5.4 近场目标把方向图测糊现象目标放在 1 米左右方向图主瓣比理论宽出一倍换到不同距离主瓣宽度还会变化。原因发射波束形成理论的成立前提是远场平面波近似。目标太近时目标到阵列每个发射天线的距离差不是简单的 d·sinθ相位补偿模型失效。解决方向图测试遵循远场条件。对于口径 D 的阵列远场距离通常按 2D²/λ 估算。12 路发射阵列口径大约 5 到 10 厘米77GHz 下远场距离在 5 到 10 米以上。所以测试时把角反射器放到 8 米外而不是放在桌面上扫。没有大场地的话可以缩短阵列口径选中间若干通道测。5.5 单帧数据信噪比太低方向图曲线毛刺现象重复测量同一组角度曲线峰值位置会抖动旁瓣起伏超过 3dB。原因每帧只取了一个 chirp 做目标幅度计算信噪比不够或者测试环境里有车辆、人员走动造成多径。解决把同一帧内 64 个 chirp 的目标距离 bin 幅度做相干平均如果允许再对多个 frame 取中位数避免个别帧里的噪声尖峰。测试场地尽量选空旷角落角反射器周围两米内不要有金属物。环境干净之后方向图曲线会平滑很多。6. 进阶用离线方向图拟合反推每路发射通道的幅度和相位误差当你已经把方向图测稳下一步可以做纯软件校准用实测方向图反推每路发射通道的复权重误差。这个技巧的价值在于它能量化地告诉你“哪一路相位偏了多少、幅度掉了多少”而不是靠肉眼猜。方法很简单远场下方向图可以写成每个通道权重和阵列导向矢量的线性叠加。在角度 θ 处理想阵列响应为 E(θ) Σ w_k·exp(j·k·2π·d·sinθ/λ)。其中 w_k 是第 k 个发射通道实际复权重包含了幅度和相位误差。我们测量多个角度的方向图数据用最小二乘估计 w_k就能得到误差分布。下面给出一段紧凑的实现import numpy as np def estimate_channel_weights(angle_meas, pattern_meas, tx_num12, d_lambda0.5): 用实测方向图反推各通道复权重 angle_meas: 实测角度度 pattern_meas: 对应角度的复幅度如果不带相位用幅度近似 theta np.deg2rad(angle_meas) # 构造导向矢量矩阵 A: shape (len(angle_meas), tx_num) index np.arange(tx_num)[None, :] A np.exp(1j * index * 2 * np.pi * d_lambda * np.sin(theta)[:, None]) # 最小二乘解 w_est, residuals, rank, s np.linalg.lstsq(A, pattern_meas, rcondNone) amp_error np.abs(w_est) phase_error np.angle(w_est) * 180 / np.pi return amp_error, phase_error使用这个函数有一点要留意如果你的pattern_meas只有幅度数据没有相位信息最小二乘只能解出各通道相对幅度相位误差无法完整还原。想反推相位误差最好在采集时让接收通道同时给出 I/Q然后以第一个通道为参考提取相对相位。或者你用多组不同波束指向角的数据做联合拟合也能提高相位误差的可辨识度。我现在的习惯是拿到一块新级联板第一件事不是直接测探测指标而是先跑一遍 0 度到正负 60 度方向图用反推出幅度和相位误差把每条通道的校准值存成一个表格之后每次配置都自动叠加这份校准值。这一步看起来多花了半天但后面所有实验数据都更可信排查问题也快很多。雷达这种设备校准做得越早后期返工越少先把发射通道的底查清楚再做探测、算法才能不被方向图上的“玄学”干扰。希望这套方法和数据代码能帮你的波束形成实践少走几段弯路。本文还有配套的精品资源点击获取
