01_简单散射体组合拟合阶段总结

UAV 角域 RCS 简单散射体组合拟合阶段报告

日期：2026-05-15

1. 研究目标

本阶段目标是用一组简单、可解释、可在 HFSS 中真实建模的散射体，复现 UAV 放大 10x 后的单站方位角域 RCS 曲线。目标曲线的主要困难在于：-90 deg、0 deg、90 deg 附近存在窄尖峰，同时非峰区域仍有连续背景起伏。

当前采用的结构是：

目标角域 RCS ≈ 三块定向平板主峰 + 二面角/三面角背景补偿

物理分工如下：

• plate：承担 -90 deg、0 deg、90 deg 三个窄主峰；
• dihedral：提供较宽的局部背景和峰肩补偿；
• trihedral：补偿更宽角域的背景起伏。

2. 为什么选择平板作为主峰单元

前期真实 HFSS/PTD 标定表明，三面角反射器主瓣过宽，容易把非峰背景抬高；平板在镜面方向响应极窄，更符合目标曲线中的尖峰特征。

简单散射体单体方向图对比如下：

散射体	俯仰角	峰值角	峰值	半功率宽度	峰值-背景P90
平板	90 deg	0 deg	41.456 dB	1.0 deg	45.283 dB
二面角	90 deg	0 deg	44.440 dB	30.0 deg	8.421 dB
三面角	90 deg	0 deg	36.695 dB	39.0 deg	8.501 dB

结论很直接：

窄尖峰 -> 平板
宽背景 -> 二面角/三面角

3. 响应库组合算法

3.1 从单个散射体到组合阵列的总体流程

优化的核心不是直接在 HFSS 里盲扫所有几何参数，而是先把每种简单散射体的角域响应做成“响应字典”，再从字典中选择少量原子组合成目标曲线。流程如下：

1. 单体标定：分别仿真 plate、dihedral、trihedral 在不同俯仰角下的方位 RCS。
2. 响应库构建：把每条单体响应曲线转成功率域，并保存为一个可平移的基函数。
3. 方位平移：把每个单体响应复制到多个候选中心角，形成大量候选响应原子。
4. 主峰锁定：根据物理判断，强制用 plate el90 负责 -90 deg、0 deg、90 deg 三个主峰。
5. 背景拟合：扣除主峰平板贡献后，只用 dihedral / trihedral 去拟合剩余背景。
6. 稀疏选择：从候选背景原子中保留贡献最大的少量原子，避免阵列过复杂。
7. 尺寸映射：把功率缩放系数映射为真实几何尺寸。
8. HFSS 验证：生成真实 PEC 几何阵列，在 HFSS/PTD 中验证相干散射结果。

对应关系可以写成：

$$\text{single scatterer responses}\to \text{shifted response atoms}\to \text{nonnegative combination}\to \text{HFSS geometry array}$$

3.2 单体响应库

对每个简单散射体进行 HFSS/PTD 单体仿真，得到方位响应曲线：

$$R_i(\varphi )={\mathrm{RCS}}_i(\varphi )[\mathrm{dB}]$$

为了进行线性组合，先转换到功率域：

$$P_i(\varphi )=10^{R_i(\varphi )/10}$$

若该单体被旋转到中心方位 c_i，则使用方位平移近似：

$$A_i(\varphi ;c_i)=P_i\left(\mathrm{wrap}(\varphi -c_i)\right)$$

其中：

$$\mathrm{wrap}(\alpha )=((\alpha +180^{\circ })\mathrm{mod}360^{\circ })-180^{\circ }$$

这样，一个单体响应曲线可以扩展成多个不同方位的候选原子：

$$\mathcal{D}=\left\{A_i(\varphi ;c_k)\mid i\in \text{scatterer library},c_k\in \text{candidate azimuths}\right\}$$

这里的 candidate azimuths 是离散方位网格，例如 -120 deg 到 120 deg 内每隔 15 deg 取一个候选中心角。优化器并不是连续搜索所有角度，而是在这个候选字典中选择组合。

3.3 主峰平板锁定

如果直接让 NNLS 自由选择所有原子，算法可能会选到“数学上能降低误差、但物理职责混乱”的组合，例如让偏俯仰的大平板或弱角反去补某个主峰。为了让结果更稳定，本阶段加入物理约束：

$${\mathcal{P}}_{\mathrm{main}}=\left\{{\text{plate}}_{90^{\circ },-90^{\circ }},{\text{plate}}_{90^{\circ },0^{\circ }},{\text{plate}}_{90^{\circ },90^{\circ }}\right\}$$

