数值线性代数第二版徐树方高立张平文上机习题第三章实验报告.pdf

WORD格式整理版 第三章上机习题 用你所熟悉的的计算机语言编制利用QR 分解求解线性方程组和线性最小二乘问题的通 用子程序，并用你编制的子程序完成下面的计算任务： （1）求解第一章上机习题中的三个线性方程组，并将所得的计算结果与前面的结果相比较， 说明各方法的优劣； （2）求一个二次多项式y=at +bt+c2 ，使得在残向量的2 范数下最小的意义下拟合表3.2 中的数据； 表 3.2 t -1 -0.75 -0.5 0 0.25 0.5 0.75 i y 1 0.8125 0.75 1 1.3125 1.75 2.3125 i （3）在房产估价的线性模型 yx ax a x a x 0 1 1 2 2 11 11 中，a,a ,,a 分别表示税、浴室数目、占地面积、车库数目、房屋数目、居室数目、房 1 2 11 龄、建筑类型、户型及壁炉数目， 代表房屋价格。现根据表3.3 和表3.4 给出的28 组数y 据，求出模型中参数的最小二乘结果。 （表3.3 和表3.4 见课本P99-100） 解 分析： （1）计算一个Householder 变换H： T  T  由于H I2ww I vv ，则计算一个Householder 变换H 等价于计算相应的 、v。  T 其中vx||x|| e, 2/(v v)。 2 1 在实际计算中， 2 2 -(x  x )  为避免出现两个相近的数出现的情形，当x 0时，令v  2 n ； 1 1 x ||x|| 1 2 v   2 T 为便于储存，将 规格化为vv/v ，相应的， 变为 2v /(v v) 1 1 为防止溢出现象，用x/||x|| 代替  （2）QR 分解： 利用Householder 变换逐步将A ,mn转化为上三角矩阵H H H A，则有 mn n n1 1 优质.参考.资料 WORD格式整理版 R   AQ ，其中QH H H ，R(1:n,:)。  0 1 2 n ~ 在实际计算中，从j 1:n，若j m，依次计算xA((j :m