QR分解解决列满秩最小二乘问题matlab程序

以下是用QR分解解决列满秩最小二乘问题的MATLAB代码：

function [x,resnorm] = lsqr_qr(A,b)
%列满秩最小二乘问题的QR分解解法
[Q,R] = qr(A,0); %QR分解
x = R\(Q'*b); %解方程
resnorm = norm(b - A*x); %计算残差范数
end

在这个程序中，我们使用MATLAB内置的`qr`函数进行QR分解，然后利用QR分解求解最小二乘问题。函数返回的$x$是最小二乘解，`resnorm`是残差范数。注意，我们在调用`qr`函数时，第二个输入参数为0，表示只对$A$进行经典的（非紧凑的）QR分解，得到的$Q$和$R$矩阵的维度与$A$相同。如果第二个输入参数为1，则进行紧凑的QR分解，得到的$Q$和$R$矩阵的维度可能会比$A$小。