calulation math.zip_highwayqpj_paintm4m_数值线性代数 上机 答案

%习题1.3 %第一问 A=10*eye(100)+diag(ones(1,99),-1)+diag(ones(1,99),1); b=round(100*rand(100,1)); %不选主元Gauss消去法 [L,U]=GaussLA(A); x1_1=Gauss(A,b,L,U); %列主元Gauss消去法 [L,U,P]=GaussCol(A); x1_2=Gauss(A,b,L,U,P); %平方根法 L=Cholesky(A); x1_3=Gauss(A,b,L,L'); %改进的平方根法 [L,D]=LDLt(A); x1_4=Gauss(A,b,L,D*L'); %第二问 A=hilb(40); b=sum(A); b=b'; %不选主元Gauss消去法 [L,U]=GaussLA(A); x2_1=Gauss(A,b,L,U); %列主元Gauss消去法 [L,U,P]=GaussCol(A); x2_2=Gauss(A,b,L,U,P); %平方根法 L=Cholesky(A); x2_3=Gauss(A,b,L,L'); %改进的平方根法 [L,D]=LDLt(A); x2_4=Gauss(A,b,L,D*L'); %解的比较 y1=1:100;y2=1:40; subplot(4,2,1);plot(y1,x1_1);title('Gauss1'); subplot(4,2,2);plot(y1,x1_2);title('PGauss1'); subplot(4,2,3);plot(y1,x1_3);title('平方根法1'); subplot(4,2,4);plot(y1,x1_4);title('改进的平方根法1'); subplot(4,2,5);plot(y2,x2_1);title('Gauss2'); subplot(4,2,6);plot(y2,x2_2);title('PGauss2'); subplot(4,2,7);plot(y2,x2_3);title('平方根法2'); subplot(4,2,8);plot(y2,x2_4);title('改进的平方根法2');