n p n = ( 1 3 )n p n = 1 3 p n−1 p n = 10 3 p n−1 − p n−20 0,10000 × 10 1 0,10000 × 10 1 0,10000 × 10 11 0,33333 × 10 0 0,33333 × 10 0 0,33333 × 10 02 0,11111 × 10 0 0,11111 × 10 0 0,11110 × 10 03 0,37037 × 10 −1 0,37036 × 10 −1 0,37000 × 10 −14 0,12346 × 10 −1 0,12345 × 10 −1 0,12230 × 10 −15 0,41152 × 10 −2 0,41150 × 10 −2 0,37660 × 10 −26 0,13717 × 10 −2 0,13717 × 10 −2 0,32300 × 10 −37 0,45725 × 10 −3 0,45723 × 10 −3 −0,26893 × 10 −28 0,15242 × 10 −3 0,15241 × 10 −3 −0,92872 × 10 −2Vidíme, že první rekurentní vztah je stabilní, kdežto druhý algoritmus je nestabilní.2.7 Kontrolní otázky a cvičeníCvičení 2.1. Určete nejmenší interval, v němž musí ležet výsledek, použijete-lipřesných hodnot místo zaokrouhlených. Předpokládá se, že všechna čísla v následujícíchvýpočtech jsou správně zaokrouhlena. (Při výpočtech zaokrouhlujte našest desetinných míst.)(i) 2,547 · 1,25,(ii)6,580,128 .Cvičení 2.2. Necht’ jsou čísla x ∗ = 0,013; y ∗ = 0,24 správně zaokrouhlena nauvedený počet míst. Spočtěte f(x,y) = xsin y pro uvedené údaje a určete početplatných míst.Cvičení 2.3.Určete nejmenší interval, v němž musí ležet výsledek operacef(x 1 ,x 2 ,x 3 ) = x 1 (x 2 + x 3 ),použijete-li přesných hodnot místo zaokrouhlených. Předpokládá se, že všechnačísla x ∗ 1 = 0,2; x∗ 2 = 0,26; x∗ 3 = 0,75 jsou správně zaokrouhlena.Cvičení 2.4. Předpokládejte, že máte n správně zaokrouhlených čísel x 1 ,... , x nna d i míst. Chcete spočítat součet těchto n čísel na d = min d i míst. Záleží na18
tom, zda všechna čísla napřed zaokrouhlíte a pak sečtete, nebo napřed sečtete apak zaokrouhlíte? Pokud ano tak proč?Cvičení 2.5. Posloupnost p n = ( 1 3 )n generujeme rekurzivně(i) p n = 5 6 p n−1 − 1 6 p n−2, p 0 = 1, p 1 = 1 3 a n ≥ 2,(ii) p n = 5 3 p n−1 − 4 9 p n−2, p 0 = 1, p 1 = 1 3 a n ≥ 2,přičemž zaokrouhlujeme čísla v normovaném tvaru na pět platných míst. Rozhodněte,zda se jedná o stabilní algoritmy.Cvičení 2.6. Pro libovolné přirozené číslo n, n ≤ 1, sestrojte algoritmus, kterýpro libovolné body x 0 ,... ,x n a x dá výstupP n+1 = (x − x 0 )(x − x 1 ) · · · (x − x n ).Cvičení 2.7. Dokažte vztahy pro odhad absolutní chyby aritmetických operací aodvod’te vztahy pro relativní chyby.2.8 VýsledkyCvičení 2.1.(i) 〈3,17039; 3,19711〉(ii) 〈51,166381; 51,646119〉Cvičení 2.2. f(x,y) = 0,309 × 10 −2 a počet platných míst je j = 1Cvičení 2.3. f(x 1 ,x 2 ,x 3 ) ∈ 〈0,1496; 0,2545〉Cvičení 2.4. Napřed sečíst a pak zaokrouhlit.Cvičení 2.5.(i) stabilní(ii) nestabilníCvičení 2.6. P := (x−x 0 ); i := k+1; While P ≠ 0 and i ≤ n Do P := P(x−x i )19
- Page 2 and 3: Matematický ústav Slezské univer
- Page 4 and 5: 6.3.2 Jacobiova a Gaussova-Seidelov
- Page 6 and 7: STRUČNÝ NÁHLED STUDIJNÍ OPORYTe
- Page 8 and 9: - juliánský kalendář byl zavede
- Page 10 and 11: Definice 2.1. Necht’ x je přesn
- Page 12 and 13: polovinu jednotky řádu poslední
- Page 14 and 15: 2.4 Chyba součtu, rozdílu, souči
- Page 16: 2.5 Dobře a špatně podmíněné
- Page 21 and 22: akterizovány hodnotami parametrů
- Page 23 and 24: 3.1 Základní tvary interpolační
- Page 25 and 26: nutné přepočítat všechny funda
- Page 27 and 28: x 0 y 0> f[x 0 ,x 1 ]x i y i f[x i
- Page 29 and 30: Důkaz. Označme g(x) = f(x) − P
- Page 31 and 32: 3.3 Interpolace na ekvidistantní s
- Page 33 and 34: Řešení. Po sestavení tabulky po
- Page 35 and 36: fundamentální polynom, který zaj
- Page 37 and 38: a jejich derivacel ′ 0 (x) = 1 6
- Page 39 and 40: Požadavky na spojitost prvních a
- Page 41 and 42: Z druhé derivace v krajním bodě
- Page 43 and 44: Řešení této soustavy jec 2 =
- Page 45 and 46: ⎧⎪⎨S(x) =⎪⎩x, −∞ < x
- Page 47 and 48: 4 METODA NEJMENŠÍCH ČTVERCŮ...
- Page 49 and 50: kde j = 0,1,... ,k, který lze jedn
- Page 51 and 52: Příklad 4.5. Nalezněte polynom p
- Page 53 and 54: Věta 4.7. Necht’ na vektorovém
- Page 55 and 56: 5 NUMERICKÉ ŘEŠENÍ NELINEÁRNÍ
- Page 57 and 58: stačí vzítl ≥ logε(1 − K)|x
- Page 59 and 60: mít shodné kořeny. Tento postup
- Page 61 and 62: x k x k+1 x k+2Obr. 5.2Předpoklád
- Page 63 and 64: (5.1). Vskutku((α = lim x k+1 = li
- Page 65 and 66: Tentokrát dělícím bodem není p
- Page 67 and 68:
Pro srovnání přidáváme hodnotu
- Page 69 and 70:
Protožex −∞ ∞ 0 − 3 22P 1
- Page 71 and 72:
Protožex −∞ ∞ 0P 1 (x) - + -
- Page 73 and 74:
Cvičení 5.8.(i) x 0 = 0, x 1 = 1,
- Page 75 and 76:
pravidlax =( D1D ,..., D ) ⊤n, (6
- Page 77 and 78:
pro n = k − 1. Pro n = k matici A
- Page 79 and 80:
(iv) Řešíme soustavu Ux = y,⎛1
- Page 81 and 82:
Příklad 6.7. Gaussovou eliminačn
- Page 83 and 84:
Věta 6.9. Necht’ ‖ · ‖ je v
- Page 85 and 86:
= H 2 x (k−2) + (H + E)h.= H k x
- Page 87 and 88:
Označíme-li ve vztahu (6.4)H =
- Page 89 and 90:
Výpočet uspořádáme do tabulky:
- Page 91 and 92:
Cvičení 6.4. Se zaokrouhlováním
- Page 93 and 94:
7 NUMERICKÉ INTEGROVÁNÍ... co by
- Page 95 and 96:
Věta 7.1. Existuje číslo η i
- Page 97 and 98:
K odvození chyby integrace využij
- Page 99 and 100:
= f ′′ (η)2= f ′′ (η)2∫
- Page 101 and 102:
Opět zavedeme substituci o označe
- Page 103 and 104:
x 2i f(x 2i ) x 2i−1 f(x 2i−1 )
- Page 105 and 106:
7.6 VýsledkyCvičení 7.1.(i) 0,30
- Page 107 and 108:
f(x,y) lipschitzovská vzhledem k y
- Page 109 and 110:
a podle (8.2) mámey i+1 = y i + hf
- Page 111 and 112:
y i+1y(x)==y ikk21xixi+h/2x i+1Pou
- Page 113 and 114:
i x i y i y(x i ) {k 1 ,k 2 ,k 2 ,k
- Page 115 and 116:
Tedy pro každé x ∈ 〈0, 0,2〉
- Page 117 and 118:
Cvičení 8.2.(i)(ii)i x i y i |y(x
- Page 119:
Literatura[1] Burden, R.L., Faires,