三块主平板的中心角固定为目标主峰位置。每块平板的初始功率系数由目标局部峰值与 1 m 单平板库峰值的比值给出：

$$w_k^{\mathrm{plate}}=\frac{{\max }_{\varphi \in {\Omega }_k}{P}_{\mathrm{target}}(\varphi )}{{\max }_{\varphi }{P}_{\mathrm{plate},1m}(\varphi )}$$

其中 \Omega_k 是第 k 个主峰附近的小角域窗口。这样做的含义是：主峰位置先由平板保证，不再让背景原子承担尖峰。

三块主平板的初始贡献为：

$$P_{\mathrm{plates}}(\varphi )=\sum _{k\in \{-90,0,90\}}{w}_k^{\mathrm{plate}}{A}_{\mathrm{plate},90^{\circ }}(\varphi ;k)$$

3.4 背景残差拟合

主峰由平板承担后，背景优化只处理剩余功率：

$$P_{\mathrm{res}}(\varphi )=\max \left(P_{\mathrm{target}}(\varphi )-P_{\mathrm{plates}}(\varphi ),0\right)$$

背景字典只允许使用二面角和三面角：

$$\mathcal{B}=\left\{A_i(\varphi ;c_k)\mid i\in \{\text{dihedral},\text{trihedral}\}\right\}$$

随后用非负最小二乘拟合背景残差：

$$P_{\mathrm{bg}}(\varphi )=\sum _j{v}_j{B}_j(\varphi )$$ $$\underset{\mathbf{v}}{\min }{\Vert W_{\mathrm{bg}}(\varphi )\left(\sum _j{v}_j{B}_j(\varphi )-P_{\mathrm{res}}(\varphi )\right)\Vert }_2$$

约束为：

$$v_j\ge 0$$

最终响应库模型为：

$$P_{\mathrm{fit}}(\varphi )=P_{\mathrm{plates}}(\varphi )+P_{\mathrm{bg}}(\varphi )$$

其中权重函数仍用于强调峰值区和斜率较大的区域：

$$W(\varphi )=\sqrt{1+2{\mathrm{peak}}_{\mathrm{norm}}(\varphi )+{\mathrm{slope}}_{\mathrm{norm}}(\varphi )}$$

3.5 稀疏原子选择

背景候选原子数量较多，如果全部保留，会导致 HFSS 阵列复杂、耦合增强、解释性变差。因此先用完整背景字典求出一组非负系数，再按每个原子的峰值贡献排序：

$$S_j=v_j\underset{\varphi }{\max }{B}_j(\varphi )$$

只保留贡献最大的若干个背景原子，再在这个小集合上重新做一次非负最小二乘。当前保留了 6 个背景原子，加上 3 个主平板，共 9 个散射体。

这个策略相当于一个简单的两步稀疏优化：

$$\text{dense NNLS}\to \text{top-}K\text{ atom selection}\to \text{re-fit selected atoms}$$

3.6 几何尺寸映射

响应库拟合得到功率缩放系数 w_i。一阶近似认为 PEC 平板/角反的峰值功率与特征尺寸四次方相关：

$$P_{\mathrm{peak}}\propto L^4$$

因此把响应库功率系数映射为 HFSS 几何尺寸：

$$L_i=L_{\mathrm{ref}}{w}_i^{1/4}$$

这只是初值映射。真实 HFSS 阵列中存在相干叠加、相位差、遮挡和耦合，因此最终必须回到 HFSS 做真实几何验证。

对于主峰平板，后续会进一步引入尺寸倍率：

$$L_k^{\mathrm{plate},\mathrm{new}}=s_k{L}_k^{\mathrm{plate}}$$

其中 s_k 通过 HFSS 尺寸扫描确定，而不是继续依赖响应库一次映射。

3.7 真实 HFSS 相干验证

响应库阶段是假设：

$$P_{\mathrm{total}}(\varphi )=\sum _i{P}_i(\varphi )$$

真实 HFSS 求解的是复场相干叠加：

$$E_{\mathrm{total}}(\varphi )=\sum _i{E}_i(\varphi )\exp \left(j{\psi }_i(\varphi )\right)$$ $$P_{\mathrm{total}}(\varphi )={\left|E_{\mathrm{total}}(\varphi )\right|}^2$$

其中传播相位近似包含：

