几何信息驱动的神经网络算子(GI-KANO):FNO+GNO+KAN 求解 3D PEC 目标电场积分方程

4031 字
20 分钟
几何信息驱动的神经网络算子(GI-KANO):FNO+GNO+KAN 求解 3D PEC 目标电场积分方程

几何信息驱动的 Kolmogorov-Arnold 神经算子(GI-KANO)#

📋 基本信息#

属性内容
🏷️ 期刊IEEE Transactions on Antennas and Propagation(SCI 一区,电磁场顶刊)
📅 发表年份2025年12月(刚发表!)
🔗 DOIhttps://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)网格密度差异——训练和推理时的网格密度可能不同。

现有神经算子的不足:

  1. 纯 FNO:只能处理规则网格(结构化数据),无法处理 3D 目标表面上非均匀分布的点云。
  2. 纯 GNO:可以处理非结构化点云,但缺少 FNO 的频域信息捕获能力和拟合精度。
  3. MLP 激活函数固定:传统神经算子使用固定的激活函数(如 GELU),限制了模型的表达能力。

核心问题:能否将 FNO 的频域学习能力、GNO 的几何适应能力、KAN 的可学习激活函数三者融合,构建一个同时处理几何多样性和物理精确性的神经算子?

📜 文献脉络 (Literature Context)#

2019 PINN [Raissi et al.] → 用 PDE 残差作为损失函数
2020 FNO [Li et al.] → 傅里叶域积分算子,处理规则网格的 PDE
2020 GNO [Li et al.] → 图神经算子,处理非结构化网格
2021 DeepONet [Lu et al.] → 分支-主干网络学习算子
2023 GINO [Li et al.] → FNO + GNO 融合,处理大规模 3D PDE
2024 KAN [Liu et al.] → 可学习激活函数替代固定激活函数
2024 PhiGRL [Shan et al.] → GNN + 物理残差学习求解 3D CFIE
2025 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 的神经积分方程形式(方法的理论根基)

