De vierkantswortel van een matrix
Voor een reëel getal
weten we wat een vierkantswortel is: dat is een getal
waarvoor
![Rendered by QuickLaTeX.com \[ b^2=a. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-26653e0cc76ef17b17308958921e62ff_l3.png)
Bij matrices stellen we precies dezelfde vraag.
Definitie
Een vierkantswortel van een matrix
is een matrix
waarvoor
![Rendered by QuickLaTeX.com \[ B^2=A. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-5ab0592dc8ecdd616c57a366f9dd26ae_l3.png)
Met andere woorden: als we

met zichzelf vermenigvuldigen, krijgen we

.
Hierbij moeten we meteen opletten dat de situatie bij matrices subtieler is dan bij gewone getallen:
- een matrix kan geen vierkantswortel hebben;
- een matrix kan meerdere vierkantswortels hebben;
- een vierkantswortel is in het algemeen niet uniek.
Daarom is het vaak beter te spreken van een vierkantswortel van een matrix dan van de vierkantswortel.
Een eerste eenvoudig voorbeeld
Neem de diagonale matrix
![Rendered by QuickLaTeX.com \[ A=\begin{pmatrix} 4&0\\ 0&9 \end{pmatrix}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-d52c995da738f4e79f6350acdcd1355d_l3.png)
Dan is
![Rendered by QuickLaTeX.com \[ B=\begin{pmatrix} 2&0\\ 0&3 \end{pmatrix} \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-28a13f28f22c8b22e54ce493a22f62fc_l3.png)
een vierkantswortel van

, want
![Rendered by QuickLaTeX.com \[ B^2= \begin{pmatrix} 2&0\\ 0&3 \end{pmatrix} \begin{pmatrix} 2&0\\ 0&3 \end{pmatrix} = \begin{pmatrix} 4&0\\ 0&9 \end{pmatrix} =A. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-b4ae9ae570c2425c7057aa252d40c598_l3.png)
Maar ook
![Rendered by QuickLaTeX.com \[ \begin{pmatrix} -2&0\\ 0&3 \end{pmatrix}, \qquad \begin{pmatrix} 2&0\\ 0&-3 \end{pmatrix}, \qquad \begin{pmatrix} -2&0\\ 0&-3 \end{pmatrix} \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-7ec79701b57692a36b1cca213d92af9f_l3.png)
zijn vierkantswortels van

.
Zelfs in dit eenvoudige geval is de vierkantswortel dus niet uniek.
De standaardmethode: diagonaliseerbare matrices
De berekening van een vierkantswortel wordt veel eenvoudiger wanneer de matrix diagonaliseerbaar is.
Stel dat
![Rendered by QuickLaTeX.com \[ A=PDP^{-1}, \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-a091bad1dba1571df5af55b22ee724a8_l3.png)
waarbij

een diagonale matrix is:
![Rendered by QuickLaTeX.com \[ D=\operatorname{diag}(\lambda_1,\lambda_2,\dots,\lambda_n). \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-c79a452e65b2b9a4d579db7155d5b19a_l3.png)
Als we voor elke diagonaalwaarde
een getal
kunnen kiezen zodat
![Rendered by QuickLaTeX.com \[ \mu_i^2=\lambda_i, \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-b92421ea05fd607b8124e5eb27403459_l3.png)
dan is
![Rendered by QuickLaTeX.com \[ \sqrt{D}=\operatorname{diag}(\mu_1,\mu_2,\dots,\mu_n) \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-c670038d7d78c1cb82283749884dda4f_l3.png)
een vierkantswortel van

.
Daaruit volgt dan onmiddellijk een vierkantswortel van
:
![Rendered by QuickLaTeX.com \[ B=P\sqrt{D}P^{-1}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-f8eac523dd03ff3e14f942d3c944cfb5_l3.png)
Inderdaad,
![Rendered by QuickLaTeX.com \[ B^2=(P\sqrt{D}P^{-1})(P\sqrt{D}P^{-1}) = P(\sqrt{D})^2P^{-1} = PDP^{-1} = A. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-dccbaff52d50f5df7283149b26bff80e_l3.png)
De strategie is dus:
- bepaal de eigenwaarden van
;
- diagonaliseer
;
- neem de vierkantswortel van de diagonale matrix;
- keer terug naar de oorspronkelijke basis.
Uitgewerkt voorbeeld
We bekijken de matrix
![Rendered by QuickLaTeX.com \[ A= \begin{pmatrix} 5&4\\ 4&5 \end{pmatrix}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-c202ecdd5208d2ab7533558a427e8424_l3.png)
We zoeken een matrix

waarvoor
![Rendered by QuickLaTeX.com \[ B^2=A. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-5ab0592dc8ecdd616c57a366f9dd26ae_l3.png)
Stap 1: eigenwaarden bepalen
Voor een matrix van de vorm
![Rendered by QuickLaTeX.com \[ \begin{pmatrix} a&b\\ b&a \end{pmatrix} \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-35d76669c93e3a23fa5f9f13d8b197c6_l3.png)
zijn de eigenwaarden
![Rendered by QuickLaTeX.com \[ a+b \quad \text{en} \quad a-b. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-eec66b8f8d39e1dd93ece2ee476794b2_l3.png)
Hier krijgen we dus
![Rendered by QuickLaTeX.com \[ \lambda_1=5+4=9, \qquad \lambda_2=5-4=1. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-7dc5471e8879f9f2a023f3975495a4b6_l3.png)
Stap 2: eigenvectoren bepalen
Bij
hoort bijvoorbeeld de eigenvector
![Rendered by QuickLaTeX.com \[ v_1=\begin{pmatrix}1\\1\end{pmatrix}, \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-d81527bc12bd8f9da3a17d1005e2c67c_l3.png)
en bij

hoort
![Rendered by QuickLaTeX.com \[ v_2=\begin{pmatrix}1\\-1\end{pmatrix}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-16f0eb374dff4a7b58b2ad55034ec5ec_l3.png)
We nemen dus
![Rendered by QuickLaTeX.com \[ P= \begin{pmatrix} 1&1\\ 1&-1 \end{pmatrix}, \qquad D= \begin{pmatrix} 9&0\\ 0&1 \end{pmatrix}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-e28c2d15b24d9b83440bf672a28c36cf_l3.png)
Dan geldt
![Rendered by QuickLaTeX.com \[ A=PDP^{-1}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-bdb4c1f0dfd9cce26387f10949b293f6_l3.png)
Stap 3: de vierkantswortel van 
Een natuurlijke keuze is
![Rendered by QuickLaTeX.com \[ \sqrt{D}= \begin{pmatrix} 3&0\\ 0&1 \end{pmatrix}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-9858afc0cd174b163ca3d0449ae0944a_l3.png)
Stap 4: teruggaan naar de oorspronkelijke basis
We berekenen
![Rendered by QuickLaTeX.com \[ B=P\sqrt{D}P^{-1}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-f8eac523dd03ff3e14f942d3c944cfb5_l3.png)
Omdat
![Rendered by QuickLaTeX.com \[ P^{-1}=\frac12 \begin{pmatrix} 1&1\\ 1&-1 \end{pmatrix}, \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-4ca4d8ef292cde5fbd0da74c7583385d_l3.png)
volgt
![Rendered by QuickLaTeX.com \[ B= \begin{pmatrix} 1&1\\ 1&-1 \end{pmatrix} \begin{pmatrix} 3&0\\ 0&1 \end{pmatrix} \frac12 \begin{pmatrix} 1&1\\ 1&-1 \end{pmatrix}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-673d26724b4297e74bded889cb7041b5_l3.png)
Eerst vinden we
![Rendered by QuickLaTeX.com \[ P\sqrt{D}= \begin{pmatrix} 3&1\\ 3&-1 \end{pmatrix}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-6f31f9a11d366c654b42a6ba45b08aad_l3.png)
Dus
![Rendered by QuickLaTeX.com \[ B= \frac12 \begin{pmatrix} 3&1\\ 3&-1 \end{pmatrix} \begin{pmatrix} 1&1\\ 1&-1 \end{pmatrix} = \frac12 \begin{pmatrix} 4&2\\ 2&4 \end{pmatrix} = \begin{pmatrix} 2&1\\ 1&2 \end{pmatrix}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-4c42675753c41e97cfa7a0dd4ed77e9f_l3.png)
Stap 5: controle
Nu controleren we:
![Rendered by QuickLaTeX.com \[ \begin{pmatrix} 2&1\\ 1&2 \end{pmatrix}^2 = \begin{pmatrix} 2&1\\ 1&2 \end{pmatrix} \begin{pmatrix} 2&1\\ 1&2 \end{pmatrix} = \begin{pmatrix} 5&4\\ 4&5 \end{pmatrix}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-c95776c2856fc79716b83e9ceb05bef6_l3.png)
Dus
![Rendered by QuickLaTeX.com \[ \begin{pmatrix} 2&1\\ 1&2 \end{pmatrix} \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-f39ddfbec31f46f9aaadd378ea8a4553_l3.png)
is inderdaad een vierkantswortel van

.
Niet elke matrix heeft een reële vierkantswortel
Neem bijvoorbeeld
![Rendered by QuickLaTeX.com \[ A= \begin{pmatrix} -1&0\\ 0&1 \end{pmatrix}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-2b92b800cf8b67b557bff23c1dc0729c_l3.png)
Als we een reële diagonale vierkantswortel zouden willen nemen, dan zouden we een reëel getal

nodig hebben met
![Rendered by QuickLaTeX.com \[ x^2=-1, \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-8348412bdc928379ee8237d34e04596a_l3.png)
en dat bestaat niet in

.
Dus niet elke reële matrix heeft een reële vierkantswortel. Vooral negatieve eigenwaarden kunnen problemen veroorzaken.
Het mooie geval: symmetrische positief-definiete matrices
Wanneer
een symmetrische matrix is met strikt positieve eigenwaarden, dan bestaat er precies één symmetrische positief-definiete vierkantswortel. Die noemt men de hoofdvierkantswortel van
.
In dat geval kan men schrijven
![Rendered by QuickLaTeX.com \[ A=QDQ^T, \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-ece2d11337a9e533bec74baff7505a5c_l3.png)
waarbij

orthogonaal is en

diagonaal met positieve diagonaalelementen. Dan is
![Rendered by QuickLaTeX.com \[ \sqrt{A}=Q\sqrt{D}Q^T. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-319ad4a956834952767a9196561beab3_l3.png)
Dit is de netjesste en meest natuurlijke situatie.
Besluit
Een vierkantswortel van een matrix
is een matrix
waarvoor
![Rendered by QuickLaTeX.com \[ B^2=A. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-5ab0592dc8ecdd616c57a366f9dd26ae_l3.png)
Voor diagonaliseerbare matrices is de berekening conceptueel eenvoudig:
![Rendered by QuickLaTeX.com \[ A=PDP^{-1} \quad \Longrightarrow \quad \sqrt{A}=P\sqrt{D}P^{-1}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-2b91f8f3c23dff6262be4fd239a5f8c7_l3.png)
Men reduceert het probleem dus tot het nemen van vierkantswortels van de eigenwaarden.
Toch moeten we opletten:
- de vierkantswortel is in het algemeen niet uniek;
- niet elke matrix heeft een reële vierkantswortel;
- voor symmetrische positief-definiete matrices bestaat er wel een natuurlijke keuze.
Het uitgewerkte voorbeeld
![Rendered by QuickLaTeX.com \[ A= \begin{pmatrix} 5&4\\ 4&5 \end{pmatrix} \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-661fb3ead22f28e8a21c3a90524cf4ed_l3.png)
laat mooi zien hoe eigenwaarden en diagonaliseerbaarheid leiden tot
![Rendered by QuickLaTeX.com \[ \sqrt{A}= \begin{pmatrix} 2&1\\ 1&2 \end{pmatrix}. \]](http://www.wiskundemagie.be/wp-content/ql-cache/quicklatex.com-fb457ab2d591306c4e5ad8d07ac360ae_l3.png)
[/latexpage]