Projeção (álgebra linear)

Em álgebra linear e análise funcional, uma projeção é uma transformação linear $P$ de um espaço vetorial em si mesmo, de modo que $P^{2}=P$, ou seja, sempre que $P$ é aplicado duas vezes a algum vetor, o resultado é o mesmo que se tivesse sido aplicado uma única vez (uma propriedade conhecida como idempotência). Embora abstrata, esta definição de "projeção" formaliza e generaliza adequadamente a ideia de projeção gráfica. Também se pode considerar o efeito de uma projeção em um objeto geométrico, examinando o efeito que a projeção tem nos pontos do objeto.

A transformação P é a projeção ortogonal na linha m

Definições

Uma projeção em um espaço vetorial $V$ é um operador linear $P\colon V\to V$ tal que $P^{2}=P$.

Quando o espaço $V$ é munido de um produto interno e é completo (ou seja, quando $V$ é um espaço de Hilbert), o conceito de ortogonalidade pode ser usado. Uma projeção $P$ em um espaço Hilbert $V$ é chamada de projeção ortogonal se satisfaz $\langle P\mathbf {x} ,\mathbf {y} \rangle =\langle \mathbf {x} ,P\mathbf {y} \rangle$ para quaisquer $\mathbf {x} ,\mathbf {y} \in V$. Uma projeção em um espaço de Hilbert que não é ortogonal é chamada de projeção oblíqua.

Matriz de projeção

• No caso de dimensões finitas, uma matriz quadrada $P$ é chamada de matriz de projeção se for igual a seu quadrado, isto é, se $P^{2}=P$.
• Uma matriz quadrada $P$ é chamada de matriz de projeção ortogonal se $P^{2}=P=P^{\mathrm {T} }$ para uma matriz real, e respectivamente $P^{2}=P=P^{*}$ para uma matriz complexa, onde $P^{\mathrm {T} }$ denota a transposição de $P$ e $P^{*}$ denota a transposição adjunta, ou hermitiana, de $P$.
• Uma matriz de projeção que não é ortogonal é chamada de matriz de projeção oblíqua.

Os autovalores de uma matriz de projeção devem ser 0 ou 1.

Exemplos

Projeção ortogonal

Por exemplo, a função que leva o ponto $(x,y,z)$ do espaço tridimensional $\mathbb {R} ^{3}$ ao ponto $(x,y,0)$ é uma projeção ortogonal no plano xy. Esta função pode ser representada pela matriz

$P={\begin{bmatrix}1&0&0\\0&1&0\\0&0&0\end{bmatrix}}.$

A ação desta matriz em um vetor arbitrário pode ser expressa por

$P{\begin{bmatrix}x\\y\\z\end{bmatrix}}={\begin{bmatrix}x\\y\\0\end{bmatrix}}.$

Para ver que $P$ é de fato uma projeção, ou seja, que vale a relação $P=P^{2}$, podemos calcular

$P^{2}{\begin{bmatrix}x\\y\\z\end{bmatrix}}=P{\begin{bmatrix}x\\y\\0\end{bmatrix}}={\begin{bmatrix}x\\y\\0\end{bmatrix}}=P{\begin{bmatrix}x\\y\\z\end{bmatrix}}$.

Observando também que $P^{\mathrm {T} }=P$, vemos que a projeção é uma projeção ortogonal.

Projeção oblíqua

Um exemplo simples de uma projeção não-ortogonal (oblíqua) pode ser dado por

$P={\begin{bmatrix}0&0\\\alpha &1\end{bmatrix}}.$

Pela multiplicação de matrizes, vê-se logo que

$P^{2}={\begin{bmatrix}0&0\\\alpha &1\end{bmatrix}}{\begin{bmatrix}0&0\\\alpha &1\end{bmatrix}}={\begin{bmatrix}0&0\\\alpha &1\end{bmatrix}}=P.$

provando que $P$ é de fato uma projeção. Porém, a projeção só será ortogonal se, e somente se, tivermos $\alpha =0$; pois somente para este caso será válida a relação $P^{\mathrm {T} }=P$.

