Java的高斯消元法

算法虽然会，但是用数组的形式进行求值是一件非常非常恶心的事情。

更恶心的是，公司里面的逻辑居然要用到方程组求解，万分无奈之下，上网找了一份Java代码。

谁知道居然是个算法数组从1开始的。

好吧，认了吧，我就把你改过来。

调试通过。

留个代码做纪念。

 packagecom;

 publicclassEquationSolver{

 /**
 *@列主元高斯消去法
 */
 staticdoublea[][];
 staticdoubleb[];
 staticdoublex[];
 staticintn;
 staticintn2;//记录换行的次数
 publicstaticvoidElimination(){//消元
 PrintA();
 for(intk=0;k
 {
 Wrap(k);
 for(inti=k+1;i
 {
 doublel=a[i][k]/a[k][k];
 a[i][k]=0.0;
 for(intj=k+1;j
 a[i][j]=a[i][j]-l*a[k][j];
 b[i]=b[i]-l*b[k];
 }
 System.out.println("第"+k+"次消元后：");
 PrintA();
 }

 }
 publicstaticvoidBack()//回代
 {
 x[n-1]=b[n-1]/a[n-1][n-1];
 for(inti=n-2;i>=0;i--)
 x[i]=(b[i]-jisuan(i))/a[i][i];
 }
 publicstaticdoublejisuan(inti){
 doublehe=0.0;
 for(intj=i;j<=n-1;j++)
 he=he+x[j]*a[i][j];
 returnhe;
 }
 publicstaticvoidWrap(intk){//换行
 doublemax=Math.abs(a[k][k]);
 intn1=k;//记住要交换的行
 for(inti=k+1;i//找到要交换的行
 {
 if(Math.abs(a[i][k])>max){
 n1=i;
 max=Math.abs(a[i][k]);
 }
 }
 if(n1!=k)
 {
 n2++;
 System.out.println("当k="+k+"时,要交换的行是："+k+"和"+n1);
 for(intj=k;j//交换a的行
 {
 doublex1;
 x1=a[k][j];
 a[k][j]=a[n1][j];
 a[n1][j]=x1;
 }
 doubleb1;//交换b的行
 b1=b[k];
 b[k]=b[n1];
 b[n1]=b1;
 System.out.println("交换后：");
 PrintA();
 }
 }
 publicstaticvoidDeterminant(){//求行列式
 doubleDM=1.0;
 for(inti=0;i
 {
 doublea2=a[i][i];
 DM=DM*a2;
 }
 doublen3=(double)n2;
 DM=DM*Math.pow(-1.0,n3);
 System.out.println("该方程组的系数行列式：detA="+DM);
 }
 publicstaticvoidPrintA(){//输出增广矩阵
 System.out.println("增广矩阵为：");
 for(inti=0;i
 {
 for(intj=0;j
 System.out.print(a[i][j]+"");
 System.out.print(b[i]+"");
 System.out.print("\n");
 }
 }
 publicstaticvoidPrint(){//输出方程的根
 System.out.println("方程组的根为：");
 for(inti=0;i
 System.out.println("x"+i+"="+x[i]);
 }
 publicstaticvoidmain(String[]args){
 //Scanneras=newScanner(System.in);
 System.out.println("输入方程组的元数：");
 //n=as.nextInt();
 n=2;
 a=newdouble[n][n];
 b=newdouble[n];
 x=newdouble[n];

 doubleinputA[][]={{1,1},{1,-1}};
 a=inputA;

 doubleinputB[]={2,0};
 b=inputB;

 doubleinputX[]={1,1};

 System.out.println("输入方程组的系数矩阵a：");
 //for(inti=1;i<=n;i++)
 //for(intj=1;j<=n;j++)
 //a[i][j]=as.nextDouble();
 System.out.println("输入方程组矩阵b：");
 //for(inti=1;i<=n;i++)
 //b[i]=as.nextDouble();
 Elimination();
 Back();
 Print();
 Determinant();
 }
 }