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í.
. některých
programech proto prvky jakobiánu nepočítají každém iteračním kroku, ale
179
.
Za odhad chyby l)-vé iterace zde považuje
i 1,2,.jsou-li znaménka e(ít) c+shodná, krok prodlouží, aniž však přesáhne
plnou délku.
Ve vstupních údajích značí počet řešených nelineárních rovnic, <0) počáteční
odhad vektoru řešení vektor vztažných absolutních hodnot jednotlivých
složek řešení, využívaných při odhadu chyby iterací. Při výpočtu stejnosměrných
nebo časových charakteristik výhodné dosazovat maximální absolutní
hodnoty jednotlivých složek, dosažené během předchozích výpočtů. [f'(x (lĚ))] f(x lk)) j
minimalizovala, popřípadě aby alespoň nepřerostla normu ||f(x(ic))||. redukci
normy využívá některá metod pro numerickou jednorozměrnou minimalizaci..
Z hlediska výpočetních operací při implementaci Newtonovy-Raphsonovy
metody každém iteračním kroku nejnáročnější výpočet prvků Jacobiho matice
a řešení příslušné soustavy linearizovaných algebraických rovnic. Dosadíme-li vek
toru velká čísla, půjde odhad absolutní chyby. příklad úplného vývojového diagramu algoritmu Newtonovy-
-Raphsonovy metody dříve uvedeným způsobem řízení délky iteračních kroků.1, pokud a{k)jy 1
Doporučená hodnota pro volbu konstantního činitele tlumení iteračních kroků y
je 0,75. Tomu bylo možné předejít vhodnější volbou počátečního bodu x(0l
Na obr.
Některé modifikace Newtonovy-Raphsonovy metody používají podstatně
důmyslnější řízení délky iteračních kroků; např., (4. Broydenově metodě při každé
iteraci činitel a(A;+1) určující délku kroku volí tak, aby norma
llf(x(k+1))|| (t) a(-k+1).
Tímto postupem lze zabránit divergenci, konvergenci však zaručit nelze. 92b Ve
druhém těchto případů však iterace místo kořeni zkonvergují lokálního
minima.29)
Pokud zvolí <'°, výraz udává odhad relativní chyby.
S využitím uvedeného postupu pro řízení délky iteračních kroků lze dosáhnout
konvergence případech oscilujících iterací naznačených obr.
V některých programech specializovaným zaměřením vychází omezen
délky iteračních kroků fyzikálních vlastností prvků analyzované soustavy; např.
v programech pro analýzu elektronických obvodů polovodičovými diodami bi-
polárními tranzistory horní hranici délky iteračních kroků volí nejvyšší pří
pustné napětí přechodu PN. Volí proto
x(k)ly, je-li a(k)/y 1
1
