对此感到困惑:我正在计算以下矩阵 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