Mathematics · Livro 5 · Bachelor Year 3

Matemática universitária — Graduação 3

Matemática universitária — Graduação 3 · Bachelor Year 3

20Subvariedades de ℝn

Esferas, toros, grupos de rotações: os habitats naturais da geometria e da mecânica não são espaços vetoriais, mas conjuntos curvos que, de perto, parecem planos. Este capítulo dá a essa frase um sentido preciso — as subvariedades de Rn\R^n — e o cálculo para trabalhar sobre elas. O alicerce é o teorema da função inversa, aqui demonstrado pelo ponto fixo de Banach; todo o resto é mudança de coordenadas: as quatro descrições equivalentes de uma subvariedade (retificação local, conjuntos de nível, gráficos, parametrizações), os espaços tangentes e a otimização com restrições por multiplicadores de Lagrange — que, como demonstração de despedida, redemonstra o teorema espectral das matrizes simétricas em três linhas de geometria. O problema de fim de semana constrói o grupo das rotações SO(3)SO(3) e seu recobrimento duplo quaterniônico: a álgebra (o Q8Q_8 do Capítulo 1, crescido) encontrando a geometria.

20.1 O teorema da função inversa

Teorema 20.1 (Teorema da função inversa)

Sejam URnU \subseteq \R^n aberto, f ⁣:URnf \colon U \to \R^n de classe C1\mathcal C^1 e aUa \in U com Df(a)Df(a) inversível. Então existem abertos VaV \ni a, Wf(a)W \ni f(a) tais que f ⁣:VWf \colon V \to W é uma bijeção com inversa C1\mathcal C^1, e

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

Se ff é Ck\mathcal C^k, então f1f^{-1} também é.

Demonstração. Normalize: substituindo ff por xDf(a)1(f(a+x)f(a))x \mapsto Df(a)^{-1}\bigl(f(a + x) - f(a)\bigr), podemos supor a=0a = 0, f(0)=0f(0) = 0, Df(0)=IDf(0) = I (o enunciado geral decorre compondo com as bijeções afins). Escreva f(x)=x+g(x)f(x) = x + g(x): Dg(0)=0Dg(0) = 0 e, pela continuidade de DgDg, escolha r>0r > 0 com Dg(x)12\vertiii{Dg(x)} \leq \frac12 em Bˉ(0,r)\bar B(0, r); a desigualdade do valor médio dá g(x)g(x)12xx\norm{g(x) - g(x')} \leq \frac12 \norm{x - x'} ali.

Bijetividade sobre uma vizinhança. Para yB(0,r2)y \in B(0, \frac r2), resolver f(x)=yf(x) = y significa encontrar um ponto fixo de Φy(x)=yg(x)\Phi_y(x) = y - g(x); Φy\Phi_y leva Bˉ(0,r)\bar B(0, r) nele mesmo (Φy(x)y+12xr\norm{\Phi_y(x)} \leq \norm y + \frac12\norm x \leq r) e é 12\frac12-lipschitziana: Banach (Teorema 7.4) dá uma única solução x=φ(y)Bˉ(0,r)x = \varphi(y) \in \bar B(0,r). Além disso, ff é injetora em 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{$*$}

Ponha W=B(0,r2)W = B(0, \frac r2) e V=f1(W)B(0,r)V = f^{-1}(W)\cap B(0, r): abertos (continuidade), com f ⁣:VWf \colon V \to W bijetora.

Continuidade e derivabilidade da inversa. ()(*) diz que φ=f1\varphi = f^{-1} é 22-lipschitziana. Fixe y0=f(x0)Wy_0 = f(x_0) \in W; a inversibilidade de A=Df(x0)A = Df(x_0) (sua distância a II é 12\leq \frac12: Neumann, Proposição 8.4) e a derivabilidade de ff dão, para y=f(x)y = f(x) perto de 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}),

usando ()(*) para converter xx02yy0\norm{x - x_0} \leq 2\norm{y - y_0}: φ\varphi é derivável em y0y_0 com a diferencial inversa. Continuidade de yDφ(y)=Df(φ(y))1y \mapsto D\varphi(y) = Df(\varphi(y))^{-1}: composição de aplicações contínuas (a inversão é contínua, Proposição 8.4): φC1\varphi \in \mathcal C^1; e repetir a mesma fórmula (bootstrapping) dá Ck\mathcal C^k.

Teorema 20.2 (Teorema da função implícita)

Seja F ⁣:URp×RqRqF \colon U \subseteq \R^p\times\R^q \to \R^q de classe C1\mathcal C^1 perto de (a,b)(a, b), com F(a,b)=0F(a,b) = 0, e suponha que a diferencial parcial DyF(a,b)L(Rq)D_yF(a,b) \in \mathcal L(\R^q) seja inversível. Então existem vizinhanças AaA \ni a, BbB \ni b e uma aplicação ψ ⁣:AB\psi \colon A \to B de classe C1\mathcal C^1 com

