在進(jìn)行數(shù)據(jù)分析的時(shí)候有時(shí)候需要求建立的矩陣的特征值,尤其是最大特征值以及其對(duì)應(yīng)的特征向量,隨著矩陣階數(shù)的擴(kuò)大,運(yùn)算量比較大,但是如果使用軟件Matlab來計(jì)算可以節(jié)省很多時(shí)間
文章源自四五設(shè)計(jì)網(wǎng)-http://www.133122.cn/39924.html
啟動(dòng)Matlab ,在命令窗口輸入要處理的矩陣A文章源自四五設(shè)計(jì)網(wǎng)-http://www.133122.cn/39924.html
A=[1,4,2,4;1/4,1,1/2,1;1/2,2,1,1/2;1/4,1,2,1]文章源自四五設(shè)計(jì)網(wǎng)-http://www.133122.cn/39924.html
輸入完成后回車軟件會(huì)按行列的形式顯示矩陣文章源自四五設(shè)計(jì)網(wǎng)-http://www.133122.cn/39924.html
順便可以檢查一下矩陣是否輸入錯(cuò)誤文章源自四五設(shè)計(jì)網(wǎng)-http://www.133122.cn/39924.html
文章源自四五設(shè)計(jì)網(wǎng)-http://www.133122.cn/39924.html
接著輸入[x,y]=eig(A) 回車文章源自四五設(shè)計(jì)網(wǎng)-http://www.133122.cn/39924.html
就可以看到矩陣的所有特征值和特征向量了文章源自四五設(shè)計(jì)網(wǎng)-http://www.133122.cn/39924.html
特征值是對(duì)角矩陣y 矩陣x的每一列對(duì)應(yīng)一個(gè)y中相應(yīng)列的特征值文章源自四五設(shè)計(jì)網(wǎng)-http://www.133122.cn/39924.html
此處注意括號(hào)必須是在英文輸入法下輸入 如果顯示紅色表示有問題需要重新輸入文章源自四五設(shè)計(jì)網(wǎng)-http://www.133122.cn/39924.html
雖然已經(jīng)計(jì)算出了矩陣的特征值和特征向量 ,但是如果只是計(jì)算這一個(gè)的話 ,完全可以觀察得到最大的特征值
接著輸入
eigenvalue=diag(y);
lamda=eigenvalue(1)
lamda表示最大的特征值 這里lamda=4.2498
求出矩陣最大特征值之后雖然可以一眼看到特征向量, 但是不方便以后大量的處理矩陣
此時(shí)再輸入
y_lamda =?x(:, 1) 即最大特征值對(duì)應(yīng)的特征向量
上面的過程只是適合進(jìn)行一次兩次這樣的計(jì)算, 如果遇到需要較多的矩陣的時(shí)候或者這樣的過程只是某些計(jì)算中的幾步需求的時(shí)候就需要把前面的計(jì)算合在一起進(jìn)行計(jì)算,同樣也可以把這段代碼放在需要的算法程序中
除了直接輸入變量結(jié)果來查看的方法外 ,還可以直接在workspace中查看變量運(yùn)算結(jié)果 ,如圖中的lamda和y_lamda的結(jié)果值


評(píng)論