用R求矩陣的特徵值和特徵向量

最近在學習多元統計分析的主成分分析時,發現需要經常計算矩陣的特徵值和特徵向量,自己就找了下用R來做計算的函數。
我們可以用sigen()函數來計算特徵對。

#創建一個矩陣
a <- matrix(c(11,sqrt(3)/2,3/2,sqrt(3)/2,21/4,5*sqrt(3)/4,3/2,5*sqrt(3)/4,31/4),
            ncol=3,
            nrow=3
            )
#使用函數eigen()計算特徵對
ev <- eigen(a)
#訪問列表val項,即特徵值
ev$val
#訪問列表vec項,即特徵向量
ev$vec

Out:
在這裏插入圖片描述
a矩陣:
在這裏插入圖片描述
ev:
在這裏插入圖片描述

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章