{(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\},

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

Demonstração. Aplique o Teorema 20.1 a Θ(x,y)=(x,F(x,y))\Theta(x, y) = (x, F(x,y)): sua diferencial em (a,b)(a,b), triangular por blocos com blocos diagonais inversíveis II e DyFD_yF, é inversível. A inversa local tem a forma Θ1(x,z)=(x,h(x,z))\Theta^{-1}(x, z) = (x, h(x, z)); ponha ψ(x)=h(x,0)\psi(x) = h(x, 0): então F(x,y)=0F(x, y) = 0 se, e somente se, Θ(x,y)=(x,0)\Theta(x,y) = (x, 0), se, e somente se, y=ψ(x)y = \psi(x), localmente. A fórmula: derive F(x,ψ(x))=0F(x, \psi(x)) = 0 pela regra da cadeia.

20.2 Subvariedades: quatro definições

Teorema 20.3 (Caracterizações equivalentes)

Sejam MRnM \subseteq \R^n, d{0,,n}d \in \{0, \dots, n\} e k1k \geq 1. São equivalentes, em cada ponto aMa \in M (e MM é uma subvariedade de dimensão dd e classe Ck\mathcal C^k se valem em todo aMa \in M):

  1. (Retificação) Existe um difeomorfismo Φ\Phi Ck\mathcal C^k de um aberto Ωa\Omega \ni a sobre um aberto ΩRn\Omega' \subseteq \R^n com

    Φ(MΩ)=Ω(Rd×{0}).\Phi(M\cap\Omega) = \Omega' \cap \bigl(\R^d\times\{0\}\bigr).
  2. (Conjunto de nível) Existe uma submersão F ⁣:ΩRndF \colon \Omega \to \R^{n-d} de classe Ck\mathcal C^k (isto é, com DF(x)DF(x) sobrejetora) em um aberto Ωa\Omega \ni a com MΩ=F1(0)M\cap\Omega = F^{-1}(0).
  3. (Gráfico) A menos de permutar coordenadas, MM é localmente o gráfico de uma aplicação ψ ⁣:ARdRnd\psi \colon A \subseteq \R^d \to \R^{n-d} de classe Ck\mathcal C^k.
  4. (Parametrização) Existe uma imersão φ ⁣:ARdRn\varphi \colon A \subseteq \R^d \to \R^n de classe Ck\mathcal C^k (com Dφ(u)D\varphi(u) injetora) com AA aberto, sendo φ\varphi um homeomorfismo de AA sobre MΩM \cap \Omega para algum aberto Ωa\Omega \ni a.

Demonstração. (1)\Rightarrow(2): F=(Φd+1,,Φn)F = (\Phi_{d+1}, \dots, \Phi_n) (as últimas coordenadas de Φ\Phi): uma submersão (DΦD\Phi é inversível). (2)\Rightarrow(3): DF(a)DF(a) é sobrejetora: algum menor q×qq \times q da jacobiana é inversível (q=ndq = n - d); após permutar coordenadas, DyF(a)D_yF(a) é inversível, e o teorema da função implícita (Teorema 20.2) exprime MM localmente como gráfico y=ψ(x)y = \psi(x). (3)\Rightarrow(4): φ(x)=(x,ψ(x))\varphi(x) = (x, \psi(x)): uma imersão (diferencial (IDψ)\bigl(\begin{smallmatrix}I\\ D\psi\end{smallmatrix}\bigr) injetora), um homeomorfismo sobre o gráfico (inversa: a projeção, contínua). (4)\Rightarrow(1): seja φ(u0)=a\varphi(u_0) = a; complete imDφ(u0)\operatorname{im}D\varphi(u_0) por um suplementar EE (dimE=nd\dim E = n - d) e defina Θ(u,v)=φ(u)+v\Theta(u, v) = \varphi(u) + v em A×EA\times E: DΘ(u0,0)D\Theta(u_0, 0) é bijetora (a imagem contém imDφ(u0)\operatorname{im}D\varphi(u_0) e EE), de modo que Θ\Theta é um difeomorfismo local (Teorema 20.1); sua inversa Φ\Phi retifica: perto de aa, os pontos de MM são exatamente os φ(u)=Θ(u,0)\varphi(u) = \Theta(u, 0) — e, para isso, a hipótese de homeomorfismo em (4) garante que MΩM\cap\Omega, para Ω\Omega pequeno, não contém outras folhas (φ(u)+v=mM\varphi(u') + v = m \in M perto de aa com v0v \neq 0 pequeno deve ser excluído: m=φ(u)m = \varphi(u'') para algum uu'' perto de u0u_0 pela propriedade de homeomorfismo, e a injetividade local de Θ\Theta força v=0v = 0). Então Φ(MΩ)=(A×{0})Φ(Ω)\Phi(M\cap\Omega) = (A\times\{0\}) \cap\Phi(\Omega), a menos de encolher.

Exemplo 20.4

A esfera Sn1={x2=1}S^{n-1} = \{\norm x^2 = 1\}: conjunto de nível da submersão F(x)=x221F(x) = \norm x_2^2 - 1 em Rn{0}\R^n\setminus\{0\} (DF(x)=2xT0DF(x) = 2x^{\mathsf T} \ne 0): uma subvariedade C\mathcal C^\infty de dimensão n1n - 1. O toro em R3\R^3: conjunto de nível de (x2+y2R)2+z2r2\bigl(\sqrt{x^2 + y^2} - R\bigr)^2 + z^2 - r^2 (0<r<R0 < r < R). O cone {x2+y2=z2}\{x^2 + y^2 = z^2\} não é uma subvariedade em 00 (Exercício 20.1). Grupos de matrizes: SLnSL_n e OnO_n são subvariedades de Mn(R)M_n(\R) (Exercícios 20.5 e 20.6) — o ponto de partida da teoria de Lie.

20.3 Espaços tangentes

Definição 20.5

Sejam MM uma dd-subvariedade e aMa \in M. O espaço tangente TaMT_aM é o conjunto dos vetores velocidade γ(0)\gamma'(0) das curvas γ ⁣:(ε,ε)M\gamma \colon \intoo{-\varepsilon}\varepsilon \to M de classe C1\mathcal C^1 com γ(0)=a\gamma(0) = a.

Proposição 20.6

TaMT_aM é um subespaço vetorial de Rn\R^n de dimensão dd, e:

  1. se M=F1(0)M = F^{-1}(0) localmente com FF submersão: TaM=kerDF(a)T_aM = \ker DF(a);
  2. se MM é parametrizada pela imersão φ\varphi (φ(u0)=a\varphi(u_0) = a): TaM=imDφ(u0)T_aM = \operatorname{im}D\varphi(u_0).

Demonstração. As curvas em MM satisfazem F(γ(t))=0F(\gamma(t)) = 0; a regra da cadeia em 00DF(a)γ(0)=0DF(a)\gamma'(0) = 0: TaMkerDF(a)T_aM \subseteq \ker DF(a). Reciprocamente, a retificação (Teorema 20.3(1)) transporta retas de Rd×{0}\R^d\times\{0\} em curvas de MM: todo vetor de um subespaço de dimensão dd é realizado; comparar dimensões (dimkerDF(a)=n(nd)=d\dim\ker DF(a) = n - (n - d) = d) força a igualdade em (1), e o mesmo argumento de transporte dá (2) (Dφ(u0)D\varphi(u_0) aplicada a retas em AA; dimensões, de novo).

Teorema 20.7 (Multiplicadores de Lagrange)

Sejam M=F1(0)M = F^{-1}(0) com F=(F1,,Fq) ⁣:ΩRqF = (F_1, \dots, F_q) \colon \Omega \to \R^q uma submersão C1\mathcal C^1, e seja f ⁣:ΩRf \colon \Omega \to \R de classe C1\mathcal C^1. Se a restrição fMf\restriction_M tem extremo local em aMa \in M, então existem reais únicos λ1,,λq\lambda_1, \dots, \lambda_q (os multiplicadores de Lagrange) com

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

Demonstração. Para toda curva γ\gamma em MM por aa: tf(γ(t))t \mapsto f(\gamma(t)) tem extremo local em 00, de modo que 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) (Proposição 20.6). Ora, kerDF(a)=imDF(a)T\ker DF(a)^\perp = \operatorname{im}DF(a)^{\mathsf T} (a identidade de posto/ortogonalidade do segundo ano, ou o Exercício 13.8 em dimensão finita: (kerT)=imT(\ker T)^\perp = \operatorname{im}T^*), que é gerado pelos gradientes Fi(a)\nabla F_i(a) — independentes, pois DF(a)DF(a) é sobrejetora: os multiplicadores existem e são únicos.

Exemplo 20.8 (O teorema espectral, geometricamente)

Seja AA uma matriz real simétrica n×nn\times n e maximize f(x)=Ax,xf(x) = \langle Ax, x\rangle na esfera Sn1S^{n-1} (compacta: o máximo é atingido, em algum v1v_1). Lagrange com F(x)=x21F(x) = \norm x^2 - 1: f=2Ax\nabla f = 2Ax e F=2x\nabla F = 2x dão Av1=λ1v1Av_1 = \lambda_1v_1 — um autovetor, com λ1=maxSn1Ax,x\lambda_1 = \max_{S^{n-1}}\langle Ax, x\rangle. Restrinja AA a v1v_1^\perp (invariante: 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) e itere: uma base ortonormal de autovetores. O teorema espectral do segundo ano, redemonstrado por pura otimização — e a sombra em dimensão infinita do mesmo argumento demonstrou o Lema 15.6.

Método 20.9

Para demonstrar que um conjunto é subvariedade: exiba-o localmente como F1(0)F^{-1}(0) com DFDF sobrejetora no conjunto (a rota mais comum), ou como gráfico. Para calcular sua dimensão e seu espaço tangente: d=nqd = n - q e Ta=kerDF(a)T_a = \ker DF(a). Para otimizar nele: Lagrange — verifique sempre antes a compacidade (ou a coercividade), de modo que exista um extremo ao qual o teorema possa se aplicar, e lembre que a equação dos multiplicadores é apenas necessária: reúna todos os pontos críticos e, depois, compare valores. Para grupos de matrizes, derive curvas na identidade para identificar espaços tangentes.

20.4 Exercícios

Exercício 20.1

(a) Verifique que os conjuntos seguintes são subvariedades C\mathcal C^\infty e dê suas dimensões: Sn1S^{n-1}; o hiperboloide {x2+y2z2=1}\{x^2 + y^2 - z^2 = 1\}; o toro do Exemplo 20.4. (b) Mostre que o cone C={x2+y2=z2}R3C = \{x^2 + y^2 = z^2\} \subseteq \R^3 não é uma 22-subvariedade em 00: determine o número de componentes conexas de (C{0})B(0,ε)\bigl(C\setminus\{0\}\bigr)\cap B(0,\varepsilon) e compare com o número que uma retificação (Teorema 20.3(1)) forçaria para um plano menos um ponto.

Solução

Solução de Exercício 20.1.

(a) Cada um é F1(0)F^{-1}(0) para uma submersão: x21\norm x^2 - 1 em Rn{0}\R^n\setminus\{0\} (gradiente 2x02x \neq 0): dimensão n1n-1; x2+y2z21x^2 + y^2 - z^2 - 1 (gradiente (2x,2y,2z)0(2x, 2y, -2z) \neq 0 no hiperboloide, em que x2+y2=1+z2>0x^2 + y^2 = 1 + z^2 > 0): dimensão 22; a função do toro G=(ρR)2+z2r2G = (\rho - R)^2 + z^2 - r^2, ρ=x2+y2\rho = \sqrt{x^2+y^2}, é C\mathcal C^\infty perto do toro (ρRr>0\rho \geq R - r > 0 ali) com G0\nabla G \neq 0 (sua componente zz é 2z2z e, onde z=0z = 0, a componente radial é 2(ρR)02(\rho - R)\ne0, pois ρR=r\abs{\rho - R} = r): dimensão 22.

(b) Para ε\varepsilon pequeno, (C{0})B(0,ε)(C\setminus\{0\})\cap B(0,\varepsilon) tem exatamente 22 componentes conexas (as folhas superior e inferior perfuradas, cada uma conexa por caminhos: ligue por círculos e semirretas). Se CC fosse uma 22-subvariedade em 00, um endireitamento daria um homeomorfismo de CΩC\cap\Omega sobre uma peça aberta de um plano, levando 00 num ponto pp; pequenas vizinhanças planas perfuradas de pp têm uma componente, e os homeomorfismos preservam o número de componentes de vizinhanças perfuradas: contradição.

Exercício 20.2

Calcule os espaços tangentes: (a) TaSn1T_aS^{n-1} para aa qualquer (resposta: aa^\perp); (b) o plano tangente ao toro do Exemplo 20.4 em um ponto arbitrário do equador externo {z=0, x2+y2=(R+r)2}\{z = 0,\ x^2 + y^2 = (R + r)^2\}; (c) a reta tangente à hélice φ(t)=(cost,sint,t)\varphi(t) = (\cos t, \sin t, t) em φ(t0)\varphi(t_0), verificando a Proposição 20.6(2).

Solução

Solução de Exercício 20.2.

(a) TaSn1=ker(2aT)=aT_aS^{n-1} = \ker\bigl(2a^{\mathsf T}\bigr) = a^\perp. (b) Em 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), de modo que o plano tangente é Vect((sinθ,cosθ,0), (0,0,1))\operatorname{Vect}\bigl((-\sin\theta, \cos\theta, 0),\ (0, 0, 1)\bigr): o plano vertical tangente ao equador externo. (c) A hélice é uma curva mergulhada com φ(t0)=(sint0,cost0,1)0\varphi'(t_0) = (-\sin t_0, \cos t_0, 1) \neq 0: a reta tangente em φ(t0)\varphi(t_0) é φ(t0)+Rφ(t0)\varphi(t_0) + \R\,\varphi'(t_0), como a Proposição 20.6(2) prescreve.

Exercício 20.3 ★★

Seja f(x,y)=(x2y2, 2xy)f(x, y) = (x^2 - y^2,\ 2xy) (isto é, zz2z \mapsto z^2). (a) Em que pontos o Teorema 20.1 se aplica? (b) Mostre que ff é local mas não globalmente inversível em R2{0}\R^2\setminus\{0\}, e exiba explicitamente as duas inversas locais definidas em uma vizinhança de (1,0)(1, 0) (os dois ramos da raiz quadrada). (c) Mesma discussão para a aplicação de coordenadas polares (r,θ)(rcosθ,rsinθ)(r, \theta) \mapsto (r\cos\theta, r\sin\theta).

Solução

Solução de Exercício 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): o teorema se aplica em todo ponto, exceto na origem. (b) f(z)=f(z)f(-z) = f(z): nunca injetora num conjunto simétrico em relação a 00; em R2{0}\R^2\setminus\{0\} ela é um difeomorfismo local em toda parte e, ainda assim, 22 para 11 globalmente. Perto de (1,0)=f(±(1,0))(1, 0) = f(\pm(1, 0)), as duas inversas são os dois ramos da raiz quadrada: em notação complexa, w±ww \mapsto \pm\sqrt w (ramo principal), isto é,

(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) Jacobiano r>0r > 0: difeomorfismo local em (0,)×R\intoo0\infty\times\R, mas θθ+2π\theta \mapsto \theta + 2\pi dá o mesmo ponto: localmente inversível (ângulo determinado a menos de 2π2\pi num semiplano), nunca globalmente.

Exercício 20.4 ★★

