理学试验直接法求解线性方程组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如如果果则则转转(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+) /*行乘数行乘数*/ 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,