R/Rakhaqr.R

Defines functions MatriksQ MatriksR

MatriksQ=function(a){
  library(pracma)#hasil kali dalam vector
  a1=matrix(x[,1],3,1)
  v1=a1
  normv1=norm(v1,"F")
  e1=v1/normv1
  a2=matrix(x[,2],3,1)
  v2=a2-(as.vector(t(a2)%*%e1))*e1
  normv2=norm(v2,"F")
  e2=v2/normv2
  a3=matrix(x[,3],3,1)
  v3=v3=a3-(as.vector(t(a3)%*%e1))*e1-(as.vector(t(a3)%*%e2))*e2
  normv3=norm(v3,"F")
  e3=v3/normv3
  a=cbind(e1,e2,e3)
}
MatriksR=function(r){
  library(pracma)#hasil kali dalam vector
  a1=matrix(x[,1],3,1)
  v1=a1
  normv1=norm(v1,"F")
  e1=v1/normv1
  a2=matrix(x[,2],3,1)
  v2=a2-(as.vector(t(a2)%*%e1))*e1
  normv2=norm(v2,"F")
  e2=v2/normv2
  a3=matrix(x[,3],3,1)
  v3=v3=a3-(as.vector(t(a3)%*%e1))*e1-(as.vector(t(a3)%*%e2))*e2
  normv3=norm(v3,"F")
  e3=v3/normv3
  r=matrix(c(dot(a1,e1),dot(a2,e1),dot(a3,e1),0,dot(a2,e2),dot(a3,e2),0,0,dot(a3,e3)),3,3,byrow=T)
}
}
rakha150/Rakhaqr documentation built on Nov. 20, 2019, 12:05 a.m.