Mathematics · Book 4 · Bachelor Year 2

Universitaire wiskunde — Bachelor jaar 2

Universitaire wiskunde — Bachelor jaar 2 · Bachelor Year 2

3Reductie van endomorfismen

Om een endomorfisme te begrijpen, vind de richtingen die het slechts rekt. Dit hoofdstuk bouwt de machinerie — eigenwaarden, karakteristieke en minimale veeltermen, het kernontbindingslemma — en haar beloningen: diagonaliseerbaarheids- en trigonaliseerbaarheidscriteria, Cayley–Hamilton, de Dunford-ontbinding, en de berekening van machten en exponentialen die Hoofdstuk 16 zal voeden. Overal is EE een eindigdimensionale KK-vectorruimte (K=RK = \R of C\C) en uL(E)u \in \mathcal{L}(E), n=dimEn = \dim E.

3.1 Eigenwaarden en eigenvectoren

Definitie 3.1

λK\lambda \in K is een eigenwaarde van uu wanneer u(x)=λxu(x) = \lambda x voor zekere x0x \neq 0 (een eigenvector); de eigenruimte is Eλ(u)=ker(uλid)E_\lambda(u) = \ker(u - \lambda\,\mathrm{id}). De verzameling van eigenwaarden is het spectrum Sp(u)\operatorname{Sp}(u). Een deelruimte FF is stabiel wanneer u(F)Fu(F) \subseteq F; eigenruimten zijn stabiel, en stabiele deelruimten laten geïnduceerde endomorfismen uFu|_F toe.

Stelling 3.2 (Onafhankelijkheid van eigenruimten)

Eigenvectoren behorend bij paarsgewijs verschillende eigenwaarden vormen een vrije familie; equivalent, is de som van de eigenruimten Eλ1++EλrE_{\lambda_1} + \dots + E_{\lambda_r} (verschillende λi\lambda_i) direct. In het bijzonder heeft uu hoogstens nn eigenwaarden.

Bewijs. Door inductie op rr. Stel x1++xr=0x_1 + \dots + x_r = 0 met xiEλix_i \in E_{\lambda_i}, de uitspraak bekend voor r1r - 1. Pas uu toe en trek λr\lambda_r maal de relatie af:

i=1r1(λiλr)xi=0,\sum_{i=1}^{r-1} (\lambda_i - \lambda_r)\, x_i = 0 ,

dus door inductie elke (λiλr)xi=0(\lambda_i - \lambda_r)x_i = 0, d.w.z. xi=0x_i = 0 voor i<ri < r, dan xr=0x_r = 0. Directe sommen van niet-nulle ruimten in een ruimte van dimensie nn hebben hoogstens nn sommanden.

De matrix A = psmallmatrix2 & 1\\ 1 & 2 psmallmatrix werkend op het vlak: de generieke vector e_1 wordt van haar lijn gestoten, maar de eigenrichtingen v_1 = (1,1) en v_2 = (1,-1) worden slechts gerekt — met 3 en met 1 (dus Av_2 = v_2: het gestippelde beeld valt samen met v_2). Diagonaliseren is de overgang naar de basis (v_1, v_2), waar A diag(3, 1) wordt.
De matrix A=(2112)A = \left(\begin{smallmatrix}2 & 1\\ 1 & 2\end{smallmatrix}\right) werkend op het vlak: de generieke vector e1e_1 wordt van haar lijn gestoten, maar de eigenrichtingen v1=(1,1)v_1 = (1,1) en v2=(1,1)v_2 = (1,-1) worden slechts gerekt — met 33 en met 11 (dus Av2=v2Av_2 = v_2: het gestippelde beeld valt samen met v2v_2). Diagonaliseren is de overgang naar de basis (v1,v2)(v_1, v_2), waar AA diag(3,1)\operatorname{diag}(3, 1) wordt.

Definitie 3.3 (Karakteristieke veelterm)

χu(X)=det(Xidu)\chi_u(X) = \det(X\,\mathrm{id} - u) — berekend in elke basis als det(XInA)\det(XI_n - A), een monische veelterm van graad nn, invariant onder gelijkvormigheid (Stelling 2.17). Haar wortels in KK zijn precies de eigenwaarden (λ\lambda eigenwaarde     uλid\iff u - \lambda\,\mathrm{id} niet injectief     χu(λ)=0\iff \chi_u(\lambda) = 0), en

χu(X)=Xn(tru)Xn1++(1)ndetu.\chi_u(X) = X^n - (\operatorname{tr} u)\, X^{n-1} + \dots + (-1)^n \det u .

De algebraïsche multipliciteit mλm_\lambda van een eigenwaarde is haar multipliciteit als wortel van χu\chi_u; de meetkundige multipliciteit is dimEλ\dim E_\lambda, en 1dimEλmλ1 \leq \dim E_\lambda \leq m_\lambda.

Bewijs van de gestelde feiten. De coëfficiëntclaims: ontwikkel det(XIA)\det(XI - A) door de permutatieformule; de identiteitspermutatie draagt i(Xaii)=Xn(aii)Xn1+\prod_i (X - a_{ii}) = X^n - (\sum a_{ii})X^{n-1} + \dots bij, en elke andere permutatie fixeert hoogstens n2n - 2 diagonaalposities, bijdragend graad n2\leq n - 2: de top twee coëfficiënten zijn zoals gesteld; X=0X = 0 geeft de constante term det(A)=(1)ndetA\det(-A) = (-1)^n\det A.

Meetkundig \leq algebraïsch: zij d=dimEλd = \dim E_\lambda en vul een basis van EλE_\lambda aan tot een basis van EE; de matrix van uu is blokbovendriehoekig met linksboven blok λId\lambda I_d, dus χu(X)=(Xλ)dχ(onderblok)(X)\chi_u(X) = (X - \lambda)^d\, \chi_{\text{(onderblok)}}(X): de multipliciteit van λ\lambda is minstens dd.

Voorbeeld 3.4 (Zelfde χ\chi, verschillende meetkunde)

De matrices

(2002)en(2102)\begin{pmatrix}2 & 0\\ 0 & 2\end{pmatrix} \qquad\text{en}\qquad \begin{pmatrix}2 & 1\\ 0 & 2\end{pmatrix}

delen de karakteristieke veelterm (X2)2(X - 2)^2, het spoor, de determinant, het spectrum — maar zijn niet gelijkvormig: de eerste heeft E2E_2 van dimensie 22 (meetkundige multipliciteit 22), de tweede van dimensie 11. De karakteristieke veelterm ziet alleen algebraïsche multipliciteiten; de eigenruimtedimensies zijn de fijnere invariant, en de minimale veelterm beslecht (X2X - 2 versus (X2)2(X - 2)^2). Moraal voor alle diagonaliseerbaarheidsdiscussies: χ\chi shortlist de kandidaten, maar kernen brengen de stemmen uit.

Definitie 3.5 (Diagonaliseerbaar, trigonaliseerbaar)

uu is diagonaliseerbaar wanneer EE een basis van eigenvectoren heeft (matrix: gelijkvormig met een diagonaalmatrix); trigonaliseerbaar wanneer haar matrix in zekere basis bovendriehoekig is.

Stelling 3.6 (Diagonaliseerbaarheidscriteria)

De volgende zijn equivalent:

  1. uu is diagonaliseerbaar;
  2. E=λSpuEλE = \bigoplus_{\lambda \in \operatorname{Sp} u} E_\lambda;
  3. χu\chi_u splitst over KK en dimEλ=mλ\dim E_\lambda = m_\lambda voor elke eigenwaarde;
  4. (voldoende, niet noodzakelijk) χu\chi_u heeft nn verschillende wortels in KK.

Bewijs. (1     \iff 2): een basis van eigenvectoren sorteert in bases van de EλE_\lambda’s, en omgekeerd geeft concatenatie van bases van de directe sommanden een basis van EE (Stelling 3.2 maakt de som direct; dimensiegelijkheid maakt haar alles).

(2     \iff 3): in de diagonale basis splitst χu=(Xλ)dimEλ\chi_u = \prod (X - \lambda)^{\dim E_\lambda} met matchende multipliciteiten. Omgekeerd, stel χu\chi_u splitst met dimEλ=mλ\dim E_\lambda = m_\lambda overal; dan heeft de directe som van de eigenruimten (direct door Stelling 3.2) dimensie

λdimEλ=λmλ=degχu=n,\sum_{\lambda}\dim E_\lambda = \sum_{\lambda} m_\lambda = \deg\chi_u = n ,

de middelste gelijkheid omdat de graad van een gesplitste veelterm de som van haar wortelmultipliciteiten is: de som is heel EE. Merk waar elke hypothese werkte: splijting vulde de graad, multipliciteitsgelijkheid vulde de dimensies.

(4 \Rightarrow 1): nn verschillende eigenwaarden geven nn onafhankelijke eigenvectoren (Stelling 3.2): een basis.

Methode 3.7 (Diagonaliseerbaarheid beslissen)

In de praktijk, test in deze volgorde — elke stap kan de klus klaren. (1) Presenteert zich een annulerende veelterm met eenvoudige gesplitste wortels (u2=idu^2 = \mathrm{id}, u2=uu^2 = u, uk=idu^k = \mathrm{id})? Zo ja: diagonaliseerbaar, geen berekening (Gevolg 3.17 hieronder). (2) Bereken χu\chi_u; als ze nn verschillende wortels in KK heeft: diagonaliseerbaar (Stelling 3.6 (4)). (3) Anders, voor elke meervoudige wortel λ\lambda alleen, vergelijk dimker(uλid)\dim\ker(u - \lambda\,\mathrm{id}) met de multipliciteit mλm_\lambda: elk tekort doodt diagonaliseerbaarheid; gelijkheid overal bewijst ze. Bereken nooit eigenruimten van enkelvoudige wortels (hun dimensie is gedwongen 11), en trigonaliseer nooit alleen om te beslissen.

Voorbeeld 3.8 (Diagonaliseren aan het werk)

A=I+J=(211121112)A = I + J = \left(\begin{smallmatrix}2 & 1 & 1\\ 1 & 2 & 1\\ 1 & 1 & 2\end{smallmatrix}\right), met JJ de all-eens-matrix: uit Sp(J)={3,0}\operatorname{Sp}(J) = \{3, 0\} (Voorbeeld 2.19), Sp(A)={4,1}\operatorname{Sp}(A) = \{4, 1\}, met eigenruimten R(1,1,1)\R(1,1,1) en het vlak {x+y+z=0}\{x + y + z = 0\}: dimensies 1+2=31 + 2 = 3, diagonaliseerbaar (Stelling 3.6 (2)). Machten zonder enige basisveranderingsmatrix: met Π=J/3\Pi = J/3 de projector op R(1,1,1)\R(1,1,1),

A=4Π+1(IΠ)Ak=4kΠ+(IΠ)=4k13J+I.A = 4\,\Pi + 1\cdot(I - \Pi) \quad\Longrightarrow\quad A^k = 4^k\,\Pi + (I - \Pi) = \frac{4^k - 1}{3}\,J + I .

(Controle k=1k = 1: 413J+I=A\frac{4-1}3 J + I = A.) Het sluitende inzicht: wanneer de eigenruimten zichtbaar zijn, berekenen spectrale projectoren machten sneller dan PDP1PDP^{-1} ooit zal — en de formule toont de dynamica: AkA^k groeit als 4k4^k langs (1,1,1)(1,1,1) en blijft staan op het orthogonale vlak.

Stelling 3.9 (Trigonalisatie)

uu is trigonaliseerbaar over KK desda χu\chi_u over KK splitst. In het bijzonder is elk endomorfisme van een C\C-vectorruimte trigonaliseerbaar.

Bewijs. (\Rightarrow) De karakteristieke veelterm van een driehoekige matrix is (Xtii)\prod(X - t_{ii}): gesplitst.