(Fólio) Sejam F(x,y)=x3+y33xyF(x,y) = x^3 + y^3 - 3xy e C=F1(0)\mathcal C = F^{-1}(0). (a) Mostre que, perto de todo ponto de C\mathcal C distinto da origem, C\mathcal C é uma 11-subvariedade, localmente um gráfico em xx ou em yy (qual, e onde?). (b) Calcule a reta tangente em (32,32)(\frac32, \frac32). (c) O que acontece na origem? (Dois ramos se cruzam: exiba duas curvas C1\mathcal C^1 em C\mathcal C por 00 com velocidades independentes e conclua que não existe retificação.)

Solução

Solução de Exercício 20.4.

(a) F=3(x2y, y2x)\nabla F = 3(x^2 - y,\ y^2 - x) se anula se, e somente se, y=x2y = x^2 e x=y2x = y^2, isto é, x4=xx^4 = x: em (0,0)(0,0) e (1,1)(1,1); só (0,0)(0,0) está em C\mathcal C (F(1,1)=1F(1,1) = -1). Logo, em C{0}\mathcal C\setminus\{0\}, FF é uma submersão: uma 11-subvariedade, localmente um gráfico y=ψ(x)y = \psi(x) onde Fy=3(y2x)0F_y = 3(y^2 - x) \neq 0 e x=χ(y)x = \chi(y) onde Fx=3(x2y)0F_x = 3(x^2 - y) \neq 0 (ao menos uma das condições vale fora da origem).

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

(c) A parametrização racional x=3t1+t3x = \frac{3t}{1 + t^3}, y=3t21+t3y = \frac{3t^2}{1+t^3} passa por 00 em t=0t = 0 com velocidade (3,0)(3, 0); trocando xyx \leftrightarrow y (a curva é simétrica, ou reparametrize por 1/t1/t), obtém-se uma segunda curva C1\mathcal C^1 por 00 com velocidade (0,3)(0, 3). Duas direções tangentes independentes são impossíveis para uma 11-subvariedade (seu espaço tangente é uma reta, Proposição 20.6): C\mathcal C não é subvariedade na origem — um autocruzamento transversal.

Exercício 20.5 ★★

Seja F(M)=MTMF(M) = M^{\mathsf T}M de Mn(R)M_n(\R) no espaço SnS_n das matrizes simétricas. (a) Mostre que DF(M)(H)=MTH+HTMDF(M)(H) = M^{\mathsf T}H + H^{\mathsf T}M e que DF(M)DF(M) é sobrejetora sobre SnS_n em todo MOnM \in O_n (dado SSnS \in S_n, tente H=12MSH = \frac12MS). (b) Conclua que On=F1(I)O_n = F^{-1}(I) é uma subvariedade C\mathcal C^\infty compacta de dimensão n(n1)2\frac{n(n-1)}2, com TIOn={H:HT=H}T_IO_n = \{H : H^{\mathsf T} = -H\}, as matrizes antissimétricas. (c) Mostre que etHOn\eu^{tH} \in O_n para toda HH antissimétrica: as direções tangentes se integram em curvas dentro do grupo.

Solução

Solução de Exercício 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. Para MOnM \in O_n e SS simétrica, H=12MSH = \frac12MSDF(M)(H)=12(S+ST)=SDF(M)(H) = \frac12(S + S^{\mathsf T}) = S: sobrejetora sobre SnS_n.

(b) On=F1(I)O_n = F^{-1}(I) com FF submersão (sobre SnS_n, de dimensão n(n+1)2\frac{n(n+1)}2) em cada um de seus pontos: uma subvariedade de dimensão n2n(n+1)2=n(n1)2n^2 - \frac{n(n+1)}2 = \frac{n(n-1)}2. Compacto: fechado (FF contínua) e limitado (as colunas são vetores unitários). Tangente em 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 (a transposição atravessa a série; as exponenciais de matrizes que comutam se multiplicam, Teorema 19.8).

Exercício 20.6 ★★

(a) Mostre que det ⁣:Mn(R)R\det \colon M_n(\R) \to \R tem diferencial Ddet(M)(H)=tr(com(M)TH)D\det(M)(H) = \operatorname{tr}\bigl(\operatorname{com}(M) ^{\mathsf T}H\bigr), não nula em todo MSLnM \in SL_n. (b) Conclua que SLn(R)SL_n(\R) é uma subvariedade de dimensão n21n^2 - 1 com TISLn={H:trH=0}T_ISL_n = \{H : \operatorname{tr}H = 0\}. (c) GLn(R)GL_n(\R) é uma subvariedade? De que dimensão?

Solução

Solução de Exercício 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) para MM inversível (expansão de det\det perto de II: o termo linear de (1+λi)\prod(1 + \lambda_i)); com 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), e a fórmula se estende a todo MM por densidade e continuidade. Em SLnSL_n, detM=1\det M = 1: Ddet(M)0D\det(M) \ne 0 (seu valor em H=MH = M é tr(I)=ndetM=n\operatorname{tr}(I)\cdot\dots = n\det M = n).

(b) SLn=det1(1)SL_n = \det^{-1}(1) com det\det submersão ali (valores em R\R): dimensão n21n^2 - 1; TISLn=kerDdet(I)={H:trH=0}T_ISL_n = \ker D\det(I) = \{H : \operatorname{tr}H = 0\}.

(c) GLnGL_n é um subconjunto aberto de Mn(R)M_n(\R) (Exercício 6.8): uma subvariedade de dimensão total n2n^2 (endireitamento: a carta identidade).

Exercício 20.7 ★★

Por multiplicadores de Lagrange: (a) encontre os extremos de f(x,y)=xyf(x,y) = xy no círculo x2+y2=1x^2 + y^2 = 1; (b) mostre que, entre todos os vetores de probabilidade (p1,,pn)(p_1, \dots, p_n) (positivos e de soma 11), a entropia pilnpi-\sum p_i\ln p_i é maximizada exatamente na distribuição uniforme; (c) encontre o ponto da elipse {x2/4+y2=1}\{x^2/4 + y^2 = 1\} mais próximo de (1,0)(1, 0) e verifique geometricamente a equação dos multiplicadores (alinhamento das normais).

Solução

Solução de Exercício 20.7.

(a) (y,x)=λ(2x,2y)(y, x) = \lambda(2x, 2y) e x2+y2=1x^2 + y^2 = 1: y=2λxy = 2\lambda x, x=2λyx = 2\lambda y dão x2=y2=12x^2 = y^2 = \frac12. Valores de xyxy: ±12\pm\frac12: máximo 12\frac12 em ±12(1,1)\pm\frac1{\sqrt2}(1,1), mínimo 12-\frac12 em ±12(1,1)\pm\frac1{\sqrt2}(1,-1) (o conjunto de restrição é compacto: os extremos existem).

(b) No interior do simplexo (pi>0p_i > 0), Lagrange para H(p)=pilnpiH(p) = -\sum p_i\ln p_i com a restrição pi=1\sum p_i = 1: lnpi1=λ-\ln p_i - 1 = \lambda para todo ii: todos os pip_i iguais, pi=1np_i = \frac1n, com H=lnnH = \ln n. O máximo no simplexo compacto é atingido; se fosse atingido no bordo (algum pi=0p_i = 0), a distribuição viveria em n1\leq n - 1 pontos e, por indução, Hln(n1)<lnnH \leq \ln(n-1) < \ln n: o ponto crítico interior é o máximo global — a ignorância uniforme maximiza a entropia.

(c) Minimize (x1)2+y2(x-1)^2 + y^2 na elipse compacta: (2(x1),2y)=λ(x2,2y)(2(x{-}1), 2y) = \lambda(\frac x2, 2y). Se y0y \neq 0: λ=1\lambda = 1, e então 2(x1)=x22(x - 1) = \frac x2x=43x = \frac43, y2=149=59y^2 = 1 - \frac49 = \frac59: distância2^2 =19+59=23= \frac19 + \frac59 = \frac23. Se y=0y = 0: x=±2x = \pm2, distâncias 11 e 33. Pontos mais próximos: (43,±53)\bigl(\frac43, \pm\frac{\sqrt5}3\bigr), à distância 2/3<1\sqrt{2/3} < 1. A equação do multiplicador diz que o segmento de (1,0)(1,0) ao ponto mais próximo é paralelo a \nabla(elipse): ele encontra a elipse ortogonalmente, como a geometria exige.

Exercício 20.8 ★★★

Escreva o Exemplo 20.8 por extenso: demonstre por indução que uma matriz real simétrica admite base ortonormal de autovetores, sendo λ1λn\lambda_1 \geq \dots \geq \lambda_n os máximos com restrição sucessivos do quociente de Rayleigh. Deduza então as fórmulas de Courant–Fischer do Exercício 15.8 em dimensão finita diretamente dessa construção.

Solução

Solução de Exercício 20.8.

Indução sobre nn; n=1n = 1 trivial. A função de Rayleigh f(x)=Ax,xf(x) = \langle Ax, x\rangle atinge seu máximo λ1\lambda_1 no compacto Sn1S^{n-1} em algum v1v_1; Lagrange (Teorema 20.7, a esfera como conjunto de nível) dá 2Av1=2λv12Av_1 = 2\lambda v_1, e λ=Av1,v1=λ1\lambda = \langle Av_1, v_1\rangle = \lambda_1. O hiperplano v1v_1^\perp é AA-invariante (simetria: Av,v1=v,Av1=0\langle Av, v_1\rangle = \langle v, Av_1\rangle = 0); a restrição é simétrica e a indução fornece uma base ortonormal de autovetores v2,,vnv_2, \dots, v_n de v1v_1^\perp com autovalores λ2λn\lambda_2 \geq \dots \geq \lambda_n, cada um o máximo de ff na esfera do complemento ortogonal restante. Courant–Fischer segue exatamente como no Exercício 15.8: expanda x=civix = \sum c_iv_i; num espaço de teste de dimensão kk, intersecte com Vect(vk,,vn)\operatorname{Vect}(v_k, \dots, v_n) (contagem de dimensões em Rn\R^n) para obter minλk\min \leq \lambda_k, e Vect(v1,,vk)\operatorname{Vect}(v_1, \dots, v_k) atinge min=λk\min = \lambda_k.

Exercício 20.9 ★★★

