IADATASDG / SENSOR SYSTEMS / TECHNICAL NOTE

Radar 原理与仿真架构

雷达器件结构、探测原理、电磁散射、几何传播、信号检测、模块流水线与标定方法。

结构与探测物理 FEATURE双路径传播FAST / WAVEFORM2026-11-01

01 / RADAR STRUCTURERadar 结构

真实毫米波雷达由波形/时钟、发射通道、天线、雷达罩、接收阵列、混频与采样、数字处理和数据接口组成。场景目标及传播介质构成外部回波链路。[1]

雷达结构:波形时钟、发射链、目标回波、接收阵列、处理和接口RADAR MODULE · 器件结构与仿真抽象对应时钟 / 波形源fc · B · Tr · CPITx 发射链PA · 天线增益Tx 天线视场 / 发射射线雷达罩 / 传播介质插损 · 雨雾 · 湿膜Rx 接收阵列到达角 · 虚拟孔径混频 / ADC仿真由频域模型抽象DSP / 检测输出R · vr · DOA · RCS场景目标 / 表面几何 · 材质 · 运动同步参考反射 / 散射 / 多径回波
结构部分作用仿真对应
波形源 / 时钟中心频率、扫频带宽、chirp 周期、CPIcenter_frequency_ghz / fmcw_* / cpi_rate_hz
发射链 / 天线 / 雷达罩发射能量、方向图、安装与罩体损失link_budget_db / antenna_pattern / rig 外参 / radome_wet_loss_db
Rx 阵列 / MIMO孔径与相位差提供角信息构型、角格、虚拟阵元规模;FAST DOA 来自主散射几何
混频 / ADC / DSP距离、速度、角度、目标检测FAST 功率格或 waveform 相干频域近似,CFAR / NMS / MTI
输出接口检测点、目标列表、诊断/信号数据PCD / NPY / JSON / 可选功率与复数 cube

雷达构型与测量能力

构型结构与测量能力代表型号
imaging_4d
高密成像雷达
大规模或级联 MIMO 阵列,测量距离、径向速度、方位角和俯仰角;以较细角格生成高密度空间点云。arbe_phoenix、ars548、huawei_12t24r
mimo_4d
稀疏 MIMO 雷达
单芯片多发多收阵列,保留方位与俯仰测量;采用较粗角格,适用于近距多目标检测。ti_awr2944、ti_iwr6843_60ghz
planar_3d
平面阵列雷达
水平线阵测量距离、径向速度与方位角;无俯仰测量,输出俯仰角和 Z 坐标为零。ars408
monopulse_mti
单脉冲动目标雷达
少通道宽视场构型,FAST 使用主散射路径几何方向估计 DOA,配合 MTI 径向速度凹口。俯仰输出由具体型号配置决定。zhwx_6130_60ghz、infineon_bgt60tr13c、hlk_ld2450_24ghz
custom
自定义参数构型
按应用配置频率、视场、量程、角格、俯仰开关与信号链参数。用户自定义配置

02 / DETECTION PRINCIPLES基本探测原理

距离来自传播时延与扫频拍频,速度来自相位/多普勒变化,角度来自接收孔径的空间相位差。静止物体也产生回波,MTI 按相对径向速度过滤检测结果。[11]

波长与往返时延[1]式 01
λ = c / fc  ;  τ = 2R / c

c=299 792 458 m/s;60 GHz 的波长约 5 mm,77 GHz 约 3.89 mm。R 为单站雷达到目标的单程距离。

FMCW 距离与径向速度[1]式 02
S = B / Tchirp  ;  fb,range ≈ 2SR / c
vr = (vtarget − vsensor) · û  ;  fD = −2vr / λ

拍频的距离项与多普勒项需分离或联合估计。这里 fD 按接收相对发射的频移约定,输出 vr 远离为正;器件混频约定可能改变频谱符号。

CPI 与分辨率[2]式 03
ΔR = c / (2B)  ;  TCPI = NsTr
Δv = λ / (2NsTr)  ;  Vmax = λ / (4Tr)

B 为带宽;Tr 为 chirp 重复周期;Ns 为慢时间样本数。CPI 决定相干观测时长,输出帧率决定扫描调度频率。