Propriedades e classificação

A transformação T é a projeção ao longo de k, sobrejetivamente em m. A imagem de t é m, e seu núcleo é k.

Idempotência

Por definição, uma projeção $P$ é idempotente (ou seja, $P^{2}=P$)

Complementariedade da imagem e do núcleo

Seja $W$ um espaço vetorial de dimensão finita e $P$ uma projeção em $W$. Suponha que os subespaços $U$ e $V$ são o imagem e o núcleo de $P$, respectivamente. Então $P$ tem as seguintes propriedades:

1. $P$ é o operador identidade em $U$ $$\forall \mathbf {x} \in U\colon P\mathbf {x} =\mathbf {x}$$
2. Temos a soma direta $W=U\oplus V$. Cada vetor $\mathbf {x} \in W$ pode ser decomposto unicamente como $\mathbf {x} =\mathbf {u} +\mathbf {v}$ com $\mathbf {u} =P\mathbf {x}$ e $\mathbf {v} =\mathbf {x} -P\mathbf {x} =\left(I-P\right)\mathbf {x}$, onde $\mathbf {u} \in U,\mathbf {v} \in V$ e $I$ é o operador identidade.

A imagem e o núcleo de uma projeção são complementares, assim como $P$ e $Q=I-P$. O operador $Q$ também é uma projeção, visto que a imagem e o núcleo de $P$ são o núcleo e imagem de $Q$, e vice-versa.

Espectro

Em espaços vetoriais de dimensão infinita, o espectro de uma projeção está contido em $\{0,1\}$, pois

$(\lambda I-P)^{-1}={\frac {1}{\lambda }}I+{\frac {1}{\lambda (\lambda -1)}}P.$

Apenas 0 ou 1 podem ser autovalores de uma projeção. Isso significa que uma projeção ortogonal $P$ é sempre uma matriz positiva semi-definida. Em geral, os autoespaços correspondentes são o núcleo e a imagem da projeção, respectivamente. A decomposição de um espaço vetorial em somas diretas não é única. Portanto, dado um subespaço $V$, pode haver muitas projeções cuja imagem (ou núcleo) é $V$.

Se a projeção é não-trivial, ela possui um polinômio mínimo $x+y$, que se fatora em raízes distintas, e assim $P$ é diagonalizável.

Produto de projeções

O produto de projeções não é, em geral, uma projeção, mesmo que ambas sejam ortogonais. Se duas projeções comutam, então seu produto é uma projeção, mas o inverso não é verdadeiro: o produto de duas projeções que não comutam pode sim ser uma projeção.

Se duas projeções ortogonais comutam, então seu produto é uma projeção ortogonal. Se o produto de duas projeções ortogonais é uma projeção ortogonal, então as duas projeções ortogonais comutam (de maneira mais geral: dois endomorfismos autoadjuntos comutam se, e somente se, seu produto é também autoadjunto).

Projeções ortogonais

Quando o espaço vetorial $W$ tem um produto interno e é completo (é um espaço de Hilbert) o conceito de ortogonalidade pode ser usado. Uma projeção ortogonal é uma projeção para a qual o intervalo $U$ e o espaço nulo $V$ são subespaços ortogonais. Assim, para cada $\mathbf {x}$ e $\mathbf {y}$ em $W$, $\langle P\mathbf {x} ,(\mathbf {y} -P\mathbf {y} )\rangle =\langle (\mathbf {x} -P\mathbf {x} ),P\mathbf {y} \rangle =0$. Equivalentemente:

$\langle \mathbf {x} ,P\mathbf {y} \rangle =\langle P\mathbf {x} ,P\mathbf {y} \rangle =\langle P\mathbf {x} ,\mathbf {y} \rangle .$

Uma projeção é ortogonal se, e somente se, for autoadjunta. Usando as propriedades autoadjunta e idempotente de $P$, para qualquer $\mathbf {x}$ e $\mathbf {y}$ em $W$ temos $P\mathbf {x} \in U$, $\mathbf {y} -P\mathbf {y} \in V$, e

