java 对矩阵的QR分解

package com.yjb;
/**
 * 把正定矩阵转化为上三角矩阵和一个正交矩阵，如：A=QR
 * */
public class ToDiagMatrix {
      private double[][] R;//上三角矩阵
      private double[][] Q;//正交矩阵
      private double[][] A;//要分解的矩阵A=QR
      private int n;
      public ToDiagMatrix（double[][] A）{
    	  this.A=A;
    	  n=A.length;
          this.R=new double[n][n];
          this.Q=new double[n][n];
          this.initR（）;
          schmidt（）;
      }
      public void initR（）{
    	  for（int i=0;i    		  for（int j=0;j    			  R[i][j]=0.0;
    		  }
    	  }
      }
      public void schmidt（）{
    	  for（int i=0;i    		  R[i][i]=this.Uii（i）;
    		  this.setQi（i）;
    		  for（int j=1;j    			  R[i][j]=this.Uij（i j）;
    		  }
    		  this.updateA（i）;
    	  }