← 返回Book 2目录
第29章

瑞利商可视化

Rayleigh Quotient Visualization

从二维到三维,揭示瑞利商与特征值的几何关系

29.1 什么是瑞利商?

瑞利商(Rayleigh Quotient)定义为:

R(x) = (xT Q x) / (xT x)

其中Q是对称矩阵,x是非零向量。分子xTQx是标量,分母xTx也是标量。

📐 二元函数形式

对于2×2矩阵Q = [[a,c],[c,b]],瑞利商是二元函数:
f(x₁,x₂) = (ax₁² + 2cx₁x₂ + bx₂²) / (x₁² + x₂²)

🎯 关键性质

• λmin ≤ R(x) ≤ λmax
• 在特征向量方向取得极值
• 与特征值分解密切相关

29.2 二元瑞利商

对于2×2对称矩阵Q,瑞利商可以看成是一个二元函数。 通过可视化等高线和3D曲面,直观看到不同矩阵Q如何影响瑞利商的形状。

正定矩阵:开口朝上的椭圆面,λmin > 0
负定矩阵:开口朝下的椭圆面,λmax < 0
半定矩阵:山谷或山脊,λmin=0 或 λmax=0
不定矩阵:双曲面(鞍形),有正有负

📈 二元瑞利商可视化

R(x) = x^T Q x / x^T x,观察不同矩阵Q如何影响瑞利商的形状

开口朝上正椭圆面(正定)

矩阵Q: [[1, 0], [0, 2]]

特征值: λ₁ = 2.00, λ₂ = 1.00

瑞利商范围: [1.00, 2.00]

🔵 正定/负定

所有特征值同号,椭圆面。 瑞利商始终为正或负。

⭕ 半正定/半负定

有零特征值,山谷或山脊。 某方向瑞利商为0。

🔄 不定

特征值有正有负,双曲面。 瑞利商可正可负。

29.3 在单位圆上看二元瑞利商

在单位圆 x₁² + x₂² = 1 上,使用极坐标表示:

x₁ = cos(θ), x₂ = sin(θ)

瑞利商成为角度θ的函数:R(θ) = xT Q x(单位圆约束下)

左图:单位圆散点,颜色映射瑞利商值
右图:R(θ)曲线,θ ∈ [0, 2π],值域恰好是[λmin, λmax]

💡 极坐标视角下,瑞利商随方向变化的规律一目了然!

29.4 在单位球上看三元瑞利商

三维情况下,在单位球面 x₁² + x₂² + x₃² = 1 上, 瑞利商 R(x) = xT Q x 直接等于二次型值。

球坐标参数化:
x₁ = sin(θ)cos(φ), x₂ = sin(θ)sin(φ), x₃ = cos(θ)
θ ∈ [0, π](极角),φ ∈ [0, 2π](方位角)
球面线框图:展示几何结构
颜色渲染:RdYlBu色谱映射瑞利商值
球面展开:(φ, θ)平面散点图、等高线图

💡 通过颜色映射球面,可以看到瑞利商在不同方向的分布,识别特征向量的大致方向!

🌐 三元瑞利商在单位球面上

在单位球面上:R(x) = x^T Q x(因为||x||=1),用颜色映射瑞利商值

正定矩阵

特征值: λ₁=1, λ₂=2, λ₃=3

瑞利商范围: [1, 3]

球面上每个点的颜色表示该方向的瑞利商值

xyz

✅ 正定/半正定

瑞利商≥0,球面颜色偏暖色。 半正定时某方向值为0(冷色)。

❌ 负定/半负定

瑞利商≤0,球面颜色偏冷色。 半负定时某方向值为0。

🔄 不定

瑞利商可正可负,球面呈现 冷暖色混合分布。

💡 几何意义:单位球面上每个点x对应一个方向。瑞利商R(x) = x^T Q x告诉我们沿该方向的"二次型能量"。 最大值λ_max和最小值λ_min分别对应Q的最大和最小特征值,出现在对应的特征向量方向。

29.5 平面上看三元瑞利商

将球面展开到θ-φ平面,类似地图的圆柱形投影。 这种展开让我们能够在平面上用散点图展示瑞利商分布。

展开映射:球面(θ,φ) → 平面矩形
• 横轴φ: [0, 2π] → [0°, 360°] (方位角)
• 纵轴θ: [0, π] → [0°, 180°] (极角)

网格散点图:θ-φ网格着色
热力图:密集采样+颜色映射
等高线图:等值线分布

29.5 平面上看三元瑞利商(球面展开)

将球面x₁²+x₂²+x₃²=1展开到(θ,φ)矩形平面

θ-φ平面热力图θφ0°90°180°0°90°180°270°360°瑞利商值2.000.50-1.00

🗺️ 球面展开原理

球坐标(θ,φ)自然映射到矩形:
• 横轴φ: 方位角 [0, 2π] → [0°, 360°]
• 纵轴θ: 极角 [0, π] → [0°, 180°]
• 类似地图的圆柱形投影

📊 可视化优势

• 平面展示便于识别周期性模式
• 等高线清晰显示极值位置
• 可以精确测量角度关系
• 补充3D球面的信息

Python核心代码:
# 球面网格
theta = np.linspace(0, np.pi, 101)
phi = np.linspace(0, 2*np.pi, 101)
Theta, Phi = np.meshgrid(theta, phi)

# 计算瑞利商
X = np.sin(Theta) * np.cos(Phi)
Y = np.sin(Theta) * np.sin(Phi)
Z = np.cos(Theta)
R = X*(Q[0,0]*X + Q[0,1]*Y + Q[0,2]*Z) + ...

# 平面展开图
plt.pcolormesh(Phi, Theta, R, cmap='RdYlBu_r')
plt.contour(Phi, Theta, R, levels=12, colors='white')

29.6 球面等高线展示三元瑞利商

在单位球面上直接绘制等高线,连接瑞利商值相等的点。 这些曲线揭示了二次型能量的"等势线"。

🎨 两种等高线方式

  • 球面白色等高线:在彩色球面上叠加白色等值线
  • 球面彩色等高线:等高线本身用颜色编码
  • 提取φ-θ平面等高线,映射回3D球面
  • 技术路线:contour → 坐标变换 → 3D投影

💡 等高线的意义

  • 连接所有具有相同瑞利商值的方向
  • 在特征向量方向,等高线最密集(梯度最大)
  • 等高线分布揭示矩阵Q的定性(正定/负定/不定)
  • 帮助识别极值点和鞍点

可视化策略:球面等高线与平面展开图互补—— 球面版保留几何直觉,平面版便于精确测量。RayleighQuotient3DExplorer组件已包含球面颜色渲染,SphereUnfoldingDemo组件提供平面等高线视图。

🎯 核心洞察与应用

几何直觉

  • 瑞利商将抽象的特征值问题转化为几何优化问题
  • 在单位球/圆上,R(x)的最大值=λmax,最小值=λmin
  • 正定矩阵→山丘,负定→盆地,不定→鞍形
  • 半定矩阵在某些方向为零(山谷/山脊底部)

实际应用

  • PCA主成分分析:最大化瑞利商找数据方差最大方向
  • 谱聚类:最小化归一化拉普拉斯矩阵的瑞利商
  • 量子力学:基态能量通过最小化哈密顿量的瑞利商估计
  • 变分原理:Rayleigh-Ritz定理的几何理解

多视角可视化:从单位圆到单位球,从等高线到色彩映射,从球面到平面展开, 不同维度的可视化互补揭示瑞利商的几何结构,让抽象的代数概念变得触手可及。