User Tools

Site Tools


04_teorie_chyb:0413_metoda_nejmensich_ctvercu

Teorie chyb

« 12. Vyrovnání měření obecně
» 14. Vyrovnání přímých měření

13. Metoda nejmenších čtverců

Úvod

Podmínka této metody obecně zní:

Ω=[v^v^]=[pvv]=vTPv=min.,\Omega {\rm =}\left[\hat{v}\hat{v}\right]=[pvv]={{\mathbf v}}^T\cdot {\mathbf P}\cdot {\mathbf v}=min.,

kde v^=v/m\hat{v}=v/m je tzv. normovaná oprava a pp tzv. váha. S těmito pojmy se nyní seznámíme podrobněji.

Váhy měření

Jsou to poměrná čísla, která kvalitativně hodnotí dosažený výsledek měření. Chceme, aby přesnější měření se ve vyrovnané veličině uplatnila více. Proto je zřejmá závislost na střední chybě. Váhy definujeme:

pi=Kmi2,p_i=\frac{K}{{\overline{m}}^2_i},

kde KK je vhodně volená konstanta. Z definice vyplývají následující vztahy:

p1: p2:p3::pn=1m12 :1m22::1mn2,pimi2=K.p_1:\ p_2:p_3:\cdots :p_n=\frac{1}{{\overline{m}}^2_1\ }:\frac{1}{{\overline{m}}^2_2}:\cdots :\frac{1}{{\overline{m}}^2_n} , p_i\cdot {\overline{m}}^2_i=K.

Slovně: součin váhy a čtverce střední chyby je v řadě měření konstantní. Když do tohoto vztahu teoreticky zavedeme takové měření, pro které p0=1p_0=1, říkáme, že tomuto měření odpovídá tzv. jednotková střední chyba m0{\overline{m}}_0. Ze vztahu přímo vyplývá význam konstanty: K=m02K={\overline{m}}^2_0. Obecný vztah pro volbu vah je pak:

pi=m02mi2.p_i=\frac{{\overline{m}}^2_0}{{\overline{m}}^2_i}.

Maticové vyjádření soustavy vah u nekorelovaných veličin:

P=(p1000p2000pn).{\mathbf P}=\left( \begin{array}{cccc} p_1 & 0 & \cdots & 0 \\ 0 & p_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & p_n \end{array} \right).

  • Poznámka: Váhy pip_i odvodíme zcela výjimečně z empirických středních chyb mim_i získaných analýzou daného souboru výsledků a to při neznalosti spolehlivých hodnot středních chyb mi{\overline{m}}_i nebo při neznalosti zákona hromadění chyb v uvažovaném případě. Praxe tedy určuje váhy dvojím způsobem. Při použití různých přístrojů, nebo metod, zvolíme váhy podle středních chyb určených pokud možno spolehlivě z dostatečně velikých výběrů. Jako konstantu KK můžeme zvolit čtverec základní střední chyby metody, nebo libovolnou číselnou hodnotu pro ulehčení výpočtů. Při stejné metodě měření určíme váhy podle zákona sčítání středních chyb. Napřed se zvolí vhodná váha jednotková tak, aby ostatní váhy nebyla čísla příliš velká nebo malá. Opravu vv budeme definovat jako rozdíl hodnoty vyrovnané a naměřené, přičemž hodnota vyrovnaná nemusí být přímo neznámá, ale může být funkcí několika neznámých. Tedy obecně:

vi=li(x1, x2, , xn)li,v_i={\overline{l}}_i\left(x_1,\ x_2,\ \cdots ,\ x_n\right)-l_i,

maticově

v=l(xT)l.{\mathbf v}=\overline{{\mathbf l}}\left({{\mathbf x}}^T\right)-{\mathbf l}.

Tento vztah používáme vždy linearizovaný.

Splnění podmínky minima funkce Ω\Omega se provádí postupem známým z matematiky:

Ω=vTPv=min.       Ωx=0\Omega ={{\mathbf v}}^T{\mathbf P \mathbf v}=min.\ {{\ \ \ \Rightarrow \ \ \ \frac{\partial \Omega }{\partial {\mathbf x}}={\mathbf 0}}} a 2Ω2x>0\frac{{\partial }^2\Omega }{{\partial }^2{\mathbf x}}>{\mathbf 0}.

Formulace vyrovnávacích úloh

Provede se řada měření (třeba i opakovaných) s různě přesnými výsledky l1l_1, …,lnl_n. Každé měření (případně skupinu opakovaných měření) pokládáme za náhodný výběr ze základního souboru možných hodnot. Dosažené výsledky představují neúplný soubor informací o neznámých skutečných hodnotách LL.

Naměřené hodnoty jsou zatíženy různými skutečnými chybami, takže platí vztah:

l + ε = L.{\mathbf l}{\ +}{\ \varepsilon }{\ =}{\ L}.

Z měření však εi{\varepsilon }_i a tudíž ani LiL_i nelze určit. Hledáme proto aproximaci LiL_i. Tu nazveme vyrovnanou hodnotou, označíme li{\overline{l}}_i a požadujeme, aby pro hodnoty lili=vi{\overline{l}}_i-l_i=v_i platila podmínka Ω=vTPv=min.\Omega ={{\mathbf v}}^T{\mathbf P \mathbf v}=min. Z této podmínky již vyplynou (podle druhu úlohy) výpočetní a kontrolní vzorce.


« 12. Vyrovnání měření obecně
» 14. Vyrovnání přímých měření

04_teorie_chyb/0413_metoda_nejmensich_ctvercu.txt · Last modified: 2024/09/17 08:15 by 127.0.0.1