詹姆斯·威尔金森
概述
詹姆斯·哈迪·威尔金森(James Hardy Wilkinson,1919–1986)是英国数学家,现代数值线性代数的奠基人。1970年 ACM 图灵奖得主,1969年当选英国皇家学会院士(FRS)。剑桥大学三一学院毕业(师从 Hardy、Littlewood),二战期间从事弹道计算,1946年加入英国国家物理实验室(NPL),与阿兰·图灵共同参与 ACE 计算机设计,并亲手设计了 Pilot ACE 的乘法运算单元。在 NPL 与早期计算机十余年的实际计算经验,使他从计算实践中提炼出后向误差分析(backward error analysis)这一革命性方法论,系统建立了数值线性代数的理论基础。其1965年专著《代数特征值问题》被公认为数值线性代数的奠基之作。
关键内容
生平
| 时间 | 事件 |
|---|---|
| 1919年9月27日 | 生于英国 |
| 剑桥大学三一学院 | 师从 Hardy、Littlewood、Besicovitch,一等荣誉学位 |
| 二战期间 | 弹道计算,用手摇计算器解12阶线性方程组 |
| 1946年 | 加入英国国家物理实验室(NPL) |
| 1946年 | 与阿兰·图灵</td> </tr> <tr> <td>1963年</td> <td>出版《代数过程中的舍入误差》(<em>Rounding Errors in Algebraic Processes</em>)</td> </tr> <tr> <td>1965年</td> <td>出版《代数特征值问题》(<em>The Algebraic Eigenvalue Problem</em>)</td> </tr> <tr> <td>1969年</td> <td>当选英国皇家学会院士(FRS)</td> </tr> <tr> <td>1970年</td> <td>荣获 ACM [[阿兰·图灵</td> </tr> <tr> <td>1971年</td> <td>与 Reinsch 合编《自动[[计算手册·线性代数卷》(EISPACK 前身) |
| 1986年10月5日 | 于英国逝世 |
后向误差分析(最重要贡献)
威尔金森的核心思想转变:将"计算结果距真实值多远"(正向误差)改为问"计算结果精确地解了哪个问题"(后向误差)。
关键洞见:在 Pilot ACE 上计算多项式零点时,他发现计算得到的零点虽然不是原多项式的精确零点,但是某个系数微小扰动后多项式的精确零点——后向误差分析由此萌生。
方法论意义: - 将算法稳定性(backward error $|\delta A|/|A|$)与问题条件性(条件数 $\kappa(A)$)清晰分离 - "后向稳定的算法"有了客观判定标准 - 彻底推翻了 Hotelling(1943)关于 Gauss 消去法误差指数增长的悲观预言
详见 后向误差分析。
QR 算法的 Wilkinson 位移
对QR算法的重要改进:对称三对角矩阵的末部 $2\times 2$ 子矩阵: $$\mu_k = a_{nn} - \frac{\operatorname{sign}(\delta)\cdot b_{n-1}^2}{|\delta| + \sqrt{\delta^2 + b_{n-1}^2}}, \quad \delta = \frac{a_{n-1,n-1} - a_{nn}}{2}$$
选取更接近末元素 $a_{nn}$ 的特征值作为位移,保证全局收敛性和三次方渐近收敛速度,在实际计算中每个特征值通常只需1–2次迭代即可隔离。
《代数特征值问题》(1965)
牛津大学出版社,662页,数值线性代数领域引用率最高的专著之一,至今未绝版。LAPACK 用户指南扉页写道:"本工作献给 Jim Wilkinson,他的思想与精神在项目的每一个环节都给予我们启迪和影响。"
主要内容: - 后向误差分析框架(全书灵魂) - 特征值条件数理论($\kappa(\lambda_i) = 1/|y_i^H x_i|$) - QR算法收敛性分析与 Wilkinson 位移 - Gauss 消去法的后向稳定性(彻底解决 Hotelling 问题) - 正交变换(Givens旋转、Householder反射)的数值优越性
软件遗产
EISPACK(1970s)→ LINPACK(1970s)→ LAPACK(1992至今)→ MATLAB、NumPy/SciPy 的核心计算引擎,均以威尔金森的理论框架和算法为基础。
SIAM 设立了两项以其名字命名的奖项(1982年"数值分析与科学计算奖",1991年"数值软件奖"),是数值计算领域最高荣誉之一。
来源
- raw/books/矩阵分析/14_wilkinson_algebraic_eigenvalue_problem_1965.md