日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

EIGENVECTORS FROM EIGENVALUES论文结论举例验证

發布時間:2023/12/20 编程问答 40 豆豆
生活随笔 收集整理的這篇文章主要介紹了 EIGENVECTORS FROM EIGENVALUES论文结论举例验证 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

陶哲軒論文原文[1]

論文總共三種證明方法來證明同一個理論:

作者博客[2]的其中一種表達形式如下:

其中:

下面舉個例子來說明下陶哲軒論文里面的上述公式到底是什么意思.
設矩陣A為:
{?141430102}\left\{ \begin{matrix} -1 & 4 & 1 \\ 4 & 3 & 0 \\ 1 & 0 & 2 \\ \end{matrix} \right\} ?????141?430?102?????

#################python代碼求解矩陣A的所有特征值和特征向量:############

import numpy as np from numpy import * # X=[ [1,2,1,1], # [3,3,1,2], # [3,5,4,3], # [5,4,5,4], # [5,6,1,5], # [6,5,2,6], # [8,7,1,2], # [9,8,3,7]] # X=np.array(X).T#這里注意,[1,2,1,1]在numpy的眼中是一列np.linalg.eig X=[[-1,4,1],[4,3,0],[1,0,2]] print("X=",X) X=matrix(X) print("------------------下面計算原始矩陣的特征值和特征向量-----------------------") eigenvalue,featurevector=np.linalg.eig(X) print("原始矩陣的特征值") print("eigenvalue=",eigenvalue) print("---------------------------") print("featurevector=\n",featurevector)

運行代碼得到特征值
eigenvalue=[-3.60201123 5.55175228 2.05025894]


得到特征向量:
[[-0.8454712 0.53175808 0.04911014]
[ 0.51225069 0.83355752 -0.20683592]
[ 0.1509228 0.14971711 0.97714231]]
例如-3.60201123對應的特征向量就是:
-0.8454712, 0.51225069, 0.1509228

####################下面分析該公式##############################
∣vi,1∣2|v_{i,1}|^2vi,1?2表示:
第i個單位化特征向量,1表示特征向量中的第一個元素的值.
已知上面的-0.8454712, 0.51225069, 0.1509228已經是單位化處理過的.

所以∣vi,1∣2|v_{i,1}|^2vi,1?2=(?0.8454712)2=0.7148215500(-0.8454712)^2=0.7148215500(?0.8454712)2=0.7148215500

##################下面計算式(1)右側的分子##################
∏k=1n?1(λi(A)?λk(M))\prod_{k=1}^{n-1}(\lambda_i(A)-\lambda_k(M))k=1n?1?(λi?(A)?λk?(M))

已知M的特征值是3和2.代入分子:

∏k=1n?1(λi(A)?λk(M))\prod_{k=1}^{n-1}(\lambda_i(A)-\lambda_k(M))k=1n?1?(λi?(A)?λk?(M))
=(-3.60201123-3)(-3.60201123-2)
=36.9845410510

注意這里的λi(A)\lambda_i(A)λi?(A)在連乘運算中是定死的,就是viv_ivi?對應的特征向量λi\lambda_iλi?

#######################下面計算式(1)右側分母######################
∏k=1,k≠in?1(λi(A)?λk(A))\prod_{k=1,k≠i}^{n-1}(\lambda_i(A)-\lambda_k(A))k=1,k?=in?1?(λi?(A)?λk?(A))
分母的意思就是"矩陣A的viv_ivi?向量對應的特征值"-“矩陣A的其他特征值”
所以:
∏k=1,k≠in?1(λi(A)?λk(A))\prod_{k=1,k≠i}^{n-1}(\lambda_i(A)-\lambda_k(A))k=1,k?=in?1?(λi?(A)?λk?(A))
=(-3.60201123-5.55175228)(-3.60201123-2.05025894)
=51.7395444308
#######################檢驗等式是否成立##################
根據上面的計算:
分子=36.9845410510
分母=51.7395444308
所以式(1)右側=36.9845410510/51.7395444308=.7148215442
∣vi,1∣2=(?0.8454712)2=.7148215442|vi,1|^2=(-0.8454712)^2=.7148215442vi,12=(?0.8454712)2=.7148215442

所以式(1)成立.

由此可知,根據式(1),
可以直接根據原矩陣的特征值子矩陣的特征值,直接得到原矩陣的向量中的任一元素的絕對值

該論文的特點:
目前難以得到特征向量里面的元素的符號,只能得到元素的絕對值或者平方.
需要事先知道原矩陣和子矩陣各自的所有特征值.

知乎[3]上更新:
徐樹方《矩陣計算的理論與方法》(1995年)Page 323th的引理3.1已經有類似結論
所以,總的來看,其實是一場誤會,這個公式在20年前就已經被研究過了.

陶哲軒本人在博客[2]的評論回復中,被人挖出來是很久以前就出版的內容以后,之后基本就不回復了了,所以看熱鬧的,散了吧.

Reference:
[1]https://arxiv.org/pdf/1908.03795.pdf
[2]https://terrytao.wordpress.com/2019/08/13/eigenvectors-from-eigenvalues/#comment-528850
[3]https://www.zhihu.com/question/355978404/answer/895583246

總結

以上是生活随笔為你收集整理的EIGENVECTORS FROM EIGENVALUES论文结论举例验证的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。