$\langle P\mathbf {x} ,\mathbf {y} -P\mathbf {y} \rangle =\langle P^{2}\mathbf {x} ,\mathbf {y} -P\mathbf {y} \rangle =\langle P\mathbf {x} ,P\left(I-P\right)\mathbf {y} \rangle =\langle P\mathbf {x} ,\left(P-P^{2}\right)\mathbf {y} \rangle =0$

onde $\langle \cdot ,\cdot \rangle$ é o produto interno associado a $W$. Portanto, $P$ e $I-P$ são projeções ortogonais. A volta, ou seja, se $P$ é ortogonal, então é auto-adjunta, segue de

$\langle \mathbf {x} ,P\mathbf {y} \rangle =\langle P\mathbf {x} ,\mathbf {y} \rangle =\langle \mathbf {x} ,P^{*}\mathbf {y} \rangle$

para qualquer $x$ e $y$ em $W$ ; portanto $P=P^{*}$.

Propriedades e casos especiais

Uma projeção ortogonal é um operador limitado. Isso ocorre porque para cada $\mathbf {v}$ no espaço vetorial temos, pela desigualdade de Cauchy-Schwarz :

$\left\|P\mathbf {v} \right\|^{2}=\langle P\mathbf {v} ,P\mathbf {v} \rangle =\langle P\mathbf {v} ,\mathbf {v} \rangle \leq \left\|P\mathbf {v} \right\|\cdot \left\|\mathbf {v} \right\|$

Assim, $\left\|P\mathbf {v} \right\|\leq \left\|\mathbf {v} \right\|$.

Para espaços vetoriais, reais ou complexos, de dimensão finita, o produto interno padrão pode ser substituído por $\langle \cdot ,\cdot \rangle$.

Projeções oblíquas

O termo projeção oblíqua é algumas vezes usado para se referir a projeções não-ortogonais. Essas projeções também são usadas para representar figuras espaciais em desenhos bidimensionais, embora não com tanta frequência quanto as projeções ortogonais. Enquanto o cálculo do ajuste em uma regressão de mínimos quadrados requer uma projeção ortogonal, o cálculo do valor ajustado de uma regressão de variáveis instrumentais requer uma projeção oblíqua.

As projeções são definidas por seu núcleo e os vetores da base usados para caracterizar sua imagem (que é o complemento do núcleo). Quando os vetores da base são ortogonais ao núcleo, a projeção é uma projeção ortogonal. Quando os vetores de base não são ortogonais ao espaço nulo, a projeção é uma projeção oblíqua. Sejam os vetores $\mathbf {u} _{1},\ldots ,\mathbf {u} _{k}$ uma base para a imagem da projeção, e monte esses vetores em uma matriz $A$, de ordem $n\times k$. A imagem e o núcleo são espaços complementares, logo o núcleo tem dimensão $n-k$. Segue que o complemento ortogonal do espaço nulo tem dimensão $k$. Seja $\mathbf {v} _{1},\ldots ,\mathbf {v} _{k}$ uma base para o complemento ortogonal do núcleo da projeção, e monte uma matriz $B$ com esses vetores. Então a projeção é definida por

$P=A\left(B^{\mathsf {T}}A\right)^{-1}B^{\mathsf {T}}.$

Esta expressão generaliza a fórmula para projeções ortogonais dada anteriormente.

Encontrando a projeção com um produto interno

Seja $V$ um espaço vetorial (neste exemplo, um plano) gerado por vetores ortogonais $\mathbf {u} _{1},\mathbf {u} _{2},\dots ,\mathbf {u} _{p}$. Seja $y$ um vetor qualquer. Podemos definir uma projeção de $\mathbf {y}$ para $V$ como

$\operatorname {proj} _{V}\mathbf {y} ={\frac {\mathbf {y} \cdot \mathbf {u} ^{i}}{\mathbf {u} ^{i}\cdot \mathbf {u} ^{i}}}\mathbf {u} ^{i}$

