The wiki page is under active construction, expect bugs.

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
statnice:bakalar:b4b01num [2025/06/07 12:45] – [3. Metoda sečen] mistrjirkastatnice:bakalar:b4b01num [2025/06/07 13:38] (current) – [Numerická integrace] mistrjirka
Line 142: Line 142:
  
 ===== Numerické metody řešení nelineárních rovnic ===== ===== Numerické metody řešení nelineárních rovnic =====
-===== Numerické metody řešení nelineární rovnice =====+
  
 **Numerické metody řešení rovnice $f(x) = 0$** slouží k nalezení kořenů, kdy analytické řešení není možné. Předpokladem je separace kořenů – nalezení intervalu $\langle a, b \rangle$ s $f(a) \cdot f(b) < 0$ a jediným kořenem. Metody se liší rychlostí konvergence, požadavky na derivace a garancemi konvergence ,  **Numerické metody řešení rovnice $f(x) = 0$** slouží k nalezení kořenů, kdy analytické řešení není možné. Předpokladem je separace kořenů – nalezení intervalu $\langle a, b \rangle$ s $f(a) \cdot f(b) < 0$ a jediným kořenem. Metody se liší rychlostí konvergence, požadavky na derivace a garancemi konvergence , 
Line 218: Line 218:
 **Visualizace**: **Visualizace**:
  
-<tikzjax> +{{:statnice:bakalar:pasted:20250607-125539.png?490}}
-\usepackage{amsmath} +
-\usepackage{pgfplots} +
-\usetikzlibrary{automata, positioning, arrows, calc, cd, intersections,arrows.meta} +
-\begin{document} +
- +
-\begin{tikzpicture}[thick,yscale=0.8] +
-% Axes +
-  \draw[-latex,name path=xaxis] (-1,0) -- (12,0) node[above]{\large $x$}; +
-  \draw[-latex] (0,-2) -- (0,8) node[right]{\large $y$}; +
- +
-  % Function plot +
-  \draw[ultra thick, orange,name path=function] +
-    plot[smooth,domain=1:9.5] +
-      (\x, {0.1*\x^2 - 1.5}) node[left]{$F(x)$}; +
- +
-  % Tangent at x^(k)=8 +
-  \node[violet,right=0.2cm] at (8,4.9) {\large tangent}; +
-  \draw[gray,thin,dotted] (8,0) -- (8,4.9) +
-    node[circle,fill,inner sep=2pt]{}; +
-  \draw[dashed, violet,name path=Tfunction] +
-    plot[smooth,domain=4.25:9.5] +
-      (\x, {1.6*\x - 7.9}); +
- +
-  % Labels for x^(k) and x^(k+1) +
-  \draw (8,0.1) -- (8,-0.1) node[below] {$x^{(k)}$}; +
-  \draw [name intersections={of=Tfunction and xaxis}] +
-    ($(intersection-1)+(0,0.1)$) -- ++(0,-0.2) +
-    node[below,fill=white] {$x^{(k+1)}$}; +
-\end{tikzpicture} +
- +
-\end{document} +
-</tikzjax>+
  
  
Line 264: Line 232:
 **Visualizace**: **Visualizace**:
  
-<tikzjax> +{{:statnice:bakalar:pasted:20250607-125807.png?450}}
-\usepackage{amsmath} +
-\usepackage{pgfplots} +
-\usetikzlibrary{automata, positioning, arrows, calc, cd, intersections,arrows.meta} +
-\begin{document}+
  
