斯图姆–刘维尔理论与本征函数展开

瑞利商(Rayleigh quotient)

/ RAY-lee /

不精确求解鼓的方程,怎样估计鼓的最低固有频率?瑞利商就是答案:由你随手猜的任一试探形状构造的一个数,它高估最低本征值,而你的猜测越好它就越接近。它把本征值问题化为一个极小化问题,后者在近似处理上往往容易得多。

对斯图姆-刘维尔算子,满足边界条件的试探函数 u 的瑞利商是 R[u] =(integral of [p (u')^2 - q u^2] dx)除以(integral of u^2 w dx)。分子在物理上是某种储能(弯曲或拉伸能加一个势能项);分母是带权范数。两个事实使它有力。其一,若代入精确本征函数 y_n,便原样得回其本征值 lambda_n。其二,该商在每个本征函数处驻定,并在所有容许 u 上被最低本征值极小化:lambda_1 等于 R[u] 的最小值。于是对任何并非真正基态的试探 u,R[u] 大于 lambda_1——一个你能算出的上界。

这是瑞利-里茨方法的引擎,也是物理与工程的常备工具。猜一个带几个自由参数的简单试探函数,把 R[u] 算成这些参数的函数,再极小化;结果是基本本征值(最低振动频率、量子力学中的基态能量、柱的屈曲载荷)的一个严格上界。即便粗糙的猜测也常落在百分之几以内,因为该商在真实模式附近驻定——形状的小误差只造成本征值的二阶误差。

对 [0, 1] 上的 y'' + lambda y = 0、y(0) = y(1) = 0,真实 lambda_1 为 pi^2 = 9.8696。试用抛物线 u = x(1 - x):R[u] =(integral of (u')^2)/(integral of u^2)=(1/3)/(1/30)= 10。仅比精确值高 1.3%,而这来自一行的猜测。

一个粗糙的试探形状仅把基本本征值高估了一个百分点左右——这正是变分方法的标志。

瑞利商的最小值只给出最低本征值,且总是从上方逼近;要够到 lambda_2,须把试探函数限制为与基态正交,且估计仍是高估、绝不会低估。

又称
Rayleigh ratioenergy quotient瑞利比