Metoda nejmenších čtverců: Porovnání verzí

Smazaný obsah Přidaný obsah
JAnDbot (diskuse | příspěvky)
m Odstraňuji šablonu {{link GA}} (vkládanou Wikidaty - skript od Amira)
Oprava interpunkce a úprava některých formulací v odstavci s odvozením metody.
Řádek 22:
\bold b \in\mathbb{R}^{n}</math>
 
takový, že <math>\bold b\not\in\mathcal{R}(\bold A)</math> (jinými slovy, neexistuje žádný vektor <math>\bold x</math> takový, aby byla nastala rovnost); symbol <math>\mathcal{R}(\bold A)</math> značí [[obor hodnot]] matice <math>\bold A</math>, tedy [[lineární obal]] jejichjejích sloupců. Metoda nejmenších čtverců hledá vektor <math>\bold x_{LS}</math> (LS je zkratkou anglického ''least squares'') minimalizujícísplňující
 
: <math>\|\bold A\bold x_{LS}-\bold b\|_2 = \min_{\bold x\in\mathbb{R}^m} \|\bold A\bold x-\bold b\|_2,</math>
 
nebo ekvivalentně, opravu pravé strany <math>\bold e_{LS}</math> minimalizujícísplňující
 
: <math>\|\bold e_{LS}\|_2 = \min_{\bold e\in\mathbb{R}^n} \|\bold e\|_2, \quad \text{kde} \quad (\bold b+\bold e)\in\mathcal{R}(\bold A).</math>
 
Zřejmě <math>\bold A \bold x_{LS} = \bold b + \bold e_{LS}</math> a oprava pravé strany <math>\bold e = \bold A \bold x - \bold b</math> je residuum odpovídající vektoru <math>\bold x</math>. Označíme-li <math>\bold e = [e_1,\ldots,e_n]^T</math> jednotlivé komponenty vektoru <math>\bold e</math>, pak z definice [[Eukleidovská norma|Eukleidovské normy]] přímo plyne, že minimalizujeme
 
: <math>\|\bold A\bold x - \bold b\|_2=\|\bold e\|_2=(\sum_{j=1}^n e_j^2)^{1/2} \qquad \Longleftrightarrow \qquad
\|\bold A\bold x - \bold b\|_2^2 = \|\bold e\|_2^2=\sum_{j=1}^n e_j^2.</math>
 
Takže minimalizujeme součet čtverců jednotlivých komponent vektoru <math>\bold e</math>, tj. odchylek jednotlivých komponent pravé strany <math>\bold b</math>. Nalezení vztahu pro vektor <math>\bold x_{LS}</math> splňující podmínku minimality, můžeme provést buď pomocí derivací (hledáním lokálního extrému) nebo přímo pomocí Pythagorovy věty.
 
'''Minimalizace kvadratického funkcionálu:''' Residuum <math>\bold e = \bold A \bold x - \bold b</math> je vektor, jehož normu minimalizovatminimalizujeme. Zřejmě
 
:<math>\arg\min \|\bold e\|_2 = \arg\min \|\bold e\|_2^2,\quad \text{kde} \quad \|\bold e\|_2^2 = \bold e^T \bold e = \sum e_j^2,</math>
Řádek 48:
 
:<math>\bold e \perp \mathcal{R}(\bold A) \qquad \Longleftrightarrow \qquad
\bold A^T\bold e = 0.</math>
 
Dosazením za residuum dostáváme
Řádek 62:
:<math>\bold x_{LS} = \left( \bold A^T \bold A \right)^{-1} \bold A^T \bold b,</math>
 
zZ předpokladu lineární nezávislosti sloupců přímo plyne, že řešení je jednoznačné.
 
Při praktickém výpočtu se ovšem normální soustava nikdy nesestavuje z důvodů numerické stability výpočtu. Použití zde odvozených vztahů pro praktický netriviální výpočet je jen obtížně ospravedlnitelné.<ref name="abjork">Åke Björck, ''Numerical Methods for Least Squares Problems'', SIAM Publications, Philadelphia PA, 1996</ref>