-\begin{tikzpicture}[thick, xscale=1, yscale=1] 
-  % Axes 
-  \draw[->] (0,0) -- (5,0) node[right]{$x$}; 
-  \draw[->] (0,-1) -- (0,4) node[above]{$f(x)$}; 
- 
-  % Function 
-  \draw[blue,smooth,domain=0.5:4.5] plot (\x,{0.1*(\x-3)^3+0.5}); 
- 
-  % Initial guesses x0=0.5, x1=1.5 
-  \coordinate (X0) at (0.5,{0.1*(0.5-3)^3+0.5}); 
-  \coordinate (X1) at (1.5,{0.1*(1.5-3)^3+0.5}); 
-  \fill (X0) circle (1.5pt) (X1) circle (1.5pt); 
- 
-  % Secant line 
-  \draw[red,name path=sec] (X0) -- (X1); 
- 
-  % Intersection X2 with f 
-  \draw[name path=fun2] plot [smooth,domain=0.5:4.5] 
-    (\x,{0.1*(\x-3)^3+0.5}); 
-  \draw[name intersections={of=sec and fun2,by=X2}]; 
-  \draw[dashed] (X2) -- (X2|-0,0) node[below]{$x_2$}; 
-  \fill (X2) circle (1.5pt); 
-\end{tikzpicture} 
- 
-\end{document} 
-</tikzjax> 
  
  
Line 308: Line 246:
 **Konvergence**: Lineární při zachování intervalu s kořenem \\ **Konvergence**: Lineární při zachování intervalu s kořenem \\
 **Visualizace**: **Visualizace**:
- 
-<tikzjax> 
-\usepackage{amsmath} 
-\usepackage{pgfplots} 
-\usetikzlibrary{automata, positioning, arrows, calc, cd, intersections,arrows.meta} 
-\begin{document} 
- 
-\begin{tikzpicture}[scale=1.0] 
-\draw[->] (0,0) -- (5,0) node[right] {$x$}; 
-\draw[blue,thick] plot[domain=0.5:4.5] (\x,{0.3*(\x-1.5)^2 - 0.7}); 
-\draw[red,dashed] (1,0) node[below] {$a_0$} -- (1,{-0.7+0.3*(1-1.5)^2}); 
-\draw[red,dashed] (4,0) node[below] {$b_0$} -- (4,{-0.7+0.3*(4-1.5)^2}); 
-\draw[orange,thick] (1,{-0.7+0.3*(1-1.5)^2}) -- (4,{-0.7+0.3*(4-1.5)^2}); 
-\draw[red,dashed] (2.2,0) node[below] {$x_0$} -- (2.2,0); 
-\draw[green!60!black,dashed] (2.2,0) -- (2.2,{-0.7+0.3*(2.2-1.5)^2}); 
-\draw[green!60!black,thick] (2.2,{-0.7+0.3*(2.2-1.5)^2}) -- (4,{-0.7+0.3*(4-1.5)^2}); 
-\draw[red,dashed] (2.8,0) node[below] {$x_1$} -- (2.8,0); 
-\end{tikzpicture} 
- 
-\end{document} 
-</tikzjax> 
  
  
 +{{:statnice:bakalar:pasted:20250607-130029.png?450}}
  
 ==== 5. Metoda prosté iterace ==== ==== 5. Metoda prosté iterace ====
Line 341: Line 259:
 **Konvergence**: Lineární při $|\phi'(x)| < 1$ v okolí kořene. Řád $p$ pro $\phi^{(k)}(x) \neq 0$ ($k = 1, \dots, p-1$) a $\phi^{(p)}(x) \neq 0$ \\ **Konvergence**: Lineární při $|\phi'(x)| < 1$ v okolí kořene. Řád $p$ pro $\phi^{(k)}(x) \neq 0$ ($k = 1, \dots, p-1$) a $\phi^{(p)}(x) \neq 0$ \\
 **Visualizace** (Cobweb diagram): **Visualizace** (Cobweb diagram):
- 
-<tikzjax> 
-\usepackage{amsmath} 
-\usepackage{pgfplots} 
-\usetikzlibrary{automata, positioning, arrows, calc, cd, intersections,arrows.meta} 
-\begin{document} 
- 
-\begin{tikzpicture}[scale=1.0] 
-\draw[->] (0,0) -- (4,0) node[right] {$x$}; 
-\draw[->] (0,0) -- (0,4) node[above] {$y$}; 
-\draw[blue,thick] plot[domain=0:3.5] (\x,{0.5*\x + 1.5}) node[right] {$y=\phi(x)$}; 
-\draw[gray,dashed] (0,0) -- (3.5,3.5) node[right] {$y=x$}; 
-\draw[red,->] (1,0) node[below] {$x_0$} -- (1,1.5); 
-\draw[red,->] (1,1.5) -- (1.5,1.5); 
-\draw[red,->] (1.5,1.5) -- (1.5,1.75); 
-\draw[red,->] (1.5,1.75) -- (1.75,1.75); 
-\draw[red,->] (1.75,1.75) -- (1.75,1.875); 
-\draw[fill=magenta] (2,2) circle (2pt) node[above] {Pevný bod}; 
-\end{tikzpicture} 
- 
-\end{document} 
-</tikzjax> 
- 
  
  
 +{{:statnice:bakalar:pasted:20250607-130215.png?450}}
 ==== Problematika separace kořenů ==== ==== Problematika separace kořenů ====
  
