matlab解线性方程组的QR分解法

解线性方程组常见的QR分解法，尤其是针对大型矩阵，比较实用

function[Hy]=holder（xk）
%建立 阵，将x向量转换为k值确定的y向量
%x、k作为已知量输入
n=length（x）;    %确认阶数n
M=max（abs（x））;
if M==0
    disp（‘M=0‘）;
    return;
else z=x/M;
end;          %规范化，防止x分量过大而溢出
s=0;
for i=k:n
    s=s+z（i）*z（i）;
end;
s=sqrt（s）;
if z（k）<0
    s=-s;
end;          %求 
p=s*（s+z（k））;  %求 
U=z;
for j=1:k-1
    U（j）=0;
end;
U（k）=s+z（k）;
H=eye（nn）-p\U*U‘;   %得 
y=zeros（n1）;
for j=1:k-1
    y（j）=x（j）
end;
y（k）=-M*s;          %得y向量

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        503  2009-11-24 09:03  数值分析上机1答案\holder.m

     文件        467  2009-11-24 09:59  数值分析上机1答案\qrholder.m

     目录          0  2010-06-27 11:12  数值分析上机1答案

----------- ---------  ---------- -----  ----

                  970                    3