Newtons metode

Frå testwiki
Hopp til navigering Hopp til søk

Newtons metode i kalkulus, òg kjend som Newton-Raphson-metoden, er ein iterativ metode for å finna nullpunkta til ein gjeven funksjon f, altså løysinga av f(x)=0.

Han må ikkje blandast saman med newtons metode i optimering, som baserer seg på å finna stasjonære punkt for ein gjeven funksjon f.

Definisjon

Gjeve eit startpunkt xn som den iterative metoden startar frå og ein funksjon f ein ønskjer å finna nullpunkta til, så er newtons algoritme i kalkulus gjeven som

xn+1=xnf(xn)Δf(xn)

Her er Δf(xn) den deriverte av f(xn). Under særskilde føresetnadar bundne av valet av startpunkt xk, så vil følgja xn+1 konvergera mot løysinga f(xx+1)=0.

I det høvet der f er ein vektorvaluert funksjon, så vert dette:

xn+1=xnJf(xn)1f(xn)

Jf er her kjend som jakobimatrisa for funksjonen f.

Utleiing av definisjon

Ein ønskjer å finna nullpunkta for ein funksjon f, altså løysinga av f(x)=0 ved ein iterativ algoritme på forma:

xk+1=xk+h

som vil konvergera mot løysinga f(xx+1)=0 gjeve eit særskilt startpunkt xk og eit steg h. Me ønskjer i denne seksjonen å motivera steget h. For eit gjeve punkt xk så ønskjer me å tilnærma funksjonen f i punktet og analysera kva for ein h for eit punkt xk+h som fører til at ein går mot eit minimum,

altså ei løysing av f(x)=0. Newtons metode går ut på å gjera dette og å tilnærma funksjonen i punktet ved hjelp av taylorpolynom.

For ein fleirvariabels funksjon f så er taylorpolynomtilnærminga av funksjonen f gjeven som:

f(xk+h)=f(xk)+Δf(xk)Th+(1/2)hTΔ2f(xk)h+

Det vert då nytta ei førsteordens taylorpolynomtilnærming av funksjonen i punktet, som er ei rett linje og er tangent til punktet xk:

f(xk+h)f(xk)+Δf(xk)Th

Me finn den h-en sånn at denne tilnærminga (tangentfunksjonen) er 0. Løyser ein dette (f(xk+h)=0 med omsyn på h), så får ein newtonsteget:

h=f(xk)Δf(xk)T

Newtonsteget i Newtons metode i kalkulus er altså utleitt frå ei førsteordens taylorpolynomtilnærming for funksjonen. For denne tilnærminga løyser ein f(xk+h)=0 med omsyn på h. Ein får då punktet xk+h som minimerer tangentfunksjonen i punktet xk.

Ut i frå h, så får me då:

xn+1=xk+h=xnf(xn)Δf(xn)T