ch5_3.java

Java环境下实现的数值分析程序

class ch5_3{
	public static void main(String[]args){
		double a[][]={{2,4},{3,-5},{10,-12},{4,11}};
		double b[]={10,-13,-26,25};
		double c0=0,c1=0,c2=0,c3=0,c4=0,c5=0;
		for(int i=0;i<4;i++){
			c0=c0+Math.pow(a[i][0],2);
			c1=c1+a[i][0]*a[i][1];
			c2=c2+a[i][0]*a[i][1];
			c3=c3+Math.pow(a[i][1],2);
			c4=c4+a[i][0]*b[i];
			c5=c5+a[i][1]*b[i];			
		}
		
		double[][] c={{c0,c1,c4},{c2,c3,c5}}; //创建增广矩阵
        double temp=0,x1,x2;
        int i,j,k,m,n,p;
        for(i=0,k=i;i<2;i++){
            double item=c[i][i];
            for(j=i+1;j<2;j++){     //找寻第i列中位于对角线及其以下元素绝对值中的最大者 
                if(item<Math.abs(c[j][i])){
                  item=c[j][i];
                  k=j;
                }
            }
            for(p=i;p<3;p++){       //将第k行与第i行互换
               temp=c[i][p];
               c[i][p]=c[k][p];
               c[k][p]=temp;
            }
            for(m=i+1;m<2;m++){     //高斯消去，将矩阵化为上三角形
        	    double key=c[m][i]/c[i][i];
                for(n=i;n<3;n++){
                    c[m][n]=c[m][n]-key*c[i][n];
                }
            }
        }
    
        x2=c[1][2]/c[1][1];      //求解同解三角矩阵
        x1=(c[0][2]-x2*c[0][1])/c[0][0];
            
        if(x1<0)System.out.println("x1="+x1); //这里的if判断仅仅是为了输出时较好看
        else System.out.println("x1= "+x1);
        if(x2<0)System.out.println("x2="+x2);
        else System.out.println("x2= "+x2);
	}
}