Line 381: Line 277:
     * Nespojitosti, periodické funkce – vyžadují speciální algoritmy (např. kombinace s derivacemi)      * Nespojitosti, periodické funkce – vyžadují speciální algoritmy (např. kombinace s derivacemi) 
  
-===== Numerické řešení soustav lineárních rovnic =====+==== Numerické řešení soustav lineárních rovnic ====
  
-**Numerické řešení soustav lineárních rovnic, možné problémy, argumenty pro použití finitních nebo iteračních metod.**+**Metody pro řešení soustav lineárních rovnic $A\vec{x} = \vec{b}$, jejich matematické formulace, problémy a kritéria volby mezi přímými a iteračními postupy.**
  
-==== 1. Finitní (přímé) metody ====+=== Přímé (finitní) metody ===
  
-  * **Příklady**: Gaussova eliminační metoda, LU dekompozice, QR rozklad.\\+Poskytují teoreticky přesné řešení v konečném počtu kroků: 1. **Gaussova eliminace**:\\ 
 +- Převádí matici $A$ na **horní trojúhelníkový tvar** pomocí ekvivalentních úprav.\\ 
 +- **Algoritmus**:\\ 
 +- Pro $k = 1$ až $n-1$:\\ 
 +$a^{(k)}_{i,j} = a^{(k-1)}_{i,j} - \frac{a^{(k-1)}_{i,k}}{a^{(k-1)}_{k,k}} \cdot a^{(k-1)}_{k,j}$\\ 
 +- **Zpětná substituce**:\\ 
 +$x_i = \frac{1}{a^{(i-1)}_{i,i}} \left( a^{(i-1)}_{i,n+1} - \sum_{j=i+1}^n a^{(i-1)}_{i,j}x_j \right)$ \\ 
 +- Výběr hlavního prvku redukuje zaokrouhlovací chyby.
  
-  * **Výhody**: +  - **LU rozklad**: 
-    * Pro malé a střední matice (např. $n < 1000$) jsou efektivní.\\+    * Rozklad $A = LU$ na **dolní ($L$)** a **horní ($U$) trojúhelníkovou matici**.\\
  
-    * Přesné řešení (do chyby zaokrouhlení) v konečném počtu kroků.\\+    * Řeší se dvě soustavy: 
 +      * $L\vec{y} = \vec{b}$ (dopředná substituce),\\
  
-  * **Problémy**: +      * $U\vec{x} = \vec{y}$ (zpětná substituce).\\
-    * **Výpočetní náročnost**: Složitost $O(n^3)$ pro metody jako LU dekompozice.\\+
  
-    * **Numerická nestabilita**: Pro nevhodné matice (např. špatně podmíněné) mohou vzniknout velké chyby zaokrouhlení.\\+    * Efektivní pro opakované výpočty s různými $\vec{b}$ 
  
-    * **Paměťové nároky**: Ukládání celé matice (i pro řídké matice).\\+=== Iterační metody ===
  
-    * **Pivotace**: Někdy nutná pro stabilitu (např. při eliminaci bez pivotace mohou nastat dělení nulou). +Generují posloupnost aproximací $\vec{x}^{(k)} \to \vec{x}$: 1. **Jacobiho metoda (JIM)**:\\ 
-  * **Příklady z materiálů**: +- Rozklad $A = D + L + U$ ($D$ diagonální, $L$ ostře dolní, $U$ ostře horní trojúhelníková).\\ 
-    * **Normální rovnice** pro řešení přeurčených soustav (např. nejmenší čtverce): $A^T A x = A^T b$.\\+- **Iterační vzorec**:\\ 
 +$\vec{x}^{(k+1)} = D^{-1} \left( \vec{b} - (L + U) \vec{x}^{(k)} \right)$ \\ 
 +- Jednoduchá implementace, paralelizovatelná, ale pomalá konvergence 
  
