Esercizio prodotti scalari

Messaggioda Pic143 » 25/01/2020, 14:45

Salve avrei bisogno di una mano per questo esercizio.
Si consideri il prodotto scalare $ phi :R^3xx R^3->R $ dato da:
$ phi (X,Y)=X^tAY. A=( ( 2 , 1 , 0 ),( 1 , 2 , 0 ),( 0 , 0 , 1 ) ) $

1) Verificare che $phi$ è definito positivo.
2) Stabilire se esiste una matrice $N$ tale che $N^t N=A$.

Il punto 1) l'ho svolto tranquillamente: ho trovato gli autovalori di A che risultano essere 1 e 3, entrambi positivi.
Mi servirebbe sapere come svolgere il punto due e il suo significato teorico.
Granzie in anticipo.
Pic143
Starting Member
Starting Member
 
Messaggio: 2 di 48
Iscritto il: 05/02/2018, 11:55

Re: Esercizio prodotti scalari

Messaggioda Bokonon » 25/01/2020, 18:57

Iniziamo dal significato teorico. A è una matrice simmetrica pertanto il teorema spettrale ci garantisce che sia sempre diagonalizzabile (con autovalori reali) e sia sempre possibile trovare una base di autovettori ortonormale tale che $A=QLambdaQ^(-1)=QLambdaQ^(T)$ dove $Lambda$ è una matrice diagonale.

Come sai, se si prende una matrice N di dimensioni a piacere (mxn), allora $N^TN$ è una matrice simmetrica quindi gode non solo delle proprietà di cui sopra ma anche del fatto che sia sempre semidefinita/definita positiva.
La dimostrazione è facile. Preso un qualsiasi vettore x (diverso da quello nullo) allora una matrice $N^TN$ è semidefinita positiva se: $x^TN^TNx=(Nx)^TNx>=0$
Dato che $Nx$ è un vettore, la sua norma al quadrato è sempre positiva o uguale a zero nel caso x appartenga al kernel di N. C.V.D.
Questa proprietà è chiave per una decomposizione estremamente importante ovvero la SVD (singular values decomposition o decomposizione a valori singolari).