(\Leftarrow) Inductie op nn. Omdat χu\chi_u splitst, heeft ze een wortel λ\lambda: kies een eigenvector e1e_1. In een basis beginnend met e1e_1 is de matrix (λ0B)\begin{pmatrix} \lambda & \ast\\ 0 & B\end{pmatrix}, en χu=(Xλ)χB\chi_u = (X - \lambda)\chi_B: χB\chi_B splitst ook. Door de inductiehypothese toegepast op de (n1)×(n1)(n-1) \times (n-1)-matrix BB is er een inverteerbare QQ met Q1BQQ^{-1}BQ bovendriehoekig; conjugatie van de hele matrix door (100Q)\begin{pmatrix}1 & 0\\ 0 & Q\end{pmatrix} driehoekigt haar.

Voorbeeld 3.10 (Met de hand trigonaliseren)

B=(3111)B = \begin{pmatrix}3 & -1\\ 1 & 1\end{pmatrix}: χB=X24X+4=(X2)2\chi_B = X^2 - 4X + 4 = (X - 2)^2, en ker(B2I)=ker(1111)\ker(B - 2I) = \ker\left(\begin{smallmatrix}1 & -1\\ 1 & -1\end{smallmatrix}\right) is de lijn opgespannen door e1=(1,1)e_1' = (1, 1): één eigenwaarde, een eendimensionale eigenruimte — niet diagonaliseerbaar, maar trigonaliseerbaar (Stelling 3.9). Vul de basis aan met e2=(1,0)e_2' = (1, 0) en bereken:

u(e1)=(2,2)=2e1,u(e2)=(3,1)=1e1+2e2,u(e_1') = (2, 2) = 2e_1', \qquad u(e_2') = (3, 1) = 1\cdot e_1' + 2\, e_2' ,

dus in de basis (e1,e2)(e_1', e_2') is de matrix T=(2102)T = \left(\begin{smallmatrix}2 & 1\\ 0 & 2\end{smallmatrix}\right). Het sluitende inzicht: de diagonaal van TT was gedwongen (beide entries moeten de dubbele eigenwaarde 22 zijn); alleen de hoekentry hing af van de keuze van e2e_2', en herschalen van e2e_2' kan er elke niet-nulle waarde van maken — de weerstandige “11” is de schaduw van het nilpotente deel dat Dunford zal isoleren.

3.2 Veeltermen van een endomorfisme

Definitie 3.11

Voor P=akXkK[X]P = \sum a_k X^k \in K[X], stel P(u)=akukL(E)P(u) = \sum a_k u^k \in \mathcal{L}(E). De afbeelding PP(u)P \mapsto P(u) is een morfisme van algebras K[X]L(E)K[X] \to \mathcal{L}(E) (Definitie 1.33); haar kern {P:P(u)=0}\{P : P(u) = 0\} is een ideaal van K[X]K[X], niet-nul (de familie (id,u,,un2)(\mathrm{id}, u, \dots, u^{n^2}) is afhankelijk in de n2n^2-dimensionale L(E)\mathcal{L}(E)), dus voortgebracht door een unieke monische veelterm μu\mu_u: de minimale veelterm (Stelling 1.26).

Propositie 3.12

  1. P(u)=0    μuPP(u) = 0 \iff \mu_u \mid P; eigenwaarden van uu zijn wortels van elke annulerende veelterm, en de wortels van μu\mu_u zijn precies de eigenwaarden.
  2. Als FF stabiel is, μuFμu\mu_{u|_F} \mid \mu_u.

Bewijs. (1) De deelbaarheid is de definitie van een voortbrenger. Als u(x)=λxu(x) = \lambda x, x0x \neq 0, dan 0=P(u)(x)=P(λ)x0 = P(u)(x) = P(\lambda)x, dus P(λ)=0P(\lambda) = 0: eigenwaarden zijn wortels van annihilatoren, in het bijzonder van μu\mu_u. Omgekeerd, als λ\lambda een wortel is, μu=(Xλ)Q\mu_u = (X - \lambda)Q met Q(u)0Q(u) \neq 0 (graad van μu\mu_u minimaal): kies yy met Q(u)(y)0Q(u)(y) \neq 0; dan (uλ)(Q(u)(y))=μu(u)(y)=0(u - \lambda)(Q(u)(y)) = \mu_u(u)(y) = 0 vertoont de eigenvector Q(u)(y)Q(u)(y).

(2) μu(uF)=μu(u)F=0\mu_u(u|_F) = \mu_u(u)|_F = 0, en pas (1) toe op uFu|_F.

Voorbeeld 3.13 (Minimale veeltermen met de hand gevonden)

De minimale veelterm wordt berekend door opeenvolgende graden te testen. Voor de all-eens-matrix JM3(R)J \in \mathcal{M}_3(\R): JλIJ \neq \lambda I (graad 11 is uit), en J2=3JJ^2 = 3J, dus

μJ=X23X=X(X3):\mu_J = X^2 - 3X = X(X - 3) :

graad 22, gesplitst, enkelvoudige wortels — JJ is diagonaliseerbaar met spectrum {0,3}\{0, 3\} (Gevolg 3.17 hieronder), bevestigend Voorbeeld 2.19 zonder één enkele determinant. Voor de swapmatrix AA van Voorbeeld 3.15: A±IA \neq \pm I en A2=IA^2 = I geven μA=X21\mu_A = X^2 - 1. In beide gevallen is het patroon hetzelfde: raad een lage-graad-identiteit uit de structuur (rang één dwingt J2=(trJ)JJ^2 = (\operatorname{tr}J)\,J; een involutie dwingt A2=IA^2 = I), controleer dan dat geen echte deler annuleert. Minimale veeltermen worden meestal gevonden, niet berekend uit χ\chi.

Stelling 3.14 (Kernontbindingslemma)

Als P=P1P2PrP = P_1 P_2 \cdots P_r met de PiP_i paarsgewijs relatief priem, dan

kerP(u)=kerP1(u)kerPr(u),\ker P(u) = \ker P_1(u) \oplus \dots \oplus \ker P_r(u),

en de projecties op de sommanden zijn veeltermen in uu.

Bewijs. Het volstaat r=2r = 2 te behandelen en te induceren. Bézout in K[X]K[X] (Stelling 1.26): UP1+VP2=1U P_1 + V P_2 = 1, dus voor elke xx,

x=U(u)P1(u)(x)=:x2+V(u)P2(u)(x)=:x1.x = \underbrace{U(u)P_1(u)(x)}_{=:\,x_2} + \underbrace{V(u)P_2(u)(x)}_{=:\,x_1}.

Als xkerP(u)x \in \ker P(u): P2(u)(x2)=U(u)P(u)(x)=0P_2(u)(x_2) = U(u)\,P(u)(x) = 0 (veeltermen in uu commuten), dus x2kerP2(u)x_2 \in \ker P_2(u), en symmetrisch x1kerP1(u)x_1 \in \ker P_1(u): de som vult kerP(u)\ker P(u); beide sommanden zitten in kerP(u)\ker P(u) (PiPP_i \mid P). Directheid: xkerP1(u)kerP2(u)x \in \ker P_1(u) \cap \ker P_2(u) geeft x=U(u)P1(u)x+V(u)P2(u)x=0x = U(u)P_1(u)x + V(u)P_2(u)x = 0. De formules voor x1,x2x_1, x_2 vertonen de projecties als V(u)P2(u)V(u)P_2(u) en U(u)P1(u)U(u)P_1(u).

Voorbeeld 3.15 (Het kernlemma met expliciete projectoren)

Zij A=(010100001)A = \left(\begin{smallmatrix}0 & 1 & 0\\ 1 & 0 & 0\\ 0 & 0 & 1\end{smallmatrix}\right) (verwissel de eerste twee coördinaten). Dan A2=IA^2 = I: de veelterm X21=(X1)(X+1)X^2 - 1 = (X - 1)(X + 1) annuleert AA, haar factoren zijn relatief priem, en Bézout is expliciet:

12(X+1)12(X1)=1.\frac{1}{2}(X + 1) - \frac12(X - 1) = 1 .

Het bewijs van Stelling 3.14 volgend, zijn de projecties op ker(AI)\ker(A - I) en ker(A+I)\ker(A + I) de veeltermen in AA

π+=A+I2=12(110110002),π=IA2=12(110110000).\pi_+ = \frac{A + I}{2} = \frac12\begin{pmatrix} 1 & 1 & 0\\ 1 & 1 & 0\\ 0 & 0 & 2\end{pmatrix}, \qquad \pi_- = \frac{I - A}{2} = \frac12\begin{pmatrix} 1 & -1 & 0\\ -1 & 1 & 0\\ 0 & 0 & 0\end{pmatrix}.

Controle: π++π=I\pi_+ + \pi_- = I, π+π=0\pi_+\pi_- = 0, π±2=π±\pi_\pm^2 = \pi_\pm, en de beelden zijn het vlak {x=y}\{x = y\} (symmetrische vectoren, eigenwaarde 11) en de lijn R(1,1,0)\R(1, -1, 0) (antisymmetrisch, eigenwaarde 1-1). Het kernlemma is geen bestaanuitspraak: Bézout-coëfficiënten zijn de projectorformules.

Voorbeeld 3.16 (Projectoren berekenen ook de exponentiaal)

Dezelfde swapmatrix, één dividend verder. Omdat A=π+πA = \pi_+ - \pi_- met algebraïsch-orthogonale projectoren (π+π=0\pi_+\pi_- = 0), gehoorzaamt elke macht Ak=π++(1)kπA^k = \pi_+ + (-1)^k\pi_-, en de exponentiaalreeks hergroepeert per projector:

etA=ktkk!(π++(1)kπ)=etπ++etπ=(coshtsinht0sinhtcosht000et).\eu^{tA} = \sum_k \frac{t^k}{k!}\bigl(\pi_+ + (-1)^k\pi_-\bigr) = \eu^{t}\,\pi_+ + \eu^{-t}\,\pi_- = \begin{pmatrix} \cosh t & \sinh t & 0\\ \sinh t & \cosh t & 0\\ 0 & 0 & \eu^{t} \end{pmatrix}.

(Controle t=0t = 0: de identiteit; afgeleide in 00: AA.) De eigenontbinding zet een matrixreeks om in twee scalaire reeksen — precies het mechanisme dat Hoofdstuk 16 op elk diagonaliseerbaar stelsel zal draaien, en de reden dat hyperbolische functies symmetrische koppelingen regeren.

Gevolg 3.17 (Diagonaliseerbaarheid via de minimale veelterm)

uu is diagonaliseerbaar     \iff μu\mu_u splitst over KK met enkelvoudige wortels     \iff zekere annulerende veelterm van uu splitst met enkelvoudige wortels.

Bewijs. Als P(u)=0P(u) = 0 met P=i(Xλi)P = \prod_{i}(X - \lambda_i) (verschillende λi\lambda_i), geeft het lemma E=kerP(u)=iker(uλi)E = \ker P(u) = \bigoplus_i \ker(u - \lambda_i): een directe som van eigenruimten, dus is uu diagonaliseerbaar (Stelling 3.6). Omgekeerd wordt een diagonaliseerbare uu gedood door λSpu(Xλ)\prod_{\lambda \in \operatorname{Sp}u}(X - \lambda) (ze doodt elke eigenruimte), die splitst met enkelvoudige wortels; en μu\mu_u deelt haar terwijl ze dezelfde wortels heeft (Propositie 3.12): μu\mu_u is precies dat product.

Voorbeeld 3.18

Projecties voldoen aan p2=pp^2 = p: geannuleerd door X(X1)X(X-1), gesplitst enkelvoudige wortels — diagonaliseerbaar met spectrum {0,1}\subseteq \{0, 1\}, en E=kerpker(pid)E = \ker p \oplus \ker(p - \mathrm{id}): de meetkundige analyse van Jaar 1, opnieuw bewezen in één regel. Symmetrieën (s2=ids^2 = \mathrm{id}, annihilator X21X^2 - 1): diagonaliseerbaar wanneer charK2\operatorname{char} K \neq 2, spectrum {±1}\subseteq \{\pm 1\}. Een endomorfisme met u3=u2u^3 = u^2 en u2uu^2 \neq u: geannuleerd door X2(X1)X^2(X - 1), niet noodzakelijk diagonaliseerbaar — het criterium detecteert het (dubbele wortel 00 moet worden getest: diagonaliseerbaar desda bovendien keru2=keru\ker u^2 = \ker u).

Voorbeeld 3.19 (Het lichaam beslist: een rotatie in R3\R^3)

Zij RR de kwartslag om de zz-as:

R=(010100001),χR=(X1)(X2+1).R = \begin{pmatrix} 0 & -1 & 0\\ 1 & 0 & 0\\ 0 & 0 & 1 \end{pmatrix}, \qquad \chi_R = (X - 1)(X^2 + 1).

Over R\R: de enige eigenwaarde is 11, met eigenruimte de as Re3\R e_3 — één lijn vaste vectoren, en geen verdere reductie: RR is noch diagonaliseerbaar noch trigonaliseerbaar in M3(R)\mathcal{M}_3(\R) (χR\chi_R splitst niet). Over C\C: drie verschillende eigenwaarden 1,i,i1, \iu, -\iu, dus is RR diagonaliseerbaar, met eigenvectoren e3e_3 en e1ie2e_1 \mp \iu e_2. De meetkunde was hoorbaar in de algebra: rotaties in het vlak hebben geen reële invariante richtingen, en de complexe eigenwaarden ±i\pm\iu van modulus 11 bewaren de hoek (±π2\pm\frac\pi2) die de reële matrix alleen kan uitdrukken door coördinaten te mengen.

Voorbeeld 3.20 (Minimaal versus karakteristiek)

Voor D=diag(2,2,3)D = \operatorname{diag}(2, 2, 3): χD=(X2)2(X3)\chi_D = (X - 2)^2(X - 3) maar μD=(X2)(X3)\mu_D = (X - 2)(X - 3), omdat (D2I)(D3I)=0(D - 2I)(D - 3I) = 0 (controleer op de canonieke basis) terwijl geen van beide factoren alleen DD doodt. Voor het shiftblok N=(0100)(3)N = \left(\begin{smallmatrix}0 & 1\\ 0 & 0\end{smallmatrix}\right) \oplus (3), d.w.z. N=(010000003)N' = \left(\begin{smallmatrix}0 & 1 & 0\\ 0 & 0 & 0\\ 0 & 0 & 3\end{smallmatrix}\right): χN=X2(X3)\chi_{N'} = X^2(X - 3) en μN=X2(X3)\mu_{N'} = X^2(X - 3) — de dubbele wortel is echt nodig omdat NN' niet diagonaliseerbaar is aan de kern-kant (Ne2=e10N'e_2 = e_1 \neq 0). Vuistregel: μ\mu en χ\chi delen hun wortels (Propositie 3.12); de multipliciteit in μ\mu meet de grootte van het grootste nilpotente blok, die in χ\chi de totale dimensie van de karakteristieke deelruimte.

Stelling 3.21 (Cayley–Hamilton)

χu(u)=0\chi_u(u) = 0; bijgevolg μuχu\mu_u \mid \chi_u, en degμun\deg \mu_u \leq n.

Bewijs. Fixeer x0x \neq 0 en zij dd maximaal met (x,u(x),,ud1(x))(x, u(x), \dots, u^{d-1}(x)) vrij; schrijf

ud(x)=a0xa1u(x)ad1ud1(x),u^d(x) = -a_0 x - a_1 u(x) - \dots - a_{d-1}u^{d-1}(x),

en stel Px=Xd+ad1Xd1++a0P_x = X^d + a_{d-1}X^{d-1} + \dots + a_0, dus Px(u)(x)=0P_x(u)(x) = 0. Vul de vrije familie aan tot een basis van EE: daarop heeft uu blokvorm (C0D)\begin{pmatrix} C & \ast\\ 0 & D\end{pmatrix} waarbij CC de begeleidende matrix van PxP_x is, waarvan de karakteristieke veelterm PxP_x is (ontwikkel det(XIC)\det(XI - C) langs de eerste kolom, door inductie op dd). Dus χu=PxχD\chi_u = P_x \cdot \chi_D, en

χu(u)(x)=χD(u)(Px(u)(x))=0.\chi_u(u)(x) = \chi_D(u)\bigl(P_x(u)(x)\bigr) = 0 .

Het argument geldt voor elke xx: χu(u)=0\chi_u(u) = 0.

Voorbeeld 3.22 (Cayley–Hamilton aan het werk)

A=(1234)A = \begin{pmatrix} 1 & 2\\ 3 & 4\end{pmatrix}: χA=X25X2\chi_A = X^2 - 5X - 2, dus A2=5A+2IA^2 = 5A + 2I. Elke macht van AA stort in tot een combinatie van II en AA:

A4=(5A+2I)2=25A2+20A+4I=145A+54I=(199290435634),A^4 = (5A + 2I)^2 = 25A^2 + 20A + 4I = 145A + 54I = \begin{pmatrix} 199 & 290\\ 435 & 634\end{pmatrix},

en de inverse komt gratis: A(A5I)=2IA(A - 5I) = 2I geeft

A1=12(A5I)=(213/21/2).A^{-1} = \tfrac12(A - 5I) = \begin{pmatrix} -2 & 1\\ 3/2 & -1/2\end{pmatrix}.

Het sluitende inzicht: Cayley–Hamilton comprimeert heel de algebra K[A]K[A] tot Vect(I,A,,An1)\operatorname{Vect}(I, A, \dots, A^{n-1})dimK[A]=degμAn\dim K[A] = \deg\mu_A \leq n, hoe groot de machten die u nodig heeft ook zijn.

Opmerking 3.23 (Veelvoorkomende valkuilen)

(i) Eigenwaarden tellen niet op: Sp(A+B)\operatorname{Sp}(A + B) is niet SpA+SpB\operatorname{Sp}A + \operatorname{Sp}B, en een som van diagonaliseerbare matrices hoeft niet diagonaliseerbaar te zijn — (1100)+(0001)=(1101)\left(\begin{smallmatrix}1 & 1\\ 0 & 0\end{smallmatrix}\right) + \left(\begin{smallmatrix}0 & 0\\ 0 & 1\end{smallmatrix}\right) = \left(\begin{smallmatrix}1 & 1\\ 0 & 1\end{smallmatrix}\right) is een som van twee diagonaliseerbare matrices (elk heeft verschillende eigenwaarden) en is niet diagonaliseerbaar; alleen commuterende families gedragen zich (Oefening 3.9). (ii) “χu\chi_u splitst” is een hypothese over het lichaam: een vlakrotatie heeft χ=X22cosθX+1\chi = X^2 - 2\cos\theta\,X + 1, gesplitst over C\C, niet over R\Rdiagonaliseerbaar in M2(C)\mathcal{M}_2(\C), niet trigonaliseerbaar in M2(R)\mathcal{M}_2(\R). (iii) De ongelijkheid loopt meetkundig \leq algebraïsch, nooit omgekeerd; alleen dimEλ1\dim E_\lambda \geq 1 testen bewijst niets over diagonaliseerbaarheid. (iv) μu\mu_u is niet χu\chi_u: gelijkheid geldt precies wanneer elke eigenwaarde een enkele blokketen heeft (bijv. begeleidende matrices, het weekendprobleem van dit hoofdstuk); χ\chi gebruiken waar μ\mu nodig is blaast elke machtsberekening op. (v) Dunfords dd en ν\nu zijn veeltermen in uu — een ontbinding u=d+νu = d' + \nu' met de juiste eigenschappen maar dννdd'\nu' \neq \nu'd' is niet Dunford en is nooit uniek.

Opmerking 3.24 (Waar dit hoofdstuk wordt gebruikt)

Reductie is het werkpaard van de rest van het boek: machten en exponentialen van matrices drijven de lineaire differentiaalstelsels van Hoofdstuk 16; de spectraalstelling van Hoofdstuk 12 is diagonaliseren orthogonaal gemaakt; genererende functies (Hoofdstuk 23) herleiden de recurrentie-asymptotiek van het weekendprobleem van dit hoofdstuk analytisch. In het Bachelor jaar 3-volume loopt hetzelfde programma in oneindige dimensie: de spectraaltheorie van compacte zelfgeadjungeerde operatoren, waar eigenwaardenvolgordes eindige spectra vervangen, en de Perron–Frobenius-theorie van positieve matrices, die verklaart waarom dominante eigenwaarden van telproblemen positief en enkelvoudig zijn.

3.3 Nilpotenten en de Dunford-ontbinding

Propositie 3.25 (Nilpotente endomorfismen)

Voor uu met χu\chi_u gesplitst zijn de volgende equivalent: un=0u^n = 0; uk=0u^k = 0 voor zekere kk; Sp(u)={0}\operatorname{Sp}(u) = \{0\}; χu=Xn\chi_u = X^n; uu is trigonaliseerbaar met nulle diagonaal. Een nilpotent endomorfisme heeft μu=X(nilpotentie-index)\mu_u = X^{\text{(nilpotentie-index)}}, en index n\leq n.

Bewijs. uk=0u^k = 0 maakt elke eigenwaarde een wortel van XkX^k: spectrum {0}\{0\} (niet-leeg wanneer χ\chi splitst — over C\C altijd). Dan χu=Xn\chi_u = X^n (alle wortels nul) en Cayley–Hamilton geeft un=0u^n = 0; trigonalisatie (Stelling 3.9) zet nullen op de diagonaal (de diagonaal draagt de eigenwaarden). Omgekeerd, zij AA strikt bovendriehoekig: aij=0a_{ij} = 0 voor jij \leq i. We tonen door inductie dat

(Ak)ij=0wanneer ji+k1,(A^k)_{ij} = 0 \qquad \text{wanneer } j \leq i + k - 1,

d.w.z. elke macht duwt de nulle regio één diagonaal hoger. Voor k=1k = 1 is dit de hypothese. Voor de stap,

(Ak+1)ij=(Ak)iaj,(A^{k+1})_{ij} = \sum_{\ell} (A^k)_{i\ell}\,a_{\ell j} ,

en elke term verdwijnt: ofwel i+k1\ell \leq i + k - 1 (de eerste factor is 00 door inductie) of i+k\ell \geq i + k, in welk geval ji+kj \leq i + k \leq \ell de tweede factor doodt. Bij k=nk = n geldt de voorwaarde ji+n1j \leq i + n - 1 voor alle i,jni, j \leq n: An=0A^n = 0. De minimale veelterm deelt XnX^n en annihilatie definieert de index.

Stelling 3.26 (Dunford-ontbinding)

Stel χu\chi_u splitst over KK (automatisch voor K=CK = \C). Dan bestaat er een uniek paar (d,ν)(d, \nu) met

u=d+ν,d diagonaliseerbaar,ν nilpotent,dν=νd,u = d + \nu, \qquad d \text{ diagonaliseerbaar}, \quad \nu \text{ nilpotent}, \quad d\nu = \nu d ,

en bovendien zijn dd en ν\nu veeltermen in uu.

Bewijs. Bestaan. Schrijf χu=i=1r(Xλi)mi\chi_u = \prod_{i=1}^{r} (X - \lambda_i)^{m_i} (verschillende λi\lambda_i) en stel Ni=ker(uλi)miN_i = \ker(u - \lambda_i)^{m_i}, de karakteristieke deelruimten. Door Cayley–Hamilton en het kernlemma (Stelling 3.14),

E=N1Nr,E = N_1 \oplus \dots \oplus N_r ,

met projecties πi\pi_i veelterm in uu; elke NiN_i is stabiel (veeltermen in uu commuten met uu). Definieer d=iλiπid = \sum_i \lambda_i \pi_i: een veelterm in uu, diagonaliseerbaar (ze werkt als λi\lambda_i op NiN_i, dus ontbindt EE in haar eigenruimten). Dan is ν=ud\nu = u - d een veelterm in uu (dus commuteert met dd), en op elke NiN_i werkt ze als uλiu - \lambda_i, met (uλi)mi=0(u - \lambda_i)^{m_i} = 0 daar: νmaxmi=0\nu^{\max m_i} = 0 op elke sommand, dus is ν\nu nilpotent.

Uniciteit. Zij u=d+νu = d' + \nu' een ander zulk paar. Omdat dd' en ν\nu' met elkaar commuten, commuten ze met u=d+νu = d' + \nu', dus met elke veelterm in uu — in het bijzonder met dd en ν\nu. Dan is ddd - d' diagonaliseerbaar (twee commuterende diagonaliseerbare afbeeldingen zijn gelijktijdig diagonaliseerbaar: Oefening 3.9) en gelijk aan νν\nu' - \nu, die nilpotent is: als νk=0\nu^k = 0 en νk=0\nu'^{k'} = 0, laat commutatie de binomiaalontwikkeling toe

(νν)k+k1=j=0k+k1(k+k1j)νj(ν)k+k1j,(\nu' - \nu)^{k + k' - 1} = \sum_{j=0}^{k+k'-1}\binom{k + k' - 1}{j} \,\nu'^{\,j}\,(-\nu)^{k + k' - 1 - j} ,

waarin elke term sterft: ofwel jkj \geq k' (eerste factor nul) of k+k1jkk + k' - 1 - j \geq k (tweede factor nul), en één van de twee geldt altijd. Een diagonaliseerbare nilpotente is nul (haar spectrum is {0}\{0\} en ze is diagonaal in zekere basis): d=dd = d', ν=ν\nu = \nu'.

Voorbeeld 3.27 (Machten en exponentialen)

A=(3111)A = \begin{pmatrix} 3 & 1\\ -1 & 1\end{pmatrix}: χA=X24X+4=(X2)2\chi_A = X^2 - 4X + 4 = (X-2)^2, enkele eigenwaarde 22, eigenruimte van dimensie 11: niet diagonaliseerbaar. Dunford: D=2ID = 2I, N=A2I=(1111)N = A - 2I = \begin{pmatrix} 1 & 1\\ -1 & -1\end{pmatrix}, N2=0N^2 = 0. Dan

Ak=(2I+N)k=2kI+k2k1N,etA=e2t(I+tN),A^k = (2I + N)^k = 2^k I + k\,2^{k-1} N , \qquad \eu^{tA} = \eu^{2t}(I + tN),

door de commuterende binomiaalstelling, resp. de exponentiaalreeks (Hoofdstuk 16) gesplitst op commuterende sommanden. Reductie zet matrixdynamica om in scalaire dynamica.

Opmerking 3.28 (Perspectieven binnen dit volume)

Reductie is een knooppunt; hier zijn de vier spaken om te volgen. In Hoofdstuk 5 zetten aangepaste normen “alle eigenwaarden van modulus <1< 1” om in “zekere operatornorm <1< 1”, waardoor spectra convergentie van machten en reeksen regeren. In Hoofdstuk 16 wordt het recept van Voorbeeld 3.27 de algemene oplossing van X=AXX' = AX: Dunford splijt etA\eu^{tA} in veelterm-maal-exponentiaal-blokken, en stabiliteit leest de reële delen van eigenwaarden af. In Hoofdstuk 12 dwingt een scalair product wat loutere lineaire algebra niet kan: symmetrische matrices worden orthogonaal diagonaliseerbaar, zonder nilpotent deel. En in Hoofdstuk 23 keren de dominante-eigenwaarde-asymptotiek van het weekendprobleem van dit hoofdstuk analytisch terug, als de kleinste singulariteit van een genererende functie — twee talen voor één groeisnelheid.

3.4 Oefeningen

Oefening 3.1

Diagonaliseer (eigenwaarden, bases van eigenruimten, inverteerbare PP):

A=(1221),B=(011101110).A = \begin{pmatrix} 1 & 2\\ 2 & 1 \end{pmatrix}, \qquad B = \begin{pmatrix} 0 & 1 & 1\\ 1 & 0 & 1\\ 1 & 1 & 0 \end{pmatrix}.
Oplossing

Oplossing van Oefening 3.1.

AA: χA=X22X3=(X3)(X+1)\chi_A = X^2 - 2X - 3 = (X - 3)(X + 1). Eigenvectoren: voor 33: (1,1)(1,1); voor 1-1: (1,1)(1,-1). Dus P=(1111)P = \begin{pmatrix} 1 & 1\\ 1 & -1\end{pmatrix} geeft P1AP=diag(3,1)P^{-1}AP = \operatorname{diag}(3, -1).

B=JIB = J - I waarbij JJ de all-eens-matrix is. JJ heeft rang 11 met Jv=3vJv = 3v voor v=(1,1,1)v = (1,1,1) en Jw=0Jw = 0 op het vlak x+y+z=0x + y + z = 0: spectrum van BB is {2,1}\{2, -1\} met eigenruimten Vect(1,1,1)\operatorname{Vect}(1,1,1) (dimensie 11) en {x+y+z=0}\{x + y + z = 0\} (dimensie 22, basis (1,1,0),(1,0,1)(1,-1,0), (1,0,-1)). PP met deze drie kolommen geeft P1BP=diag(2,1,1)P^{-1}BP = \operatorname{diag}(2, -1, -1).

Oefening 3.2

Toon dat C=(1101)C = \begin{pmatrix} 1 & 1\\ 0 & 1\end{pmatrix} niet diagonaliseerbaar is, tweemaal: via eigenruimten, en via de minimale veelterm.

Oplossing

Oplossing van Oefening 3.2.

Eigenruimten: χC=(X1)2\chi_C = (X-1)^2, enkele eigenwaarde 11; ker(CI)=ker(0100)\ker(C - I) = \ker\begin{pmatrix} 0&1\\ 0&0\end{pmatrix} is de lijn Vect(e1)\operatorname{Vect}(e_1): dimensie 1<2=m11 < 2 = m_1, dus niet diagonaliseerbaar (Stelling 3.6).

Minimale veelterm: μC\mu_C deelt (X1)2(X-1)^2 en CIC \neq I, dus μC=(X1)2\mu_C = (X-1)^2: een dubbele wortel, dus niet diagonaliseerbaar (Gevolg 3.17).

Oefening 3.3

Zij uu voldoen aan u25u+6id=0u^2 - 5u + 6\,\mathrm{id} = 0. Bewijs dat uu diagonaliseerbaar is, bepaal de mogelijke spectra, en bereken uku^k als combinatie van id\mathrm{id} en uu.

Oplossing

Oplossing van Oefening 3.3.

X25X+6=(X2)(X3)X^2 - 5X + 6 = (X-2)(X-3): gesplitst met enkelvoudige wortels, dus is uu diagonaliseerbaar (Gevolg 3.17), met Sp(u){2,3}\operatorname{Sp}(u) \subseteq \{2, 3\}. Mogelijke spectra: {2}\{2\} (u=2idu = 2\,\mathrm{id}), {3}\{3\} (u=3idu = 3\,\mathrm{id}), of {2,3}\{2, 3\}.

Machten: zoek uk=akid+bkuu^k = a_k\,\mathrm{id} + b_k\,u. Op de eigenruimten luidt dit 2k=ak+2bk2^k = a_k + 2b_k en 3k=ak+3bk3^k = a_k + 3b_k: oplossend, bk=3k2kb_k = 3^k - 2^k, ak=32k23ka_k = 3\cdot2^k - 2\cdot 3^k:

uk=(32k23k)id+(3k2k)u.u^k = (3\cdot 2^k - 2\cdot 3^k)\,\mathrm{id} + (3^k - 2^k)\, u .

(Geldig voor alle drie de spectra: de identiteiten gelden eigenwaarde-gewijs.)

Oefening 3.4 ★★

Zij uu diagonaliseerbaar en FF een stabiele deelruimte. Bewijs dat uFu|_F diagonaliseerbaar is (beperk een annulerende veelterm met enkelvoudige gesplitste wortels).

Oplossing

Oplossing van Oefening 3.4.

uu diagonaliseerbaar: P=λ(Xλ)P = \prod_{\lambda}(X - \lambda) over het spectrum annuleert uu, splitst, enkelvoudige wortels. Dan P(uF)=P(u)F=0P(u|_F) = P(u)|_F = 0: de restrictie wordt geannuleerd door een gesplitste veelterm met enkelvoudige wortels, dus diagonaliseerbaar (Gevolg 3.17).

Oefening 3.5 ★★

(Fibonacci) Zij A=(1110)A = \begin{pmatrix} 1 & 1\\ 1 & 0\end{pmatrix}. Diagonaliseer AA over R\R, en leid de formule van Binet af voor de Fibonacci-rij (F0=0F_0 = 0, F1=1F_1 = 1, Fn+1=Fn+Fn1F_{n+1} = F_n + F_{n-1}):

Fn=φnψn5,φ=1+52, ψ=152.F_n = \frac{\varphi^n - \psi^n}{\sqrt 5}, \qquad \varphi = \frac{1 + \sqrt5}{2},\ \psi = \frac{1 - \sqrt5}{2}.
Oplossing

Oplossing van Oefening 3.5.

χA=X2X1\chi_A = X^2 - X - 1, wortels φ\varphi en ψ\psi (verschillend): diagonaliseerbaar, met eigenvectoren (φ,1)(\varphi, 1) en (ψ,1)(\psi, 1). De recurrentie geeft (Fn+1Fn)=An(10)\begin{pmatrix} F_{n+1}\\ F_n \end{pmatrix} = A^n \begin{pmatrix}1\\ 0\end{pmatrix}. Ontbind (1,0)(1, 0) op de eigenvectoren: (1,0)=1φψ((φ,1)(ψ,1))(1,0) = \frac{1}{\varphi - \psi}\bigl((\varphi, 1) - (\psi, 1)\bigr) met φψ=5\varphi - \psi = \sqrt5. AnA^n toepassen vermenigvuldigt elke eigencomponent met de nn-de macht van haar eigenwaarde; de tweede coördinaat lezend:

Fn=φnψn5.F_n = \frac{\varphi^n - \psi^n}{\sqrt 5} .

(Controle: n=1n = 1 geeft φψ5=1\frac{\varphi - \psi}{\sqrt5} = 1.)

Oefening 3.6 ★★

Zij uL(E)u \in \mathcal{L}(E) met u2u^2 diagonaliseerbaar en uu inverteerbaar (K=CK = \C). Bewijs dat uu diagonaliseerbaar is. Geef een tegenvoorbeeld wanneer uu niet inverteerbaar is.

Oplossing

Oplossing van Oefening 3.6.

Zij P=i(Xμi)P = \prod_i (X - \mu_i) u2u^2 annuleert, gesplitst met enkelvoudige wortels μi\mu_i (het spectrum van u2u^2). Omdat uu inverteerbaar is, is 00 geen eigenwaarde van u2u^2 (detu2=(detu)20\det u^2 = (\det u)^2 \neq 0), dus alle μi0\mu_i \neq 0. Dan

Q(X)=i(X2μi)=i(Xμi)(X+μi)Q(X) = \prod_i (X^2 - \mu_i) = \prod_i (X - \sqrt{\mu_i})(X + \sqrt{\mu_i})

annuleert uu:   Q(u)=i(u2μiid)=P(u2)=0\;Q(u) = \prod_i (u^2 - \mu_i\,\mathrm{id}) = P(u^2) = 0. Haar wortels ±μi\pm \sqrt{\mu_i} (complexe vierkantswortels) zijn paarsgewijs verschillend omdat de μi\mu_i verschillend en niet-nul zijn (μi=μj\sqrt{\mu_i} = -\sqrt{\mu_j} zou μi=μj\mu_i = \mu_j geven). Gesplitst + enkelvoudige wortels: uu is diagonaliseerbaar.

Tegenvoorbeeld zonder inverteerbaarheid: u=(0100)u = \begin{pmatrix} 0 & 1\\ 0 & 0\end{pmatrix}: u2=0u^2 = 0 is diagonaliseerbaar, uu is het niet.

Oefening 3.7 ★★

Bereken de Dunford-ontbinding, AkA^k, en etA\eu^{tA} voor

A=(210021002).A = \begin{pmatrix} 2 & 1 & 0\\ 0 & 2 & 1\\ 0 & 0 & 2 \end{pmatrix}.
Oplossing

Oplossing van Oefening 3.7.

A=2I+NA = 2I + N met NN de shift (Ne2=e1N e_2 = e_1, Ne3=e2Ne_3 = e_2), N3=0N^3 = 0, N2=E13N^2 = E_{13}: dit is de Dunford-ontbinding (2I2I diagonaal, NN nilpotent, ze commuten; uniciteit maakt haar de unieke). Binomiaal met commuterende termen:

Ak=2kI+k2k1N+(k2)2k2N2=(2kk2k1(k2)2k202kk2k1002k),A^k = 2^k I + k 2^{k-1} N + \binom k2 2^{k-2} N^2 = \begin{pmatrix} 2^k & k2^{k-1} & \binom k2 2^{k-2}\\ 0 & 2^k & k2^{k-1}\\ 0 & 0 & 2^k \end{pmatrix},
etA=e2t(I+tN+t22N2)=e2t(1tt2/201t001).\eu^{tA} = \eu^{2t}\Bigl(I + tN + \frac{t^2}{2}N^2\Bigr) = \eu^{2t}\begin{pmatrix} 1 & t & t^2/2\\ 0 & 1 & t\\ 0 & 0 & 1 \end{pmatrix}.

Oefening 3.8 ★★

Zij AMn(C)A \in \mathcal{M}_n(\C) met Ak=IA^k = I voor zekere k1k \geq 1. Bewijs dat AA diagonaliseerbaar is en haar eigenwaarden kk-de eenheidswortels zijn. Leid af dat een eindige-orde inverteerbare complexe matrix gelijkvormig met een driehoekige matrix met eenheidsdiagonaal de identiteit is.

Oplossing

Oplossing van Oefening 3.8.

Xk1X^k - 1 annuleert AA en splitst over C\C met de kk verschillende wortels e2iπj/k\eu^{2\iu\pi j/k}: AA is diagonaliseerbaar (Gevolg 3.17) en haar eigenwaarden, wortels van Xk1X^k - 1, zijn kk-de eenheidswortels.

Als bovendien AA gelijkvormig is met een driehoekige matrix met eenheidsdiagonaal: alle eigenwaarden gelijk 11, en AA, diagonaliseerbaar met enige eigenwaarde 11, is PIP1=IP\,I\,P^{-1} = I.

Oefening 3.9 ★★★

(Gelijktijdige diagonaliseerbaarheid) Zij u,vu, v diagonaliseerbaar en commuterend. Bewijs dat ze gelijktijdig diagonaliseerbaar zijn: zekere basis diagonaliseert beide. (Elke eigenruimte van uu is vv-stabiel; diagonaliseer de restricties van vv daar, met Oefening 3.4.)

Oplossing

Oplossing van Oefening 3.9.

Schrijf E=λEλ(u)E = \bigoplus_\lambda E_\lambda(u) (Stelling 3.6). Elke Eλ(u)E_\lambda(u) is vv-stabiel: voor xEλx \in E_\lambda, u(v(x))=v(u(x))=λv(x)u(v(x)) = v(u(x)) = \lambda v(x). De restrictie van vv tot Eλ(u)E_\lambda(u) is diagonaliseerbaar (Oefening 3.4): kies een basis van Eλ(u)E_\lambda(u) bestaande uit vv-eigenvectoren. Concatenatie van deze bases over alle λ\lambda geeft een basis van EE waarvan de vectoren eigenvectoren zijn van beide uu (door lidmaatschap van Eλ(u)E_\lambda(u)) en vv (per constructie).

Oefening 3.10 ★★★

Zij uL(Cn)u \in \mathcal{L}(\C^n). Bewijs dat uu diagonaliseerbaar is desda elke uu-stabiele deelruimte een uu-stabiele complementaire deelruimte heeft. (Voor \Leftarrow: pas de eigenschap toe op F=λEλ(u)F = \sum_\lambda E_\lambda(u), de som van alle eigenruimten; als een stabiel complement GG niet-nul was, zou trigonaliseren van uGu|_G een eigenvector van uu in GG produceren — in strijd met GF={0}G \cap F = \{0\}.)

Oplossing

Oplossing van Oefening 3.10.

(\Rightarrow) Zij uu diagonaliseerbaar en FF stabiel. Dan is uFu|_F diagonaliseerbaar (Oefening 3.4): FF heeft een basis van eigenvectoren, die zich, binnen elke globale eigenruimte EλE_\lambda, uitbreidt tot een basis van EλE_\lambda (onvolledige-basisstelling binnen EλE_\lambda, startend vanuit het deel van de basis van FF dat daar ligt — merk op F=λ(FEλ)F = \bigoplus_\lambda (F \cap E_\lambda) omdat uFu|_F diagonaliseerbaar is). De toegevoegde vectoren spannen een stabiel complement op (elk ligt in zekere EλE_\lambda, dus is hun opspansel uu-stabiel).

(\Leftarrow) Zij F=λEλ(u)F = \sum_\lambda E_\lambda(u) (een stabiele deelruimte) en GG een stabiel complement. Als G{0}G \neq \{0\}: χuG\chi_{u|_G} splitst over C\C, dus heeft uGu|_G een eigenvector xGx \in G (Stelling 3.9 of direct bestaan van een wortel); maar elke eigenvector van uu ligt in FF, dus xFG={0}x \in F \cap G = \{0\}: contradictie. Dus G={0}G = \{0\} en E=FE = F: de eigenruimten vullen EE, d.w.z. uu is diagonaliseerbaar.

Oefening 3.11 ★★★

(Spectrale straal à la Gelfand-lite, 2×22\times2-smaak van de analyse die komt) Zij AM2(C)A \in \mathcal{M}_2(\C) met beide eigenwaarden van modulus <1< 1. Bewijs dat Ak0A^k \to 0 entrygewijs als kk \to \infty. (Trigonaliseer: A=P(T)P1A = P(T)P^{-1} met TT bovendriehoekig; bereken TkT^k expliciet — onderscheid gelijke en verschillende eigenwaarden — en begrens.)

Oplossing

Oplossing van Oefening 3.11.

Trigonaliseer: A=PTP1A = PTP^{-1}, T=(λc0μ)T = \begin{pmatrix} \lambda & c\\ 0 & \mu\end{pmatrix}, λ,μ<1\abs\lambda, \abs\mu < 1. Dan Ak=PTkP1A^k = PT^kP^{-1}, en het volstaat dat Tk0T^k \to 0.

Verschillende eigenwaarden: inductie geeft

Tk=(λkcλkμkλμ0μk),T^k = \begin{pmatrix} \lambda^k & c\,\dfrac{\lambda^k - \mu^k}{\lambda - \mu}\\[4pt] 0 & \mu^k \end{pmatrix},

en elke entry neigt naar 00 (λk,μk0\abs{\lambda}^k, \abs\mu^k \to 0).

Gelijke eigenwaarden (μ=λ\mu = \lambda): T=λI+cE12T = \lambda I + cE_{12} en Tk=λkI+kλk1cE12T^k = \lambda^k I + k\lambda^{k-1}cE_{12}; de entry kλk10k\lambda^{k-1} \to 0 omdat λ<1\abs\lambda < 1 (geometrisch verslaat polynomiaal). In beide gevallen Tk0T^k \to 0 entrygewijs, dus Ak=PTkP10A^k = PT^kP^{-1} \to 0 (matrixvermenigvuldiging met vaste P,P1P, P^{-1} is continu in de entries — elke entry van het product is een vaste lineaire combinatie).

Oefening 3.12 ★★

Zij uL(Cn)u \in \mathcal{L}(\C^n) met rku=1\operatorname{rk} u = 1 (n2n \geq 2). Toon dat χu=Xn1(Xtru)\chi_u = X^{n-1}(X - \operatorname{tr} u), en dat uu diagonaliseerbaar is desda tru0\operatorname{tr} u \neq 0. (Herinner uit Oefening 2.5 dat u2=(tru)uu^2 = (\operatorname{tr} u)\,u.)

Oplossing

Oplossing van Oefening 3.12.

keru\ker u heeft dimensie n1n - 1 (rang–nulheid), dus is 00 een eigenwaarde van meetkundige multipliciteit n1n - 1, en is χu\chi_u deelbaar door Xn1X^{n-1} (Definitie 3.3: meetkundig \leq algebraïsch). Schrijf χu=Xn1(Xα)\chi_u = X^{n-1}(X - \alpha); de coëfficiënt van Xn1X^{n-1} is tru-\operatorname{tr} u, dus α=tru\alpha = \operatorname{tr} u: χu=Xn1(Xtru)\chi_u = X^{n-1}(X - \operatorname{tr} u).

Als tru0\operatorname{tr} u \neq 0: de eigenwaarde tru\operatorname{tr} u is een wortel van χu\chi_u, dus draagt ze een eigenvector; de eigenruimten voor 00 en tru\operatorname{tr} u hebben dimensies n1n - 1 en 1\geq 1, sommerend tot n\geq n: ze vullen EE, en uu is diagonaliseerbaar (Stelling 3.6). Als tru=0\operatorname{tr} u = 0: door Oefening 2.5, u2=(tru)u=0u^2 = (\operatorname{tr} u)u = 0 met u0u \neq 0: uu is een niet-nulle nilpotente, en een diagonaliseerbare nilpotente is nul (Propositie 3.25): niet diagonaliseerbaar.

3.5 Probleem: Lineaire recurrenties en begeleidende matrices

Een lineaire recurrentie un+k=ak1un+k1++a0unu_{n+k} = a_{k-1}u_{n+k-1} + \dots + a_0 u_n is een matrixmacht in vermomming, en reductie zet ze om in gesloten formules, groeisnelheden, en foutschattingen. Dit weekendprobleem ontwikkelt het woordenboek — begeleidende matrices aan de ene kant, de shiftoperator op de rijruimte aan de andere — bewijst de fundamentele stelling van lineaire recurrenties (de algemene oplossing is iQi(n)λin\sum_i Q_i(n)\lambda_i^n over de wortels van de karakteristieke veelterm), en besteedt de dividenden aan Diophantische benadering van 2\sqrt2, aan het tellen van wandelingen en woorden, en aan een ring van gekoppelde rijen die alleen gelijktijdige diagonaliseerbaarheid kan ontwarren.

Probleem 3.1

Weekendprobleem — de fundamentele stelling van lineaire recurrenties

Fixeer k1k \geq 1, scalaren a0,,ak1Ca_0, \dots, a_{k-1} \in \C met a00a_0 \neq 0, de monische veelterm P=Xkak1Xk1a1Xa0P = X^k - a_{k-1}X^{k-1} - \dots - a_1 X - a_0, en de recurrentie

(R) ⁣:un+k=ak1un+k1++a1un+1+a0un(n0).(\mathcal R)\colon\quad u_{n+k} = a_{k-1}u_{n+k-1} + \dots + a_1 u_{n+1} + a_0 u_n \qquad (n \geq 0).

De begeleidende matrix van PP is

C=(0101a0a1ak1)Mk(C).C = \begin{pmatrix} 0 & 1 & & \\ & \ddots & \ddots & \\ & & 0 & 1\\ a_0 & a_1 & \cdots & a_{k-1} \end{pmatrix} \in \mathcal{M}_k(\C).

Deel I — Het begeleidende woordenboek.

  1. Toon dat een rij (un)(u_n) aan (R)(\mathcal R) voldoet desda de vectoren vn=(un,un+1,,un+k1)Tv_n = (u_n, u_{n+1}, \dots, u_{n+k-1})^{\mathsf T} voldoen aan vn+1=Cvnv_{n+1} = Cv_n, dus vn=Cnv0v_n = C^n v_0.
  2. Bewijs dat χC=P\chi_C = P (ontwikkel det(XIC)\det(XI - C) langs de eerste kolom en induceer op kk), dan dat μC=P\mu_C = P eveneens (ga over op CTC^{\mathsf T}, waarvoor e1e_1 cyclisch is, en merk op dat een matrix en haar getransponeerde dezelfde minimale veelterm hebben).
  3. Toon dat voor elke wortel λ\lambda van PP de vector (1,λ,,λk1)T(1, \lambda, \dots, \lambda^{k-1})^{\mathsf T} de eigenruimte van CC voor λ\lambda opspant; leid af dat elke eigenruimte van CC dimensie 11 heeft, en dat CC diagonaliseerbaar is desda PP kk verschillende wortels heeft.
  4. Stel PP heeft verschillende wortels λ1,,λk\lambda_1, \dots, \lambda_k. Toon dat de meetkundige rijen (λin)n(\lambda_i^n)_n een basis vormen van de oplossingsruimte van (R)(\mathcal R), dus is elke oplossing un=iciλinu_n = \sum_i c_i\lambda_i^n voor unieke constanten cic_i.
  5. Los volledig op: un+2=un+1+6unu_{n+2} = u_{n+1} + 6u_n, u0=1u_0 = 1, u1=8u_1 = 8.

Deel II — De shiftoperator en de fundamentele stelling. Zij S\mathcal{S} de C\C-vectorruimte van alle complexe rijen en SL(S)S \in \mathcal{L}(\mathcal{S}) de shift, S((un)n)=(un+1)nS\bigl((u_n)_n\bigr) = (u_{n+1})_n.

  1. Toon dat de oplossingsverzameling van (R)(\mathcal R) kerP(S)\ker P(S) is, en dat ze precies dimensie kk heeft (beeld een oplossing af op haar beginwaarden).
  2. Leg uit waarom het kernontbindingslemma (Stelling 3.14) van toepassing is op SS op de oneindigdimensionale S\mathcal{S} zonder enige verandering, en schrijf de resulterende ontbinding van kerP(S)\ker P(S) voor P=i=1r(Xλi)miP = \prod_{i=1}^{r}(X - \lambda_i)^{m_i} (verschillende λi\lambda_i, alle niet-nul omdat a00a_0 \neq 0).
  3. Voor λ0\lambda \neq 0 en m1m \geq 1, toon

    ker(Sλid)m={(Q(n)λn)n:QCm1[X]},\ker\,(S - \lambda\,\mathrm{id})^m = \bigl\{\,\bigl(Q(n)\,\lambda^n\bigr)_n : Q \in \C_{m-1}[X]\,\bigr\},

    van dimensie mm. (Bereken (Sλ)(Q(n)λn)=λn+1(ΔQ)(n)(S - \lambda)\bigl(Q(n)\lambda^n\bigr) = \lambda^{n+1}(\Delta Q)(n) met ΔQ=Q(X+1)Q(X)\Delta Q = Q(X + 1) - Q(X), en gebruik dat Δ\Delta de graad verlaagt; voor de dimensie, begrens ze met mm via beginwaarden.)

  4. (De fundamentele stelling van lineaire recurrenties) Concludeer: als P=i=1r(Xλi)miP = \prod_{i=1}^{r}(X - \lambda_i)^{m_i} met de λi\lambda_i verschillend en niet-nul, zijn de oplossingen van (R)(\mathcal R) precies de rijen

    un=i=1rQi(n)λin,QiCmi1[X],u_n = \sum_{i=1}^{r} Q_i(n)\,\lambda_i^n, \qquad Q_i \in \C_{m_i - 1}[X],

    met uniek bepaalde veeltermen QiQ_i.

  5. Los volledig op: un+2=4un+14unu_{n+2} = 4u_{n+1} - 4u_n, u0=1u_0 = 1, u1=0u_1 = 0, en controleer het antwoord op u2u_2.

Deel III — Dominante wortels en Diophantische dividenden.

  1. Stel de wortels zijn enkelvoudig met λ1>λi\abs{\lambda_1} > \abs{\lambda_i} voor i2i \geq 2, en un=iciλinu_n = \sum_i c_i \lambda_i^n met c10c_1 \neq 0. Toon unc1λ1nu_n \sim c_1\lambda_1^n en un+1/unλ1u_{n+1}/u_n \to \lambda_1.
  2. (Pell) Definieer an+1=an+2bna_{n+1} = a_n + 2b_n, bn+1=an+bnb_{n+1} = a_n + b_n, a0=b0=1a_0 = b_0 = 1. Toon dat q(a,b)=a22b2q(a, b) = a^2 - 2b^2 voldoet aan q(an+1,bn+1)=q(an,bn)q(a_{n+1}, b_{n+1}) = -q(a_n, b_n), dus an22bn2=(1)n+1a_n^2 - 2b_n^2 = (-1)^{n+1}; relateer dit aan de determinant van M=(1211)M = \left(\begin{smallmatrix}1 & 2\\ 1 & 1\end{smallmatrix}\right).
  3. Leid de foutschatting af

    anbn2=1bn(an+2bn)12bn2,\Bigl|\frac{a_n}{b_n} - \sqrt2\Bigr| = \frac{1}{b_n\,(a_n + \sqrt2\,b_n)} \leq \frac{1}{2b_n^2},

    en toon dat ze geometrisch vervalt met ratio 3223 - 2\sqrt2 (vind de eigenwaarden van MM en de groei van bnb_n).

  4. (Algemene groei) Uit vraag 9, bewijs: (a) als elke wortel λiρ\abs{\lambda_i} \leq \rho voldoet, dan unCnm1ρn\abs{u_n} \leq C\,n^{m-1}\rho^n met m=maximim = \max_i m_i; (b) als er een unieke wortel λ1\lambda_1 van maximale modulus is en Q10Q_1 \neq 0, dan un+1/unλ1u_{n+1}/u_n \to \lambda_1 — controleer op de oplossing van vraag 10.

Deel IV — Wandelen en woorden tellen. Voor een eindige graaf met knopenverzameling {1,,N}\{1, \dots, N\} heeft de adjacentiematrix AA Aij=1A_{ij} = 1 als ijij een rand is, anders 00.

  1. Bewijs dat (An)ij(A^n)_{ij} het aantal wandelingen van lengte nn van ii naar jj is (rijen van nn randen, elke stap langs een rand).
  2. (De driehoek) Voor de complete graaf op 33 knopen, A=JIA = J - I: met het spectrum van JJ (Voorbeeld 2.19), toon

    (An)ii=2n+2(1)n3,(An)ij=2n(1)n3(ij),(A^n)_{ii} = \frac{2^n + 2(-1)^n}{3}, \qquad (A^n)_{ij} = \frac{2^n - (-1)^n}{3} \quad (i \neq j),

    en controleer beide bij n=2n = 2 door wandelingen op te sommen.

  3. (Woorden zonder 1111) Zij wnw_n het aantal binaire woorden van lengte nn zonder twee opeenvolgende 11’en. Codeer woorden door hun laatste letter om een transfermatrix te krijgen, toon wn+2=wn+1+wnw_{n+2} = w_{n+1} + w_n, leid af wn=Fn+2w_n = F_{n+2} (Fibonacci, Oefening 3.5), en geef de groeisnelheid limwn+1/wn\lim w_{n+1}/w_n.
  4. (Het pad) Voor de padgraaf 1231 - 2 - 3, toon dat de eigenwaarden van AA 2,0,2\sqrt2, 0, -\sqrt2 zijn met eigenvectoren (1,±2,1)(1, \pm\sqrt2, 1) en (1,0,1)(1, 0, -1), en leid af dat het aantal wandelingen van lengte nn van eind tot eind ((2)n+(2)n)/4\bigl((\sqrt2)^n + (-\sqrt2)^n\bigr)/4 is: nul voor oneven nn, en 2n/212^{\,n/2 - 1} voor even nn. Controleer bij n=4n = 4.
  5. (Spoorformule) Toon dat het totale aantal gesloten wandelingen van lengte nn (alle startpunten) tr(An)=iλin\operatorname{tr}(A^n) = \sum_i \lambda_i^n is, en verifieer op de driehoek.

Deel V — Een ring van rijen: gelijktijdige diagonaliseerbaarheid. Fixeer k3k \geq 3, zij ω=e2iπ/k\omega = \eu^{2\iu\pi/k}, en zij WMk(C)W \in \mathcal{M}_k(\C) de cyclische shift: Wei=ei+1W e_i = e_{i+1} (indices mod kk, kolommen geïndexeerd 0,,k10, \dots, k-1).

  1. Toon dat WTW^{\mathsf T} de begeleidende matrix van Xk1X^k - 1 is, leid af χW=μW=Xk1\chi_W = \mu_W = X^k - 1, en dat WW diagonaliseerbaar is met de kk enkelvoudige eigenwaarden ωj\omega^j en eigenvectoren fj=(1,ωj,ω2j,,ω(k1)j)Tf_j = (1, \omega^{-j}, \omega^{-2j}, \dots, \omega^{-(k-1)j})^{\mathsf T}.
  2. Een circulante matrix is C=c0I+c1W++ck1Wk1C = c_0 I + c_1 W + \dots + c_{k-1}W^{k-1}. Toon dat alle circulanten commuten, dat de basis (f0,,fk1)(f_0, \dots, f_{k-1}) alle ervan gelijktijdig diagonaliseert, en dat de eigenwaarden van CC c^(ωj)=mcmωjm\widehat c(\omega^j) = \sum_m c_m \omega^{jm} zijn, j=0,,k1j = 0, \dots, k-1.
  3. Leid af detC=j=0k1c^(ωj)\det C = \prod_{j=0}^{k-1} \widehat c(\omega^j), en controleer dat k=3k = 3 de factorisatie van Oefening 2.8 herwint.
  4. (Het kettinggemiddelde) Zij x(n+1)=Mx(n)x^{(n+1)} = Mx^{(n)} met M=12(W+W1)M = \frac12(W + W^{-1}): elk van kk getallen in een ring gerangschikt wordt vervangen door het gemiddelde van zijn twee buren. Toon dat de eigenwaarden van MM cos(2πj/k)\cos(2\pi j/k) zijn, en dat de coëfficiënt van x(0)x^{(0)} op f0f_0 het gemiddelde 1kmxm(0)\frac1k\sum_m x^{(0)}_m is (sommeer de coördinaten van de fjf_j).
  5. Concludeer: voor oneven kk convergeert x(n)x^{(n)} naar de constante vector waarvan de waarde het gemiddelde van de beginwaarden is; voor k=4k = 4, vertoon de eigenwaarde verantwoordelijk voor niet-convergentie en de exacte obstructie (een alternerend-gemiddelde- coëfficiënt die moet verdwijnen).
  6. (Synthese) In één zin elk: hoe de begeleidende matrix analyse van (R)(\mathcal R) omzet in reductie; waar het kernontbindingslemma geen eindige dimensie nodig had; waarom dominante eigenwaarden groeisnelheden en Diophantische fout regeren; waarom machten van de adjacentiematrix wandelingen tellen; en wat commuterende matrices kopen. Noem de twee toppen: de fundamentele stelling van lineaire recurrenties, en — voor de positieve matrices van Deel IV, in het Bachelor jaar 3-volume — de stelling van Perron–Frobenius.
Oplossing

Oplossing van Probleem 3.1.

1. De eerste k1k - 1 coördinaten van CvnCv_n zijn un+1,,un+k1u_{n+1}, \dots, u_{n+k-1} (de superdiagonaalverschuivingen), en de laatste is a0un++ak1un+k1a_0 u_n + \dots + a_{k-1}u_{n+k-1}. Dus vn+1=Cvnv_{n+1} = Cv_n geldt voor alle nn desda de laatste coördinaten voor alle nn matchen, d.w.z. desda (R)(\mathcal R) geldt. Itererend, vn=Cnv0v_n = C^nv_0.

2. Ontwikkel Dk(X)=det(XIkC)D_k(X) = \det(XI_k - C) langs de eerste kolom: de twee niet-nulle entries zijn XX (positie (1,1)(1,1)) en a0-a_0 (positie (k,1)(k,1)). De eerste minor is Dk1D_{k-1}-vormig voor de coëfficiënten a1,,ak1a_1, \dots, a_{k-1}; de tweede minor is bovendriehoekig met diagonaal 1-1: determinant (1)k1(-1)^{k-1}, met teken (1)k+1(-1)^{k+1} van de positie. Inductie op kk (basis k=1k = 1: Xa0X - a_0) geeft

Dk(X)=X(Xk1ak1Xk2a1)a0=P(X).D_k(X) = X\bigl(X^{k-1} - a_{k-1}X^{k-2} - \dots - a_1\bigr) - a_0 = P(X).

Voor μC\mu_C: omdat Q(CT)=Q(C)TQ(C^{\mathsf T}) = Q(C)^{\mathsf T} voor elke veelterm, hebben CC en CTC^{\mathsf T} dezelfde annihilatoren, dus dezelfde minimale veelterm. Voor CTC^{\mathsf T}: de kolommen lezen CTe1=e2C^{\mathsf T}e_1 = e_2, …, CTek1=ekC^{\mathsf T}e_{k-1} = e_k, dus is (e1,CTe1,,(CT)k1e1)(e_1, C^{\mathsf T}e_1, \dots, (C^{\mathsf T})^{k-1}e_1) de canonieke basis: vrij. Een veelterm Q0Q \neq 0 van graad <k< k heeft dan Q(CT)e10Q(C^{\mathsf T})e_1 \neq 0 (ze is een niet-triviale combinatie van basisvectoren): degμk\deg\mu \geq k. Omdat μχ=P\mu \mid \chi = P met degP=k\deg P = k: μC=P\mu_C = P.

3. Voor v=(1,λ,,λk1)Tv = (1, \lambda, \dots, \lambda^{k-1})^{\mathsf T}: rijen 11 tot k1k-1 van CvCv geven λ,λ2,,λk1\lambda, \lambda^2, \dots, \lambda^{k-1}, d.w.z. λ\lambda maal de eerste k1k - 1 entries van vv; de laatste rij geeft mamλm=λkP(λ)=λk=λλk1\sum_m a_m\lambda^m = \lambda^k - P(\lambda) = \lambda^k = \lambda\cdot\lambda^{k-1}. Dus Cv=λvCv = \lambda v. Omgekeerd luiden de vergelijkingen (Cx)i=λxi(Cx)_i = \lambda x_i voor i<ki < k xi+1=λxix_{i+1} = \lambda x_i: elke eigenvector is proportioneel met vv — elke eigenruimte heeft precies dimensie 11. Diagonaliseerbaar desda de eigenruimtedimensies tot kk sommeren (Stelling 3.6) desda er kk verschillende eigenwaarden zijn desda PP kk verschillende wortels heeft (de eigenwaarden zijn de wortels van χC=P\chi_C = P).

4. Elke (λin)n(\lambda_i^n)_n lost (R)(\mathcal R) op: λin+k=λinλik=λinmamλim\lambda_i^{n+k} = \lambda_i^n\,\lambda_i^k = \lambda_i^n\sum_m a_m\lambda_i^m. Vrijheid: een verdwijnende combinatie iciλin=0\sum_i c_i\lambda_i^n = 0 voor n=0,,k1n = 0, \dots, k-1 is een Vandermondestelsel (Oefening 2.11) in de cic_i: alle ci=0c_i = 0. De oplossingsruimte heeft dimensie kk (vraag 6, waarvan het bewijs elementair en onafhankelijk is): kk vrije oplossingen vormen een basis, en coördinaten zijn uniek.

5. P=X2X6=(X3)(X+2)P = X^2 - X - 6 = (X - 3)(X + 2): algemene oplossing un=A3n+B(2)nu_n = A\,3^n + B(-2)^n. Beginvoorwaarden: A+B=1A + B = 1, 3A2B=83A - 2B = 8: A=2A = 2, B=1B = -1:

un=23n(2)n.u_n = 2\cdot 3^n - (-2)^n .

(Controle: u2=u1+6u0=14u_2 = u_1 + 6u_0 = 14 en 294=142\cdot9 - 4 = 14.)

6. P(S)((un))P(S)\bigl((u_n)\bigr) is de rij nun+kak1un+k1a0unn \mapsto u_{n+k} - a_{k-1}u_{n+k-1} - \dots - a_0u_n: ze verdwijnt desda (R)(\mathcal R) geldt, dus is de oplossingsverzameling kerP(S)\ker P(S), een deelruimte. De afbeelding kerP(S)Ck\ker P(S) \to \C^k, u(u0,,uk1)u \mapsto (u_0, \dots, u_{k-1}), is lineair, injectief (de recurrentie bepaalt uk,uk+1,u_{k}, u_{k+1}, \dots uit de eerste kk waarden, door inductie) en surjectief (definieer unu_n recursief uit willekeurige beginwaarden): dimensie kk.

7. Het bewijs van Stelling 3.14 gebruikt alleen: de Bézout-identiteit in C[X]\C[X], en het feit dat veeltermen in een vast endomorfisme commuten. Geen van beide noemt de dimensie van de omgevingsruimte: het lemma geldt letterlijk voor SL(S)S \in \mathcal{L}(\mathcal{S}). Dus

kerP(S)=i=1rker(Sλiid)mi.\ker P(S) = \bigoplus_{i=1}^{r} \ker\,(S - \lambda_i\,\mathrm{id})^{m_i}.

8. Voor QC[X]Q \in \C[X]: (Sλ)(Q(n)λn)n(S - \lambda)\bigl(Q(n)\lambda^n\bigr)_n heeft nn-de term Q(n+1)λn+1λQ(n)λn=λn+1(ΔQ)(n)Q(n{+}1)\lambda^{n+1} - \lambda Q(n)\lambda^n = \lambda^{n+1}(\Delta Q)(n), met ΔQ=Q(X+1)Q(X)\Delta Q = Q(X{+}1) - Q(X) van graad degQ1\deg Q - 1 (de leidende termen cancelleren). Itererend, (Sλ)m(Q(n)λn)=(λn+m(ΔmQ)(n))n(S - \lambda)^m\bigl(Q(n)\lambda^n\bigr) = \bigl(\lambda^{n+m}(\Delta^m Q)(n)\bigr)_n, en ΔmQ=0\Delta^m Q = 0 wanneer degQm1\deg Q \leq m - 1: de rechterverzameling is bevat in de kern. Ze is een deelruimte van dimensie mm: de rijen (njλn)n(n^j\lambda^n)_n, 0j<m0 \leq j < m, zijn vrij, omdat jcjnjλn=0\sum_j c_j n^j\lambda^n = 0 voor alle nn (delen door λn0\lambda^n \neq 0) de veelterm jcjXj\sum_j c_jX^j dwingt te verdwijnen in elke nNn \in \N, dus nul te zijn. Omgekeerd dimker(Sλ)mm\dim\ker(S - \lambda)^m \leq m: ontwikkelend (Sλ)m=j(mj)(λ)mjSj(S - \lambda)^m = \sum_j \binom mj(-\lambda)^{m-j}S^j, is de vergelijking (Sλ)mu=0(S - \lambda)^m u = 0 een lineaire recurrentie van orde mm (leidende coëfficiënt 11), dus is uu bepaald door u0,,um1u_0, \dots, u_{m-1} zoals in vraag 6. Dimensiegelijkheid besluit.

9. Combineer vragen 7 en 8: elke oplossing ontbindt uniek als som van elementen van de ker(Sλi)mi\ker(S - \lambda_i)^{m_i}, d.w.z. un=iQi(n)λinu_n = \sum_i Q_i(n)\lambda_i^n met degQimi1\deg Q_i \leq m_i - 1; de QiQ_i zijn uniek omdat de ontbinding direct is en, binnen elke sommand, de coëfficiënten van QiQ_i coördinaten zijn in de basis (njλin)j(n^j\lambda_i^n)_j (vraag 8). Saniteitscontrole op dimensies: imi=k\sum_i m_i = k.

10. P=X24X+4=(X2)2P = X^2 - 4X + 4 = (X - 2)^2: oplossingen (a+bn)2n(a + bn)2^n. Beginwaarden: a=1a = 1, 2(a+b)=02(a + b) = 0, dus b=1b = -1:

un=(1n)2n.u_n = (1 - n)\,2^n .

Controle: u2=4u14u0=4u_2 = 4u_1 - 4u_0 = -4, en (12)4=4(1 - 2)\cdot4 = -4.

11. Schrijf un=λ1n(c1+i2ci(λi/λ1)n)u_n = \lambda_1^n\bigl(c_1 + \sum_{i\geq2} c_i(\lambda_i/\lambda_1)^n\bigr); elke ratio heeft modulus <1< 1, dus neigt de haak naar c10c_1 \neq 0: unc1λ1nu_n \sim c_1\lambda_1^n. In het bijzonder un0u_n \neq 0 voor grote nn, en

un+1un=λ1c1+o(1)c1+o(1)λ1.\frac{u_{n+1}}{u_n} = \lambda_1\,\frac{c_1 + o(1)}{c_1 + o(1)} \longrightarrow \lambda_1 .

12. Bereken:

q(an+1,bn+1)=(an+2bn)22(an+bn)2=an2+2bn2=q(an,bn).q(a_{n+1}, b_{n+1}) = (a_n + 2b_n)^2 - 2(a_n + b_n)^2 = -a_n^2 + 2b_n^2 = -q(a_n, b_n).

Met q(a0,b0)=12=1q(a_0, b_0) = 1 - 2 = -1: an22bn2=(1)n+1a_n^2 - 2b_n^2 = (-1)^{n+1}. Structureel: q(a,b)=(a2b)(a+2b)q(a, b) = (a - \sqrt2\,b)(a + \sqrt2\,b) en de lineaire afbeelding MM vermenigvuldigt de factor a+2ba + \sqrt2 b met 1+21 + \sqrt2 en de factor a2ba - \sqrt2 b met 121 - \sqrt2 (bereken: an+1+2bn+1=(1+2)(an+2bn)a_{n+1} + \sqrt2 b_{n+1} = (1 + \sqrt2)(a_n + \sqrt2 b_n)); het product wordt vermenigvuldigd met (1+2)(12)=1=detM(1 + \sqrt2)(1 - \sqrt2) = -1 = \det M bij elke stap.

13. Omdat an22bn2=(an2bn)(an+2bn)=(1)n+1a_n^2 - 2b_n^2 = (a_n - \sqrt2 b_n)(a_n + \sqrt2 b_n) = (-1)^{n+1},

anbn2=an22bn2bn(an+2bn)=1bn(an+2bn)12bn2,\Bigl|\frac{a_n}{b_n} - \sqrt2\Bigr| = \frac{\abs{a_n^2 - 2b_n^2}}{b_n(a_n + \sqrt2 b_n)} = \frac{1}{b_n(a_n + \sqrt2 b_n)} \leq \frac1{2b_n^2},

gebruikend anbn1a_n \geq b_n \geq 1 (inductie: beide stijgen) dus an+2bn(1+2)bn2bna_n + \sqrt2 b_n \geq (1 + \sqrt2)b_n \geq 2b_n. Eigenwaarden van MM: χM=X22X1\chi_M = X^2 - 2X - 1, wortels 1±21 \pm \sqrt2; omdat (a0,b0)(a_0, b_0) een niet-nulle component op de dominante eigenvector heeft (alle entries positief), bnc(1+2)nb_n \sim c(1 + \sqrt2)^n met c>0c > 0 (vraag 11). Dus is de fout (1+2)2n=(3+22)n\asymp (1 + \sqrt2)^{-2n} = (3 + 2\sqrt2)^{-n}: geometrisch verval met ratio 1/(3+22)=3220.1721/(3 + 2\sqrt2) = 3 - 2\sqrt2 \approx 0.172.

14. (a) Uit vraag 9: uniQi(n)λin(iQi(n))ρn\abs{u_n} \leq \sum_i \abs{Q_i(n)}\abs{\lambda_i}^n \leq \bigl(\sum_i \abs{Q_i(n)}\bigr)\rho^n, en elke Qi(n)Cinmi1Cinm1\abs{Q_i(n)} \leq C_i n^{m_i - 1} \leq C_i n^{m-1} voor n1n \geq 1: sommeer de constanten. (b) Zij ρ=maxi2λi<λ1\rho' = \max_{i \geq 2}\abs{\lambda_i} < \abs{\lambda_1} en d=degQ1d = \deg Q_1, leidende coëfficiënt c0c \neq 0. Dan un=Q1(n)λ1n+Rnu_n = Q_1(n)\lambda_1^n + R_n met RnCnm1ρn\abs{R_n} \leq Cn^{m-1}\rho'^n, en

RnQ1(n)λ1n=O(nm1d(ρ/λ1)n)0\frac{R_n}{Q_1(n)\lambda_1^n} = O\Bigl(n^{m-1-d} \bigl(\rho'/\abs{\lambda_1}\bigr)^n\Bigr) \longrightarrow 0

(geometrisch verslaat polynomiaal). Dus

unQ1(n)λ1ncndλ1n,un+1unλ1(omdat Q1(n+1)/Q1(n)1).u_n \sim Q_1(n)\,\lambda_1^n \sim c\,n^d\lambda_1^n, \qquad \frac{u_{n+1}}{u_n} \longrightarrow \lambda_1 \quad\text{(omdat } Q_1(n{+}1)/Q_1(n) \to 1\text{)}.

Controle op vraag 10: voor un=(1n)2nu_n = (1-n)2^n is de ratio

(n)2n+1(1n)2n=2n1n2=λ1.\frac{(-n)2^{n+1}}{(1-n)2^n} = 2\,\frac{-n}{1-n} \longrightarrow 2 = \lambda_1 .

15. Inductie op nn. Voor n=1n = 1 telt AijA_{ij} wandelingen van lengte 11. Stap: een wandeling van lengte n+1n + 1 van ii naar jj is een wandeling van lengte nn van ii naar zekere knoop \ell gevolgd door een rand j\ell j:

#{wandelingen}=(An)iAj=(An+1)ij.\#\{\text{wandelingen}\} = \sum_{\ell} (A^n)_{i\ell}A_{\ell j} = (A^{n+1})_{ij}.

16. J=3ΠJ = 3\Pi waarbij Π=J/3\Pi = J/3 de projectie is op Vect(1,1,1)\operatorname{Vect}(1,1,1) langs het vlak x+y+z=0x + y + z = 0 (Π2=Π\Pi^2 = \Pi omdat J2=3JJ^2 = 3J). Dan A=JI=2Π(IΠ)A = J - I = 2\Pi - (I - \Pi), en omdat Π\Pi en IΠI - \Pi complementaire projecties zijn,

An=2nΠ+(1)n(IΠ),d.w.z.(An)ij=2n3+(1)n(δij13),A^n = 2^n\,\Pi + (-1)^n (I - \Pi), \qquad\text{d.w.z.}\qquad (A^n)_{ij} = \frac{2^n}3 + (-1)^n\Bigl(\delta_{ij} - \frac13\Bigr),

wat de twee weergegeven formules geeft. Bij n=2n = 2: diagonaal (4+2)/3=2(4 + 2)/3 = 2 (wandelingen iii \to \ell \to i voor de twee buren \ell); buitendiagonaal (41)/3=1(4 - 1)/3 = 1 (de enkele wandeling iji \to \ell \to j door de derde knoop).

17. Zij wn(0),wn(1)w_n^{(0)}, w_n^{(1)} toegelaten woorden van lengte nn eindigend op 00, resp. 11 tellen. Een letter toevoegen: een 00 mag op alles volgen, een 11 alleen op een 00:

(wn+1(0)wn+1(1))=(1110)(wn(0)wn(1)).\begin{pmatrix} w_{n+1}^{(0)}\\ w_{n+1}^{(1)}\end{pmatrix} = \begin{pmatrix} 1 & 1\\ 1 & 0\end{pmatrix} \begin{pmatrix} w_n^{(0)}\\ w_n^{(1)}\end{pmatrix}.

Sommerend, wn+2=wn+1+wnw_{n+2} = w_{n+1} + w_n (of: conditioneer op de eerste letter). Met w1=2w_1 = 2, w2=3w_2 = 3: wn=Fn+2w_n = F_{n+2} door inductie (F3=2F_3 = 2, F4=3F_4 = 3, zelfde recurrentie). Groei: de wortels van X2X1X^2 - X - 1 zijn φ>ψ\varphi > \abs\psi (Oefening 3.5), en de φ\varphi-component is niet-nul (de wnw_n zijn positief en ψn0\psi^n \to 0), dus geeft vraag 11 wn+1/wnφ=1+52w_{n+1}/w_n \to \varphi = \frac{1 + \sqrt5}2.

18. A=(010101010)A = \left(\begin{smallmatrix} 0&1&0\\ 1&0&1\\ 0&1&0\end{smallmatrix}\right). Controle:

A(1,±2,1)T=(±2,2,±2)T=±2(1,±2,1)T,A(1,0,1)T=0:A(1, \pm\sqrt2, 1)^{\mathsf T} = (\pm\sqrt2, 2, \pm\sqrt2)^{\mathsf T} = \pm\sqrt2\,(1, \pm\sqrt2, 1)^{\mathsf T}, \qquad A(1, 0, -1)^{\mathsf T} = 0 :

eigenwaarden 2,2,0\sqrt2, -\sqrt2, 0 (=2cosπ4,2cos3π4,2cosπ2= 2\cos\frac\pi4, 2\cos\frac{3\pi}4, 2\cos\frac\pi2). Ontbind e1e_1 op de eigenbasis en lees de derde coördinaat, of gebruik symmetrie: met v±=(1,±2,1)v_\pm = (1, \pm\sqrt2, 1), v0=(1,0,1)v_0 = (1, 0, -1), controleert men e1=14v++14v+12v0e_1 = \frac14 v_+ + \frac14 v_- + \frac12 v_0, dus voor n1n \geq 1

(An)13=(14(2)nv++14(2)nv+0) ⁣3=(2)n+(2)n4,(A^n)_{13} = \Bigl(\tfrac14(\sqrt2)^n v_+ + \tfrac14(-\sqrt2)^n v_- + 0\Bigr)_{\!3} = \frac{(\sqrt2)^n + (-\sqrt2)^n}{4},

nul voor oneven nn (bipartiete graaf: einden liggen op even afstand), en 22n/2/4=2n/212\cdot 2^{n/2}/4 = 2^{n/2 - 1} voor even nn. Bij n=4n = 4: 21=22^{1} = 2, matchend de twee wandelingen 121231\,2\,1\,2\,3 en 123231\,2\,3\,2\,3.

19. Gesloten wandelingen van lengte nn vanuit ii zijn (An)ii(A^n)_{ii}; sommeren over ii geeft tr(An)\operatorname{tr}(A^n). Trigonaliserend AA (over C\C), is AnA^n driehoekig met diagonaal λin\lambda_i^n: tr(An)=iλin\operatorname{tr}(A^n) = \sum_i\lambda_i^n. Driehoek: tr(An)=32n+2(1)n3=2n+2(1)n=2n+(1)n+(1)n\operatorname{tr}(A^n) = 3\,\frac{2^n + 2(-1)^n}3 = 2^n + 2(-1)^n = 2^n + (-1)^n + (-1)^n: het spectrum {2,1,1}\{2, -1, -1\}, consistent met vraag 16.

20. De kolommen van WTW^{\mathsf T}: WTei=ei1W^{\mathsf T}e_i = e_{i-1} voor i1i \geq 1 en WTe0=ek1W^{\mathsf T}e_0 = e_{k-1}; herlabelen in de volgorde e0,e1,e_0, e_1, \dots is dit precies de begeleidende matrix van Xk1X^k - 1 (a0=1a_0 = 1, andere am=0a_m = 0). Vraag 2: χW=χWT=Xk1=μW\chi_{W} = \chi_{W^{\mathsf T}} = X^k - 1 = \mu_{W}. De wortels ωj\omega^j (j=0,,k1j = 0, \dots, k-1) zijn de kk verschillende kk-de eenheidswortels: WW is diagonaliseerbaar (vraag 3, of Oefening 3.8: Wk=IW^k = I). Eigenvectoren: Wfj=mωjmem+1=mωj(m1)em=ωjfjWf_j = \sum_m \omega^{-jm}e_{m+1} = \sum_{m'}\omega^{-j(m'-1)}e_{m'} = \omega^j f_j.

21. Circulanten zijn veeltermen in WW, en veeltermen in een vaste matrix commuten met elkaar. Elke fjf_j is een eigenvector van elke macht: Wmfj=ωjmfjW^m f_j = \omega^{jm}f_j, dus

Cfj=mcmωjmfj=c^(ωj)fj:Cf_j = \sum_m c_m\omega^{jm} f_j = \widehat c(\omega^j)\,f_j :

de basis (f0,,fk1)(f_0, \dots, f_{k-1}) (vrij: Vandermonde in de verschillende ωj\omega^{-j}, Oefening 2.11) diagonaliseert elke circulant tegelijk, met de gestelde eigenwaarden.

22. De determinant is het product van de eigenwaarden (diagonaliseer): detC=jc^(ωj)\det C = \prod_{j}\widehat c(\omega^j). Voor k=3k = 3, c0=ac_0 = a, c1=bc_1 = b, c2=cc_2 = c en ω=j=e2iπ/3\omega = j = \eu^{2\iu\pi/3}:

detC=(a+b+c)(a+bj+cj2)(a+bj2+cj4),\det C = (a + b + c)(a + bj + cj^2)(a + bj^2 + cj^4),

en j4=jj^4 = j: precies de factorisatie van Oefening 2.8.

23. M=12(W+W1)M = \frac12(W + W^{-1}) is een circulant (W1=Wk1W^{-1} = W^{k-1}), met eigenwaarden 12(ωj+ωj)=cos2πjk\frac12(\omega^j + \omega^{-j}) = \cos\frac{2\pi j}k op dezelfde basis fjf_j. Coördinaten: schrijf x(0)=jαjfjx^{(0)} = \sum_j \alpha_j f_j. De coördinaten van fjf_j sommeren tot mωjm\sum_m \omega^{-jm}, wat kk is voor j=0j = 0 en 00 anders (geometrische som met ratio ωj1\omega^{-j} \neq 1). Sommerend de coördinaten van x(0)x^{(0)}: mxm(0)=α0k\sum_m x^{(0)}_m = \alpha_0\,k, dus α0=1kmxm(0)\alpha_0 = \frac1k\sum_m x^{(0)}_m, het gemiddelde.

24. x(n)=Mnx(0)=jαjcosn(2πjk)fjx^{(n)} = M^nx^{(0)} = \sum_j \alpha_j\cos^n\bigl(\tfrac{2\pi j}k\bigr)f_j. Voor oneven kk, cos(2πj/k)<1\abs{\cos(2\pi j/k)} < 1 voor elke j0j \neq 0 (de hoek is nooit 00 of π\pi), dus neigen alle termen behalve j=0j = 0 naar 00: x(n)α0f0x^{(n)} \to \alpha_0 f_0, de constante vector gelijk aan het gemiddelde — middelen op een oneven ring egaliseert. Voor k=4k = 4 zijn de eigenwaarden 1,0,1,01, 0, -1, 0: de j=2j = 2-term α2(1)nf2\alpha_2(-1)^nf_2 met f2=(1,1,1,1)Tf_2 = (1, -1, 1, -1)^{\mathsf T} oscilleert eeuwig. De obstructie is het alternerende gemiddelde: de coördinaten van x(0)x^{(0)} met (1)m(-1)^m vermenigvuldigen en sommeren, dezelfde geometrische-som-berekening geeft m(1)mxm(0)=4α2\sum_m (-1)^mx^{(0)}_m = 4\alpha_2: het proces convergeert desda x0(0)x1(0)+x2(0)x3(0)=0x^{(0)}_0 - x^{(0)}_1 + x^{(0)}_2 - x^{(0)}_3 = 0, en convergeert dan naar het gemiddelde.

25. De begeleidende matrix zet een scalaire recurrentie van orde kk om in een eerste-orde vectorrecurrentie, zodat gesloten formules uitspraken over CnC^n worden — het thuisveld van reductie (vragen 1–5). Het kernontbindingslemma is pure veeltermalgebra (Bézout plus commutatie), dus splijt ze kerP(S)\ker P(S) hoewel S\mathcal{S} oneindigdimensionaal is (vragen 7–9). Dominante eigenwaarden regeren groei omdat elke andere bijdrage na normalisatie geometrisch verwaarloosbaar is — wat ook is waarom de Pell-fout vervalt met het kwadraat van de dominante wortel (vragen 11–14). Machten van de adjacentiematrix tellen wandelingen omdat matrixvermenigvuldiging over tussenliggende knopen sommeert, zodat spectra gesloten wandelingen tellen (vragen 15–19). Commuterende matrices delen een eigenbasis, en één Fourierbasis diagonaliseert dan heel de circulante algebra in één slag (vragen 20–24). Toppen: de fundamentele stelling van lineaire recurrenties (vraag 9); en voor niet-negatieve matrices is de reden dat dominante wortels zoals φ\varphi of 1+21 + \sqrt2 automatisch reëel, positief en enkelvoudig zijn de stelling van Perron–Frobenius, bewezen in het Bachelor jaar 3-volume.