接收阵列测角[2]式 04
Δφ ≈ (2πd / λ) sin α  ;  Δα ≈ λ / (Nd)
d = λ/2  ⇒  Δα ≈ 2/N (rad)

阵元间相位差用于真实器件 DOA;角孔径公式为近轴尺度近似。当前 waveform 用虚拟阵列规模推导角格,FAST 用路径几何角而非逐阵元相位差求解。

分辨率、格宽与精度不同。 增大带宽改善可分离距离;延长 CPI 改善速度格宽;增大孔径改善角分辨尺度。测量精度还取决于 SNR、估计方法和系统误差。

03 / PHYSICAL FEATURES物理仿真特性与公式

各模块使用路径几何、运动状态和非视觉物性计算回波。下表对应各项仿真特性、实现模块与公式。

Feature实现模块 / 抽象物理作用公式
几何传播CPU LineTrace / RTX Inline RayQuery闭合路程、遮挡、到达角[3]05
刚体多普勒状态观测器 / 命中点运动学线运动+旋转+雷达运动[4]06
频变电磁材质39 类预设 / PM 自定义覆写10–120 GHz,幂律 / Debye 先验[8]07
复介质与薄层CPU 散射核 / GPU ShaderFresnel TE/TM / Airy 相干干涉[5]07、08
镜面 / 漫反射表面粗糙度散射核近法向镜面叶瓣+粗糙漫散射[6]09
逆反射材质 retro 属性增强接近返回方向的散射贡献[6]10
射线密度积分立体角微元 / 截面上界Δσ = σ⁰ΔA;按射线功率汇总[7]10
多跳镜面续传SBR + Return Visibility逐跳反射权重与返回可见性[3]05、10
链路 / 天线LinkBudget / AntennaPatternR⁻⁴、方向增益、归一化噪声参考[7]11
FAST 检测近似4D power matrix功率质心;主散射 metadata[9]12
waveform 相干近似FMCWProcessor物理格宽、路径相位、干涉与噪声[10]03、13
CFAR / NMS检测器 / 构型策略局部背景自适应阈值+峰值合并[9]14
MTI / 6130monopulse_mti 构型径向速度凹口,主散射几何 DOA[11]15
天气与湿膜独立 Radar Weather传播等效损失、湿膜物性、检测级响应[14]16–19
测量 / 数据合同Noise / Encoding / WriterISO 坐标、复现种子、同帧目标列表[15]15

传播、材质与回波功率

闭合路径与运动

闭合传播路径[3]式 05
L = Σj=0…K ‖pj+1 − pj‖  ;  Req = L/2
τ = L/c  ;  Φ = 2πL/λ = 4πReq/λ

p0 为 Tx,pK+1 为 Rx;K 为命中次数。多径距离取总闭合路程的一半,到达角取最终散射点到接收机的方向。

刚体命中点速度与多径 Doppler[4]式 06
vj = vcm,j + ωj × (pj − pcm,j)
dL/dt = Σ êj · (vj+1 − vj)  ;  veq = ½ dL/dt

线速度、旋转表面点速度和雷达自身运动共同作用;静止物体仍能反射,是否通过 MTI 取决于相对径向速度。

发射射线逐跳沿镜面反射方向续传。多跳回波需要接收机返回方向的可见性;当前多跳准直判断采用反射方向与返回方向点积大于 0.95 的几何门限。到达角不直接照搬发射角。[3]

复介质、薄层与粗糙表面

复介质 Fresnel 反射[5]式 07
ε = ε′(1 + j tanδ)  ;  q = √(ε − sin²θ)
rTE = (cosθ − q)/(cosθ + q)
rTM = (εcosθ − q)/(εcosθ + q)
Γ = (|rTE|² + |rTM|²)/2

θ 为入射角;Γ 是功率反射率。复数符号采用实现约定;高导电金属采用 PEC 近似,Γ=1。

相干薄层反射[5]式 08
Z = exp(−4πdqi/λ) · exp(j4πdqr/λ)
rslab = r(1 − Z)/(1 − r²Z)

