Skip to content

预解算子与预解集

Resolvent operator · Resolvent set

由复参数移位后的算子逆定义预解集、谱与对扰动敏感性的尺度。

形式陈述

X 是复 Banach 空间,AB(X)有界线性算子。若对 zC,算子

zIA:XX

是双射,则有界逆定理保证其逆也是有界线性算子。所有这样的 z 构成 A 的预解集

ρ(A)={zC:zIA 可逆且 (zIA)1B(X)}.

预解算子是定义在 ρ(A) 上的算子值函数

R(z,A)=(zIA)1.

补集

σ(A)=Cρ(A)

称为 A 的谱。对有界算子,σ(A) 是非空紧集且包含于闭圆盘 {|z|A}。预解集是开集,zR(z,A) 在其每个连通分支上是算子范数意义的解析函数。

两个预解参数满足预解恒等式

R(z,A)R(w,A)=(wz)R(z,A)R(w,A).

由此也可见同一算子的预解彼此交换。若 (zw)R(w,A)<1,Neumann 级数给出

R(z,A)=R(w,A)k=0((zw)R(w,A))k,

从而在 w 附近仍可逆,并量化预解集的开放性。

直觉

zIA 询问的是:给定输出 y,方程 (zIA)x=y 是否对每个 y 有唯一、连续依赖于 y 的解。预解算子把右端项送回解,R(z,A) 就是该求解问题的最坏放大倍数。靠近谱时,某些方向几乎被 zIA 压扁,逆映射必须把它们猛烈放大。

谱不仅记录严格不可逆的位置;预解范数还描绘“离不可逆有多近”。对正规矩阵,这个大小完全由 z 到谱的距离决定。对非正规矩阵,特征向量接近平行会使远离特征值的区域也出现巨大放大,这正是伪谱比特征值图提供更多信息的原因。

例子与计算

取对角矩阵 A=diag(1,3)。当 z1,3 时,

R(z,A)=((z1)100(z3)1).

在 Euclidean 算子范数下,

R(z,A)2=max{1|z1|,1|z3|}=1dist(z,{1,3}).

这里谱恰为两个特征值,预解范数的等高线是围绕它们的圆弧。

非正规例子取

A=(0M00),M>0.

它只有特征值 0,但对 z0

R(z,A)=(z1Mz20z1).

即使 |z|=1,当 M 很大时预解范数也至少为 M 量级,远大于 1/dist(z,σ(A))=1。只看特征值会漏掉这种强烈的受迫响应与扰动敏感性。

边界与失败情形

有限维中,不可逆等价于行列式为零,所以谱等于特征值集合。无穷维中谱还可能包含没有特征向量的点。例如右移算子在某些序列空间上的谱包含整个单位圆盘,而点谱只占其中一部分;“找不到特征向量”不能推出参数在预解集中。

对无界算子,定义还必须要求 zIA 从算子定义域 D(A) 双射到 X,且逆作为 XX 的算子有界。不能忽略稠密定义域、闭性等条件,把无界微分算子当成普通矩阵处理。本页的正式陈述限于全空间上有界的 A

预解范数大并不表示 z 已经属于谱;它表示求逆病态。给定 ε>0,常用

σε(A)=σ(A){zρ(A):R(z,A)>ε1}

定义 ε-伪谱。阈值与所选算子范数有关,比较结果时必须声明范数。

正负号约定也会造成表面差异。有些文献定义 (AzI)1;它与本页的 R(z,A) 相差负号,预解集和范数不变,但预解恒等式的符号应随约定一致。

推论与应用

|z|>A 时,

R(z,A)=1z(IAz)1=1zk=0Akzk,

所以大圆盘外必属于预解集,并有 R(z,A)1/(|z|A)。这个级数把算子幂、谱半径与复分析联系起来。

预解的围道积分可定义谱投影和解析函数演算。在线性动力系统中,Laplace 变换后的解通过 (zIA)1 表示;在数值分析中,预解估计控制非正规演化的瞬态增长以及离散方案对扰动的响应。

参考资料
  • Walter Rudin, Functional Analysis, 2nd ed., McGraw-Hill, 1991,Ch. 10。
  • John B. Conway, A Course in Functional Analysis, 2nd ed., Springer, 1990,Ch. VII。
  • Lloyd N. Trefethen and Mark Embree, Spectra and Pseudospectra, Princeton University Press, 2005,Ch. 2。