美文网首页
由本征值计算本征向量

由本征值计算本征向量

作者: 小李飞刀s3 | 来源:发表于2020-04-06 11:28 被阅读0次

最近需要求解任意本征值的本征矢量,目前有两种方法计算本征值,一种是传统的方法,另外一种是陶哲轩[1]去年提出的一种方法。当然后来发现这种方法不是他首次提出,但是确实是一种新的思路。本次笔记就是为了展示如何用他的结论。例子也是他的论文找的。

陶哲轩

计算本征值

计算本征值很容易,例如对于矩阵M,

MX=\lambda X\implies (M-\lambda I)X=0\implies \det(M-\lambda I)=0

方法简述

传统方法

(M-\lambda I)X=0
将计算的本征值带回原矩阵求一个代数方程即可。

陶哲轩新方法

|v_{i,j}|^2 \prod_{k=1,k\neq i}^{n}\left(λ_i(A)-λ_k(A)\right)=\prod_{k=1}^{n-1}\left(λ_i(A)-λ_k(M_j)\right)

其中

  • v_{i,j} 是特征值\lambda_{i} 对应特征向量的第j个元素
  • \lambda_{i}(A)是矩A的第i个特征向量(本征向量、本征矢)
  • M_{j}为矩阵A的第j个余子式,\lambda_{k}(M_{j})是该余子式的第k个特征值。

这种方法的好处是直接带公式就可以,缺点是需要多计算几个本征值。

示例

我们令

M=\left[\begin{matrix} 1 &1 &-1\\ 1 &3 &1\\ -1 &1 &3 \end{matrix}\right]
其本征值可以通过解方程得出,有三个,

\lambda_{1}=0,\lambda_{2}=3,\lambda_{3}=4

传统方法

传统方法,求解本征矢,对于\lambda=0,我们有

\left[\begin{matrix} 1-0 &1 &-1\\ 1 &3-0 &1\\ -1 &1 &3-0 \end{matrix}\right]\cdot\left[\begin{matrix} x\\ y\\ z \end{matrix}\right]=0

v_{1}=\frac{1}{\sqrt{6}}\left[\begin{matrix} 2\\ -1\\ 1 \end{matrix}\right]
同样的做法,我们可以知道

v_{2}=\frac{1}{\sqrt{3}}\left[\begin{matrix} 1\\ 1\\ -1 \end{matrix}\right]
v_{3}=\frac{1}{\sqrt{2}}\left[\begin{matrix} 0\\ 1\\ 1 \end{matrix}\right]

陶哲轩方法

按照陶哲轩的方法,首先将不同的余子式表示出来

M_{1}=\left[\begin{matrix} 3 &1\\ 1 &3 \end{matrix}\right],\quad \lambda_{11}=4,\quad \lambda_{12}=2
M_{2}=\left[\begin{matrix} 1 &-1\\ -1 &3 \end{matrix}\right],\quad \lambda_{21}=2+\sqrt{2},\quad \lambda_{22}=2-\sqrt{2}
M_{1}=\left[\begin{matrix} 1 &1 \\ 1 &3 \\ \end{matrix}\right],\quad \lambda_{21}=2+\sqrt{2},\quad \lambda_{22}=2-\sqrt{2}

这样就可以依次求出系数了,

|v_{1,1}|^2=\frac{(0-2)(0-4)}{(0-3)(0-4)}=\frac{2}{3}
|v_{1,2}|^2=\frac{(0-2-\sqrt{2})(0-2+\sqrt{2})}{(0-3)(0-4)}=\frac{1}{6}
|v_{1,3}|^2=\frac{(0-2-\sqrt{2})(0-2+\sqrt{2})}{(0-3)(0-4)}=\frac{1}{6}

|v_{2,1}|^2=\frac{(3-2)(3-4)}{(3-0)(3-4)}=\frac{1}{3}
|v_{2,2}|^2=\frac{(3-2-\sqrt{2})(3-2+\sqrt{2})}{(3-0)(3-4)}=\frac{1}{3}
|v_{2,3}|^2=\frac{(3-2-\sqrt{2})(3-2+\sqrt{2})}{(3-0)(3-4)}=\frac{1}{3}

|v_{3,1}|^2=\frac{(4-2)(4-4)}{(4-0)(4-3)}=0
|v_{3,2}|^2=\frac{(4-2-\sqrt{2})(4-2+\sqrt{2})}{(4-0)(4-3)}\frac{1}{2}
|v_{3,3}|^2=\frac{(4-2-\sqrt{2})(4-2+\sqrt{2})}{(4-0)(4-3)}=\frac{1}{2}

可见该方法的不足是不能分辨正负。


  1. https://arxiv.org/abs/1908.03795

相关文章

  • 由本征值计算本征向量

    最近需要求解任意本征值的本征矢量,目前有两种方法计算本征值,一种是传统的方法,另外一种是陶哲轩[1]去年提出的一种...

  • 人工智能学习笔记-Day15

    本征值计算 本征向量可以作为非常好的基,同时本征值表示线性变换在此方向上的放缩。 若一个矩阵将一个向量变换成0向量...

  • 第五章 本征值与本征向量

    title: 第五章 本征值与本征向量category: 笔记date: 2019/10/11mathjax: t...

  • 2.21 本征向量和本征值

    https://www.youtube.com/watch?v=5sHCGsEYtLE&list=PL65jGfV...

  • 厄米变换

    特征:1.本征值为实数。2.不同本征值的本征矢彼此正交。3.本征矢构成完备集。

  • 我理解的本征值,本征态和表象

    我理解的量子力学中,波函数就是一个系统的各种可能的状态的叠加,比如描述一个系统,只考虑电子的自旋这个参数的话,只有...

  • 本征矢

    在一个复矢量线性空间中,对每一线性变换都存在一个“特殊”的向量,这个向量接受变换后等于原向量乘以一个标量。这种向量...

  • Arduino教学拓展——二极管内部结构及其原理

    在讲二极管之前先引入本征半导体和杂质半导体的概念 一、本征半导体和杂质半导体 (1)本征半导体 什么是半导体? 导...

  • Fermi黄金规则

    Fermi黄金规则(Fermi's golden rules)是一个描述量子体系从一个能量本征态到另一个能量本征态...

  • 关于inversed eigenvalue problems

    概述一下问题:使用等谱流方法(isospectral flow method)求解逆谱问题。即给定一个矩阵的本征值...

网友评论

      本文标题:由本征值计算本征向量

      本文链接:https://www.haomeiwen.com/subject/uivvphtx.html