对 TE/TM 分别计算,再平均功率。d 与 λ 使用同一长度单位;厚度控制层内传播相位与吸收。

镜面 / 漫散射与粗糙度[6]式 09
g = 4πσhcosθ/λ  ;  ρs = exp(−g²)
m = clamp(0.08 + 0.12σh,mm, 0.08, 0.60)
σ⁰ = Γ[ρscos⁴θ · exp(−tan²θ/m²) + (1−ρs)cos²θ]

σh 为 RMS 粗糙度;m 为项目散射核的坡度尺度。前项为近法向镜面叶瓣,后项为粗糙漫散射;数值实现含指数截断。

非视觉 PM 资产中的物性独立于渲染材质的 Metallic / Roughness。预设物性按频率合成,用户覆写通过 CPU 材质解析或 GPU 段表映射使用。含水介质采用 Debye 类型色散先验,常规介质按频率幂律缩放。[8]

频变材质先验[8]式 21
ε′(f)=ε′77(f/77)b  ;  σcond(f)=σ77(f/77)d
tanδ(f)=tanδ77(f/77)d−1−b
εDebye(ω)=ε∞+(εs−ε∞)/(1+jωτD)

幂律用于常规介质,Debye 模型用于含水介质;τD 是介质弛豫时间,区别于传播时延 τ。

面积积分、镜面续传与链路

射线微元与镜面续传[7]式 10
ΔA ≈ R²ΔΩ  ;  Δσ = σ⁰ΔA
dout = din − 2(din·n̂)n̂
Wspec = ∏ Γbρs,b

ΔΩ 为单射线立体角;CPU 面积微元还受目标包围盒截面上界约束。多跳续传累乘反射与相干镜面权重;返回回波需可见性和方向约束。

雷达方程与归一化接收功率[7]式 11
Pr = PtGtGrλ²σ / [(4π)³R⁴Lsys]
Pnorm = Δσ · Gant · 10(Clink − 40log₁₀R − Lweather)/10

上行为单站点目标参考雷达方程;下行为项目归一化功率。Clink 即 link_budget_db,汇总发射/接收/噪声参考等常数。

方向图可选 uniform 或 cosine_falloff;后者按离轴角余弦及指数控制方向增益。RCS 使用 m² 的线性量参与功率计算,输出时换算 dBsm。固定目标/阈值条件下 +3 dB 约两倍功率,+12 dB 约两倍理想距离。[7]

FAST、waveform 与检测链

FAST · 快速检测级近似

不保留逐路径相干相位,按功率累加。适合几何/检测级数据与高效扫描;6130 的主散射方向由几何近似提供。

datacube_mode: detection

waveform · 相干频域近似

携带路径传播相位并合成 I/Q,体现相干叠加与相消;带宽、慢时间和阵列参数决定信号格宽。

datacube_mode: waveform

FAST 功率格与质心[9]式 12
Pcell = Σ Pi  ;  x̄ = Σ(Pixi)/ΣPi
x ∈ {R, v, α, ε}

路径进入 (range, velocity, azimuth, elevation) 检测格。功率加权质心输出连续距离、速度和角度;主贡献路径传递对象与材质 metadata。

相干 waveform 与复数噪声[10]式 13
zk = Σ aki√Pi exp(jΦi) + nk
Pk = |zk|² = Ik² + Qk²
n ~ CN(0, N)  ;  E|n|² = N

aki 为离散主旁瓣振幅核;实现采用 0.60 主格功率权重、4×0.09 距离/速度旁格和 2×0.02 方位旁格的核,边界处可能截断。相干合成保留路径干涉。

CA-CFAR 与检测级 SNR[9]式 14
N̂ = (ΣtrainingP)/Ntraining + Nthermal
T = ACFARN̂  ;  Pcluster > T
SNRdB = 10log₁₀(Pcluster/N̂)

训练窗排除 CUT / guard 区域;Acfar 即线性 cfar_threshold。阈值随训练窗背景功率变化,倍率控制检测灵敏度。

RCS 反解与 MTI[11]式 15
σ̂ = Pcluster / [GantGpath(R)]
RCSdBsm = 10log₁₀σ̂
MTI:vnotch > 0 且 |vr| < vnotch 时剔除

