Fattorizzazione PA = LU

Messaggioda davicos » 27/05/2017, 16:13

Salve a tutti,
circa il problema: "Data una matrice $A$ si risolva il sistema lineare $Ax = b$ mediante fattorizzazione $PA = LU$, con $b$ definito in modo tale che la corrispondente soluzione $x$ coincida con il vettore unitario".

Il programma è:
Codice:

A = [11 2 3 4; 0 8 2 3; 0 0 4 0;1 0 0 5];
b = sum(A,2);
[L,U,P] = lu(A);
y = L\(P*b);
x = U\y



Non riesco a capire perchè viene calcolato $b$ in quella maniera se so già che la soluzione $x$ è il vettore unitario.
Non capisco cosa si debba risolvere. Come incognita ho la $b$ che la posso calcolare con $b = A*x$ e mi ritrovo il medesimo valore di $b$ calcolato come da programma.
Alla fine del programma $x$ è un vettore unitario, ma questo lo sapevo già dal testo. In definitiva non capisco il senso di tutto ciò e non mi è chiaro il significato del comando 'sum'.
Se poteste darmi delucidazioni in merito ve ne sarei grato.
Grazie.
davicos
Average Member
Average Member
 
Messaggio: 310 di 862
Iscritto il: 16/09/2015, 11:24

Re: Fattorizzazione PA = LU

Messaggioda Raptorista » 29/05/2017, 08:54

La soluzione ti viene data solo per poter verificare che tu abbia fatto bene XD
Tu devi calcolare \(b = Ax\), poi dimenticarti di aver avuto \(x\) e risolvere il sistema lineare. Quando hai fatto confronti la tua \(x\) trovata con quella esatta e vedi se hai fatto bene o male.
Un matematico ha scritto:... come mia nonna che vuole da anni il sistema per vincere al lotto e crede che io, in quanto matematico, sia fallito perché non glielo trovo


Immagine
Avatar utente
Raptorista
Moderatore
Moderatore
 
Messaggio: 4388 di 9616
Iscritto il: 28/09/2008, 19:58

Re: Fattorizzazione PA = LU

Messaggioda davicos » 29/05/2017, 14:55

Bene.
Grazie.
davicos
Average Member
Average Member
 
Messaggio: 311 di 862
Iscritto il: 16/09/2015, 11:24


Torna a Analisi Numerica e Ricerca Operativa

Chi c’è in linea

Visitano il forum: Nessuno e 1 ospite