Sto codificando un algoritmo di decomposizione QR in MATLAB, giusto per assicurarsi di avere i meccanismi corretti. Ecco il codice per la funzione principale: function [Q,R] = QRgivens(A)
n = le
sto cercando di calcolare in R una matrice di proiezione P di un arbitrario N x J matrice S: P = S (S'S)^-1 S'
ho state cercando di eseguire questo con la seguente funzione: P <- function(S){
o