• No results found

Oppsummering 6/4 Retningsfelt

N/A
N/A
Protected

Academic year: 2022

Share "Oppsummering 6/4 Retningsfelt"

Copied!
15
0
0

Laster.... (Se fulltekst nå)

Fulltekst

(1)

12.04.2016

1

Oppsummering 6/4

Retningsfelt

En måte å presentere differensialligninger Små rette linjer (piler) fordelt utover planet angir hvor raskt løsningen endrer seg.

(2)

12.04.2016

2

Numerisk løsning av 𝑦

= 𝐹(𝑥, 𝑦) og startverdi 𝑦 𝑥

0

= 𝑦

0

ved Eulers metode

Lineariserer i startpunktet:

𝑦 = 𝑦0+ 𝑦 𝑥0 𝑥 − 𝑥0

Dette gir 𝑦 𝑥1 ≈ 𝑦0+ 𝐹(𝑥0, 𝑦0)(𝑥1− 𝑥0)dersom 𝑥1− 𝑥0er liten.

Visualisering av Eulers metode i retningsfelt:

Generelt, når vi har beregnet 𝑦1, 𝑦2, … , 𝑦𝑛kan vi bestemme 𝑦𝑛+1: 𝑦𝑛+1= 𝑦𝑛+ ℎ𝐹(𝑥𝑛, 𝑦𝑛)

Separable differensialligninger

En differensialligning på formen 𝑞 𝑦 𝑦 = 𝑝(𝑥)

kalles separabel. Den kan skrives om til 𝑞 𝑦 𝑑𝑦 = 𝑝 𝑥 𝑑𝑥

slik at løsningen bestemmes ved antiderivasjon.

(3)

12.04.2016

3

Oppgave

Oppgave: Løs

𝑒𝑥𝑦𝑦 = 1 + 𝑦2.

Løsning: Differensialligningen kan skrives på formen

𝑦

1+𝑦2𝑦 = 𝑒−𝑥. Dvs den er separabel slik at

𝑦

1+𝑦2𝑑𝑦 = 𝑒−𝑥𝑑𝑥.

Oppgave

Deloppgave: Regn ut

𝑦

1+𝑦2𝑑𝑦og 𝑒−𝑥𝑑𝑥

Løsning: Substitusjon med 𝑢 = 1 + 𝑦2gir 𝑦

1 + 𝑦2𝑑𝑦 = 𝑦 𝑢

𝑑𝑢 2𝑦=1

2 1 𝑢𝑑𝑢 =1

2ln 𝑢 + 𝐶1=1

2ln 1 + 𝑦2 + 𝐶1 Substitusjon med 𝑢 = −𝑥gir

𝑒−𝑥𝑑𝑥 = 𝑒𝑢𝑑𝑢

−1= −𝑒𝑢+ 𝐶2= −𝑒−𝑥+ 𝐶2. Dermed, med 𝐶3= 𝐶2− 𝐶1,er

1

2ln 1 + 𝑦2 = − 𝑒−𝑥+ 𝐶3 slik at

1 + 𝑦2= 𝑒−2𝑒−𝑥+2𝐶3 = 𝐶𝑒−2𝑒−𝑥, 𝐶 = 𝑒2𝐶3, 𝑦 𝑥 = ± 𝐶𝑒−2𝑒−𝑥 − 1

(4)

12.04.2016

4

Førsteordens lineær differensialligning

Er på formen

𝑦 + 𝑓 𝑥 𝑦 = 𝑔(𝑥)

Løses ved først å multiplisere med integrerende faktor ℎ(𝑥)

ℎ 𝑥 𝑦′ + ℎ(𝑥)𝑓 𝑥 𝑦 = ℎ(𝑥)𝑔(𝑥) og så observere at venstre side er den deriverte av produktet ℎ 𝑥 𝑦 dersom ℎ = ℎ𝑓(𝑥)(separabel differensialligning med løsning ℎ 𝑥 = 𝑒 𝑓 𝑥 𝑑𝑥).

Førsteordens lineær differensialligning

Altså:

ℎ 𝑥 𝑦 = ℎ 𝑥 𝑔 𝑥 , Som løses ved antiderivasjon:

ℎ 𝑥 𝑦 = ℎ 𝑥 𝑔 𝑥 𝑑𝑥 slik at

𝑦 = 1

ℎ(𝑥) ℎ 𝑥 𝑔 𝑥 𝑑𝑥

(5)

12.04.2016

5

Tips til oblig: bruk av Matlab til å beregne integralet

𝑎𝑏

𝑓 𝑥 𝑑𝑥

𝑎 𝑏

𝑓 𝑥 𝑑𝑥

Simpson.m:

Tips til oblig: bruk av Matlab til å løse

differensialligningen 𝑦 = 𝐹 𝑥, 𝑦 , 𝑦 𝑎 = 𝑐 på (𝑎, 𝑏)

𝑦

= 𝐹 𝑥, 𝑦 , 𝑦 𝑎 = 𝑐

Eulerm.m:

(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)

Referanser

RELATERTE DOKUMENTER

Linjene for Trondheim og Nord-Jæren i figur 6 er så å si identiske med tilsvarende linjer i figur 4, Mens linjene for omegnskommunene er om lag 25% lavere enn tilsvarende linjer

Det ble ikke testet på spor mellom 20-50 mm, men det er ikke noe som tyder på at dette ville gitt resultater som endrer inntrykket av at løsningen med slapselameller påmontert

Direkte regulering av hedgefond er vanskelig fordi det ikke finnes en entydig definisjon av hva et hedgefond er, aktivitetene til hedgefond endrer seg raskt, og mange hedgefond

Kan vi finne en formel for hvor mange felter vi maksimalt kan dele planet i ved hjelp av n

La oss (ikke helt realistisk) anta at luft strømmer laminært rundt metallkula, og at luftmotstanden (friksjonskraften) kan skrives p˚ a formen f = − bv, der v er kulas hastighet, og b

Lys gr˚ a angir en flate med konstant negativ verdi av orbitalen, mørk gr˚ a angir en flate med tilsvarende konstant positiv verdi.. (Til høyre: z-aksen inn i planet, som i figuren p˚

Disse faktorene forventes å påvirke andelen kvinner i nasjonalforsamlinger, men ikke i domstolene direkte, så det utledes ingen hypotese om sammenhengen mellom disse og

Den må ikke være en film av en fungerende prototype i en reell brukskontekst: den kan presentere en modell som demonstrerer at tekniske løsningen fungerer og hvordan, og der