Cos'è ElectroYou | Login Iscriviti

ElectroYou - la comunità dei professionisti del mondo elettrico

Metodo minimi quadrati polinomio approssimazione

Analisi, geometria, algebra, topologia...

Moderatori: Foto UtentePietroBaima, Foto UtenteIanero

0
voti

[1] Metodo minimi quadrati polinomio approssimazione

Messaggioda Foto Utentejupiter » 11 lug 2014, 17:32

Ciao a tutti

perché il sistema lineare Ac = y con cui voglio trovare le incognite (i coefficienti) del polinomio approssimante, che utilizza il metodo dei minimi quadrati, diviene A^{^{t}} Ac = A^{^{t}} y ?

Grazie :-)
Avatar utente
Foto Utentejupiter
10 3
 
Messaggi: 34
Iscritto il: 21 dic 2008, 13:49

1
voti

[2] Re: Metodo minimi quadrati polinomio approssimazione

Messaggioda Foto UtenteIanero » 11 lug 2014, 17:51

La soluzione esatta ovviamente non esiste per i coefficienti di un polinomio approssimante.
La soluzione più "vicina", la migliore, si ottiene proiettando il vettore dei termini noti sullo spazio delle colonne della matrice del sistema.
Riesci a farti un'idea del perché?

Prova a iniziare da qui per dimostrare quella relazione, sfruttando le proprietà dei complementi ortogonali.

Stasera continuiamo se non ci sono news da parte tua :D
:shock:
Avatar utente
Foto UtenteIanero
8.069 5 8 11
Master EY
Master EY
 
Messaggi: 4320
Iscritto il: 21 mar 2012, 15:47

1
voti

[3] Re: Metodo minimi quadrati polinomio approssimazione

Messaggioda Foto Utentejupiter » 11 lug 2014, 18:34

se osservo la matrice A ho un sistema n x m mentre se applico la trasposizione il sistema diventa n x n....

:-k
Avatar utente
Foto Utentejupiter
10 3
 
Messaggi: 34
Iscritto il: 21 dic 2008, 13:49

0
voti

[4] Re: Metodo minimi quadrati polinomio approssimazione

Messaggioda Foto UtenteIanero » 12 lug 2014, 11:24

AX=B

Affinchè B sia soluzione, deve appartenere allo spazio delle colonne di A.



Come vedi B nel nostro caso non appartiene allo spazio delle colonne della matrice del sistema, perché non esiste nessun X tale che AX=B (se esistesse B sarebbe nel iperpiano \text{col}(A)).

Adesso è evidente che la soluzione migliore è rappresentata dal vettore proiezione di B su \text{col}(A) perché questo vettore che chiamiamo AZ è il meno "distante" a B tra tutti i vettori possibili appartenenti a \text{col}(A) (è un teorema facilmente dimostrabile).

Allora al lavoro.

AZ=\text{proj}_{\text{col}(A)}B

B-AZ=B-\text{proj}_{\text{col}(A)}B\perp AY

Dove Y è un generico vettore.

AY\bullet \left (B-AZ  \right )=0

\left (AY  \right )^T \left (B-AZ  \right )=0

Y^T A^T \left (B-AZ  \right )=0

Y \bullet  A^T \left (B-AZ  \right )=0

Essendo Y generico vettore di \mathbb{R}^n otteniamo:

A^T \left (B-AZ  \right )=0

A^T B- A^T AZ  \right )=0

\left (A^T A  \right )Z  \right )=\left (A^T B  \right )




se osservo la matrice A ho un sistema n x m mentre se applico la trasposizione il sistema diventa n x n


Ottima osservazione!
Prima il sistema era sovradeterminato e senza soluzione.
Dopo invece... :D
:shock:
Avatar utente
Foto UtenteIanero
8.069 5 8 11
Master EY
Master EY
 
Messaggi: 4320
Iscritto il: 21 mar 2012, 15:47


Torna a Matematica generale

Chi c’è in linea

Visitano il forum: Nessuno e 22 ospiti