(Desigualdade de Hadamard) Para MGLn(R)M \in GL_n(\R) com colunas c1,,cnc_1, \dots, c_n:

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

com igualdade se, e somente se, as colunas são ortogonais. (Reduza a colunas de norma 11 por reescalamento; maximize det\det no produto compacto de esferas (Sn1)n(S^{n-1})^n; em um maximizador, Lagrange em cada coluna separadamente dá cidet=λici\nabla_{c_i}\det = \lambda_ic_i, e cidet\nabla_{c_i}\det é a ii-ésima coluna de com(M)\operatorname{com}(M): deduza que MTMM^{\mathsf T}M é diagonal, logo =I= I, logo det=±1\det = \pm1.) Leitura geométrica: o volume de um paralelepípedo é no máximo o produto dos comprimentos de suas arestas.

Solução

Solução de Exercício 20.9.

Reescalar cada coluna para norma unitária divide det\abs{\det} por ci\prod\norm{c_i}: basta demonstrar detM1\abs{\det M} \leq 1 quando todas as colunas são unitárias, com igualdade se, e somente se, MOnM \in O_n. A função det\det é contínua no compacto (Sn1)n(S^{n-1})^n: ela atinge um máximo mdetI=1>0m \geq \det I = 1 > 0 em alguma MM. Fixando todas as colunas exceto a ii-ésima, det\det é linear em cic_i, com gradiente a ii-ésima coluna de com(M)\operatorname{com}(M); Lagrange na ii-ésima esfera: com(M)i=λici\operatorname{com}(M)_{\cdot i} = \lambda_i c_i. A identidade MTcom(M)=det(M)IM^{\mathsf T}\operatorname{com}(M) = \det(M)\,I se lê cj,com(M)i=det(M)δij\langle c_j, \operatorname{com}(M)_{\cdot i}\rangle = \det(M)\,\delta_{ij}, isto é, λicj,ci=det(M)δij\lambda_i\langle c_j, c_i\rangle = \det(M)\delta_{ij}; tomando j=ij = i: λi=detM=m0\lambda_i = \det M = m \neq 0, e então jij \neq ici,cj=0\langle c_i, c_j\rangle = 0: as colunas são ortonormais, MOnM \in O_n, m=detM=1m = \abs{\det M} = 1. Logo detci\abs{\det} \leq \prod\norm{c_i} sempre, com igualdade exatamente para colunas ortogonais (reescale de volta): o volume de um paralelepípedo é máximo, para comprimentos de aresta dados, quando as arestas são perpendiculares.

Exercício 20.10 ★★

Perto de quais de seus pontos o círculo S1S^1 é um gráfico y=ψ(x)y = \psi(x)? E um gráfico x=χ(y)x = \chi(y)? Verifique explicitamente a caracterização por gráficos (Teorema 20.3(3)) em (1,0)(1, 0) e explique em uma frase por que alguma permutação de coordenadas é sempre suficiente, mas nenhuma única serve sempre.

Solução

Solução de Exercício 20.10.

y=±1x2y = \pm\sqrt{1 - x^2} serve perto de todo ponto com y0y \neq 0; x=±1y2x = \pm\sqrt{1 - y^2} perto de todo ponto com x0x \neq 0; em (1,0)(1, 0): o gráfico x=1y2x = \sqrt{1 - y^2} sobre y(1,1)y \in \intoo{-1}1, que é o Teorema 20.3(3) com as coordenadas trocadas. Alguma permutação sempre serve porque a reta tangente, sendo unidimensional, não pode ser simultaneamente vertical e horizontal — mas ela pode ser qualquer uma das duas, de modo que nenhuma escolha fixa da coordenada “dependente” serve em todo ponto.

Exercício 20.11 ★★

(O grupo ortogonal como subvariedade, quantitativamente) (a) Mostre que On={M:MTM=I}O_n = \{M : M^{\mathsf T}M = I\} é compacto: limitado (cada coluna é um vetor unitário, de modo que Mn\norm M \leq \sqrt n para a norma euclidiana de matrizes) e fechado. (b) Mostre que seu espaço tangente em II é o espaço das matrizes antissimétricas, de dimensão n(n1)2\frac{n(n-1)}2, e em um AOnA \in O_n geral: TAOn={AK:KT=K}T_AO_n = \{AK : K^{\mathsf T} = -K\}. (c) Deduza que a aplicação tAexp(tK)t \mapsto A\exp(tK) é, para cada KK antissimétrica, uma curva em OnO_n por AA com velocidade AKAK (verifique exp(tK)On\exp(tK) \in O_n usando exp(X)T=exp(XT)\exp(X)^{\mathsf T} = \exp(X^{\mathsf T}) e exp(X)exp(X)=I\exp(-X)\exp(X) = I): todo vetor tangente é realizado por uma curva explícita, sem precisar do teorema da função implícita.

Solução

Solução de Exercício 20.11.

(a) A aplicação que o define, F(M)=MTMIF(M) = M^{\mathsf T}M - I, é contínua: On=F1(0)O_n = F^{-1}(0) é fechado; as colunas de uma matriz ortogonal são vetores unitários, de modo que a norma euclidiana (de Frobenius) vale exatamente n\sqrt n: limitado. Compacto por Heine–Borel em Mn(R)Rn2M_n(\R) \cong \R^{n^2}.

(b) OnO_n é o conjunto de nível F=0F = 0 estudado no capítulo: DF(A)H=ATH+HTADF(A)H = A^{\mathsf T}H + H^{\mathsf T}A, sobrejetora sobre as matrizes simétricas em cada AOnA \in O_n (dada SS simétrica, tome H=12ASH = \frac12AS), de modo que OnO_n é uma subvariedade de dimensão n2n(n+1)2=n(n1)2n^2 - \frac{n(n+1)}2 = \frac{n(n-1)}2 com

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

em A=IA = I essas são as matrizes antissimétricas.

(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 (as duas matrizes ±tK\pm tK comutam, de modo que o produto das exponenciais é a exponencial da soma): exp(tK)On\exp(tK) \in O_n, e γ(t)=Aexp(tK)\gamma(t) = A\exp(tK) é uma curva em OnO_n com γ(0)=A\gamma(0) = A, γ(0)=AK\gamma'(0) = AK. Quando KK percorre as matrizes antissimétricas, AKAK varre TAOnT_AO_n: a exponencial realiza todo o espaço tangente por curvas explícitas — o atalho de grupos de Lie que o Problema 20.1 explora para SO(3)SO(3).

Exercício 20.12 ★★

(Pontos críticos da distância) Sejam MRnM \subseteq \R^n uma subvariedade e pMp \notin M. Mostre que, se x0Mx_0 \in M minimiza a distância a pp (um tal ponto existe quando MM é fechada e não vazia — por quê?), então

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

(derive tγ(t)p2t \mapsto \norm{\gamma(t) - p}^2 ao longo de curvas em MM). Deduza: o ponto mais próximo em uma esfera está na semirreta que passa pelo centro; e use a condição para calcular a distância de p=(2,0)p = (2, 0) à parábola y=x2y = x^2 (reduza a uma cúbica e resolva-a numericamente com três algarismos).

Solução

Solução de Exercício 20.12.

Existência: intersecte MM com uma bola fechada grande em torno de pp para obter um compacto não vazio; a distância, contínua, atinge ali seu mínimo, e os pontos fora da bola estão mais longe. Condição de primeira ordem: para uma curva γ\gamma em MM com γ(0)=x0\gamma(0) = x_0, a função h(t)=γ(t)p2h(t) = \norm{\gamma(t) - p}^2 é derivável com um mínimo em 00:

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

e γ(0)\gamma'(0) varre Tx0MT_{x_0}M: px0Tx0Mp - x_0 \perp T_{x_0}M. Esfera S(c,r)S(c, r): o espaço tangente em x0x_0 é (x0c)(x_0 - c)^\perp, de modo que px0x0cp - x_0 \parallel x_0 - c: x0x_0 está na reta por cc e pp, à distância rr de cc — o ponto do raio, como a geometria insiste. Parábola: em x0=(x,x2)x_0 = (x, x^2) a tangente é gerada por (1,2x)(1, 2x); a ortogonalidade a px0=(2x,x2)p - x_0 = (2 - x, -x^2) se lê

(2x)2x3=0,isto eˊ2x3+x2=0,(2 - x) - 2x^3 = 0, \qquad\text{isto é}\qquad 2x^3 + x - 2 = 0,

com única raiz real (x2x3+xx \mapsto 2x^3 + x é estritamente crescente) x0.835x \approx 0.835; então x0(0.835,0.698)x_0 \approx (0.835, 0.698) e d(p,M)=(20.835)2+0.69821.358d(p, M) = \sqrt{(2 - 0.835)^2 + 0.698^2} \approx 1.358.

20.5 Problema: SO(3)SO(3) e os quatérnios

Problema 20.1

Problema de fim de semana — rotações, o grupo S3S^3 e o recobrimento duplo

Os quatérnios H={t+xi+yj+zk}\mathbb H = \{t + x\mathrm i + y\mathrm j + z\mathrm k\} — a álgebra cujo grupo de unidades contém o Q8Q_8 do Problema 1.1 — parametrizam as rotações tridimensionais duas vezes: a aplicação “conjugar por um quatérnio unitário” é um morfismo sobrejetor S3SO(3)S^3 \to SO(3) de núcleo {±1}\{\pm1\}. Construímos tudo. Recorde/defina: a multiplicação é R\R-bilinear com i2=j2=k2=ijk=1\mathrm i^2 = \mathrm j^2 = \mathrm k^2 = \mathrm{ijk} = -1; o conjugado de q=t+xi+yj+zkq = t + x\mathrm i + y\mathrm j + z\mathrm k é 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.

Parte I — A álgebra H\mathbb H e o grupo S3S^3.

  1. Verifique que H\mathbb H é uma R\R-álgebra associativa de centro R\R, que pq=qˉpˉ\overline{pq} = \bar q\,\bar p e que N(pq)=N(p)N(q)N(pq) = N(p)N(q) (uma rota limpa: represente qq pela matriz complexa 2×22\times2 (αββˉαˉ)\bigl(\begin{smallmatrix}\alpha & \beta\\ -\bar\beta & \bar\alpha\end{smallmatrix}\bigr), q=α+βjq = \alpha + \beta\mathrm j, e use det\det).
  2. Deduza que todo q0q \neq 0 é inversível (q1=qˉ/N(q)q^{-1} = \bar q/N(q)): H\mathbb H é um corpo (não comutativo), e S3={N(q)=1}S^3 = \{N(q) = 1\} é um grupo — e uma 33-subvariedade compacta de R4\R^4 (Exemplo 20.4).

Parte II — O morfismo de rotação. Identifique R3\R^3 com os quatérnios puros P={xi+yj+zk}P = \{x\mathrm i + y\mathrm j + z\mathrm k\} e, para qS3q \in S^3, defina ρq(v)=qvqˉ\rho_q(v) = q\,v\,\bar q.

  1. Mostre que ρq\rho_q leva PP em PP (os quatérnios puros são aqueles com vˉ=v\bar v = -v), é R\R-linear, preserva a norma, e que ρ ⁣:qρq\rho \colon q \mapsto \rho_q é um morfismo de grupos S3O(3)S^3 \to O(3).
  2. Calcule o núcleo: ρq=id\rho_q = \mathrm{id} se, e somente se, qq comuta com i,j,k\mathrm i, \mathrm j, \mathrm k, se, e somente se, qRS3={±1}q \in \R\cap S^3 = \{\pm1\}.
  3. Escreva q=cosθ2+sinθ2uq = \cos\frac\theta2 + \sin\frac\theta2\,u com uPu \in P, N(u)=1N(u) = 1 (por que isso é sempre possível para qS3q \in S^3?). Mostre que ρq\rho_q fixa uu e, no plano uPu^\perp\cap P, age como a rotação de ângulo θ\theta (calcule ρq(w)\rho_q(w) para wuw \perp u usando uw=wuuw = -wu para unitários puros ortogonais — demonstre essa identidade a partir da tabela de multiplicação, ou de uw+wu=2u,wuw + wu = -2\langle u, w\rangle).
  4. Conclua: imρSO(3)\operatorname{im}\rho \subseteq SO(3) (cada ρq\rho_q é uma rotação de eixo e ângulo como calculados — determinante +1+1 por continuidade de qdetρqq \mapsto \det\rho_q no conexo S3S^3, ou diretamente), e ρ\rho é sobrejetor sobre SO(3)SO(3): toda rotação de R3\R^3 tem eixo (demonstre: uma matriz ortogonal real 3×33\times3 com det=1\det = 1 tem autovalor 11 — considere o polinômio característico) e, portanto, é algum ρq\rho_q. Em resumo:

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

Parte III — SO(3)SO(3) como subvariedade; Rodrigues.

  1. Mostre que SO(3)SO(3) é uma subvariedade compacta de dimensão 33 de M3(R)M_3(\R) com TISO(3)=T_ISO(3) = as matrizes antissimétricas (Exercício 20.5; a condição sobre o determinante seleciona uma reunião de componentes).
  2. Para a matriz antissimétrica AuA_u associada a uR3u \in \R^3 (Auv=uvA_uv = u\wedge v, o produto vetorial), demonstre a fórmula de 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)

    (a partir de Au3=AuA_u^3 = -A_u: separe a série exponencial segundo as potências de AuA_u), e identifique-a como a rotação de eixo uu e ângulo θ\theta. Deduza que exp\exp leva as matrizes antissimétricas sobre SO(3)SO(3).

  3. Relacione as duas parametrizações: mostre que tρq(t)t \mapsto \rho_{q(t)} com q(t)=cost2+sint2uq(t) = \cos\frac t2 + \sin\frac t2\,u é um grupo a um parâmetro de rotações cuja derivada em t=0t = 0 é AuA_u — as exponenciais quaterniônica e matricial contam a mesma história, a meia velocidade e a velocidade plena, respectivamente.

