Numerisk løysing av likningar ved iterasjon.
Numerisk løysing av likningar
Mange likningar kan ikkje løysast eksakt med algebraiske metodar. Til dømes:
-
-
-
Likninga bruker cosinus, som høyrer til R2. Døme 2 og oppgåve 3 a) er difor fordjuping. Vi bruker òg nokre deriverte som du lærer å finne i kapittel 4: og .
Newtons metode gir oss ein måte å finne tilnærma løysingar med høg presisjon.
Ideen: Start med ei gjetting, og forbetre ho iterativt ved å følgje tangentlinja.
1. Vel ein startverdi nær løysinga
2. Iterer med formelen:
Geometrisk tolking:
- Teikn tangenten til i punktet
- Finn kvar tangenten kryssar -aksen
- Dette gir
Prosessen blir gjenteken til ønskt nøyaktigheit er nådd.
Utleiing av formelen
Tangentlinja til i har likning:
Vi finn kvar tangenten kryssar -aksen ():
Dette gir oss .
Bruk Newtons metode til å finne ved å løyse .
-
-
- Iterasjonsformel:
Startverdi:
Iterasjonar:
| 0 | 1 |
| 1 | |
| 2 | |
| 3 |
Svar: (etter 3 iterasjonar)
Finn ei løysing av med Newtons metode.
Derivert:
Iterasjonsformel:
Startverdi: (frå grafisk inspeksjon)
Iterasjonar:
| 0 | 0,5 | 0,378 |
| 1 | 0,7552 | |
| 2 | 0,7391 |
Svar: er løysinga av .
- Startverdien er nær løysinga
- nær løysinga
- er kontinuerleg og deriverbar nær løysinga
Når kan det gå gale?
- (divisjon med null)
- Dårleg startverdi (kan gje divergens eller feil løysing)
- Fleire løysingar (kan lande på "feil" løysing)
Konvergensfart:
Når metoden fungerer, doblast vanlegvis talet på korrekte siffer for kvar iterasjon (kvadratisk konvergens).
Bruk Newtons metode med 3 iterasjonar.
Finn ved å løyse med .
Finn ved å løyse med .
Løys med .
Likninga har éi reell løysing.
Vis at løysinga ligg mellom 2 og 3.
Bruk Newtons metode med til å finne løysinga med 4 desimalars nøyaktigheit.
Gjør dette i GeoGebra
Utforsk Newtons metode grafisk
Skriv inn én linje om gangen i inntastingsfeltet:
f(x) = x^3 - 2x - 5 x0 = 3 A = (x0, f(x0)) tang1 = Tangent(A, f) B = Intersect(tang1, xAxis) x1 = x(B) C = (x1, f(x1)) tang2 = Tangent(C, f) D = Intersect(tang2, xAxis)
Med språket satt til norsk forstår GeoGebra kommandonavnene både på norsk og engelsk.
Prøv å løyse (kubikkrot) med Newtons metode. Bruk .
-
-
Iterasjonsformel:
Iterasjonar:
| 0 | |
| 1 | |
| 2 | |
| 3 | |
| 4 |
Iterata divergerer — verdiane veks og sprett mellom positiv og negativ.
Årsak:
Nullpunktet er , men der er . Den deriverte er ubegrensa nær nullpunktet, og tangentlinjene peikar lenger og lenger bort. Newtons metode konvergerer ikkje.
Svar: Newtons metode divergerer her, og vi finn ikkje nullpunktet frå .
I praksis bruker vi eit stoppkriterium for å avgjere når vi har god nok tilnærming:
1. Endring mellom iterasjonar: (t.d. )
2. Funksjonsverdi nær null:
Vanlegvis er 5–10 iterasjonar nok til å oppnå meir enn 10 desimalars nøyaktigheit, takka vere den kvadratiske konvergensen.
Bruk Newtons metode på ulike problem.
(Fordjuping, R2-stoff) Finn nullpunktet til med . Utfør 3 iterasjonar.
Likninga har ei løysing nær . Bruk Newtons metode til å finne ho.
Forklar kvifor er ein dårleg startverdi for å løyse med Newtons metode.
Skriv ein algoritme (pseudokode eller Python) for Newtons metode.
Skriv ein funksjon newton(f, f_derivert, x0, toleranse, maks_iter) som returnerer eit tilnærma nullpunkt. Funksjonen skal iterere til eller maks tal på iterasjonar er nådd.
Test algoritmen på for å finne . Kor mange iterasjonar trengst for 10 desimalars nøyaktigheit?
Halveringsmetoden
Newtons metode krev den deriverte og ein god startverdi. Halveringsmetoden krev berre at er kontinuerleg og skifter forteikn i eit intervall.
Er og av motsett forteikn, er , og har eit nullpunkt i . Vi reknar ut midtpunktet og beheld den halvdelen der forteiknsskiftet er:
- Er , ligg nullpunktet i , og vi set .
- Elles ligg nullpunktet i , og vi set .
Kvar runde halverer intervallet. Startar vi med lengda , er lengda etter rundar, og nullpunktet ligg alltid i det intervallet vi har att.
La .
a) Vis at har eit nullpunkt i .
b) Utfør tre rundar med halveringsmetoden for hand.
c) Programmet nedanfor bruker halveringsmetoden. Forklar kva det gjer, og kva det skriv ut.
from math import exp
def f(x):
return exp(x) - 4
a = 1
b = 2
while b - a > 0.001:
m = (a + b) / 2
if f(a) * f(m) < 0:
b = m
else:
a = m
print(a, b)a) og .
er kontinuerleg og skifter forteikn, så har eit nullpunkt i .
Svar: , så det finst eit nullpunkt i .
b)
| Runde | |||
|---|---|---|---|
| 1 | 1 | 2 | 1,5 |
| 2 | 1 | 1,5 | 1,25 |
| 3 | 1,25 | 1,5 | 1,375 |
| Runde | Nytt intervall | ||
| 1 | |||
| 2 | |||
| 3 |
Svar: Etter tre rundar ligg nullpunktet i .
c)
f er funksjonen, a og b er endepunkta. Løkka køyrer så lenge intervallet er lengre enn 0,001. I kvar runde reknast midtpunktet m ut, og if-testen beheld den halvdelen der skifter forteikn.Lengda startar på 1 og halverast kvar runde. Etter 10 rundar er ho , og løkka stoppar.
Det eksakte nullpunktet er .
Svar: Programmet skriv ut endepunkta i det siste intervallet: 1.3857421875 1.38671875. Nullpunktet ligg mellom dei.
La .
Vis at har eit nullpunkt i .
Utfør tre rundar med halveringsmetoden. I kva intervall ligg nullpunktet etterpå?
Programmet nedanfor bruker halveringsmetoden.
def f(x):
return x**2 - 3
a = 1
b = 2
while b - a > 0.2:
m = (a + b) / 2
print(m)
if f(a) * f(m) < 0:
b = m
else:
a = mKva skriv programmet ut?
Kva tal nærmar midtpunkta seg viss vi let løkka gå lenger? Grunn svaret.
Kor mange rundar køyrer løkka viss vi byter 0.2 med 0.001?
Vi vil løyse likninga med programmet frå eksempel 4.
Endre programmet slik at det finn løysinga som ligg i . Kva skriv det ut?
Kva skriv programmet ut viss vi startar med a = -1 og b = 0? Er det ei løysing av likninga?
Legg til ein kontroll slik at programmet skriv ut ei melding når det ikkje er forteiknsskifte i startintervallet.
| Newtons metode | |
|---|---|
| Fart | Kvadratisk konvergens (rask) |
| Krav | , god startverdi |
| Pålitelegheit | Kan divergere |
| Halveringsmetoden | |
| Fart | Lineær konvergens (langsam) |
| Krav | Berre eit forteiknsskifte: |
| Pålitelegheit | Fungerer alltid når vi har forteiknsskifte |
I praksis: Halveringsmetoden gir eit grovt intervall rundt nullpunktet, og Newtons metode med startverdi i det intervallet gir raskt mange riktige desimalar.
Oppsummering
Newtons metode:
Framgangsmåte:
1. Omform likninga til
2. Finn
3. Vel ein god startverdi
4. Iterer til ønskt nøyaktigheit
Halveringsmetoden: Halver intervallet og behald halvdelen der skifter forteikn. Langsam, men verkar alltid når er kontinuerleg og .
Fordelar med Newtons metode:
- Rask konvergens (doblar presisjon per steg)
- Fungerer for dei fleste funksjonar
Ulemper med Newtons metode:
- Krev god startverdi
- Krev at
Dette kapitlet er skrevet av Anthropics toppmodeller (Claude Opus og Claude Fable) og er foreløpig ikke manuelt gjennomgått — kvalitetskontrollen gjøres av uavhengige KI-agenter, og innmeldte feil rettes fortløpende. Funnet en feil? Meld fra, så retter vi den. Les mer om hvordan innholdet lages.
