昨天才回答過這個問題..你可以再搜搜的Gauss消去法的分析。其包括兩個過程：消去過程：把方程組係數矩陣A化為同解的上三角矩陣；回代過程：按相反的順序，從xn至x1逐個求解上三角方程組。%高斯消去法的MATLAB程式function x=gauss(a,b); %編寫高斯消去法函式%a表示方程組的係數矩陣，b表示方程組的值%X表示最終的輸出結果，即方程組的解n=length(b); %計算方程組的維數%下面的程式在不斷的消去，直到變成a變成上三角矩陣未知for k=1:n-1for i=k+1:na(i,k)=a(i,k)/a(k,k);for j=k+1:na(i,j)=a(i,j)-a(i,k)*a(k,j);endb(i)=b(i)-a(i,k)*b(k);endend%表示高斯消去法的迴帶過程x=zeros(n,1);x(n)=b(n)/a(n,n);for k=n-1:-1:1s=b(k);for j=k+1:ns=s-a(k,j)*x(j);endx(k)=s/a(k,k);end例項驗證： %呼叫編好的消去法函式A=[1,2,3;2,2,3;-1,-3,10];B=[0,3,2];gauss(A,B)ans =3.0000-1.55170.0345A=[1,2,3;2,2,3;-1,-3,10];B=[0,3,2];x=gauss(A,B)x =3.0000-1.55170.0345A*x %反代求解進行比較ans =0.00003.00002.0000
昨天才回答過這個問題..你可以再搜搜的Gauss消去法的分析。其包括兩個過程：消去過程：把方程組係數矩陣A化為同解的上三角矩陣；回代過程：按相反的順序，從xn至x1逐個求解上三角方程組。%高斯消去法的MATLAB程式function x=gauss(a,b); %編寫高斯消去法函式%a表示方程組的係數矩陣，b表示方程組的值%X表示最終的輸出結果，即方程組的解n=length(b); %計算方程組的維數%下面的程式在不斷的消去，直到變成a變成上三角矩陣未知for k=1:n-1for i=k+1:na(i,k)=a(i,k)/a(k,k);for j=k+1:na(i,j)=a(i,j)-a(i,k)*a(k,j);endb(i)=b(i)-a(i,k)*b(k);endend%表示高斯消去法的迴帶過程x=zeros(n,1);x(n)=b(n)/a(n,n);for k=n-1:-1:1s=b(k);for j=k+1:ns=s-a(k,j)*x(j);endx(k)=s/a(k,k);end例項驗證： %呼叫編好的消去法函式A=[1,2,3;2,2,3;-1,-3,10];B=[0,3,2];gauss(A,B)ans =3.0000-1.55170.0345A=[1,2,3;2,2,3;-1,-3,10];B=[0,3,2];x=gauss(A,B)x =3.0000-1.55170.0345A*x %反代求解進行比較ans =0.00003.00002.0000