草稿:泛函连接理论

泛函连接理论(TFC)是为执行泛函插值而开发的数学框架。该框架提供了一种将约束优化问题转换为等效的无约束问题的方法。通过利用这种转换,TFC可以被有效地应用于微分方程求解问题。但是泛函插值到底意味着什么呢?

从函数插值到泛函插值

在介绍TFC问题之前,首先考虑一个涉及n个约束且具有普遍性的插值问题,例如微分方程的变值问题(BVP)求解。这些约束可能是一致的,也可能是不一致的。例如,在域   的问题中,约束   不一致,因为它们在   处产生不同的值。

如果 n 个约束一致,则可以通过选择 n 个线性无关的支撑函数(例如单项式   )来构造插入这些约束的函数。所选的支撑函数(support function)集可能与给定的约束一致,也可能不一致。例如,约束    与支撑函数   不一致。如果支撑函数与约束一致,则可以求解泛函插值问题,从而产生一个插值函数——一个满足所有约束的函数。选择不同的支撑函数将导致不同的插值。当解决了一个插值问题并确定了一个初始插值函数时,原则上,可以通过使用与约束一致的每一组不同的线性独立支撑函数进行插值过程来生成所有可能的插值函数。但是,这种方法是不切实际的,因为可能的支撑函数集的数量是无限的。

随着泛函连接理论(TFC)的发展,这个问题得以解决。泛函连接理论是德克萨斯A&M大学的D. Mortari教授在其著作中引入的用于执行泛函插值的分析框架[1][2]。该方法构造了一个约束泛函(constrained functional) ,其中自由函数   无论取何种表达式,  都满足给定约束。因此约束泛函   提供了所有可能的插值的完整表示。通过改变  ,可以生成整个插值集,包括那些不连续或部分定义的插值。

 
从函数插值到泛函插值的流程图

图1流程图说明了从函数插值到泛函插值的过程。通过泛函连接理论求解得到的泛函封装了本身满足施加约束的函数的子空间,从而将搜索空间缩小到满足约束的区域。通过,约束优化问题可以重新表述为无约束问题,从而通过更简单、更稳健、更准确、更高效和可靠的方法解决它们。因此,TFC框架将约束问题转化为不受约束的问题,从而显著简化了求解过程。

论文[2]展示了如何解决涉及点、导数、积分以及这些的任意线性组合的单变量约束。然后,该理论被扩展以适应无限和多元约束,并应用于求解常微分方程、偏微分方程和积分-微分方程。TFC 的单变量版本可以用以下两种形式之一表示。

 

式中,n 表示线性约束的数量;  是自由函数; 是 n 个自定义的线性无关的支撑函数;  是系数泛函,  是开关泛函(在各自的约束下计算时取值1,在其他约束下取值0),   是用自由函数表示约束的投影泛函。这些元素的详细说明和性质在参考文献[3][4][5]中提供。

一个例子

为了说明TFC如何泛化插值,在下文中将介绍一个算例。考虑约束条件    ,满足这些约束的插值函数是

 

可以很容易地验证   满足约束条件。在此基础上,构造函数 。

 

同理,不论 取为何种形式, 的导数在   处和   处为0。因此,通过将    相加,得到约束泛函   如下。

 

对于约束泛函   ,无论   取为何种形式,只要     处有定义,必然满足约束条件。由于   的灵活性,约束泛函   就可以在指定的约束以外任意取值。重要的是,这种灵活性不仅限于本例中选择的特定约束。相反,它普遍适用于任何一组约束。这种普遍性说明了 TFC 如何进行泛函插值:它构造一个满足给定约束的函数,同时通过改变   使其在约束之外的地方足够自由。本质上,这个例子表明约束泛函捕获了所有满足给定约束的可能函数,展示了 TFC 在处理各种插值问题方面的强大功能和通用性。

泛函连接理论的应用

自文献[2]发表以来,它在许多领域得到了广泛应用。这些应用包括它在剪切型和混合导数问题中的应用 [6]、分数算子的分析[7][8] 、在弯曲空间中确定边值问题的测地线[9] 以及对延拓法的贡献[10] 。此外,TFC 已应用于间接最优控制[11]、刚性化学动力学建模[12]和流行病学动力学研究。它还在非线性规划 和结构力学等领域显示出潜力。

特别地,TFC 在神经网络中表现出了非凡的效率 ,尤其是在提高准确性和解决高维问题方面。TFC 通过有效消除优化过程中的约束,显著提高了物理信息神经网络 (PINN) 的性能,这是传统神经网络通常难以解决的。此功能显著提高了计算效率和准确性,从而能够更轻松地解决复杂问题。

