基于MATLAB的數值分析编程上机作业.docVIP

PAGE 基于MATLAB的数值分析编程上机作业（1） 班级： 学号： 姓名： 日期： 页数： 基于MATLAB的数值分析编程上机作业（1） PAGE 3 一、给定向量x≠0，计算初等反射阵H。 1、基本原理： 若的分量不全为零，则由 确定的镜面反射阵H使得；当时，由 有 算法1： 输入x 将x规范化， 如果M＝0，则转出停机 否则 计算，如果，则 4、 5、计算 6、 7、 8、按要求输出 结束 2、程序代码： 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。 %使用举例： % [p,u]=holder2(x) %Define variables: % x － 输入的n维向量； % n － n维向量x的维数； % M － M是向量x的无穷范数，即x中绝对值最大的一项的绝对值； %