A triangolare inferiore

Se la matrice $A$ è triangolare inferiore si presenta nella forma


\begin{displaymath}
A= \left(
\begin{array}{cccc}
a_{11} & 0 & \cdots & 0 \\ ...
... \\
a_{n1} & a_{n2} & \cdots & a_{nn}
\end{array}
\right)
\end{displaymath}

caratterizzata da $a_{ij}=0 \mbox{ } \forall j>i$. Anche in questo caso il determinante è facilmente calcolabile, infatti si ottiene ancora come il prodotto degli elementi diagonali della matrice, allora

\begin{displaymath}det(A) \neq 0 \quad \Leftrightarrow \quad \prod_{i=1}^n a_{ii...
...ftrightarrow
\quad a_{ii} \neq 0 \mbox{ } \forall i=1,\cdots,n\end{displaymath}

Proviamo con un esempio di una matrice triangolare inferiore $3
\times 3$ a vedere come funziona il modo di risoluzione:


\begin{displaymath}
\begin{array}{cccc}
\par
\end{displaymath}

(
\begin{array}{ccc}
a_{11} & 0 & 0 \\
a_{21} & a_{22} & 0 \\
a_{31} & a_{32} & a_{33}
\end{array}
)

\begin{displaymath}
\par
&
\par
\end{displaymath}

(
\begin{array}{c}
x_1\\
x_2 \\
x_3
\end{array}
)

\begin{displaymath}
\par
&
=
&
\par
\end{displaymath}

(
\begin{array}{c}
b_1\\
b_2 \\
b_3
\end{array}
)

\begin{displaymath}
\par
\end{array}
\end{displaymath}


\begin{displaymath}\Downarrow\end{displaymath}


\begin{displaymath}
\left\{
\begin{array}{ccccccc}
a_{11}x_1 & & & & & = & b_...
...2}x_2 & + & a_{33}x_3 & = & b_3
\end{array}
\right.\mbox{.}
\end{displaymath}

Si vede che la prima componente del vettore soluzione $x$, $x_1$, è già disponibile, per cui possiamo trovare

\begin{displaymath}x_1=\frac{b_1}{a_{11}}\end{displaymath}

ed il fatto che gli elementi diagonali della matrice siano diversi da zero rende l'espressione ben formata.

Noto $x_1$, possiamo adesso ricavare $x_2$ nel seguente modo

\begin{displaymath}x_2=\frac{b_2 -
a_{21}\overbrace{\frac{b_1}{a_{11}}}^{x_1}}{a_{22}}\end{displaymath}

ed adesso siamo in grado di trovare l'ultima componente del vettore soluzione come


\begin{displaymath}x_3=\frac{b_3 - a_{31}x_1 - a_{32}x_2}{a_{33}}\mbox{.}\end{displaymath}

Quello che abbiamo cercato di fare con questo semplice esempio era mostrare come si poteva ottenere il vettore $x$: tramite un algoritmo iterativo si ottengono le componenti del vettore soluzione tramite la relazione

\begin{displaymath}x_i=\frac{b_i - \sum_{j=1}^{i-1} a_{ij}x_j}{a_{ii}}\end{displaymath}

e partendo da $i=1$ si arriva alla soluzione desiderata.



Subsections
Morpheus 2004-01-04