下载后可任意编辑一、给定向量 x≠0，计算初等反射阵 Hk。1．程序功能： 给定向量 x≠0，计算初等反射阵 Hk。2．基本原理： 若的重量不全为零，则由确定的镜面反射阵 H 使得；当时，由有算法：(1)输入 x，若 x 为零向量，则报错(2)将 x 法律规范化，假如 M＝0，则报错同时转出停机否则(3)计算，假如，则(4)(5)计算(6)下载后可任意编辑(7)(8)按要求输出,结束3．变量说明：x－输入的 n 维向量；n－n 维向量 x 的维数；M－M 是向量 x 的无穷范数，即 x 中绝对值最大的一项的绝对值；p －Householder 初等变换阵的系数 ρ；u －Householder 初等变换阵的向量 Us－向量 x 的二范数；x－输入的 n 维向量；n－n 维向量 x 的维数；p －Householder 初等变换阵的系数 ρ；u －Householder 初等变换阵的向量 Uk－数 k，H*x=y，使得 y 的第 k+1 项到最后项全为零；4．程序代码：(1)function [p,u]=holder2(x)%HOLDER2 给定向量 x≠0,计算 Householder 初等变换阵的 p,u%程序功能：函数 holder2 给定向量 x≠0,计算 Householder 初等变换阵的 p,u；%输入：n 维向量 x；%输出：[p,u]。p 是 Householder 初等变换阵的系数 ρ，% u 是 Householder 初等变换阵的向量 U。n=length(x); % 得到 n 维向量 x 的维数；p=1;u=0; % 初始化 p,u；M=max(abs(x)); % 得到向量 x 的无穷范数，即 x 中绝对值最大的一项的绝对值；if M==0 % 假如 x=0，提示出错,程序终止; disp('Error: M=0'); return;else x=x/M; % 法律规范化end;s=norm(x); % 求 x 的二范数if x(1)<0 % 首项为负，s 值要变号 s=-s;endu=x; % 除首项外，其余各项 x,u 相同u(1)=s+x(1); % 计算 u 的首项p=s*u(1); % 计算 pif n==1 u=0; end % 若 x 是 1×1 维向量，则 u=0(2)下载后可任意编辑function H=holderk(x,k)%HOLDERK 给定向量 x≠0,数 k，计算初等反射阵 Hk，使 HkX=Y，其中 Y 的第 k+1 项到最后项全为零；%程序功能：函数 holderk 给定向量 x≠0，数 k，计算初等反射阵 Hk，使 HkX=Y，%程序功能：函数 holder2 给定向量 x≠0,计算 Householder 初等变换阵的 p,u；%输入：n 维向量 x,数 k；%输出：H。H 是 Householder 初等变换阵，H*x=y，使得 y 的第 k+1 项到最后项全为零；%引用函数：holder2； n=length(x); % 得到 n 维向...