Parte IV — O recobrimento duplo, sentido.

  1. Mostre que o caminho q(t)=cost2+sint2kq(t) = \cos\frac t2 + \sin\frac t2\,\mathrm k, t[0,2π]t \in \intcc0{2\pi}, é um laço em SO(3)SO(3) (sua imagem ρq(t)\rho_{q(t)} volta à identidade) cujo levantamento quaterniônico não é um laço: q(2π)=q(0)q(2\pi) = -q(0). Continuar até t=4πt = 4\pi fecha o levantamento. Explique em um parágrafo curto o que isso diz: uma rotação de 2π2\pi não pode ser desfeita continuamente, mas uma de 4π4\pi pode (o truque do cinto), porque os laços de SO(3)SO(3) são detectados em seu recobrimento duplo S3S^3.
  2. Deduza também o dividendo prático: compor rotações = multiplicar quatérnios (o equivalente a 44 multiplicações de dados em vez de 99, sem deriva da ortogonalidade) — verifique na composição de duas rotações de um quarto de volta em torno de i\mathrm i e de j\mathrm j: calcule o eixo e o ângulo do produto.

Parte V — A matriz explícita: Euler–Rodrigues. Escreva q=a+bi+cj+dkS3q = a + b\mathrm i + c\mathrm j + d\mathrm k \in S^3, de modo que a2+b2+c2+d2=1a^2 + b^2 + c^2 + d^2 = 1.

  1. Calcule ρq(i)\rho_q(\mathrm i) por extenso a partir da tabela de multiplicação; obtenha em seguida ρq(j)\rho_q(\mathrm j) e ρq(k)\rho_q(\mathrm k) pela substituição cíclica 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) (justifique-a: ciclar i,j,k\mathrm i, \mathrm j, \mathrm k se estende a um automorfismo de H\mathbb H, pois as relações de definição são ciclicamente simétricas). Conclua que a matriz de ρq\rho_q na base (i,j,k)(\mathrm i, \mathrm j, \mathrm k) é a matriz de Euler–Rodrigues

    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. (Lendo uma rotação de trás para diante) Mostre que

    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,

    na notação das questões 5 e 8. Deduza um algoritmo que recupera ±q\pm q a partir de uma matriz de rotação RR: o ângulo pelo traço; o eixo pela parte antissimétrica quando 0<θ<π0 < \theta < \pi; e, quando θ=π\theta = \pi, demonstre e use a identidade R+I=2uuTR + I = 2\,uu^{\mathsf T}.

  3. Avalie RqR_q para o produto q=12(1+i+j+k)q = \frac12(1 + \mathrm i + \mathrm j + \mathrm k) da questão 11: aparece uma matriz de permutação. Identifique a rotação e reconcilie com o eixo e o ângulo encontrados na questão 11.

Parte VI — Dentro de S3S^3: SU(2)SU(2), classes de conjugação, exponenciais.

  1. Mostre que a representação matricial da questão 1 (chame-a Φ\Phi) se restringe a um isomorfismo de grupos de S3S^3 sobre o grupo especial unitário

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

    (para a sobrejetividade, escreva as equações U1=UU^{-1} = U^* e detU=1\det U = 1 para uma matriz complexa 2×22\times2 geral).

  2. Mostre que a parte real é um invariante de conjugação em S3S^3Re(pqpˉ)=Req\operatorname{Re}(pq\bar p) = \operatorname{Re}q para todo pS3p \in S^3 — e, reciprocamente, que dois quatérnios unitários de mesma parte real são conjugados em S3S^3 (reduza a levar um eixo unitário puro sobre outro, o que a Parte II fornece). Descreva geometricamente as classes de conjugação de S3S^3; traduza para SU(2)SU(2) (conjuntos de nível do traço); e projete por ρ\rho: duas rotações são conjugadas em SO(3)SO(3) se, e somente se, têm o mesmo ângulo θ[0,π]\theta \in \intcc0\pi.
  3. Defina exp\exp em H\mathbb H pela série exponencial; verifique a convergência absoluta, usando pq=pq\abs{pq} = \abs p\,\abs q para q=N(q)\abs q = \sqrt{N(q)}. Mostre, para uu puro unitário e θR\theta \in \R,

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

    deduza que exp\exp leva o hiperplano PP sobre S3S^3, e verifique que ρexp(su)=e2sAu\rho_{\exp(su)} = \eu^{2sA_u}: de novo o fenômeno do meio ângulo da questão 9.

  4. Para quatérnios puros v,wv, w, demonstre a regra do produto vw=v,w+vwvw = -\langle v, w\rangle + v\wedge w, e daí a identidade do comutador vwwv=2vwvw - wv = 2\,v\wedge w; demonstre também [Av,Aw]=Avw[A_v, A_w] = A_{v\wedge w} para as matrizes da questão 8. Conclua que a derivada de ρ\rho em 11 ao longo das curvas texp(tv)t \mapsto \exp(tv) é o isomorfismo linear v2Avv \mapsto 2A_v de PP sobre as matrizes antissimétricas, e que ele transporta o comutador de quatérnios ao comutador de matrizes.

