[数值分析] 利用QR分解求解最小二乘法的原理

当方程阶数较大的时,法方程的条件数变得非常大,需要使用QR分解来求解。
Ax=b

    [Q,R] = qr(A,0); % 
    x  = R\(Q'*b); % 

原理:
在这里插入图片描述
因为 c 2 2 2 \|c_2\|_2^2为常数, 所以上述方程可以在 R x = c 1 Rx=c_1 的时候取得极值,是为最小二乘法的解。
x = R 1 ( Q T b ) x=R^{-1} (Q^Tb) (利用精简QR分解得到的R)

猜你喜欢

转载自blog.csdn.net/huangzhywin/article/details/89205696