当前位置: 首页 > news >正文

谱分解与Cholesky分解

要真正理解一般正规算子的谱分解是困难的,你几乎是绕不开“算子代数”的,我现在说一下思路:

1:有界的正规算子

设 TT 是希尔伯特空间上的有界线性算子,它可以生成一个交换的 C∗C^* 算子代数 ATA_T ,而且它的极大理想可以和 σ(T)\sigma(T) 构成的对应,因此我们可以有一个从 ATA_T 到 C(σ(T))C(\sigma(T)) 的Gelfand map Γ\Gamma ,而且这个映射是 ∗* 代数同构的,也就是说,反过来,我们随便给一个连续函数 u∈C(σ(T))u\in C(\sigma(T)) ,那么 Γ−1(u)\Gamma^{-1} (u) 是一个有界线性算子,特别的,我们可以有 Γ−1(1)=I\Gamma^{-1}(1)=I , Γ−1(z)=T\Gamma^{-1}(z)=T , 假设 σ(T)={λ1,⋯,λn}\sigma(T)=\{\lambda_1,\cdots,\lambda_n\} ,也就是说谱刚刚好是有限个点谱, 如果是一个矩阵,那么它的谱肯定是点谱。那么如果这个时候我们考虑函数 χλi\chi_{\lambda_i} ,也就是每一个点上的示性函数,在这个离散的情况下,这个函数是连续函数,所以我们可以构造 Ei=Γ−1(χλi)E_i=\Gamma^{-1}(\chi_{\lambda_i}) ,不难证明 EiE_i 是正交投影算子,而且

Tx:=∑λiEixTx:=\sum \lambda_i E_i x .

发现核心问题了吗?对于一般的 σ(T)\sigma(T),示性函数不是连续函数,也就是说Gelfand map没用了,可是投影算子是对应示性函数的,事实上 设 PP 是投影算子,那么 P2=PP^2=P ,如果它对应一个 σ(T)\sigma(T) 的函数 uu ,那么 u2(z)=u(z)u^2(z)=u(z) ,所以在每个点这个函数只能是0或者1.

为了处理一般的算子,我们需要引入 W∗W^* -算子代数 WTW_T ,本质上是 ATA_T在一个弱拓扑 下的闭包,

最大的优点是可以让 Γ\Gamma 可以覆盖示性函数从而得到

Tx:=∫λdEλTx:=\int \lambda d E_\lambda (这里也从离散和变为一般的积分)

“一般的对称算子“

首先根据上面的结论,对于任何unitary算子 UU ,因为它的谱肯定在一个圆环上,我们可以找到一个谱分解

U=∫02πeiθdFθU=\int_0^{2\pi} e^{i\theta} d F_\theta ,

然后对于任何无界的对称算子 AA ,Cayley 变换 B=A−iIA+iIB=\frac{A-i I}{A+iI} 是一个unitary算子,也就是说

B=∫02πeiθdFθB=\int_0^{2\pi} e^{i\theta} d F_\theta ,如我们引入定义 λ=−cot⁡(θ/2),E(λ):=F(θ)\lambda=-\cot(\theta/2), E(\lambda):=F(\theta) ,那么变换可以把一个圆环拉成实数轴,这个时候

U=∫Rλ−iλ+idE(λ)U=\int_{\mathbb{R}}\frac{\lambda -i}{\lambda+i} dE(\lambda) ,

如果我们定义算子

C=∫RλdE(λ)C=\int_{\mathbb{R}} \lambda dE(\lambda) ,

可以证明 C=AC=A ,这里主要需要一些Cayley变换和对称算子的一些性质。

很惭愧我泛函学得不怎么样,没法像dhchen那样做深入浅出的分析。对我来说,谱分解就是有限维线性代数中“复正规矩阵都可以酉对角化”这个结论的无穷维推广;在无穷维Hilbert空间中,对角矩阵的类似物就是乘法算子(把每一个函数乘以一个固定的函数f给出的算子),那么谱分解定理可以表述为:对于Hilbert space H上的normal operator A(我实在不喜欢 正常算子 这个翻译。。),存在某个紧Hausdorff拓扑空间X (我不记得这个X是不是就是A的spectrum),以及unitary equivalence ϕ:H→L2(X)\phi:H\to L^2(X) , 乘法算子 F:L2(X)→L2(X), g↦gfF:L^2(X)\to L^2(X),\ g\mapsto gf for some f∈C(X)f \in C(X) ,使得 A=ϕ−1∘F∘ϕA=\phi^{-1}\circ F\circ \phi .

这是我本科的泛函老师给出的一种比较简单的解释,实际上就是说无穷维的normal operator也是酉等价于“无穷维对角矩阵”(也就是乘法算子)的,只是你要把空间换一下。当然也可以基于谱测度或者算子代数去理解,像dhchen解释的那样,我不怎么懂泛函就不多说了。


http://www.mrgr.cn/news/39232.html

相关文章:

  • 《后端程序猿 · Spring事务失效场景》
  • g2o 中定义好的常用边的类型
  • 修改qtcreator创建项目的模板
  • DialMAT:跨模态特征提取与对抗训练的结合
  • Arduino Uno板实现流水灯的设计方案
  • 六、动画系统
  • 9.29 LeetCode 3304、3300、3301
  • B - 1D Keyboard题解
  • k8s架构,从clusterIP到光电半导体,再从clusterIP到企业管理
  • SD2.0 Specification之CPRM(Content Protection for Recordable Media)
  • 「OC」多线程的学习——NSThread
  • Redis: RDB与AOF的选择和容灾备份以及Redis数据持久化的优化方案
  • ADC的原理
  • 二进制字符串的和
  • spring boot项目对接人大金仓
  • CSS 的color-adjust 属性,是给打印机节省墨水的属性
  • 计数排序(counting sort)
  • 文件传输工具 | 闪电藤 v2.5.5 绿色版
  • MFC工控项目实例之十八手动测试界面输入信号实时检测
  • 算法:852.山脉数组的峰顶索引