Tôi đang mã hóa thuật toán phân tách QR trong MATLAB, chỉ để đảm bảo rằng tôi có cơ chế chính xác. Dưới đây là mã cho các chức năng chính: function [Q,R] = QRgivens(A)
n = length(A(:,1));
Tôi đang cố gắng để tính toán trong R một ma trận chiếu P của một ma trận J N x tùy ý S: P = S (S'S)^-1 S'
Tôi đã đã cố gắng để thực hiện điều này với các chức năng sau: P <- function(S){
outpu