凭一篇博士论文留名数学史｜纪念杨大卫百年诞辰

luyuanhong

凭一篇博士论文留名数学史｜纪念杨大卫百年诞辰

今年是美国计算数学家 David M. Young Jr.（杨大卫）诞辰一百周年。他在博士学位论文中提出的“逐次超松弛迭代法”，在计算机求解大型线性方程组方面发挥着重要作用，成为留名数学史的杰出工作。

撰文 | 丁玖（美国南密西西比大学数学系教授）

光看标题，读者可能以为杨大卫是一位华人学者。其实他是纯粹的美国人，全名为 David M. Young Jr. (1923 年 10 月 20 日 - 2008 年 12 月 21 日)。他不像德国“大卫”—— David Hilbert (1862-1943) 那样蜚声国际数学界，我甚至查不到他全名的标准中文翻译；再加上他的姓在西方人中罕见地只有一个元音，我就顺势给他起了个响亮的中文名字：“杨大卫”。

杨大卫生前在计算数学界是个响当当的人物。他毕生致力于求解线性方程组的迭代法，最伟大的成就是发明了 SOR 方法，即逐次超松弛迭代法（successive over-relaxation method）。更了不起的是，这项与时俱进的发明脱胎于他在 1950 年完成的博士论文。按照抽象代数名家和教学名师丁石孙 (1927-2019) 校长所言，数学博士论文“有百分之九十大几的比例对所在学科没有影响”。但是，杨大卫的博士论文影响了整个线性迭代法。1984 年，我在南京大学读到了丁校长的采访记录，感受到他希望研究生们多读书以拓宽知识面，而不囿于一篇学位论文的苦心。

杨大卫于 1944 年在纽约州的韦伯海军建筑学院获得学士学位，为海军工作到二战结束，其后进入哈佛大学数学系读研究生，分别于 1947 年和 1950 年获得硕士和博士学位。他的博士论文导师是加勒特·伯克霍夫 (Garrett Birkhoff，1911-1996)。加勒特和父亲乔治 (George David Birkhoff，1884-1944) 是少见的数学家父子，而且都是哈佛教授，父亲在动力系统和遍历理论领域有开创性贡献，儿子的成就主要集中在格论 (lattice theory)。

格论是抽象代数的一个分支，起始于布尔代数，研究的对象是有序集，小伯克霍夫在老伯克霍夫任教的哈佛大学数学系本科毕业后，去了剑桥大学，先学数学物理，后转向抽象代数。在游学欧洲的那年，他从几大数学名家那里获益良多。1933 年后，他一直待在哈佛。期间发表的一系列论文，构成了他 1940 年出版的经典著作《格论》，其第三版至今还在印刷。二战中，他的兴趣转向应用数学，尤其是流体力学，研究成果包含在他于 1950 年出版的专著《流体动力学》中。由于和冯·诺伊曼是好朋友，他在三十五岁后对计算的爱好与日俱增。

这就是 1948 年杨大卫投入伯克霍夫教授门下攻读博士时的背景。导师给了弟子一个题目：研究数值求解泊松偏微分方程。正是对这一博士学位论题的研究推动了 SOR 方法的问世。之后若干年，在杨大卫研究工作的基础上，伯克霍夫与 1950 年进入哈佛数学系读研究生的瓦尔加 (Richard Steven Varga，1928-2022) 合作，研究了用于微分方程和正算子的迭代法。瓦尔加于 1962 年出版的《矩阵迭代分析》成了这一领域的经典著作。

杨大卫在他 1990 年撰写的综述性文章《迭代方法的历史回顾》（ A historical review of iterative methods）中回忆，当他找伯克霍夫教授选博士论文题目时，心中的头选是纯数学领域的李代数，但是教授建议他做计算数学的“松弛方法”，并递给他几篇文献，其中包括论文《拉普拉斯和泊松方程的数值解》（Numerical solution of Laplace's and Poisson's equations）和英国数学家索斯韦尔爵士 (Sir Richard V. Southwell，1888-1970) 的著作《理论物理中的松弛方法》（ Relaxation Methods in Theoretical Physics）。

所谓“松弛方法”，从广义上讲是指获得偏微分方程近似数值解的过程，但从狭义上讲，它仅指用于求解线性代数方程组Ax = b的迭代过程。在此过程中，“松弛”技术与方程组的残量 r = b - Ax 密切相关，目的是加快收敛。

伯克霍夫对杨大卫的博士论文研究帮助很大，他不仅提供了指导和鼓励、给出建议的参考文献，还仔细阅读了论文草稿。据杨大卫回忆， “逐次超松弛方法”的英文全称就是伯克霍夫建议的——“我相信这是一个很好的选择”。这样的导师堪称良师益友。

杨大卫的博士论文研究开局并不顺利，甚至来访的索斯韦尔爵士也留下一句令人沮丧的评述“任何机械化松弛方法的尝试都是浪费时间（Any attempt to mechanize relaxation methods would be a waste of time）。”但他毫不气馁，继续工作。不久他有了一个好发现：对某些线性方程组，高斯-赛德尔迭代矩阵的特征值是雅可比迭代矩阵特征值的平方。得益于从阅读相关文章获得的灵感，杨大卫有了突破性的进展：对于他通过引进超松弛因子 ω 而设计出的 SOR 方法，如果方程组系数矩阵 A 一致有序 (consistently ordered)，则其迭代矩阵的特征值与雅可比方法迭代矩阵的特征值有个关键性的关系。他所数值求解的椭圆型偏微分方程在区域网格按从左到右和向上的自然顺序编号，使用标准五点差分格式，会产生一致有序的矩阵。



上世纪四十年代末，杨大卫研究迭代方法伊始，对于在新生的电子计算机上使用迭代法求解大型问题的想法，有人表示怀疑。但自从他开创性的博士论文问世后，迭代法已被广泛应用于科学和工程中，并衍生出许多新的变种。

在科学计算的广阔领域，杨大卫的功绩永远不会被人遗忘，今年 12 月 21 日是他十五周年忌辰，谨以此短文感谢他，并简单介绍他留名数学史的博士论文成果—— SOR 方法。

写于 2023 年 11 月 23 日星期四，感恩节

美国哈蒂斯堡夏日山庄

原创 丁玖 返朴 2023-12-23 11:42 发表于北京