RCS 由归一化链路反解;MTI 阈值边界保留。默认 6130 凹口为 0.05 m/s,关闭凹口可保留满足检测阈值的静态回波。

标准构型采用二维 CFAR 和峰值处理;monopulse_mti 使用专用的动目标/主散射方向策略。NMS 聚合邻近响应并保留局部峰值。[9]

恒定天气、湿膜与检测响应

默认晴天和恒定天气状态。改变场景视觉雨雾不会自动改变 Radar;需要在雷达实例显式指定天气参数。物理材质层用于表达湿面响应,检测级层用于表达等效增益、杂波和统计退化;离线工具可只运行后者以生成数据,不必重复 UE 场景几何追踪。

天气等效链路损失[12]式 16
Lweather,dB = 2γR/1000 + Lradome + Lclutter

R 为 m,γ 为单程 dB/km;雷达罩项和杂波项为等效 dB 损失。天气独立于 UDS 视觉效果,默认晴天且状态恒定;随机化是可选应用。

湿膜厚度与连续材质调制[13]式 17
nup = clamp(n̂·ẑworld, 0, 1)
hw = hmaxWnup  ;  fw = hw/(σh,dry+hw)
σh,wet = σh,dry(1−fw)
√ε′wet = (1−fw)√ε′dry + fw√ε′water
tanδwet = (1−fw)tanδdry + fwtanδwater

W∈[0,1],hmax=3 mm;式中长度统一为 mm。当前连续混合控制粗糙度、介电实部和损耗,高导电金属只执行粗糙度调制;不需要为每种基底枚举天气标签。

天气检测级响应[14]式 18
ΔRCSphys = { Gscos⁴θ − Droadnupsin²θ,K=1
(K−1)Gmp,K≥2 }
ΔRCS = ΔRCSphys + wᵀφ
φ = [1, R/100, cosθ, K−1, (1−cosθ)(K−1), clamp(RCS,−40,40)/20]

Gs、Droad、Gmp 和残差权重为可配置 dB 系数。路径散射物性和检测级响应是不同层;离线合成可使用检测级参数模型而无需重跑场景追踪。

天气随机项与检测扰动[14]式 19
Nclutter ~ Poisson(λclutter)
pkeep = (1−pbase) / [1+exp(−(SNR−SNRmin)/s)]
σangle = σradome√[10 / max(10SNR/10,1)]

clutter_density 为每次扫描期望杂波数量;s>0 时启用近阈值软丢点,否则只有基本概率和硬阈值。角扰动与 SNR 相关;固定天气参数下,测量噪声与杂波由随机种子控制。

04 / MODULE PIPELINE仿真模块 pipeline

几何传播提供“路径”,材质和链路提供“能量”,信号链生成“响应”,检测器形成“目标”,输出层保证坐标和数据合同。

仿真 pipeline:配置、扫描、传播、路径回波、FAST 或 waveform、检测、天气测量和输出配置 / 扫描位姿 · 材质 · 射线CPU / RTX 传播命中 · 多跳 · 可见性路径回波记录Req · vr · DOA · PFAST功率格 / 主路径信息waveform相干 I/Q / 复数噪声CFAR / NMS / MTIRCS · 目标几何信息天气 / 测量 / 输出PCD · JSON · cube

配置与型号解析[2]

输入 capture JSON、预设、安装外参、频率、显式内存预算;实例字段覆写预设。

FIACaptureConfigLoader / FIARadarConfig

公式:03

频变物性与天气状态[8]

构建非视觉材质表;PM 自定义参数覆盖相应属性。默认晴天、天气状态恒定;按需配置独立天气。

MaterialLibrary / MaterialClassifier / WeatherModel

公式:07、17

扫描射线与状态观测[4]

按视场与 rays_per_deg 建立射线,获取传感器位姿及目标线/角速度,保存扫描时刻。

Sensor / StateObserver / Rig

公式:06

几何追踪与多跳续传[3]

CPU 查询碰撞,RTX 查询 TLAS;解析命中法线、材质和对象,镜面反射续传,并检查返回可见性。