-    * **QR rozklad** pro řešení soustav s LN sloupci: $Ax = b \Rightarrow Q R x = b \Rightarrow x = R^{-1} Q^T b$.+  - **Gaussova-Seidelova metoda (GSM)**: 
 +    * Využívá již aktualizované hodnoty v aktuální iteraci:\\ 
 +$\vec{x}^{(k+1)} = (D + L)^{-1} \left( \vec{b} - U \vec{x}^{(k)} \right)$ \\
  
-==== 2. Iterační metody ====+    * Rychlejší konvergence než JIM, ale sekvenční výpočet  
 +  - **Superrelaxační metoda (SOR)**: 
 +    * Zavádí relaxační parametr $\omega$ pro urychlení konvergence:\\ 
 +$\vec{x}^{(k+1)} = (D + \omega L)^{-1} \left[ (1-\omega)D \vec{x}^{(k)} - \omega U \vec{x}^{(k)} \right] + \omega (D + \omega L)^{-1} \vec{b}$ \\
  
-  * **Příklady**: Jacobiho metoda, Gaussova-Seidelova metoda, metoda konjugovaných gradientů (pro symetrické pozitivně definitní matice).\\+    * Pro $\omega = 1$ přechází na GSM. Optimální $\omega$ zrychluje konvergenci, ale špatná volba způsobí divergenci 
  
-  * **Výhody**: +=== Problémy a kritéria volby metod ===
-    * Efektivní pro **velké řídké matice** (např. $n > 10^4$), protože využívají jen operace s nenulovými prvky.\\ +
- +
-    * Nižší paměťové nároky (neukládá se celá matice, jen její struktura).\\+
  
   * **Problémy**:   * **Problémy**:
-    * **Konvergence** není zaručena: Závisí na vlastnostech matice (např. diagonální dominance, pozitivní definitnost).\\+    * **Špatná podmíněnost**: Malé změny v $A$ nebo $\vec{b}$ vedou k velkým změnám řešení (kritérium: velké $\|A^{-1}\|$) \\
  
-    * **Pomalá konvergence** bez předpodmínky (preconditioner).\\+    * **Zaokrouhlovací chyby**: Akumulují se v přímých metodách, zejména bez výběru hlavního prvku \\
  
-    * **Citlivost na volbu počátečního odhadu**. +    * **Konvergence iterací**: Zajištěna pouze pokud $\rho(B) < 1$ (spektrální poloměr iterační matice) \\
-  * **Příklady z materiálů**: +
-    * **Nejmenší čtverce** pro přeurčené soustavy: Minimalizace $||Ax - b||_2$ pomocí QR nebo SVD.\\+
  
-    * **Pseudoinverze** pro matice bez plné hodnosti (např. SVD rozklad).+    * **Řídké matice**: Přímé metody ztrácejí na efektivitě kvůli “zaplnění” (fill-in), iterační metody ji zachovávají  
 +  * **Volba metody**:
  
-==== 3. Porovnání a výběr metody ====+^**Parametr**             ^**Přímé metody**               ^**Iterační metody**                       ^ 
 +|**Velikost matice**      |Malé až střední systémy        |Velké systémy ($n > 10^3$)                | 
 +|**Struktura matice**     |Plné nebo řídké                |Řídké (např. diskretizace PDE)            | 
 +|**Přesnost**             |Vysoká (teoreticky přesné)     |Přibližná (kontrola rezidua $\|\vec{r}\|$)| 
 +|**Vícenásobná $\vec{b}$**|Efektivní (jednorázový rozklad)|Výhodné s dobrým počátečním odhadem       | 
 +|**Paměťová náročnost**   |Vyšší (např. LU: $O(n^2)$)     |Nižší (např. JIM: $O(n)$)                 |
  
