Mathematics · Livro 3 · Bachelor Year 1

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

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

18Espaços Vetoriais

A álgebra linear começa aqui: os axiomas dos espaços vetoriais isolam o que R2\R^2, R3\R^3, os espaços de polinômios e os espaços de funções têm em comum — pode-se somar e multiplicar por escalares. Dois capítulos constroem a teoria (o Capítulo 19 acrescenta a dimensão); a linguagem que eles instalam — espaço gerado, família livre, base, soma direta — é o pão de cada dia de todos os capítulos seguintes. Ao longo do texto, KK designa R\R ou C\C (escalares).

18.1 Definição e exemplos

Definição 18.1 (Espaço vetorial)

Um espaço vetorial sobre KK é um conjunto EE munido de uma adição que faz de (E,+)(E, +) um grupo abeliano (o zero escrito 0E0_E ou 00), e de uma multiplicação por escalares K×EEK \times E \to E tal que, para todos λ,μK\lambda, \mu \in K e x,yEx, y \in E:

λ(x+y)=λx+λy,(λ+μ)x=λx+μx,λ(μx)=(λμ)x,1x=x.\lambda(x + y) = \lambda x + \lambda y,\quad (\lambda + \mu) x = \lambda x + \mu x,\quad \lambda(\mu x) = (\lambda\mu) x,\quad 1\,x = x .

Consequências: 0x=0E0\,x = 0_E, λ0E=0E\lambda\,0_E = 0_E, (1)x=x(-1)x = -x, e λx=0E    λ=0\lambda x = 0_E \implies \lambda = 0 ou x=0Ex = 0_E (multiplique por λ1\lambda^{-1}).

Demonstração das consequências. Para 0x=0E0\,x = 0_E: de (0+0)x=0x+0x(0 + 0)x = 0x + 0x e (0+0)x=0x(0+0)x = 0x, cancele 0x0x no grupo (E,+)(E, +). Para λ0E\lambda\,0_E: o mesmo truque em λ(0E+0E)\lambda(0_E + 0_E). Para (1)x(-1)x: some xx,

x+(1)x=1x+(1)x=(1+(1))x=0x=0E,x + (-1)x = 1\,x + (-1)x = \bigl(1 + (-1)\bigr)x = 0\,x = 0_E ,

logo (1)x(-1)x é o inverso aditivo de xx. Por fim, se λx=0E\lambda x = 0_E com λ0\lambda \neq 0: multiplique por λ1\lambda^{-1} (os escalares formam um corpo) e use os dois axiomas λ1(λx)=(λ1λ)x=1x=x\lambda^{-1}(\lambda x) = (\lambda^{-1}\lambda)x = 1x = x junto com λ10E=0E\lambda^{-1}0_E = 0_E: x=0Ex = 0_E. Por menores que sejam, estas quatro regras são usadas silenciosamente em cada página que se segue — e a última é exatamente onde os corpos são necessários: sobre os escalares Z\Z, o “espaço” Z/2Z\Z/2\Z a violaria com 2x=02\,x = 0.

Exemplo 18.2

KnK^n (operações coordenada a coordenada); os polinômios K[X]K[X]; as funções F(A,K)\mathcal{F}(A, K) de um conjunto qualquer AA em KK (operações ponto a ponto) — contendo as funções contínuas, as sequências F(N,R)\mathcal{F}(\N, \R), etc.; C\C como espaço vetorial sobre R\R. Em cada caso os axiomas são herdados dos de KK.

Definição 18.3 (Subespaço)

FEF \subseteq E é um subespaço quando 0EF0_E \in F e FF é estável por adição e por multiplicação por escalares — equivalentemente:

Fex,yF, λK,x+λyF.F \neq \emptyset \qquad\text{e}\qquad \forall x, y \in F,\ \forall \lambda \in K,\quad x + \lambda y \in F .

Um subespaço é ele próprio um espaço vetorial. Toda interseção de subespaços é um subespaço; uma união quase nunca é (mesma demonstração do Exercício 7.6).

Exemplo 18.4

Em F(R,R)\mathcal{F}(\R, \R): as funções contínuas, as deriváveis, os polinômios de grau n\leq n (escrito Kn[X]K_n[X] dentro de K[X]K[X]), as soluções de uma equação diferencial linear homogênea (o Teorema 5.10 dizia exatamente isso). Não exemplos: {f:f(0)=1}\{f : f(0) = 1\} (sem o zero), grau exatamente nn (não estável por adição).

Exemplo 18.5 (Subespaço ou não: quatro veredictos, argumentados)

No espaço das sequências reais:

  • {u:u limitada}\{u : u \text{ limitada}\} é um subespaço: 00 é limitada, e se unM\abs{u_n} \leq M, vnM\abs{v_n} \leq M', então un+λvnM+λM\abs{u_n + \lambda v_n} \leq M + \abs\lambda M'.
  • {u:un1}\{u : u_n \to 1\} não é: falta a sequência nula (e a soma de dois membros tende a 22).
  • {u:u monotone}\{u : u \text{ monotone}\} não é: un=nu_n = n e vn=n+(1)nv_n = -n + (-1)^n são monótonas, e a sua soma (1)n(-1)^n não é; a estabilidade por adição é o axioma que falha, embora o conjunto contenha 00 e todos os múltiplos escalares dos seus membros.
  • {u:un+1=un2}\{u : u_{n+1} = u_n^2\} não é: contém 00, mas 2u2u escapa assim que uu é um membro não nulo (2un+1(2un)22u_{n+1} \neq (2u_n)^2 em geral) — a não linearidade está no quadrado.

A ordem de trabalho é sempre a mesma: teste 00 primeiro (o mais barato), depois a estabilidade — e, para refutar, um único par explícito de contraexemplos vale mais do que qualquer dose de dúvida.

18.2 Espaço gerado, somas, somas diretas

Definição 18.6 (Combinações lineares, espaço gerado)

Uma combinação linear da família (x1,,xp)(x_1, \dots, x_p) de vetores de EE é qualquer λ1x1++λpxp\lambda_1 x_1 + \dots + \lambda_p x_p (λiK\lambda_i \in K). O conjunto de todas elas é o espaço gerado Vect(x1,,xp)\operatorname{Vect}(x_1, \dots, x_p): é um subespaço, o menor que contém a família.

Demonstração das duas afirmações. Estabilidade: uma soma de duas combinações lineares λixi+μixi=(λi+μi)xi\sum\lambda_i x_i + \sum\mu_i x_i = \sum(\lambda_i + \mu_i)x_i é ainda uma, e um múltiplo escalar μλixi=(μλi)xi\mu\sum\lambda_i x_i = \sum(\mu\lambda_i) x_i também; a combinação nula mostra que 00 pertence: o espaço gerado é um subespaço. Minimalidade: seja HH um subespaço qualquer contendo x1,,xpx_1, \dots, x_p. Pela estabilidade por multiplicação por escalares, cada λixiH\lambda_i x_i \in H, e pela estabilidade por adição a sua soma pertence a HH: toda combinação linear pertence a HH, isto é, Vect(x1,,xp)H\operatorname{Vect}(x_1, \dots, x_p) \subseteq H. Assim o espaço gerado está contido em todo subespaço que contém a família: é o menor deles.

Definição 18.7 (Soma, soma direta)

Para subespaços F,GF, G de EE:

F+G={u+v:uF, vG}F + G = \{\,u + v : u \in F,\ v \in G\,\}

é um subespaço (o menor que contém FGF \cup G). A soma é direta, escrita FGF \oplus G, quando todo elemento de F+GF + G se decompõe de modo único como u+vu + v; equivalentemente (veja abaixo) quando FG={0}F \cap G = \{0\}. Quando E=FGE = F \oplus G, os subespaços são suplementares em EE.

Exemplo 18.8 (Uma soma de duas retas)

Em R3\R^3, sejam F=Vect((1,0,1))F = \operatorname{Vect}\bigl((1,0,1)\bigr) e G=Vect((0,1,1))G = \operatorname{Vect}\bigl((0,1,1)\bigr). A sua soma é

F+G={a(1,0,1)+b(0,1,1)}={(a, b, a+b)}={(x,y,z):z=x+y},F + G = \{\,a(1,0,1) + b(0,1,1)\,\} = \{(a,\ b,\ a + b)\} = \{(x, y, z) : z = x + y\},

o plano pela origem que contém as duas retas. Ele é estritamente maior do que a união FGF \cup G (a mera cruz das duas retas): o vetor (1,1,2)=(1,0,1)+(0,1,1)(1, 1, 2) = (1,0,1) + (0,1,1) está na soma mas em nenhuma das duas retas. E FG={0}F \cap G = \{0\} (um vetor comum exigiria a(1,0,1)=b(0,1,1)a(1,0,1) = b(0,1,1), cujas duas primeiras coordenadas forçam a=b=0a = b = 0): a soma é direta, e FGF \oplus G é exatamente aquele plano.

Proposição 18.9

F+GF + G é direta se e somente se FG={0}F \cap G = \{0\}.

Demonstração. Se algum w0w \neq 0 pertence a FGF \cap G: w=w+0=0+ww = w + 0 = 0 + w são duas decomposições de ww. Reciprocamente, se u+v=u+vu + v = u' + v' com u,uFu, u' \in F, v,vGv, v' \in G, então uu=vvu - u' = v' - v pertence a FG={0}F \cap G = \{0\}: as decomposições são únicas.

