C++高斯全主元法计算矩阵.docVIP

C++高斯全主元法计算矩阵 //解线性方程组 #includeiostream #includeiomanip #includestdlib.h using namespace std; const int Number=15; double a[Number][Number],b[Number]; int A_y[Number]; int lenth,copy_lenth; double a_sum; char * x; void input(); void guass_all(); int Doolittle_check(double a[][Number],double b[Number]); void xiaoqu_u_l(); void calculate_u_l(); double calculate_A(int n,int m); double quanpailie_A(); void exchange(int m,int i); void exchange_lie(int j); void exchange_hang(int m,int n); void gauss_all_xiaoqu(); void gauss_calculate(); void exchange_a_lie(int m,int n); void exchange_x(int m,int n); //主函数 void main() { int flag=1; input(); //输入方程 while(flag) { guass_all(); break; } system(pause);/*保持当前画面*/ } void input() { int i,j; cout请输入方程的个数:; cinlenth; if(lenthNumber) { coutIt is too big.\n; return; } x=new char[lenth]; for(i=0;ilenth;i++) x[i]=a+i; //输入方程矩阵 //提示如何输入 cout请在每个方程里输入lenth系数和一个常数:\n; cout例：\n方程:a; for(i=1;ilenth;i++) { cout+i+1x[i]; } cout=10\n; cout应输入:; for(i=0;ilenth;i++) couti+1 ; cout10\n; //输入每个方程 for(i=0;ilenth;i++) { cout输入方程i+1:; for(j=0;jlenth;j++) cina[i][j]; cinb[i]; } } double calculate_A(int n,int m) //计算行列式 { int i; if(n==1) { a_sum+= quanpailie_A(); } else{for(i=0;in;i++) { exchange(m,m+i); calculate_A(n-1,m+1); exchange(m,m+i); } } return a_sum; } double quanpailie_A() //计算行列式中一种全排列的值 { int i,j,l; double sum=0,p; for(i=0,l=0;ilenth;i++) for(j=0;A_y[j]!=ijlenth;j++) if(A_y[j]i) l++; for(p=1,i=0;ilenth;i++) p*=a[i][A_y[i]]; sum+=p*((l%2==0)?(1):(-1)); return sum; } void guass_all() { int i,j; gauss_all_xiaoqu(); for(i=0;ilenth;i++) { for(j=0;jlenth;j++) coutsetw(10)setprecision(5)a[i][j]; coutsetw(10)b[i]endl; } if(a[lenth-1][lenth-1]!=0) { cout系数行列式不为零,方程有唯一的解：\n; gauss_c