几何信息驱动的神经网络算子(GI-KANO):FNO+GNO+KAN 求解 3D PEC 目标电场积分方程
几何信息驱动的 Kolmogorov-Arnold 神经算子(GI-KANO)
📋 基本信息
| 属性 | 内容 |
|---|---|
| 🏷️ 期刊 | IEEE Transactions on Antennas and Propagation(SCI 一区,电磁场顶刊) |
| 📅 发表年份 | 2025年12月(刚发表!) |
| 🔗 DOI | https://doi.org/10.1109/TAP.2025.3611802 |
| 👤 作者 | Yang Liao(廖洋,复旦硕士生), Ji-Ping Deng, Feng Wang(王峰/通讯作者,复旦副教授) |
| 🎯 方向 | 神经算子(Neural Operator)——FNO + GNO + KAN 三层融合用于电磁建模 |
| 🏷️ 类型 | 方法 |
📝 摘要
In this article, a geometry-informed neural operator (GINO) with learnable activation functions is proposed. It incorporates a data-driven approach for solving the electric field integral equation (EFIE) for 3-D PEC targets. The geometry-informed Kolmogorov-Arnold neural operator (GI-KANO) consists of Fourier neural operators (FNOs), graph neural operators (GNOs), and Kolmogorov-Arnold networks (KANs). FNO can learn the function mappings of query points regularly distributed in space, and GNO can map query points to the surface points of a target as the decoder. Moreover, the kernel functions of the GNO are represented by KAN. The activation functions of KAN are learnable. GI-KANO can solve the EFIE in 3-D space… Owing to the super-resolution capabilities of neural operators, models trained on coarse point clouds can achieve comparable results on dense point clouds during testing.
中文翻译:
本文提出了一种具有可学习激活函数的几何信息神经算子(GINO)。它采用数据驱动方法求解三维完美电导体目标的电场积分方程(EFIE)。几何信息Kolmogorov-Arnold神经算子(GI-KANO)由傅里叶神经算子(FNO)、图神经算子(GNO)和Kolmogorov-Arnold网络(KAN)组成。FNO可以学习空间中规则分布的查询点的函数映射,GNO可以作为解码器将查询点映射到目标的表面点。此外,GNO的核函数由KAN表示,KAN的激活函数是可学习的。GI-KANO可以在三维空间中求解EFIE。三维PEC目标的形状由点云表示,每个点代表每个离散网格单元的中心。利用神经算子的超分辨率能力,在粗点云上训练的模型可以在测试时对密点云达到可比的结果。
🎯 问题与动机 (Problem & Motivation)
三维电磁散射问题的 EFIE 求解面临三重重挑战:(1)计算效率——MoM 需要 O(N²) 填充阻抗矩阵 + O(N³) 求解;(2)几何多样性——不同目标形状各异,传统 ML 方法(CNN)要求固定输入尺寸,无法直接处理任意形状的点云;(3)网格密度差异——训练和推理时的网格密度可能不同。
现有神经算子的不足:
- 纯 FNO:只能处理规则网格(结构化数据),无法处理 3D 目标表面上非均匀分布的点云。
- 纯 GNO:可以处理非结构化点云,但缺少 FNO 的频域信息捕获能力和拟合精度。
- MLP 激活函数固定:传统神经算子使用固定的激活函数(如 GELU),限制了模型的表达能力。
核心问题:能否将 FNO 的频域学习能力、GNO 的几何适应能力、KAN 的可学习激活函数三者融合,构建一个同时处理几何多样性和物理精确性的神经算子?
📜 文献脉络 (Literature Context)
2019 PINN [Raissi et al.] → 用 PDE 残差作为损失函数2020 FNO [Li et al.] → 傅里叶域积分算子,处理规则网格的 PDE2020 GNO [Li et al.] → 图神经算子,处理非结构化网格2021 DeepONet [Lu et al.] → 分支-主干网络学习算子2023 GINO [Li et al.] → FNO + GNO 融合,处理大规模 3D PDE2024 KAN [Liu et al.] → 可学习激活函数替代固定激活函数2024 PhiGRL [Shan et al.] → GNN + 物理残差学习求解 3D CFIE2025 GI-KANO(本文)→ FNO + GNO + KAN 三层融合!首个将 KAN 引入 EM 神经算子本文的突破:首次将 FNO、GNO、KAN 三者融合为单一神经算子架构,专门针对 3D EM 散射问题。其中 GINO 骨架来自 Li et al. (2023),但将那篇中的 GNO 核函数从 MLP 替换为 KAN(FastKAN 实现),实现了可学习激活函数。
🔬 方法详解 (Method Deep-Dive)
核心公式
公式1:EFIE 的神经积分方程形式(方法的理论根基)
公式:
逐行解读:
J_S^k= 第 k 次迭代的表面电流(未知函数)W= 前一次迭代的权重系数(可训练)K_1, K_2= 可学习的核函数(由 KAN 表示),分别作用于未变换和变换后的电流T(J_S^k)= 变换网络(由 KAN 表示),在进入 FNO 前对电流做特征变换∫_s ... ds'= 核积分算子——FNO 在傅里叶域高效计算C - n̂×E_i= 常数项和入射场(作为偏置输入)直观理解: 传统 MoM 的每一步迭代是做固定的矩阵-向量乘 Z·u_k。GI-KANO 将其替换为两个并行的可学习积分算子——一个处理原始电流,一个处理变换后的电流——然后求和。关键在于:积分算子中的核函数 K_1, K_2 和变换 T 都由 KAN 表示,它们的激活函数在训练中自学习,而非预设。
公式2:FNO 的傅里叶层(频域学习的核心)
公式:
逐行解读:
v_t= 第 t 层的特征函数F, F^{-1}= 傅里叶变换和逆变换F(K)·F(v_t)= 在频域做逐元素乘法(替代空域的积分运算!)W= 线性变换(可训练的权重矩阵)σ= 非线性激活函数直观理解: FNO 的巧妙之处——空域中的全局积分算子 ∫K(x,y)v(y)dy 在傅里叶域变成了简单的乘法 F(K)·F(v)。这意味着 FNO 可以”看到”整个空间域,而不仅仅是局部邻域。在 GI-KANO 中,FNO 负责处理空间中规则分布的 query points,学习全局的场-源映射关系。
公式3:GNO 的离散图卷积(几何解码器)
公式:
逐行解读:
N(x)= 点 x 的邻居集合(以 x 为中心、半径 r 的球体内的 query points)K(x, y)= 核函数(由 KAN 表示,而非 MLP——这是本文的核心创新!)1/|N(x)|= 平均聚合(mean aggregation)W= 线性变换直观理解: GNO 充当 FNO 的”解码器”——FNO 输出的是规则 query points 上的解,GNO 将这些 query points 映射回目标表面的实际采样点。每个表面点的值由它邻域内 query points 的加权平均决定,权重由 KAN 核函数学习。这使模型能处理任意形状的目标。
算法流程
输入:目标点云(表面离散点的坐标)+ 入射波参数(θ, φ)
第1步:最远点采样 → 统一采样点数(用于批训练)
第2步:FNO 编码器 2a. 在 48×48×48 均匀网格生成 query points 2b. query points + 入射波编码 → 两个并行的 FNO 分支 - 变换分支:KAN(T) 变换 → FNO - 未变换分支:升维 → FNO 2c. FNO 在傅里叶域计算积分算子 → 输出 query points 上的电流
第3步:GNO + KAN 解码器 3a. 对每个目标表面点,以半径 r 搜索 query point 邻居 3b. KAN 核函数 K(x,y) 计算邻居权重 → 加权聚合 → 输出表面电流 3c. KAN 的激活函数为径向基函数(FastKAN,替代 B-spline 加速)
第4步:输出 → 6 通道(J_x, J_y, J_z 的实部和虚部)
训练:loss = ||J_MoM - J_model||_2 / ||J_MoM||_2(相对 L2 误差)GI-KANO 架构示意(三模块融合)
┌─────────────────────────────────────────────────────┐│ GI-KANO 架构 ││ ││ 入射波 (θ,φ) ──→ 编码 ──┐ ││ ├──→ FNO(傅里叶层 ×N) ││ Query Points (48³) ──────┘ │ ││ │ ↓ ││ │ ┌── 变换分支(KAN(T)) ││ │ │ ↓ ││ └──────────────┤ FNO ││ │ ↓ ││ ├── 未变换分支(升维) ││ │ ↓ ││ │ FNO ││ └──────┬──────┘ ││ ↓ ││ ┌── GNO 解码器 ──┐ ││ │ KAN核函数K(x,y)│ ← 可学习激活函数 ││ │ 邻居聚合 │ ││ └───────┬────────┘ ││ ↓ ││ 目标点云 (表面点) ──→ 输出 6ch 电流 (J_re, J_im) │└─────────────────────────────────────────────────────┘关键创新
-
FNO + GNO 双层融合(GINO 骨架):FNO 处理规则 query points(利用傅里叶变换的全局感受野和计算效率)→ GNO 将 query points 映射到任意形状的表面点(利用图的灵活性处理几何多样性)。两者互补:FNO 负责”函数学习”,GNO 负责”几何适配”。
-
KAN 替代 MLP 作为 GNO 核函数:传统 GINO [40] 用 MLP(固定 GELU 激活)作为核函数 → GI-KANO 用 KAN(可学习径向基激活函数)。实验证明 GI-KANO 在所有测试案例上优于 GINO-MLP(误差降低 2-17%),证明可学习激活函数提升了拟合能力。
-
点云 + 神经算子的超分辨率:训练时用 500 个粗采样点,推理时用 2000 个细采样点——误差几乎不增加(表 IV)。这是神经算子的独特优势:由于在连续函数空间上定义,解不依赖于离散分辨率。对大规模问题意味着训练成本可降低 ~10%。
-
端到端跳过矩阵构造:GI-KANO 直接学习”入射波 + 几何 → 表面电流”的映射,完全跳过 MoM 的阻抗矩阵构造(Z 矩阵)。推理时间仅 0.068-0.077 秒(vs MoM 的 2.8-17.2 秒),加速 40-220 倍。
📊 实验与验证 (Experiments & Results)
实验1:未知入射角的泛化
| 模型 | 几何 | 训练误差 | 测试误差 |
|---|---|---|---|
| GINO-MLP | 椭球 | 0.0559 | 0.0597 |
| GI-KANO | 椭球 | 0.0487 | 0.0498 |
| GINO-MLP | 锥台 | 0.0477 | 0.0479 |
| GI-KANO | 锥台 | 0.0300 | 0.0302 |
| GINO-MLP | 四棱柱 | 0.0350 | 0.0355 |
| GI-KANO | 四棱柱 | 0.0265 | 0.0270 |
GI-KANO 在所有几何上全面优于 GINO-MLP,KAN 的贡献使误差降低 2-37%。
实验2:未知几何的泛化
| 模型 | 椭球 | 锥台 | 四棱柱 |
|---|---|---|---|
| GNO | 0.530 | 0.734 | 0.605 |
| GINO-MLP | 0.155 | 0.164 | 0.189 |
| GI-KANO | 0.153 | 0.158 | 0.159 |
纯 GNO 几乎不可用(误差 >50%),FNO 的加入使其降至 ~16%,KAN 进一步改进。
实验3:超分辨率(粗→细)
| 训练点数 | 测试点数 | 锥台误差 | 四棱柱误差 | 椭球误差 |
|---|---|---|---|---|
| 500 | 2000 | 0.1650 | 0.1795 | 0.1602 |
| 1000 | 2000 | 0.1615 | 0.1917 | 0.1575 |
| 2000 | 2000 | 0.1749 | 0.1610 | 0.1703 |
500 点训练的模型在 2000 点测试上误差几乎不变!这验证了神经算子的超分辨率特性。
实验4:计算时间对比
| 方法 | 网格单元 | 矩阵构造(s) | 求解(s) | 总时间(s) |
|---|---|---|---|---|
| MoM | 1640 | 2.467 | 0.345 | 2.812 |
| GI-KANO | 1640 | 0 | 0.068 | 0.068(41×加速) |
| MoM | 2056 | 3.932 | 0.678 | 4.610 |
| GI-KANO | 2056 | 0 | 0.073 | 0.073(63×加速) |
| MoM | 3844 | 13.141 | 4.011 | 17.152 |
| GI-KANO | 3844 | 0 | 0.077 | 0.077(223×加速) |
加速比随未知数增加而增大——GI-KANO 的推理时间几乎不随网格密度增长(0.068→0.077s),这是相对于 MoM 的根本优势。
实验5:导弹目标
| 模型 | 测试误差 |
|---|---|
| GINO-MLP | 0.1884 |
| GI-KANO | 0.0798 |
在复杂导弹几何上 GI-KANO 仍保持良好性能,且 KAN 的优势更加显著(误差降低 58%)。
实验评价
- ✅ 实验设计系统:覆盖角度泛化、几何泛化、超分辨率、效率对比、真实复杂目标
- ✅ 消融实验完整:GNO vs GINO-MLP vs GI-KANO 三级对比
- ✅ 超分辨率验证严谨:三种训练密度 + 三种几何
- ⚠️ 仅与自己的 GINO-MLP 对比,未与 PhiGRL 等同期方法对比
- ⚠️ 1GHz 单一频率,频率泛化能力未测试
🤔 批判性思考 (Critical Analysis)
优点
- 架构设计精巧:FNO(频域全局学习)+ GNO(空域几何适配)+ KAN(可学习激活函数)三层融合,各自负责不同维度的学习任务,互补性强。这不是简单的模块堆砌,而是基于神经积分方程理论的系统设计。
- 超分辨率特性极具工程价值:训练时 500 点、推理时 2000 点误差几乎不变——意味着可以用粗网格训练、细网格推理,大幅降低训练资源消耗。这对你的 FPRFNO 处理大尺寸飞行器网格的场景非常实用。
- 跳过矩阵构造是根本性加速:传统 MoM 的 O(N²) 矩阵构造是耗时大户,GI-KANO 端到端跳过这一步,推理时间几乎 O(1) 级别(0.068-0.077s),而非 MoM 的 O(N²)。
- FastKAN 加速 KAN 训练:原生 KAN 使用 B-spline 激活函数训练极慢,本文使用 FastKAN(径向基函数实现),使 KAN 的训练效率可行。
局限
- 纯数据驱动,无物理信息嵌入:GI-KANO 完全依赖 MoM 生成的数据集训练,没有像 PhiGRL 那样将 CFIE 残差作为输入或损失函数。论文在结论中明确承认这是未来方向(“Further research on embedding physical information in models is necessary”)。
- 训练成本高:1600 样本训练 300 epoch 需要 13-15 小时(RTX 3090)。对于你的飞行器应用,数据生成(MoM 求解)本身就昂贵,再加上训练时间,端到端成本不低。
- 形状泛化有限:未知几何的误差(~15%)远高于已知几何(~3%)。这意味着需要用足够多样化的几何覆盖训练集才能保证泛化——对于飞行器这种特定形状,可能需要大量变体数据。
- 点云采样策略简单:最远点采样 + 固定半径邻居搜索。对于飞行器这种有细长部件(机翼、尾翼)的几何,固定半径可能同时包含太多或太少邻居。
- 实验平台差异:MoM 跑在 CPU 上(FEKO),GI-KANO 跑在 GPU 上——“加速 223 倍”中有一部分来自 CPU vs GPU 本身的速度差,不完全公平。
对我研究的启发(重点!)
- GI-KANO 是 FPRFNO 最直接的参考架构:你的 FPRFNO 也可以用 FNO 处理频域信息 + GNO 处理几何 + KAN 增强非线性表达。特别是 GNO 解码器的设计——如果你的飞行器网格是点云格式,可以直接借鉴 GNO 的邻居聚合方式将频域解映射到表面点。
- 超分辨率对多保真框架的启发:你论文中的”多保真”(multi-fidelity)可以用 GI-KANO 的超分辨率来实现——低精度 MoM 生成粗网格数据用于训练,高精度推理直接用细网格。或者反过来:MoM 的粗网格 + FNO 的超分辨率 = 等效于细网格精度。
- 物理信息缺失是你 FPRFNO 的差异化优势:GI-KANO 的最弱点是”纯数据驱动、无物理嵌入”。你的 FPRFNO 如果有”频率投影残差”作为物理约束(类似于 PhiGRL 的 CFIE 残差输入),那就在这个维度上超越了 GI-KANO。这可能是你的论文创新点的核心定位。
- KAN 可以考虑但非必须:KAN 的改进在本文中是边际的(3-17%),但引入了训练速度的代价。对于你的 FPRFNO,可以把 KAN 作为”锦上添花”的扩展实验,而非核心创新点。
可追问的问题
- GI-KANO 跳过矩阵构造后如何处理 EFIE 的奇异性(singularity)?MoM 中这部分需要精细的奇异积分处理。
- 如果把 PhiGRL 的物理残差输入(R_k = b - Z·u_k)嵌入 GI-KANO 的迭代结构中,能否结合两者的优势——GI-KANO 的架构效率 + PhiGRL 的物理保真?
- 超分辨率的上限在哪?500→2000 点可行,500→20000 点会怎样?
💻 可复现性 (Reproducibility)
| 维度 | 状态 | 说明 |
|---|---|---|
| 代码开源 | ❌ 未公开 | 论文未提供代码仓库 |
| 数据公开 | ❌ 未公开 | MoM 数据由 FEKO 生成,未公开 |
| 文档质量 | 良好 | 架构参数在 Fig.4 中详细标注(层数、通道数、邻居半径),训练超参数完整 |
| 复现难度 | 中-高 | 需要实现 FNO + GNO + KAN 三个模块 + MoM 数据生成流水线 |
🔗 关联笔记
方法相关
- Fourier_Neural_Operator — FNO:傅里叶神经算子,频域学习(GI-KANO 的核心编码器)
- Graph_Neural_Operator — GNO:图神经算子,几何解码器
- KAN_Learnable_Activation — KAN:可学习激活函数(GI-KANO 中替代 MLP 核函数)
- Operator_Learning — 算子学习:学习函数空间之间的映射
- Super_Resolution — 超分辨率:粗网格训练 → 细网格推理
应用相关
- Point_Cloud_Geometry — 点云几何表示(3D 目标的输入格式)
- Integral_Equation_Solver — 积分方程求解器(GI-KANO 的目标应用)
引用论文
- Li_2020_FNO — FNO 原始论文(GI-KANO 的 FNO 模块来源)
- Li_2023_GINO — GINO 原始论文(GI-KANO 的架构骨架来源)
- Liu_2024_KAN — KAN 原始论文(GI-KANO 的可学习激活函数来源)
- 2024-IEEE_TAP-物理信息图残差学习求解CFIE — PhiGRL(同一应用场景的对比方法,物理信息学习路线)
文章分享
如果这篇文章对你有帮助,欢迎分享给更多人!