Mathematics · Book 5 · Bachelor Year 3

Universitaire wiskunde — Bachelor jaar 3

Universitaire wiskunde — Bachelor jaar 3 · Bachelor Year 3

20Deelvariëteiten van ℝn

Bollen, tori, rotatiegroepen: de natuurlijke habitats van meetkunde en mechanica zijn geen vectorruimten maar gekromde verzamelingen die van dichtbij vlak lijken. Dit hoofdstuk geeft die zin een precieze betekenis — deelvariëteiten van Rn\R^n — en de calculus om erop te werken. De grondslag is de stelling van de inverse functie, hier bewezen via het Banach-vastpunt; al het overige is coördinatenwisseling: de vier equivalente beschrijvingen van een deelvariëteit (lokale rechttrekking, niveaucurven, grafieken, parametriseringen), raakruimten, en geoptimaliseerd onder nevenvoorwaarden via multiplicatoren van Lagrange — die, als afscheidende demonstratie, de spectraalstelling voor symmetrische matrices in drie regels meetkunde herbewezen. Het weekendprobleem bouwt de rotatiegroep SO(3)SO(3) en haar quaternionische dubbeldekking: algebra (Q8Q_8 van Hoofdstuk 1, volwassen) ontmoet meetkunde.

20.1 De stelling van de inverse functie

Stelling 20.1 (Stelling van de inverse functie)

Zij URnU \subseteq \R^n open, f ⁣:URnf \colon U \to \R^n van klasse C1\mathcal C^1, en aUa \in U met Df(a)Df(a) inverteerbaar. Dan bestaan open verzamelingen VaV \ni a, Wf(a)W \ni f(a) zodat f ⁣:VWf \colon V \to W een bijectie is met C1\mathcal C^1-inverse, en

D(f1)(y)=(Df(f1(y)))1(yW).D(f^{-1})(y) = \bigl(Df(f^{-1}(y))\bigr)^{-1} \qquad (y \in W).

Als ff van klasse Ck\mathcal C^k is, zo ook f1f^{-1}.

