R语言 丨Gram-Schimth正交化代码

d<-1e-10
a<-matrix(c(1,d,0,0,1,0,d,0,1,0,0,d),ncol=3)
v<-matrix(c(1,d,0,0,1,0,d,0,1,0,0,d),ncol=3)
r<-matrix(c(0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0),ncol=4)
q<-matrix(c(1,d,0,0,0,0,0,0,0,0,0,0,0,0,0,0),ncol=4)
#is.matrix(q[,c(1)])
#as.matrix(q[,c(1)])
#t(r[1,2]%*%q[,1])
#q[,c(1)]
#as.matrix(v[,2])
n=3
#r[1,2]#说明R矩阵的下标是从1开始计数的
#a<-matrix(c(1,0,0,0,0,1,0,0,0,0,1,0),ncol=3)
for (j in 1:n){
  v[,j]<-as.matrix(a[,j])#结果是列向量！
  for (i in 1:j-1) {
    r[i,j]<-t(q[,i])%*%v[,j]
    #print(r[i,j])
    #print(q[,i])
    v[,j]<-v[,j]-t(r[i,j]%*%t(as.matrix(q[,i])))#报错显示长度不对！
  }
  #print(sqrt(t(v[,j])%*%v[,j]))
  r[j,j]<-as.numeric(sqrt(t(v[,j])%*%v[,j]))
  #print(r[j,j])
  q[,j]<-v[,j]/r[j,j]
  print(q[,j])
}
q
r