example2_4.java

清华大学2002年出版的《科学与工程数值计算算法Java》配套源码

import java.applet.*;
import java.awt.*;
public class Example2_4 extends Applet
{
public boolean Mdzinv(int n , double[][] a)
   { int i,j,k,m;
     double  w,g;
     double[] b=new double[n];
    // 循环重新编号求解
    for(k=0;k<n;k++)
   { w = a[0][0];
    //求解失败
    if (Math.abs(w)==0) return false;
    m=n-k-1;
    for(i=1;i<n;i++)
    { g = a[i][0];
      b[i] = g/w;
      if (i <= m)b[i] = -b[i];
      for(j=1;j<=i;j++)a[i-1][j-1]= a[i][j] + g * b[j];
    }
      a[n-1][n-1]= 1.0/ w;
      for(i=1;i<n;i++) a[n-1][i-1]= b[i];

   }
    for(i=0;i<n-1;i++)
      for(j=i+1;j<n;j++)
        a[i][j]=a[j][i];
//求解成功
return true;
}

public void Mrcheng(double[][] a,double[][] b,double[][]c,int m,int n,int l)
  {double[][] d=new double[m][l];
   int i,j,k;
   for(i=0;i<m;i++)
     for(j=0;j<l;j++)
      {d[i][j]=0;
       for(k=0;k<n;k++)d[i][j]+=a[i][k]*b[k][j];}
  for(i=0;i<m;i++)
   for(j=0;j<l;j++)
    c[i][j]=d[i][j];
}

public void paint(Graphics g)
{double[][] a={{5,7,6,5},{7,10,8,7},{6,8,10,9},{5,7,9,10}};
 double[][] b=new double[4][4];
  int i,j;
  for(i=0;i<4;i++)
   for(j=0;j<4;j++)
    b[i][j]=a[i][j];
  g.drawString("A=",10,20);
  for(i=0;i<4;i++)
  g.drawString(""+a[i][0]+"  "+a[i][1]+"  "+a[i][2]+"  "+a[i][3],10,40+20*i);
  if(Mdzinv(4,a))
  {g.drawString("A的逆阵B为",10,120);
    for(i=0;i<4;i++)
    g.drawString(""+a[i][0]+"  "+a[i][1]+"  "+a[i][2]+"  "+a[i][3],10,140+20*i);
   Mrcheng(a,b,a,4,4,4);
   g.drawString("A*B=",10,220);
   for(i=0;i<4;i++)
   g.drawString(""+a[i][0]+"  "+a[i][1]+"  "+a[i][2]+"  "+a[i][3],10,240+20*i);
  }
  else
   g.drawString("求解失败,原矩阵不是对称正定阵",10,120);
 }
}