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
R语言 丨Gram-Schimth正交化代码
最新推荐文章于 2024-09-14 08:25:16 发布