In deze serie voorbeelden hebben we te maken met matrices waarvan de eigenwaarden reële waarden blijken te zijn. Oftewel, de eigenwaarden en de eigenvectoren bevinden zich in $\mathbb{R}^n$.
Stel, de volgende matrix is gegeven: \begin{equation*} \mathbf{A}= \begin{pmatrix} 8 & 0 & -5 \\ 9 & 3 & -6 \\ 10 & 0 & -7 \end{pmatrix}. \end{equation*} De opdracht is om de eigenwaarden en de bijbehorende eigenvectoren te berekenen. 1. Karakteristieke vergelijking Stel eerst de karakteristieke vergelijking op en los deze vervolgens op. De gevonden oplossingen zijn de eigenwaarden van matrix $ \mathbf{A} $. Als $ \mathbf{I} $ de eenheidsmatrix (ook wel identiteitsmatrix) is van $ \mathbf{A} $ en $ \lambda $ de te vinden eigenwaarde(n), dan is de karakteristieke vergelijking: \begin{equation*} \det(\mathbf{A}-\lambda \mathbf{I})=0. \end{equation*} Als we die opschrijven in matrixvorm, verkrijgen we \begin{equation} \label{eq:characteristic1} \begin{vmatrix} 8-\lambda & 0 & -5 \\ 9 & 3-\lambda & -6 \\ 10 & 0 & -7-\lambda \end{vmatrix} =0. \end{equation} Kies de makkelijkste rij of kolom om de determinant van de matrix van vergelijking \eqref{eq:characteristic1} te berekenen. Oftewel, hier gaan we langs de middelste kolom naar beneden. Deze kolom bevat namelijk de meeste nullen, wat de lengte van de karakteristieke vergelijking danig zal verkorten: \begin{align*} &-0\begin{vmatrix}9 & -6 \\ 10 & -7-\lambda \end{vmatrix} \\ &\quad + (3-\lambda)\begin{vmatrix}8-\lambda & -5 \\ 10 & -7-\lambda \end{vmatrix} \\ &\qquad- 0\begin{vmatrix}8-\lambda & -5 \\ 9 & -6 \end{vmatrix} = 0, \end{align*} en dus \begin{equation} (3-\lambda)\begin{vmatrix}8-\lambda & -5 \\ 10 & -7-\lambda \end{vmatrix} = 0. \end{equation} Als we dit verder uitwerken, verkrijgen we \begin{align*} (3-\lambda)[(8-\lambda)(-7-\lambda)-(-5)(10)] &= 0, \\ (3-\lambda)[-56-8\lambda+7\lambda+\lambda^2+50] &= 0, \\ (3-\lambda)[\lambda^2-\lambda-6] &= 0, \end{align*} waardoor we een hanteerbaardere vorm van de karakteristieke vergelijking hebben gevonden: \begin{equation} (3-\lambda)[(\lambda+2)(\lambda-3)] = 0, \end{equation} waarvan de oplossingen van $ \lambda $ vervolgens makkelijk af te lezen zijn: \begin{equation} \therefore \lambda = 3 \vee \lambda = -2 \vee \lambda = 3. \end{equation} We kunnen ten slotte nog twee controles uitvoeren. De som van de gevonden waarden moet gelijk zijn aan de spoor, of trace, van $\mathbf{A}$, oftewel, de som van de hoofddiagonaal van $\mathbf{A}$: \begin{equation*} \text{tr }\mathbf{A}=8+3-7=4. \end{equation*} En inderdaad is ook de som van de gevonden waarden van $\lambda$ gelijk aan 4. Ten slotte kunnen we nog controleren of het product van de gevonden waarden gelijk is aan $\det\mathbf{A}$: \begin{align*} \det\mathbf{A} &= -0\begin{vmatrix}9&-6\\10&-7\end{vmatrix} \\ &\qquad +3\begin{vmatrix}8&-5\\10&-7\end{vmatrix} \\ &\qquad\quad -0\begin{vmatrix}8&-5\\9&-6\end{vmatrix} \\ &= 3\begin{vmatrix}8&-5\\10&-7\end{vmatrix} \\ &= 3((8)(-7)-(-5)(10)) \\ &= -18. \end{align*} En inderdaad is ook het product van de gevonden waarden van $\lambda$ gelijk aan $3\cdot-2\cdot3=-18$. 2. Benoem de eigenwaarden De eigenwaarden van matrix $ \mathbf{A} $ zijn dus $ \lambda = -2 $ en $ \lambda = 3$. 3. Eigenvectorvergelijkingen We herschrijven nu de karakteristieke vergelijking in matrixvorm naar een stelsel van drie lineaire vergelijkingen. Omdat het de bedoeling is dat we een (of meer) eigenvector(en) $ \mathbf{v} $ gaan vinden, stel dan dat \begin{equation} \label{eq:v01} \mathbf{v} = \begin{pmatrix} x_1 \\ x_2 \\ x_3 \end{pmatrix} \end{equation} en dat \begin{equation} (\mathbf{A}-\lambda\mathbf{I})\mathbf{v}=\mathbf{0}. \end{equation} In dat geval schrijven we \begin{equation} \label{eq:A-lambda I times v = 0} \begin{pmatrix} 8-\lambda & 0 & -5 \\ 9 & 3-\lambda & -6 \\ 10 & 0 & -7-\lambda \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \\ x_3 \end{pmatrix} = \mathbf{0}, \end{equation} wat we vervolgens als een stelsel van drie lineaire vergelijkingen kunnen opstellen: \begin{equation*} \left\{ \begin{matrix}[r] (8-\lambda)x_1 + 0x_2 – 5x_3 = 0, \\ 9x_1 + (3-\lambda)x_2 – 6x_3 = 0, \\ 10x_1 + 0x_2 + (-7-\lambda)x_3 = 0. \end{matrix} \right. \end{equation*} Als we dit nog enigszins vereenvoudigen, hebben we de volgende eigenvectorvergelijkingen verkregen: \begin{equation} \label{eq:eigenvectorvergelijkingen01} \left\{ \begin{matrix}[r] (8-\lambda)x_1 – 5x_3 &= 0, \\ 9x_1 + (3-\lambda)x_2 – 6x_3 &= 0, \\ 10x_1 + (-7-\lambda)x_3 &= 0. \end{matrix} \right. \end{equation} 4. Substitueer elke verkregen eigenwaarde $\boldsymbol{\lambda}$ in de eigenvectorvergelijkingen 4.1. Eigenwaarde $\boldsymbol{\lambda = -2}$ Laten we beginnen met eigenwaarde $ \lambda = -2 $. Als we dit in eigenvectorvergelijkingen \eqref{eq:eigenvectorvergelijkingen01} substitueren, verkrijgen we \begin{align*} (8-(-2))x_1 – 5x_3 &= 0, \\ 9x_1 + (3-(-2))x_2 – 6x_3 &= 0, \\ 10x_1 + (-7-(-2))x_3 &= 0. \end{align*} Dat kunnen we vereenvoudigen tot \begin{align*} 10x_1 &= 5x_3, \\ 9x_1 + 5x_2 – 6x_3 &= 0, \\ 10x_1 &= 5x_3. \end{align*} Van de eerste en de derde vergelijking blijft over: $ 2x_1 = x_3 $. Laten we $x_3$ in de tweede vergelijking substitueren. \begin{align*} 9x_1 + 5x_2 -6(2x_1) &= 0 \\ \therefore 3x_1 &= 5x_2. \end{align*} Met andere woorden, als $ x_1 = 5 $, dan is $ x_2 = 3 $ en $ x_3 = 10 $. We kunnen nu dus de waarden van $ \mathbf{v} $ in \eqref{eq:v01} invullen: \begin{equation} \mathbf{v} = \begin{pmatrix} x_1 \\ x_2 \\ x_3 \end{pmatrix} = \begin{pmatrix} 5 \\ 3 \\ 10 \end{pmatrix}. \end{equation} Met andere woorden, een eigenvector bij de eigenwaarde $ \lambda = -2 $ is $ \begin{pmatrix}5 & 3 & 10\end{pmatrix}^T $. (Noot: we schrijven hier met opzet een eigenvector, aangezien bijvoorbeeld de eigenvector $ \begin{pmatrix}30 & 18 & 60\end{pmatrix}^T $ ook een eigenvector is bij deze eigenwaarde. Zolang $ 2x_1 = x_3 $ en $ 3x_1 = 5x_2 $, met andere woorden, zolang de verhoudingen tussen $ x_1$, $ x_2 $ en $ x_3 $ maar onveranderd blijven, is het een eigenvector die bij deze eigenwaarde hoort. De conventie is echter om de laagst mogelijke gehele waarden te hanteren.) Dit kunnen we nog controleren door het inwendig product van $\mathbf{A}$ met de eigenvector te berekenen. Als het goed is, verkrijgen we het inwendige product met dezelfde eigenvector als resultaat, oftewel, $ \mathbf{Av}=\lambda \mathbf{v} $. Als we dit uitschrijven in matrixnotatie (waarbij we voor de duidelijkheid aangeven welk deel welke variabele is), verkrijgen we inderdaad \begin{align*} & \underbrace{\begin{pmatrix} 8 & 0 & -5 \\ 9 & 3 & -6 \\ 10 & 0 & -7 \end{pmatrix}}_{\mathbf{A}} \underbrace{\begin{pmatrix} 5 \\ 3 \\ 10 \end{pmatrix}}_{\mathbf{v}} \\ &= \begin{pmatrix} 8\cdot5 + 0\cdot3 + -5\cdot10 \\ 9\cdot5 + 3\cdot3 + -6\cdot10 \\ 10\cdot5 + 0\cdot3 + -7\cdot10 \\ \end{pmatrix} \\ &= \begin{pmatrix} -10 \\ -6 \\ -20 \end{pmatrix} = \underbrace{-2}_{\lambda} \underbrace{ \begin{pmatrix} 5 \\ 3 \\ 10 \end{pmatrix}}_{\mathbf{v}}. \end{align*} 4.2. Eigenwaarde $ \boldsymbol{\lambda = 3} $ Op dezelfde manier vervangen we $ \lambda = 3 $ in de eigenvectorvergelijkingen van \eqref{eq:eigenvectorvergelijkingen01}. We schrijven dan de volgende drie lineaire vergelijkingen op: \begin{align*} (8-3)x_1 – 5x_3 &= 0, \\ 9x_1 + (3-3)x_2 – 6x_3 &= 0, \\ 10x_1 + (-7-3)x_3 &= 0. \end{align*} Die kunnen we verder uitwerken: \begin{align*} 5x_1 &= 5x_3, \\ 9x_1 &= 6x_3, \\ 10x_1 &= 10x_3. \end{align*} We zien dat alleen $ x_1=x_3=0 $ een oplossing kan vormen voor dit stelsel van gelijktijdige vergelijkingen. En omdat $ x_2 $ altijd wordt vermenigvuldigd met 0 (die komt dan ook niet langer voor in dit stelsel), kan $ x_2 $ iedere waarde aannemen. Een eigenvector die bij eigenwaarde $ \lambda = 3 $ hoort, is dan ook simpelweg $ \begin{pmatrix}0 & 1 & 0\end{pmatrix}^T $. Als we dit checken met behulp van het inwendig product van matrix $\mathbf{A}$ met $ \begin{pmatrix}0 & 1 & 0\end{pmatrix}^T $, verkrijgen we inderdaad $\lambda\mathbf{v}$: \begin{align*} &\begin{pmatrix} 8 & 0 & -5 \\ 9 & 3 & -6 \\ 10 & 0 & -7 \end{pmatrix} \begin{pmatrix} 0 \\ 1 \\ 0 \end{pmatrix} \\ &= \begin{pmatrix} 8\cdot0 + 0\cdot1 + -5\cdot0 \\ 9\cdot0 + 3\cdot1 + -6\cdot0 \\ 10\cdot0 + 0\cdot1 + -7\cdot0 \\ \end{pmatrix} \\ &= \begin{pmatrix} 0 \\ 3 \\ 0 \end{pmatrix} = 3 \begin{pmatrix} 0 \\ 1 \\ 0 \end{pmatrix}, \end{align*} dus dat klopt.
