理學試驗直接法求解線性方程組PPT學習教案

1、會計學1 理學試驗直接法求解線性方程組理學試驗直接法求解線性方程組 第1頁/共7頁 1 2 1 2 2 : . ( , )max( , ) .( , ), , ( , )( , )/( , ), ( )( )/( , ) k i n kN A p kA i k A p kTOL pk A k jA k jA k k b kb kA k k 對對于于 選選主主元元 如如果果則則停停止止 2.32.3如如果果則則轉(zhuǎn)轉(zhuǎn)(2.4),(2.4),否否則則 A(k,k:N)A(p,k:N),b(k)b(p)A(k,k:N)A(p,k:N),b(k)b(p) 2.42.4對對于于j=k+1:Nj=k+1:N

2、 第2頁/共7頁 2 4 111 . ( ,:)( ,:)( , )*( ,:) ( )( )( , )* ( ) A i kNA i kNA i kA k kN b ib iA i kb k 對對于于i=k+1:Ni=k+1:N 3 1 . - ( )( )( ,i+1:)* (:)b ib iA iNb iN 回回代代求求解解 對對于于i=Ni=N： 1:11:1 1()k kk a 對對于于一一些些特特殊殊類類型型的的方方程程組組，可可以以保保證證不不會會很很小小， 因因而而不不需需要要選選主主元元 第3頁/共7頁 #include #include #include #include

3、#define N 4 #define EPS 1e-5 void Guass(double AN+1) int i,j,k=0,I,t; double w=0,p,b; for(k=0;kN-1;k+) /*選主元選主元*/ p=Akk; I=k; for(i=k+1;ifabs(p) p=Aik; I=i; if(fabs(p)EPS) printf(the number is error! i =%d, k=%dn,i,k); exit(0); /*換行換行*/ if(I!=k) for(j=0;jN+1;j+) b=Akj; Akj=AIj; AIj=b; 第4頁/共7頁 for(i=

4、k+1;iN;i+) /*行乘數(shù)行乘數(shù)*/ Aik=Aik/Akk; for(j=k+1;jN+1;j+) Aij=Aij-Aik*Akj; /*輸出變換后的矩陣輸出變換后的矩陣*/ printf(The Matrix:n); for(i=0;iN;i+) for(j=0;jN+1;j+) printf(A%d%d=%8.4f, i,j, Aij); printf(n); if(fabs(AN-1N-1)=0;k-) w=0; for(j=k+1;jN;j+) w=w+Akj*AjN; AkN=AkN-w; AkN=AkN/Akk; 第5頁/共7頁 /*方程組的解方程組的解*/ printf(the root:n); for(i=0;iN;i+) printf(x%d=%fn,i,AiN) ; void main() double ANN+1=/1, -1, 2, -1, -16,2, -2, 3,