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 ...
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