From 6b30e9e8cb288fee8d3d989ad31480a0a40f855b Mon Sep 17 00:00:00 2001 From: Giovanni Di Grezia Date: Wed, 15 Oct 2014 15:38:30 +0200 Subject: [PATCH] Aggiunto metodi di fattorizzazione di matrice --- samples/LU.m | 31 +++++++++++++++++++ ...e_complete_liner_equation_system_resolve.m | 2 +- 2 files changed, 32 insertions(+), 1 deletion(-) create mode 100644 samples/LU.m diff --git a/samples/LU.m b/samples/LU.m new file mode 100644 index 0000000..de51582 --- /dev/null +++ b/samples/LU.m @@ -0,0 +1,31 @@ + +%se A è la matrice originale +% e A1 la matrice dopo la riduzione di gauss naif +% allora A = L * U + +L=tril(A1,-1) + eye (size(A1)); +U=triu(A1); + +A = L * U; + + + + + +%se A è la matrice originale +% e A1 la matrice dopo la riduzione di gauss pivot +% allora P A = L * U + +L=tril(A1,-1) + eye (size(A1)); +U=triu(A1); +% P matrice di permutazione in base al pivot + +P * A = L * U; + + + +%comando matlab per fattorizzare A +[L,U,P] = lu(A); + + + diff --git a/samples/sample_complete_liner_equation_system_resolve.m b/samples/sample_complete_liner_equation_system_resolve.m index 0372e5c..734f338 100644 --- a/samples/sample_complete_liner_equation_system_resolve.m +++ b/samples/sample_complete_liner_equation_system_resolve.m @@ -6,6 +6,6 @@ U=[2,4,-2;1,-1,5;4,1,-2] % matrice di input b=[6;-3;-10] %termini noti -[x,y,z] = convert_matrix_to_triangular_matrix_gauss_pivoting(U,b) +[x,y] = convert_matrix_to_triangular_matrix_gauss_naif(U,b) [x1] = linear_system_resolver_triangular_matrix(x,y)