Parte VII — Estrutura global.

  1. (Não há seção contínua) Suponha s ⁣:SO(3)S3s \colon SO(3) \to S^3 contínua com ρs=id\rho \circ s = \operatorname{id}. Para o laço R(t)=ρq(t)R(t) = \rho_{q(t)} da questão 10, ponha ε(t)=s(R(t))q(t)1\varepsilon(t) = s(R(t))\,q(t)^{-1} para t[0,2π]t \in \intcc0{2\pi}. Mostre que ε\varepsilon é contínua com valores em {±1}\{\pm1\} e obtenha uma contradição: não existe escolha global contínua de um quatérnio unitário que represente cada rotação.
  2. (O modelo da bola) Sejam BˉR3\bar B \subseteq \R^3 a bola fechada de raio π\pi e E(v)=eAvE(v) = \eu^{A_v}, com E(0)=IE(0) = I. Mostre que EE leva Bˉ\bar B sobre SO(3)SO(3), é injetora na bola aberta e, na esfera de fronteira, identifica exatamente os antípodas: E(πu)=E(πu)=2uuTIE(\pi u) = E(-\pi u) = 2uu^{\mathsf T} - I, sem nenhuma outra coincidência. Assim, SO(3)SO(3) é a bola com os pontos antípodas da fronteira colados — o espaço projetivo RP3\mathbb{RP}^3 — e um diâmetro se torna o laço não contrátil da questão 10.
  3. Mostre que ρpρqρp1=ρpqpˉ\rho_p\rho_q\rho_p^{-1} = \rho_{pq\bar p}; que as involuções de SO(3)SO(3) (os RIR \neq I com R2=IR^2 = I) são exatamente as meias-voltas ρw\rho_w com ww quatérnio puro unitário; e que o centro de SO(3)SO(3) é trivial.
  4. Mostre que toda rotação é produto de duas meias-voltas: para q=cosθ2+sinθ2uq = \cos\frac\theta2 + \sin\frac\theta2\,u, escolha um puro unitário wuw \perp u, verifique que w=qww' = qw é de novo um quatérnio puro unitário, e confirme ρq=ρwρw\rho_q = \rho_{w'}\rho_w. Onde ficam os dois eixos, e que ângulo eles formam?
  5. Conclua o resumo topológico: SO(3)SO(3) é compacto e conexo por caminhos (dê duas demonstrações: imagem contínua de S3S^3 por ρ\rho; imagem de exp\exp), ao passo que O(3)O(3) tem exatamente duas componentes conexas, cada uma homeomorfa a SO(3)SO(3).
  6. (Uma composição, de três maneiras) Sejam R1R_1 a rotação de π2\frac\pi2 em torno do eixo zz e R2R_2 a rotação de π2\frac\pi2 em torno do eixo xx. Calcule o eixo e o ângulo de R2R1R_2R_1: (i) multiplicando as duas matrizes 3×33\times3 e usando o traço e a parte antissimétrica (Parte V); (ii) multiplicando os quatérnios unitários correspondentes q2q1q_2q_1. Verifique que as duas respostas coincidem: ângulo 2π3\frac{2\pi}3, eixo 13(1,1,1)\frac1{\sqrt3}(1, -1, 1).
  7. (A transformada de Cayley) Para KK antissimétrica, mostre que I+KI + K é inversível e

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

    com 1-1 nunca sendo autovalor de C(K)C(K); mostre que KC(K)K \mapsto C(K) é uma bijeção das matrizes antissimétricas sobre {RSO(n):1SpR}\{R \in SO(n) : -1 \notin \operatorname{Sp}R\}, com inversa R(IR)(I+R)1R \mapsto (I - R)(I + R)^{-1}. (Uma carta racional de SO(n)SO(n), companheira da exp\exp transcendente do Exercício 20.11.)

Solução

Solução de Problema 20.1.

1. Aplique 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) com α=t+ix\alpha = t + \iu x, β=y+iz\beta = y + \iu z: verifica-se que 1,i,j,k1, \mathrm i, \mathrm j, \mathrm k vão para 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), cujos produtos reproduzem a tabela dos quatérnios: a aplicação é um morfismo injetor de álgebras, de modo que H\mathbb H herda a associatividade; N(q)=α2+β2=detN(q) = \abs\alpha^2 + \abs\beta^2 = \det é multiplicativa, e a conjugação corresponde à transposta da adjunta, dando pq=qˉpˉ\overline{pq} = \bar q\bar p. Centro: comutar com i\mathrm i força y=z=0y = z = 0, e com j\mathrm j força x=0x = 0: R\R.

2. qqˉ=N(q)q\bar q = N(q): para q0q \neq 0, q1=qˉ/N(q)q^{-1} = \bar q/N(q): uma álgebra de divisão. Em S3S^3: N(pq)=1N(pq) = 1 e N(q1)=1N(q^{-1}) = 1: um grupo; e S3R4S^3 \subseteq \R^4 é a esfera unitária: uma 33-subvariedade compacta.

3. vv é puro se, e somente se, vˉ=v\bar v = -v; então qvqˉ=qvˉqˉ=qvqˉ\overline{qv\bar q} = q\bar v\bar q = -qv\bar q: ρq\rho_q preserva PP. A linearidade é clara; N(qvqˉ)=N(q)N(v)N(q)=N(v)N(qv\bar q) = N(q)N(v)N(q) = N(v): uma isometria de (P,N)(R3,2)(P, N) \cong (\R^3, \norm\cdot^2): ρqO(3)\rho_q \in O(3). E ρ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)): um morfismo.

4. ρq=id\rho_q = \mathrm{id} se, e somente se, qv=vqqv = vq para todo vv puro, se, e somente se, qq comuta com i,j,k\mathrm i, \mathrm j, \mathrm k, se, e somente se, qq é central (questão 1): qRS3={±1}q \in \R\cap S^3 = \{\pm1\}.

5. Escreva q=t+pq = t + p (tRt \in \R, pp puro): 1=N(q)=t2+N(p)1 = N(q) = t^2 + N(p), de modo que t=cosθ2t = \cos\frac\theta2 e p=sinθ2up = \sin\frac\theta2\,u com N(u)=1N(u) = 1 para algum θ\theta (se p=0p = 0, q=±1q = \pm1 age trivialmente). Como u2=N(u)=1u^2 = -N(u) = -1, qq e uu comutam, e ρq(u)=quqˉ=uqqˉ=u\rho_q(u) = qu\bar q = uq\bar q = u: o eixo. Para unitários puros wuw \perp u: a regra do produto vw=v,w+vwvw = -\langle v, w\rangle + v\wedge w (expanda em coordenadas a partir da tabela) dá uw=uw=wuuw = u\wedge w = -wu. Então

ρ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),

usando uwu=u2w=wuwu = -u^2w = w e as fórmulas do arco duplo: a rotação de ângulo θ\theta no plano orientado (w,uw)(w, u\wedge w).

6. Cada ρq\rho_q é uma rotação em torno de uu de θ\theta: na base ortonormal (u,w,uw)(u, w, u\wedge w), sua matriz tem determinante +1+1: imρSO(3)\operatorname{im}\rho \subseteq SO(3). Sobrejetividade: uma matriz RSO(3)R \in SO(3) tem 11 como autovalor, já que

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),

de modo que det(RI)=0\det(R - I) = 0. Tome um autovetor unitário uu; RR preserva uu^\perp e ali se restringe a uma rotação de algum ângulo θ\theta (ortogonal planar, de determinante 11): R=ρqR = \rho_q para q=cosθ2+sinθ2uq = \cos\frac\theta2 + \sin\frac\theta2\,u. Com a questão 4 e o primeiro teorema de isomorfismo (Teorema 1.3): SO(3)S3/{±1}SO(3) \cong S^3/\{\pm1\}.

7. O3O_3 é uma subvariedade compacta de dimensão 33 (Exercício 20.5); det\det é contínuo nela, com valores em {±1}\{\pm1\}, de modo que SO(3)=O3{det=1}SO(3) = O_3\cap\{\det = 1\} é aberto e fechado em O3O_3: uma união de componentes conexas, logo ela própria uma 33-subvariedade compacta, com o mesmo espaço tangente em II: as matrizes antissimétricas.

8. Au2v=u(uv)=u,vuvA_u^2v = u\wedge(u\wedge v) = \langle u, v\rangle u - v (uu unitário), de modo que Au3v=u(u,vuv)=uvA_u^3v = u\wedge(\langle u,v\rangle u - v) = -u\wedge v: Au3=AuA_u^3 = -A_u. Separando a série exponencial pelos restos das potências módulo a relação 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 .

Em uu: Auu=0A_uu = 0: fixo. Em 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: a rotação de eixo uu e ângulo θ\theta — Rodrigues. Toda rotação tem essa forma (questão 6): exp\exp é sobrejetora sobre SO(3)SO(3) a partir das matrizes antissimétricas.

9. Com q(t)=cost2+sint2uq(t) = \cos\frac t2 + \sin\frac t2\,u: a questão 5 mostra que ρq(t)\rho_{q(t)} é a rotação de eixo uu e ângulo tt, isto é, ρq(t)=etAu\rho_{q(t)} = \eu^{tA_u}, cuja derivada em t=0t = 0 é AuA_u. O quatérnio corre à metade do ângulo — o traço analítico do recobrimento duplo.

10. ρq(t)\rho_{q(t)} é a rotação em torno de k\mathrm k de ângulo tt: em t=2πt = 2\pi ela retorna à identidade — um laço em SO(3)SO(3). Seu levantamento satisfaz q(2π)=cosπ=1=q(0)q(2\pi) = \cos\pi = -1 = -q(0): o caminho levantado não é fechado; só em t=4πt = 4\pi é que qq retorna a 11. Interpretação: o laço das rotações completas não é contrátil em SO(3)SO(3) — seu levantamento termina na outra folha do recobrimento —, ao passo que o laço duplo é; um corpo preso ao ambiente por tiras (o truque do cinto) retorna a um estado sem torção após 4π4\pi, mas não após 2π2\pi. Os grupos de rotação lembram a paridade das voltas completas; S3S^3, sendo simplesmente conexo, é onde essa memória vive.

11. Quartos de volta: 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. Produto (aplicando primeiro a volta de j\mathrm j):

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),

de norma 11, com cosθ2=12\cos\frac\theta2 = \frac12: θ=2π3\theta = \frac{2\pi}3, e eixo u=i+j+k3u = \frac{\mathrm i + \mathrm j + \mathrm k}{\sqrt3} (a parte pura normalizada). Dois quartos de volta sucessivos em torno de eixos ortogonais produzem uma rotação de 120120^\circ em torno da diagonal principal do cubo — o equivalente a quatro multiplicações reais de contabilidade, com a ortogonalidade preservada exatamente: eis por que os softwares de voo e os motores gráficos compõem rotações via quatérnios.

12. Da tabela, ji=k\mathrm{ji} = -\mathrm k e ki=j\mathrm{ki} = \mathrm j, de modo que

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 .

Multiplicando por qˉ=abicjdk\bar q = a - b\mathrm i - c\mathrm j - d\mathrm k com a regra escalar–vetor (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, em que p1=(a,d,c)p_1 = (a, d, -c) e p2=(b,c,d)p_2 = (-b, -c, -d): a parte escalar é ab(ab)=0-ab - (-ab) = 0 (pura, como deve ser), e a parte vetorial é

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):

