Visste du att...

Naturvetenskap.org varje månad donerar 10 % av sina annonsintäkter till välgörenhet?

Eulers stegmetod

Skriven av Johan Asplund
Senast uppdaterad söndag, 24 juli 2011 13:07

Eulers stegmetod är en numerisk metod för att approximalt bestämma ett värde av en funktion om man får givet en differentialekvation som funktionen uppfyller, och ett startvärde. Detta gör man genom att utnyttja linjära approximationer,

\nor y^\pri(x_1) \ap y(x_0) + y^\pri(x_0) \cd h \ .

En mer konkret förklaring är att om man har en differentialekvation och vet att \nor y(0.5) = 1 så kan man approximalt bestämma t.ex. \nor y(3) med hjälp av Eulers stegmetod. Det hela går ut på att man först bestämmer hur långa steg i \nor x-led man gör. Om man väljer att ta steg med längden 0.5 så kan man inleda med att beräkna lutningen i punkten \nor (0.5;1), som var startpunkten (\nor A_0). Sedan går man 1 steg med längden 0.5 i \nor x-led i den riktning som lutningen går. Sedan beräknar man den nya lutningen i punkten man nu befinner sig (\nor A_1) och fortsätter så, tills man är vid \nor y(3).

Eulers stegmetod

Figuren illustrerar en röd lösningskurva till en differentialekvation,
och gröna linjesegment som approximerar värden på lösningskurvan med hjälp av Eulers stegmetod.

 

Exempel

Vi har givet att \nor xy^\pri - y = x^2, \ y(1) = 0, och vi vill beräkna \nor y(2). Vi kan inte lösa den typen av differentialekvationer med de metoder vi känner till, så vi får använda Eulers stegmetod istället. Vi sätter steglängden till \nor h = 0.5, eftersom approximationen blir bättre, ju kortare steg man tar. Vi vet att startpunkten är i \nor (1,0). Vi skriver om differentialekvationen som

\nor y^\pri(x) = \fr{x^2+y}{x} \ .

Detta ger lutningen

\nor y^\pri(1) = \fr{1^2+0}{1} = 1 \ .

Vi går i \nor x-led med steglängden 0.5 med lutningen 1. Då får vi \nor y(1.5).

\nor y(1.5) \ap y(1) + y^\pri(1) \cd h = 0 + 1 \cd 0.5 = 0.5 \ .

Nu har vi en ny punkt. \nor (1.5;0.5). Lutningen i denna punkten är

\nor y^\pri(1.5) = \fr{1.5^2+0.5}{1.5} = 1.83\ .

Vi tar ytterligare ett steg och \nor y(2) kan beräknas.

\nor y(2) \ap y(1.5) + y^\pri(1.5) \cd h = 0.5 + 1.83 \cd 0.5 = 1.415 \ .

Det visar sig att lösningen till differentialekvationen är

\nor y = x^2-x \ , 

vilket gör att det verkliga värdet för \nor y(2) är 2. Vi hade alltså ett fel på 0.6 ungefär. Ju kortare steg man tar, ju bättre approximation får man.

Eulers stegmetod

Röd kurva: Lösningen till differentialekvationen.

Gröna linjer: Vårat resultat med Eulers stegmetod. \nor h = 0.5.

Blåa linjer: Eulers stegmetod med \nor h = 0.25.