Gram Schmidt z R
Oto kod MATLAB do wykonywania Gram Schmidt na stronie 1http://web.mit.edu/18.06/www/Essays/gramschmidtmat.pdf
Próbuję godzinami wykonać to z R, ponieważ nie mam MATLABA Oto moja R
f=function(x){
m=nrow(x);
n=ncol(x);
Q=matrix(0,m,n);
R=matrix(0,n,n);
for(j in 1:n){
v=x[,j,drop=FALSE];
for(i in 1:j-1){
R[i,j]=t(Q[,i,drop=FALSE])%*%x[,j,drop=FALSE];
v=v-R[i,j]%*%Q[,i,drop=FALSE]
}
R[j,j]=max(svd(v)$d);
Q[,j,,drop=FALSE]=v/R[j,j]}
return(list(Q,R))}
Ciągle mówi, że w jednym z nich występują błędy:
v=v-R[i,j]%*%Q[,i,drop=FALSE]
lub
R[j,j]=max(svd(v)$d);
Co robię źle, tłumacząc kod MATLAB na R ???