0

对此感到困惑:我正在计算以下矩阵 X 的 QR 分解,我希望通过计算乘积 QR 来返回 X。好吧,该产品给出了“几乎”X,但是第 1 列和第 2 列倒置了……这怎么可能?

> X
     [,1]  [,2] [,3]
[1,] 1-0i -1+1i 1-1i
[2,] 1-1i -1+2i 1-1i
[3,] 1-1i -1+1i 1+0i
> 
> QRX <- qr(X)
> X <- qr.X(QRX)
> Q <- qr.Q(QRX)
> R <- qr.R(QRX)
> 
> Q%*%R
      [,1] [,2] [,3]
[1,] -1+1i 1-0i 1-1i
[2,] -1+2i 1-1i 1-1i
[3,] -1+1i 1-1i 1+0i
4

0 回答 0