Newtons metod

Newtons metod, eller Newton–Raphsons metod (efter Isaac Newton och Joseph Raphson) är en numerisk metod för att approximera nollställen till en funktion. Man använder alltså en numerisk metod för att hitta en rot till en ekvation, vilken går ut på att man väljer en punkt på kurvan som man räknar ut tangenten för. Det x-värde vid vilket tangenten skär x-axeln används sedan för att räkna ut en ny tangent i en iterativ process till dess att önskad noggrannhet uppnåtts.

Tangenten till en funktion ${\displaystyle f(x)}$ i punkten ${\displaystyle x_{0}}$ har enligt enpunktsformeln ekvationen

${\displaystyle y=f'(x_{0})(x-x_{0})+f(x_{0}).}$

Den skär x-axeln då y = 0, dvs:

${\displaystyle {\begin{aligned}0&=f'(x_{0})(x-x_{0})+f(x_{0})\\x-x_{0}&=-{\frac {f(x_{0})}{f'(x_{0})}}\\x&=x_{0}-{\frac {f(x_{0})}{f'(x_{0})}}.\end{aligned}}}$

Iterationsformeln blir alltså

${\displaystyle x_{n+1}=x_{n}-{\frac {f(x_{n})}{f'(x_{n})}}.}$

Idén är att steg för steg beräkna bättre och bättre approximationer till en rot till en ekvation ${\displaystyle f(x)=0}$. Vi börjar med en approximation ${\displaystyle x_{0}}$. Tangenten tillhörande funktionen f i punkten ${\displaystyle (x_{0},f(x_{0}))}$ skär x-axeln i en punkt (förutsatt att ${\displaystyle f'(x)\neq 0}$) som betecknas ${\displaystyle x_{1}}$. Man bestämmer denna punkt genom formeln

${\displaystyle x_{1}=x_{0}-{\frac {f(x_{0})}{f'(x_{0})}},}$

där ${\displaystyle f'(x_{0})}$ är värdet av derivatan till f i ${\displaystyle x_{0}}$, och så itererar man förloppet med ${\displaystyle x_{1}\;}$ som startpunkt, och så vidare. Den allmänna formeln blir då

${\displaystyle x_{n+1}=x_{n}-{\frac {f(x_{n})}{f'(x_{n})}}}$

Talföljden ${\displaystyle x_{0},x_{1},x_{2},\dots }$ konvergerar mot en rot r förutsatt att ${\displaystyle x_{0}}$ är tillräckligt nära den rot som ska approximeras.

Blir derivatans värden svårberäknade kan man approximera dem med formeln:

${\displaystyle f'(x_{n})={\frac {f(x_{n}+h)-f(x_{n})}{h}}}$

${\displaystyle \lim _{h\,\to \ 0}{\frac {f(x_{n}+h)-f(x_{n})}{h}},}$

eller mer exakt approximation med

${\displaystyle f'(x_{n})={\frac {f(x_{n}+h)-f(x_{n}-h)}{2h}}}$

${\displaystyle \lim _{h\,\to \ 0}{\frac {f(x_{n}+h)-f(x_{n}-h)}{2h}}.}$

Kvadratroten ur ett tal

Hur hittar man kvadratroten ur ett tal?. Det finns åtskilliga metoder för att hitta rötter och Newton Raphsons metod är en.

T.ex. om man önskar hitta kvadratroten ur 1395, så är det ekvivalent med att:

${\displaystyle \,x^{2}=1395}$

Funktionen i Newton Raphsons metod blir då,

${\displaystyle \,f(x)=x^{2}-1395}$

med derivatan,

${\displaystyle f'(x)=2x.\,}$

Med en inledande gissning 12, så blir ordningsföljden enligt Newton Raphsons metod:

${\displaystyle {\begin{matrix}x_{1}&=&x_{0}-{\dfrac {f(x_{0})}{f'(x_{0})}}&=&12-{\dfrac {12^{2}-1395}{2\cdot 12}}&=&64.125\quad \quad \quad {}\\x_{2}&=&x_{1}-{\dfrac {f(x_{1})}{f'(x_{1})}}&=&64.125-{\dfrac {64.125^{2}-1395}{2\cdot 64.125}}&=&42.93969298\\x_{3}&=&\vdots &=&\vdots &=&{\underline {37}}.71355835\\x_{4}&=&\vdots &=&\vdots &=&{\underline {37.3}}5145405\\x_{5}&=&\vdots &=&\vdots &=&{\underline {37.349698}}84\\x_{6}&=&\vdots &=&\vdots &=&{\underline {37.34969879}}\end{matrix}}}$

Dvs

${\displaystyle {\sqrt {1395}}={\underline {37,34969879}}}$

Där de korrekta siffrorna är understrukna. Vi ser här att bara med några få iterationer så får vi fram en lösning som stämmer överens på många decimaler.

Väljer vi att ${\displaystyle x_{0}\;\to \;0}$ i samma funktion som i exemplet ovan

${\displaystyle x_{1}=x_{0}-{\frac {f(x_{0})}{f^{\prime }(x_{0})}}\Leftrightarrow x_{1}=0-{\dfrac {0^{2}-1395}{2\cdot 0}}\;\to \;\infty }$ Som vi ser så går funktion mot oändligheten då x går mot noll.

Väljer vi att ${\displaystyle x_{0}\;\to \ \infty }$

${\displaystyle x_{1}=x_{0}-{\frac {f(x_{0})}{f^{\prime }(x_{0})}}\Leftrightarrow x_{1}=x_{0}-{\dfrac {x_{0}^{2}-1395}{2\cdot x_{0}}}\Leftrightarrow x_{1}={\dfrac {2x_{0}^{2}-x_{0}^{2}+1395}{2x_{0}}}\Leftrightarrow x_{1}={\dfrac {\not x_{0}\cdot (x_{0}+{\dfrac {1395}{x_{0}}})}{\not x_{0}\cdot 2}}\Leftrightarrow x_{1}={\dfrac {(x_{0}+{\dfrac {1395}{x_{0}}})}{2}}\;\to \;\infty \;}$

Här ser vi att funktionen ${\displaystyle \to \;\infty }$ för både noll och x. Detta implicerar då att Newton-Raphsons metod endast fungerar om den inledande gissningen är större än noll dvs ${\displaystyle x\;>0}$ och mindre än oändligheten dvs ${\displaystyle x<\infty }$ Så vårt intervall ${\displaystyle I=[a\;b]}$ borde då vara ${\displaystyle I=]0\;\infty [}$

Redan på Babylons tid visste man hur man kunde approximera rötter och rotberäkningar förekom också i de tidiga kulturerna i Egypten, Kina, Indien (Aryabhata) och Grekland (Heron). I Indien finns metoder för att beräkna närmevärden till kvadratrötter beskrivna på 500-talet f.Kr., i bland annat Baudhayanasutran.

Men det var Newton och Raphson som använde sig av analys för att generalisera denna urgamla metod för att hitta rötterna till en godtycklig ekvation. Newtons metod publicerades först år 1685 i boken A Treatise of Algebra both Historical and Practical av John Wallis. 5 år senare publicerade Joseph Raphson en förenklad version i avhandlingen Analysis aequationum universalis. Där Raphson visade att det är algebraisk metod som är begränsad till polynom. Newtons metod har också beskrivits av Isaac Newton år 1669 i sin bok De analysi per aequationes numero terminorum infinitas (som publicerades 1711 av William Jones) och i De metodis fluxionum et serierum infinitarum(författad 1671, översatt och publicerad som Method of Fluxions år 1736 av John Colson). Men det ska nämnas att den beskrivning som Newton gav i de ovannämnda böckerna skiljer sig mycket från den beskrivning som har angivits ovan.

• Forsling & Neymark, Matematisk Analys En variabel, 228-229, ISBN 91-47-05188-4
• Författarna och Liber AB, Tabeller och formler för NV- och TE- Programmen femte upplagan, S 56, ISBN 978-91-47-01746-1
• Engelska wikipedias beskrivning av Newton Raphsons metod http://en.wikipedia.org/wiki/Newton%27s_method

• Wikimedia Commons har media som rör Newtons metod.
  Bilder & media