Kniha je úvodem do metod praktického modelování, analýzy, návrhu a optimalizace elektrotechnických zařízeni na číslicovém počítači. Výklad je doprovázen jednoduchými názornými příklady řešených úloh z různých odvětví elektrotechniky.Kniha je určena inženýrům a technikům, kteří se zabývají moderním návrhem elektrotechnických zařízení.
2.6. praxi proto nejčastěji používají integrační
metody řádem 6. Komplikovaný vzorec
(6. Ukázali jsme již,
že přímá zpětná Eulerova metoda představuje aproximaci řešení prvními dvěma
členy Taylorova rozvoje lichoběžníková metoda aproximuje řešení rozdílem
prvních dvou členů jeho Taylorových rozvojů dvou bodech. Metody vyšších řádů
založené Taylorově rozvoji však vyžadují znalost vyšších derivací funkce f(x, t),
které obvykle nejsou dispozici analytickém tvaru nutné proto počítat
numericky. Jednokrokové metody vyšších řádů
Místní zbytková chyba integračních metod řádu jdoucím nule zmenšuje
rychlostí 0(hr+1).42), tak (6. Jednou
z možností, které zde nabízejí, využití Taylorova rozvoje.
6.můžeme dosadit jak (6.44)
kde
k, f(*„,í)
/ h\
k J
( h\
k f(*„ k2’ 2)
k f(xn hk3, h)
Jelikož jde metodu čtvrtého řádu, dalo usuzovat, metoda dovoluje
volbu dlouhého integračního kroku. praxi nejoblíbenější metoda Rungeho
a Kutty čtvrtého řádu, která dána vzorcem
x„+i ~{ki 2k3 k4) (6. Navíc nesmíme zapomenout, řádem metody roste počet aritmetic
kých operací prováděných jednom integračním kroku, tím roste strojový čas
a zvětšují chyby zaokrouhlení. Bohužel tomu tak není. Proto metody vyšších řádů zpravidla dovolují při určité zadané
přesnosti použít delší integrační krok než metody nižších řádů. Jelikož numerické derivování obtížné nepřesné, integrační metody
vyšších řádů vycházející přímo Taylorova rozvoje praxi neujaly.44) nedovoluje odvodit jednoduchý vztah pro odhad místní zbytkové chyby
metody, proto krok obvykle nutné volit opatrnosti mnohem kratší, než by
310
.43) opět buď nebo hodnotu
vypočítanou pomocí některého prediktoru.
Z Taylorova rozvoje však byly odvozeny metody Rungeho Kutty, nichž
se potřeba vyšších derivací funkce f(x, obchází jejím vyhodnocováním bodech,
které leží uvnitř integračního kroku. ovšem platí
pouze předpokladu, délka kroku není omezena hlediska numerické ne
stability metody, neboť metody vyšších řádů obvykle bývají méně stabilní než metody
nižších řádů.
Integrační metody vyšších řádů lze konstruovat různými způsoby
