18.01.2015 Views

Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

7.3. Metodo <strong>di</strong> eliminazione <strong>di</strong> Gauss<br />

coefficiente (cioè per − 3 2 ) e poi sottraiamo la terza equazione dalla seconda moltiplicata per − 3 2 :<br />

3<br />

2 x 2 + 4x 3 = 15 −<br />

− 3 2 (−x 2 − 2x 3 = −8) ⇐⇒ 3 2 x 2 + 3x 3 = 12 =<br />

x 3 = 3<br />

– Sostituiamo questa equazione alla terza del sistema, ricavando il sistema equivalente<br />

2x 1 + x 2 + 2x 3 = 10<br />

−x 2 − 2x 3 = −8<br />

x 3 = 3<br />

Con tutte le trasformazioni effettuate, abbiamo trasformato il sistema <strong>di</strong> partenza in uno equivalente, che<br />

si può risolvere facilmente me<strong>di</strong>ante sostituzioni all’in<strong>di</strong>etro. Dall’ultima equazione abbiamo x 3 = 3. Sostituendo<br />

questo valore nella seconda equazione otteniamo −x 2 − 6 = −8 da cui x 2 = 2. Infine, sostituendo i<br />

valori <strong>di</strong> x 3 e x 2 nella prima equazione abbiamo 2x 1 + 2 + 6 = 10 da cui x 1 = 1.<br />

7.3.3 Eliminazione <strong>di</strong> Gauss: caso generale<br />

Sia dato un sistema <strong>di</strong> n equazioni, in cui la matrice A è piena (o densa, cioè abbia quasi tutti gli elementi<br />

non nulli). Applichiamo trasformazioni elementari per riga in modo da ridurre il sistema <strong>di</strong> partenza in uno<br />

equivalente <strong>di</strong> forma triangolare superiore, che potremo risolvere me<strong>di</strong>ante sostituzione all’in<strong>di</strong>etro.<br />

La soluzione del problema Ax = b, infatti, non cambia se moltiplichiamo una riga per una costante, se<br />

sottraiamo il multiplo <strong>di</strong> una riga da un’altra riga o se facciamo scambi <strong>di</strong> righe, come abbiamo detto all’inizio<br />

della Sezione 7.2.<br />

Supponiamo, per il momento, che tutti gli elementi della <strong>di</strong>agonale principale <strong>di</strong> A siano non nulli.<br />

G Al primo passo vogliamo eliminare gli elementi della prima colonna al <strong>di</strong> sotto della <strong>di</strong>agonale<br />

principale:<br />

– sottraiamo la prima equazione moltiplicata per a 21<br />

a 11<br />

dalla seconda equazione:<br />

a 21 x 1 + a 22 x 2 + a 23 x 3 + ... + a 2n x n = b 2 −<br />

a 21<br />

(a 11 x 1 + a 12 x 2 + a 13 x 3 + ... + a 1n x n ) = a 21<br />

b 1 =<br />

a 11 a 11<br />

(a 22 − a 21<br />

a 11<br />

a 12 )x 2 + (a 23 − a 21<br />

a 11<br />

a 13 )x 3 + ... + (a 2n − a 21<br />

a 11<br />

a 1n )x n = b 2 − a 21<br />

a 11<br />

b 1<br />

– sottraiamo la prima equazione moltiplicata per a 31<br />

a 11<br />

dalla terza equazione.<br />

– ...<br />

– sottraiamo la prima equazione moltiplicata per a n1<br />

a 11<br />

dalla n-sima equazione.<br />

Come risultato <strong>di</strong> questa operazione avremo una nuova matrice con gli elementi della prima<br />

colonna, eccetto quello <strong>di</strong> posto 11, tutti uguali a zero.<br />

⎛<br />

⎞⎛<br />

⎞ ⎛ ⎞<br />

a 11 a 12 ... a 1n x 1 b 1<br />

0 a (1)<br />

22<br />

... a (1)<br />

2n<br />

x 2<br />

⎜<br />

⎝<br />

.<br />

. ...<br />

⎟⎜<br />

. ⎠⎝<br />

⎟<br />

. ⎠ = b (1)<br />

2 ⎜ ...<br />

⎟<br />

⎝ ⎠<br />

0 a (1)<br />

n2<br />

... a nn<br />

(1) x n b n<br />

(1)<br />

101

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!