Till KTH:s startsida Till KTH:s startsida

Laboration 1 (Matlab)

I den här labben ska ni lösa en godtycklig ekvation y=f(x) med newton-raphsons metod.
Newton raphsons metod erbjuder en möjlighet att lösa ekvationer som saknar analytisk lösning, den är en iterativ numerisk ekvationslösare. Metoden består i att:

1 - ange ett närmevärde (x0) på en lösning
2 - beräkna funktionsvärdet f(x0)
3 - beräkna derivatan f'(x0)
4 - bilda korrektionstermen f(x0)/f'(x0)
5 - x1 = x0 - korrektionstermen
6 - upprepa från steg 2 med det nya x-värdet, gör så tills relativfelet < 1e-3


Kommentarer till ovanstående

- Lägg funktionen i en funktionsfil
- Rita upp funktionen och titta ungefär var en lösning finns, detta blir ert x0.
- Derivatan beräknar ni analytiskt eller enligt derivatans definition med ett litet h-värde (säg 1e-6)
- upprepningen av steg 2-6 görs med en slinga med avbrottsvilkor enligt relativfelet
- relativfelet kan uttryckas som abs((x1-x0)/x1), abs är absolutbeloppet
- Detalj: Skrivsättet 1e-6 är ett sätt att skriva 10^(-6) i Matlab.