-^**Kritérium**        ^**Finitní metody**                       ^**Iterační metody**                    ^ +=== Odhady chyb a optimalizace ===
-|**Rozměr matice**    |Malé a střední ($n < 1000$)              |Velké ($n > 10^4$)                     | +
-|**Dostupnost paměti**|Vyžadují větší paměť                     |Efektivní pro řídké matice             | +
-|**Přesnost**         |Vysoká (do chyby zaokrouhlení)           |Závislá na konvergenci a iteracích     | +
-|**Použití**          |Řešení přesných soustav, nejmenší čtverce|Řídké matice, paralelizovatelné systémy|+
  
-==== 4. Typické problémy v numerice ====+  * **Reziduální odhad**: $\vec{r} = \vec{b} - A\vec{x}_c$ umožňuje zpřesnění řešení \\
  
-  * **Špatně podmíněné matice**: Malá nebo nulová determinant, velké chyby v řešení.\\+  * **Složitost**: 
 +    * Přímé metody: $O(n^3)$ pro GEM, $O(n^2)$ pro zpětný chod.\\
  
-  * **Numerická nestabilita**: Např. v Gaussově eliminaci bez pivotace.\\+    * Iterační metody: $O(n^2)$ na iteraci (pro husté matice), $O(n)$ pro řídké 
  
-  * **Chyby zaokrouhlení**: Akumulace v přímých metodách pro velké $n$.\\ 
  
-  * **Řídké matice**: Iterační metody výrazně převyšují přímé metody.+===== Numerická integrace =====
  
-==== 5. Ukázkové příklady z materiálů ====+**Numerická integrace slouží k přibližnému výpočtu určitého integrálu, když analytické řešení není možné. Zahrnuje diskretizaci intervalu, výpočet ploch pod křivkou pomocí jednoduchých geometrických tvarů a řízení chyb.**
  
-  * **Přeurčená soustava** $Ax = b$: Řešeno metodou nejmenších čtverců pomocí normálních rovnic nebo QR.\\+==== Metody numerické integrace ====
  
-  * **Příklad řídké matice**: Iterační metoda (např. CG) pro řešení $Ax = b$, kde $A$ má jen $O(n)$ nenulových prvků.+  - **Metoda levých obdélníků** 
 +    * **Myšlenka:** Funkce se aproximuje hodnotou v **levém konci** každého podintervalu.\\
  
-==== 6. Závěr ====+    * **Matematicky:** Pro interval $[a, b]$ rozdělený na $n$ podintervalů o šířce $h = \frac{b-a}{n}$:\\ 
 +$$ 
 +I_L = h \sum_{i=0}^{n-1} f(x_i), \quad x_i = a + i \cdot h. 
 +$$\\
  
-  * **Finitní metody** jsou vhodné pro malé matice a přesné řešení.\\+    * **Řád metody:** 1 (přesná pro polynomy stupně 0, např. konstantní funkce). 
 +  - **Metoda středních obdélníků** 
 +    * **Myšlenka:** Funkce se aproximuje hodnotou ve **středu** každého podintervalu.\\
  
-  * **Iterační metody** dominují u velkých řídkých matic a při náročnosti na paměť.\\ +    * **Matematicky:**\\
- +
-  * Výběr závisí na rozměru matice, její struktuře a požadované přesnosti. +
-===== Numerická integrace ===== +
- +
-**Numerická integrace - princip, možné problémy, volba metody, problematika odhadu chyb.** +
- +
-==== Princip numerické integrace ==== +
- +
-Numerická integrace je postup přibližného výpočtu **definitního integrálu**\\+
 $$ $$
-\int_{a}^{b} f(x) \, dx+I_S = h \sum_{i=0}^{n-1} f\left(x_i + \frac{h}{2}\right).
 $$\\ $$\\
-pro případy, kdy analytická řešení (primitivní funkce) nelze získat nebo je výpočet složitý. Základní přístup spočívá v **discretizaci** intervalu $[a, b]$ a aproximaci integrálu pomocí **polynomů**, **geometrických útvarů** nebo **kvadraturních vzorců**. 
  
-==== Hlavní metody numerické integrace ====+    * **Řád metody:** 2 (přesná pro lineární polynomy, např. $x^1$). 
 +  - **Lichoběžníková metoda (Trapezoid)** 
 +    * **Myšlenka:** Funkce se aproximuje **lineárně** mezi konci podintervalu (plocha lichoběžníku).\\
  