Método 18.10 (Demonstrar que E=FGE = F \oplus G)

Duas coisas a verificar, cada uma com o seu lance de abertura padrão.

  1. Interseção trivial. Tome xFGx \in F \cap G, escreva as duas condições de pertinência e esprema até x=0x = 0. (Nunca argumente por desenho: cf. as armadilhas abaixo.)
  2. A soma é tudo. Tome um xEx \in E arbitrário e produza a decomposição x=f+gx = f + g — seja adivinhando ff a partir do alvo (ff deve satisfazer a propriedade que define FF, o que em geral dita a sua fórmula), seja resolvendo o sistema linear que exprime xx contra geradores de FF e de GG.

Quando a fórmula da decomposição é adivinhada, a unicidade é automática pelo passo 1; quando só a existência é duvidosa, é no passo 2 que mora o trabalho. Os dois exemplos abaixo executam o método: para as funções pares e ímpares, a fórmula de ff é forçada avaliando a identidade desejada em xx e em x-x; para os polinômios que se anulam num ponto, avaliando em aa.

Exemplo 18.11

Em F(R,R)\mathcal{F}(\R, \R), as funções pares P\mathcal{P} e as funções ímpares I\mathcal{I} são suplementares: toda ff se escreve

f(x)=f(x)+f(x)2par+f(x)f(x)2ıˊmpar,f(x) = \underbrace{\frac{f(x) + f(-x)}{2}}_{\text{par}} + \underbrace{\frac{f(x) - f(-x)}{2}}_{\text{ímpar}},

e uma função ao mesmo tempo par e ímpar é nula. (Aplicado a exp\exp, isto é o par (cosh,sinh)(\cosh, \sinh) do Capítulo 4.)

Exemplo 18.12 (Um par suplementar em Kn[X]K_n[X])

Fixe aKa \in K e ponha F={PKn[X]:P(a)=0}F = \{P \in K_n[X] : P(a) = 0\}, G=Vect(1)G = \operatorname{Vect}(1) (as constantes). Então Kn[X]=FGK_n[X] = F \oplus G. De fato FGF \cap G consiste nas constantes que se anulam em aa, isto é, {0}\{0\}; e todo PP se decompõe como

P=(PP(a))F+P(a)G.P = \underbrace{\bigl(P - P(a)\bigr)}_{\in F} + \underbrace{P(a)}_{\in G} .

Vale a pena memorizar a decomposição: subtrair o valor num ponto é a maneira padrão de projetar sobre as “funções que se anulam em aa”. Note que FF é um subespaço grande e GG um pequeno; um par suplementar não precisa ser equilibrado em sentido algum.

Exemplo 18.13 (Um subespaço suplementar nunca é único)

Em R2\R^2, seja F=Vect((1,0))F = \operatorname{Vect}\bigl((1,0)\bigr) (o eixo xx). Tanto G=Vect((0,1))G = \operatorname{Vect}\bigl((0,1)\bigr) quanto G=Vect((1,1))G' = \operatorname{Vect}\bigl((1,1)\bigr) são suplementares de FF: cada um encontra FF apenas em 00, e cada par soma R2\R^2. As decomposições de um mesmo vetor diferem:

