1

私はマトリックスの世界に不慣れです。この基本的な質問で申し訳ありませんが、私は理解できませんでした：

私は4つのマトリックスを持っています（1つは不明です）。

マトリックスX

x <- c(44.412, 0.238, -0.027, 93.128, 0.238, 0.427, -0.193, 0.673, 0.027, 
     -0.193, 0.094, -0.428, 93.128, 0.673, -0.428, 224.099)

X <- matrix(x, ncol = 4 )

行列 B : 解決する必要がある、1 X 4 (列 x 行)、b1、b2、b3、b4 値

マトリックスG

g <- c(33.575, 0.080, -0.006, 68.123, 0.080, 0.238, -0.033, 0.468, -0.006, 
-0.033, 0.084, -0.764, 68.123, 0.468, -0.764, 205.144)

G <- matrix(g, ncol = 4)

マトリックス A

a <- c(1, 1, 1, 1) # one this case but can be any value 
A <- matrix(a, ncol = 1)

解決：

B = inv(X) G A  # inv(X) is inverse of the X matrix multiplied by G and A

これを適切に解決する方法、特に逆行列を解決する方法がわかりませんでした。あなたの助けに感謝。

4

2 に答える 2

4

ニックとベンはどちらも教師であり、他の人の宿題をするよりもさらに多くのことをしていると思いますが、完全な解決策への道は本当に明白だったので、そうしないことはあまり意味がありませんでした次のステップ：

B = solve(X) %*% G %*% A 
> B
             [,1]
[1,] -2.622000509
[2,]  7.566857261
[3,] 17.691911600
[4,]  2.318762273

反転のQRメソッドは、2番目の引数として単位行列を指定することで呼び出すことができます。

> qr.solve(G, diag(1,4))
                [,1]             [,2]          [,3]             [,4]
[1,]  0.098084556856 -0.0087200426695 -0.3027373205 -0.0336789016478
[2,] -0.008720042669  4.4473233701790  1.7395207242 -0.0007717410073
[3,] -0.302737320546  1.7395207241703 13.9161591761  0.1483895429511
[4,] -0.033678901648 -0.0007717410073  0.1483895430  0.0166129089935
于 2011-11-09T13:07:50.207 に答える
3

より計算上安定した解決策は、qrではなくを使用することですsolve。

method1 <- solve(X) %*% G %*% A
method2 <- qr.coef(qr(X), G) %*% A
stopifnot(isTRUE(all.equal(method1, method2)))

の例を参照してください?qr。

于 2011-11-09T15:24:06.400 に答える