ExecuteTracesCPU / RTViewExtension / RayQueryCS

公式:05、10

表面响应与路径回波[7]

求散射截面、方向增益、链路功率与路径 Doppler,得到 Req、到达角、功率、物性和 ID。

FPathReturn / 射线回读记录

公式:06–11

FAST / waveform 分支[10]

FAST 将功率非相干累加到 4D 格;waveform 以路径相位合成复数频域响应并加噪。

DetectionProcessor / FMCWProcessor

公式:12、13

检测、峰值与 MTI[9]

CFAR 估计背景,NMS/构型策略聚合峰值,反解 RCS;对径向速度执行 MTI。非 CFAR 路径按直接阈值检测。

RunCFAR / ConvertDetectionsToOutput

公式:14、15

天气响应与测量噪声[14]

按目标信息应用天气检测级修正、杂波/丢点和角扰动,再施加距离、速度和角测量噪声。

WeatherModel / ApplyMeasurementNoise

公式:18、19

坐标、编码与异步写盘[15]

转换到 ISO 输出,记录发射扫描时间与位姿;点云与目标列表配套。cube 在信号链阶段生成,独立输出。

Encode* / FEncodedPayload / Writer

公式:15

数据输出与生命周期

数据生成阶段 / 语义合同
检测点与目标 JSON检测、天气/测量处理后的点与对应目标信息点云与目标列表配套保存,目标列表携带对象、材质和路径几何信息
功率 cubeFAST 或 waveform 信号链响应float32;(NEl, NAz, NV, NR),R 为最快维
复数 cubewaveform I/Q,或 FAST 检测功率赋相位的合成表示complex64;不是硬件 ADC 原始采样,类型与 shape 不替代语义说明
异步 GPU 扫描派发 / GPU 回读 / 编码 / 写盘分离回读绑定派发位姿、帧号和仿真时间;停止采集需排空尾部扫描
双雷达合并按外参把两台传感器检测转换到公共坐标仅变换与拼接;保留各来源径向速度,不自动去重、关联或跟踪

传感器输出为 ISO 8855 右手系(X 前、Y 左、Z 上);UE 局部坐标到输出的换算为 pISO,m=0.01(xUE,−yUE,zUE)。安装 JSON 的位置为 cm,物理核心按 SI 单位使用。[15]

05 / CALIBRATION标定原理

标定使用实测数据估计安装外参、系统链路、材质和天气参数。各类参数使用对应的参考目标与观测量。[16]

外参 / 系统链路

已知对应点求刚体变换 TBS;已知 RCS、距离和 SNR 求系统链路常数。未观测/退化项目保留原配置。

材质 / 天气参数

多入射角单跳观测拟合电磁响应;晴天/天气数据拟合等效退化项。使用物性先验和约束稳定求解,再以独立复采验证应用效果。

标定的先验约束优化[16]式 20
θ̂ = arg minθ { Σ ρ[wi‖yi−hi(θ)‖²] + ‖Λ(θ−θprior)‖² }

外参、链路、材质和天气使用各自观测与目标函数。先验正则化约束参数偏离参考物性的幅度,观测残差衡量模型与实测数据的拟合程度。

观测准备、频率与角度覆盖、参数回写和复采流程见 Radar 标定章节;使用操作见 应用手册。