na qual índices repetidos são somados (seguindo a convenção de Einstein). O vetor $\mathbf {y}$ pode ser escrito como uma soma ortogonal de modo que $\mathbf {y} =\operatorname {proj} _{V}\mathbf {y} +\mathbf {z}$. O operador $\operatorname {proj} _{V}\mathbf {y}$ é às vezes denotado como ${\hat {\mathbf {y} }}$. Há um teorema em álgebra linear que afirma que este $\mathbf {z}$ é a distância mais curta de $\mathbf {y}$ para $V$, e é frequentemente usado em áreas como aprendizado de máquina.

y está sendo projetado no espaço vetorial V

Formas canônicas

Qualquer projeção $P=P^{2}$ em um espaço vetorial de dimensão $d$ é uma matriz diagonalizável, uma vez que seu polinômio mínimo divide $x^{2}-x$, que se separa em fatores lineares distintos. Portanto, existe uma base na qual $P$ tem a forma

$P=I_{r}\oplus 0_{d-r}$

Onde $r$ é a ordem de $P$. Aqui $I_{r}$ é a matriz de identidade de ordem $r$, e $0_{d-r}$ é a matriz nula de ordem $d-r$. Se o espaço vetorial é complexo e munido de um produto interno, então existe uma base ortonormal na qual a matriz de P é dada por

$P={\begin{bmatrix}1&\sigma _{1}\\0&0\end{bmatrix}}\oplus \cdots \oplus {\begin{bmatrix}1&\sigma _{k}\\0&0\end{bmatrix}}\oplus I_{m}\oplus 0_{s}.$

onde $\sigma _{1}\geq \sigma _{2}\geq \dots \geq \sigma _{k}>0$. Os números inteiros $k,s,m$ e os números reais $\sigma _{i}$ estão determinados de forma única. Observe que $2k+s+m=d$. O fator $I_{m}\oplus 0_{s}$ corresponde ao subespaço invariante máximo no qual $P$ atua como uma projeção ortogonal (de modo que o próprio P seja ortogonal se, e somente se, $k=0$) e os blocos $\sigma _{i}$ correspondem às componentes oblíquas.

Projeções em espaços vetoriais normados

Quando o espaço vetorial subjacente $X$ é um espaço vetorial normado (não necessariamente de dimensão finita), questões analíticas, irrelevantes no caso de dimensão finita, precisam ser consideradas. Suponhamos então que $X$ é um espaço de Banach.

Muitos dos resultados algébricos discutidos acima sobrevivem à passagem para este contexto. Uma determinada decomposição de soma direta de $X$ em subespaços complementares ainda especifica uma projeção, e vice-versa. Se $X$ é a soma direta $X=U\oplus V$, então o operador definido por $P(u+v)=u$ ainda é uma projeção com imagem $U$ e núcleo $V$. Também é claro que $P^{2}=P$. Por outro lado, se $P$ é a projeção em $X$, i.e., $P^{2}=P$, então é facilmente verificado que $(1-P)^{2}=(1-P)$. Em outras palavras, $1-P$ também é uma projeção. A relação $P^{2}=P$ implica $1=P+(1-P)$ e $X$ é a soma direta $\operatorname {rg} (P)\oplus \operatorname {rg} (1-P)$.

No entanto, diferentemente do caso de dimensão finita, as projeções não precisam ser contínuas, em geral. Se um subespaço $U$ de $X$ não é fechado na topologia normal, então a projeção em $U$ não é contínua. Em outras palavras, a imagem de uma projeção contínua $P$ deve ser um subespaço fechado. Além disso, o núcleo de uma projeção contínua (na verdade, de um operador linear contínuo qualquer) é fechado. Portanto, uma projeção contínua $P$ leva a uma decomposição de $X$ em dois subespaços fechados complementares: $X=\operatorname {rg} (P)\oplus \ker(P)=\ker(1-P)\oplus \ker(P)$.

