计算方法上机练习题.docxVIP

计算方法上机练习题信息44 肖卓昊 2140502094习题2-13：程序：A=[0.4721 0.2352 -0.2613 0.8421; 0.2352 0.7411 -0.0463 0.1569; -0.2613 -0.0463 -0.8955 0.1748; 0.8421 0.1569 0.1748 0.9841];b=[-0.2317;0.3219;0.6217;0.9838];% 不选主元的高斯消去法x=Gauss(A,b);x=x;disp(不选主元的解x=)x% 选主元的高斯消去法x=Gauss_pivot(A,b);disp(选主元的解x=)xfunction [x]=Gauss_pivot(A,b)[n,n]=size(A);x=zeros(n,1);C=zeros(1,n);B=[A b];for i=1:n-1 [max1,j]=max(abs(B(i:n,i))); C=B(i,:); B(i,:)=B(j+i-1,:); B(j+i-1,:)=C; if A(i,i)==0 break end for k=i+1:n; temp=B(k,i)/B(i,i); B(k,i+1:n+1)=B(k,i+1:n+1)-temp*B(i,i+1:n+1); end endx(n)=B(n,n+1)/B(n,n);for k=n-1:-1:1 x(k)=(B(k,n+1)-B(k,k+1:n)*x(k+1:n))/B(k,k);endfor s=1:1:n-1; for i=s+1:1:n c=B(i,s)/B(s,s); B(i,:)=B(i,:)-c*B(s,:); endenddisp(选主元的高斯消去后的增广矩阵A_b=)Bfunction [x]=Gauss(A,b)A_b=[A,b];[m,n]=size(A_b);for k=1:1:m-1; for i=k+1:1:m c=A_b(i,k)/A_b(k,k); A_b(i,:)=A_b(i,:)-c*A_b(k,:); endenddisp(不选主元的高斯消去后的增广矩阵A_b=)A_bx(m)=A_b(m,n)/A_b(m,m);for k=m-1:-1:1 x(k)=(A_b(k,n)-A_b(k,k+1:m)*x(k+1:m))/A_b(k,k);end运行结果：（1）不选主元A_b = 0.4721 0.2352 -0.2613 0.8421 -0.2317 0 0.6239 0.0839 -0.2626 0.4373 0 0 -1.0514 0.6762 0.4347 0 0 0 -0.1936 1.8607解x = 14.2241 -2.4575 -6.5933 -9.6090选主元A_b= 0.8421 0.1569 0.1748 0.9841 0.9838 0.0000 0.6535 -0.1439 -0.3928 -0.2277 -0.0000 0 -0.7754 0.8166 1.2499 0.0000 0 0.0000 -0.8413 -2.3816解x= 14.2241 -2.4575 -6.5933 -9.6090（2）不选主元A_b= 0.5398 0.7161 -0.5554 -0.2982 0.2058 0 -0.0050 0.8974 -0.3351 -0.2507 0 0 -300.6974 112.9461 84.0065 0 0 0 -0.0835 -0.0866解x= 0.5162 0.4152 0.1100 1.0365选主元A_b= 0.6465 -0.8187 -0.1872 0.1291 0.1070 0 2.4125 -0.4412 -0.6364 -0.0066 0 0 1.5249 0.1200 -0.