与谱方法的区别

乍一看,TFC 和谱方法在解决约束优化问题的方法上似乎很相似。但是,它们之间有两个基本区别:

1. 解的表示:谱法将解表示为基函数的和,而 TFC 将自由函数表示为基函数的和。这种区别使 TFC 能够在分析上满足约束条件,而频谱方法将约束视为附加数据,以取决于残差的准确度对其进行近似。

2. 变值问题中的计算方法:在线性变值问题中,两种方法的计算策略差异很大。谱方法通常采用迭代技术(例如打靶法)将变值问题重新表述为初值问题。相反,TFC 通过线性最小二乘技术直接解决这些问题,无需迭代过程。

这两种方法都可以使用Galerkin方法(确保残差向量与所选基函数正交)或Collocation方法(最小化残差向量的范数)来执行优化。

与拉格朗日乘子法的区别

在存在约束的优化问题中。该方法引入了其他变量,称为拉格朗日乘子,必须计算这些变量才能强制执行约束。在某些情形下,拉格朗日乘子法的计算较为困难。相比之下,TFC 不添加新变量,并且能够推导约束泛函,从而简化了问题的求解。但是,需要注意的是,拉格朗日乘子方法具有处理不等式约束的优势,而TFC目前缺乏这种能力。

这两种方法的一个共同局限性是它们容易产生局部最优解而不是全局最优解,尤其是在非凸问题的求解中。因此,可能需要补充验证程序或替代方法来评估和确认所获得解决方案的质量和全局有效性。总之,虽然 TFC 并不能完全取代拉格朗日乘子方法,但在乘数计算变得过于复杂或不可行的情况下,只要约束仅限于等式约束,它就可以作为一个强大的替代方案。

  1. ^ Daniele, Mortari. Least-Squares Solution of Linear Differential Equations. Mathematics. 2017, 5 (4). doi:10.3390/math5040048. 
  2. ^ 2.0 2.1 2.2 Daniele, Mortari. The Theory of Connections: Connecting Points. Mathematics. 2017, 5 (4). doi:10.3390/math5040057. 
  3. ^ Hunter, Johnston. The Theory of Functional Connections: A Journey from Theory to Application, PhD thesis, Texas A&M University, August 2021
  4. ^ Carl, Leake; Hunter, Johnston; Daniele, Mortari. The Theory of Functional Connections: A Functional Interpolation. Framework with Applications, Lulu, Morrisville NC, 2022.
  5. ^ Carl, Leake. The Multivariate Theory of Functional Connections: An n-Dimensional Constraint Embedding Technique Applied to Partial Differential Equations, PhD thesis, Texas A&M University, August 2021.
  6. ^ Daniele, Mortari. Theory of Functional Connections Subject to Shear-Type and Mixed Derivatives. Mathematics. 2022, 10 (24). doi:10.3390/math10244692. 
  7. ^ Daniele, Mortari. Representation of Fractional Operators Using the Theory of Functional Connections. Mathematics 2023. 2023, 11 (23). doi:10.3390/math11234772. 
  8. ^ Daniele, Mortari; Roberto, Garrappa; Luigi, Nicolò. Theory of functional connections extended to fractional operators. Mathematics. 2023, 11 (7). doi:10.3390/math11071721. 
  9. ^ Daniele, Mortari. Using the Theory of Functional Connections to Solve Boundary Value Geodesic Problems. Mathematical and Computational Applications. 2022, 27 (4). doi:10.3390/mca27040064. 
  10. ^ Wang, Yang; Topputo, Francesco. A TFC-based homotopy continuation algorithm with application to dynamics and control problems. Journal of Computational and Applied Mathematics. 2022, 401. doi:10.1016/j.cam.2021.113777. 
  11. ^ Andrea, D’Ambrosio; Enrico, Schiassi; Hunter, Johnston; Fabio, Curti; Daniele, Mortari; Roberto, Furfaro. Time-energy optimal landing on planetary bodies via theory of functional connections. Advances in Space Research. 2022, 69 (12): 4198-4220. doi:10.1016/j.asr.2022.04.009. 
  12. ^ Mario, De Florio; Enrico, Schiassi; Roberto, Furfaro. Physics-informed neural networks and functional interpolation for stiff chemical kinetics. Chaos: An Interdisciplinary Journal of Nonlinear Science. 2022, 32 (6). doi:10.1063/5.0086649.