Bewijs. Normaliseer: vervang ff door xDf(a)1(f(a+x)f(a))x \mapsto Df(a)^{-1}\bigl(f(a + x) - f(a)\bigr), zodat we a=0a = 0, f(0)=0f(0) = 0, Df(0)=IDf(0) = I mogen aannemen (de algemene uitspraak volgt door samenstellen met de affiene bijecties). Schrijf f(x)=x+g(x)f(x) = x + g(x): Dg(0)=0Dg(0) = 0, en kies door continuïteit van DgDg een r>0r > 0 met Dg(x)12\vertiii{Dg(x)} \leq \frac12 op Bˉ(0,r)\bar B(0, r); de middelwaardeongelijkheid geeft g(x)g(x)12xx\norm{g(x) - g(x')} \leq \frac12 \norm{x - x'} daar.

Bijectiviteit op een omgeving. Voor yB(0,r2)y \in B(0, \frac r2) betekent f(x)=yf(x) = y oplossen het vinden van een vast punt van Φy(x)=yg(x)\Phi_y(x) = y - g(x); Φy\Phi_y beeldt Bˉ(0,r)\bar B(0, r) in zichzelf af (Φy(x)y+12xr\norm{\Phi_y(x)} \leq \norm y + \frac12\norm x \leq r) en is 12\frac12-Lipschitz: Banach (Stelling 7.4) levert een unieke oplossing x=φ(y)Bˉ(0,r)x = \varphi(y) \in \bar B(0,r). Bovendien is ff injectief op Bˉ(0,r)\bar B(0,r):

f(x)f(x)xxg(x)g(x)12xx.()\norm{f(x) - f(x')} \geq \norm{x - x'} - \norm{g(x) - g(x')} \geq \tfrac12\norm{x - x'} . \tag{$*$}

Stel W=B(0,r2)W = B(0, \frac r2) en V=f1(W)B(0,r)V = f^{-1}(W)\cap B(0, r): open (continuïteit), met f ⁣:VWf \colon V \to W bijectief.

Continuïteit en differentieerbaarheid van de inverse. ()(*) zegt dat φ=f1\varphi = f^{-1} 22-Lipschitz is. Fixeer y0=f(x0)Wy_0 = f(x_0) \in W; inverteerbaarheid van A=Df(x0)A = Df(x_0) (haar afstand tot II is 12\leq \frac12: Neumann, Propositie 8.4) en differentieerbaarheid van ff geven, voor y=f(x)y = f(x) nabij y0y_0:

φ(y)φ(y0)A1(yy0)=A1(f(x)f(x0)A(xx0))=A1o(xx0)=o(yy0),\varphi(y) - \varphi(y_0) - A^{-1}(y - y_0) = -A^{-1}\bigl(f(x) - f(x_0) - A(x - x_0)\bigr) = -A^{-1}\,o(\norm{x - x_0}) = o(\norm{y - y_0}),

met ()(*) om xx02yy0\norm{x - x_0} \leq 2\norm{y - y_0} om te zetten: φ\varphi is differentieerbaar in y0y_0 met de inverse differentiaal. Continuïteit van yDφ(y)=Df(φ(y))1y \mapsto D\varphi(y) = Df(\varphi(y))^{-1}: samenstelling van continue afbeeldingen (inversie is continu, Propositie 8.4): φC1\varphi \in \mathcal C^1; bootstrappen van dezelfde formule geeft Ck\mathcal C^k.

Stelling 20.2 (Stelling van de impliciete functie)

Zij F ⁣:URp×RqRqF \colon U \subseteq \R^p\times\R^q \to \R^q van klasse C1\mathcal C^1 nabij (a,b)(a, b), F(a,b)=0F(a,b) = 0, en veronderstel dat de partiële differentiaal DyF(a,b)L(Rq)D_yF(a,b) \in \mathcal L(\R^q) inverteerbaar is. Dan bestaan omgevingen AaA \ni a, BbB \ni b en een C1\mathcal C^1-afbeelding ψ ⁣:AB\psi \colon A \to B met

{(x,y)A×B:F(x,y)=0}={(x,ψ(x)):xA},\bigl\{(x, y)\in A\times B : F(x,y) = 0\bigr\} = \{(x, \psi(x)) : x \in A\},

en Dψ(x)=DyF(x,ψ(x))1DxF(x,ψ(x))D\psi(x) = -D_yF(x, \psi(x))^{-1}\,D_xF(x, \psi(x)).

Bewijs. Pas Stelling 20.1 toe op Θ(x,y)=(x,F(x,y))\Theta(x, y) = (x, F(x,y)): haar differentiaal in (a,b)(a,b), blok-driehoekig met inverteerbare diagonale blokken II en DyFD_yF, is inverteerbaar. De lokale inverse heeft de vorm Θ1(x,z)=(x,h(x,z))\Theta^{-1}(x, z) = (x, h(x, z)); stel ψ(x)=h(x,0)\psi(x) = h(x, 0): dan F(x,y)=0F(x, y) = 0 desda Θ(x,y)=(x,0)\Theta(x,y) = (x, 0) desda y=ψ(x)y = \psi(x), lokaal. De formule: differentieer F(x,ψ(x))=0F(x, \psi(x)) = 0 met de kettingregel.

20.2 Deelvariëteiten: vier definities

Stelling 20.3 (Equivalente karakteriseringen)

Zij MRnM \subseteq \R^n, d{0,,n}d \in \{0, \dots, n\}, en k1k \geq 1. De volgende zijn equivalent, voor elk punt aMa \in M (en MM is een dd-dimensionale deelvariëteit van klasse Ck\mathcal C^k als zij in elke aMa \in M gelden):

  1. (Rechttrekking) Er bestaat een Ck\mathcal C^k-diffeomorfisme Φ\Phi van een open Ωa\Omega \ni a op een open ΩRn\Omega' \subseteq \R^n met

    Φ(MΩ)=Ω(Rd×{0}).\Phi(M\cap\Omega) = \Omega' \cap \bigl(\R^d\times\{0\}\bigr).
  2. (Niveaucurve) Er bestaat een Ck\mathcal C^k-submersie F ⁣:ΩRndF \colon \Omega \to \R^{n-d} (d.w.z. DF(x)DF(x) surjectief) op een open Ωa\Omega \ni a met MΩ=F1(0)M\cap\Omega = F^{-1}(0).
  3. (Grafiek) Op permutatie van coördinaten na is MM lokaal de grafiek van een Ck\mathcal C^k-afbeelding ψ ⁣:ARdRnd\psi \colon A \subseteq \R^d \to \R^{n-d}.
  4. (Parametrisering) Er bestaat een Ck\mathcal C^k-immersie φ ⁣:ARdRn\varphi \colon A \subseteq \R^d \to \R^n (Dφ(u)D\varphi(u) injectief) met AA open, φ\varphi een homeomorfisme van AA op MΩM \cap \Omega voor enige open Ωa\Omega \ni a.

Bewijs. (1)\Rightarrow(2): F=(Φd+1,,Φn)F = (\Phi_{d+1}, \dots, \Phi_n) (laatste coördinaten van Φ\Phi): een submersie (DΦD\Phi inverteerbaar). (2)\Rightarrow(3): DF(a)DF(a) surjectief: enige q×qq \times q-minor van de Jacobiaan is inverteerbaar (q=ndq = n - d); na permutatie van coördinaten is DyF(a)D_yF(a) inverteerbaar, en de stelling van de impliciete functie (Stelling 20.2) drukt MM lokaal uit als grafiek y=ψ(x)y = \psi(x). (3)\Rightarrow(4): φ(x)=(x,ψ(x))\varphi(x) = (x, \psi(x)): een immersie (differentiaal (IDψ)\bigl(\begin{smallmatrix}I\\ D\psi\end{smallmatrix}\bigr) injectief), een homeomorfisme op de grafiek (inverse: de projectie, continu). (4)\Rightarrow(1): zij φ(u0)=a\varphi(u_0) = a; vul imDφ(u0)\operatorname{im}D\varphi(u_0) aan met een supplement EE (dimE=nd\dim E = n - d) en definieer Θ(u,v)=φ(u)+v\Theta(u, v) = \varphi(u) + v op A×EA\times E: DΘ(u0,0)D\Theta(u_0, 0) is bijectief (beeld bevat imDφ(u0)\operatorname{im}D\varphi(u_0) en EE), dus is Θ\Theta een lokaal diffeomorfisme (Stelling 20.1); haar inverse Φ\Phi trekt recht: nabij aa zijn punten van MM precies de φ(u)=Θ(u,0)\varphi(u) = \Theta(u, 0) — hiervoor garandeert de homeomorfisme-hypothese in (4) dat MΩM\cap\Omega, voor kleine Ω\Omega, geen andere bladen bevat (φ(u)+v=mM\varphi(u') + v = m \in M dicht bij aa met kleine v0v \neq 0 moet worden uitgesloten: m=φ(u)m = \varphi(u'') voor enige uu'' nabij u0u_0 door de homeomorfisme-eigenschap, en lokale injectiviteit van Θ\Theta dwingt v=0v = 0). Dan Φ(MΩ)=(A×{0})Φ(Ω)\Phi(M\cap\Omega) = (A\times\{0\}) \cap\Phi(\Omega) na inkrimpen.

Voorbeeld 20.4

De sfeer Sn1={x2=1}S^{n-1} = \{\norm x^2 = 1\}: niveaucurve van de submersie F(x)=x221F(x) = \norm x_2^2 - 1 op Rn{0}\R^n\setminus\{0\} (DF(x)=2xT0DF(x) = 2x^{\mathsf T} \ne 0): een C\mathcal C^\infty-deelvariëteit van dimensie n1n - 1. De torus in R3\R^3: niveaucurve van (x2+y2R)2+z2r2\bigl(\sqrt{x^2 + y^2} - R\bigr)^2 + z^2 - r^2 (0<r<R0 < r < R). De kegel {x2+y2=z2}\{x^2 + y^2 = z^2\} is geen deelvariëteit in 00 (Oefening 20.1). Matrixgroepen: SLnSL_n en OnO_n zijn deelvariëteiten van Mn(R)M_n(\R) (Oefeningen 20.5 en 20.6) — het startpunt van de Lie-theorie.

20.3 Raakruimten

Definitie 20.5

Zij MM een dd-deelvariëteit en aMa \in M. De raakruimte TaMT_aM is de verzameling snelheidsvectoren γ(0)\gamma'(0) van C1\mathcal C^1-krommen γ ⁣:(ε,ε)M\gamma \colon \intoo{-\varepsilon}\varepsilon \to M met γ(0)=a\gamma(0) = a.

Propositie 20.6

TaMT_aM is een dd-dimensionale vectordeelruimte van Rn\R^n, en:

  1. als M=F1(0)M = F^{-1}(0) lokaal met FF een submersie: TaM=kerDF(a)T_aM = \ker DF(a);
  2. als MM geparametriseerd wordt door de immersie φ\varphi (φ(u0)=a\varphi(u_0) = a): TaM=imDφ(u0)T_aM = \operatorname{im}D\varphi(u_0).

Bewijs. Krommen in MM voldoen aan F(γ(t))=0F(\gamma(t)) = 0; de kettingregel in 00 geeft DF(a)γ(0)=0DF(a)\gamma'(0) = 0: TaMkerDF(a)T_aM \subseteq \ker DF(a). Omgekeerd transporteert rechttrekking (Stelling 20.3(1)) lijnen van Rd×{0}\R^d\times\{0\} naar krommen van MM: elke vector van een dd-dimensionale deelruimte wordt gerealiseerd; dimensies vergelijken (dimkerDF(a)=n(nd)=d\dim\ker DF(a) = n - (n - d) = d) dwingt gelijkheid in (1), en hetzelfde transportargument geeft (2) (Dφ(u0)D\varphi(u_0) toegepast op rechte lijnen in AA; dimensies opnieuw).

Stelling 20.7 (Multiplicatoren van Lagrange)

Zij M=F1(0)M = F^{-1}(0) met F=(F1,,Fq) ⁣:ΩRqF = (F_1, \dots, F_q) \colon \Omega \to \R^q een C1\mathcal C^1-submersie, en zij f ⁣:ΩRf \colon \Omega \to \R van klasse C1\mathcal C^1. Als de restrictie fMf\restriction_M een lokaal extremum heeft in aMa \in M, dan bestaan unieke reële λ1,,λq\lambda_1, \dots, \lambda_q (multiplicatoren van Lagrange) met

f(a)=λ1F1(a)++λqFq(a).\nabla f(a) = \lambda_1\nabla F_1(a) + \dots + \lambda_q\nabla F_q(a) .

Bewijs. Voor elke kromme γ\gamma in MM door aa: tf(γ(t))t \mapsto f(\gamma(t)) heeft een lokaal extremum in 00, dus 0= ⁣d ⁣dtf(γ(t))0=f(a),γ(0)0 = \frac{\dd}{\dd t}f(\gamma(t))\big|_0 = \langle\nabla f(a), \gamma'(0)\rangle: f(a)TaM=kerDF(a)\nabla f(a) \perp T_aM = \ker DF(a) (Propositie 20.6). Nu is kerDF(a)=imDF(a)T\ker DF(a)^\perp = \operatorname{im}DF(a)^{\mathsf T} (de rang/orthogonaliteitsidentiteit van Jaar 2, of Oefening 13.8 in eindige dimensie: (kerT)=imT(\ker T)^\perp = \operatorname{im}T^*), die wordt opgespannen door de gradiënten Fi(a)\nabla F_i(a) — onafhankelijk, omdat DF(a)DF(a) surjectief is: de multiplicatoren bestaan en zijn uniek.

Voorbeeld 20.8 (De spectraalstelling, meetkundig)

Zij AA een reële symmetrische n×nn\times n-matrix en maximaliseer f(x)=Ax,xf(x) = \langle Ax, x\rangle op de sfeer Sn1S^{n-1} (compact: het maximum wordt bereikt, in enige v1v_1). Lagrange met F(x)=x21F(x) = \norm x^2 - 1: f=2Ax\nabla f = 2Ax en F=2x\nabla F = 2x geven Av1=λ1v1Av_1 = \lambda_1v_1 — een eigenvector, met λ1=maxSn1Ax,x\lambda_1 = \max_{S^{n-1}}\langle Ax, x\rangle. Restricteer AA tot v1v_1^\perp (invariant: Av,v1=v,Av1=λ1v,v1=0\langle Av, v_1\rangle = \langle v, Av_1\rangle = \lambda_1\langle v, v_1\rangle = 0) en itereer: een orthonormal basis van eigenvectoren. De spectraalstelling van Jaar 2, herbewezen door pure optimalisatie — en de oneindigdimensionale schaduw van hetzelfde argument bewees Lemma 15.6.

Methode 20.9

Om te bewijzen dat een verzameling een deelvariëteit is: toon haar lokaal als F1(0)F^{-1}(0) met DFDF surjectief op de verzameling (de meest voorkomende route), of als grafiek. Om dimensie en raakruimte te berekenen: d=nqd = n - q en Ta=kerDF(a)T_a = \ker DF(a). Om erop te optimaliseren: Lagrange — controleer altijd eerst compactheid (of coërciviteit), zodat een extremum bestaat waarop de stelling kan worden toegepast, en onthoud dat de multiplicatorvergelijking slechts noodzakelijk is: verzamel alle kritieke punten, vergelijk dan waarden. Voor matrixgroepen, differentieer krommen in de identiteit om raakruimten te identificeren.

20.4 Oefeningen

Oefening 20.1

(a) Verifieer dat de volgende C\mathcal C^\infty-deelvariëteiten zijn en geef hun dimensies: Sn1S^{n-1}; de hyperboloïde {x2+y2z2=1}\{x^2 + y^2 - z^2 = 1\}; de torus van Voorbeeld 20.4. (b) Toon aan dat de kegel C={x2+y2=z2}R3C = \{x^2 + y^2 = z^2\} \subseteq \R^3 geen 22-deelvariëteit is in 00: bepaal het aantal samenhangscomponenten van (C{0})B(0,ε)\bigl(C\setminus\{0\}\bigr)\cap B(0,\varepsilon), en vergelijk met het aantal dat een rechttrekking (Stelling 20.3(1)) zou dwingen voor een vlak minus een punt.

Oplossing

Oplossing van Oefening 20.1.

(a) Elk is F1(0)F^{-1}(0) voor een submersie: x21\norm x^2 - 1 op Rn{0}\R^n\setminus\{0\} (gradiënt 2x02x \neq 0): dimensie n1n-1; x2+y2z21x^2 + y^2 - z^2 - 1 (gradiënt (2x,2y,2z)0(2x, 2y, -2z) \neq 0 op de hyperboloïde, waar x2+y2=1+z2>0x^2 + y^2 = 1 + z^2 > 0): dimensie 22; de torusfunctie G=(ρR)2+z2r2G = (\rho - R)^2 + z^2 - r^2, ρ=x2+y2\rho = \sqrt{x^2+y^2}, is C\mathcal C^\infty nabij de torus (ρRr>0\rho \geq R - r > 0 daar) met G0\nabla G \neq 0 (haar zz-component is 2z2z, en waar z=0z = 0 is de radiale component 2(ρR)02(\rho - R)\ne0 want ρR=r\abs{\rho - R} = r): dimensie 22.

(b) Voor kleine ε\varepsilon heeft (C{0})B(0,ε)(C\setminus\{0\})\cap B(0,\varepsilon) precies 22 samenhangscomponenten (bovenste en onderste doorboorde nappen, elk padsamenhangend: verbind via cirkels en stralen). Als CC een 22-deelvariëteit was in 00, zou een rechttrekking een homeomorfisme geven van CΩC\cap\Omega op een open stuk van een vlak dat 00 naar een punt pp stuurt; kleine doorboorde vlakomgevingen van pp hebben één component, en homeomorfismen bewaren het aantal componenten van doorboorde omgevingen: tegenspraak.

Oefening 20.2

Bereken de raakruimten: (a) TaSn1T_aS^{n-1} voor willekeurige aa (antwoord: aa^\perp); (b) het raakvlak aan de torus van Voorbeeld 20.4 in een willekeurig punt van de buitenste evenaar {z=0, x2+y2=(R+r)2}\{z = 0,\ x^2 + y^2 = (R + r)^2\}; (c) de raaklijn aan de helix φ(t)=(cost,sint,t)\varphi(t) = (\cos t, \sin t, t) in φ(t0)\varphi(t_0), met controle van Propositie 20.6(2).

Oplossing

Oplossing van Oefening 20.2.

(a) TaSn1=ker(2aT)=aT_aS^{n-1} = \ker\bigl(2a^{\mathsf T}\bigr) = a^\perp. (b) In p=((R+r)cosθ,(R+r)sinθ,0)p = ((R+r)\cos\theta, (R+r)\sin\theta, 0): G=(2rcosθ,2rsinθ,0)\nabla G = (2r\cos\theta, 2r\sin\theta, 0), dus is het raakvlak Vect((sinθ,cosθ,0), (0,0,1))\operatorname{Vect}\bigl((-\sin\theta, \cos\theta, 0),\ (0, 0, 1)\bigr): het verticale vlak raak aan de buitenste evenaar. (c) De helix is een ingebedde kromme met φ(t0)=(sint0,cost0,1)0\varphi'(t_0) = (-\sin t_0, \cos t_0, 1) \neq 0: de raaklijn in φ(t0)\varphi(t_0) is φ(t0)+Rφ(t0)\varphi(t_0) + \R\,\varphi'(t_0), zoals Propositie 20.6(2) voorschrijft.

Oefening 20.3 ★★

Zij f(x,y)=(x2y2, 2xy)f(x, y) = (x^2 - y^2,\ 2xy) (d.w.z. zz2z \mapsto z^2). (a) In welke punten past Stelling 20.1 toe? (b) Toon aan dat ff lokaal maar niet globaal inverteerbaar is op R2{0}\R^2\setminus\{0\}, en geef expliciet de twee lokale inversen gedefinieerd op een omgeving van (1,0)(1, 0) (de twee vierkantsworteltakken). (c) Zelfde discussie voor de poolcoördinatenafbeelding (r,θ)(rcosθ,rsinθ)(r, \theta) \mapsto (r\cos\theta, r\sin\theta).

Oplossing

Oplossing van Oefening 20.3.

(a) Df(x,y)=(2x2y2y2x)Df(x,y) = \bigl(\begin{smallmatrix}2x & -2y\\ 2y & 2x\end{smallmatrix}\bigr), det=4(x2+y2)\det = 4(x^2 + y^2): de stelling past toe in elk punt behalve de oorsprong. (b) f(z)=f(z)f(-z) = f(z): nooit injectief op een verzameling symmetrisch om 00; op R2{0}\R^2\setminus\{0\} is zij overal een lokaal diffeomorfisme maar globaal 22-op-11. Nabij (1,0)=f(±(1,0))(1, 0) = f(\pm(1, 0)) zijn de twee inversen de twee vierkantsworteltakken: in complexe notatie w±ww \mapsto \pm\sqrt w (principale tak), d.w.z.

(u,v)±(u+u2+v22, v2(u+u2+v2)/2).(u, v) \longmapsto \pm\Bigl(\sqrt{\tfrac{u + \sqrt{u^2+v^2}}{2}},\ \frac{v}{2\sqrt{(u + \sqrt{u^2+v^2})/2}}\Bigr).

(c) Jacobiaan r>0r > 0: lokaal diffeomorfisme op (0,)×R\intoo0\infty\times\R, maar θθ+2π\theta \mapsto \theta + 2\pi geeft hetzelfde punt: lokaal inverteerbaar (hoek bepaald tot op 2π2\pi op een halfvlak), nooit globaal.

Oefening 20.4 ★★

(Folium) Zij F(x,y)=x3+y33xyF(x,y) = x^3 + y^3 - 3xy en C=F1(0)\mathcal C = F^{-1}(0). (a) Toon aan dat nabij elk punt van C\mathcal C behalve de oorsprong C\mathcal C een 11-deelvariëteit is, lokaal een grafiek in xx of in yy (welke, waar?). (b) Bereken de raaklijn in (32,32)(\frac32, \frac32). (c) Wat gebeurt er in de oorsprong? (Twee takken kruisen: toon twee C1\mathcal C^1-krommen in C\mathcal C door 00 met onafhankelijke snelheden, en concludeer dat geen rechttrekking bestaat.)

Oplossing

Oplossing van Oefening 20.4.

(a) F=3(x2y, y2x)\nabla F = 3(x^2 - y,\ y^2 - x) verdwijnt desda y=x2y = x^2 en x=y2x = y^2, d.w.z. x4=xx^4 = x: in (0,0)(0,0) en (1,1)(1,1); alleen (0,0)(0,0) ligt op C\mathcal C (F(1,1)=1F(1,1) = -1). Dus is FF op C{0}\mathcal C\setminus\{0\} een submersie: een 11-deelvariëteit, lokaal een grafiek y=ψ(x)y = \psi(x) waar Fy=3(y2x)0F_y = 3(y^2 - x) \neq 0 en x=χ(y)x = \chi(y) waar Fx=3(x2y)0F_x = 3(x^2 - y) \neq 0 (minstens één geldt buiten de oorsprong).

(b) In (32,32)(\frac32, \frac32): F=3(9432)(1,1)=94(1,1)\nabla F = 3(\frac94 - \frac32)(1, 1) = \frac94(1,1): raaklijn x+y=3x + y = 3.

(c) De rationale parametrisering x=3t1+t3x = \frac{3t}{1 + t^3}, y=3t21+t3y = \frac{3t^2}{1+t^3} gaat door 00 in t=0t = 0 met snelheid (3,0)(3, 0); xyx \leftrightarrow y omwisselen (de curve is symmetrisch, of herparametriseer met 1/t1/t) geeft een tweede C1\mathcal C^1-kromme door 00 met snelheid (0,3)(0, 3). Twee onafhankelijke raakrichtingen zijn onmogelijk voor een 11-deelvariëteit (haar raakruimte is een lijn, Propositie 20.6): C\mathcal C is geen deelvariëteit in de oorsprong — een transversale zelfkruising.

Oefening 20.5 ★★

Zij F(M)=MTMF(M) = M^{\mathsf T}M van Mn(R)M_n(\R) naar de ruimte SnS_n van symmetrische matrices. (a) Toon aan DF(M)(H)=MTH+HTMDF(M)(H) = M^{\mathsf T}H + H^{\mathsf T}M en dat DF(M)DF(M) surjectief is op SnS_n in elke MOnM \in O_n (gegeven SSnS \in S_n, probeer H=12MSH = \frac12MS). (b) Concludeer dat On=F1(I)O_n = F^{-1}(I) een compacte C\mathcal C^\infty-deelvariëteit is van dimensie n(n1)2\frac{n(n-1)}2, met TIOn={H:HT=H}T_IO_n = \{H : H^{\mathsf T} = -H\}, de antisymmetrische matrices. (c) Toon aan dat etHOn\eu^{tH} \in O_n voor elke antisymmetrische HH: de raakrichtingen integreren tot krommen in de groep.

Oplossing

Oplossing van Oefening 20.5.

(a) F(M+H)=MTM+MTH+HTM+HTHF(M + H) = M^{\mathsf T}M + M^{\mathsf T}H + H^{\mathsf T}M + H^{\mathsf T}H: DF(M)(H)=MTH+HTMDF(M)(H) = M^{\mathsf T}H + H^{\mathsf T}M. Voor MOnM \in O_n en SS symmetrisch geeft H=12MSH = \frac12MS DF(M)(H)=12(S+ST)=SDF(M)(H) = \frac12(S + S^{\mathsf T}) = S: surjectief op SnS_n.

(b) On=F1(I)O_n = F^{-1}(I) met FF een submersie (op SnS_n, van dimensie n(n+1)2\frac{n(n+1)}2) in elk van haar punten: een deelvariëteit van dimensie n2n(n+1)2=n(n1)2n^2 - \frac{n(n+1)}2 = \frac{n(n-1)}2. Compact: gesloten (FF continu) en begrensd (kolommen zijn eenheidsvectoren). Raakruimte in II: kerDF(I)={H:H+HT=0}\ker DF(I) = \{H : H + H^{\mathsf T} = 0\}.

(c) (etH)TetH=etHTetH=etHetH=I\bigl(\eu^{tH}\bigr)^{\mathsf T}\eu^{tH} = \eu^{tH^{\mathsf T}}\eu^{tH} = \eu^{-tH}\eu^{tH} = I (transponeren gaat door de reeks; de exponentialen van commuterende matrices vermenigvuldigen, Stelling 19.8).

Oefening 20.6 ★★

(a) Toon aan dat det ⁣:Mn(R)R\det \colon M_n(\R) \to \R differentiaal Ddet(M)(H)=tr(com(M)TH)D\det(M)(H) = \operatorname{tr}\bigl(\operatorname{com}(M) ^{\mathsf T}H\bigr) heeft, niet-nul in elke MSLnM \in SL_n. (b) Concludeer dat SLn(R)SL_n(\R) een deelvariëteit is van dimensie n21n^2 - 1 met TISLn={H:trH=0}T_ISL_n = \{H : \operatorname{tr}H = 0\}. (c) Is GLn(R)GL_n(\R) een deelvariëteit? Van welke dimensie?

Oplossing

Oplossing van Oefening 20.6.

(a) det(M+H)=detMdet(I+M1H)=detM(1+tr(M1H)+O(H2))\det(M + H) = \det M\,\det(I + M^{-1}H) = \det M\bigl(1 + \operatorname{tr}(M^{-1}H) + O(\norm H^2)\bigr) voor inverteerbare MM (ontwikkeling van det\det nabij II: de lineaire term van (1+λi)\prod(1 + \lambda_i)); met detMM1=com(M)T\det M\cdot M^{-1} = \operatorname{com}(M)^{\mathsf T}: Ddet(M)(H)=tr(com(M)TH)D\det(M)(H) = \operatorname{tr}\bigl(\operatorname{com} (M)^{\mathsf T}H\bigr), en de formule strekt zich tot alle MM uit door densiteit en continuïteit. Op SLnSL_n, detM=1\det M = 1: Ddet(M)0D\det(M) \ne 0 (haar waarde op H=MH = M is tr(I)=ndetM=n\operatorname{tr}(I)\cdot\dots = n\det M = n).

(b) SLn=det1(1)SL_n = \det^{-1}(1) met det\det daar een submersie (waarden in R\R): dimensie n21n^2 - 1; TISLn=kerDdet(I)={H:trH=0}T_ISL_n = \ker D\det(I) = \{H : \operatorname{tr}H = 0\}.

(c) GLnGL_n is een open deelverzameling van Mn(R)M_n(\R) (Oefening 6.8): een deelvariëteit van volle dimensie n2n^2 (rechttrekking: de identiteitskaart).

Oefening 20.7 ★★

Via multiplicatoren van Lagrange: (a) vind de extrema van f(x,y)=xyf(x,y) = xy op de cirkel x2+y2=1x^2 + y^2 = 1; (b) toon aan dat onder alle kansvectoren (p1,,pn)(p_1, \dots, p_n) (positief, sommerend tot 11) de entropie pilnpi-\sum p_i\ln p_i precies in de uniforme verdeling wordt gemaximaliseerd; (c) vind het punt van de ellips {x2/4+y2=1}\{x^2/4 + y^2 = 1\} dichtst bij (1,0)(1, 0), en controleer de multiplicatorvergelijking meetkundig (normaal-uitlijning).

Oplossing

Oplossing van Oefening 20.7.

(a) (y,x)=λ(2x,2y)(y, x) = \lambda(2x, 2y) en x2+y2=1x^2 + y^2 = 1: y=2λxy = 2\lambda x, x=2λyx = 2\lambda y geven x2=y2=12x^2 = y^2 = \frac12. Waarden van xyxy: ±12\pm\frac12: maximum 12\frac12 in ±12(1,1)\pm\frac1{\sqrt2}(1,1), minimum 12-\frac12 in ±12(1,1)\pm\frac1{\sqrt2}(1,-1) (de constrainteverzameling is compact: extrema bestaan).

(b) Op het inwendige van de simplex (pi>0p_i > 0), Lagrange voor H(p)=pilnpiH(p) = -\sum p_i\ln p_i met constrainte pi=1\sum p_i = 1: lnpi1=λ-\ln p_i - 1 = \lambda voor alle ii: alle pip_i gelijk, pi=1np_i = \frac1n, met H=lnnH = \ln n. Het maximum over de compacte simplex wordt bereikt; als het op de rand werd bereikt (enige pi=0p_i = 0), leeft de verdeling op n1\leq n - 1 punten en door inductie Hln(n1)<lnnH \leq \ln(n-1) < \ln n: het inwendige kritieke punt is het globale maximum — uniforme onwetendheid maximaliseert entropie.

(c) Minimaliseer (x1)2+y2(x-1)^2 + y^2 op de compacte ellips: (2(x1),2y)=λ(x2,2y)(2(x{-}1), 2y) = \lambda(\frac x2, 2y). Als y0y \neq 0: λ=1\lambda = 1, dan 2(x1)=x22(x - 1) = \frac x2 geeft x=43x = \frac43, y2=149=59y^2 = 1 - \frac49 = \frac59: afstand2^2 =19+59=23= \frac19 + \frac59 = \frac23. Als y=0y = 0: x=±2x = \pm2, afstanden 11 en 33. Dichtstbijzijnde punten: (43,±53)\bigl(\frac43, \pm\frac{\sqrt5}3\bigr), op afstand 2/3<1\sqrt{2/3} < 1. De multiplicatorvergelijking zegt dat het segment van (1,0)(1,0) naar het dichtstbijzijnde punt parallel is aan \nabla(ellips): het snijdt de ellips orthogonaal, zoals de meetkunde eist.

Oefening 20.8 ★★★

Schrijf Voorbeeld 20.8 volledig uit: bewijs door inductie dat een reële symmetrische matrix een orthonormal basis van eigenvectoren toelaat, met λ1λn\lambda_1 \geq \dots \geq \lambda_n de opeenvolgende geconstrueerde maxima van het Rayleigh-quotiënt. Leid dan de formules van Courant–Fischer van Oefening 15.8 in eindige dimensie rechtstreeks uit deze constructie af.

Oplossing

Oplossing van Oefening 20.8.

Inductie op nn; n=1n = 1 triviaal. De Rayleigh-functie f(x)=Ax,xf(x) = \langle Ax, x\rangle bereikt haar maximum λ1\lambda_1 op de compacte Sn1S^{n-1} in enige v1v_1; Lagrange (Stelling 20.7, sfeer als niveaucurve) geeft 2Av1=2λv12Av_1 = 2\lambda v_1, en λ=Av1,v1=λ1\lambda = \langle Av_1, v_1\rangle = \lambda_1. Het hypervlak v1v_1^\perp is AA-invariant (symmetrie: Av,v1=v,Av1=0\langle Av, v_1\rangle = \langle v, Av_1\rangle = 0); de restrictie is symmetrisch, en inductie levert een orthonormale eigenbasis v2,,vnv_2, \dots, v_n van v1v_1^\perp met eigenwaarden λ2λn\lambda_2 \geq \dots \geq \lambda_n, elk het maximum van ff op de sfeer van het resterende orthocomplement. Courant–Fischer volgt precies zoals in Oefening 15.8: ontwikkel x=civix = \sum c_iv_i; op een kk-dimensionale testruimte snijd met Vect(vk,,vn)\operatorname{Vect}(v_k, \dots, v_n) (dimensietelling in Rn\R^n) om minλk\min \leq \lambda_k te krijgen, en Vect(v1,,vk)\operatorname{Vect}(v_1, \dots, v_k) bereikt min=λk\min = \lambda_k.

Oefening 20.9 ★★★

(Ongelijkheid van Hadamard) Voor MGLn(R)M \in GL_n(\R) met kolommen c1,,cnc_1, \dots, c_n:

detM    i=1nci2,\abs{\det M} \;\leq\; \prod_{i=1}^n\norm{c_i}_2 ,

met gelijkheid desda de kolommen orthogonaal zijn. (Reduceer tot kolommen van norm 11 door schaling; maximaliseer det\det op het compacte product van bollen (Sn1)n(S^{n-1})^n; in een maximizer geeft Lagrange in elke kolom apart cidet=λici\nabla_{c_i}\det = \lambda_ic_i, en cidet\nabla_{c_i}\det is de ii-de kolom van com(M)\operatorname{com}(M): leid af dat MTMM^{\mathsf T}M diagonaal is, dus =I= I, dus det=±1\det = \pm1.) Meetkundige lezing: het volume van een parallellepipedum is hoogstens het product van zijn randlengten.

Oplossing

Oplossing van Oefening 20.9.

Schalen van elke kolom naar eenheidsnorm deelt det\abs{\det} door ci\prod\norm{c_i}: het volstaat detM1\abs{\det M} \leq 1 te bewijzen wanneer alle kolommen eenheid zijn, met gelijkheid desda MOnM \in O_n. De functie det\det is continu op de compacte (Sn1)n(S^{n-1})^n: zij bereikt een maximum mdetI=1>0m \geq \det I = 1 > 0 in enige MM. Alle kolommen behalve de ii-de fixeren, is det\det lineair in cic_i met gradiënt de ii-de kolom van com(M)\operatorname{com}(M); Lagrange op de ii-de sfeer: com(M)i=λici\operatorname{com}(M)_{\cdot i} = \lambda_i c_i. De identiteit MTcom(M)=det(M)IM^{\mathsf T}\operatorname{com}(M) = \det(M)\,I leest cj,com(M)i=det(M)δij\langle c_j, \operatorname{com}(M)_{\cdot i}\rangle = \det(M)\,\delta_{ij}, d.w.z. λicj,ci=det(M)δij\lambda_i\langle c_j, c_i\rangle = \det(M)\delta_{ij}; met j=ij = i: λi=detM=m0\lambda_i = \det M = m \neq 0, en dan jij \neq i geeft ci,cj=0\langle c_i, c_j\rangle = 0: de kolommen zijn orthonormaal, MOnM \in O_n, m=detM=1m = \abs{\det M} = 1. Dus detci\abs{\det} \leq \prod\norm{c_i} altijd, met gelijkheid precies voor orthogonale kolommen (herschaal terug): het volume van een parallellepipedum is, bij gegeven randlengten, het grootst wanneer de randen loodrecht staan.

Oefening 20.10 ★★

Nabij welke van haar punten is de cirkel S1S^1 een grafiek y=ψ(x)y = \psi(x)? Een grafiek x=χ(y)x = \chi(y)? Verifieer de grafiekkarakterisering (Stelling 20.3(3)) expliciet in (1,0)(1, 0), en leg in één zin uit waarom enige coördinatenpermutatie altijd volstaat maar geen enkele altijd werkt.

Oplossing

Oplossing van Oefening 20.10.

y=±1x2y = \pm\sqrt{1 - x^2} werkt nabij elk punt met y0y \neq 0; x=±1y2x = \pm\sqrt{1 - y^2} nabij elk punt met x0x \neq 0; in (1,0)(1, 0): de grafiek x=1y2x = \sqrt{1 - y^2} over y(1,1)y \in \intoo{-1}1, wat Stelling 20.3(3) is met de coördinaten omgewisseld. Enige permutatie werkt altijd omdat de raaklijn, eendimensionaal, niet tegelijk verticaal en horizontaal kan zijn — maar wel een van beide, dus dient geen vaste keuze van “afhankelijke” coördinaat in elk punt.

Oefening 20.11 ★★

(De orthogonale groep als deelvariëteit, kwantitatief) (a) Toon aan dat On={M:MTM=I}O_n = \{M : M^{\mathsf T}M = I\} compact is: begrensd (elke kolom is een eenheidsvector, dus Mn\norm M \leq \sqrt n voor de Euclidische matrixnorm) en gesloten. (b) Toon aan dat haar raakruimte in II de ruimte van antisymmetrische matrices is, van dimensie n(n1)2\frac{n(n-1)}2, en in een algemene AOnA \in O_n: TAOn={AK:KT=K}T_AO_n = \{AK : K^{\mathsf T} = -K\}. (c) Leid af dat de afbeelding tAexp(tK)t \mapsto A\exp(tK) voor elke antisymmetrische KK een kromme is in OnO_n door AA met snelheid AKAK (verifieer exp(tK)On\exp(tK) \in O_n met exp(X)T=exp(XT)\exp(X)^{\mathsf T} = \exp(X^{\mathsf T}) en exp(X)exp(X)=I\exp(-X)\exp(X) = I): elke raakvector wordt gerealiseerd door een expliciete kromme, zonder de stelling van de impliciete functie.

Oplossing

Oplossing van Oefening 20.11.

(a) De definiërende afbeelding F(M)=MTMIF(M) = M^{\mathsf T}M - I is continu: On=F1(0)O_n = F^{-1}(0) is gesloten; kolommen van een orthogonale matrix zijn eenheidsvectoren, dus is de Euclidische (Frobenius-)norm precies n\sqrt n: begrensd. Compact door Heine–Borel in Mn(R)Rn2M_n(\R) \cong \R^{n^2}.

(b) OnO_n is de niveaucurve F=0F = 0 bestudeerd in het hoofdstuk: DF(A)H=ATH+HTADF(A)H = A^{\mathsf T}H + H^{\mathsf T}A, surjectief op symmetrische matrices in elke AOnA \in O_n (gegeven symmetrische SS, neem H=12ASH = \frac12AS), dus is OnO_n een deelvariëteit van dimensie n2n(n+1)2=n(n1)2n^2 - \frac{n(n+1)}2 = \frac{n(n-1)}2 met

TAOn=kerDF(A)={H:ATH antisymmetrisch}={AK:KT=K};T_AO_n = \ker DF(A) = \{H : A^{\mathsf T}H \text{ antisymmetrisch}\} = \{AK : K^{\mathsf T} = -K\} ;

in A=IA = I zijn dit de antisymmetrische matrices.

(c) exp(tK)Texp(tK)=exp(tKT)exp(tK)=exp(tK)exp(tK)=I\exp(tK)^{\mathsf T}\exp(tK) = \exp(tK^{\mathsf T}) \exp(tK) = \exp(-tK)\exp(tK) = I (de twee matrices ±tK\pm tK commuteren, dus is het product van exponentialen de exponentiaal van de som): exp(tK)On\exp(tK) \in O_n, en γ(t)=Aexp(tK)\gamma(t) = A\exp(tK) is een kromme in OnO_n met γ(0)=A\gamma(0) = A, γ(0)=AK\gamma'(0) = AK. Als KK over de antisymmetrische matrices loopt, veegt AKAK TAOnT_AO_n: de exponentiaal realiseert de hele raakruimte door expliciete krommen — de Lie-groep-snelweg die Probleem 20.1 voor SO(3)SO(3) uitbuit.

Oefening 20.12 ★★

(Kritieke punten van de afstand) Zij MRnM \subseteq \R^n een deelvariëteit en pMp \notin M. Toon aan dat als x0Mx_0 \in M de afstand tot pp minimaliseert (zo’n punt bestaat wanneer MM gesloten en niet-leeg is — waarom?), dan

px0    Tx0Mp - x_0 \;\perp\; T_{x_0}M

(differentieer tγ(t)p2t \mapsto \norm{\gamma(t) - p}^2 langs krommen in MM). Leid af: het dichtstbijzijnde punt op een sfeer ligt op de straal door het middelpunt; en gebruik de voorwaarde om de afstand van p=(2,0)p = (2, 0) tot de parabool y=x2y = x^2 te berekenen (reduceer tot een derdegraads en los die numeriek op tot drie cijfers).

Oplossing

Oplossing van Oefening 20.12.

Bestaan: snijd MM met een grote gesloten bal rond pp om een niet-lege compacte te krijgen; de continue afstand bereikt daar haar minimum, en punten buiten de bal zijn verder. Eerste-orde-voorwaarde: voor een kromme γ\gamma in MM met γ(0)=x0\gamma(0) = x_0 is de functie h(t)=γ(t)p2h(t) = \norm{\gamma(t) - p}^2 differentieerbaar met een minimum in 00:

0=h(0)=2γ(0), x0p,0 = h'(0) = 2\,\langle\gamma'(0),\ x_0 - p\rangle,

en γ(0)\gamma'(0) veegt Tx0MT_{x_0}M: px0Tx0Mp - x_0 \perp T_{x_0}M. Sfeer S(c,r)S(c, r): de raakruimte in x0x_0 is (x0c)(x_0 - c)^\perp, dus px0x0cp - x_0 \parallel x_0 - c: x0x_0 ligt op de lijn door cc en pp, op afstand rr van cc — het straalpunt, zoals de meetkunde eist. Parabool: in x0=(x,x2)x_0 = (x, x^2) wordt de raaklijn opgespannen door (1,2x)(1, 2x); orthogonaliteit op px0=(2x,x2)p - x_0 = (2 - x, -x^2) leest

(2x)2x3=0,d.w.z.2x3+x2=0,(2 - x) - 2x^3 = 0, \qquad\text{d.w.z.}\qquad 2x^3 + x - 2 = 0,

met unieke reële wortel (x2x3+xx \mapsto 2x^3 + x is strikt stijgend) x0,835x \approx 0{,}835; dan x0(0,835,0,698)x_0 \approx (0{,}835, 0{,}698) en d(p,M)=(20,835)2+0,69821,358d(p, M) = \sqrt{(2 - 0{,}835)^2 + 0{,}698^2} \approx 1{,}358.

20.5 Probleem: SO(3)SO(3) en de quaternionen

Probleem 20.1

Weekendprobleem — rotaties, de groep S3S^3, en de dubbeldekking

De quaternionen H={t+xi+yj+zk}\mathbb H = \{t + x\mathrm i + y\mathrm j + z\mathrm k\} — de algebra waarvan de eenheidsgroep Q8Q_8 van Probleem 1.1 bevat — parametriseren driedimensionale rotaties tweevoudig: de afbeelding “conjugeren met een eenheidsquaternion” is een surjectief morfisme S3SO(3)S^3 \to SO(3) met kern {±1}\{\pm1\}. We bouwen alles. Herinner/definieer: vermenigvuldiging is R\R-bilineair met i2=j2=k2=ijk=1\mathrm i^2 = \mathrm j^2 = \mathrm k^2 = \mathrm{ijk} = -1; de geconjugeerde van q=t+xi+yj+zkq = t + x\mathrm i + y\mathrm j + z\mathrm k is qˉ=txiyjzk\bar q = t - x\mathrm i - y\mathrm j - z\mathrm k; N(q)=qqˉ=t2+x2+y2+z2N(q) = q\bar q = t^2 + x^2 + y^2 + z^2.

Deel I — De algebra H\mathbb H en de groep S3S^3.

  1. Verifieer dat H\mathbb H een associatieve R\R-algebra is met centrum R\R, dat pq=qˉpˉ\overline{pq} = \bar q\,\bar p, en dat N(pq)=N(p)N(q)N(pq) = N(p)N(q) (één schone route: representeer qq als de 2×22\times2 complexe matrix (αββˉαˉ)\bigl(\begin{smallmatrix}\alpha & \beta\\ -\bar\beta & \bar\alpha\end{smallmatrix}\bigr), q=α+βjq = \alpha + \beta\mathrm j, en gebruik det\det).
  2. Leid af dat elke q0q \neq 0 inverteerbaar is (q1=qˉ/N(q)q^{-1} = \bar q/N(q)): H\mathbb H is een (niet-commutatief) lichaam, en S3={N(q)=1}S^3 = \{N(q) = 1\} is een groep — en een compacte 33-deelvariëteit van R4\R^4 (Voorbeeld 20.4).

Deel II — Het rotatiemorfisme. Identificeer R3\R^3 met de pure quaternionen P={xi+yj+zk}P = \{x\mathrm i + y\mathrm j + z\mathrm k\}, en voor qS3q \in S^3 definieer ρq(v)=qvqˉ\rho_q(v) = q\,v\,\bar q.

  1. Toon aan dat ρq\rho_q PP naar PP stuurt (pure quaternionen zijn die met vˉ=v\bar v = -v), R\R-lineair is, de norm bewaart, en dat ρ ⁣:qρq\rho \colon q \mapsto \rho_q een groeps-morfisme S3O(3)S^3 \to O(3) is.
  2. Bereken de kern: ρq=id\rho_q = \mathrm{id} desda qq commuteert met i,j,k\mathrm i, \mathrm j, \mathrm k desda qRS3={±1}q \in \R\cap S^3 = \{\pm1\}.
  3. Schrijf q=cosθ2+sinθ2uq = \cos\frac\theta2 + \sin\frac\theta2\,u met uPu \in P, N(u)=1N(u) = 1 (waarom is dit altijd mogelijk voor qS3q \in S^3?). Toon aan dat ρq\rho_q uu vastzet en, op het vlak uPu^\perp\cap P, werkt als de rotatie over hoek θ\theta (bereken ρq(w)\rho_q(w) voor wuw \perp u met uw=wuuw = -wu voor orthogonale pure eenheden — bewijs deze identiteit uit de vermenigvuldigingstafel, of uit uw+wu=2u,wuw + wu = -2\langle u, w\rangle).
  4. Concludeer: imρSO(3)\operatorname{im}\rho \subseteq SO(3) (elke ρq\rho_q is een rotatie met as en hoek zoals berekend — determinant +1+1 door continuïteit van qdetρqq \mapsto \det\rho_q op de samenhangende S3S^3, of rechtstreeks), en ρ\rho is op SO(3)SO(3): elke rotatie van R3\R^3 heeft een as (bewijs: een reële 3×33\times3 orthogonale matrix met det=1\det = 1 heeft eigenwaarde 11 — beschouw het karakteristieke polynoom) en is dus enige ρq\rho_q. Samenvatting:

    SO(3)    S3/{±1}.SO(3) \;\cong\; S^3/\{\pm 1\} .

Deel III — SO(3)SO(3) als deelvariëteit; Rodrigues.

  1. Toon aan dat SO(3)SO(3) een compacte 33-dimensionale deelvariëteit is van M3(R)M_3(\R) met TISO(3)=T_ISO(3) = antisymmetrische matrices (Oefening 20.5; de determinantvoorwaarde selecteert een unie van componenten).
  2. Voor de antisymmetrische matrix AuA_u geassocieerd met uR3u \in \R^3 (Auv=uvA_uv = u\wedge v, het kruisproduct), bewijs de formule van Rodrigues:

    eθAu=I+sinθAu+(1cosθ)Au2(u=1)\eu^{\theta A_u} = I + \sin\theta\,A_u + (1 - \cos\theta)\,A_u^2 \qquad (\norm u = 1)

    (uit Au3=AuA_u^3 = -A_u: splits de exponentiaalreeks langs de machten van AuA_u), en identificeer haar als de rotatie om as uu over hoek θ\theta. Leid af dat exp\exp de antisymmetrische matrices op SO(3)SO(3) afbeeldt.

  3. Relateer de twee parametriseringen: toon aan dat tρq(t)t \mapsto \rho_{q(t)} met q(t)=cost2+sint2uq(t) = \cos\frac t2 + \sin\frac t2\,u een éénparameter-groep van rotaties is waarvan de afgeleide in t=0t = 0 AuA_u is — de quaternionische en matrix-exponentialen vertellen hetzelfde verhaal op half en vol tempo respectievelijk.

Deel IV — De dubbeldekking, gevoeld.

  1. Toon aan dat het pad q(t)=cost2+sint2kq(t) = \cos\frac t2 + \sin\frac t2\,\mathrm k, t[0,2π]t \in \intcc0{2\pi}, een lus is in SO(3)SO(3) (haar beeld ρq(t)\rho_{q(t)} keert terug naar de identiteit) waarvan de quaternionische lift geen lus is: q(2π)=q(0)q(2\pi) = -q(0). Doorgaan tot t=4πt = 4\pi sluit de lift. Leg in een korte alinea uit wat dit zegt: een 2π2\pi-rotatie is niet continu ongedaan te maken terwijl een 4π4\pi-rotatie dat wel is (de riemtruc), omdat lussen van SO(3)SO(3) worden gedetecteerd in haar dubbeldekking S3S^3.
  2. Leid ook het praktische dividend af: samenstelling van rotaties == vermenigvuldiging van quaternionen (44 vermenigvuldigingen aan data in plaats van 99, geen drift van orthogonaliteit) — verifieer op de samenstelling van twee kwartdraaien om i\mathrm i en j\mathrm j: bereken as en hoek van het product.

Deel V — De expliciete matrix: Euler–Rodrigues. Schrijf q=a+bi+cj+dkS3q = a + b\mathrm i + c\mathrm j + d\mathrm k \in S^3, zodat a2+b2+c2+d2=1a^2 + b^2 + c^2 + d^2 = 1.

  1. Bereken ρq(i)\rho_q(\mathrm i) volledig uit de vermenigvuldigingstafel; verkrijg dan ρq(j)\rho_q(\mathrm j) en ρq(k)\rho_q(\mathrm k) door de cyclische substitutie ijki\mathrm i \to \mathrm j \to \mathrm k \to \mathrm i, (b,c,d)(c,d,b)(b, c, d) \to (c, d, b) (rechtvaardig: cycleren van i,j,k\mathrm i, \mathrm j, \mathrm k zet zich voort tot een automorfisme van H\mathbb H, omdat de definiërende relaties cyclisch symmetrisch zijn). Concludeer dat de matrix van ρq\rho_q in de basis (i,j,k)(\mathrm i, \mathrm j, \mathrm k) de Euler–Rodrigues-matrix is

    Rq=(a2+b2c2d22(bcad)2(bd+ac)2(bc+ad)a2b2+c2d22(cdab)2(bdac)2(cd+ab)a2b2c2+d2).R_q = \begin{pmatrix} a^2 + b^2 - c^2 - d^2 & 2(bc - ad) & 2(bd + ac)\\ 2(bc + ad) & a^2 - b^2 + c^2 - d^2 & 2(cd - ab)\\ 2(bd - ac) & 2(cd + ab) & a^2 - b^2 - c^2 + d^2 \end{pmatrix}.
  2. (Een rotatie achterstevoren lezen) Toon aan dat

    trRq=4a21=1+2cosθ,12(RqRqT)=sinθAu,\operatorname{tr}R_q = 4a^2 - 1 = 1 + 2\cos\theta, \qquad \tfrac12\bigl(R_q - R_q^{\mathsf T}\bigr) = \sin\theta\,A_u,

    in de notatie van vragen 5 en 8. Leid een algoritme af dat ±q\pm q uit een rotatiematrix RR terugwint: de hoek uit het spoor; de as uit het antisymmetrische deel wanneer 0<θ<π0 < \theta < \pi; en, wanneer θ=π\theta = \pi, bewijs en gebruik de identiteit R+I=2uuTR + I = 2\,uu^{\mathsf T}.

  3. Evalueer RqR_q voor het product q=12(1+i+j+k)q = \frac12(1 + \mathrm i + \mathrm j + \mathrm k) van vraag 11: een permutatiematrix verschijnt. Identificeer de rotatie en verenig met de as en hoek gevonden in vraag 11.

Deel VI — Binnen S3S^3: SU(2)SU(2), conjugatieklassen, exponentialen.

  1. Toon aan dat de matrixrepresentatie van vraag 1 (noem haar Φ\Phi) restrictie tot een groepsisomorfisme van S3S^3 op de speciale unitaire groep

    SU(2)={UM2(C):UU=I, detU=1}SU(2) = \bigl\{U \in M_2(\C) : U^*U = I,\ \det U = 1\bigr\}

    (voor surjectiviteit, schrijf de vergelijkingen U1=UU^{-1} = U^* en detU=1\det U = 1 uit voor een algemene 2×22\times2 complexe matrix).

  2. Toon aan dat het reële deel conjugatie-invariant is op S3S^3Re(pqpˉ)=Req\operatorname{Re}(pq\bar p) = \operatorname{Re}q voor alle pS3p \in S^3 — en, omgekeerd, dat twee eenheidsquaternionen met hetzelfde reële deel geconjugeerd zijn in S3S^3 (reduceer tot het verplaatsen van één pure eenheidsas op een andere, wat Deel II levert). Beschrijf de conjugatieklassen van S3S^3 meetkundig; vertaal naar SU(2)SU(2) (niveaucurven van het spoor); en projecteer door ρ\rho: twee rotaties zijn geconjugeerd in SO(3)SO(3) dan en slechts dan als zij dezelfde hoek θ[0,π]\theta \in \intcc0\pi hebben.
  3. Definieer exp\exp op H\mathbb H via de exponentiaalreeks; controleer absolute convergentie, met pq=pq\abs{pq} = \abs p\,\abs q voor q=N(q)\abs q = \sqrt{N(q)}. Toon, voor een pure eenheid uu en θR\theta \in \R,

    exp(θu)=cosθ+sinθu,\exp(\theta u) = \cos\theta + \sin\theta\,u ,

    leid af dat exp\exp het hypervlak PP op S3S^3 afbeeldt, en controleer dat ρexp(su)=e2sAu\rho_{\exp(su)} = \eu^{2sA_u}: het halfhoek-fenomeen van vraag 9 opnieuw.

  4. Voor pure quaternionen v,wv, w bewijs de productregel vw=v,w+vwvw = -\langle v, w\rangle + v\wedge w, dus de commutatoridentiteit vwwv=2vwvw - wv = 2\,v\wedge w; bewijs ook [Av,Aw]=Avw[A_v, A_w] = A_{v\wedge w} voor de matrices van vraag 8. Concludeer dat de afgeleide van ρ\rho in 11 langs de krommen texp(tv)t \mapsto \exp(tv) het lineaire isomorfisme v2Avv \mapsto 2A_v is van PP op de antisymmetrische matrices, en dat het de quaternioncommutator naar de matrixcommutator transporteert.

Deel VII — Globale structuur.

  1. (Geen continue sectie) Stel s ⁣:SO(3)S3s \colon SO(3) \to S^3 is continu met ρs=id\rho \circ s = \operatorname{id}. Voor de lus R(t)=ρq(t)R(t) = \rho_{q(t)} van vraag 10, stel ε(t)=s(R(t))q(t)1\varepsilon(t) = s(R(t))\,q(t)^{-1} voor t[0,2π]t \in \intcc0{2\pi}. Toon aan dat ε\varepsilon continu is met waarden in {±1}\{\pm1\}, en leid een tegenspraak af: er is geen continue globale keuze van een eenheidsquaternion die elke rotatie voorstelt.
  2. (Het balmodel) Zij BˉR3\bar B \subseteq \R^3 de gesloten bal van straal π\pi en E(v)=eAvE(v) = \eu^{A_v}, met E(0)=IE(0) = I. Toon aan dat EE Bˉ\bar B op SO(3)SO(3) afbeeldt, injectief is op de open bal, en op de randbol precies antipoden identificeert: E(πu)=E(πu)=2uuTIE(\pi u) = E(-\pi u) = 2uu^{\mathsf T} - I, zonder andere samenvallen. Dus is SO(3)SO(3) de bal met antipodale randpunten gelijmd — de projectieve ruimte RP3\mathbb{RP}^3 — en een diameter wordt de niet-contractiele lus van vraag 10.
  3. Toon aan dat ρpρqρp1=ρpqpˉ\rho_p\rho_q\rho_p^{-1} = \rho_{pq\bar p}; dat de involuties van SO(3)SO(3) (de RIR \neq I met R2=IR^2 = I) precies de halve-draaien ρw\rho_w zijn met ww een pure eenheidsquaternion; en dat het centrum van SO(3)SO(3) triviaal is.
  4. Toon aan dat elke rotatie een product is van twee halfdraaien: voor q=cosθ2+sinθ2uq = \cos\frac\theta2 + \sin\frac\theta2\,u, kies een pure eenheid wuw \perp u, controleer dat w=qww' = qw opnieuw een pure eenheidsquaternion is, en verifieer ρq=ρwρw\rho_q = \rho_{w'}\rho_w. Waar liggen de twee assen, en welke hoek maken zij?
  5. Concludeer de topologische samenvatting: SO(3)SO(3) is compact en padsamenhangend (geef twee bewijzen: continu beeld van S3S^3 onder ρ\rho; beeld van exp\exp), terwijl O(3)O(3) precies twee samenhangscomponenten heeft, elk homeomorf met SO(3)SO(3).
  6. (Een samenstelling, drie manieren) Zij R1R_1 de rotatie over π2\frac\pi2 om de zz-as en R2R_2 de rotatie over π2\frac\pi2 om de xx-as. Bereken as en hoek van R2R1R_2R_1: (i) door de twee 3×33\times3-matrices te vermenigvuldigen en spoor/antisymmetrisch deel te gebruiken (Deel V); (ii) door de bijbehorende eenheidsquaternionen q2q1q_2q_1 te vermenigvuldigen. Controleer dat de twee antwoorden overeenkomen: hoek 2π3\frac{2\pi}3, as 13(1,1,1)\frac1{\sqrt3}(1, -1, 1).
  7. (De Cayley-transformatie) Voor KK antisymmetrisch, toon aan dat I+KI + K inverteerbaar is en

    C(K)=(IK)(I+K)1SO(n),C(K) = (I - K)(I + K)^{-1} \in SO(n),

    met 1-1 nooit een eigenwaarde van C(K)C(K); toon aan dat KC(K)K \mapsto C(K) een bijectie is van antisymmetrische matrices op {RSO(n):1SpR}\{R \in SO(n) : -1 \notin \operatorname{Sp}R\}, met inverse R(IR)(I+R)1R \mapsto (I - R)(I + R)^{-1}. (Een rationale kaart van SO(n)SO(n), metgezel van de transcendentale exp\exp van Oefening 20.11.)

Oplossing

Oplossing van Probleem 20.1.

1. Beeld q=t+xi+yj+zk(αββˉαˉ)q = t + x\mathrm i + y\mathrm j + z\mathrm k \mapsto \bigl(\begin{smallmatrix}\alpha & \beta\\ -\bar\beta & \bar\alpha\end{smallmatrix}\bigr) met α=t+ix\alpha = t + \iu x, β=y+iz\beta = y + \iu z: men controleert dat 1,i,j,k1, \mathrm i, \mathrm j, \mathrm k naar II, (i00i)\bigl(\begin{smallmatrix}\iu & 0\\ 0 & -\iu\end{smallmatrix}\bigr), (0110)\bigl(\begin{smallmatrix}0 & 1\\ -1 & 0\end{smallmatrix}\bigr), (0ii0)\bigl(\begin{smallmatrix}0 & \iu\\ \iu & 0\end{smallmatrix}\bigr) gaan, waarvan de producten de quaternionentafel reproduceren: de afbeelding is een injectief algebramorfisme, dus erft H\mathbb H associativiteit; N(q)=α2+β2=detN(q) = \abs\alpha^2 + \abs\beta^2 = \det is multiplicatief, en conjugatie correspondeert met de adjunct-transpose, wat pq=qˉpˉ\overline{pq} = \bar q\bar p geeft. Centrum: commuteren met i\mathrm i dwingt y=z=0y = z = 0, met j\mathrm j dwingt x=0x = 0: R\R.

2. qqˉ=N(q)q\bar q = N(q): voor q0q \neq 0, q1=qˉ/N(q)q^{-1} = \bar q/N(q): een delingsalgebra. Op S3S^3: N(pq)=1N(pq) = 1 en N(q1)=1N(q^{-1}) = 1: een groep; en S3R4S^3 \subseteq \R^4 is de eenheidssfeer: een compacte 33-deelvariëteit.

3. vv is puur desda vˉ=v\bar v = -v; dan qvqˉ=qvˉqˉ=qvqˉ\overline{qv\bar q} = q\bar v\bar q = -qv\bar q: ρq\rho_q bewaart PP. Lineariteit is duidelijk; N(qvqˉ)=N(q)N(v)N(q)=N(v)N(qv\bar q) = N(q)N(v)N(q) = N(v): een isometrie van (P,N)(R3,2)(P, N) \cong (\R^3, \norm\cdot^2): ρqO(3)\rho_q \in O(3). En ρpq(v)=pqvpq=p(qvqˉ)pˉ=ρp(ρq(v))\rho_{pq}(v) = pqv\overline{pq} = p(qv\bar q)\bar p = \rho_p(\rho_q(v)): een morfisme.

4. ρq=id\rho_q = \mathrm{id} desda qv=vqqv = vq voor alle pure vv, desda qq commuteert met i,j,k\mathrm i, \mathrm j, \mathrm k, desda qq centraal is (vraag 1): qRS3={±1}q \in \R\cap S^3 = \{\pm1\}.

5. Schrijf q=t+pq = t + p (tRt \in \R, pp puur): 1=N(q)=t2+N(p)1 = N(q) = t^2 + N(p), dus t=cosθ2t = \cos\frac\theta2 en p=sinθ2up = \sin\frac\theta2\,u met N(u)=1N(u) = 1 voor enige θ\theta (als p=0p = 0, werkt q=±1q = \pm1 triviaal). Omdat u2=N(u)=1u^2 = -N(u) = -1, commuteren qq en uu, en ρq(u)=quqˉ=uqqˉ=u\rho_q(u) = qu\bar q = uq\bar q = u: de as. Voor pure eenheden wuw \perp u: de productregel vw=v,w+vwvw = -\langle v, w\rangle + v\wedge w (coördinaatsgewijs uit de tafel) geeft uw=uw=wuuw = u\wedge w = -wu. Dan

ρq(w)=(cosθ2+sinθ2u)w(cosθ2sinθ2u)=cosθw+sinθ(uw),\rho_q(w) = \bigl(\cos\tfrac\theta2 + \sin\tfrac\theta2u\bigr)\,w\,\bigl(\cos\tfrac\theta2 - \sin\tfrac\theta2u\bigr) = \cos\theta\,w + \sin\theta\,(u\wedge w),

met uwu=u2w=wuwu = -u^2w = w en de dubbele-hoekformules: de rotatie over hoek θ\theta in het georiënteerde vlak (w,uw)(w, u\wedge w).

6. Elke ρq\rho_q is een rotatie om uu over θ\theta: in de orthonormale basis (u,w,uw)(u, w, u\wedge w) heeft haar matrix determinant +1+1: imρSO(3)\operatorname{im}\rho \subseteq SO(3). Surjectiviteit: een matrix RSO(3)R \in SO(3) heeft 11 als eigenwaarde, want

det(RI)=detRdet(IRT)=det(IR)=(1)3det(RI),\det(R - I) = \det R\,\det(I - R^{\mathsf T}) = \det(I - R) = (-1)^3\det(R - I),

dus det(RI)=0\det(R - I) = 0. Neem een eenheidseigenvector uu; RR bewaart uu^\perp en restrictie daar is een rotatie over enige hoek θ\theta (vlak orthogonaal, determinant 11): R=ρqR = \rho_q voor q=cosθ2+sinθ2uq = \cos\frac\theta2 + \sin\frac\theta2\,u. Met vraag 4 en de eerste isomorfismestelling (Stelling 1.3): SO(3)S3/{±1}SO(3) \cong S^3/\{\pm1\}.

7. O3O_3 is een compacte 33-dimensionale deelvariëteit (Oefening 20.5); det\det is continu erop met waarden in {±1}\{\pm1\}, dus is SO(3)=O3{det=1}SO(3) = O_3\cap\{\det = 1\} open en gesloten in O3O_3: een unie van samenhangscomponenten, dus zelf een compacte 33-deelvariëteit, met dezelfde raakruimte in II: de antisymmetrische matrices.

8. Au2v=u(uv)=u,vuvA_u^2v = u\wedge(u\wedge v) = \langle u, v\rangle u - v (eenheid uu), dus Au3v=u(u,vuv)=uvA_u^3v = u\wedge(\langle u,v\rangle u - v) = -u\wedge v: Au3=AuA_u^3 = -A_u. Splits de exponentiaalreeks naar resten van machten mod de relatie A3=AA^3 = -A:

eθAu=I+(θθ33!+)Au+(θ22!θ44!+)Au2=I+sinθAu+(1cosθ)Au2.\eu^{\theta A_u} = I + \Bigl(\theta - \frac{\theta^3}{3!} + \cdots\Bigr)A_u + \Bigl(\frac{\theta^2}{2!} - \frac{\theta^4}{4!} + \cdots\Bigr)A_u^2 = I + \sin\theta\,A_u + (1 - \cos\theta)\,A_u^2 .

Op uu: Auu=0A_uu = 0: vast. Op wuw \perp u: eθAuw=w+sinθuw+(1cosθ)(w)=cosθw+sinθuw\eu^{\theta A_u}w = w + \sin\theta\,u\wedge w + (1 - \cos\theta)(-w) = \cos\theta\,w + \sin\theta\,u\wedge w: de rotatie om as uu, hoek θ\theta — Rodrigues. Elke rotatie heeft deze vorm (vraag 6): exp\exp is op SO(3)SO(3) vanaf de antisymmetrische matrices.

9. Met q(t)=cost2+sint2uq(t) = \cos\frac t2 + \sin\frac t2\,u: vraag 5 toont dat ρq(t)\rho_{q(t)} de rotatie is om as uu over hoek tt, d.w.z. ρq(t)=etAu\rho_{q(t)} = \eu^{tA_u}, waarvan de afgeleide in t=0t = 0 AuA_u is. Het quaternion loopt op halve hoek — het analytische spoor van de dubbeldekking.

10. ρq(t)\rho_{q(t)} is de rotatie om k\mathrm k over hoek tt: in t=2πt = 2\pi keert zij terug naar de identiteit — een lus in SO(3)SO(3). Haar lift voldoet aan q(2π)=cosπ=1=q(0)q(2\pi) = \cos\pi = -1 = -q(0): het gelifte pad is niet gesloten; pas in t=4πt = 4\pi keert qq terug naar 11. Interpretatie: de lus van volle rotaties is niet contractiel in SO(3)SO(3) — haar lift eindigt op het andere blad van de dekking — terwijl de dubbele lus dat wel is; een lichaam verbonden met zijn omgeving door riemen (de riemtruc) keert terug naar een onverdraaide toestand na 4π4\pi maar niet na 2π2\pi. Rotatiegroepen onthouden de pariteit van volle slagen; S3S^3, enkelvoudig samenhangend, is waar dat geheugen leeft.

11. Kwartdraaien: qi=cosπ4+sinπ4iq_{\mathrm i} = \cos\frac\pi4 + \sin\frac\pi4\,\mathrm i, qj=cosπ4+sinπ4jq_{\mathrm j} = \cos\frac\pi4 + \sin\frac\pi4\,\mathrm j. Product (eerst de j\mathrm j-draai):

qiqj=12(1+i)(1+j)=12(1+i+j+k),q_{\mathrm i}q_{\mathrm j} = \tfrac12(1 + \mathrm i)(1 + \mathrm j) = \tfrac12\bigl(1 + \mathrm i + \mathrm j + \mathrm k\bigr),

van norm 11, met cosθ2=12\cos\frac\theta2 = \frac12: θ=2π3\theta = \frac{2\pi}3, en as u=i+j+k3u = \frac{\mathrm i + \mathrm j + \mathrm k}{\sqrt3} (het pure deel genormaliseerd). Twee opeenvolgende kwartdraaien om orthogonale assen maken een 120120^\circ-rotatie om de hoofddiagonaal van de kubus — vier reële vermenigvuldigingen aan boekhouding, orthogonaliteit exact bewaard: waarom vluchtsoftware en grafische engines rotaties via quaternionen samenstellen.

12. Uit de tafel, ji=k\mathrm{ji} = -\mathrm k en ki=j\mathrm{ki} = \mathrm j, dus

qi=aib+c(ji)+d(ki)=b+ai+djck.q\,\mathrm i = a\mathrm i - b + c(\mathrm{ji}) + d(\mathrm{ki}) = -b + a\mathrm i + d\mathrm j - c\mathrm k .

Vermenigvuldigen met qˉ=abicjdk\bar q = a - b\mathrm i - c\mathrm j - d\mathrm k met de scalaire–vectorregel (t1+p1)(t2+p2)=t1t2p1,p2+t1p2+t2p1+p1p2(t_1 + p_1)(t_2 + p_2) = t_1t_2 - \langle p_1, p_2\rangle + t_1p_2 + t_2p_1 + p_1\wedge p_2, waar p1=(a,d,c)p_1 = (a, d, -c) en p2=(b,c,d)p_2 = (-b, -c, -d): het scalaire deel is ab(ab)=0-ab - (-ab) = 0 (puur, zoals het moet), en het vectordeel is

b(b,c,d)+a(a,d,c)+(c2d2, bc+ad, bdac)=(a2+b2c2d2, 2(bc+ad), 2(bdac)):b(b, c, d) + a(a, d, -c) + (-c^2 - d^2,\ bc + ad,\ bd - ac) = \bigl(a^2 + b^2 - c^2 - d^2,\ 2(bc + ad),\ 2(bd - ac)\bigr):

de eerste kolom van RqR_q. De cyclische afbeelding σ(i)=j\sigma(\mathrm i) = \mathrm j, σ(j)=k\sigma(\mathrm j) = \mathrm k, σ(k)=i\sigma(\mathrm k) = \mathrm i bewaart de relaties i2=j2=k2=ijk=1\mathrm i^2 = \mathrm j^2 = \mathrm k^2 = \mathrm{ijk} = -1 (het woord ijk\mathrm{ijk} is cyclisch invariant op de relatie ijk=jki\mathrm{ijk} = \mathrm{jki}, die in elke ring geldt: conjugeren van ijk=1\mathrm{ijk} = -1 met de inverteerbare i\mathrm i), dus zet σ\sigma zich voort tot een R\R-algebra-automorfisme, en σ(ρq(v))=ρσ(q)(σ(v))\sigma(\rho_q(v)) = \rho_{\sigma(q)}(\sigma(v)). Ontwikkelend is het beeld van j\mathrm j de eerste-kolom-formule na de substitutie (b,c,d)(c,d,b)(b, c, d) \to (c, d, b) met de basis herlabeld ijki\mathrm i \to \mathrm j \to \mathrm k \to \mathrm i, precies de getoonde tweede kolom; één slag meer geeft de derde.

13. Het spoor sommerend, trRq=3a2(b2+c2+d2)=4a21\operatorname{tr}R_q = 3a^2 - (b^2 + c^2 + d^2) = 4a^2 - 1 (eenheidsnorm), en met a=cosθ2a = \cos\frac\theta2: 4cos2θ21=1+2cosθ4\cos^2\frac\theta2 - 1 = 1 + 2\cos\theta. Antisymmetrisch deel: de drie onafhankelijke entries van RqRqTR_q - R_q^{\mathsf T} zijn 4ab,4ac,4ad4ab, 4ac, 4ad (in posities (3,2),(1,3),(2,1)(3,2), (1,3), (2,1)), dus 12(RqRqT)=Am\frac12(R_q - R_q^{\mathsf T}) = A_m met m=2a(b,c,d)=2cosθ2sinθ2u=sinθum = 2a\,(b, c, d) = 2\cos\frac\theta2\sin\frac\theta2\,u = \sin\theta\,u. Algoritme: θ=arccostrR12[0,π]\theta = \arccos\frac{\operatorname{tr}R - 1}{2} \in \intcc0\pi; als 0<θ<π0 < \theta < \pi, lees uu af uit RRT2sinθ\frac{R - R^{\mathsf T}}{2\sin\theta} en stel q=±(cosθ2+sinθ2u)q = \pm(\cos\frac\theta2 + \sin\frac\theta2 u); als θ=0\theta = 0, q=±1q = \pm1. Voor θ=π\theta = \pi: a=0a = 0, en Rodrigues (vraag 8) geeft R=I+2Au2=I+2(uuTI)=2uuTIR = I + 2A_u^2 = I + 2(uu^{\mathsf T} - I) = 2uu^{\mathsf T} - I, d.w.z. R+I=2uuTR + I = 2uu^{\mathsf T}; elke niet-nul kolom van R+IR + I, genormaliseerd, is ±u\pm u, en q=±uq = \pm u.

14. Met a=b=c=d=12a = b = c = d = \frac12: alle diagonale entries verdwijnen, 2(bcad)=02(bc - ad) = 0, 2(bd+ac)=12(bd + ac) = 1, 2(bc+ad)=12(bc + ad) = 1, 2(cdab)=02(cd - ab) = 0, 2(bdac)=02(bd - ac) = 0, 2(cd+ab)=12(cd + ab) = 1:

Rq=(001100010),R_q = \begin{pmatrix} 0 & 0 & 1\\ 1 & 0 & 0\\ 0 & 1 & 0 \end{pmatrix},

de cyclische permutatie e1e2e3e1e_1 \to e_2 \to e_3 \to e_1. Haar spoor is 0=1+2cosθ0 = 1 + 2\cos\theta, dus θ=2π3\theta = \frac{2\pi}3, en zij fixeert (1,1,1)(1,1,1): de rotatie over 120120^\circ om de hoofddiagonaal — precies het antwoord van vraag 11, nu zichtbaar als de matrix die de coördinaatassen cyclisch verwisselt.

15. Op de basis controleert men Φ(qˉ)=Φ(q)\Phi(\bar q) = \Phi(q)^* (de matrix van qˉ\bar q heeft α=αˉ\alpha' = \bar\alpha, β=β\beta' = -\beta, wat de geconjugeerde transpose is van (αββˉαˉ)\bigl(\begin{smallmatrix}\alpha & \beta\\ -\bar\beta & \bar\alpha\end{smallmatrix}\bigr)). Dus Φ(q)Φ(q)=Φ(qˉq)=N(q)I\Phi(q)^*\Phi(q) = \Phi(\bar qq) = N(q)I en detΦ(q)=α2+β2=N(q)\det\Phi(q) = \abs\alpha^2 + \abs\beta^2 = N(q): voor qS3q \in S^3, Φ(q)SU(2)\Phi(q) \in SU(2), en Φ\Phi is een injectief morfisme (vraag 1). Surjectiviteit: zij U=(αβγδ)U = \bigl(\begin{smallmatrix}\alpha & \beta\\ \gamma & \delta\end{smallmatrix}\bigr) met detU=1\det U = 1; dan U1=(δβγα)U^{-1} = \bigl(\begin{smallmatrix}\delta & -\beta\\ -\gamma & \alpha\end{smallmatrix}\bigr), en U1=U=(αˉγˉβˉδˉ)U^{-1} = U^* = \bigl(\begin{smallmatrix}\bar\alpha & \bar\gamma\\ \bar\beta & \bar\delta\end{smallmatrix}\bigr) dwingt δ=αˉ\delta = \bar\alpha, γ=βˉ\gamma = -\bar\beta, en dan 1=detU=α2+β21 = \det U = \abs\alpha^2 + \abs\beta^2: U=Φ(q)U = \Phi(q) voor het eenheidsquaternion qq met coördinaten α=a+ib\alpha = a + \iu b, β=c+id\beta = c + \iu d. Dus S3SU(2)S^3 \cong SU(2).

16. Reële scalaren zijn centraal en N(p)=1N(p) = 1 geeft pqpˉ=pqˉpˉ\overline{pq\bar p} = p\bar q\bar p, dus pqpˉ+pqpˉ=p(q+qˉ)pˉ=q+qˉpq\bar p + \overline{pq\bar p} = p(q + \bar q)\bar p = q + \bar q: het reële deel is invariant. Omgekeerd, zij Req=Req=a\operatorname{Re}q = \operatorname{Re}q' = a; dan hebben de pure delen dezelfde norm 1a2=s\sqrt{1 - a^2} = s. Als s=0s = 0, q=q=±1q = q' = \pm1. Als s>0s > 0, schrijf q=a+suq = a + su, q=a+suq' = a + su' met u,uu, u' pure eenheden; vraag 6 levert een rotatie die uu' naar uu draagt, d.w.z. pS3p \in S^3 met ρp(u)=u\rho_p(u') = u, en dan pqpˉ=a+sρp(u)=qpq'\bar p = a + s\rho_p(u') = q. De klassen van S3S^3 zijn dus {1}\{1\}, {1}\{-1\}, en voor elke a(1,1)a \in \intoo{-1}1 de 22-sfeer {a+su:u pure eenheid}\{a + su : u \text{ pure eenheid}\} van straal ss. Onder Φ\Phi, trΦ(q)=α+αˉ=2Req\operatorname{tr}\Phi(q) = \alpha + \bar\alpha = 2\operatorname{Re}q: de klassen van SU(2)SU(2) zijn de niveaucurven van het spoor. Projecterend: als q=pqpˉq' = pq\bar p dan ρq=ρpρqρp1\rho_{q'} = \rho_p\rho_q\rho_p^{-1}; omgekeerd dwingt ρq=ρpρqρp1=ρpqpˉ\rho_{q'} = \rho_p\rho_q\rho_p^{-1} = \rho_{pq\bar p} q=±pqpˉq' = \pm pq\bar p (kern), dus Req=±Req\operatorname{Re}q' = \pm\operatorname{Re}q, d.w.z. cosθ2=Req=Req=cosθ2\cos\frac{\theta'}2 = \abs{\operatorname{Re}q'} = \abs{\operatorname{Re}q} = \cos\frac\theta2 voor de hoeken in [0,π]\intcc0\pi: geconjugeerde rotaties hebben gelijke hoeken. Omgekeerd laten gelijke hoeken representanten toe met hetzelfde niet-negatieve reële deel, geconjugeerd door het bovenstaande: in SO(3)SO(3) is de conjugatieklasse van een rotatie precies haar hoek.

17. NN is multiplicatief, dus is \abs\cdot een multiplicatieve norm op HR4\mathbb H \cong \R^4 en qk=qk\abs{q^k} = \abs q^k: de reeks qk/k!\sum q^k/k! convergeert absoluut in de eindigdimensionale (dus volledige) ruimte, gedomineerd door qk/k!=eq\sum\abs q^k/k! = \eu^{\abs q}. Voor een pure eenheid uu: u2=1u^2 = -1, dus (θu)2m=(1)mθ2m(\theta u)^{2m} = (-1)^m\theta^{2m} en (θu)2m+1=(1)mθ2m+1u(\theta u)^{2m+1} = (-1)^m\theta^{2m+1}u; de reeks splitsend,

exp(θu)=m(1)mθ2m(2m)!+um(1)mθ2m+1(2m+1)!=cosθ+sinθu.\exp(\theta u) = \sum_m\frac{(-1)^m\theta^{2m}}{(2m)!} + u\sum_m\frac{(-1)^m\theta^{2m+1}}{(2m+1)!} = \cos\theta + \sin\theta\,u .

Elke qS3q \in S^3 is cosα+sinαu\cos\alpha + \sin\alpha\,u met α[0,π]\alpha \in \intcc0\pi (vraag 5): q=exp(αu)q = \exp(\alpha u), dus exp(P)=S3\exp(P) = S^3. Ten slotte exp(su)=coss+sinsu=q(2s)\exp(su) = \cos s + \sin s\,u = q(2s) in de notatie van vraag 9, en ρq(t)=etAu\rho_{q(t)} = \eu^{tA_u} daar: ρexp(su)=e2sAu\rho_{\exp(su)} = \eu^{2sA_u}.

18. vwvw coördinaatsgewijs ontwikkelen met de tafel: de producten ii=1\mathrm i\cdot\mathrm i = -1, … geven de scalaire (v1w1+v2w2+v3w3)-(v_1w_1 + v_2w_2 + v_3w_3), en de gemengde producten (ij=k\mathrm{ij} = \mathrm k, ji=k\mathrm{ji} = -\mathrm k, …) geven de vector (v2w3v3w2, v3w1v1w3, v1w2v2w1)(v_2w_3 - v_3w_2,\ v_3w_1 - v_1w_3,\ v_1w_2 - v_2w_1): vw=v,w+vwvw = -\langle v, w\rangle + v\wedge w. Het omgekeerde product aftrekken: vwwv=2vwvw - wv = 2\,v\wedge w (de scalaire delen heffen elkaar op, de kruisproducten tellen op). Voor de matrices, met a(bc)=ba,cca,ba\wedge(b\wedge c) = b\langle a, c\rangle - c\langle a, b\rangle:

[Av,Aw]x=v(wx)w(vx)=wv,xvw,x=(vw)x=Avwx.[A_v, A_w]x = v\wedge(w\wedge x) - w\wedge(v\wedge x) = w\langle v, x\rangle - v\langle w, x\rangle = (v\wedge w)\wedge x = A_{v\wedge w}x .

Afgeleide: exp(tv)=exp(tv)\overline{\exp(tv)} = \exp(-tv) (conjugatie is continu en negeert pure quaternionen), dus

 ⁣d ⁣dtt=0exp(tv)xexp(tv)=vxxv=2vx=2Avx:\frac{\dd}{\dd t}\Bigr|_{t=0}\exp(tv)\,x\,\exp(-tv) = vx - xv = 2\,v\wedge x = 2A_vx :

de differentiaal is v2Avv \mapsto 2A_v, een lineaire bijectie van PP op de antisymmetrische matrices, en [2Av,2Aw]=4Avw=2A2vw=2A[v,w][2A_v, 2A_w] = 4A_{v\wedge w} = 2A_{2v\wedge w} = 2A_{[v,w]} toont dat zij de quaternioncommutator naar de matrixcommutator draagt.

19. ρ\rho toepassend: ρ(ε(t))=ρ(s(R(t)))ρ(q(t))1=R(t)R(t)1=id\rho(\varepsilon(t)) = \rho(s(R(t)))\,\rho(q(t))^{-1} = R(t)R(t)^{-1} = \operatorname{id}, dus ε(t)kerρ={±1}\varepsilon(t) \in \ker\rho = \{\pm1\} (vraag 4). Als product van de continue afbeeldingen ts(R(t))t \mapsto s(R(t)) en tq(t)1=q(t)ˉt \mapsto q(t)^{-1} = \bar{q(t)} is ε\varepsilon continu op het samenhangende interval [0,2π]\intcc0{2\pi} met waarden in het discrete paar {±1}\{\pm1\}: zij is constant, zeg ε(t)ε\varepsilon(t) \equiv \varepsilon. Maar R(0)=R(2π)=IR(0) = R(2\pi) = I, dus s(R(0))=s(R(2π))s(R(0)) = s(R(2\pi)), terwijl s(R(0))=εq(0)=εs(R(0)) = \varepsilon\,q(0) = \varepsilon en s(R(2π))=εq(2π)=εs(R(2\pi)) = \varepsilon\,q(2\pi) = -\varepsilon: tegenspraak. Geen continue sectie bestaat: de tekenambigüiteit ±q\pm q is globaal, geen gebrek van een bepaalde formule.

20. Op: elke RSO(3)R \in SO(3) is eθAu\eu^{\theta A_u} voor enige eenheid uu en θ[0,2π]\theta \in \intcc0{2\pi} (vragen 6 en 8); als θ>π\theta > \pi, geeft Rodrigues eθAu=e(2πθ)Au\eu^{\theta A_u} = \eu^{(2\pi - \theta)A_{-u}} (beide gelijk I+sinθAu+(1cosθ)Au2I + \sin\theta A_u + (1 - \cos\theta)A_u^2, en Au=AuA_{-u} = -A_u met sin(2πθ)=sinθ\sin(2\pi - \theta) = -\sin\theta, cos(2πθ)=cosθ\cos(2\pi - \theta) = \cos\theta), dus R=E(v)R = E(v) met vπ\norm v \leq \pi. Injectief binnenin: als E(v)=E(v)IE(v) = E(v') \neq I met v,v<π\norm v, \norm{v'} < \pi, wint vraag 13 dezelfde hoek θ=v=v(0,π)\theta = \norm v = \norm{v'} \in \intoo0\pi uit het spoor en, omdat sinθ0\sin\theta \neq 0, dezelfde as uit het antisymmetrische deel: v=vv = v'; en E(v)=IE(v) = I dwingt θ{0}\theta \in \{0\} op de open bal. Rand: E(πu)=I+2Au2=2uuTIE(\pi u) = I + 2A_u^2 = 2uu^{\mathsf T} - I hangt van uu alleen af via uuTuu^{\mathsf T}, vandaar E(πu)=E(πu)E(\pi u) = E(-\pi u); omgekeerd geeft 2uuTI=2uuTI2uu^{\mathsf T} - I = 2u'u'^{\mathsf T} - I toegepast op uu u=u,uuu = \langle u', u\rangle u', dus u=±uu' = \pm u. Inwendig en rand botsen nooit (spoor >1> -1 versus =1= -1). Dus induceert EE een continue bijectie van de bal-met-antipodale-lijming — compact — op SO(3)SO(3): een homeomorfisme, en SO(3)RP3SO(3) \cong \mathbb{RP}^3. Een diameter van πu\pi u naar πu-\pi u heeft gelijmde eindpunten: het is een lus in SO(3)SO(3), en haar EE-beschrijving past bij de familie rotaties om uu van vraag 10 die een volle slag vegen.

21. ρ\rho is een morfisme en ρp1=ρp1=ρpˉ\rho_p^{-1} = \rho_{p^{-1}} = \rho_{\bar p}, dus ρpρqρp1=ρpqpˉ\rho_p\rho_q\rho_p^{-1} = \rho_{pq\bar p}; volgens vraag 16 bewaart conjugeren van een rotatie haar hoek en draait haar as door ρp\rho_p. Involuties: ρq2=ρq2=id\rho_q^2 = \rho_{q^2} = \operatorname{id} desda q2=±1q^2 = \pm1. Als q2=1q^2 = 1 dan (q1)(q+1)=q21=0(q - 1)(q + 1) = q^2 - 1 = 0 (centrale scalaren, dus is deze factorisatie geldig) en q=±1q = \pm1 in het delingslichaam H\mathbb H, wat ρq=I\rho_q = I geeft, uitgesloten; q2=1q^2 = -1 met q=a+suq = a + su geeft a2s2+2asu=1a^2 - s^2 + 2as\,u = -1, dus a=0a = 0, s=1s = 1: qq is een pure eenheid ww, en ρw\rho_w is de halfdraai om ww (hoek π\pi, vraag 5). Centrum: als ρq\rho_q commuteert met elke ρp\rho_p, dan ρpqpˉ=ρq\rho_{pq\bar p} = \rho_q, dus pqpˉ=ε(p)qpq\bar p = \varepsilon(p)\,q met ε(p){±1}\varepsilon(p) \in \{\pm1\}; pε(p)=(pqpˉ)q1p \mapsto \varepsilon(p) = (pq\bar p)q^{-1} is continu op de samenhangende S3S^3 en gelijk aan 11 in p=1p = 1, dus 1\equiv 1: qq commuteert met heel S3S^3, dus met heel H\mathbb H (herschaal), dus qRS3={±1}q \in \R \cap S^3 = \{\pm1\} (vraag 1) en ρq=I\rho_q = I: het centrum is triviaal.

22. Omdat uwu \perp w pure eenheden zijn, is uw=uwuw = u\wedge w puur (vraag 18), dus

w=qw=cosθ2w+sinθ2uww' = qw = \cos\tfrac\theta2\,w + \sin\tfrac\theta2\,u\wedge w

is puur, van norm qw=1\abs q\abs w = 1. Dan ww=qww=qw2=qw'w = qw\cdot w = qw^2 = -q, en

ρwρw=ρww=ρq=ρq.\rho_{w'}\rho_w = \rho_{w'w} = \rho_{-q} = \rho_q .

Beide assen ww en w=cosθ2w+sinθ2(uw)w' = \cos\frac\theta2 w + \sin\frac\theta2(u\wedge w) liggen in het vlak uu^\perp orthogonaal op de rotatieas, en w,w=cosθ2\langle w', w\rangle = \cos\frac\theta2: zij maken de halfhoek θ2\frac\theta2. Dit is de klassieke generatie: twee halfdraaien om assen die hoek θ2\frac\theta2 maken, stellen samen tot de rotatie over hoek θ\theta om hun gemeenschappelijke loodlijn.

23. Compactheid is vraag 7. Padsamenhang: SO(3)=ρ(S3)SO(3) = \rho(S^3) is het continue beeld van de padsamenhangende sfeer; alternatief, voor R=eAR = \eu^{A} met AA antisymmetrisch (vraag 8), is tetAt \mapsto \eu^{tA} een pad in SO(3)SO(3) van II naar RR (orthogonaal want (etA)T=etA(\eu^{tA})^{\mathsf T} = \eu^{-tA}, determinant 11 door continuïteit vanaf t=0t = 0). Voor O(3)O(3): det\det is continu op {±1}\{\pm1\}, dus is O(3)O(3) onsamenhangend, O(3)=SO(3)DSO(3)O(3) = SO(3) \sqcup D\,SO(3) voor vaste DD met detD=1\det D = -1 (b.v. D=ID = -I), en linkse vermenigvuldiging met DD is een homeomorfisme: precies twee componenten, elk een kopie van SO(3)SO(3). De twee-op-één ρ\rho, sectie-vrij door vraag 19, is dus een eerlijke dubbeldekking van een samenhangende compacte groep door de enkelvoudig samenhangende S3S^3 — de meetkunde achter de riemtruc.

24. (i) Matrices:

R1=(010100001),R2=(100001010),R2R1=(010001100).R_1 = \begin{pmatrix} 0 & -1 & 0\\ 1 & 0 & 0\\ 0 & 0 & 1 \end{pmatrix}, \quad R_2 = \begin{pmatrix} 1 & 0 & 0\\ 0 & 0 & -1\\ 0 & 1 & 0 \end{pmatrix}, \quad R_2R_1 = \begin{pmatrix} 0 & -1 & 0\\ 0 & 0 & -1\\ 1 & 0 & 0\end{pmatrix}.

Spoor 0=1+2cosθ0 = 1 + 2\cos\theta geeft cosθ=12\cos\theta = -\frac12: θ=2π3\theta = \frac{2\pi}3. Antisymmetrisch deel RRT2\frac{R - R^{\mathsf T}}2 heeft entries die sinθ(v3,v2,v1)\sin\theta\,(v_3, -v_2, v_1)-gewijs de as coderen: hier RRT2=12(011101110)\frac{R - R^{\mathsf T}}2 = \frac12\bigl(\begin{smallmatrix}0 & -1 & -1\\ 1 & 0 & -1\\ 1 & 1 & 0\end{smallmatrix}\bigr), wat (woordenboek van Deel V voor AvA_v) vsinθ=12(1,1,1)v\sin\theta = \frac12(1, -1, 1) leest; met sin2π3=32\sin\frac{2\pi}3 = \frac{\sqrt3}2: v=13(1,1,1)v = \frac1{\sqrt3}(1, -1, 1). (ii) Quaternionen: q1=cosπ4+sinπ4k=22(1+k)q_1 = \cos\frac\pi4 + \sin\frac\pi4\,k = \frac{\sqrt2}2(1 + k), q2=22(1+i)q_2 = \frac{\sqrt2}2(1 + i), en

q2q1=12(1+i)(1+k)=12(1+k+i+ik)=12(1+ij+k)q_2q_1 = \tfrac12(1 + i)(1 + k) = \tfrac12(1 + k + i + ik) = \tfrac12\bigl(1 + i - j + k\bigr)

(ik=jik = -j). Dus cosθ2=12\cos\frac\theta2 = \frac12: θ=2π3\theta = \frac{2\pi}3, en het vectordeel 12(ij+k)\frac12(i - j + k) heeft richting 13(1,1,1)\frac1{\sqrt3}(1, -1, 1) — hetzelfde antwoord, met de quaternionroute die één regel vermenigvuldiging vraagt in plaats van een matrixproduct: de praktische reden waarom vluchtsoftware attitudes in S3S^3 samenstelt.

25. I+KI + K inverteerbaar: (I+K)v=0(I + K)v = 0 geeft 0=v,v+v,Kv=v20 = \langle v, v\rangle + \langle v, Kv\rangle = \norm v^2 (antisymmetrie doodt de tweede term): v=0v = 0. Orthogonaliteit van C=C(K)C = C(K): met (I±K)T=IK(I \pm K)^{\mathsf T} = I \mp K en het feit dat alle vier matrices I±KI \pm K, (I±K)1(I \pm K)^{-1} commuteren (polynomiale uitdrukkingen in KK, plus limieten):

CTC=(I+K)T(IK)T(IK)(I+K)1=(IK)1(I+K)(IK)(I+K)1=I.C^{\mathsf T}C = (I + K)^{-\mathsf T}(I - K)^{\mathsf T} (I - K)(I + K)^{-1} = (I - K)^{-1}(I + K)(I - K)(I + K)^{-1} = I .

Determinant: det(IK)=det((IK)T)=det(I+K)\det(I - K) = \det\bigl((I - K)^{\mathsf T}\bigr) = \det(I + K), dus detC=1\det C = 1: CSO(n)C \in SO(n). Geen eigenwaarde 1-1: Cv=vCv = -v betekent (IK)w=(I+K)w(I - K)w = -(I + K)w voor w=(I+K)1vw = (I + K)^{-1}v, d.w.z. 2w=02w = 0: v=0v = 0. Inversie: uit C(I+K)=IKC(I + K) = I - K, los K(I+C)=ICK(I + C) = I - C op; omdat 1SpC-1 \notin \operatorname{Sp}C, is I+CI + C inverteerbaar en K=(IC)(I+C)1K = (I - C)(I + C)^{-1}, wat antisymmetrisch is wanneer CC orthogonaal is zonder eigenwaarde 1-1 (transponeer de uitdrukking en gebruik CT=C1C^{\mathsf T} = C^{-1}: KT=(IC1)(I+C1)1=(CI)(C+I)1=KK^{\mathsf T} = (I - C^{-1})(I + C^{-1})^{-1} = (C - I)(C + I)^{-1} = -K). De twee afbeeldingen zijn elkaars inverse per constructie: een globale rationale parametrisering van het dichte open stuk van SO(n)SO(n) dat eigenwaarde 1-1 vermijdt — geen reeks, geen trigonometrie, en in dimensie 33 is het de halfhoek-substitutie K=tanθ2AvK = \tan\frac\theta2\,A_v in vermomming.