(2, 1.5)=(2,0)F+(0,1.5)G=(0.5, 0)F+(1.5, 1.5)G.(2,\ 1.5) = \underbrace{(2, 0)}_{\in F} + \underbrace{(0, 1.5)}_{\in G} = \underbrace{(0.5,\ 0)}_{\in F} + \underbrace{(1.5,\ 1.5)}_{\in G'} .

De fato, toda reta distinta do próprio FF é um suplementar de FF em R2\R^2: os suplementares são abundantes, e falar em “o” suplementar não faz sentido enquanto uma estrutura adicional (um produto interno, Capítulo 23) não singularizar um deles.

Duas decomposições do mesmo ponto de ℝ2 ao longo de F (o eixo x): com o suplementar G (queda vertical) e com o suplementar G' (queda oblíqua). As componentes em F diferem: uma projeção depende da direção da descida.
Duas decomposições do mesmo ponto de R2\R^2 ao longo de FF (o eixo xx): com o suplementar GG (queda vertical) e com o suplementar GG' (queda oblíqua). As componentes em FF diferem: uma projeção depende da direção da descida.

18.3 Famílias livres, famílias geradoras, bases

Definição 18.14

Uma família (x1,,xp)(x_1, \dots, x_p) de vetores de EE é:

  • geradora (de EE) quando Vect(x1,,xp)=E\operatorname{Vect}(x_1,\dots,x_p) = E;
  • livre (os seus vetores linearmente independentes) quando

    λ1x1++λpxp=0    λ1==λp=0;\lambda_1 x_1 + \dots + \lambda_p x_p = 0 \implies \lambda_1 = \dots = \lambda_p = 0 ;

    caso contrário, linearmente dependente;

  • uma base quando é livre e geradora.

Proposição 18.15 (Coordenadas)

(e1,,en)(e_1, \dots, e_n) é uma base de EE se e somente se todo xEx \in E é de modo único uma combinação x=λ1e1++λnenx = \lambda_1 e_1 + \dots + \lambda_n e_n; os escalares λi\lambda_i são as coordenadas de xx na base.

Demonstração. Geradora == existência da decomposição. Unicidade == liberdade: duas decomposições do mesmo xx diferem por uma combinação igual a 00; a liberdade força todos os seus coeficientes — as diferenças de coordenadas — a se anularem. Reciprocamente, uma combinação nula não trivial dá as duas decomposições 0=λiei=0ei0 = \sum \lambda_i e_i = \sum 0\,e_i.

Exemplo 18.16

A base canônica de KnK^n: ei=(0,,1,,0)e_i = (0, \dots, 1, \dots, 0) (11 na posição ii). Os monômios (1,X,X2,,Xn)(1, X, X^2, \dots, X^n): uma base de Kn[X]K_n[X] (liberdade: uma combinação nula é o polinômio nulo, logo todos os coeficientes se anulam, Definição 8.1). Em C\C sobre R\R: a base (1,i)(1, \iu).

Observação 18.17 (As coordenadas são um trabalho de equipe)

A primeira coordenada de xx numa base (e1,,en)(e_1, \dots, e_n) depende de todos os vetores da base, e não apenas de e1e_1. Em R2\R^2: o vetor (3,1)(3, 1) tem primeira coordenada 33 na base canônica, mas primeira coordenada 22 na base ((1,0),(1,1))\bigl((1,0), (1,1)\bigr) — resolva (3,1)=a(1,0)+b(1,1)(3,1) = a(1,0) + b(1,1): b=1b = 1, a=2a = 2. Trocar um vetor da base embaralha todas as coordenadas; o Capítulo 21 empacotará esse embaralhamento na matriz de mudança de base.

Exemplo 18.18 (Testar uma candidata a base, do início ao fim)

Será F=(1+X, 1+X2, X+X2)\mathcal{F} = (1 + X,\ 1 + X^2,\ X + X^2) uma base de R2[X]\R_2[X]? Escreva u1,u2,u3u_1, u_2, u_3 para os três polinômios. Liberdade: uma combinação nula au1+bu2+cu3=0a\,u_1 + b\,u_2 + c\,u_3 = 0 dá, coeficiente a coeficiente,

a+b=0,a+c=0,b+c=0;a + b = 0, \qquad a + c = 0, \qquad b + c = 0 ;

subtraindo as duas primeiras, b=cb = c, e então a terceira dá 2b=02b = 0: a=b=c=0a = b = c = 0, livre. Geradora: em vez de resolver três sistemas, note a combinação simétrica

u1+u2u3=(1+X)+(1+X2)(X+X2)=2,u_1 + u_2 - u_3 = (1 + X) + (1 + X^2) - (X + X^2) = 2 ,

logo 1=12(u1+u2u3)1 = \frac12(u_1 + u_2 - u_3); então

X=u11=12(u1u2+u3),X2=u21=12(u1+u2+u3).X = u_1 - 1 = \tfrac12\bigl(u_1 - u_2 + u_3\bigr), \qquad X^2 = u_2 - 1 = \tfrac12\bigl(-u_1 + u_2 + u_3\bigr).

Os monômios estão no espaço gerado, logo tudo está: F\mathcal{F} é uma base. Como bônus, juntar as três fórmulas exibidas dá as coordenadas de qualquer P=α+βX+γX2P = \alpha + \beta X + \gamma X^2:

P=α+βγ2u1+αβ+γ2u2+α+β+γ2u3.P = \frac{\alpha + \beta - \gamma}{2}\,u_1 + \frac{\alpha - \beta + \gamma}{2}\,u_2 + \frac{-\alpha + \beta + \gamma}{2}\,u_3 .

(Verificação de coerência com P=XP = X: coordenadas (12,12,12)\bigl(\frac12, -\frac12, \frac12\bigr), como encontrado acima.) Duas lições: a simetria numa família costuma esconder uma combinação atalho; e assim que a dimensão estiver disponível (Capítulo 19), toda a metade “geradora” deste trabalho virá sem custo algum — três vetores livres de um espaço de dimensão 33 formam sempre uma base.

Proposição 18.19 (Critérios úteis de liberdade)

  1. Uma família de polinômios não nulos de graus dois a dois distintos é livre.
  2. Acrescentar um vetor a uma família livre mantém-na livre se e somente se o vetor está fora do espaço gerado pela família.
  3. Toda subfamília de uma família livre é livre; toda família que contém uma família geradora é geradora.

Demonstração. (1) Numa combinação nula, olhe para o maior grau presente: o seu coeficiente deve se anular (nada cancela esse grau), e desça em cascata.

(2) Se xVect(x1,,xp)x \in \operatorname{Vect}(x_1, \dots, x_p), a relação xλixi=0x - \sum\lambda_i x_i = 0 é não trivial. Reciprocamente, uma combinação nula não trivial de (x1,,xp,x)(x_1, \dots, x_p, x) deve envolver xx com coeficiente não nulo (senão contradiz a liberdade da família pequena), e isolar xx o coloca no espaço gerado.

(3) Subfamília: uma combinação nula da subfamília é uma da família inteira com os coeficientes ausentes postos a 00; a liberdade da família grande mata todos eles. Superfamília: todo vetor de EE já é uma combinação da parte geradora; dê aos vetores extras o coeficiente 00.

Exemplo 18.20 (O princípio da escada)

Sejam P0,P1,,PnKn[X]P_0, P_1, \dots, P_n \in K_n[X] com degPk=k\deg P_k = k para cada kk (uma “escada” de graus). Então (P0,,Pn)(P_0, \dots, P_n) é uma base de Kn[X]K_n[X]. A liberdade é a Proposição 18.19 (1). Para a propriedade geradora, argumente por descida finita sobre o grau: seja QKn[X]Q \in K_n[X], Q0Q \neq 0, de grau dd, com coeficiente dominante aa, e seja b0b \neq 0 o coeficiente dominante de PdP_d. Então QabPdQ - \frac ab P_d tem grau <d< d (os termos de topo cancelam-se); substituindo QQ por essa diferença e iterando, após no máximo n+1n + 1 passos chega-se ao polinômio nulo, e desfazer as subtrações exprime QQ como combinação dos PkP_k. Duas escadas já encontradas: as potências transladadas ((Xa)k)0kn\bigl((X-a)^k\bigr)_{0 \leq k \leq n} (Exercício 18.4), e os produtos de Newton ((Xx0)(Xx1)(Xxk1))0kn\bigl((X - x_0)(X - x_1)\cdots(X - x_{k-1})\bigr)_{0 \leq k \leq n}, postos a trabalhar no problema de fim de semana.

Exemplo 18.21 (Liberdade em espaços de funções)

Em F(R,R)\mathcal{F}(\R,\R), a família (ea1x,,eapx)(\eu^{a_1 x}, \dots, \eu^{a_p x}) com a1<<apa_1 < \dots < a_p é livre: divida uma combinação nula por eapx\eu^{a_p x} e faça x+x \to +\infty; o último coeficiente morre, e desce-se em cascata (o Exercício 18.8 detalha isto e as suas variantes). A liberdade de funções demonstra-se avaliando: em pontos bem escolhidos, no infinito, ou depois de derivar.

Exemplo 18.22 (Uma relação escondida encolhe um espaço gerado)

Em F(R,R)\mathcal{F}(\R, \R), o que é Vect(1, cos2, sin2)\operatorname{Vect}(1,\ \cos^2,\ \sin^2)? A identidade cos2+sin2=1\cos^2 + \sin^2 = 1 é uma combinação nula não trivial

11+(1)cos2+(1)sin2=0:1\cdot\mathbf{1} + (-1)\cos^2 + (-1)\sin^2 = 0 :

a família é linearmente dependente, e o espaço gerado já é gerado apenas por (1,cos2)(1, \cos^2) (sin2=1cos2\sin^2 = 1 - \cos^2). Essa família menor é livre: a+bcos2x=0a + b\cos^2 x = 0 para todo xx dá, em x=0x = 0 e x=π2x = \frac\pi2: a+b=0a + b = 0 e a=0a = 0. Logo o espaço gerado é um plano dentro do espaço de funções — e ele contém também cos2x=2cos2x1\cos 2x = 2\cos^2 x - 1: famílias de funções trigonométricas com ar de independentes colapsam rotineiramente sob identidades, razão pela qual a liberdade deve ser demonstrada, nunca suposta a partir do tamanho da lista.

Observação 18.23 (Armadilhas comuns)

Quatro armadilhas clássicas. Dois a dois não basta: em R2\R^2, os vetores (1,0)(1,0), (0,1)(0,1), (1,1)(1,1) são dois a dois não proporcionais, e no entanto linearmente dependentes — a liberdade é uma propriedade da família inteira, testada por uma única combinação global, nunca duas a duas. O vetor nulo envenena tudo: toda família que contém 00 é linearmente dependente (10=01\cdot 0 = 0 é uma relação não trivial), por mais inocentes que sejam os outros vetores. União não é soma: FGF \cup G quase nunca é um subespaço (Definição 18.3); o menor subespaço que contém ambos é F+GF + G, em geral muito maior do que a união — em R2\R^2, duas retas distintas têm por união uma cruz e por soma o plano inteiro. Direta exige interseção trivial, não disjunção: dois subespaços nunca são disjuntos (ambos contêm 00); a condição correta é FG={0}F \cap G = \{0\}, e ela deve ser demonstrada, não lida num desenho — cf. o Exemplo 18.13, em que muitos GG diferentes servem. A liberdade depende dos escalares: o par (1,i)(1, \iu) é livre em C\C visto como espaço vetorial sobre R\R, mas linearmente dependente em C\C visto como espaço vetorial sobre C\C (i1+(1)i=0\iu\cdot 1 + (-1)\cdot\iu = 0). Saiba sempre qual corpo está agindo antes de declarar livre uma família — o problema de fim de semana do Capítulo 19 transforma exatamente essa sensibilidade em demonstrações de irracionalidade.

Observação 18.24 (Para onde vai esta linguagem)

Tudo o que vem depois deste capítulo fala a linguagem instalada aqui. O Capítulo 19 conta vetores de uma base e transforma “livre” e “geradora” em desigualdades sobre um único inteiro, a dimensão. O Capítulo 20 estuda as aplicações compatíveis com as duas operações; as somas diretas tornam-se lá projetores. O Capítulo 21 codifica vetores pelas suas coordenadas numa base — a Proposição 18.15 é a licença para essa codificação — e o Capítulo 23 acrescenta comprimentos e ângulos sobre a estrutura linear. No volume do segundo ano de graduação os mesmos axiomas, palavra por palavra, valem sobre corpos arbitrários e em dimensão infinita; nada neste capítulo usou a finitude em lugar algum.

Observação 18.25 (Três fios a seguir ao longo do Livro 3)

Observe três ideias específicas deste capítulo crescerem. O princípio da escada (Exemplo 18.20) reaparece como a base de Newton no problema de fim de semana deste capítulo, como a base binomial (Bk)(B_k) ali, e como o truque do alternante polinomial no problema de fim de semana do Capítulo 22: um lema, três dividendos sem determinantes. A avaliação como teste de liberdade (Exemplo 18.21) torna-se o isomorfismo de interpolação do Capítulo 20, depois o critério de Vandermonde do Capítulo 22, depois o teste de Gram do Capítulo 23: o mesmo reflexo, afiado três vezes. As somas diretas (Definição 18.7) tornam-se projetores no Capítulo 20, decomposições ortogonais E=FFE = F \oplus F^\perp no Capítulo 23, e a decomposição explicado-mais-resíduo dos mínimos quadrados no problema de fim de semana do Capítulo 25. Muito pouco deste livro não é, no fundo, uma destas três ideias vestindo roupa nova.

18.4 Exercícios

Exercício 18.1

Quais dos seguintes são subespaços?

  1. {(x,y,z)R3:x+2yz=0}\{(x, y, z) \in \R^3 : x + 2y - z = 0\};
  2. {(x,y,z)R3:x+2yz=1}\{(x, y, z) \in \R^3 : x + 2y - z = 1\};
  3. {(x,y)R2:xy0}\{(x, y) \in \R^2 : xy \geq 0\};
  4. {PR[X]:P(1)=0}\{P \in \R[X] : P(1) = 0\};
  5. {fF(R,R):f limitada}\{f \in \mathcal{F}(\R,\R) : f \text{ limitada}\}.
Solução

Solução de Exercício 18.1.

  1. Sim: contém 00, e a equação que o define é linear (estável por x+λyx + \lambda y).
  2. Não: não contém (0,0,0)(0,0,0).
  3. Não: (1,0)(1, 0) e (0,1)(0, -1) pertencem (xy=0xy = 0), mas a sua soma (1,1)(1, -1) não (xy=1<0xy = -1 < 0).
  4. Sim: 00 se anula em 11; (P+λQ)(1)=P(1)+λQ(1)=0(P + \lambda Q)(1) = P(1) + \lambda Q(1) = 0.
  5. Sim: a função nula é limitada; se fM\abs f \leq M e gM\abs g \leq M', então f+λgM+λM\abs{f + \lambda g} \leq M + \abs\lambda M'.

Exercício 18.2

Em R3\R^3, (1,2,1)(1, 2, 1) pertence a Vect((1,0,1),(1,1,0))\operatorname{Vect}\bigl((1,0,1),\, (1,1,0)\bigr)? E (2,1,1)(2, 1, 1)? Descreva Vect((1,0,1),(1,1,0))\operatorname{Vect}\bigl((1,0,1),(1,1,0)\bigr) por uma equação.

Solução

Solução de Exercício 18.2.

(1,2,1)=a(1,0,1)+b(1,1,0)(1,2,1) = a(1,0,1) + b(1,1,0) exige a+b=1a + b = 1, b=2b = 2, a=1a = 1: incompatível (a+b=31a + b = 3 \neq 1): não está no espaço gerado. (2,1,1)=a(1,0,1)+b(1,1,0)(2,1,1) = a(1,0,1) + b(1,1,0): b=1b = 1, a=1a = 1, a+b=2a + b = 2: compatível, logo (2,1,1)=(1,0,1)+(1,1,0)(2,1,1) = (1,0,1) + (1,1,0), está no espaço gerado.

Equação: (x,y,z)=(a+b,b,a)(x, y, z) = (a + b, b, a) significa x=y+zx = y + z: o espaço gerado é o plano {xyz=0}\{x - y - z = 0\}.

Exercício 18.3

Decida a liberdade em R3\R^3:   ((1,1,0),(1,0,1),(0,1,1))\;\bigl((1,1,0), (1,0,1), (0,1,1)\bigr);   ((1,2,3),(2,4,6))\;\bigl((1,2,3), (2,4,6)\bigr);   ((1,0,0),(1,1,0),(1,1,1),(0,1,1))\;\bigl((1,0,0), (1,1,0), (1,1,1), (0,1,1)\bigr).

Solução

Solução de Exercício 18.3.

Primeira família: λ(1,1,0)+μ(1,0,1)+ν(0,1,1)=0\lambda(1,1,0) + \mu(1,0,1) + \nu(0,1,1) = 0λ+μ=0\lambda + \mu = 0, λ+ν=0\lambda + \nu = 0, μ+ν=0\mu + \nu = 0: somando, 2(λ+μ+ν)=02(\lambda + \mu + \nu) = 0, e subtraindo cada equação, λ=μ=ν=0\lambda = \mu = \nu = 0: livre.

Segunda: (2,4,6)=2(1,2,3)(2,4,6) = 2(1,2,3): linearmente dependente.

Terceira: quatro vetores em R3\R^3 — necessariamente dependentes assim que a dimensão estiver disponível (Capítulo 19); diretamente: (0,1,1)=(1,0,0)+0(1,1,0)+(1,1,1)(0,1,1) = -(1,0,0) + 0\cdot(1,1,0) + (1,1,1), de fato (1,0,0)+(1,1,1)=(0,1,1)(-1,0,0) + (1,1,1) = (0,1,1): uma relação não trivial.

Exercício 18.4

Demonstre que (1,X1,(X1)2,(X1)3)(1, X - 1, (X-1)^2, (X-1)^3) é uma base de R3[X]\R_3[X], e dê as coordenadas de X3X^3 nela. (Taylor em 11!)

Solução

Solução de Exercício 18.4.

Os polinômios 1,(X1),(X1)2,(X1)31, (X-1), (X-1)^2, (X-1)^3 têm graus distintos 0,1,2,30, 1, 2, 3: livres (Proposição 18.19 (1)); quatro vetores livres e geradores (todo PR3[X]P \in \R_3[X] se desenvolve em potências de X1X - 1, por exemplo por Taylor para polinômios, cf. a demonstração da Proposição 8.11): uma base. Para X3X^3, Taylor em 11: P=X3P = X^3, P(1)=1P(1) = 1, P(1)=3P'(1) = 3, P(1)=6P''(1) = 6, P(1)=6P'''(1) = 6:

X3=1+3(X1)+3(X1)2+(X1)3,X^3 = 1 + 3(X - 1) + 3(X-1)^2 + (X-1)^3 ,

coordenadas (1,3,3,1)(1, 3, 3, 1) (a linha de Pascal, como se espera de X3=((X1)+1)3X^3 = ((X-1)+1)^3).

Exercício 18.5 ★★

Em R4\R^4, sejam F={(x,y,z,t):x=y=z}F = \{(x,y,z,t) : x = y = z\} e G={(x,y,z,t):x=t=0}G = \{(x,y,z,t) : x = t = 0\}. Demonstre que FG=R4F \oplus G = \R^4, e decomponha (1,2,3,4)(1,2,3,4) em consequência.

Solução

Solução de Exercício 18.5.

FGF \cap G: as condições x=y=zx = y = z e x=t=0x = t = 0 dão juntas x=0x = 0, donde y=z=0y = z = 0, e t=0t = 0: a interseção é {0}\{0\}. Soma: dado (x,y,z,t)(x,y,z,t), procure (a,a,a,b)F(a,a,a,b) \in F e (0,c,d,0)G(0,c,d,0) \in G cuja soma seja ele: a=xa = x, b=tb = t, c=yxc = y - x, d=zxd = z - x: sempre possível. Logo R4=FG\R^4 = F \oplus G, e

(1,2,3,4)=(1,1,1,4)+(0,1,2,0).(1,2,3,4) = (1,1,1,4) + (0,1,2,0) .

Exercício 18.6 ★★

No espaço das sequências, seja FF o conjunto das sequências convergentes e G=Vect(u)G = \operatorname{Vect}(u), onde un=(1)nu_n = (-1)^n. Demonstre que FG={0}F \cap G = \{0\}. Será F+GF + G o espaço inteiro das sequências?

Solução

Solução de Exercício 18.6.

Um elemento de GG é λu\lambda u; se ele converge, então (como λun=λ(1)n\lambda u_n = \lambda(-1)^n tem os dois limites de subsequências ±λ\pm \lambda) necessariamente λ=0\lambda = 0: FG={0}F \cap G = \{0\}.

F+GF + G não é tudo: consiste nas sequências da forma cn+λ(1)nc_n + \lambda(-1)^n com (cn)(c_n) convergente. A sequência vn=nv_n = n não é dessa forma (vnλ(1)nv_n - \lambda(-1)^n é ilimitada, nunca convergente). Logo FGF \oplus G \subsetneq (espaço de todas as sequências).

Exercício 18.7 ★★

Sejam F,G,HF, G, H subespaços de EE. Demonstre que

F(G+(FH))=(FG)+(FH),F \cap (G + (F \cap H)) = (F \cap G) + (F \cap H),

e mostre por um exemplo em R2\R^2 que a distributividade irrestrita F(G+H)=(FG)+(FH)F \cap (G + H) = (F\cap G) + (F \cap H) é falsa.

Solução

Solução de Exercício 18.7.

(\supseteq) Tanto FGF \cap G quanto FHF \cap H estão em FF, e a sua soma está em G+(FH)G + (F \cap H): a inclusão segue porque o lado esquerdo é um subespaço que contém as duas partes — concretamente, um elemento g+hg + h com gFGg \in F\cap G, hFHh \in F \cap H está em FF (soma de dois elementos de FF) e em G+(FH)G + (F \cap H).

(\subseteq) Seja xFx \in F com x=g+hx = g + h, gGg \in G, hFHh \in F \cap H. Então g=xhFg = x - h \in F (diferença de elementos de FF), logo gFGg \in F \cap G, e x=g+h(FG)+(FH)x = g + h \in (F \cap G) + (F \cap H).

Contraexemplo à distributividade plena em R2\R^2: F=Vect(1,1)F = \operatorname{Vect}(1,1), G=Vect(1,0)G = \operatorname{Vect}(1,0), H=Vect(0,1)H = \operatorname{Vect}(0,1). Então G+H=R2G + H = \R^2, logo F(G+H)=FF \cap (G+H) = F, ao passo que FG=FH={0}F \cap G = F \cap H = \{0\}: o lado direito é {0}F\{0\} \neq F.

Exercício 18.8 ★★★

Demonstre que as seguintes famílias de F(R,R)\mathcal{F}(\R, \R) são livres:

  1. (ea1x,,eapx)(\eu^{a_1 x}, \dots, \eu^{a_p x}) para a1<<apa_1 < \dots < a_p;
  2. (cosx,sinx,cos2x,sin2x)(\cos x, \sin x, \cos 2x, \sin 2x);
  3. (xxa1,,xxap)(x \mapsto \abs{x - a_1}, \dots, x \mapsto \abs{x - a_p}) para aia_i distintos (a derivabilidade falha em exatamente um ponto por função).
Solução

Solução de Exercício 18.8.

  1. Suponha iλieaix=0\sum_{i} \lambda_i \eu^{a_i x} = 0 para todo xx. Multiplique por eapx\eu^{-a_p x}: λp+i<pλie(aiap)xλp\lambda_p + \sum_{i < p} \lambda_i \eu^{(a_i - a_p)x} \to \lambda_p quando x+x \to +\infty (cada expoente aiap<0a_i - a_p < 0). O lado esquerdo é identicamente 00, logo λp=0\lambda_p = 0; itere para baixo.
  2. Seja acosx+bsinx+ccos2x+dsin2x=0a\cos x + b \sin x + c \cos 2x + d\sin 2x = 0 para todo xx. Avalie em x=0x = 0: a+c=0a + c = 0; em x=πx = \pi: a+c=0-a + c = 0; logo a=c=0a = c = 0, e a relação reduz-se a bsinx+dsin2x=0b\sin x + d \sin 2x = 0. Avalie em x=π2x = \frac\pi2: b=0b = 0; depois em x=π4x = \frac\pi4: d=0d = 0.
  3. Suponha λixai=0\sum \lambda_i \abs{x - a_i} = 0 para todo xx. A função ijλixai\sum_{i \neq j} \lambda_i\abs{x - a_i} é derivável em aja_j (cada termo o é, longe do seu próprio bico), logo λjxaj-\lambda_j \abs{x - a_j}, a diferença delas, também deve ser derivável em aja_j — o que força λj=0\lambda_j = 0 (\abs{\,\cdot\,} tem um bico). Isto vale para todo jj.

Exercício 18.9 ★★★

Seja EE um espaço vetorial sobre KK e F,G,HF, G, H subespaços com F+G=F+HF + G = F + H, FG=FHF \cap G = F \cap H e GHG \subseteq H. Demonstre que G=HG = H. Dê um contraexemplo sem a hipótese GHG \subseteq H.

Solução

Solução de Exercício 18.9.

Seja hHh \in H. Como hHF+H=F+Gh \in H \subseteq F + H = F + G, escreva h=f+gh = f + g com fFf \in F, gGg \in G. Então f=hgHf = h - g \in H (os dois termos em HH, usando GHG \subseteq H), logo fFH=FGGf \in F \cap H = F \cap G \subseteq G, e h=f+gGh = f + g \in G. Portanto HGH \subseteq G, e com a hipótese GHG \subseteq H: igualdade.

Contraexemplo sem GHG \subseteq H: em R2\R^2, tome F=Vect(1,0)F = \operatorname{Vect}(1,0), G=Vect(0,1)G = \operatorname{Vect}(0,1), H=Vect(1,1)H = \operatorname{Vect}(1,1): então F+G=F+H=R2F + G = F + H = \R^2 e FG=FH={0}F \cap G = F \cap H = \{0\}, mas GHG \neq H.

Exercício 18.10 ★★

Em R[X]\R[X], seja P\mathcal P o conjunto dos polinômios pares (P(X)=P(X)P(-X) = P(X)) e I\mathcal I o conjunto dos ímpares (P(X)=P(X)P(-X) = -P(X)). Demonstre que R[X]=PI\R[X] = \mathcal P \oplus \mathcal I, e mostre que P=Vect(1,X2,X4,)\mathcal P = \operatorname{Vect}(1, X^2, X^4, \dots), isto é, que os polinômios pares são exatamente os polinômios em X2X^2.

Solução

Solução de Exercício 18.10.

Os dois conjuntos são subespaços (as condições que os definem são lineares e valem para 00). Decomposição: para PR[X]P \in \R[X],

P(X)=P(X)+P(X)2P+P(X)P(X)2I,P(X) = \underbrace{\frac{P(X) + P(-X)}{2}}_{\in\,\mathcal P} + \underbrace{\frac{P(X) - P(-X)}{2}}_{\in\,\mathcal I},

e um polinômio ao mesmo tempo par e ímpar satisfaz P=PP = -P, logo P=0P = 0: a soma é direta e vale R[X]\R[X].

Agora seja P=kakXkP = \sum_k a_k X^k par. Então P(X)P(X)=2k ıˊmparakXkP(X) - P(-X) = 2 \sum_{k \text{ ímpar}} a_k X^k é o polinômio nulo, logo todo coeficiente de grau ímpar se anula (Definição 8.1): PVect(1,X2,X4,)P \in \operatorname{Vect}(1, X^2, X^4, \dots), isto é, P=Q(X2)P = Q(X^2) para um polinômio QQ. Reciprocamente, todo polinômio em X2X^2 é par.

Exercício 18.11 ★★

Seja (x1,x2,x3)(x_1, x_2, x_3) uma família livre de um espaço vetorial real EE. Demonstre que (x1+x2, x2+x3, x3+x1)(x_1 + x_2,\ x_2 + x_3,\ x_3 + x_1) é livre. A família análoga de quatro vetores (x1+x2, x2+x3, x3+x4, x4+x1)(x_1 + x_2,\ x_2 + x_3,\ x_3 + x_4,\ x_4 + x_1) é livre quando (x1,x2,x3,x4)(x_1, x_2, x_3, x_4) o é?

Solução

Solução de Exercício 18.11.

Suponha a(x1+x2)+b(x2+x3)+c(x3+x1)=0a(x_1 + x_2) + b(x_2 + x_3) + c(x_3 + x_1) = 0. Reagrupando na família livre (x1,x2,x3)(x_1, x_2, x_3):

(a+c)x1+(a+b)x2+(b+c)x3=0    a+c=a+b=b+c=0.(a + c)\,x_1 + (a + b)\,x_2 + (b + c)\,x_3 = 0 \implies a + c = a + b = b + c = 0 .

Subtrair as duas primeiras equações dá c=bc = b; a terceira então dá 2b=02b = 0, logo b=c=0b = c = 0, e depois a=0a = 0: a família é livre.

Para quatro vetores a família análoga é sempre linearmente dependente:

(x1+x2)(x2+x3)+(x3+x4)(x4+x1)=0(x_1 + x_2) - (x_2 + x_3) + (x_3 + x_4) - (x_4 + x_1) = 0

é uma combinação nula não trivial (coeficientes 1,1,1,11, -1, 1, -1), qualquer que seja (x1,x2,x3,x4)(x_1, x_2, x_3, x_4). É a paridade do comprimento do ciclo que decide.

Exercício 18.12 ★★★

Seja EE um espaço vetorial sobre R\R (ou C\C) e F1,,FkF_1, \dots, F_k subespaços próprios de EE (cada FiEF_i \neq E).

  1. Trate diretamente o caso k=2k = 2: se F1⊈F2F_1 \not\subseteq F_2 e F2⊈F1F_2 \not\subseteq F_1, tome xF1F2x \in F_1 \setminus F_2 e yF2F1y \in F_2 \setminus F_1 e localize x+yx + y.
  2. Demonstre em geral que EF1FkE \neq F_1 \cup \dots \cup F_k: um espaço vetorial sobre um corpo infinito nunca é uma união finita de subespaços próprios. (Tome kk mínimo, escolha xF1x \in F_1 fora dos outros FiF_i, escolha yF1y \notin F_1, e siga a reta ty+txt \mapsto y + tx.)
Solução

Solução de Exercício 18.12.

  1. Se F1F2F_1 \subseteq F_2 ou F2F1F_2 \subseteq F_1, a união é um dos dois, logo própria. Caso contrário, escolha xF1F2x \in F_1 \setminus F_2 e yF2F1y \in F_2 \setminus F_1, e considere x+yx + y. Se x+yF1x + y \in F_1, então y=(x+y)xF1y = (x + y) - x \in F_1: contradição. Se x+yF2x + y \in F_2, então xF2x \in F_2: contradição. Logo x+yF1F2x + y \notin F_1 \cup F_2, e EF1F2E \neq F_1 \cup F_2.
  2. Suponha, por absurdo, que E=F1FkE = F_1 \cup \dots \cup F_k, com kk escolhido mínimo entre todas essas coberturas. A minimalidade proíbe F1F2FkF_1 \subseteq F_2 \cup \dots \cup F_k (senão descarte F1F_1), logo existe xF1x \in F_1 com xFix \notin F_i para todo i2i \geq 2. Como F1F_1 é próprio, escolha yF1y \notin F_1. Para cada escalar tt, o vetor y+txy + t x está em algum FiF_i. Ele nunca está em F1F_1: caso contrário y=(y+tx)txF1y = (y + tx) - tx \in F_1 (pois xF1x \in F_1). O corpo é infinito, logo escolha kk escalares distintos t1,,tkt_1, \dots, t_k: os kk vetores y+tjxy + t_j x caem nos k1k - 1 subespaços F2,,FkF_2, \dots, F_k, e dois deles, digamos y+txy + t x e y+txy + t' x com ttt \neq t', estão no mesmo FiF_i (i2i \geq 2). Então a sua diferença (tt)xFi(t - t')x \in F_i, logo xFix \in F_i: contradição. Portanto não existe cobertura finita por subespaços próprios.

18.5 Problema: interpolação, três bases para um espaço

Problema 18.1

Fixe n+1n + 1 pontos distintos x0,x1,,xnx_0, x_1, \dots, x_n de R\R. Este problema revisita a interpolação de Lagrange (o Teorema 8.23) com os olhos deste capítulo: o espaço Rn[X]\R_n[X] carrega três bases naturais — a de Lagrange, a de Newton e (para pontos igualmente espaçados) a base binomial — e cada base torna fácil uma pergunta. O caminho termina num teorema aritmético genuíno: a caracterização de Pólya dos polinômios que levam Z\Z em Z\Z.

Parte I — A base de Lagrange. Para 0in0 \leq i \leq n ponha

Li  =  jiXxjxixj    Rn[X].L_i \;=\; \prod_{j \neq i} \frac{X - x_j}{x_i - x_j} \;\in\; \R_n[X].
  1. Verifique que degLi=n\deg L_i = n e que Li(xj)=1L_i(x_j) = 1 se j=ij = i, e 00 se jij \neq i.
  2. Demonstre que a família (L0,,Ln)(L_0, \dots, L_n) é livre.
  3. Demonstre que para todo PRn[X]P \in \R_n[X],

    P  =  i=0nP(xi)Li,P \;=\; \sum_{i=0}^{n} P(x_i)\, L_i ,

    e conclua que (L0,,Ln)(L_0, \dots, L_n) é uma base de Rn[X]\R_n[X]. (Considere a diferença dos dois lados e conte as suas raízes, Corolário 8.8.)

  4. Deduza o teorema da interpolação: para quaisquer valores y0,,ynRy_0, \dots, y_n \in \R existe um único PRn[X]P \in \R_n[X] com P(xi)=yiP(x_i) = y_i para todo ii. Na base de Lagrange, quais são as coordenadas de um polinômio PP?
  5. Demonstre as identidades

    i=0nLi=1e, para 0kn,i=0nxikLi=Xk.\sum_{i=0}^{n} L_i = 1 \qquad\text{e, para } 0 \leq k \leq n,\qquad \sum_{i=0}^{n} x_i^{k}\, L_i = X^{k} .

Parte II — A base de Newton e as diferenças divididas. Ponha N0=1N_0 = 1 e Nk=(Xx0)(Xx1)(Xxk1)N_k = (X - x_0)(X - x_1) \cdots (X - x_{k-1}) para 1kn1 \leq k \leq n. Para uma função ff definida nos nós, defina as diferenças divididas por f[xi]=f(xi)f[x_i] = f(x_i) e

f[xi,,xi+k]  =  f[xi+1,,xi+k]f[xi,,xi+k1]xi+kxi.f[x_i, \dots, x_{i+k}] \;=\; \frac{f[x_{i+1}, \dots, x_{i+k}] - f[x_i, \dots, x_{i+k-1}]} {x_{i+k} - x_i} .
  1. Demonstre que (N0,N1,,Nn)(N_0, N_1, \dots, N_n) é uma base de Rn[X]\R_n[X].
  2. Calcule f[x0,x1]f[x_0, x_1] e f[x0,x1,x2]f[x_0, x_1, x_2] em termos dos valores de ff, e depois calcule todas as diferenças divididas de f(x)=x2f(x) = x^2 em três nós arbitrários.
  3. (Lema de Aitken) Seja RR interpolando ff em x0,,xn1x_0, \dots, x_{n-1} e QQ interpolando ff em x1,,xnx_1, \dots, x_n, ambos de grau n1\leq n - 1. Demonstre que

    S  =  (Xx0)Q(Xxn)Rxnx0S \;=\; \frac{(X - x_0)\,Q - (X - x_n)\,R}{x_n - x_0}

    interpola ff em x0,x1,,xnx_0, x_1, \dots, x_n.

  4. Deduza, por indução sobre o número de nós, que o coeficiente de XkX^{k} no interpolante de ff em x0,,xkx_0, \dots, x_k é exatamente f[x0,,xk]f[x_0, \dots, x_k].
  5. Demonstre a fórmula de interpolação de Newton: o interpolante de ff em x0,,xnx_0, \dots, x_n é

    P  =  k=0nf[x0,,xk]Nk,P \;=\; \sum_{k=0}^{n} f[x_0, \dots, x_k]\, N_k ,

    e deduza a fórmula fechada

    f[x0,,xk]  =  i=0kf(xi)ji,jk(xixj),f[x_0, \dots, x_k] \;=\; \sum_{i=0}^{k} \frac{f(x_i)}{\prod_{j \neq i,\, j \leq k} (x_i - x_j)} ,

    que mostra que f[x0,,xk]f[x_0, \dots, x_k] não depende da ordenação dos nós.

Parte III — Nós igualmente espaçados: o operador de diferenças. A partir daqui os nós são 0,1,2,0, 1, 2, \dots e, para um polinômio PP, pomos

ΔP(X)=P(X+1)P(X),Bk=X(X1)(Xk+1)k!(B0=1).\Delta P(X) = P(X + 1) - P(X), \qquad B_k = \frac{X(X-1)\cdots(X-k+1)}{k!} \quad (B_0 = 1).
  1. Mostre que se degP=m1\deg P = m \geq 1 com coeficiente dominante aa, então degΔP=m1\deg \Delta P = m - 1 com coeficiente dominante mam\,a, e que Δ\Delta mata as constantes.
  2. Mostre que (B0,B1,,Bn)(B_0, B_1, \dots, B_n) é uma base de Rn[X]\R_n[X] e que ΔBk=Bk1\Delta B_k = B_{k-1} para k1k \geq 1.
  3. (Fórmula das diferenças progressivas de Newton) Demonstre que todo PRn[X]P \in \R_n[X] satisfaz

    P  =  k=0n(ΔkP)(0)Bk.P \;=\; \sum_{k=0}^{n} \bigl(\Delta^{k} P\bigr)(0)\, B_k .
  4. Demonstre que para todo k0k \geq 0,

    (ΔkP)(0)  =  j=0k(1)kj(kj)P(j).\bigl(\Delta^{k} P\bigr)(0) \;=\; \sum_{j=0}^{k} (-1)^{k-j} \binom{k}{j} P(j) .
  5. Mostre que se degP=n\deg P = n com coeficiente dominante ana_n, então ΔnP\Delta^{n} P é a constante n!ann!\,a_n e Δn+1P=0\Delta^{n+1} P = 0.

Parte IV — Polinômios a valores inteiros. Um polinômio PR[X]P \in \R[X] é a valores inteiros quando P(m)ZP(m) \in \Z para todo mZm \in \Z.

  1. Demonstre que cada BkB_k é a valores inteiros. (Trate mkm \geq k, 0m<k0 \leq m < k e m<0m < 0 separadamente; para m=q<0m = -q < 0, mostre que Bk(q)=(1)k(q+k1k)B_k(-q) = (-1)^k \binom{q + k - 1}{k}.)
  2. Demonstre a caracterização de Pólya: PRn[X]P \in \R_n[X] é a valores inteiros se e somente se as suas coordenadas na base (B0,,Bn)(B_0, \dots, B_n) são inteiras.
  3. Deduza: se PRn[X]P \in \R_n[X] assume valores inteiros em n+1n + 1 inteiros consecutivos a,a+1,,a+na, a+1, \dots, a+n, então PP é a valores inteiros. (Translade: aplique o estudo a Q(X)=P(X+a)Q(X) = P(X + a).)
  4. Deduza da questão 16 que um produto de kk inteiros consecutivos é sempre divisível por k!k!.
  5. Seja P=X(X+1)(2X+1)6P = \dfrac{X(X+1)(2X+1)}{6}. Calcule a sua tabela de Newton em 0,1,2,30, 1, 2, 3, escreva PP na base (Bk)(B_k), e conclua que PP é a valores inteiros embora nenhum dos seus coeficientes monomiais seja inteiro. Verifique ΔP=(X+1)2\Delta P = (X+1)^2 e deduza P(m)=12+22++m2P(m) = 1^2 + 2^2 + \dots + m^2 para mNm \in \N.

Parte V — Dividendos.

  1. Seja PRn[X]P \in \R_n[X] interpolando os valores 2i2^i em i=0,1,,ni = 0, 1, \dots, n. Mostre que P=B0+B1++BnP = B_0 + B_1 + \dots + B_n e que P(n+1)=2n+11P(n + 1) = 2^{n+1} - 1: o “padrão de duplicação” quebra sempre no ponto imediatamente seguinte.
  2. (Primitiva discreta) Demonstre que para todos os inteiros m1m \geq 1 e k0k \geq 0,

    j=0m1Bk(j)  =  Bk+1(m),\sum_{j=0}^{m-1} B_k(j) \;=\; B_{k+1}(m),

    isto é, a identidade do taco de hóquei j=km1(jk)=(mk+1)\sum_{j=k}^{m-1} \binom{j}{k} = \binom{m}{k+1}.

  3. Expanda X2X^2 e X3X^3 na base (Bk)(B_k) e deduza fórmulas fechadas para j=0m1j2\sum_{j=0}^{m-1} j^2 e j=0m1j3\sum_{j=0}^{m-1} j^3; recupere a identidade de Nicômaco 13++m3=(1++m)21^3 + \dots + m^3 = (1 + \dots + m)^2.
  4. Tome n=2n = 2 e os nós 0,1,20, 1, 2. Escreva as coordenadas de X2X^2 nas três bases deste problema: a base monomial, a base de Lagrange, a base de Newton. Confira as três respostas contra as questões 4 e 9.
  5. Síntese. Em quatro frases: qual conceito de espaço vetorial torna a questão 4 automática; por que a base de Newton calcula as coordenadas recursivamente enquanto a base de Lagrange as lê instantaneamente; qual critério de liberdade as duas bases partilham; e em que sentido preciso o teorema de Pólya diz que a integralidade de um polinômio é uma propriedade das suas coordenadas na base certa.
Solução

Solução de Problema 18.1.

1. LiL_i é um produto de nn fatores lineares dividido por uma constante não nula (os xix_i são distintos), logo degLi=n\deg L_i = n. Avaliando em xjx_j com jij \neq i: o fator XxjX - x_j do numerador se anula, logo Li(xj)=0L_i(x_j) = 0. Em xix_i, numerador e denominador coincidem: Li(xi)=1L_i(x_i) = 1.

2. Suponha iλiLi=0\sum_i \lambda_i L_i = 0. Avalie em xjx_j: todos os termos morrem exceto λjLj(xj)=λj\lambda_j L_j(x_j) = \lambda_j, logo λj=0\lambda_j = 0 para todo jj: a família é livre.

3. Seja D=PiP(xi)LiD = P - \sum_i P(x_i) L_i. Então degDn\deg D \leq n e, pela questão 1, D(xj)=P(xj)P(xj)=0D(x_j) = P(x_j) - P(x_j) = 0 nos n+1n + 1 pontos distintos x0,,xnx_0, \dots, x_n. Um polinômio não nulo de grau n\leq n tem no máximo nn raízes (Corolário 8.8), logo D=0D = 0. Assim todo PRn[X]P \in \R_n[X] é uma combinação dos LiL_i: a família é geradora e, com a questão 2, uma base.

4. Dados y0,,yny_0, \dots, y_n, o polinômio P=iyiLiP = \sum_i y_i L_i tem grau n\leq n e interpola. Unicidade: um PP interpolador tem, pela questão 3, coordenadas (P(x0),,P(xn))=(y0,,yn)(P(x_0), \dots, P(x_n)) = (y_0, \dots, y_n) na base (Li)(L_i), e as coordenadas numa base são únicas (Proposição 18.15). As coordenadas de PP na base de Lagrange são os seus valores nos nós — é este todo o interesse da base.

5. Aplique a questão 3 a P=XkP = X^k (0kn0 \leq k \leq n):

Xk=i=0nxikLi,X^k = \sum_{i=0}^{n} x_i^{k} L_i ,

e k=0k = 0iLi=1\sum_i L_i = 1.

6. degNk=k\deg N_k = k exatamente: a família (N0,,Nn)(N_0, \dots, N_n) é uma escada de graus em Rn[X]\R_n[X], logo uma base pelo Exemplo 18.20 (liberdade pela Proposição 18.19 (1), geração por descida finita sobre o grau).

7. Da recorrência,

f[x0,x1]=f(x1)f(x0)x1x0,f[x0,x1,x2]=f[x1,x2]f[x0,x1]x2x0.f[x_0, x_1] = \frac{f(x_1) - f(x_0)}{x_1 - x_0}, \qquad f[x_0, x_1, x_2] = \frac{f[x_1, x_2] - f[x_0, x_1]}{x_2 - x_0}.

Para f(x)=x2f(x) = x^2:

f[x0,x1]=x12x02x1x0=x0+x1,f[x_0, x_1] = \frac{x_1^2 - x_0^2}{x_1 - x_0} = x_0 + x_1,

e então

f[x0,x1,x2]=(x1+x2)(x0+x1)x2x0=x2x0x2x0=1.f[x_0, x_1, x_2] = \frac{(x_1 + x_2) - (x_0 + x_1)}{x_2 - x_0} = \frac{x_2 - x_0}{x_2 - x_0} = 1 .

8. degSn\deg S \leq n, pois Q,RQ, R têm grau n1\leq n - 1. Em x0x_0: S(x0)=(x0xn)R(x0)xnx0=R(x0)=f(x0)S(x_0) = \frac{-(x_0 - x_n) R(x_0)}{x_n - x_0} = R(x_0) = f(x_0). Em xnx_n: S(xn)=(xnx0)Q(xn)xnx0=Q(xn)=f(xn)S(x_n) = \frac{(x_n - x_0) Q(x_n)}{x_n - x_0} = Q(x_n) = f(x_n). Num nó interior xix_i (1in11 \leq i \leq n-1), tanto QQ quanto RR tomam o valor f(xi)f(x_i), logo

S(xi)=(xix0)(xixn)xnx0f(xi)=f(xi).S(x_i) = \frac{(x_i - x_0) - (x_i - x_n)}{x_n - x_0}\, f(x_i) = f(x_i) .

9. Indução sobre o número de nós. Um nó: o interpolante é a constante f(x0)=f[x0]f(x_0) = f[x_0]. Suponha a afirmação para kk nós e seja SS interpolando em x0,,xkx_0, \dots, x_k; pela unicidade (questão 4), SS é dado pelo lema de Aitken a partir de RR (nós x0,,xk1x_0, \dots, x_{k-1}) e QQ (nós x1,,xkx_1, \dots, x_k). O coeficiente de XkX^{k} em SS é

[Xk1]Q[Xk1]Rxkx0=f[x1,,xk]f[x0,,xk1]xkx0=f[x0,,xk]\frac{[X^{k-1}]\,Q - [X^{k-1}]\,R}{x_k - x_0} = \frac{f[x_1, \dots, x_k] - f[x_0, \dots, x_{k-1}]}{x_k - x_0} = f[x_0, \dots, x_k]

pela hipótese de indução e pela recorrência que as define.

10. Seja PkP_k o interpolante de ff em x0,,xkx_0, \dots, x_k. A diferença PkPk1P_k - P_{k-1} tem grau k\leq k e anula-se em x0,,xk1x_0, \dots, x_{k-1}, logo, pelo teorema do fator aplicado kk vezes (Teorema 8.7), ela vale cNkc\,N_k para uma constante cc; comparando os coeficientes de XkX^{k} e usando a questão 9, c=f[x0,,xk]c = f[x_0, \dots, x_k]. Telescopando a partir de P0=f(x0)N0P_0 = f(x_0) N_0 obtém-se a fórmula de Newton. Para a forma fechada, escreva Pk=ikf(xi)LiP_k = \sum_{i \leq k} f(x_i) L_i (Lagrange, nos nós x0,,xkx_0, \dots, x_k) e leia o coeficiente de XkX^{k}: cada LiL_i contribui com 1ji(xixj)\frac{1}{\prod_{j \neq i}(x_i - x_j)}, donde

f[x0,,xk]=i=0kf(xi)ji,jk(xixj).f[x_0, \dots, x_k] = \sum_{i=0}^{k} \frac{f(x_i)}{\prod_{j \neq i,\, j \leq k}(x_i - x_j)} .

O lado direito é invariante por qualquer permutação dos nós, logo a diferença dividida não depende da ordenação deles.

11. Se P=aXm+(graus inferiores)P = a X^m + (\text{graus inferiores}), o teorema binomial dá

ΔP=a((X+1)mXm)+=amXm1+(graus inferiores),\Delta P = a\bigl((X+1)^m - X^m\bigr) + \dots = a\,m\,X^{m-1} + (\text{graus inferiores}),

pois (X+1)mXm=mXm1+(X+1)^m - X^m = m X^{m-1} + \dots e a parte de grau inferior de PP contribui com grau m2\leq m - 2 depois de Δ\Delta (ou seja, com termos de grau m2\leq m-2). Logo degΔP=m1\deg \Delta P = m - 1 com coeficiente dominante mam a. Uma constante ccΔc=cc=0\Delta c = c - c = 0.

12. degBk=k\deg B_k = k: escada, logo uma base de Rn[X]\R_n[X] (Exemplo 18.20). Para ΔBk\Delta B_k (k1k \geq 1), fatore o produto comum:

k!ΔBk=(X+1)X(Xk+2)X(X1)(Xk+1)=X(X1)(Xk+2)[(X+1)(Xk+1)]=kX(X1)(Xk+2),\begin{align*} k!\,\Delta B_k &= (X+1)X\cdots(X-k+2) - X(X-1)\cdots(X-k+1) \\ &= X(X-1)\cdots(X-k+2)\,\bigl[(X+1) - (X-k+1)\bigr] \\ &= k\,X(X-1)\cdots(X-k+2), \end{align*}

logo ΔBk=X(X1)(Xk+2)(k1)!=Bk1\Delta B_k = \frac{X(X-1)\cdots(X-k+2)}{(k-1)!} = B_{k-1}.

13. Escreva P=k=0nckBkP = \sum_{k=0}^{n} c_k B_k (base, questão 12). Aplique Δj\Delta^{j}: pela questão 12, ΔjP=kjckBkj\Delta^{j} P = \sum_{k \geq j} c_k B_{k-j}. Avalie em 00: B0(0)=1B_0(0) = 1 e Bm(0)=0B_m(0) = 0 para m1m \geq 1 (o fator XX se anula), logo (ΔjP)(0)=cj\bigl(\Delta^{j}P\bigr)(0) = c_j. Esta é a fórmula das diferenças progressivas.

14. Indução sobre kk. Para k=0k = 0 a identidade lê-se P(0)=P(0)P(0) = P(0). Suponha-a para kk e aplique-a a ΔP\Delta P:

(Δk+1P)(0)=j=0k(1)kj(kj)(P(j+1)P(j)).\bigl(\Delta^{k+1} P\bigr)(0) = \sum_{j=0}^{k} (-1)^{k-j}\binom kj \bigl(P(j+1) - P(j)\bigr).

Junte o coeficiente de P(i)P(i): ele é (1)ki+1(ki1)(1)0(-1)^{k-i+1}\binom k{i-1}\cdot(-1)^{0} na primeira soma (deslocada) e (1)ki(ki)-(-1)^{k-i}\binom ki na segunda — juntos

(1)k+1i((ki1)+(ki))=(1)k+1i(k+1i)(-1)^{k+1-i}\Bigl(\binom k{i-1} + \binom ki\Bigr) = (-1)^{k+1-i}\binom{k+1}i

pela regra de Pascal, que é a identidade na ordem k+1k + 1.

15. Iterando a questão 11 a partir do grau nn, com coeficiente dominante ana_n: depois de um Δ\Delta, grau n1n-1 e coeficiente dominante nann a_n; depois de dois, n(n1)ann(n-1)a_n; depois de nn passos, grau 00 e valor n(n1)1an=n!ann(n-1)\cdots 1\, a_n = n!\,a_n, uma constante. Mais um Δ\Delta a mata: Δn+1P=0\Delta^{n+1}P = 0.

16. Se mkm \geq k: Bk(m)=(mk)NB_k(m) = \binom mk \in \N. Se 0m<k0 \leq m < k: um fator de m(m1)(mk+1)m(m-1)\cdots(m-k+1) é nulo, logo Bk(m)=0B_k(m) = 0. Se m=qm = -q com q1q \geq 1:

Bk(q)=(q)(q1)(qk+1)k!=(1)kq(q+1)(q+k1)k!=(1)k(q+k1k),B_k(-q) = \frac{(-q)(-q-1)\cdots(-q-k+1)}{k!} = (-1)^k\,\frac{q(q+1)\cdots(q+k-1)}{k!} = (-1)^k \binom{q+k-1}{k},

um inteiro. Logo todo BkB_k leva Z\Z em Z\Z.

17. (\Leftarrow) Se P=kckBkP = \sum_k c_k B_k com ckZc_k \in \Z, então para mZm \in \Z, P(m)=kckBk(m)ZP(m) = \sum_k c_k B_k(m) \in \Z pela questão 16. (\Rightarrow) Se PP é a valores inteiros, as suas coordenadas são ck=(ΔkP)(0)=j=0k(1)kj(kj)P(j)c_k = \bigl(\Delta^k P\bigr)(0) = \sum_{j=0}^k (-1)^{k-j}\binom kj P(j) (questões 13 e 14), uma combinação inteira dos inteiros P(0),,P(k)P(0), \dots, P(k). Esta é a caracterização de Pólya dos polinômios a valores inteiros.

18. Ponha Q(X)=P(X+a)Q(X) = P(X + a), um polinômio de grau n\leq n com Q(0),Q(1),,Q(n)ZQ(0), Q(1), \dots, Q(n) \in \Z. As suas coordenadas em (Bk)kn(B_k)_{k \leq n} são ck=jk(1)kj(kj)Q(j)Zc_k = \sum_{j \leq k}(-1)^{k-j}\binom kj Q(j) \in \Z (a questão 14 usa apenas os valores em 0,,kn0, \dots, k \leq n). Pela questão 17 (\Leftarrow), QQ é a valores inteiros em todo o Z\Z, e portanto P(X)=Q(Xa)P(X) = Q(X - a) também.

19. Um produto de kk inteiros consecutivos é m(m1)(mk+1)=k!Bk(m)m(m-1) \cdots(m-k+1) = k!\,B_k(m) para algum mZm \in \Z, e Bk(m)ZB_k(m) \in \Z pela questão 16: o produto é divisível por k!k!.

20. Valores de P=X(X+1)(2X+1)6P = \frac{X(X+1)(2X+1)}{6} em 0,1,2,30,1,2,3: 0,1,5,140, 1, 5, 14. Tabela de diferenças: linha Δ\Delta: 1,4,91, 4, 9; linha Δ2\Delta^2: 3,53, 5; linha Δ3\Delta^3: 22. Logo, pela questão 13,

P=0B0+1B1+3B2+2B3,P = 0\cdot B_0 + 1\cdot B_1 + 3\,B_2 + 2\,B_3 ,

com coordenadas inteiras: PP é a valores inteiros (questão 17), ao passo que os seus coeficientes monomiais 13,12,16\frac13, \frac12, \frac16 não são inteiros. Cálculo direto:

ΔP=(X+1)(X+2)(2X+3)X(X+1)(2X+1)6=(X+1)[(X+2)(2X+3)X(2X+1)]6=(X+1)(6X+6)6=(X+1)2.\begin{align*} \Delta P &= \frac{(X+1)(X+2)(2X+3) - X(X+1)(2X+1)}{6} \\ &= \frac{(X+1)\bigl[(X+2)(2X+3) - X(2X+1)\bigr]}{6} = \frac{(X+1)(6X+6)}{6} = (X+1)^2 . \end{align*}

Telescopando P(m)=j=0m1ΔP(j)=j=1mj2P(m) = \sum_{j=0}^{m-1}\Delta P(j) = \sum_{j=1}^{m} j^2 (com P(0)=0P(0) = 0): a fórmula da soma dos quadrados.

21. Os valores 2i2^i em i=0,,ni = 0, \dots, n têm tabela de diferenças constantemente igual a 11 na borda esquerda: Δk\Delta^k da sequência (2i)(2^i) é de novo (2i)(2^i) (pois 2i+12i=2i2^{i+1} - 2^i = 2^i), logo (ΔkP)(0)=20=1\bigl(\Delta^k P\bigr)(0) = 2^0 = 1 para todo knk \leq n, e P=B0+B1++BnP = B_0 + B_1 + \dots + B_n pela questão 13. Então

P(n+1)=k=0n(n+1k)=2n+1(n+1n+1)=2n+112n+1:P(n+1) = \sum_{k=0}^{n}\binom{n+1}{k} = 2^{n+1} - \binom{n+1}{n+1} = 2^{n+1} - 1 \neq 2^{n+1}:

o padrão quebra no primeiro ponto não controlado.

22. Pela questão 12, Bk=ΔBk+1B_k = \Delta B_{k+1}, logo

j=0m1Bk(j)=j=0m1(Bk+1(j+1)Bk+1(j))=Bk+1(m)Bk+1(0)=Bk+1(m).\sum_{j=0}^{m-1} B_k(j) = \sum_{j=0}^{m-1}\bigl(B_{k+1}(j+1) - B_{k+1}(j)\bigr) = B_{k+1}(m) - B_{k+1}(0) = B_{k+1}(m).

Para j<kj < k os termos Bk(j)B_k(j) se anulam, logo a soma começa de fato em j=kj = k: j=km1(jk)=(mk+1)\sum_{j=k}^{m-1}\binom jk = \binom m{k+1}, a identidade do taco de hóquei.

23. As tabelas de diferenças (ou o desenvolvimento direto) dão

X2=B1+2B2,X3=B1+6B2+6B3X^2 = B_1 + 2 B_2, \qquad X^3 = B_1 + 6 B_2 + 6 B_3

(verificação: B1+2B2=X+X(X1)=X2B_1 + 2B_2 = X + X(X-1) = X^2; em X=1,2,3X = 1, 2, 3 o segundo dá 1,8,271, 8, 27). A questão 22 fornece então

j=0m1j2=B2(m)+2B3(m)=(m2)+2(m3)=m(m1)(2m1)6,\sum_{j=0}^{m-1} j^2 = B_2(m) + 2B_3(m) = \binom m2 + 2\binom m3 = \frac{m(m-1)(2m-1)}{6},
j=0m1j3=B2(m)+6B3(m)+6B4(m)=(m2)+6(m3)+6(m4).\sum_{j=0}^{m-1} j^3 = B_2(m) + 6B_3(m) + 6B_4(m) = \binom m2 + 6\binom m3 + 6\binom m4 .

Desenvolvendo a última expressão: (m2)+6(m3)+6(m4)=m(m1)2[1+2(m2)+(m2)(m3)2]=m2(m1)24=(m2)2\binom m2 + 6\binom m3 + 6\binom m4 = \frac{m(m-1)}{2}\bigl[1 + 2(m-2) + \frac{(m-2)(m-3)}{2}\bigr] = \frac{m^2(m-1)^2}{4} = \binom m2^2. Substituindo mm por m+1m + 1: 13++m3=(m(m+1)2)2=(1++m)21^3 + \dots + m^3 = \bigl(\frac{m(m+1)}2\bigr)^2 = (1 + \dots + m)^2, a identidade de Nicômaco.

24. Nós 0,1,20, 1, 2, polinômio X2X^2. Base monomial (1,X,X2)(1, X, X^2): coordenadas (0,0,1)(0, 0, 1). Base de Lagrange: as coordenadas são os valores (0,1,4)(0, 1, 4) (questão 4). Base de Newton (1,X,X(X1))(1, X, X(X-1)): diferenças divididas f[0]=0f[0] = 0, f[0,1]=1f[0,1] = 1, f[0,1,2]=312=1f[0,1,2] = \frac{3 - 1}{2} = 1 (questão 9), logo coordenadas (0,1,1)(0, 1, 1) — de fato X+X(X1)=X2X + X(X-1) = X^2. Três bases, três vetores de coordenadas, um só polinômio.

25. (i) A questão 4 é automática porque (Li)(L_i) é uma base: existência e unicidade da interpolação são exatamente existência e unicidade das coordenadas. (ii) A base de Newton é uma escada, de modo que as coordenadas se calculam por divisões sucessivas — cada novo nó acrescenta um termo sem perturbar os anteriores — ao passo que as coordenadas de PP na base de Lagrange são os valores P(xi)P(x_i), disponíveis sem cálculo algum. (iii) As duas bases são livres pelos mesmos dois critérios da Proposição 18.19: graus distintos para Newton, avaliação nos nós para Lagrange. (iv) O teorema de Pólya diz que “P(Z)ZP(\Z) \subseteq \Z”, uma propriedade dos valores, equivale à integralidade das coordenadas na base (Bk)(B_k) — a aritmética de um polinômio só se torna visível na base adaptada à pergunta.