$${\psi }_i(\varphi )\approx -\frac{4\pi }{\lambda }\left(x_i\cos \varphi +y_i\sin \varphi \right)+{\psi }_{i0}$$

所以响应库算法用于给出可解释的设计初值，真实效果以 HFSS 阵列仿真为准。

4. 当前阵列设计

当前阵列由 9 个简单散射体组成：三块平板负责主峰，六个二面角/三面角负责背景。

角色	散射体	方位角	俯仰角	尺寸
主峰	平板	-90 deg	90 deg	0.820 m
主峰	平板	0 deg	90 deg	0.614 m
主峰	平板	90 deg	90 deg	0.815 m
背景	二面角	0 deg	90 deg	0.143 m
背景	三面角	30 deg	105 deg	0.179 m
背景	三面角	-30 deg	105 deg	0.169 m
背景	二面角	90 deg	75 deg	0.676 m
背景	二面角	-90 deg	75 deg	0.656 m
背景	三面角	-45 deg	75 deg	0.144 m

阵列三维示意图：

阵列俯视图：

5. 当前 HFSS 验证结果

真实 HFSS/PTD 阵列仿真结果如下。曲线对齐只用于比较形状，不改变“无对齐”指标。

无对齐对比：

整体指标：

方案	MAE	RMSE	最大误差	相似度
三平板主峰 + 背景散射体，均值对齐	5.260 dB	6.695 dB	17.175 dB	90.620 %
三平板主峰 + 背景散射体，无对齐	6.891 dB	8.334 dB	22.139 dB	87.711 %
响应库非相干拟合	5.322 dB	7.377 dB	23.811 dB	90.509 %
纯三面角真实阵列	约 9.68 dB	约 12.11 dB	约 31.00 dB	较低

三主峰检查，均值对齐后：

主峰	目标峰位	HFSS 峰位	角度误差	峰值误差
左峰	-90 deg	-90 deg	0 deg	-15.960 dB
中峰	0 deg	0 deg	0 deg	-6.549 dB
右峰	90 deg	90 deg	0 deg	-16.526 dB

6. 结果分析

当前方案已经解决了最关键的结构问题：

三个主峰的位置全部对齐

这说明“平板负责主峰”的物理建模是正确的。相比纯三面角阵列，当前方案的整体误差明显下降，且背景不再被宽主瓣严重抬高。

当前主要问题变成：

主峰高度不足，尤其 -90 deg 和 +90 deg 两侧峰偏低约 16 dB

这说明尺寸映射公式给出的主平板尺寸偏保守。原因可能包括：

1. 真实有限平板峰值并不严格按 L^4 缩放；
2. 阵列位置导致复场相干干涉；
3. 背景散射体改变局部场分布；
4. 目标 UAV 的尖峰并非理想单平板镜面散射；
5. HFSS/PTD 模板对大尺寸平板和小尺寸补偿体的数值处理存在实际偏差。

因此下一步不应再大改结构，而应做主平板尺寸校准。

7. 当前结论

阶段性结论：

1. 平板是当前最适合承担窄主峰的简单散射体。
2. 二面角/三面角更适合补偿背景和宽项。
3. “三平板主峰 + 背景散射体”已经在真实 HFSS 中取得当前最好结果。
4. 当前误差主要来自主峰幅度不足，而不是峰位错误。
5. 下一步重点是主平板尺寸扫描，然后再重新微调背景项。

8. 下一步计划

固定背景散射体不变，先扫描三块主平板尺寸：

两侧平板倍率: 1.5, 2.0, 2.5, 3.0
中心平板倍率: 1.2, 1.5, 2.0

目标：

• 拉高 -90 deg 和 +90 deg 两侧主峰；
• 适度拉高 0 deg 中心峰；
• 检查平板放大后是否引入额外旁瓣或背景抬升；
• 找到主峰高度和背景控制之间的折中点。

完成主峰尺寸校准后，再重新拟合背景残差：

$$P_{\mathrm{residual}}(\varphi )=\max \left(P_{\mathrm{target}}(\varphi )-P_{\mathrm{plates},\mathrm{HFSS}}(\varphi ),0\right)$$

之后继续用二面角/三面角响应库拟合：

$$P_{\mathrm{background}}(\varphi )=\sum _j{v}_j{B}_j(\varphi )$$

最终形成闭环：

响应库初值 -> HFSS 真实阵列 -> 主峰校准 -> 背景残差拟合 -> HFSS 再验证