Metoda navadne iteracije: Razlika med redakcijama
Izbrisana vsebina Dodana vsebina
Ena od pomembnejših metod za numerično reševanje enačb |
m dp/negibna točka glede na Brouwerjev izrek o negibni točki/decimalne vejice |
||
Vrstica 1:
'''Metóda navádne iterácije''' je v [[matematika|matematiki]] preprosta in učinkovita [[numerična metoda]] za določanje
: <math> x=g(x) \,\! . </math>
== Potek postopka ==▼
▲==Potek postopka==
Enačbo oblike ''x'' = ''g''(''x'') rešimo po metodi navadne iteracije z naslednjim postopkom:
*
*
* Na enak način izračunamo še nadaljnje približke: <math>x_2=g(x_1),~ x_3=g(x_2),~\ldots,~x_n=g(x_{n-1})\,\!</math>.
* Postopek končamo, ko smo zadovoljni z natančnostjo dobljenega približka.
Tako dobljeno zaporedje približkov [[konvergenca|konvergira]] k fiksni točki, če je v neki okolici fiksne točke [[odvod]] funkcije ''g'' po absolutni vrednosti manjši od 1 in če začetni približek leži v tej okolici.
== Zgled ==
Rešimo enačbo <math>x= \cos x\,\!</math> (kjer je ''x'' [[kot]] v [[radian]]ih). Ker odvod funkcije po absolutni vrednosti nikoli ni večji od 1, je izbira začetnega približka precej poljubna. Izberimo npr. začetni približek 1. Po večkratnem izvajanju funkcije [[kosinus]] dobimo naslednje približke:
:<math>x_0=1\,\!</math>▼
: <math>
: <math>
: <math>
: <math>
: <math>
Vidimo, da se razlika med zaporednima približkoma manjša, kar je znak, da zaporedje konvergira. Po večjem številu korakov dobimo približek:
:<math>x_{20}=0.73918\,\!</math>▼
:<math>x_{21}=0.73901\,\!</math>▼
Iz tega sklepamo, da so prve tri decimalke pravilne rešitve enake <math>x=0.739\,\!</math>. Če želimo izračunati [[rešitev enačbe|rešitev]] na več decimalk natančno, postopek nadaljujemo.▼
==Drugi načini uporabe==▼
===Iskanje ničel funkcije===▼
▲Iz tega sklepamo, da so prve tri decimalke pravilne rešitve enake <math>x=0
▲== Drugi načini uporabe ==
▲=== Iskanje ničel funkcije ===
Metodo navadne iteracije lahko uporabimo za iskanje [[ničla funkcije|ničel]] poljubne [[funkcija|funkcije]] ''f'', če enačbo ''f''(''x'') = 0 najprej preoblikujemo v obliko ''x'' = ''g''(''x''). To je ponavadi možno na več načinov, vendar pa vsa preoblikovanja niso enako uspešna.
Vrstica 32 ⟶ 41:
<li>
Polinom lahko preoblikujemo po naslednejm postopku:
:<math>0=x^3-3x+5\,\!</math>▼
: <math>
: <math>
: <math> x=\frac{x^3+5}{3} </math>
V dobljeno iteracijsko funkcijo <math>g(x)=\frac{x^3+5}{3}</math> vstavljamo različne začetne približke in opazimo, da dobljeno zaporedje nikoli ne konvergira.
</li>
<li>
Isti polinom lahko preoblikujemo tudi drugače:
▲:<math>x^3-3x+5=0\,\!</math>
: <math> x^3
: <math> x
: <math> x=\sqrt[3]{3x-5} </math>
Tako dobljena iteracijska funkcija :<math>g(x)=\sqrt[3]{3x-5}</math> je uspešnejša. Pri različnih izbirah začetneg približka dobimo že po nekaj korakih ničlo na več decimalk natančno:
: <math> x=-2
</li>
</ul>
=== Reševanje enačbe oblike ''h''(''x'') = ''k''(''x'') ===
Splošno enačbo oblike ''h''(''x'') = ''k''(''x'') lahko preoblikujemo z uporabo [[inverzna funkcija|inverza]] funkcije ''h'' ali ''k'' (če obstajata). Tako dobimo dve obliki:
▲:<math>x=h^{-1}(k(x))\,\!</math>
: <math>x=
Praviloma vsaj ena od teh dveh metod (ob izbiri primernega začetnega približka) vodi k rešitvi.
== Glej tudi ==
* [[metoda pospešene iteracije]] (Newtonova ali tangentna metoda)
[[Kategorija:Numerična analiza]]
|