a primeira coluna de RqR_q. A aplicação cíclica σ(i)=j\sigma(\mathrm i) = \mathrm j, σ(j)=k\sigma(\mathrm j) = \mathrm k, σ(k)=i\sigma(\mathrm k) = \mathrm i preserva as relações i2=j2=k2=ijk=1\mathrm i^2 = \mathrm j^2 = \mathrm k^2 = \mathrm{ijk} = -1 (a palavra ijk\mathrm{ijk} é ciclicamente invariante a menos da relação ijk=jki\mathrm{ijk} = \mathrm{jki}, que vale em qualquer anel: conjugue ijk=1\mathrm{ijk} = -1 pelo inversível i\mathrm i), de modo que σ\sigma se estende a um automorfismo de R\R-álgebras, e σ(ρq(v))=ρσ(q)(σ(v))\sigma(\rho_q(v)) = \rho_{\sigma(q)}(\sigma(v)). Desenrolando, a imagem de j\mathrm j é a fórmula da primeira coluna após a substituição (b,c,d)(c,d,b)(b, c, d) \to (c, d, b) com a base reetiquetada ijki\mathrm i \to \mathrm j \to \mathrm k \to \mathrm i, que é exatamente a segunda coluna exibida; mais uma volta dá a terceira.

13. Somando a diagonal, trRq=3a2(b2+c2+d2)=4a21\operatorname{tr}R_q = 3a^2 - (b^2 + c^2 + d^2) = 4a^2 - 1 (norma unitária) e, com a=cosθ2a = \cos\frac\theta2: 4cos2θ21=1+2cosθ4\cos^2\frac\theta2 - 1 = 1 + 2\cos\theta. Parte antissimétrica: as três entradas independentes de RqRqTR_q - R_q^{\mathsf T} são 4ab,4ac,4ad4ab, 4ac, 4ad (nas posições (3,2),(1,3),(2,1)(3,2), (1,3), (2,1)), de modo que 12(RqRqT)=Am\frac12(R_q - R_q^{\mathsf T}) = A_m com 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. Algoritmo: θ=arccostrR12[0,π]\theta = \arccos\frac{\operatorname{tr}R - 1}{2} \in \intcc0\pi; se 0<θ<π0 < \theta < \pi, leia uu em RRT2sinθ\frac{R - R^{\mathsf T}}{2\sin\theta} e ponha q=±(cosθ2+sinθ2u)q = \pm(\cos\frac\theta2 + \sin\frac\theta2 u); se θ=0\theta = 0, q=±1q = \pm1. Para θ=π\theta = \pi: a=0a = 0, e Rodrigues (questão 8) dá R=I+2Au2=I+2(uuTI)=2uuTIR = I + 2A_u^2 = I + 2(uu^{\mathsf T} - I) = 2uu^{\mathsf T} - I, isto é, R+I=2uuTR + I = 2uu^{\mathsf T}; qualquer coluna não nula de R+IR + I, normalizada, é ±u\pm u, e q=±uq = \pm u.

14. Com a=b=c=d=12a = b = c = d = \frac12: todas as entradas diagonais se anulam, 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},

a permutação cíclica e1e2e3e1e_1 \to e_2 \to e_3 \to e_1. Seu traço é 0=1+2cosθ0 = 1 + 2\cos\theta, de modo que θ=2π3\theta = \frac{2\pi}3, e ela fixa (1,1,1)(1,1,1): a rotação de 120120^\circ em torno da diagonal principal — precisamente a resposta da questão 11, agora visível como a matriz que cicla os eixos coordenados.

15. Na base, verifica-se Φ(qˉ)=Φ(q)\Phi(\bar q) = \Phi(q)^* (a matriz de qˉ\bar q tem α=αˉ\alpha' = \bar\alpha, β=β\beta' = -\beta, que é a transposta conjugada de (αββˉαˉ)\bigl(\begin{smallmatrix}\alpha & \beta\\ -\bar\beta & \bar\alpha\end{smallmatrix}\bigr)). Logo Φ(q)Φ(q)=Φ(qˉq)=N(q)I\Phi(q)^*\Phi(q) = \Phi(\bar qq) = N(q)I e detΦ(q)=α2+β2=N(q)\det\Phi(q) = \abs\alpha^2 + \abs\beta^2 = N(q): para qS3q \in S^3, Φ(q)SU(2)\Phi(q) \in SU(2), e Φ\Phi é um morfismo injetor (questão 1). Sobrejetividade: seja U=(αβγδ)U = \bigl(\begin{smallmatrix}\alpha & \beta\\ \gamma & \delta\end{smallmatrix}\bigr) com detU=1\det U = 1; então U1=(δβγα)U^{-1} = \bigl(\begin{smallmatrix}\delta & -\beta\\ -\gamma & \alpha\end{smallmatrix}\bigr), e U1=U=(αˉγˉβˉδˉ)U^{-1} = U^* = \bigl(\begin{smallmatrix}\bar\alpha & \bar\gamma\\ \bar\beta & \bar\delta\end{smallmatrix}\bigr) força δ=αˉ\delta = \bar\alpha, γ=βˉ\gamma = -\bar\beta, e então 1=detU=α2+β21 = \det U = \abs\alpha^2 + \abs\beta^2: U=Φ(q)U = \Phi(q) para o quatérnio unitário qq de coordenadas α=a+ib\alpha = a + \iu b, β=c+id\beta = c + \iu d. Logo S3SU(2)S^3 \cong SU(2).

