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

????..A..H.???????162??2010-06-27?11:11??數值分析上機1答案\~$一.doc

?????目錄??????????0??2010-06-27?11:12??數值分析上機1答案

?????文件??????68096??2009-11-24?10:29??數值分析上機1答案\一.doc

-----------?---------??----------?-----??----

????????????????69228????????????????????5