// fattorizzazione LU
function [L,U] = LU1 (A)
  [m,n] = size (A)
  if n ~= m
    error ( 'La matrice non è quadrata')
  end
  L = eye ( n,n)
  for k = 1 : ( n-1)
    for i = k +1 : n
      L (i,k) = A (i,k)/ A(k,k)
      for j = k +1  : n
        A(i,j) = A(i,j) - L(i,k)*A(k,j)
      end
    end
  end
  U = triu (A)
endfunction

Citazione:

---

Messaggio inserito da fozzy04

Citazione:

---

Messaggio inserito da miss smile

Salve ragazzi, come da oggetto vorrei sapere se qualcuno è riuscito ad implementare in Scilab l'algoritmo di eliminazione di Gauss "normale" (cioè senza pivot e senza tecnica nel massimo pivot parziale).
Questo è il codice che ho scritto io e non riesco a capire perchè non funzioni!!!

function U=gauss(A,b)
  [m,n]=size(A)
  if m~=n
    error('La matrice dei coefficienti A deve essere quadrata!')
  elseif n~=length(b)
    error('Il vettore dei termini noti deve avere lo stesso numero di righe della matrice A!')
  end
  A1=[A b]
  M=zeros(n,n)
  for k=1:n
    for i=k+1:n
      M(i,k)=-A1(i,k)/A(k,k)
      for j=k+1:n
        A1(i,j)=A(i,j)+(M(i,k)*A(k,j))
      end
    end
  end
  U=A1
  endfunction

Ringrazio in anticipo chi mi fosse d'aiuto.

---

for j=1:n

---

Citazione:

---

Messaggio inserito da fozzy04

A proposito... c'è qualcuno che mi sa dire come si fa a realizzare una libreria in Scilab, così la smetto ogni volta di caricarmi a mano tutte le function che mi servono??? Ho provato ad usare la "genlib" ma credo che voglia il codice compilato...

---

function [x]=linearNN(A,b)
  // RISOLUZIONE DI UN SISTEMA LINEARE QUADRATO NxN
  // MEDIANTE FATTORIZZAZIONE LU CLASSICA
  //
  //  INPUT : A -> MATRICE QUADRATA DEI COEFFICIENTI
  //          b -> VETTORE COLONNA DEI TERMINI NOTI
  //  OUTPUT: x -> VETTORE COLONNA DELLE SOLUZIONI DEL SISTEMA
  //
  //                  STRATEGIA RISOLUTIVA
  //  Per risolvere il sistema Ax=b si eseguono i seguenti passi:
  //  1 - si fattorizza A come prodotto L*U quindi avremo L*U*x=b
  //  2 - si pone U*x=y  
  //  3 - si pone L*y=b
  //  4 - nell'equazione al punto 3 si ricava il vettore y con 
  //      il metodo di sostituzione in avanti
  //  5 - nell'equazione al punto 2 si ricava il vettore x con 
  //      il metodo di sostituzione all'indietro
  
 
  // RECUPERO LA DIMENSIONE DELLA MATRICE DEI COEFFICIENTI
  [m,n]= size(A);
  
  // RECUPERO LA LUNGHEZZA DEL VETTORE DEI TERMINI NOTI
  [p,q]= size(b);
  
  // VERIFICO LE DIMENSIONI DELL'INPUT
  if (m<>n)  |  (p<>m)  |  (q>1)
    error('ERRORE: DIMENSIONI NON CONGRUE!')
  else
    printf('\n *******************************************')
    printf('\n *  RISOLUZIONE DI UN SISTEMA LINEARE NxN  *')
    printf('\n *  MEDIANTE FATTORIZZAZIONE L*U CLASSICA  *')
    printf('\n ******************************************* \n')
  end

  // FATTORIZZAZIONE LU CLASSICA (SENZA PIVOTING)
  [L,U] = fattorizzaLU(A)
  
  // RISOLUZIONE DEL SISTEMA TRIANGOLARE INFERIORE
  y = trianginf(L,b)
  
  // RISOLUZIONE DEL SISTEMA TRIANGOLARE SUPERIORE
  x = triangsup(U,y)

endfunction

Citazione:

---

Messaggio inserito da miss smile

Salve ragazzi, come da oggetto vorrei sapere se qualcuno è riuscito ad implementare in Scilab l'algoritmo di eliminazione di Gauss "normale" (cioè senza pivot e senza tecnica nel massimo pivot parziale).
Questo è il codice che ho scritto io e non riesco a capire perchè non funzioni!!!

function U=gauss(A,b)
  [m,n]=size(A)
  if m~=n
    error('La matrice dei coefficienti A deve essere quadrata!')
  elseif n~=length(b)
    error('Il vettore dei termini noti deve avere lo stesso numero di righe della matrice A!')
  end
  A1=[A b]
  M=zeros(n,n)
  for k=1:n
    for i=k+1:n
      M(i,k)=-A1(i,k)/A(k,k)
      for j=k+1:n
        A1(i,j)=A(i,j)+(M(i,k)*A(k,j))
      end
    end
  end
  U=A1
  endfunction

Ringrazio in anticipo chi mi fosse d'aiuto.

---

Citazione:

---

Messaggio inserito da miss smile

non va cmq!

---

Citazione:

---

Messaggio inserito da miss smile

A1 = [A b] funziona, ho provato. Come seconda cosa, U la devo ottenere fattorizzando A, non semplicemente facendone la sua triangolare superiore...o no?

---