Gruppuppgifter 1 MMA132, Numeriska metoder, distans

Gruppuppgifter 1
MMA132, Numeriska metoder, distans
Uppgifter märkta med ? redovisas
1. Läs om felkalkyl i enkla fall sidan 1.2-1.3. Givet
a = 1,23, Ea = 0,005
c = 0,00438 ± 0,5 · 10−5
b = 23,71, Eb = 0,003
d = 20,104 ± 0,2
Hur många korrekta decimaler och signifikanta siffror har närmevärdena? Hur många korrekta
decimaler och signifikanta siffror har närmevärdena?
2. Antag att alla givna siffor är korrekta för nedanstående värden:
a = 1,435171
c = 4,7556 · 10−5
b = 23,6212
d = 13,1
Skriv a, b, c och d med fyra korrekta siffror.
3. Låt x = 1,2 ± 0,3 och y = 8,2 ± 0,1 Beräkna z med absoluta och relativa felgränser i följande
fall:
a) z = x + y
b) z = 0,7x
c) z = 3x + 7y
4. Låt x = 7,2 ± 0,3 och y = 4,7 ± 0,1 Beräkna z med absoluta och relativa felgränser i följande
fall:
a) z = 0,4xy
b) z =
c) z =
?5.
6.
y
2x
x3 y 2
1000
(a) Låt x = 1,414 och y = 0,09125 vara värden givna med 4 korrekta siffror.
Ange z1 = x + y och z2 = x · y med rätt antal korrekta siffror.
(b) Låt x = 31,415 och y = 0,027182 vara värden givna med 5 korrekta siffror.
Ange z1 = x + y och z2 = x · y med rätt antal korrekta siffror.
√
a) Låt y = x + sin(x).Bestäm ett närmevärde på y om x ≈ 2,222 (3 korrekta decimaler).
Ge en bra felgräns för y och ange antalet korrekta decimaler och siffror i y.
b) Låt y = cos(t). Bestäm ett närmevärde på y om t = 0,02 ± 0,01. Ge en bra felgräns för y.
Tips: Använd 3 termer i Taylor-utvecklingen.
7. Låt y = 5e−x sin(t). Bestäm ett närmevärde på y om x ≈ 0,34 (2 korrekta decimaler) och
t ≈ 1,3 (1 korrekt decimal). Ge en bra felgräns för y och ange antalet korrekta decimaler och
siffror i y.
?8.
Låt f (x) = x3 − 15,1x2 + 10. Fyll i tredje kolumnen nedanstående tabell (det räcker med 4
korrekta siffror)
x
9,4
9,6
10,4
10,6
x̃
9
10
10
11
ef (x) = f (x̃) − f (x)
Ef ≈ |f 0 (x̃)|Ex
Notera att |x̃ − x| < 0,5 = Ex .
Uppskatta med hjälp av fortplantningsformeln Ef ≈ |f 0 (x̃)|Ex och fyll i fjärde kolumnen i
ovanstående tabell.
Ser du några problem med denna uppskattning?
Hur skulle du göra en bättre uppskattning?
9. Läs om kancellation och utskiftning sidan 1.5.
√
√
a) Beräkna 800,0001 − 800 på räknedosan. Beräkna också värdet av det ekvivalenta ut0,0001 √
trycket √800,0001+
. Varför blir resultaten olika? Vilket resultat är mest noggrant?
800
b) Låt x = 0,001 (exakt). Beräkna f (x) = 2 sin(x) − sin(2x) på räknedosan. Beräkna också
värdet av det ekvivalenta uttrycket f (x) = 4 sin(x) sin2 (0,5x). Varför blir resultaten olika?
Vilket resultat är mest noggrant?
10. Låt f (x) = π − | cos(x)| − x. Bestäm nollstället f (x) = 0 i intervallet [0, 3,2] med två signifikanta siffror med intervallhalveringsmetoden.
11. Antag att ekvationen 2 − x − ex = 0 har en lösning nära x = 0,45. Gör två iterationer med
Newton-Raphsons metod för att hitta lösningen och ange en feluppskattning för svaret.
12. Antag att ekvationen x = 5(1 − e−x ) har en lösning nära x = 5. Gör två iterationer med
sekantmetoden för att hitta lösningen och ange en feluppskattning för svaret.
?13.
Bestäm en lösning till ekvationen e 2 = e2 + x med minst 3 korrekta decimaler med NewtonRaphsons metod. Tips: Nollstället ligger mellan 0 och 2.
?14.
Bestäm en lösning till ekvationen 7e−y = 5 med minst 3 korrekta decimaler med sekantmetoden. Tips: Funktionen har två nollställen. Ett mellan 0 och 1 och ett mellan 0 och -1.
x
x
2
15. Beräkna hur många iterationer med intervallhalveringsformeln som skulle behövas för att hitta
ett nollställe för funktionen f (x) = 3x2 + sin(x − 11) − ecos(ln(x)) i intervallet [0,2] med 4
signifikanta decimaler.
Tips: Om din miniräknare inte har log2 så gäller log2 (a) = log(a)
log(2) där log är logaritmen avseende
på valfri bas.
16.
a) Gör minst tre iterationer för att beräkna roten till ekvationen x = 0,2 cos(x) med fixpunktsmetoden och gör feluppskattning (metodoberoende).
b) Beräkna en rot till ekvationen x = 3 cos(x) med minst 2 korrekta decimaler med fixpunktsmetoden (ekvationen måste först skrivas om för att få konvergens).
17. Läs om Newton-Raphsons modifierade metod sidan 2.7. Lös x = 1 − 0,8 sin(x) med denna
metod. Använd f 0 (a) ≈ f 0 (x0 ).
Lösningar och svar
1.
a
b
c
d
korrekta
decimaler
2
2
5
0
signifikanta
siffror
3
4
3
2
2. Om vi skriver ut de olika värdena med 4 korrekta siffror får vi
b = 2,362 · 101
a = 1,435
c = 4,756 · 10−5
det är inte möjligt att skriva d med fyra korrekta siffor då vi endast har tre korrekta siffror ifrån
början.
3. Låt x = 1,2 ± 0,3 och y = 8,2 ± 0,1.
a) Med z = x + y får vi approximationen z̃ = x̃ + ỹ = 1,2 + 8,2 = 9,4 och den absoluta
felgränsen
Ez = Ex + Ey = 0,3 + 0,1 = 0,4.
Den relativa felgränsen blir
Rz =
Ez
0,4
=
= 0,043.
z̃
9,4
b) Med z = 0,7x får vi approximationen z̃ = 0,7x̃ = 0,7 · 1,2 = 0,84 och den absoluta
felgränsen
Ez = 0,7Ex = 0,7 · 0,3 = 0,21.
Den relativa felgränsen blir
Rz =
Ez
0,21
=
= 0,25.
z̃
0,84
c) Med z = 3x + 7y får vi approximationen z̃ = 3x̃ + 7ỹ = 3 · 1,2 + 7 · 8,2 = 61 och den
absoluta felgränsen
Ez = 3Ex + 7Ey = 3 · 0,3 + 7 · 0,1 = 1,6.
Den relativa felgränsen blir
Ez
1,6
=
= 0,026.
z̃
61
Observera att här har vi avrundat uppåt för säkerhets skull.
Rz =
4. Låt x = 7,2 ± 0,3 och y = 4,7 ± 0,1.
a) Med z = 0,4xy får vi approximationen z̃ = 0,4x̃ỹ = 0,4 · 7,2 · 4,7 = 13,536 och den
relativa felgränsen
Rz = Rx + Ry =
Ex Ey
0,3 0,1
+
=
+
= 0,063.
x̃
ỹ
7,2 4,7
Observera att här har vi avrundat uppåt för säkerhets skull. Den absoluta felgränsen blir
Ez = Rz z̃ = 0,063 · 13,536 = 0,86.
Observera att här har vi avrundat uppåt för säkerhets skull.
b) Med z =
y
får vi approximationen z̃ =
2x
Rz = Rx + Ry =
ỹ
2x̃
=
4,7
2·7,2
= 0,3264 och den relativa felgränsen
Ex Ey
0,3 0,1
+
=
+
= 0,063.
x̃
ỹ
7,2 4,7
Observera att här har vi avrundat uppåt för säkerhets skull. Den absoluta felgränsen blir
Ez = Rz z̃ = 0,063 · 0,3264 = 0,021.
Observera att här har vi avrundat uppåt för säkerhets skull.
x3 y 3
7,23 · 4,72
c) Med z =
får vi approximationen z̃ =
= 8,245 och den relativa fel1000
1000
gränsen
0,3
0,1
Ex Ey
+
=3·
+2·
= 0,168.
Rz = 3Rx + 2Ry =
x̃
ỹ
7,2
4,7
Den absoluta felgränsen blir
Ez = Rz z̃ = 0,168 · 8,245 = 1,39.
Observera att här har vi avrundat uppåt för säkerhets skull.
6.
a) Enligt felfortplantningsformeln: Ey ≈ |y 0 |Ex ≈ 0,2707 · 0,0005 ≈ 0,00014 ≈ 0,0002,
y = 2,2860 ± 0,0002. 3 korrekta decimaler, 4 korrekta siffror.
00
b) Enligt felfortplantningsformeln med två termer: Ey ≈ |y 0 |Et + |y2 | Et2 ≈ 0,02 · 0,01 + 1 ·
0,012
2 ≈ 0,00025 ≈ 0,0003, y = 0,9998 ± 0,0003. 3 korrekta decimaler, 3 korrekta siffror.
7. Löses enklast genom att använda den flerdimensionella varianten av felfortplantningsformeln
(sidan 1.4 i Egnesund)
n
X
Ef ≈
Ei |fxi (x)|
i=1
på y = f (x,t) = 5e−x sin(t) med x = (x̃,t̃) = (0,34, 1,3)
Ey ≈ |fx (x̃,t̃)|Ex +|ft (x̃,t̃)|Et = |−5e−0,34 sin(1,3)|·0,005+|5e−0,34 cos(1,3)|·0,05 ≈ 0,0647
Alternativ lösning: y = 5e−x sin(t) = a(x)b(t)
Kan då konstatera att Ry = Ra + Rb ≈
Ea
|ã|
+
Eb
|b̃|
Med Ea ≈ |a0 (x̃)|Ex och Eb ≈ |b0 (t̃)|Et fås
Ry = Ra + Rb ≈
| − 5e−0,34 |0,005 | cos(1,3)|0,05
+
= 0,005 + | tan 1,3| · 0,05 ≈ 0,1851
|5e−0,34 |
| sin(1,3)|
9. Resultaten kan bli olika på olika räknedosor. MATLAB ger:
√
√
a) 800,0001− 800 ≈ 28,28427301522880−28,28427124746190 ≈ 0,000001767766896421108
(ca 14 decimaler, dvs 9 korrekta siffror).
0,0001 √
√
≈ 0,000001767766897723655 (ca 16 korrekta siffror, mest noggrant).
800,0001+ 800
b) f (0,001) = 2 sin(0,001) − sin(2 · 0,001) ≈ 9,999997501071444 · 10−10 ,
f (0,001) = 4 sin(0,001) sin2 (0,5 · 0,001) ≈ 9,999997500000248 · 10−10 (mest noggrant).
Det är alltid uttrycket utan subtraktion som är mest noggrannt eftersom det inte kan ske någon
cancellation.
10. x ≈ 2.4
11. x2 = 0.442854, feluppskattning: Ex =
12. x2 = 4.96511, feluppskattning: Ex =
Ez
|f 0 (x2 )|
Ez
|f 0 (x2 )|
=
=
|f (x2 )|
|f 0 (x2 )|
|f (x2 )|
|f 0 (x2 )|
≤ 0.5 · 10−5
≤ 0.5 · 10−5
15. x ≈ 0.56898
16. Här behöver vi inte bry oss om den specifika funktionen, utan vi kan direkt använda oss av
formeln Exn = 2b−a
n+1 .
Exn =
b−a
b−a
⇒ 2n+1 =
⇒ n = log2
2n+1
Exn
b−a
Exn
−1
Då vi har 4 signifikanta decimaler så är Exn = 0,5 · 10−4 .
b−a
n = log2
− 1 = 14,29
Exn
Då vi bara kan göra ett helt antal iterationer så avrundar vi svaret uppåt n = 15.
17.
a) Till exempel x0 = 0,2, x1 ≈ 0,1960133, x2 ≈ 0,19617016, x3 ≈ 0,19616405.
f (x) = x − 0,2 cos(x), Ea ≈
|f (x3 )|
≈ 2,3 · 10−7
|f 0 (x3 )|
b) Ekvationen har flera rötter. Vi väljer den rot som ligger i intervallet [0,π]. I detta intervall
är ekvationen ekvivalent med x = G(x) = arccos( x3 ) som ger konvergens för fixpunktsmetoden. Låt t ex
x0 = 1, x1 = 1.23095, x2 = 1,14799, x3 = 1,17812
Rätt rot: x ≈ 1,17012.