chaoding.java.bak

数值分析之计算方法与实验5：曲线拟合 by java

//解超定方程组
class Chaoding{
	static int row=0,line=0;
	public static void main(String[] args){
		double[][] ab={{2,4,10},{3,-5,-13},{10,-12,-26},{4,11,25}}; 
		
		line=ab[0].length;row=line-1;
		double[][] h=new double[row][line];
		for(int i=0;i<row;i++){
			for(int j=0;j<row;j++){
				for(int k=0;k<ab.length;k++)h[i][j]+=ab[k][i]*ab[k][j];
			}
		}
		for(int i=0;i<row;i++){
			for(int k=0;k<ab.length;k++)h[i][line-1]+=ab[k][i]*ab[k][line-1];
		}
			
		double[] x=directLU(h);
		for(int i=0;i<row;i++)System.out.print("x"+(i+1)+"="+x[i]+"  ");
		}
	static double[] directLU(double[][] ab){
		double[][] u=new double[row][line];
		double[][] l=new double[row][row];
		for(int j=0;j<line;j++)u[0][j]=ab[0][j];
		for(int i=1;i<row;i++)l[i][0]=ab[i][0]/ab[0][0];
		for(int k=1;k<row;k++){
			for(int j=k;j<line;j++){
				u[k][j]=ab[k][j];
				for(int q=0;q<=k-1;q++){
					u[k][j]=u[k][j]-u[q][j]*l[k][q];
				}
			}
			for(int i=k+1;i<row;i++){
				l[i][k]=u[k][i]/u[k][k];
			}
		}
		double[] x=new double[row];
		x[row-1]=u[row-1][line-1]/u[row-1][line-2];
		for(int i=row-2;i>=0;i--){
			x[i]=u[i][line-1];
			
			for(int j=i+1;j<=line-2;j++){
				x[i]=x[i]-u[i][j]*x[j];
			}
			x[i]=x[i]/u[i][i];
		}
		return x;	
	}
}