16. Os escalares reais são centrais e N(p)=1N(p) = 1pqpˉ=pqˉpˉ\overline{pq\bar p} = p\bar q\bar p, de modo que pqpˉ+pqpˉ=p(q+qˉ)pˉ=q+qˉpq\bar p + \overline{pq\bar p} = p(q + \bar q)\bar p = q + \bar q: a parte real é invariante. Reciprocamente, seja Req=Req=a\operatorname{Re}q = \operatorname{Re}q' = a; então as partes puras têm a mesma norma 1a2=s\sqrt{1 - a^2} = s. Se s=0s = 0, q=q=±1q = q' = \pm1. Se s>0s > 0, escreva q=a+suq = a + su, q=a+suq' = a + su' com u,uu, u' puros unitários; a questão 6 fornece uma rotação que leva uu' em uu, isto é, pS3p \in S^3 com ρp(u)=u\rho_p(u') = u, e então pqpˉ=a+sρp(u)=qpq'\bar p = a + s\rho_p(u') = q. As classes de S3S^3 são, portanto, {1}\{1\}, {1}\{-1\} e, para cada a(1,1)a \in \intoo{-1}1, a 22-esfera {a+su:u puro unitaˊrio}\{a + su : u \text{ puro unitário}\} de raio ss. Sob Φ\Phi, trΦ(q)=α+αˉ=2Req\operatorname{tr}\Phi(q) = \alpha + \bar\alpha = 2\operatorname{Re}q: as classes de SU(2)SU(2) são os conjuntos de nível do traço. Projetando: se q=pqpˉq' = pq\bar p, então ρq=ρpρqρp1\rho_{q'} = \rho_p\rho_q\rho_p^{-1}; reciprocamente, ρq=ρpρqρp1=ρpqpˉ\rho_{q'} = \rho_p\rho_q\rho_p^{-1} = \rho_{pq\bar p} força q=±pqpˉq' = \pm pq\bar p (núcleo), de modo que Req=±Req\operatorname{Re}q' = \pm\operatorname{Re}q, isto é, cosθ2=Req=Req=cosθ2\cos\frac{\theta'}2 = \abs{\operatorname{Re}q'} = \abs{\operatorname{Re}q} = \cos\frac\theta2 para os ângulos em [0,π]\intcc0\pi: rotações conjugadas têm ângulos iguais. Reciprocamente, ângulos iguais permitem representantes com a mesma parte real não negativa, conjugados pelo que precede: em SO(3)SO(3), a classe de conjugação de uma rotação é exatamente o seu ângulo.

17. NN é multiplicativa, de modo que \abs\cdot é uma norma multiplicativa em HR4\mathbb H \cong \R^4 e qk=qk\abs{q^k} = \abs q^k: a série qk/k!\sum q^k/k! converge absolutamente no espaço de dimensão finita (logo completo), dominada por qk/k!=eq\sum\abs q^k/k! = \eu^{\abs q}. Para um uu puro unitário: u2=1u^2 = -1, de modo que (θu)2m=(1)mθ2m(\theta u)^{2m} = (-1)^m\theta^{2m} e (θu)2m+1=(1)mθ2m+1u(\theta u)^{2m+1} = (-1)^m\theta^{2m+1}u; separando a série,

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 .

Qualquer qS3q \in S^3 é cosα+sinαu\cos\alpha + \sin\alpha\,u com α[0,π]\alpha \in \intcc0\pi (questão 5): q=exp(αu)q = \exp(\alpha u), de modo que exp(P)=S3\exp(P) = S^3. Por fim, exp(su)=coss+sinsu=q(2s)\exp(su) = \cos s + \sin s\,u = q(2s) na notação da questão 9, e ρq(t)=etAu\rho_{q(t)} = \eu^{tA_u} ali: ρexp(su)=e2sAu\rho_{\exp(su)} = \eu^{2sA_u}.

18. Expandindo vwvw coordenada a coordenada com a tabela: os produtos ii=1\mathrm i\cdot\mathrm i = -1, … dão o escalar (v1w1+v2w2+v3w3)-(v_1w_1 + v_2w_2 + v_3w_3), e os produtos mistos (ij=k\mathrm{ij} = \mathrm k, ji=k\mathrm{ji} = -\mathrm k, …) dão o vetor (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. Subtraindo o produto invertido: vwwv=2vwvw - wv = 2\,v\wedge w (as partes escalares se cancelam, os produtos vetoriais se somam). Para as matrizes, com 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 .

Derivada: exp(tv)=exp(tv)\overline{\exp(tv)} = \exp(-tv) (a conjugação é contínua e nega os quatérnios puros), de modo que

 ⁣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 :

a diferencial é v2Avv \mapsto 2A_v, uma bijeção linear de PP sobre as matrizes antissimétricas, e [2Av,2Aw]=4Avw=2A2vw=2A[v,w][2A_v, 2A_w] = 4A_{v\wedge w} = 2A_{2v\wedge w} = 2A_{[v,w]} mostra que ela leva o comutador de quatérnios no comutador de matrizes.

19. Aplicando ρ\rho: ρ(ε(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}, de modo que ε(t)kerρ={±1}\varepsilon(t) \in \ker\rho = \{\pm1\} (questão 4). Como produto das aplicações contínuas ts(R(t))t \mapsto s(R(t)) e tq(t)1=q(t)ˉt \mapsto q(t)^{-1} = \bar{q(t)}, ε\varepsilon é contínua no conexo intervalo [0,2π]\intcc0{2\pi} com valores no par discreto {±1}\{\pm1\}: ela é constante, digamos ε(t)ε\varepsilon(t) \equiv \varepsilon. Mas R(0)=R(2π)=IR(0) = R(2\pi) = I, de modo que s(R(0))=s(R(2π))s(R(0)) = s(R(2\pi)), ao passo que s(R(0))=εq(0)=εs(R(0)) = \varepsilon\,q(0) = \varepsilon e s(R(2π))=εq(2π)=εs(R(2\pi)) = \varepsilon\,q(2\pi) = -\varepsilon: contradição. Não existe seção contínua: a ambiguidade de sinal ±q\pm q é global, não um defeito de uma fórmula particular.

20. Sobrejetora: todo RSO(3)R \in SO(3) é eθAu\eu^{\theta A_u} para algum uu unitário e θ[0,2π]\theta \in \intcc0{2\pi} (questões 6 e 8); se θ>π\theta > \pi, Rodrigues dá eθAu=e(2πθ)Au\eu^{\theta A_u} = \eu^{(2\pi - \theta)A_{-u}} (ambos iguais a I+sinθAu+(1cosθ)Au2I + \sin\theta A_u + (1 - \cos\theta)A_u^2, e Au=AuA_{-u} = -A_u com sin(2πθ)=sinθ\sin(2\pi - \theta) = -\sin\theta, cos(2πθ)=cosθ\cos(2\pi - \theta) = \cos\theta), de modo que R=E(v)R = E(v) com vπ\norm v \leq \pi. Injetora no interior: se E(v)=E(v)IE(v) = E(v') \neq I com v,v<π\norm v, \norm{v'} < \pi, a questão 13 recupera o mesmo ângulo θ=v=v(0,π)\theta = \norm v = \norm{v'} \in \intoo0\pi a partir do traço e, como sinθ0\sin\theta \neq 0, o mesmo eixo a partir da parte antissimétrica: v=vv = v'; e E(v)=IE(v) = I força θ{0}\theta \in \{0\} na bola aberta. Bordo: E(πu)=I+2Au2=2uuTIE(\pi u) = I + 2A_u^2 = 2uu^{\mathsf T} - I depende de uu apenas através de uuTuu^{\mathsf T}, donde E(πu)=E(πu)E(\pi u) = E(-\pi u); reciprocamente, 2uuTI=2uuTI2uu^{\mathsf T} - I = 2u'u'^{\mathsf T} - I aplicado a uuu=u,uuu = \langle u', u\rangle u', de modo que u=±uu' = \pm u. Interior e bordo nunca colidem (traço >1> -1 contra =1= -1). Logo EE induz uma bijeção contínua da bola com colagem antipodal — compacta — sobre SO(3)SO(3): um homeomorfismo, e SO(3)RP3SO(3) \cong \mathbb{RP}^3. Um diâmetro de πu\pi u a πu-\pi u tem os extremos colados: ele é um laço em SO(3)SO(3), e sua descrição por EE coincide com a família da questão 10 de rotações em torno de uu varrendo uma volta completa.

21. ρ\rho é um morfismo e ρp1=ρp1=ρpˉ\rho_p^{-1} = \rho_{p^{-1}} = \rho_{\bar p}, de modo que ρpρqρp1=ρpqpˉ\rho_p\rho_q\rho_p^{-1} = \rho_{pq\bar p}; pela questão 16, conjugar uma rotação preserva seu ângulo e gira seu eixo por ρp\rho_p. Involuções: ρq2=ρq2=id\rho_q^2 = \rho_{q^2} = \operatorname{id} se, e somente se, q2=±1q^2 = \pm1. Se q2=1q^2 = 1, então (q1)(q+1)=q21=0(q - 1)(q + 1) = q^2 - 1 = 0 (escalares centrais, de modo que essa fatoração é válida) e q=±1q = \pm1 no anel de divisão H\mathbb H, dando ρq=I\rho_q = I, excluído; q2=1q^2 = -1 com q=a+suq = a + sua2s2+2asu=1a^2 - s^2 + 2as\,u = -1, de modo que a=0a = 0, s=1s = 1: qq é um puro unitário ww, e ρw\rho_w é a meia-volta em torno de ww (ângulo π\pi, questão 5). Centro: se ρq\rho_q comuta com todo ρp\rho_p, então ρpqpˉ=ρq\rho_{pq\bar p} = \rho_q, de modo que pqpˉ=ε(p)qpq\bar p = \varepsilon(p)\,q com ε(p){±1}\varepsilon(p) \in \{\pm1\}; pε(p)=(pqpˉ)q1p \mapsto \varepsilon(p) = (pq\bar p)q^{-1} é contínua no conexo S3S^3 e vale 11 em p=1p = 1, logo 1\equiv 1: qq comuta com todo S3S^3, logo com todo H\mathbb H (reescale), de modo que qRS3={±1}q \in \R \cap S^3 = \{\pm1\} (questão 1) e ρq=I\rho_q = I: o centro é trivial.

22. Como uwu \perp w são puros unitários, uw=uwuw = u\wedge w é puro (questão 18), de modo que

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

é puro, de norma qw=1\abs q\abs w = 1. Então ww=qww=qw2=qw'w = qw\cdot w = qw^2 = -q, e

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

Ambos os eixos ww e w=cosθ2w+sinθ2(uw)w' = \cos\frac\theta2 w + \sin\frac\theta2(u\wedge w) estão no plano uu^\perp ortogonal ao eixo da rotação, e w,w=cosθ2\langle w', w\rangle = \cos\frac\theta2: eles formam o semiângulo θ2\frac\theta2. Essa é a geração clássica: duas meias-voltas em torno de eixos que se encontram sob o ângulo θ2\frac\theta2 compõem a rotação de ângulo θ\theta em torno de sua perpendicular comum.

23. A compacidade é a questão 7. Conexidade por caminhos: SO(3)=ρ(S3)SO(3) = \rho(S^3) é a imagem contínua da esfera, conexa por caminhos; alternativamente, para R=eAR = \eu^{A} com AA antissimétrica (questão 8), tetAt \mapsto \eu^{tA} é um caminho em SO(3)SO(3) de II a RR (ortogonal, pois (etA)T=etA(\eu^{tA})^{\mathsf T} = \eu^{-tA}, de determinante 11 por continuidade a partir de t=0t = 0). Para O(3)O(3): det\det é contínuo sobre {±1}\{\pm1\}, de modo que O(3)O(3) é desconexo, O(3)=SO(3)DSO(3)O(3) = SO(3) \sqcup D\,SO(3) para qualquer DD fixo com detD=1\det D = -1 (por exemplo, D=ID = -I), e a multiplicação à esquerda por DD é um homeomorfismo: exatamente duas componentes, cada uma uma cópia de SO(3)SO(3). O ρ\rho dois para um, sem seção pela questão 19, é assim um honesto recobrimento duplo de um grupo compacto conexo pelo simplesmente conexo S3S^3 — a geometria por trás do truque do cinto.

24. (i) Matrizes:

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}.

O traço 0=1+2cosθ0 = 1 + 2\cos\thetacosθ=12\cos\theta = -\frac12: θ=2π3\theta = \frac{2\pi}3. A parte antissimétrica RRT2\frac{R - R^{\mathsf T}}2 tem entradas que codificam, na escala de sinθ(v3,v2,v1)\sin\theta\,(v_3, -v_2, v_1), o eixo: aqui 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), o que se lê (dicionário AvA_v da Parte V) vsinθ=12(1,1,1)v\sin\theta = \frac12(1, -1, 1); com sin2π3=32\sin\frac{2\pi}3 = \frac{\sqrt3}2: v=13(1,1,1)v = \frac1{\sqrt3}(1, -1, 1). (ii) Quatérnios: 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), e

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). Logo cosθ2=12\cos\frac\theta2 = \frac12: θ=2π3\theta = \frac{2\pi}3, e a parte vetorial 12(ij+k)\frac12(i - j + k) tem direção 13(1,1,1)\frac1{\sqrt3}(1, -1, 1) — a mesma resposta, com a rota dos quatérnios exigindo uma linha de multiplicação em vez de um produto de matrizes: a razão prática pela qual os softwares de voo compõem atitudes em S3S^3.

25. I+KI + K inversível: (I+K)v=0(I + K)v = 00=v,v+v,Kv=v20 = \langle v, v\rangle + \langle v, Kv\rangle = \norm v^2 (a antissimetria mata o segundo termo): v=0v = 0. Ortogonalidade de C=C(K)C = C(K): usando (I±K)T=IK(I \pm K)^{\mathsf T} = I \mp K e o fato de que as quatro matrizes I±KI \pm K, (I±K)1(I \pm K)^{-1} comutam (expressões polinomiais em KK, mais limites):

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 .

Determinante: det(IK)=det((IK)T)=det(I+K)\det(I - K) = \det\bigl((I - K)^{\mathsf T}\bigr) = \det(I + K), de modo que detC=1\det C = 1: CSO(n)C \in SO(n). Sem autovalor 1-1: Cv=vCv = -v significa (IK)w=(I+K)w(I - K)w = -(I + K)w para w=(I+K)1vw = (I + K)^{-1}v, isto é, 2w=02w = 0: v=0v = 0. Inversão: de C(I+K)=IKC(I + K) = I - K, resolva K(I+C)=ICK(I + C) = I - C; como 1SpC-1 \notin \operatorname{Sp}C, I+CI + C é inversível e K=(IC)(I+C)1K = (I - C)(I + C)^{-1}, que é antissimétrica sempre que CC é ortogonal sem autovalor 1-1 (transponha a expressão e use 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). As duas aplicações são mutuamente inversas por construção: uma parametrização racional global da peça aberta densa de SO(n)SO(n) que evita o autovalor 1-1 — sem séries, sem trigonometria e, em dimensão 33, é a substituição do arco metade K=tanθ2AvK = \tan\frac\theta2\,A_v disfarçada.