如何從矢量創建矩陣

-3

我有一個矩陣mydata。我把它的本徵如何從矢量創建矩陣

ev <- eigen(mydata) 

       [,1]     [,2]     [,3]     [,4]     [,5]   [,6] 
[1,] -0.44852453+0i 0.11670674-0.43102606i 0.11670674+0.43102606i 0.02678808-0.00239041i 0.02678808+0.00239041i 0.013711630+0i 
[2,] -0.40701659+0i 0.17265390+0.01329737i 0.17265390-0.01329737i 0.32050237+0.31677316i 0.32050237-0.31677316i -0.946287625+0i 
[3,] -0.65381971+0i 0.35662142+0.31807798i 0.35662142-0.31807798i -0.83423480+0.00000000i -0.83423480+0.00000000i -0.265152784+0i 
[4,] -0.37545331+0i -0.72987261+0.00000000i -0.72987261+0.00000000i -0.00637960-0.10181135i -0.00637960+0.10181135i 0.004667737+0i 
[5,] -0.05763717+0i 0.02473111+0.01282302i 0.02473111-0.01282302i -0.00818839+0.03180396i -0.00818839-0.03180396i 0.122647730+0i 
[6,] -0.24778175+0i

萬一結果沒有排序，我得到對應的最大特徵值從而最大特徵向量：

k <- which(abs(ev$values)==max(abs(ev$values))) 
G<-ev$vectors[,k] 

[1] -0.44852453+0i -0.40701659+0i -0.65381971+0i -0.37545331+0i -0.05763717+0i -0.24778175+0i

調用這個最大特徵向量

G = [克，克 ...克 ñ ] Ť

因此，在這種情況下，G = [-0.44852453 + 0I -0.40701659 + 0I -0.65381971 + 0I -0.37545331 + 0I -0.05763717 + 0I -0.24778175 + 0I，其中n = 6

我的問題是，我想從這個特徵向量構建一個比率矩陣，以便B = {bij = gi/gj}。 R中有這樣的函數嗎？如果不是什麼是最簡潔的方式來獲得沒有循環的矩陣？ {抱歉，我不知道如何在這個界面中顯示數學符號。應該看，矩陣B =條目b（子ij）= g（子i）/ g（子j），g是G的成員}

來源

2014-12-26 Ivan

我認爲這與特徵值無關，想要從矢量創建矩陣B，其中 = G[i]/G[j]用於矢量G。

這可以通過outer函數實現，該函數生成外部產品，但可以指定函數。

G <- 1:4 
outer(G, G, FUN='/') 
##  [,1] [,2]  [,3] [,4] 
## [1,] 1 0.5 0.3333333 0.25 
## [2,] 2 1.0 0.6666667 0.50 
## [3,] 3 1.5 1.0000000 0.75 
## [4,] 4 2.0 1.3333333 1.00

來源

2014-12-26 19:51:27

對不起，這是不正確的。例如，矩陣是6x6。對角線條目清楚地全部爲1. [1,2] = -0.44852453/-0.40701659 = 1.101980954。 [2,1] = -0.40701659/-0.44852453 = 0.9074567 [5,4] = -0.05763717/-0.37545331 = 0.153513549等等。 – Ivan

這個結果並沒有使用你的輸入向量，所以當然這些值是不正確的！試着將'outer'應用到你的實際向量上。 –

對不起，duuh！謝謝。 – Ivan

如何從矢量創建矩陣

回答

相關問題