D <- matrix(
c(2,1,3,1,1,1,0,-3,3,-2,2,-6),
nrow=3
)
print(D)
print(svd(D))
round(svd(D)$d, 5)
J <- matrix(
c(1,2,3,4,2,5,6,7,3,6,8,9,4,7,9,10),
nrow=4
)
print(J)
eigen(J)$values
W = 1/16 * matrix(
c(11,-5,-3,1,4,2,
-5,11,-3,1,4,2,
-3,-3,3,-1,-4,-2,
1,1,-1,11,-4,6,
4,4,-4,-4,8,0,
2,2,-2,6,0,4),
nrow=6
)
print(W)
print(W%*%W)
identical(W%*%W, W)
(11+11+3+11+8+4)/16
J <- matrix(
c(1,2,3,4,2,5,6,7,3,6,8,9,4,7,9,10),
nrow=4
)
print(J)
B = diag(c(1,2,3,4))
print(B)
print(solve(B) %*% J)
v <- eigen(solve(B) %*% J)
print(v)
a <- sum(v$vectors[,1]^2 * c(1,2,3,4))
x <- v$vectors[,1] / sqrt(a)
print(abs(x))