SturdyCobble's Study Note

(LDU Decomposition만을 진행하며, 예외 처리 등은 없는 코드입니다. 직접 작성해 오류나 비효율적인 부분이 존재할 수 있습니다.) function [L,D,U] = decomLDU(A) % Decompose the given matrix A as LDU form i = 1; j = 1; [R, C] = size(A); % square matrix A if (R == C) N = R; % since A is N by N. Maybe redundant. L = eye(N); D = eye(N); while i

MATLAB에는 rref라는 함수로 위 Reduced Row Echelon Form을 구할 수 있습니다. 이번 글에서는 개략적으로 이 함수의 동작원리를 살펴보겠습니다. 자세한 사항은 https://en.wikipedia.org/wiki/Gaussian_elimination#Pseudocode 의 코드를 참조하시면 편합니다. (rref함수의 코드는 type rref로 구할 수 있습니다.) function A = getRef(A) % 1 [R,C] = size(A); i = 1; j = 1; while i getRref(A) A = 15 16 17 18 19 7 8 9 10 11 11 12 13 14 15 1 2 3 4 5 A = 1.0000 1.0667 1.1333 1.2000 1.2667 0 0.5333..