符号与单位速查
fc / B:中心频率 / 带宽,Hz
λ / R / L:波长 / 单程距离 / 闭合路径长,m
Tr / TCPI:重复周期 / 相干处理时长,s
vr / veq:径向 / 等效径向速度,m/s
α / ε / θ:方位 / 俯仰 / 入射角
σ / σ⁰:RCS(m²)/ 散射截面系数
σh / d:RMS 粗糙度 / 薄层厚度,材质核用 mm
Γ / ρs:功率反射率 / 相干镜面因子
ε′ / tanδ:介电常数实部 / 损耗角正切
K / ΔΩ:命中跳数 / 射线立体角(sr)
Clink / γ:链路常数(dB)/ 比衰减(dB/km)
I / Q / z:正交复响应分量 / 复数响应
附注
  1. 电磁波与 FMCW 测距。基本关系 λ=c/f、τ=2R/c、拍频 fᵦ≈2SR/c(忽略或分离多普勒项)。用于说明真实器件探测原理;仿真 waveform 采用频域路径合成,不逐 ADC 样点实现接收机。
    实现锚点:IARadarFMCW.cpp · 外部原理 / 标准参考
  2. 距离 / 速度 / 阵列分辨率。FRadarFMCWProcessor::Configure 推导 ΔR、Δv、Vmax 与半波长阵列角格。λ/(Nd) 是孔径尺度近似,分辨率不等于精度。
    实现锚点:IARadarFMCW.cpp
  3. 射线传播与闭合多径。GPU Inline RayQuery 与 CPU LineTrace;返回接收机可见性、多跳准直与到达方向。
    实现锚点:IARadarSensor.cpp
  4. 刚体运动学与路径 Doppler。命中点速度 v=vcm+ω×r,闭合路径 v_eq=½dL/dt;输出径向速度采用远离为正。
    实现锚点:IARadarSensor.cpp
  5. Fresnel 与薄层干涉。ComputeSurfaceReflection:复数 TE/TM Fresnel 系数,Airy 薄层反射;高导电金属采用 Γ=1 的 PEC 近似。
    实现锚点:IARadarRayTrace.usf · 外部原理 / 标准参考
  6. 粗糙度、镜面与漫散射。ComputeCoherentSpecularFactor / ComputeBackscatterCoefficient:瑞利粗糙度因子及项目双分量散射核。
    实现锚点:IARadarRayTrace.usf
  7. 射线微元与链路预算。ComputeDifferentialRCSM2 / ComputePathGainLinear:面积微元与归一化接收功率;link_budget_db 汇总系统常数。
    实现锚点:IARadarSensor.cpp
  8. 频变非视觉物性。BuildTableForFrequency:39 类材质,77 GHz 锚点、幂律与含水介质 Debye 色散。外部标准用于先验,不替代逐材质实测。
    实现锚点:IARadarMaterialLibrary.h · 外部原理 / 标准参考
  9. FAST 功率格与 CFAR。FRadarDetectionProcessor:4D 功率/质心累加、range–velocity 训练窗、积分图、CFAR 与构型相关 NMS。
    实现锚点:IARadarDetection.cpp
  10. 频域相干合成与复高斯噪声。BuildDatacube / InjectNoise:路径相位、离散主旁瓣振幅核、I/Q 加噪。不是实际硬件 ADC 原始采样。
    实现锚点:IARadarFMCW.cpp
  11. MTI 与 6130 构型。ConvertDetectionsToOutput 与 monopulse 分箱策略:|vr|实现锚点:IARadarSensor.cpp
  12. 大气 / 雷达罩 / 等效杂波损失。ComputeWeatherLossDb:2γRkm+Lradome+Lclutter。当前频率缩放为项目先验,不是完整雨强/极化求解器。
    实现锚点:IARadarWeatherModel.cpp · 外部原理 / 标准参考
  13. 湿膜连续物性模型。ModulateBounceEMForWetness:法线投影厚度、粗糙峰淹没比例、折射率实部混合和损耗线性插值。hmax=3 mm 为模型常数。
    实现锚点:IARadarWeatherModel.cpp
  14. 天气检测级响应与随机项。ComputeTargetRcsAdjustmentDb / ApplyPostDetectionWeatherEffects:单跳/多跳增益、六维残差、Poisson 杂波、Bernoulli 丢点、SNR 自适应角扰动。
    实现锚点:IARadarWeatherModel.cpp
  15. 坐标、数据输出与扫描生命周期。FIARadarSensor 的 ISO 坐标、异步回读发射位姿、PCD/NPY/JSON 和 cube 编码。文件格式不能代替物理语义。
    实现锚点:IARadarSensor.cpp
  16. 外参、系统与材质标定。IATools registration / em_inversion / weather:刚体配准、已知 RCS/SNR 链路估计、带物性先验的材料拟合、天气参数拟合。优化使用先验是设计选择。
    实现锚点:calibration