O inverso também vale, com uma hipótese adicional. Suponha ser $U$ um subespaço fechado de $X$. Se existe um subespaço fechado $V$ tal que $X=U\oplus V$, então a projeção $P$ com imagem $U$ e núcleo $V$ é contínua. Isso segue do teorema do gráfico fechado. Suponha que $x_{n}\to x$ e $Px_{n}\to y$. É preciso mostrar que $Px=y$. Já que $U$ é fechado e $\{Px_{n}\}\subset U$, $y$ encontra-se em $U$, isto é, $Py=y$. Além disso, $x_{n}-Px_{n}=(I-P)x_{n}\to x-y$. Por $V$ ser fechado e $\{(I-P)x_{n}\}\subset V$, temos que $x-y\in V$, i.e., $P(x-y)=Px-Py=Px-y=0$, o que demonstra a afirmação.

O argumento acima faz uso da hipótese de que ambos $U$ e $V$ são fechados. Em geral, dado um subespaço fechado $U$, não é necessário que exista um subespaço fechado complementar $V$, embora para espaços de Hilbert isso sempre possa ser feito tomando o complemento ortogonal. Para espaços de Banach, um subespaço unidimensional sempre tem um subespaço complementar fechado. Esta é uma consequência imediata do teorema de Hahn-Banach. Seja $U$ ser o subespaço gerado por $u$. Por Hahn-Banach, existe um funcional linear limitado $\varphi$ tal que φ(u) = 1. O operador $P(x)=\varphi (x)u$ satisfaz $P^{2}=P$, logo é uma projeção.A limitação de $\varphi$ implica a continuidade de $P$ e portanto $\ker(P)=\operatorname {rg} (I-P)$ é um subespaço complementar fechado de $U$.

Aplicações e outras considerações

As projeções (ortogonais ou não) desempenham um papel importante em algoritmos para certos problemas de álgebra linear:

• Decomposição QR (ver transformação de Householder e decomposição de Gram – Schmidt);
• Decomposição de valor singular
• Redução para a forma de Hessenberg (a primeira etapa em muitos algoritmos de autovalor)
• Regressão linear
• Elementos projetivos de álgebras de matriz são usados na construção de certos K-grupos na K-teoria de operadores

Como mencionado acima, as projeções são um caso especial de operadores idempotentes. Analiticamente, as projeções ortogonais são generalizações não-comutativas de funções características. Operadores idempotentes são usados na classificação, por exemplo, de álgebras semisimples, enquanto a teoria da medida tem seu início em funções características de conjuntos mensuráveis. Portanto, como se pode imaginar, as projeções são frequentemente encontradas no contexto de álgebras de operadores. Em particular, uma álgebra de von Neumann é gerada por sua rede completa de projeções.

Generalizações

De maneira mais geral, dado um mapa entre espaços vetoriais normados $T\colon V\to W,$ pode-se exigir analogamente que este mapa seja uma isometria no complemento ortogonal do núcleo: que $(\ker T)^{\perp }\to W$ ser uma isometria; em particular, ela deve ser sobrejetora. O caso de uma projeção ortogonal ocorre quando $W$ é um subespaço de $V$. Na geometria Riemanniana, isso é usado na definição de uma submersão Riemanniana.

Ver também

• Matriz de centralização, que é um exemplo de matriz de projeção.
• Ortogonalização
• Subespaço invariável
• Propriedades do traço
• O algoritmo de projeção de Dykstra para calcular a projeção em uma interseção de conjuntos

Referências

• Banerjee, Sudipto; Roy, Anindya (2014), Linear Algebra and Matrix Analysis for Statistics, ISBN 978-1420095388, Texts in Statistical Science 1st ed. , Chapman and Hall/CRC 
• Dunford, N.; Schwartz, J. T. (1958). Linear Operators, Part I: General Theory. [S.l.]: Interscience 
• Meyer, Carl D. (2000). Matrix Analysis and Applied Linear Algebra. [S.l.]: Society for Industrial and Applied Mathematics. ISBN 978-0-89871-454-8 
• MIT Linear Algebra Lecture on Projection Matrices , do MIT OpenCourseWare
• Linear Algebra 15d: The Projection Transformation, de Pavel Grinfeld.
• Tutorial de projeções geométricas planas - um tutorial simples de se seguir que explica os diferentes tipos de projeções geométricas planas.