Corsi di laurea › Ing. Telecomunicazioni › Algebra Lineare e Geometria › 6. Prodotto scalare e ortogonalità
Esercizio 82(Im f)⊥ = Ker f per una matrice simmetrica In questa pagina 4 (a) Im f \operatorname{Im} f Im f , ( Im f ) ⊥ (\operatorname{Im} f)^\perp ( Im f ) ⊥ e Ker f \operatorname{Ker} f Ker f (b) Base ortogonale di Im f \operatorname{Im} f Im f (c) v = u + w v = u + w v = u + w con u ∈ Ker f u \in \operatorname{Ker} f u ∈ Ker f , w ∈ Im f w \in \operatorname{Im} f w ∈ Im f Errori comuni
Testo (Appello del 2 luglio 2021, esercizio 3). Sia f : R 4 → R 4 f: \mathbb{R}^4 \to \mathbb{R}^4 f : R 4 → R 4 la funzione lineare la cui matrice, rispetto alla base canonica, è
A = ( 4 − 2 2 0 − 2 2 − 3 1 2 − 3 5 − 2 0 1 − 2 1 ) . A = \begin{pmatrix} 4 & -2 & 2 & 0 \\ -2 & 2 & -3 & 1 \\ 2 & -3 & 5 & -2 \\ 0 & 1 & -2 & 1 \end{pmatrix}. A = 4 − 2 2 0 − 2 2 − 3 1 2 − 3 5 − 2 0 1 − 2 1 .
(a) Trovare una base di Im f \operatorname{Im} f Im f e una base di ( Im f ) ⊥ (\operatorname{Im} f)^\perp ( Im f ) ⊥ e verificare che ( Im f ) ⊥ = Ker f (\operatorname{Im} f)^\perp = \operatorname{Ker} f ( Im f ) ⊥ = Ker f .
(b) Usando il procedimento di Gram-Schmidt trovare una base ortogonale di Im f \operatorname{Im} f Im f .
(c) Dato il vettore v = ( 1 , 5 , − 3 , 1 ) v = (1, 5, -3, 1) v = ( 1 , 5 , − 3 , 1 ) , trovare u ∈ Ker f u \in \operatorname{Ker} f u ∈ Ker f e w ∈ Im f w \in \operatorname{Im} f w ∈ Im f tali che v = u + w v = u + w v = u + w .
Teoria usata: Complemento ortogonale e proiezioni ortogonaliL'ortogonale U⊥ di un sottospazio è un sottospazio di dimensione n − dim U, e R^n = U ⊕ U⊥; ogni vettore si scompone in proiezione su U più componente ortogonale; la proiezione è il punto di U più vicino e si calcola con un sistema o con la matrice di proiezione A(AᵀA)⁻¹Aᵀ.Complemento ortogonale e proiezioni ortogonali → , Basi ortonormali e Gram-SchmidtUna base ortonormale è fatta di vettori di norma 1 a due a due ortogonali: le coordinate si calcolano con prodotti scalari e le proiezioni con una formula diretta. Il procedimento di Gram-Schmidt trasforma una base qualsiasi in una base ortogonale (e poi ortonormale) dello stesso sottospazio.Basi ortonormali e Gram-Schmidt → , Nucleo e immagineIl nucleo (vettori mandati in 0) e l'immagine (vettori raggiunti) di una funzione lineare sono sottospazi; f è iniettiva se e solo se Ker f = {0}; dim Ker f + dim Im f = dim V (nullità + rango); l'antimmagine di un vettore è una soluzione particolare più il nucleo.Nucleo e immagine → , Eliminazione di GaussCon tre operazioni elementari sulle righe (scambio, moltiplicazione per uno scalare non nullo, somma di un multiplo di un'altra riga) ogni matrice si riduce a scala senza cambiare il rango; serve a calcolare ranghi, risolvere sistemi, trovare relazioni di dipendenza e matrici che riducono a scala.Eliminazione di Gauss → .
(a) Im f \operatorname{Im} f Im f , ( Im f ) ⊥ (\operatorname{Im} f)^\perp ( Im f ) ⊥ e Ker f \operatorname{Ker} f Ker f
Rango con Gauss
A → R 3 → R 3 − 1 2 R 1 R 2 → R 2 + 1 2 R 1 ( 4 − 2 2 0 0 1 − 2 1 0 − 2 4 − 2 0 1 − 2 1 ) → R 4 → R 4 − R 2 R 3 → R 3 + 2 R 2 ( 4 − 2 2 0 0 1 − 2 1 0 0 0 0 0 0 0 0 ) A \xrightarrow[\;R_3 \to R_3 - \frac12 R_1\;]{R_2 \to R_2 + \frac12 R_1} \begin{pmatrix} 4 & -2 & 2 & 0 \\ 0 & 1 & -2 & 1 \\ 0 & -2 & 4 & -2 \\ 0 & 1 & -2 & 1 \end{pmatrix} \xrightarrow[\;R_4 \to R_4 - R_2\;]{R_3 \to R_3 + 2R_2} \begin{pmatrix} 4 & -2 & 2 & 0 \\ 0 & 1 & -2 & 1 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{pmatrix} A R 2 → R 2 + 2 1 R 1 R 3 → R 3 − 2 1 R 1 4 0 0 0 − 2 1 − 2 1 2 − 2 4 − 2 0 1 − 2 1 R 3 → R 3 + 2 R 2 R 4 → R 4 − R 2 4 0 0 0 − 2 1 0 0 2 − 2 0 0 0 1 0 0
Due pivot (nella colonna 1 e nella colonna 2): rango ( A ) = 2 \operatorname{rango}(A) = 2 rango ( A ) = 2 , quindi dim Im f = 2 \dim \operatorname{Im} f = 2 dim Im f = 2 e, per il teorema nullità + rangodim Ker f + dim Im f = dim del dominio = 4.Nucleo e immagine → , dim Ker f = 4 − 2 = 2 \dim \operatorname{Ker} f = 4 - 2 = 2 dim Ker f = 4 − 2 = 2 .
Base di Im f \operatorname{Im} f Im f
L'immagine è generata dalle colonne di A A A ; le colonne dove stanno i pivot (la prima e la seconda) sono indipendenti e formano una base:
v 1 = ( 4 , − 2 , 2 , 0 ) , v 2 = ( − 2 , 2 , − 3 , 1 ) . v_1 = (4, -2, 2, 0), \qquad v_2 = (-2, 2, -3, 1). v 1 = ( 4 , − 2 , 2 , 0 ) , v 2 = ( − 2 , 2 , − 3 , 1 ) .
Base di ( Im f ) ⊥ (\operatorname{Im} f)^\perp ( Im f ) ⊥
x ∈ ( Im f ) ⊥ x \in (\operatorname{Im} f)^\perp x ∈ ( Im f ) ⊥ se e solo se x ⋅ v 1 = 0 x \cdot v_1 = 0 x ⋅ v 1 = 0 e x ⋅ v 2 = 0 x \cdot v_2 = 0 x ⋅ v 2 = 0 :
{ 4 x 1 − 2 x 2 + 2 x 3 = 0 − 2 x 1 + 2 x 2 − 3 x 3 + x 4 = 0 \begin{cases} 4x_1 - 2x_2 + 2x_3 = 0 \\ -2x_1 + 2x_2 - 3x_3 + x_4 = 0 \end{cases} { 4 x 1 − 2 x 2 + 2 x 3 = 0 − 2 x 1 + 2 x 2 − 3 x 3 + x 4 = 0
Dalla prima (divisa per 2 2 2 ): x 2 = 2 x 1 + x 3 x_2 = 2x_1 + x_3 x 2 = 2 x 1 + x 3 . Nella seconda: − 2 x 1 + 4 x 1 + 2 x 3 − 3 x 3 + x 4 = 0 -2x_1 + 4x_1 + 2x_3 - 3x_3 + x_4 = 0 − 2 x 1 + 4 x 1 + 2 x 3 − 3 x 3 + x 4 = 0 , cioè x 4 = − 2 x 1 + x 3 x_4 = -2x_1 + x_3 x 4 = − 2 x 1 + x 3 . Liberi x 1 , x 3 x_1, x_3 x 1 , x 3 :
( x 1 , x 3 ) = ( 1 , 0 ) (x_1, x_3) = (1, 0) ( x 1 , x 3 ) = ( 1 , 0 ) : w 1 = ( 1 , 2 , 0 , − 2 ) w_1 = (1, 2, 0, -2) w 1 = ( 1 , 2 , 0 , − 2 ) ;
( x 1 , x 3 ) = ( 0 , 1 ) (x_1, x_3) = (0, 1) ( x 1 , x 3 ) = ( 0 , 1 ) : w 2 = ( 0 , 1 , 1 , 1 ) w_2 = (0, 1, 1, 1) w 2 = ( 0 , 1 , 1 , 1 ) .
Verifica di ( Im f ) ⊥ = Ker f (\operatorname{Im} f)^\perp = \operatorname{Ker} f ( Im f ) ⊥ = Ker f
Due sottospazi della stessa dimensione, uno contenuto nell'altro, coincidono (sottospazio di uguale dimensioneSe W ⊆ U e dim W = dim U (finita), allora W = U.Dimensione → ). Basta allora mostrare ( Im f ) ⊥ ⊆ Ker f (\operatorname{Im} f)^\perp \subseteq \operatorname{Ker} f ( Im f ) ⊥ ⊆ Ker f , cioè A w 1 = 0 Aw_1 = 0 A w 1 = 0 e A w 2 = 0 Aw_2 = 0 A w 2 = 0 :
A w 1 = ( 4 − 4 + 0 + 0 − 2 + 4 + 0 − 2 2 − 6 + 0 + 4 0 + 2 + 0 − 2 ) = 0 , A w 2 = ( 0 − 2 + 2 + 0 0 + 2 − 3 + 1 0 − 3 + 5 − 2 0 + 1 − 2 + 1 ) = 0. A w_1 = \begin{pmatrix} 4 - 4 + 0 + 0 \\ -2 + 4 + 0 - 2 \\ 2 - 6 + 0 + 4 \\ 0 + 2 + 0 - 2 \end{pmatrix} = 0, \qquad A w_2 = \begin{pmatrix} 0 - 2 + 2 + 0 \\ 0 + 2 - 3 + 1 \\ 0 - 3 + 5 - 2 \\ 0 + 1 - 2 + 1 \end{pmatrix} = 0 . A w 1 = 4 − 4 + 0 + 0 − 2 + 4 + 0 − 2 2 − 6 + 0 + 4 0 + 2 + 0 − 2 = 0 , A w 2 = 0 − 2 + 2 + 0 0 + 2 − 3 + 1 0 − 3 + 5 − 2 0 + 1 − 2 + 1 = 0.
Quindi ( Im f ) ⊥ ⊆ Ker f (\operatorname{Im} f)^\perp \subseteq \operatorname{Ker} f ( Im f ) ⊥ ⊆ Ker f e, siccome entrambi hanno dimensione 2 2 2 , ( Im f ) ⊥ = Ker f (\operatorname{Im} f)^\perp = \operatorname{Ker} f ( Im f ) ⊥ = Ker f ✓.
Perché succede? A A A è simmetrica (A T = A A^T = A A T = A ). Per una matrice qualunque vale Ker A = ( spazio delle righe ) ⊥ \operatorname{Ker} A = (\text{spazio delle righe})^\perp Ker A = ( spazio delle righe ) ⊥ : infatti A x = 0 Ax = 0 A x = 0 significa che x x x è ortogonale a ogni riga. Se A A A è simmetrica le righe coincidono con le colonne, che generano Im f \operatorname{Im} f Im f : da qui Ker f = ( Im f ) ⊥ \operatorname{Ker} f = (\operatorname{Im} f)^\perp Ker f = ( Im f ) ⊥ . È un'anticipazione del teorema spettraleUna funzione lineare è simmetrica se f(v)·w = v·f(w); in una base ortonormale ha matrice simmetrica. Teorema spettrale: f è simmetrica se e solo se esiste una base ortonormale di autovettori, cioè A simmetrica ⇔ PᵀAP diagonale con P ortogonale. Applicato alle forme quadratiche, permette di scriverle come somma di quadrati con gli autovalori come coefficienti.Teorema spettrale e forme quadratiche → .
(b) Base ortogonale di Im f \operatorname{Im} f Im f
v 1 ⋅ v 2 = − 8 − 4 − 6 + 0 = − 18 ≠ 0 v_1 \cdot v_2 = -8 - 4 - 6 + 0 = -18 \ne 0 v 1 ⋅ v 2 = − 8 − 4 − 6 + 0 = − 18 = 0 . Gram-Schmidt: v 1 ′ = v 1 v_1' = v_1 v 1 ′ = v 1 , v 2 ′ = v 2 + α v 1 v_2' = v_2 + \alpha v_1 v 2 ′ = v 2 + α v 1 con
α = − v 1 ⋅ v 2 v 1 ⋅ v 1 = − − 18 16 + 4 + 4 + 0 = 18 24 = 3 4 , \alpha = -\frac{v_1 \cdot v_2}{v_1 \cdot v_1} = -\frac{-18}{16 + 4 + 4 + 0} = \frac{18}{24} = \frac34, α = − v 1 ⋅ v 1 v 1 ⋅ v 2 = − 16 + 4 + 4 + 0 − 18 = 24 18 = 4 3 ,
v 2 ′ = ( − 2 , 2 , − 3 , 1 ) + 3 4 ( 4 , − 2 , 2 , 0 ) = ( 1 , 1 2 , − 3 2 , 1 ) . v_2' = (-2, 2, -3, 1) + \frac34 (4, -2, 2, 0) = \left(1, \frac12, -\frac32, 1\right). v 2 ′ = ( − 2 , 2 , − 3 , 1 ) + 4 3 ( 4 , − 2 , 2 , 0 ) = ( 1 , 2 1 , − 2 3 , 1 ) .
Moltiplicando per 2 2 2 (non cambia l'ortogonalità): ( 2 , 1 , − 3 , 2 ) (2, 1, -3, 2) ( 2 , 1 , − 3 , 2 ) .
Controllo: v 1 ⋅ ( 2 , 1 , − 3 , 2 ) = 8 − 2 − 6 + 0 = 0 v_1 \cdot (2, 1, -3, 2) = 8 - 2 - 6 + 0 = 0 v 1 ⋅ ( 2 , 1 , − 3 , 2 ) = 8 − 2 − 6 + 0 = 0 ✓.
Risultato (b): { ( 4 , − 2 , 2 , 0 ) , ( 2 , 1 , − 3 , 2 ) } \{(4, -2, 2, 0),\ (2, 1, -3, 2)\} {( 4 , − 2 , 2 , 0 ) , ( 2 , 1 , − 3 , 2 )} .
(c) v = u + w v = u + w v = u + w con u ∈ Ker f u \in \operatorname{Ker} f u ∈ Ker f , w ∈ Im f w \in \operatorname{Im} f w ∈ Im f
Siccome Ker f = ( Im f ) ⊥ \operatorname{Ker} f = (\operatorname{Im} f)^\perp Ker f = ( Im f ) ⊥ , scrivere v = u + w v = u + w v = u + w con w ∈ Im f w \in \operatorname{Im} f w ∈ Im f e u ∈ ( Im f ) ⊥ u \in (\operatorname{Im} f)^\perp u ∈ ( Im f ) ⊥ significa proprio scomporre v v v in proiezione ortogonale su Im f \operatorname{Im} f Im f più componente ortogonalePer ogni sottospazio W si ha ℝⁿ = W ⊕ W⊥, quindi ogni v si scrive in un solo modo come w + w' con w ∈ W, w' ∈ W⊥.Complemento ortogonale e proiezioni ortogonali → . La decomposizione è unica.
Con la base ortogonale di (b), v 1 ′ = ( 4 , − 2 , 2 , 0 ) v_1' = (4, -2, 2, 0) v 1 ′ = ( 4 , − 2 , 2 , 0 ) e v 2 ′ ′ = ( 2 , 1 , − 3 , 2 ) v_2'' = (2, 1, -3, 2) v 2 ′′ = ( 2 , 1 , − 3 , 2 ) :
v ⋅ v 1 ′ = 4 − 10 − 6 + 0 = − 12 v \cdot v_1' = 4 - 10 - 6 + 0 = -12 v ⋅ v 1 ′ = 4 − 10 − 6 + 0 = − 12 , v 1 ′ ⋅ v 1 ′ = 24 \ v_1' \cdot v_1' = 24 v 1 ′ ⋅ v 1 ′ = 24 : coefficiente − 1 2 -\frac12 − 2 1 ;
v ⋅ v 2 ′ ′ = 2 + 5 + 9 + 2 = 18 v \cdot v_2'' = 2 + 5 + 9 + 2 = 18 v ⋅ v 2 ′′ = 2 + 5 + 9 + 2 = 18 , v 2 ′ ′ ⋅ v 2 ′ ′ = 4 + 1 + 9 + 4 = 18 \ v_2'' \cdot v_2'' = 4 + 1 + 9 + 4 = 18 v 2 ′′ ⋅ v 2 ′′ = 4 + 1 + 9 + 4 = 18 : coefficiente 1 1 1 .
w = − 1 2 ( 4 , − 2 , 2 , 0 ) + ( 2 , 1 , − 3 , 2 ) = ( − 2 , 1 , − 1 , 0 ) + ( 2 , 1 , − 3 , 2 ) = ( 0 , 2 , − 4 , 2 ) , w = -\frac12 (4, -2, 2, 0) + (2, 1, -3, 2) = (-2, 1, -1, 0) + (2, 1, -3, 2) = (0, 2, -4, 2), w = − 2 1 ( 4 , − 2 , 2 , 0 ) + ( 2 , 1 , − 3 , 2 ) = ( − 2 , 1 , − 1 , 0 ) + ( 2 , 1 , − 3 , 2 ) = ( 0 , 2 , − 4 , 2 ) ,
u = v − w = ( 1 , 5 , − 3 , 1 ) − ( 0 , 2 , − 4 , 2 ) = ( 1 , 3 , 1 , − 1 ) . u = v - w = (1, 5, -3, 1) - (0, 2, -4, 2) = (1, 3, 1, -1). u = v − w = ( 1 , 5 , − 3 , 1 ) − ( 0 , 2 , − 4 , 2 ) = ( 1 , 3 , 1 , − 1 ) .
Controlli:
u ∈ Ker f u \in \operatorname{Ker} f u ∈ Ker f : A u = ( 4 − 6 + 2 + 0 , − 2 + 6 − 3 − 1 , 2 − 9 + 5 + 2 , 0 + 3 − 2 − 1 ) = ( 0 , 0 , 0 , 0 ) Au = (4 - 6 + 2 + 0,\ -2 + 6 - 3 - 1,\ 2 - 9 + 5 + 2,\ 0 + 3 - 2 - 1) = (0, 0, 0, 0) A u = ( 4 − 6 + 2 + 0 , − 2 + 6 − 3 − 1 , 2 − 9 + 5 + 2 , 0 + 3 − 2 − 1 ) = ( 0 , 0 , 0 , 0 ) ✓;
w ∈ Im f w \in \operatorname{Im} f w ∈ Im f : w = v 1 + 2 v 2 = ( 4 − 4 , − 2 + 4 , 2 − 6 , 0 + 2 ) = ( 0 , 2 , − 4 , 2 ) w = v_1 + 2v_2 = (4 - 4,\ -2 + 4,\ 2 - 6,\ 0 + 2) = (0, 2, -4, 2) w = v 1 + 2 v 2 = ( 4 − 4 , − 2 + 4 , 2 − 6 , 0 + 2 ) = ( 0 , 2 , − 4 , 2 ) ✓.
Risultato (c): u = ( 1 , 3 , 1 , − 1 ) u = (1, 3, 1, -1) u = ( 1 , 3 , 1 , − 1 ) , w = ( 0 , 2 , − 4 , 2 ) w = (0, 2, -4, 2) w = ( 0 , 2 , − 4 , 2 ) .
Metodo alternativo (quello della soluzione ufficiale). Si pone w = α v 1 + β v 2 w = \alpha v_1 + \beta v_2 w = α v 1 + β v 2 , u = v − w u = v - w u = v − w , e si impone A u = 0 Au = 0 A u = 0 (usando la forma a scala di A A A , che ha le stesse soluzioni e meno righe). Si ottiene un sistema in α , β \alpha, \beta α , β con soluzione α = 1 \alpha = 1 α = 1 , β = 2 \beta = 2 β = 2 . Questo metodo non richiede la base ortogonale e funziona anche quando Ker f \operatorname{Ker} f Ker f e Im f \operatorname{Im} f Im f non sono ortogonali (purché siano in somma diretta).
Errori comuni
Prendere come base dell'immagine le righe non nulle della forma a scala : le righe della forma a scala generano lo spazio delle righe , non quello delle colonne. Per Im f \operatorname{Im} f Im f si prendono le colonne di A A A nelle posizioni dei pivot.
Credere che ( Im f ) ⊥ = Ker f (\operatorname{Im} f)^\perp = \operatorname{Ker} f ( Im f ) ⊥ = Ker f valga sempre : vale per le matrici simmetriche; in generale ( Im f ) ⊥ = Ker ( A T ) (\operatorname{Im} f)^\perp = \operatorname{Ker}(A^T) ( Im f ) ⊥ = Ker ( A T ) .
Precedente Esercizio 81 · base ortogonale, U⊥, proiezione e un sottospazio impossibile Successiva Esercizio 83 · U da un sistema, U⊥, proiezioni su U e U⊥