Come trovare N (e qua l'ho dedotto da me, visto che non uso libri).
Hai già dimostrato che A è definita positiva (quindi esiste un N che soddisfa i requisiti).
Quindi $det(A)>0$.
Quindi cerchiamo una matrice N di dimensioni (mx3) tale che $N^TN=A$ (che è appunto 3x3)
$det(N^TN)=det(A)>0$ e questo accade se e solo se $m>=3$: tradotto le colonne di N devono essere indipendenti e affinchè ciò accada devono essere composte da vettori di almeno 3 componenti.

$A=N^TN=QLambdaQ^(T) rArr Q^TN^TNQ=Lambda rArr (NQ)^TNQ=Lambda$
Ponendo $NQ=K$ abbiamo che $N=KQ^T$

E ora mi sono chiesto quale matrice K possa risolvere $K^TK=Lambda$.
La risposta è per $K=K^T$, quindi abbiamo $K^2=Lambda rArr K=Sigma$

Come scoprirai quando studierai la decomposizione SVD, $Sigma$ è la matrice diagionale dei valori singolari, che non sono altro che la radice quadrata degli autovalori.

Ora, non hai completato l'esercizio trovando anche $Q$. Quando l'avrai fatto, troverai K e di conseguenza N.
Una soluzione è $ N=( ( -1/sqrt(2) , 1/sqrt(2) , 0 ),( 0 , 0 , 1 ),( sqrt(3/2) , sqrt(3/2) , 0 ) ) $
Avatar utente
Bokonon
Cannot live without
Cannot live without
 
Messaggio: 1987 di 5942
Iscritto il: 25/05/2018, 20:22

Re: Esercizio prodotti scalari

Messaggioda Pic143 » 26/01/2020, 12:00

Sergio ha scritto:
Pic143 ha scritto:Mi servirebbe sapere come svolgere il punto due e il suo significato teorico.

In effetti il punto due è "strano" perché $A=N^TN$ è la scomposizione di Cholesky, legata alle scomposizioni LU e LDL, e anche alla SVD (v. sopra).
Possibile che non hai studiato nessuno di questi argomenti?


No in effetti non ho studiato nessuna di queste. Abbiamo fatto esclusivamente la riduzione di Gauss.
Pic143
Starting Member
Starting Member
 
Messaggio: 3 di 48
Iscritto il: 05/02/2018, 11:55

Re: Esercizio prodotti scalari

Messaggioda Pic143 » 26/01/2020, 12:08

Bokonon ha scritto:Iniziamo dal significato teorico. A è una matrice simmetrica pertanto il teorema spettrale ci garantisce che sia sempre diagonalizzabile (con autovalori reali) e sia sempre possibile trovare una base di autovettori ortonormale tale che $A=QLambdaQ^(-1)=QLambdaQ^(T)$ dove $Lambda$ è una matrice diagonale.

Come sai, se si prende una matrice N di dimensioni a piacere (mxn), allora $N^TN$ è una matrice simmetrica quindi gode non solo delle proprietà di cui sopra ma anche del fatto che sia sempre semidefinita/definita positiva.
La dimostrazione è facile. Preso un qualsiasi vettore x (diverso da quello nullo) allora una matrice $N^TN$ è semidefinita positiva se: $x^TN^TNx=(Nx)^TNx>=0$
Dato che $Nx$ è un vettore, la sua norma al quadrato è sempre positiva o uguale a zero nel caso x appartenga al kernel di N. C.V.D.
Questa proprietà è chiave per una decomposizione estremamente importante ovvero la SVD (singular values decomposition o decomposizione a valori singolari).

Come trovare N (e qua l'ho dedotto da me, visto che non uso libri).
Hai già dimostrato che A è definita positiva (quindi esiste un N che soddisfa i requisiti).
Quindi $det(A)>0$.
Quindi cerchiamo una matrice N di dimensioni (mx3) tale che $N^TN=A$ (che è appunto 3x3)
$det(N^TN)=det(A)>0$ e questo accade se e solo se $m>=3$: tradotto le colonne di N devono essere indipendenti e affinchè ciò accada devono essere composte da vettori di almeno 3 componenti.

$A=N^TN=QLambdaQ^(T) rArr Q^TN^TNQ=Lambda rArr (NQ)^TNQ=Lambda$
Ponendo $NQ=K$ abbiamo che $N=KQ^T$

E ora mi sono chiesto quale matrice K possa risolvere $K^TK=Lambda$.
La risposta è per $K=K^T$, quindi abbiamo $K^2=Lambda rArr K=Sigma$

Come scoprirai quando studierai la decomposizione SVD, $Sigma$ è la matrice diagionale dei valori singolari, che non sono altro che la radice quadrata degli autovalori.

Ora, non hai completato l'esercizio trovando anche $Q$. Quando l'avrai fatto, troverai K e di conseguenza N.
Una soluzione è $ N=( ( -1/sqrt(2) , 1/sqrt(2) , 0 ),( 0 , 0 , 1 ),( sqrt(3/2) , sqrt(3/2) , 0 ) ) $


Grazie mille. L'esercizio chiede solo di stabilirne l'esistenza quindi non mi devo preoccupare di argomenti non trattati. La dimostrazione della sola esistenza è alla mia portata. Grazie ancora.
Pic143
Starting Member
Starting Member
 
Messaggio: 4 di 48
Iscritto il: 05/02/2018, 11:55

Re: Esercizio prodotti scalari

Messaggioda Bokonon » 26/01/2020, 20:43

Prego
Avatar utente
Bokonon
Cannot live without
Cannot live without
 
Messaggio: 2000 di 5942
Iscritto il: 25/05/2018, 20:22


Torna a Geometria e algebra lineare

Chi c’è in linea

Visitano il forum: Nessuno e 1 ospite