6.4 高频渐近、混合算法与仿真工程

6.4 高频渐近、混合算法与仿真工程

前三节的方法都在"严格求解"的旗帜下作业，但一艘万吨舰船在十吉赫下的散射，严格方法连起步都困难——未知量以亿计。本节请出值班室的效率组：几何光学与物理光学放弃逐波长求解，承认"波长足够短时波的行为像光"，用渐近近似直接给出答案；弹跳射线追踪把光学近似工程化；混合算法让严格方法管小细节、渐近方法管大场面（旧版混合算法与多尺度内容并入）；最后交代仿真软件的工程学——网格、收敛、并行与验证（旧版软件架构内容并入）。学完本节，你将拥有一张完整的"问题到方法"选型地图，并能识破仿真结果里最常见的三类假货。

一、高频渐近：用光学思维换速度

几何光学（GO）只用直射、反射、绕射（经 GTD 修正）三条规则追迹射线，完全不剖分场——计算量与频率无关。物理光学（PO）再进一步：假设照明区的表面感应电流就是局部入射场决定的那一个（消除暗区电流），对表面电流做辐射积分得到散射场。物理光学在镜面方向与后向散射主瓣内精度可观，是雷达截面快速估算的主力；它的盲区恰是 5.2 节绕射理论的辖区——棱边与爬行波需要 PTD（物理光学绕射修正）补账。

弹跳射线（SBR）把光学近似装配成流水线：从发射口面发出上万条射线，每条按几何光学在目标表面弹跳、追踪幅度相位与极化，末端用物理光学积分合成远场。频率升高，射线数不必增加——渐近方法越高频越便宜，与严格方法的账单走势正好相反。整机隐身评估、汽车电磁兼容、射电望远镜面型公差分析，都是 SBR 的日常工位（5.3 节隐身预算表正是用它们算出来的）。

# 平板物理光学区 RCS 快速估算与方向图
import numpy as np

A = 1.0                                  # 一平方米金属平板
lam = 0.03                               # 10 GHz
theta = np.radians(np.linspace(-60, 60, 241))
k = 2*np.pi/lam
ax, ay = 1.0, 1.0                        # 方板边长
sinc_x = np.sinc(ax*np.sin(theta)/lam)   # np.sinc 自带 pi 归一
sigma = 4*np.pi*(A**2/lam**2)*(np.cos(theta)**2)*sinc_x**2
peak_db = 10*np.log10(4*np.pi*A**2/lam**2)
print(f"法向峰值 RCS {peak_db:.1f} dBsm（平板公式 4πA平方除λ平方）")
half = theta[np.where(10*np.log10(sigma/10**(peak_db/10)) >= -3)[0]]
print(f"主瓣半功率宽度约 {np.degrees(half[-1]-half[0]):.2f} 度——镜面闪光的致命窄角")
​

一平方米平板的法向回波高达 41 dBsm，但主瓣只有几度宽——5.3 节"外形斜置把闪光调离威胁方向"的数值依据，就是这条 sinc 平方曲线。

二、混合算法与多尺度交接

真实工程问题往往尺度混杂：卫星天线（厘米级细节）装在两米平台上，基站阵列贴着高楼（百米级环境）。单一方法顾此失彼，混合算法的纪律于是成了值班室的日常：MoM 管天线细节、PO 管平台大面、UTD 管棱边绕射，交界处用等效原理交换等效电流（6.3 节 FEM-MoM 混合同理）；区域分解让 FEM 或 FDTD 在局部子域各自加密，子域间以边界交换交接。交接纪律三句话：交界面选在场变化平缓处（别切在缝隙上）；两侧网格密度匹配（短板效应）；每级近似都要有一次严格方法校验（ PO 结果定期拿 MoM 抽查）。多尺度问题的另一半功夫在降阶模型：把反复调用的子结构（连接器、过孔）预先算好等效模型存库，主仿真调用查表——工程仿真的快，一半来自肯建库。

三、仿真工程学：从算法到生产力

大型电磁软件是把算法装配成整车：前处理（几何清理、剖分）、求解器（多算法混编）、后处理（近远场变换、参数提取）。工程纪律比软件名牌更重要，三条铁律值得刻在值班室门上。其一，网格即精度：剖分准则（每波长若干网格、金属边缘加密）、共形网格与阶梯网格的差别、边界（PML）的放置距离，每一项都直接写进误差预算。其二，收敛验证是道德：加密网格重跑、换算法抽查、与解析特例对表（平板、球、标准增益喇叭有现成解析或基准数据），三选一也必须做——仿真软件从不保证你设的对，它只保证按你设的算。其三，并行与账单：矩量法内存墙、FEM 的稀疏求解器并行度、FDTD 的空间剖分并行，选型时把集群规格一起考虑；一次整机仿真的电费与机时，是方法选型表里的真实条目。

四、假结果鉴定：三类最常见的仿真赝品

值班室出过最贵的冤案，不是算不出来，而是把赝品签收出库。三类高发假货值得列入交接台账。第一类：未收敛的漂亮方向图——网格加密一档后副瓣电平挪了三分贝以上，说明第一次的结果还在离散误差里漂着；赝品往往图面光洁、数值圆整，恰因误差互相抵消得好看。第二类：吸收边界反射造出的假谐振——腔体或天线 Q 值被 PML 的回波垫高，谐振频率出现成对微小劈裂；把 PML 拉远十格重跑，谐振点搬家即露馅。第三类：端口与馈电的假阻抗——端口模式设错（高阶模混入）或归一化阻抗与系统不符，S 参数曲线整体平移，方向图却可能毫发无损，于是最容易被签收。鉴定纪律共通：每一个关键数都要有第二来源——解析特例、另一算法、或实测基准，三者至少占一样。

基准考题不是入门练习，而是每次大版本网格策略变更后的例行晨检——半小时的对表，换来的是整机仿真数字敢写进报告的底气。

本节要点回顾

• 渐近交换：几何光学与物理光学以光学近似换速度，计算量与频率无关、越高频越便宜。
• PO 的边界：镜面与主瓣区可信，棱边绕射要 PTD 补账，暗区电流假设是其软肋。
• SBR 流水线：射线弹跳加物理光学积分，整机级隐身与兼容分析的主力引擎。
• 混合纪律：交界面选平缓处、网格密度匹配、近似结果定期用严格方法抽查。
• 仿真三铁律：网格即精度、收敛验证是道德、并行账单进选型——软件只保证按你设的算。

值班室六日班次到此排满。下一章换换口味：电波的出勤路上还有天气与意外——噪声、随机介质、非线性互调、混沌与随机共振，统计无线电物理的登场。