公式: JSk+1=WJSk+sJSkK1(r,r)ds+sT(JSk)K2(r,r)ds+Cn^×Ei\mathbf{J}_S^{k+1} = \mathbf{W}\mathbf{J}_S^k + \int_s \mathbf{J}_S^k K_1(\mathbf{r},\mathbf{r}') ds' + \int_s T(\mathbf{J}_S^k) K_2(\mathbf{r},\mathbf{r}') ds' + \mathbf{C} - \widehat{n} \times \mathbf{E}_i

逐行解读:

  • 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 的傅里叶层(频域学习的核心)

公式: vt+1=σ(Wvt+F1(F(K)F(vt)))v_{t+1} = \sigma\left(W v_t + \mathcal{F}^{-1}\left(\mathcal{F}(K) \cdot \mathcal{F}(v_t)\right)\right)

逐行解读:

  • 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 的离散图卷积(几何解码器)

公式: vt+1(x)σ(Wvt(x)+1N(x)yN(x)K(x,y)vt(x))v_{t+1}(x) \approx \sigma\left(W v_t(x) + \frac{1}{|N(x)|} \sum_{y \in N(x)} K(x, y) v_t(x)\right)

逐行解读:

  • 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) │
└─────────────────────────────────────────────────────┘

关键创新#

  1. FNO + GNO 双层融合(GINO 骨架):FNO 处理规则 query points(利用傅里叶变换的全局感受野和计算效率)→ GNO 将 query points 映射到任意形状的表面点(利用图的灵活性处理几何多样性)。两者互补:FNO 负责”函数学习”,GNO 负责”几何适配”。

  2. KAN 替代 MLP 作为 GNO 核函数:传统 GINO [40] 用 MLP(固定 GELU 激活)作为核函数 → GI-KANO 用 KAN(可学习径向基激活函数)。实验证明 GI-KANO 在所有测试案例上优于 GINO-MLP(误差降低 2-17%),证明可学习激活函数提升了拟合能力。

  3. 点云 + 神经算子的超分辨率:训练时用 500 个粗采样点,推理时用 2000 个细采样点——误差几乎不增加(表 IV)。这是神经算子的独特优势:由于在连续函数空间上定义,解不依赖于离散分辨率。对大规模问题意味着训练成本可降低 ~10%。

  4. 端到端跳过矩阵构造:GI-KANO 直接学习”入射波 + 几何 → 表面电流”的映射,完全跳过 MoM 的阻抗矩阵构造(Z 矩阵)。推理时间仅 0.068-0.077 秒(vs MoM 的 2.8-17.2 秒),加速 40-220 倍。

📊 实验与验证 (Experiments & Results)#

实验1:未知入射角的泛化#

模型几何训练误差测试误差
GINO-MLP椭球0.05590.0597
GI-KANO椭球0.04870.0498
GINO-MLP锥台0.04770.0479
GI-KANO锥台0.03000.0302
GINO-MLP四棱柱0.03500.0355
GI-KANO四棱柱0.02650.0270

GI-KANO 在所有几何上全面优于 GINO-MLP,KAN 的贡献使误差降低 2-37%。

实验2:未知几何的泛化#

模型椭球锥台四棱柱
GNO0.5300.7340.605
GINO-MLP0.1550.1640.189
GI-KANO0.1530.1580.159

纯 GNO 几乎不可用(误差 >50%),FNO 的加入使其降至 ~16%,KAN 进一步改进。

实验3:超分辨率(粗→细)#

训练点数测试点数锥台误差四棱柱误差椭球误差
50020000.16500.17950.1602
100020000.16150.19170.1575
200020000.17490.16100.1703

500 点训练的模型在 2000 点测试上误差几乎不变!这验证了神经算子的超分辨率特性。

实验4:计算时间对比#

方法网格单元矩阵构造(s)求解(s)总时间(s)
MoM16402.4670.3452.812
GI-KANO164000.0680.068(41×加速)
MoM20563.9320.6784.610
GI-KANO205600.0730.073(63×加速)
MoM384413.1414.01117.152
GI-KANO384400.0770.077(223×加速)

加速比随未知数增加而增大——GI-KANO 的推理时间几乎不随网格密度增长(0.068→0.077s),这是相对于 MoM 的根本优势。

实验5:导弹目标#

模型测试误差
GINO-MLP0.1884
GI-KANO0.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(同一应用场景的对比方法,物理信息学习路线)

文章分享

如果这篇文章对你有帮助,欢迎分享给更多人!

几何信息驱动的神经网络算子(GI-KANO):FNO+GNO+KAN 求解 3D PEC 目标电场积分方程
https://sciml.com.cn/posts/gi-kano-neural-operator-efie/
作者
星飞帆
发布于
2026-07-06
许可协议
CC BY-NC-SA 4.0
相关文章 智能推荐
1
神经算子:从 DeepONet 到 FNO
论文精读 神经算子(Neural Operators)的核心思想——学习无限维函数空间之间的映射,覆盖 DeepONet 和 Fourier Neural Operator 两种主流架构。
2
利用物理信息图残差学习求解组合场积分方程用于三维PEC目标电磁散射
论文精读 本研究提出物理信息图残差学习(PhiGRL),作为一种有效且鲁棒的深度学习方法用于三维电磁建模。PhiGRL从物理信息监督残差学习(PhiSRL)扩展而来,通过应用图神经网络(GNN)预测修正量,模拟不动点迭代法的计算过程以迭代修改候选解直至收敛。GNN的应用使PhiGRL能够自适应处理三维电...
3
基于分块的克雷洛夫子空间基函数加速压缩感知矩量法求解双站散射
论文精读 为提高矩量法中结合压缩感知和Krylov子空间求解电磁散射问题时构造基函数的效率,本文提出了一种新方法,采用分块技术构造Krylov子空间基函数。首先将目标划分为若干小块,每块进行扩展以保证电流的连续性。然后在每块上计算Krylov子空间基函数。进一步通过奇异值分解技术增强基函数的正交性以获得...
4
基于物理启发式机器学习与实验设计的导电目标高效 RCS 预测
论文精读 本文提出了一种结合机器学习和实验设计的混合方法,用于高效准确地预测导电目标随入射角变化的单站雷达散射截面(RCS)。该方法称为物理光学启发支持向量回归(POI-SVR),其核函数的设计灵感来源于物理光学法。引入均匀设计(UD)和均匀设计采样(UDS)来获得具有高度代表性的训练样本。对简单和复杂...
5
基于动态计算方法的倾转旋翼机雷达散射截面特性研究
论文精读 为研究倾转旋翼机的雷达散射截面(RCS)特性,本文提出了一种考虑旋翼旋转和短舱倾转的动态计算方法。采用物理光学法和物理绕射理论处理目标的瞬时电磁散射。分析了飞机在直升机模式、固定翼模式和过渡模式下的RCS。结果表明,在固定翼模式下,桨叶对水平面内迎面入射波的偏转效应较弱。直升机模式改善了旋翼在...
随机文章 随机推荐
Profile Image of the Author
星飞帆
西北工业大学24级研究生 | SciML · PINNs · 神经算子 · PDE 约束学习
欢迎
欢迎来到星飞帆的 SciML 研究笔记:这里记录论文阅读、数学推导、代码复现、实验踩坑,也分享研究生阶段可能用得上的资料。
分类
标签
站点统计
文章
21
分类
4
标签
69
总字数
70,675
运行时长
0
最后活动
0 天前

文章目录