作者embed (暱稱是我)
看板Statistics
標題[問題] R中的非調和陣列
時間Thu Oct 30 17:25:21 2008
自己寫了一個雙重迴圈
程式碼如下(A是一個5*5方陣)
orthonormal <- function(A){
W <- matrix(nrow=5, ncol=5)
V <- diag(x=1, nrow=5)
for (i in 1:dim(V)[1]){
for (j in 1:i-1){
V[,i] <- A[,i] - (prod(V[,j],A[,i])/norm2(V[,j]))*V[,j]
} }
W <- V
show(W)
}
但是執行的結果卻顯示"錯誤在prod(V[, j], A[, i])/norm2(V[, j]) : 非調和陣列"
卻不知道怎麼下手解決...
依照平常計算gram-shmidt的公式 這樣子的表達式是可以理解的
究竟該怎麼解決非調和陣列的問題?
謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.4.234
1F：推 bugle:當i=1時.... 10/30 22:00
2F：→ clickhere:Gram-Schmidt 用的應該是 V[,j] %*% A[,i] 而非 prod 10/31 01:07