-=== 1. Obdélníková metoda (Rectangle Rule) === +    * **Matematicky:**\\
- +
-  * **Princip**: Rozdělí interval $[a, b]$ na $n$ podintervalů, v každém z nich integrand nahradí **konstantou** (hodnotou v levém, pravém nebo středním bodě).\\ +
- +
-  * **Formule**:\\+
 $$ $$
-\int_{a}^{b} f(x) \, dx \approx \frac{b - a}{n} \sum_{i=1}^{n} f(x_i)+I_T = \frac{h}{2} \left[ f(a) + 2 \sum_{i=1}^{n-1} f(x_i) + f(b) \right].
 $$\\ $$\\
-kde $x_i$ je bod v $i$-tém podintervalu (např. střed = **střední metoda**).\\ 
- 
-  * **Řád metody**: pro levé a pravé obdelníky: 1 (chyba je $O(h)$, kde $h = \frac{b-a}{n}$)., pro střední je řád 2 (chyba je $O(h^2)$) 
- 
-  * **Výhody**: Jednoduchost, malá početní náročnost.\\ 
- 
-  * **Nevýhody**: Malá přesnost, nevhodné pro funkcí s rychlým zrychlením. 
- 
-=== 2. Lichoběžníková metoda (Trapezoidal Rule) === 
  
-  * **Princip**: Interval rozdělí na $n$ podintervalů, které aproximuje **lichoběžníky** (lineární interpolace).\\+    * **Řád metody:** 2 (přesná pro lineární polynomy). 
 +  - **Simpsonova metoda** 
 +    * **Myšlenka:** Funkce se aproximuje **kvadratickou** parabolou přes dvojice podintervalů.\\
  
-  * **Formule**:\\+    * **Matematicky (pro sudé $n$):**\\
 $$ $$
-\int_{a}^{b} f(x) \, dx \approx \frac{h}{2} \left[ f(a) + 2\sum_{i=1}^{n-1} f(x_i) + f(b) \right], \quad h = \frac{b - a}{n}+I_S = \frac{h}{3} \left[ f(a) + 4 \sum_{i=1,3,\dots}^{n-1} f(x_i) + 2 \sum_{i=2,4,\dots}^{n-2} f(x_i) + f(b) \right].
 $$\\ $$\\
  
-  * **Řád metody**: 2 (chyba $O(h^2)$).\\+    * **Řád metody:** 4 (přesná pro kubické polynomy, např. $x^3$). 
 +  - **Gaussova kvadratura** 
 +    * **Myšlenka:** Volí **optimální body a váhy** v intervalu $[-1, 1]$ pro maximální přesnost. Transformuje se na $[a, b]$.\\
  
-  * **Výhody**: Lépe přesná než obdélníková, stabilní pro hladké funkce.\\ +    * **Matematicky:** Pro $n$ bodů:\\
- +
-  * **Nevýhody**: Chyba se zvětšuje pro funkce s velkými druhými derivacemi. +
- +
-=== 3. Simpsonova metoda (Simpson’s Rule) === +
- +
-  * **Princip**: Užívá **parabolické interpolace** (kvadratické polynomy) na párových podintervalech.\\ +
- +
-  * **Formule**:\\+
 $$ $$
-\int_{a}^{b} f(x) \, dx \approx \frac{h}{3} \left[ f(a) + 4\sum_{i=1}^{n/2} f(x_{2i-1}) + 2\sum_{i=1}^{(n/2)-1} f(x_{2i}) + f(b) \right]+I_G = \frac{b-a}{2} \sum_{i=1}^n w_i f\left( \frac{b-a}{2} t_i + \frac{a+b}{2} \right),
 $$\\ $$\\
-(vyžaduje sudé $n$).\\+kde $t_i$ jsou kořeny Legendreových polynomů, $w_i$ odpovídající váhy.\\
  
-  * **Řád metody**: 4 (chyba $O(h^4)$).\\ +    * **Řád metody:** $2n-1$ (pro $n$ bodů přesná pro polynomy stupně $\leq 2n-1$).
- +
-  * **Výhody**: Vysoká přesnost pro hladké funkce.\\ +
- +
-  * **Nevýhody**: Závislost na počtu podintervalů, složitější výpočet -> velmi pomalé.+
  
 === 4. Newton-Cotesovy vzorce === === 4. Newton-Cotesovy vzorce ===
