作者ru486good (金仔)
看板NCCU09_STGRA
标题[讨论] 统模第二题
时间Mon May 24 21:15:43 2010
power method那题
这题卡关好久啊
蔡丞庭的方法总让我觉得怪怪的...
幸好有Google大神
让我可以顺利找到剩下的 那三个 害羞的 eigen vectors
大家可以参考下面这个网站
http://www.miislita.com/information-retrieval-tutorial/matrix-tutorial-3-
eigenvalues-eigenvectors.html
#power-method
不会缩网址 将就一下
看The Power Method (Vector Iteration) 、The Deflation Method这两个章节
应该就OK了
###########code 参考看看############
#Power Method
p.m = function(A,v=c(rep(1,nrow(A))),tol=10^-4){
w = c(rep(0,nrow(A)))
iter = 0
while( mean(abs(w-v)) >= tol){
w = v
v = A%*%v
b = sum(v^2)
v = v/sqrt(b)
iter = iter+1
}
list(value=sqrt(b),vector=v,iter=iter)
}
a1 = p.m(A) ; a1
## Deflation Method
Z = a1$value*a1$vector%*%t(a1$vector)
A2 = A - Z
a2 = p.m(A2) ; a2
#
Z2 = a2$value*a2$vector%*%t(a2$vector)
A3 = A2 - Z2
a3 = p.m(A3) ; a3
#
Z3 = a3$value*a3$vector%*%t(a3$vector)
A4 = A3 - Z3
a4 = p.m(A4) ; a4
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.119.232.236
1F：推 farte::) 金厉害 05/24 21:42