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

PAGE \* MERGEFORMAT PAGE \* MERGEFORMAT 1 PAGE \* MERGEFORMAT PAGE \* MERGEFORMAT 1 第三章上机习题 用你所熟悉的的计算机语言编制利用QR分解求解线性方程组和线性最小二乘问题的通用子程序，并用你编制的子程序完成下面的计算任务： 求解第一章上机习题中的三个线性方程组，并将所得的计算结果与前面的结果相比较，说明各方法的优劣； 求一个二次多项式，使得在残向量的2范数下最小的意义下拟合表3.2中的数据； 表 3.2 ti -1 -0.75 -0.5 0 0.25 0.5 0.75 yi 1 0.8125 0.75 1 1.3125 1.75 2.3125 在房产估价的线性模型 中，分别表示税、浴室数目、占地面积、车库数目、房屋数目、居室数目、房龄、建筑类型、户型及壁炉数目，代表房屋价格。现根据表3.3和表3.4给出的28组数据，求出模型中参数的最小二乘结果。 （表3.3和表3.4见课本P99-100） 解 分析： 计算一个Householder变换H： 由于，则计算一个Householder变换H等价于计算相应的。其中。 在实际计算中， 为避免出现两个相近的数出现的情形，当时，令； 为便于储存，将规格化为，相应的，变为 为防止溢出现象，用代替 QR分解： 利用Householder变换逐步将转化为上三角矩阵，则有 ，其中，。 在实际计算中，从，若，依次计算对应的 即对应的，，将储存到，储存到，迭代结束后再次计算，有，（时） 求解线性方程组或最小二乘问题的步骤为 = 1 \* roman \* MERGEFORMAT i 计算的QR分解； = 2 \* roman \* MERGEFORMAT ii 计算，其中 = 3 \* roman \* MERGEFORMAT iii 利用回代法求解上三角方程组 （4）对第一章第一个线性方程组，由于R的结果最后一行为零，故使用前代法时不计最后一行，而用运行结果计算。 运算matlab程序为 1 计算Householder变换 [v,belta]=house(x) function [v,belta]=house(x) n=length(x); x=x/norm(x,inf); sigma=x(2:n)*x(2:n); v=zeros(n,1); v(2:n,1)=x(2:n); if sigma==0 belta=0; else alpha=sqrt(x(1)^2+sigma); if x(1)=0 v(1)=x(1)-alpha; else v(1)=-sigma/(x(1)+alpha); end belta=2*v(1)^2/(sigma+v(1)^2); v=v/v(1,1); end end 2 计算的QR分解 [Q,R]=QRfenjie(A) function [Q,R]=QRfenjie(A) [m,n]=size(A); Q=eye(m); for j=1:n if jm [v,belta]=house(A(j:m,j)); H=eye(m-j+1)-belta*v*v; A(j:m,j:n)=H*A(j:m,j:n); d(j)=belta; A(j+1:m,j)=v(2:m-j+1); end end R=triu(A(1:n,:)); for j=1:n if jm H=eye(m); temp=[1;A(j+1:m,j)]; H(j:m,j:m)=H(j:m,j:m)-d(j)*temp*temp; Q=Q*H; end end end 3 解下三角形方程组的前代法 x=qiandaifa(L,b) function x=qiandaifa(L,b) n=length(b); for j=1:n-1 b(j)=b(j)/L(j,j); b(j+1:n)=b(j+1:n)-b(j)*L(j+1:n,j); end b(n)=b(n)/L(n,n); x=b; end 4 求解第一章上机习题中的三个线性方程组 ex3_1 clear;clc; %第一题 A=6*eye(84)+diag(8*ones(1,83),-1)+diag(ones(1,