Line 530: Line 414:
  
 ==== Problémy a optimalizace ==== ==== Problémy a optimalizace ====
 +  * **Problémy:** Numerické chyby (zaokrouhlovací, metody), singularity, oscilace funkce.\\
  
-=== 1. Odhad chyby a kompozitní metody ===+=== Metoda polovičního kroku a odhad chyby ===
  
-  * **Kompozice**: Rozdělení $[a, b]$ na $n$ podintervalů a aplikace metody na každý.\\ +  * **Princip:** Integrál $I_h$ se spočte s krokem $h$, pak s $h/2$ pro jemnější dělení. Chyba pro metodu řádu $p$ je:\\
- +
-  * **Chyba trapezoidální metody**:\\+
 $$ $$
-E \approx -\frac{(b - a)^3}{12n^2} f''(\xi), \quad \xi \in [a,b]+E_h \approx \frac{I_h - I_{h/2}}{2^p - 1}.
 $$\\ $$\\
  
-  * **Chyba Simpsonovy metody**:\\ +  * **Vylepšení integrálu:**\\
-$$ +
-E \approx -\frac{(b - a)^5}{180n^4} f^{(4)}(\xi)+
 $$ $$
 +I_{\text{lepšené}} = I_{h/2} + \frac{I_{h/2} - I_h}{2^p - 1}.
 +$$\\
  
-=== 2. Richardsonova extrapolace ===+  * **Souvislost s Richardsonovou extrapolací:** Tento postup je jejím speciálním případem. Kombinací výsledků pro různé $h$ eliminuje vedoucí člen chyby 
  
-Sloučení dvou odhadů s různými kroky $h$ a $h/2$ pro odstranění vedení v chybě.\\ +=== Řád metody ===
-- **Příklad**: U trapezoidální metody:\\ +
-$$ +
-  I(h/2) = I(h) + \frac{E(h)}{4} \Rightarrow \text{vyšší řád} O(h^4) +
-  $$+
  
-=== 3. Adaptivní kvadratura ===+  * **Význam:** Udává, do jakého stupně polynomu metoda počítá integrál přesně. Např. metoda řádu 2 přesně integruje polynomy stupně $\leq 1$ (tj. $x^{1}$ a nižší).\\ 
 + 
 +  * **Chyba metody:** Pro krok $h$ a řád $p$ je globální chyba $O(h^p)$. 
 + 
 + 
 +==== Řád chyby ==== 
 + 
 +  * **Definice:** Řád chyby $p$ znamená, že chyba metody klesá jako $h^p$ při zmenšování kroku $h$. 
 +  * **Příklad:** Pro metodu s řádem 2, zmenšení kroku $h$ na polovinu sníží chybu 4krát, protože $E \propto h^2$. 
 + 
 +=== Adaptivní kvadratura ===
  
 Metody jako **Gauss-Kronrod** automaticky rozdělují intervaly podle lokálního chybového odhadu (lepší efektivita pro funkcí s různými charakteristikami). Metody jako **Gauss-Kronrod** automaticky rozdělují intervaly podle lokálního chybového odhadu (lepší efektivita pro funkcí s různými charakteristikami).
Line 579: Line 468:
 ==== Kritéria volby metody ==== ==== Kritéria volby metody ====
  
-  - **Hladkost funkce**: Pro funkce s velkými derivacemi preferovat metody s vysokým řádem (např. Simpsonova).\\ +  - **Hladkost funkce**: Pro funkce s velkými derivacemi preferovat metody s vysokým řádem (např. Simpsonova). 
- +  - **Interval délky**: Pro velké intervaly nebo singulární body použít adaptivní metody nebo substituci.
-  - **Interval délky**: Pro velké intervaly nebo singulární body použít adaptivní metody nebo substituci.\\ +
   - **Výpočetní náročnost**: Obdélníková je nejrychlejší, ale malá přesnost.   - **Výpočetní náročnost**: Obdélníková je nejrychlejší, ale malá přesnost.
  

Trace:

Navigation

Playground

QR Code
QR Code statnice:bakalar:b4b01num (generated for current page)