This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
1.3 Guldins regler och Tyngdpunkt Guldins två huvudregler är:
𝑅𝑜𝑡𝑎𝑡𝑖𝑜𝑛𝑑𝑎𝜋𝑒𝑎 = 𝑇𝑦𝑛𝑔𝑑𝑝𝑢𝑛𝑘𝑡𝑒𝑛𝑑 𝑣ä𝑔 ∙ 𝐾𝑢𝜋𝑣𝑙ä𝑛𝑔𝑑𝑒𝑛
𝑅𝑜𝑡𝑎𝑡𝑖𝑜𝑛𝑑𝑣𝑜𝑙𝑦𝑚 = 𝑇𝑦𝑛𝑔𝑑𝑝𝑢𝑛𝑘𝑡𝑒𝑛𝑑 𝑣ä𝑔 ∙ 𝐴𝜋𝑒𝑎𝑛
Reglerna är väldigt användbara framförallt vid rotation som inte är runt en lodrät eller horisontal axel.
För tillämpning av Guldins regler gällande rotationsarea, se avsnittet om rotationsarea för funktionskurvor respektive funktioner i polär form. Nedan beskrivs endast rotationsvolymer enligt Guldins regler.
Rotationsvolym
Rotationsvolym runt x-axeln Teorin om rotation kring x-axeln utgår från att man delar upp funktionen i ”diskar” som man låter rotera i en cirkel. Tyngdpunktens väg blir här en cirkel vars radie är halva ”diskens” höjd.
𝑅𝑜𝑡𝑎𝑡𝑖𝑜𝑛𝑑𝑣𝑜𝑙𝑦𝑚 = 𝑇𝑦𝑛𝑔𝑑𝑝𝑢𝑛𝑘𝑡𝑒𝑛𝑑 𝑣ä𝑔 ∙ 𝐴𝜋𝑒𝑎𝑛 ⇔
𝑑𝑑 = 2𝜋𝑓(𝑥)
2∙ 𝑓(𝑥)𝑑𝑥 ⇒ 𝑑 = �𝑑𝑣 = 𝜋�𝑓(𝑥)2
𝑏
𝑎
𝑑𝑥
Rotationsvolym runt y-axeln
𝑅𝑜𝑡𝑎𝑡𝑖𝑜𝑛𝑑𝑣𝑜𝑙𝑦𝑚 = 𝑇𝑦𝑛𝑔𝑑𝑝𝑢𝑛𝑘𝑡𝑒𝑛𝑑 𝑣ä𝑔 ∙ 𝐴𝜋𝑒𝑎𝑛 ⇔
𝑑𝑑 = 2𝜋𝑥 ∙ 𝑓(𝑥)𝑑𝑥 ⇒ 𝑑 = �𝑑𝑣 = 2𝜋�𝑥𝑓(𝑥)𝑏
𝑎
𝑑𝑥
𝑦
𝑥 𝑋𝑇 = 2𝜋 ∙
𝑓(𝑥)2
𝒚 = 𝒇(𝒙)
𝜋 =𝑓(𝑥)
2
𝑦
𝑥
𝑌𝑇 = 2𝜋𝑥
𝒚 = 𝒇(𝒙)
𝜋 = 𝑥
Sida 9 av 44
Rotationsvolym runt linje För rotationsvolym kring en rät linje (som inte är horisontell eller vertikal) finns det ingen bestämd formel. För att beräkna rotationsvolymen måste man ta fram uttryck för:
1. Areaelement 2. Tyngdpunktens väg.
I exemplet nedan visas hur man tar fram respektive uttryck.
1 – Areaelement Areaelementet är den röda rektangeln vars bas utgörs av skillnaden i x (𝑑𝑥) och höjd som utgörs av skillnaden mellan funktionerna (𝑦2 − 𝑦1).
2 – Tyngdpunktens väg Tyngdpunktens väg innebär den cirkel som roterar runt 𝑦 = 𝑥 som har en radie som ligger på halva areaelementet.
Generell metod: Generellt sätt kan man lösa rotationsproblem med följande metod.
1. Gör en skiss av figuren 2. Ta fram uttryck för ett areaelement. 3. Bestäm uttryck för tyngdpunktens väg. 4. Beräkna integralen.
Exempel: Uppgift 6 från Tenta 2011-06-09 TATA42
Området givet av 4(𝑥 − 1)2 + (𝑦 − 4)2 ≤ 4 roteras ett varv kring linjen 𝑦 = 𝑥. Beräkna rotationskroppens volym.
Uttrycket 4(𝑥 − 1)2 + (𝑦 − 4)2 ≤ 4 beskriver en ellips (med halvaxellängderna 1 respektive 2) med centrum i (1, 4). För att beräkna rotationsvolymen är det smidigast att använda Guldins regel:
𝑅𝑜𝑡𝑎𝑡𝑖𝑜𝑛𝑑𝑣𝑜𝑙𝑦𝑚 = 𝑇𝑦𝑛𝑔𝑑𝑝𝑢𝑛𝑘𝑡𝑒𝑛𝑑 𝑣ä𝑔 ∙ 𝐴𝜋𝑒𝑎𝑛
Area ellips: 𝐴 = 𝑎𝑏𝜋 där a och b är ellipsens halvaxellängder.
Uppgiften skiljer sig från ”typuppgiften” då arean för ellipsen kan beräknas direkt. Därför behöver man heller inte ställa upp någon integral för att beräkna volymen, utan det räcker med att multiplicera kroppens totala area med tyngdpunkten.
Räknelagar för Ordo Syftet med Ordo är att ”samla ihop” felet i approximationen. 𝑂(𝑥2) innebär att felet är en begränsad funktion där den störst betydande termen är 𝑥2.
Lite mer drastiskt kan man se Ordo-funktionen som ett svart hål där man slänger allt som saknar större betydelse.
Då Maclaurinutvecklingar gäller för 𝑥 𝑛ä𝜋𝑎 0 så innebär det i regel att t.ex. 𝑥3 > 𝑥4, eller allmänt: 𝑥𝑚 > 𝑥𝑛, då 𝑚 < 𝑛. Detta förklarar en del regler som intuitivt kan tyckas märkvärdiga.
Generell metod – Gränsvärden 1. Gemensamma nämnare 2. Utveckla nämnare (så lite som möjligt – Lägg märke till graden.) 3. Utveckla täljaren (till samma grad som nämnaren) 4. Jämför nämnare och täljare - Förkorta
Exempel: Uppgift 3a från Tenta 2011-08-25 TATA42
Bestäm konstanten 𝑎 så att gränsvärdet lim𝑥→0ln(1+𝑎𝑥)−2𝑥𝑥 arctan𝑥
existerar ändligt samt beräkna
gränsvärdet.
1. Gemensamma nämnare. Uttrycket är skrivet med gemensam nämnare. 2. Utveckla nämnaren. Använder standardutveckling för arctan𝑥.
𝑥→ ∞ 𝑑å 𝑥 → 0 måste (𝑎 − 2) = 0 för att ett ändligt
gränsvärde ska existera. (𝑎 − 2) = 0 ⇔ 𝑎 = 2
lim𝑥→0
2 − 2𝑥 − 22
2 + 𝑂(𝑥)1 + 𝑂(𝑥2) = lim
𝑥→0
−2 + 𝑂(𝑥)1 + 𝑂(𝑥2) = −2
(Obs! 𝑂(𝑥2) = 𝑏(𝑥) ∙ 𝑥2 där 𝑏(𝑥) är en begränsad funktion och 𝑥2 → 0, alltså går 𝑂(𝑥2) → 0. )
Svar: −2
Sida 13 av 44
2.2 Tillämpning av Lagranges form på resttermen Om man ska beräkna felet vid en approximation går det inte att längre att använda Ordo som restfunktion (man kan aldrig räkna ut vad Ordo blir). Vid problem där man måste beräkna felet ska (måste) man använda Lagranges form på resttermen.
Maclaurinutveckling med Lagranges form på resttermen:
𝑓(𝑥) = 𝑓(0) + 𝑓′(0) ∙ 𝑥 + 𝑓′′(𝑥) ∙𝑥2
2!+. . . +
𝒇(𝒏+𝟏)(𝝃)(𝒏 + 𝟏)!
∙ 𝒙𝒏+𝟏
Resttermen motsvarar alltså felet I approximationen.
Exempel: Från föreläsning 5, TATA42
Approximera integralen
�𝑒𝑥2 − 1
𝑥
1/2
0
𝑑𝑥
så att felet blir mindre än 1/1000.
1. Maclaurinutvecklar funktionen (för komplex för att ta fram primitiv funktion). Obs! Man får chansa på hur långt man vill utveckla (och se hur stort felet blir).
Standard: 𝑒𝑡 = 1 + 𝑡 + 𝑡2
2+ 𝑒𝜉𝑡3
3!, 𝑑ä𝜋 0 ≤ 𝜉 ≤ 𝑡
Här: 𝑒𝑥2 = 1 + 𝑥2 + 𝑥4
2+ 𝑒𝜉𝑥6
3! 𝑑ä𝜋 0 ≤ 𝜉 ≤ 𝑥2
𝑒𝑥2 − 1𝑥
=1 + 𝑥2 + 𝑥4
2 + 𝑒𝜉𝑥63! − 1
𝑥= 𝑥 +
𝑥3
2+𝑒𝜉𝑥5
6
2. Beräknar primitiva funktionen. � 𝑥 +
𝑥3
2+𝑒𝜉𝑥5
6
1/2
0
𝑑𝑥
3. Delar upp intergranden i två delar: 1) Approximationen och 2) resten/felet.
� 𝑥 +𝑥3
2+𝑒𝜉𝑥5
6
1/2
0
𝑑𝑥 = � 𝑥 +𝑥3
2
1/2
0
𝑑𝑥 + �𝑒𝜉𝑥5
6
1/2
0
𝑑𝑥
4. Beräknar resten (det är inte lönt att gå vidare om resten skulle vara större än 1/1000).
�𝑒𝜉𝑥5
6
1/2
0
𝑑𝑥 = �𝑒𝜉𝑥6
36�0
1/2
=𝑒𝜉
36∙ �
12�6− 0 =
𝑒𝜉
2304
Eftersom 0 ≤ 𝜉 ≤ 𝑥2 och 0 ≤ 𝑥 ≤ 12 innebär det att största möjliga värde för 𝑒𝜉 är
𝑒(1/2)2 = 𝑒1/4 < 2 Detta ger oss följande jämförelse:
𝑒𝜉
2304<
22304
=1
1152<
11000
Felet är alltså OK enligt kravet! 5. Beräknar approximationen. � 𝑥 +
𝑥3
2
1/2
0
𝑑𝑥 = �𝑥2
2+𝑥4
8�0
1/2
= �18
+1
128� =
17128
Sida 14 av 44
Exempel: Uppgift 6, från Tenta 2011-05-28 TATA42
Bestäm ett närmevärde för längden av kurvan 𝑦 = ln 𝑥 , 5 ≤ 𝑥 ≤ 10, så att felet är mindre än 11000
.
1. Ställ upp uttryck för kurvlängden. 𝑑 = ��1 + 𝑓′(𝑥)2
𝑏
𝑎
𝑑𝑥
𝑓(𝑥) = ln𝑥 ,𝑓′(𝑥) =1𝑥
⇒ 𝑑 = � �1 + �1𝑥�210
5
𝑑𝑥 = � �1 +1𝑥2
10
5
𝑑𝑥
2. Maclaurinutvecklar uttrycket. Använder Lagranges form på resttermen för att kunna uppskatta felet.
𝑓(𝑥) = 𝑓(0) + 𝑓′(0) ∙ 𝑥 + 𝑓′′(𝑥) ∙𝑥2
2!+. . . +
𝒇(𝒏+𝟏)(𝝃)(𝒏 + 𝟏)!
∙ 𝒙𝒏+𝟏
Vi testar med ordning 2:
𝑓(𝑥) = 𝑓(0) + 𝑓′(0) ∙ 𝑥 +𝒇(𝟐)(𝝃)
(𝟐)!∙ 𝒙𝟐
Substitution: 𝟏𝒙𝟐
= 𝒕 Beräknar derivator: 𝑓(𝑡) = √1 + 𝑡 ⇒ 𝑓(0) = 1
𝑓′(𝑡) =1
2√1 + 𝑡 ⇒ 𝑓′(0) =
12
𝑓′′(𝑡) = −1
4(1 + 𝑡)3/2 ⇒ 𝑓′′(𝜉) = −1
4(1 + 𝜉)3/2
Sätter in i Maclaurinutvecklingen:
√1 + 𝑡 = 1 +12𝑡 −
18(1 + 𝜉)3/2 𝑡
2 ⇔
�1 +1𝑥2
= 1 +1
2𝑥2−
1
8(1 + 𝜉)32∙
1𝑥4
3. Ställer upp integralen. Delar upp den i approximation och rest/fel.
𝑑 = � �1 + �1𝑥�210
5
𝑑𝑥 = � 1 +1
2𝑥2−
18(1 + 𝜉)3/2 ∙
1𝑥4
10
5
𝑑𝑥 =
� 1 +1
2𝑥2
10
5
𝑑𝑥 − �1
8(1 + 𝜉)3/2 ∙1𝑥4
10
5
𝑑𝑥
Approximation | Rest/fel 4. Uppskattar felet.
�1
8(1 + 𝜉)3/2 ∙1𝑥4
10
5
𝑑𝑥 =1
8(1 + 𝜉)3/2 �1𝑥4
10
5
𝑑𝑥
𝜉 är konstant och 0 ≤ 𝜉 ≤ 1/𝑥2, vilket gör att vi kan ställa upp jämförelsen:
18(1 + 𝜉)3/2 �
1𝑥4
10
5
𝑑𝑥 ≤18�
1𝑥4
10
5
𝑑𝑥 =18 �−
13𝑥3�5
10=
18�−
13 ∙ 103
− −1
3 ∙ 53� =
=1
24�
1125
−1
1000� <
11000
Felet är alltså OK! 5. Beräknar approximationen. � 1 +
12𝑥2
10
5
𝑑𝑥 = �𝑥 −1
2𝑥�5
10= �10 −
120� − �5 −
110� = 5 −
120
+1
10=
=100 − 1 + 2
20=
10120
Sida 15 av 44
2.3 Uppskatta värden Med hjälp av Maclaurinutvecklingar kan man uppskatta värde på olika uttryck. Obs! Om man vill uppskatta felet får man inte använda ordoform på resttermen, utan föreslagsvis Lastranges restform.
Exempel: Från föreläsning 5, TATA42
Beräkna värdet av √66 med ett fel mindre än 11000
.
Uträkning Jämförelse med formler
√66 = √64 + 2 = �64 �1 +1
32� = 8 �1 +
132�12
Modell: 𝑓(𝑥) = (1 + 𝑥)1/2
𝑓(𝑥) = 1 +12𝑥 −
14 (1 + 𝜉)−
32
2!∙ 𝑥2 ⇔
𝑓(𝑥) = 1 +12𝑥 −
18
(1 + 𝜉)−32 ∙ 𝑥2
√66 = 8𝑓 �1
32� = 8�1 +
12∙
132
−18
(1 + 𝜉)−32 ∙ �
132�2
� =
= 8�1 +1
64−
18
(1 + 𝜉)−32 ∙ �
132�2
� =
= 8 +18− (1 + 𝜉)−
32 ∙ �
132�2
Där 8 + 1
8 motsvarar approximationen och
(1 + 𝜉)−32 ∙ �
132�2
motsvarar felet, och 0 ≤ 𝜉 ≤ 1/32
(1 + 𝜉)−32 ∙ �
132�2
=1
(1 + 𝜉)32∙
1322
≤1
322=
11024
<1
1000
Felet är alltså inom det tillåtna intervallet.
Maclaurinutveckling med Lastranges restform:
𝑓(𝑥) = 𝑓(0) + 𝑓′(0) ∙ 𝑥 + 𝑓′′(𝑥) ∙𝑥2
2!+ 𝜋(𝑥)
𝑑ä𝜋
𝜋(𝑥) =𝑓(𝑛+1)(𝜉)(𝑛 + 1)!
∙ 𝑥𝑛+1 , 𝑜𝑐ℎ 0 ≤ 𝜉 ≤ 𝑥
Derivator:
𝑓(𝑥) = (1 + 𝑥)1/2, 𝑓(0) = 1
𝑓′(𝑥) =12
(1 + 𝑥)−1/2, 𝑓′(0) =12
𝑓′′(𝑥) = −14
(1 + 𝑥)−3/2
Svar: √66 ≈ 8 + 18
Anmärkning: Man hade kunnat utveckla med standardutveckling, dock behöver man ändå ta fram derivatan för att kunna tillämpa Lastranges restform…
Sida 16 av 44
Exempel: Uppgift 4a, från Tenta 2011-06-09 TATA42
Bestäm ett rationellt tal som approximerar 𝑒−2 med ett fel som är mindre än 110
.
Uträkning Jämförelse med formler Ansats: 𝑓(𝑥) = 𝑒𝑥 Testar att utveckla till ordning 2.
(I facit till uppgiften har de utvecklat 𝑒−2 = (𝑒−1)2 och fått fram att en tillräckligt god approximation är 1/9, men metoden ovan borde vara fullt godkänd. 𝑒−2 ≈ 0,135 vilket ger att alla svar mellan 0,035 och 0,235 borde räknas som godkänt.)
Sida 17 av 44
3. Differentialekvationer
3.1 Teori
Första ordningen – Integrerande faktor När man tar fram lösning till en differentialekvation av första ordningen använder man sig av något som kallas integrerande faktor. Genom att förlänga båda led med den integrerande faktor kan man ta fram lösningen: 𝑦.
Anmärkning. Integrerande faktor är definierat för första ordningens differentialekvationer. Dock kan man använda metoden för differentialekvationer av högre ordningen om de är karaktäristiska. För att se sambandet kan man göra en variabelsubstitution.
Vi ser att vi i 𝑑𝐿 har en funktion bestående av endast 𝑥 samt 𝑦′. 𝐻𝐿 består av en funktion av endast 𝑦 (och en konstant). Differentialekvationen är alltså separabel.
Byter ut 𝑦′ = 𝑑𝑦𝑑𝑥
och flyttar om.
𝑦′ =𝑑𝑦𝑑𝑥
⇒ 𝑥𝑑𝑦𝑑𝑥
= 𝑦2 + 1 ⇒ 1
𝑦2 + 1𝑑𝑦 =
1𝑥𝑑𝑥
Integrerar båda leden. � 1
𝑦2 + 1𝑑𝑦 = �
1𝑥𝑑𝑥 ⇔ arctan𝑦 = ln𝑥 + 𝐶
Använder villkoret för att bestämma konstanten.
Villkor: 𝑦(1) = 1 ⇒ arctan 1 = ln 1 + 𝐶 ⇒ = 𝜋4
⇒
Svar: 𝑦 = tan �ln 𝑥 +𝜋4�
Sida 19 av 44
Homogen lösning En homogen lösning innebär en lösning till differentialekvationen då 𝐻𝐿 = 0. Man ställer upp en karaktäristisk ekvation och beräknar dess rötter.
Differentialekvationer av andra ordningen Rötterna till den karaktäristiska ekvationen ger följande lösningar:
3.2 Exempel Exempel: Uppgift 4 från Tenta 2011-03-19 TATA42
Bestäm alla lösningar till differentialekvationen 𝑦′′′ − 3𝑦′′ + 2𝑦′ = 2𝑥 + 6𝑥2.
Konstaterar främst att det finns en homogen lösning samt en partikulärlösning som kommer att vara ett polynom. Vi tar först fram den honomgena lösningen och därefter partikulärlösningen.
Beskrivning Uträkning 1. Homogen lösning 𝑦′′′ − 3𝑦′′ + 2𝑦′ = 0 har den karaktäristiska ekvationen:
𝑥𝑦′′ + (𝑥 − 1)𝑦′ = 𝑒−𝑥 Ekvationen innehåller första och andra ordningens derivator. Genom att göra en variabelsubstitution kan man forma om differentialekvationen till en första ordningens ekvation och använda sig av integrerande faktor. 𝑦′ = 𝑧, 𝑦′′ = 𝑧′ ⇒ 𝑥𝑧′ + (𝑥 − 1)𝑧 = 𝑒−𝑥
Delar på 𝑥 för att kunna isolera den integrerande faktorn.
𝑧′ +(𝑥 − 1)𝑧
𝑥=𝑒−𝑥
𝑥 ⇔ 𝑧′ + �1 −
1𝑥� 𝑧 =
𝑒−𝑥
𝑥
𝑔 = �1 −1𝑥� , 𝐺 = ��1 −
1𝑥�𝑑𝑥 = 𝑥 − ln 𝑥
Integrerande faktor: 𝑒𝑥−ln𝑥 = 𝑒𝑥 ∙ 𝑒− ln𝑥 =𝑒𝑥
𝑒ln𝑥=𝑒𝑥
𝑥
Multiplicerar båda leden med den integrerande faktorn.
𝑧′ ∙𝑒𝑥
𝑥+ �1 −
1𝑥� 𝑧 ∙
𝑒𝑥
𝑥=𝑒−𝑥
𝑥∙𝑒𝑥
𝑥 ⇔
�𝑧 ∙𝑒𝑥
𝑥 �′
=1𝑥2
Integrerar båda leden och löser därefter ut 𝑧. ��𝑧 ∙
𝑒𝑥
𝑥 �′
= �1𝑥2
⇔ 𝑧 ∙𝑒𝑥
𝑥= −
1𝑥
+ 𝐶1 ⇔ 𝑧 = −1𝑒𝑥
+𝐶1𝑥𝑒𝑥
⇔
𝑧 = −𝑒−𝑥 + 𝐶1𝑥𝑒−𝑥
Tar fram den primitiva funktionen till 𝑧. Obs! Här måste man ta lägga till yterliggare en konstant.
Anmärkning: Variabelsubstitutionen i inledningen är inte nödvändig, utan gör det endast enklare att se sambandet. Man hade lika gärna kunnat gå direkt till:
𝑥𝑦′′ + (𝑥 − 1)𝑦′ = 𝑒−𝑥 ⇔ 𝑦′′ + �1 −1𝑥�𝑦′ =
𝑒−𝑥
𝑥⇔ �𝑦′ ∙
𝑒𝑥
𝑥 �′
=1𝑥2
⇔ 𝑦′ ∙𝑒𝑥
𝑥= �
1𝑥2𝑑𝑥 ⇔
𝑦′ = 𝑥𝑒−𝑥 �1𝑥2𝑑𝑥 ⇔ 𝒚 = ��𝒙𝒆−𝒙 �
𝟏𝒙𝟐𝒅𝒙�𝒅𝒙
Sida 24 av 44
Obs! Specialfall Ett märkligt fall är när den homogena lösningen till en differentialekvation sammanfaller med partikulärlösningen.
T.ex. differentialekvationen 𝒚′′ + 𝟗𝒚 = 𝐜𝐨𝐬𝟑𝒙.
Den homogena lösningen beräknas: 𝜋2 + 9 = 0 ⇔ 𝜋 = ±𝑖3 ⇒ 𝑦ℎ = 𝐶1 cos 3𝑥 + 𝐶2 sin 3𝑥
Vanligtvis skulle man nu göra en ansats 𝒚 = 𝒅𝐜𝐨𝐬𝒌𝒙 + 𝒅𝐬𝐢𝐧𝒌𝒙 för att ta fram partikulärlösningen, men detta hade gett samma svar som den homogena lösningen. För att beräkna partikulärlösningen här är det bästa sättet att göra en komplex ansats.
Den komplexa exponentialfunktionen definieras: 𝒚 = 𝒆𝒊𝒂𝒙 = 𝐜𝐨𝐬𝒂𝒙 + 𝒊 𝐬𝐢𝐧𝒂𝒙. Här kan man utnyttja att funktionen har en reell del samt en imaginär del:
• Re�𝑒𝑖3𝑥� = cos 3𝑥 • Im�𝑒𝑖3𝑥� = sin 3𝑥.
Vi sätter upp en ny ekvation: 𝒖′′ + 𝟗𝒖 = 𝒆𝟑𝒊𝒙. Om 𝒖𝒑 står för partikulärlösningen så kan man med resonemanget ovan sluta sig till att 𝒚𝒑 = 𝑹𝒆(𝒖𝒑)
Nu löser man differentialekvationen 𝑢′′ + 9𝑢 = 𝑒3𝑖𝑥 med en vanlig ansats: 𝒖 = 𝒛𝒆𝟑𝒊𝒙
2. Tar fram partikulärlösning. Obs! Här hade man vanligtvis gjort en trigonometrisk ansats. Eftersom att den homogena lösningen är en trigonometrisk ansats får man istället göra en komplex ansats.
𝑦′′ + 4𝑦 = sin 2𝑥
𝑒2𝑖𝑥 = cos 2𝑥 + 𝑖 sin 2𝑥 ⇔ 𝐼𝑚�𝑒2𝑖𝑥� = sin 2𝑥
Ställer upp en stödekvation:
𝑢′′ + 4𝑢 = 𝑒2𝑖𝑥
Ansats: 𝒖 = 𝒛𝒆𝟐𝒊𝒙
𝑢′′ = (𝑧′′ + 𝑖4𝑧′ − 4𝑧)𝑒2𝑖𝑥 Sätter in i ekvationen:
Det man ska kunna: • Bestämma om en serie är konvergent eller inte • Beräkna värdet av serie • Bestämma fel vid approximation
I avsnitt fyra beskrivs endast hur man bestämmer om en serie är konvergent eller divergent. Hur man beräknar värdet av serier och bestämmer fel vid approximationer beskrivs i avsnitt 5 (potensserier) då metoderna i princip är samma.
4.1 Teori
Serier och integraler Vid bedömning av konvergens och divergens kan integraler och serier behandlas snarlikt.
Sats: Integralkriteriet.
𝑎) 𝑂𝑚 � 𝑓(𝑥)∞
1
𝑑𝑥 ä𝜋 𝑘𝑜𝑛𝑣𝑒𝜋𝑔𝑒𝑛𝑡, 𝑑å ä𝜋 �𝑓(𝑘)∞
𝑘=1
𝑘𝑜𝑛𝑣𝑒𝜋𝑔𝑒𝑛𝑡.
𝑏) 𝑂𝑚 � 𝑓(𝑥)∞
1
𝑑𝑥 ä𝜋 𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡, 𝑑å ä𝜋 �𝑓(𝑘)∞
𝑘=1
𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡.
Begrepp:
𝐷𝑒𝑓𝑖𝑛𝑖𝑡𝑖𝑜𝑛 𝑑𝑒𝜋𝑖𝑒:�𝑎𝑘
∞
𝑘=1
= 𝑎1 + 𝑎2 + 𝑎3 + 𝑎4+. ..
Konvergens: När delar från olika håll närmar sig varandra. En konvergent serie erhåller alltså ett värde när antalet termer går mot oändligheten
Divergens: Isärgående; spridning åt olika håll. En divergent serie får inte något ändligt värde när serien går mot oändligheten.
OBS! Innan man påbörjar någon beräkning måste man ha klart för sig om serien är positiv, varierande eller alternerande.
Typ av serie Definition Exempel Positiv serie
�𝑎𝑘
∞
𝑘=1
= 𝑎1 + 𝑎2 + 𝑎3+. . .𝑑ä𝜋 𝑎𝑘 ≥ 0
�1𝑘
∞
𝑘=1
= 1 + 2 + 3 + 4. ..
”Varierande” serie
Har både positiva och negativa termer. �
sin𝑘𝑘2
∞
𝑘=1
Alternerande serie �𝑎𝑘
∞
𝑘=1
𝑑ä𝜋 𝑎𝑘 > 0 𝑓ö𝜋 𝑢𝑑𝑑𝑎 𝑘 = 1,3,5, …
𝑜𝑐ℎ 𝑎𝑘 < 0 𝑓ö𝜋 𝑗ä𝑚𝑛𝑎 𝑘 = 2,4,6, … (eller tvärtom)
�(−1)𝑘−1
𝑘
∞
𝑘=1
= 1 −12
+13−
14
+ ⋯
Sida 27 av 44
Divergenstest (Gäller för alla typer av serier) För en konvergent serie måste lim𝑛→∞ 𝑎𝑛 = 0. Om inte växer serien hela tiden och går då mot oändligheten. Kriteriet innebär ett enkelt sätt att testa om en serie är divergent.
𝑂𝑚 𝑙𝑖𝑚 𝑛→∞
𝑎𝑛 ≠ 0 ä𝜋 𝑑𝑒𝜋𝑖𝑒𝑛 �𝑎𝑘
∞
𝑘=𝑛
𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡.
Geometriska serien (Kan vara positiv eller alternerande)
�𝑎𝑞𝑛∞
𝑛=0
= 𝑎 + 𝑎𝑞1 + 𝑎𝑞2 + 𝑎𝑞3+. . .�=
𝑎1 − 𝑞
𝑜𝑚 |𝑞| < 1
𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡 𝑜𝑚 |𝑞| ≥ 1
Förhållandet ter sig ganska självklart om man jämför med en exponentialfunktion:
𝑦 = 𝑎𝑞𝑥 → ∞ 𝑑å 𝑥 → ∞ 𝑜𝑚 |𝑞| > 1
𝑦 = 𝑎𝑞𝑥 → 0 𝑑å 𝑥 → ∞ 𝑜𝑚 |𝑞| < 1
Jämförelsesatser (Endast för positiva serier) Jämförelsesats I
𝐹ö𝜋�𝑎𝑘
∞
𝑘=1
𝑜𝑐ℎ �𝑏𝑘
∞
𝑘=1
𝑑ä𝜋 0 ≤ 𝑎𝑘 ≤ 𝑏𝑘 𝑔ä𝑙𝑙𝑒𝜋
𝒂) 𝑂𝑚 �𝑏𝑘
∞
𝑘=1
ä𝜋 𝑘𝑜𝑛𝑣𝑒𝜋𝑔𝑒𝑛𝑡 ⇒ �𝑎𝑘
∞
𝑘=1
𝑘𝑜𝑛𝑣𝑒𝜋𝑔𝑒𝑛𝑡 �𝑜𝑐ℎ �𝑎𝑘
∞
𝑘=1
≤ �𝑏𝑘
∞
𝑘=1
�
𝒃) 𝑂𝑚 �𝑏𝑘
∞
𝑘=1
ä𝜋 𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡 ⇒ �𝑎𝑘
∞
𝑘=1
𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡
Jämförelsesats II
𝐹ö𝜋�𝑎𝑘
∞
𝑘=1
𝑜𝑐ℎ �𝑏𝑘
∞
𝑘=1
𝑑ä𝜋 𝑎𝑘 ≥ 0, 𝑏𝑘 ≥ 0 𝑜𝑐ℎ 𝑎𝑘 = 𝑐𝑘 ∙ 𝑏𝑘
𝑑ä𝜋 𝑐𝑘ä𝜋 𝑒𝑡𝑡 ä𝑛𝑑𝑙𝑖𝑔𝑡 𝑔𝜋ä𝑛𝑑𝑣ä𝜋𝑑𝑒 𝑔ä𝑙𝑙𝑒𝜋 ∶
𝒂) 𝑂𝑚 �𝑏𝑘
∞
𝑘=1
ä𝜋 𝑘𝑜𝑛𝑣𝑒𝜋𝑔𝑒𝑛𝑡 ⇒ �𝑎𝑘
∞
𝑘=1
𝑘𝑜𝑛𝑣𝑒𝜋𝑔𝑒𝑛𝑡
𝒃) 𝑂𝑚 �𝑏𝑘
∞
𝑘=1
ä𝜋 𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡 ⇒ �𝑎𝑘
∞
𝑘=1
𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡
𝑦
𝑥
𝑦
𝑥
Sida 28 av 44
Jämförelseserier:
𝐀)�1𝑥𝑎
𝑑𝑥∞
1
ä𝜋 �𝑘𝑜𝑛𝑣𝑒𝜋𝑔𝑒𝑛𝑡 om 𝑎 > 1𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡 om 𝑎 ≤ 1
𝐁)�1𝑥𝑎
𝑑𝑥1
0
ä𝜋 �𝑘𝑜𝑛𝑣𝑒𝜋𝑔𝑒𝑛𝑡 om 𝑎 < 1𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡 om 𝑎 ≥ 1
Absolutkonvergens (Säger ingenting om serien är positiv… Bra om serien är varierande)
𝑂𝑚 �𝑎𝑘
∞
𝑘=1
ä𝜋 𝑒𝑛 𝑑𝑒𝜋𝑖𝑒 𝑔ä𝑙𝑙𝑒𝜋:
𝒂) �|𝑎𝑘|∞
𝑘=1
𝑘𝑜𝑛𝑣𝑒𝜋𝑔𝑒𝑛𝑡 ⇒ �𝑎𝑘
∞
𝑘=1
𝑘𝑜𝑛𝑣𝑒𝜋𝑔𝑒𝑛𝑡
𝒃) ��𝑎𝑘
∞
𝑘=1
� ≤ �|𝑎𝑘|∞
𝑘=1
Obs! Bara för att en serie inte är absolutkonvergent behöver den inte vara divergent…
Rotkriteriet & Kvotkriteriet (Gäller för alla typer av serier)
𝐹ö𝜋 𝑑𝑒𝜋𝑖𝑒𝑛 �𝑎𝑘
∞
𝑘=1
𝑔ä𝑙𝑙𝑒𝜋 𝑎𝑡𝑡 𝑜𝑚
𝑎) �|𝑎𝑘|𝑘 → 𝑄 𝑑å 𝑘 → ∞
𝑏)|𝑎𝑘+1|
|𝑎𝑘| → 𝑄 𝑑å 𝑘 → ∞
𝑑ä𝜋 𝑄 ä𝜋 𝑒𝑡𝑡 𝜋𝑒𝑒𝑙𝑙𝑡 𝑡𝑎𝑙 ≥ 0 𝑔ä𝑙𝑙𝑒𝜋
𝑄 < 1 ⇒ 𝑆𝑒𝜋𝑖𝑒𝑛 ä𝜋 𝑎𝑏𝑑𝑜𝑙𝑢𝑡𝑘𝑜𝑛𝑣𝑒𝜋𝑔𝑒𝑛𝑡
𝑄 > 1, (𝑄 = ∞) ⇒ 𝑆𝑒𝜋𝑖𝑒𝑛 ä𝜋 𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡
Obs! Om 𝑄 = 1 säger testerna ingenting.
Testerna fyller samma funktion och fungerar på likadant sätt. Vilket man väljer beror av hur serien ser ut.
Liebniz Kriterium (Endast alternerande serier) Liebniz kriterium innebär att en serie är konvergent om:
Kriterium Samband • Serien är alternerande … • Termernas belopp avtar (för stora 𝑛) Obs! Kan visas på två sätt:
|𝑎𝑛+1| < |𝑎𝑛| 𝑓′(𝑥) < 0
• Termerna går mot noll lim𝑛→∞
𝑎𝑛 = 0
Obs! Bara för att kriterierna inte uppfylls måste inte serien vara divergent…
Sida 29 av 44
4.2 Bestämma om en serie är konvergent eller divergent Strategi Man kan inte använda alla verktyg på alla serier. Nedan beskrivs vilka verktyg man kan använda till vilka serier (det vanligaste verktyget till respektive typ är fetmarkerat).
Metod Bevisar Divergenstestet Bevisar om en serie är divergent.
Kan inte visa att en serie är konvergent. Jämförelsesatser Kan bevisa både konvergens respektive divergens. Absolutkonvergent En serie som en absolutkonvergent är konvergent.
En serie som inte är absolutkonvergent behöver inte vara divergent. Liebniz Kriterium Kan bevisa att en serie är konvergent.
Serien behöver inte vara divergent bara för att kriterierna inte uppfylls. Rot- & Kvotkriteriet Kan visa både divergens respektive konvergens.
Om 𝑄 = 1 säger testet ingenting.
Exempel Exempel: Uppgift 5b från Tenta 2012-01-14 TATA42
Avgör om �(−1)𝑛
𝑛 − √𝑛
∞
𝑛=2
är konvergent.
Vi ser att serien är alternerande, då den innehåller termen (−1)𝑛. Första testet blir därför att titta på Liebniz kriterium.
Liebniz Kriterium
1. Serien är alternerande.
2. Termerna går mot noll.
𝐥𝐢𝐦𝒏→∞
(−1)𝑛
𝑛 − √𝑛= 𝐥𝐢𝐦
𝒏→∞ 𝟏𝒏∙
(−𝟏)𝒏
𝟏 − 𝟏√𝒏
= 𝟎
3. Termernas belopp avtar (för stora 𝑛)
|𝑎𝑛+1| < |𝑎𝑛| ⇔ �(−1)𝑛+1
𝑛 + 1 − √𝑛 + 1� < �
(−1)𝑛
𝑛 − √𝑛� ⇔
1𝑛 + 1 − √𝑛 + 1
<1
𝑛 − √𝑛⇔
𝑛 − √𝑛 < 𝑛 + 1 − √𝑛 + 1 ⇔ √𝑛 + 1 < 1 + √𝑛
Svar: Konvergent enligt Liebniz.
Sida 30 av 44
Exempel: Uppgift 5a från Tenta 2012-01-14 TATA42
Avgör om �𝑥2
𝑥2 + 𝑥
∞
1
𝑑𝑥 är konvergent.
Analogt med teorin för serier så gäller att divergenstestet även för integraler: Om lim𝑥→∞ 𝑓(𝑥) ≠ 0 så är integralen divergent (annars hade den ju gått mot oändligheten…).
𝐃𝐢𝐯𝐞𝐫𝐠𝐞𝐧𝐬𝐭𝐞𝐬𝐭: lim𝑥→∞
𝑥2
𝑥2 + 𝑥= lim
𝑥→∞
𝑥2
𝑥2(1 + 1/𝑥)= lim
𝑥→∞
11 + 1/𝑥
= 1 ≠ 0 ⇒ 𝐷𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡
𝐀𝐥𝐭𝐞𝐫𝐧𝐚𝐭𝐢𝐯 𝐥ö𝐬𝐧𝐢𝐧𝐠:�𝑥2
𝑥2 + 𝑥
∞
1
𝑑𝑥 = �1
1 + 1/𝑥
∞
1
𝑑𝑥 ≥ �1
2𝑥𝑑𝑥
∞
1
, 𝑑𝑜𝑚 ä𝜋 𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡. (𝐽ä𝑚𝑓ö𝜋𝑒𝑙𝑑𝑒𝑑𝑎𝑡𝑑 𝐼)
𝐀𝐥𝐭𝐞𝐫𝐧𝐚𝐭𝐢𝐯 𝐥ö𝐬𝐧𝐢𝐧𝐠:�𝑥2
𝑥2 + 𝑥
∞
1
𝑑𝑥 ≥ �𝑥2
𝑥2 + 𝑥2
∞
1
𝑑𝑥 = �12𝑑𝑥
∞
1
, 𝑑𝑜𝑚 ä𝜋 𝑑𝑖𝑣𝑒𝜋𝑔𝑒𝑛𝑡.
Sida 31 av 44
Exempel: Uppgift 6 från Tenta 2011-08-25 TATA42
För vilka 𝑎 ∈ 𝑹 är den generaliserade integralen � cos(𝑒𝑎𝑥)𝑑𝑥∞
0
konvergent?
Integralen är generaliserad vid 𝑥 → ∞. Den inre funktionen 𝑡 = 𝑒𝑎𝑥 gör det naturligt att dela upp undersökningen i fall där 𝑎 är större, lika med eller mindre än noll (eftersom 𝑡 = 𝑒𝑎𝑥 då är växande, konstant respektive avtagande).
Obs! Satsen säger ingenting om hur serien beter sig då |𝒙| = 𝑹, detta måste undersökas separat.
Derivata och primitiv funktion För en potensserie med konvergensradien R och |𝑥| < 𝑅 (den är konvergent) gäller följande:
𝑓(𝑥) = �𝑐𝑘
∞
𝑘=0
𝑥𝑘 = 𝑐0 + 𝑐1𝑥 + 𝑐2𝑥2 + 𝑐3𝑥3+. ..
𝑓′(𝑥) = �𝑘 ∙ 𝑐𝑘
∞
𝑘=1
𝑥𝑘−1 = 𝑐1 + 2𝑐2𝑥 + 3𝑐3𝑥2+. ..
𝐹(𝑥) = �𝑐𝑘𝑥𝑘+1
𝑘 + 1
∞
𝑘=0
= 𝑐0𝑥 +𝑐1𝑥2
2+𝑐2𝑥3
3+𝑐3𝑥4
4+. ..
Alla serierna har samma konvergensradie: R.
Sida 33 av 44
Värdet av en potensserie Med hjälp av derivata och primitiva funktioner kan man beräkna värdet av potensserier. Man utnyttjar formeln för beräkning av en geometrisk summa och deriverar/tar fram primitiv funktion implicit.
För konvergenta serier där |𝒙| < 𝟏 gäller:
𝑓(𝑥) = �𝑥𝑘∞
𝑘=0
=1
1 − 𝑥
𝑓′(𝑥) = �𝑘 ∙ 𝑥𝑘−1∞
𝑘=1
=1
(1 − 𝑥)2
𝐹(𝑥) = �𝑥𝑘+1
𝑘 + 1
∞
𝑘=0
= �𝑥𝑘
𝑘
∞
𝑘=1
= − 𝑙𝑛(1 − 𝑥)
Maclaurinserier De vanliga Maclaurinutvecklingarna kan skrivas som serier. Sambanden kan användas för att bland annat beräkna värdet av en serie.
𝑒𝑥 = �𝑥𝑛
𝑛!
∞
𝑘=0
𝑑𝑖𝑛 𝑥 = �(−1)𝑘𝑥2𝑘+1
(2𝑘 + 1)!
∞
𝑘=0
𝑐𝑜𝑑 𝑥 = �(−1)𝑘𝑥2𝑘
(2𝑘)!
∞
𝑘=0
𝑙𝑛(1 + 𝑥) = �(−1)𝑘−1𝑥𝑘
𝑘
∞
𝑘=0
, 𝑓ö𝜋 − 1 < 𝑥 ≤ 1
(1 + 𝑥)𝑎 = ��𝑎𝑘�𝑥𝑘
∞
𝑘=0
𝑎𝜋𝑐𝑡𝑎𝑛 𝑥 = �(−1)𝑘𝑥2𝑘+1
2𝑘 + 𝑥
∞
𝑘=0
, 𝑓ö𝜋 − 1 ≤ 𝑥 ≤ 1
Sida 34 av 44
5.2 Konvergensradie Exempel: Uppgift 6a från Tenta 2012-01-14 TATA42
Avgör för vilka 𝑥 ∈ 𝐑 serien �𝑛𝑥𝑛
7𝑛
∞
𝑛=1
är konvergent.
Ställer upp kvotkriteriet för att ta fram ett värde på konvergensradien.
𝐊𝐯𝐨𝐭𝐤𝐫𝐢𝐭𝐞𝐫𝐢𝐞𝐭: |𝑎𝑘+1|
|𝑎𝑘| → 𝑄 𝑑å 𝑘 → ∞
�𝑛𝑥𝑛+1
7𝑛+1� �
𝑛𝑥𝑛
7𝑛 �� =
𝑛|𝑥𝑛𝑥|7𝑛 ∙ 7
∙7𝑛
𝑛|𝑥𝑛| =|𝑥|7
= 𝑄 < 1 ⇔ |𝑥| < 7
Serien är alltså absolutkonvergent för |𝑥| < 7, och således divergent för |𝑥| > 7. Dock måste man ta reda på vad som händer i ändpunkterna: 𝑥 = ±7
𝒙 = 𝟕 ⇒ �𝑛7𝑛
7𝑛
∞
𝑛=1
= �𝑛∞
𝑛=1
= 1,2,3, …
Växande och positiv: Divergent!
𝒙 = −𝟕 ⇒ �𝑛(−7)𝑛
7𝑛
∞
𝑛=1
= �𝑛(−1)𝑛∞
𝑛=1
= −1, 2,−3, 4 …
Divergent!
Svar: Serien är konvergent för −7 < 𝑥 < 7
Exempel: Uppgift 4b från Tenta 2012-01-14 TATA42
För vilka 𝑥 ∈ 𝐑 serien �4𝑛+1𝑥2𝑛
𝑛
∞
𝑛=1
är konvergent. ?
Ställer upp kvotkriteriet för att ta fram ett värde på konvergensradien.
𝐊𝐯𝐨𝐭𝐤𝐫𝐢𝐭𝐞𝐫𝐢𝐞𝐭: |𝑎𝑘+1|
|𝑎𝑘| → 𝑄 𝑑å 𝑘 → ∞
�4𝑛+1+1𝑥2(𝑛+1)
𝑛 + 1� �
4𝑛+1𝑥2𝑛
𝑛�� =
4𝑛 ∙ 42 ∙ 𝑥2𝑛 ∙ 𝑥2
𝑛 + 1∙
𝑛4𝑛 ∙ 4 ∙ 𝑥2𝑛
=
=𝑛4𝑥2
𝑛 + 1=
4𝑥2
1 + 1/𝑛→ 4𝑥2 = 𝑄 < 1 ⇔ 4𝑥2 < 1 ⇔ |𝑥| <
12
Således absolutkonvergent för |𝑥| < 1
2, och divergent för
|𝑥| > 12. Kontrollerar nu
ändpunkterna: 𝑥 = ±1/2
𝒙 = ±𝟏𝟐⇒ �
4𝑛+1 �± 12�
2𝑛
𝑛
∞
𝑛=1
= �4𝑛 ∙ 4𝑛
∞
𝑛=1
∙1
4𝑛= �
4𝑛
∞
𝑛=1
som är divergent.
Svar: Serien är konvergent för −1/2 < 𝑥 < 1/2
Sida 35 av 44
5.3 Beräkna värde av potensserie Exempel: Uppgift 6b från Tenta 2012-01-14 TATA42
Beräkna �𝑛
4𝑛−1
∞
𝑛=1
Vi har 𝑛 både som koefficient och exponent vilket leder tankarna till derivata. Målet nu är att matcha serien med en formel för att därefter kunnat beräkna värdet.
Uträkning Jämförelse med formler
�𝑛
4𝑛−1
∞
𝑛=1
= �𝑛 ∙ �14�𝑛−1∞
𝑛=1
Serien stämmer överrens med formeln för 𝑓′(𝑥) för en allmän potensserie.
�𝑛 ∙ �14�𝑛−1∞
𝑛=1
=1
�1 − 14�
2 =19
16=
169
𝑓(𝑥) = �𝑥𝑘∞
𝑘=0
=1
1 − 𝑥
𝑓′(𝑥) = �𝑘 ∙ 𝑥𝑘−1∞
𝑘=1
=1
(1 − 𝑥)2
𝐹(𝑥) = �𝑥𝑘+1
𝑘 + 1
∞
𝑘=0
= �𝑥𝑘
𝑘
∞
𝑘=1
= − 𝑙𝑛(1 − 𝑥)
Svar: 16/9
Exempel: Från föreläsning 11, TATA42
Beräkna �1𝑛5𝑛
∞
𝑛=1
I serien finns 𝑛 som nämnare och som exponent. Tankarna går då till primitiva funktioner av potensserier.
Uträkning Jämförelse med formler
�1𝑛5𝑛
∞
𝑛=1
= �1𝑛∙ �
15�𝑛∞
𝑛=1
Serien stämmer överrens med formeln för 𝐹(𝑥) för en allmän potensserie.
�1𝑛∙ �
15�𝑛∞
𝑛=1
= − ln �1 −15� = − ln �
45� = ln �
54�
𝑓(𝑥) = �𝑥𝑘∞
𝑘=0
=1
1 − 𝑥
𝑓′(𝑥) = �𝑘 ∙ 𝑥𝑘−1∞
𝑘=1
=1
(1 − 𝑥)2
𝐹(𝑥) = �𝑥𝑘+1
𝑘 + 1
∞
𝑘=0
= �𝑥𝑘
𝑘
∞
𝑘=1
= − 𝑙𝑛(1 − 𝑥)
Svar: ln �54�
Sida 36 av 44
Exempel: Uppgift 5c från Tenta 2011-06-09 TATA42
Bestäm summan av serien �𝑛 + 2
3𝑛
∞
𝑛=0
Vi skriver om serien och delar upp den:
�𝑛 + 2
3𝑛
∞
𝑛=0
= �(𝑛 + 2) ∙ �13�𝑛∞
𝑛=0
= �𝑛�13�𝑛
+ 2 �13�𝑛∞
𝑛=0
= �𝑛�13�𝑛∞
𝑛=0
+ 2 ��13�𝑛∞
𝑛=0
Nu kan beräkna värdet av serierna var för sig.
2 ��13�𝑛∞
𝑛=0
= 2 ∙1
1 − �13�
= 2 ∙123
= 3
�𝑛�13�𝑛∞
𝑛=0
När 𝑛 finns både som koefficient och som exponent ser man att uttrycket påminner om en derivata. Nu gäller det att ”forma om” uttrycket så att det stämmer helt med uttrycken för funktion/derivata.
Uträkning Jämförelse med formler
�𝑛�13�𝑛∞
𝑛=0
Jämfört med den deriverade formen av en potensserie ser vi att det skiljer sig två saker:
a) Serien börjar vid 𝑛 = 0 istället för 𝑘 = 1 b) Exponenten är 𝑛 istället för 𝑘 − 1
Vi vill nu forma om uttrycket så det stämmer överrens.
�𝑛�13�𝑛∞
𝑛=0
= 0 �13�0
+ 1 �13�1
+ 2 �13�2
+. . . =
�13� �1 + 2 �
13� + 3 �
13�2
+. . .� =13∙ � 𝑛 �
13�𝑛−1∞
𝑛=1
13∙ � 𝑛�
13�𝑛−1∞
𝑛=1
=13∙
1
�1 − 13�
2 =13∙
149
=34
𝑓(𝑥) = �𝑥𝑘∞
𝑘=0
=1
1 − 𝑥
𝑓′(𝑥) = �𝑘 ∙ 𝑥𝑘−1∞
𝑘=1
=1
(1 − 𝑥)2
𝐹(𝑥) = �𝑥𝑘+1
𝑘 + 1
∞
𝑘=0
= �𝑥𝑘
𝑘
∞
𝑘=1
= − 𝑙𝑛(1 − 𝑥)
𝐒𝐯𝐚𝐫: �𝑛�13�𝑛∞
𝑛=0
+ 2 ��13�𝑛∞
𝑛=0
=34
+ 3 =𝟏𝟓𝟒
Sida 37 av 44
5.4 Differentialekvationer och potensserier Genom att se en potensserie som en funktion kan man lösa differentialekvationer som innehåller produkter av 𝑥-termer och funktioner eller derivator.
Lösningsmetoden beskrivs först med ett ganska enkelt exempel (som man hade kunnat lösa mycket smidigare) och appliceras därefter på ett svårt exempel, där ansättning som en potensserie är enda (bästa) tillvägagångssättet.
Exempel: Uppgift 10.22 från boken (Forsling)
Bestäm alla lösningar i en potensserieform till differentialekvationen 𝑦′′ − 𝑦 = 𝑥 med villkoret att 𝑦 = 1 och 𝑦′ = 0 för 𝑥 = 0.
1. Ansätter 𝑦 som en potensserie. 𝑦 = �𝑐𝑘
∞
𝑘=0
𝑥𝑘 = 𝑐0 + 𝑐1𝑥 + 𝑐2𝑥2 + 𝑐3𝑥3 + 𝑐4𝑥4. ..
2. Bestämmer 𝑦′′ 𝑦′ = �𝑘𝑐𝑘
∞
𝑘=1
𝑥𝑘−1 = 𝑐1 + 2𝑐2𝑥1 + 3𝑐3𝑥2 + 4𝑐4𝑥3. ..
𝑦′′ = �𝑘(𝑘 − 1)𝑐𝑘𝑥𝑘−2∞
𝑘=2
= 2𝑐2 + 6𝑐3𝑥 + 12𝑐4𝑥2+. ..
3. Anpassar 𝑦′′ genom att göra ändringar i termerna (vi vill att summan ska börja vid 𝑘 = 0).
𝑦′′ = �𝑘(𝑘 − 1)𝑐𝑘𝑥𝑘−2∞
𝑘=2
= �(𝑘 + 2)(𝑘 + 1)𝑐𝑘+2𝑥𝑘∞
𝑘=0
Tänk: Summorna ska börja med samma värde. 𝑘 = 2 insatt i det vänstra uttrycket ska ge samma som 𝑘 = 0 insatt i det högra.
4. Sätter in uttrycken i differentialekvationen.
𝑦′′ − 𝑦 = 𝑥 ⇔
�(𝑘 + 2)(𝑘 + 1)𝑐𝑘𝑥𝑘∞
𝑘=0
−� 𝑐𝑘
∞
𝑘=0
𝑥𝑘 = 𝑥 ⇔
��(𝑘 + 2)(𝑘 + 1)𝑐𝑘+2𝑥𝑘 − 𝑐𝑘𝑥𝑘�∞
𝑘=0
5. Jämför termer mot termer. Detta görs smidigast genom att utveckla summan i 𝑑𝐿 och jämföra med 𝐻𝐿.
Villkor: 𝑦 = 1 och 𝑦′ = 0 för 𝑥 = 0 Sätter in i serierna:
𝑦(𝟎) = �𝑐𝑘
∞
𝑘=0
𝟎𝑘 = 𝑐0 + 𝑐1𝟎+ 𝑐2𝟎2 + … = 𝑐0 = 𝟏
𝑦′(𝟎) = �𝑘𝑐𝑘
∞
𝑘=1
𝟎𝑘−1 = 𝑐1 + 2𝑐2𝟎1 + 3𝑐3𝟎2+. . . = 𝑐1 = 𝟎
7. Bestämmer uttryck för 𝑐𝑘 med hjälp av förhållandena i steg 5 och värdena i steg 6. �
2𝑐2 − 𝑐0 = 06𝑐3 − 𝑐1 = 1
(𝑘 + 2)(𝑘 + 1)𝑐𝑘+2 − 𝑐𝑘 = 0⇔
⎩⎨
⎧𝑐2 = 1/2𝑐3 = 1/6
𝑐𝑘+2 =𝑐𝑘
�(𝑘 + 1)(𝑘 + 2)�
Sida 38 av 44
𝑘 = 2 ⇒ 𝑐𝟒 =𝑐2
3 ∙ 4=
12∙
13 ∙ 4
=1𝟒!
𝑘 = 3 ⇒ 𝑐𝟓 =𝑐3
4 ∙ 5=
16∙
14 ∙ 5
=1𝟓!
𝑐𝒌 =1𝒌!
8. Skriver upp ett uttryck för funktionen. 𝑦 = 𝑐0 + 𝑐1𝑥 + 𝑐2𝑥2 + 𝑐3𝑥3 + 𝑐4𝑥4+. . ., 𝑑ä𝜋 �
𝑐0 = 1𝑐1 = 0𝑐𝑘 = 1
𝑘!
⇒ 𝑦 = 1 +𝑥2
2!+𝑥3
3!+𝑥4
4!+. ..
9. Skriver som en serie. 𝒚 = 1 +
𝑥2
2!+𝑥3
3!+𝑥4
4!+. . . = 1 + �
𝑥𝑘
𝑘
∞
𝑘=2
𝐒𝐯𝐚𝐫: 𝒚 = 1 + �𝑥𝑘
𝑘
∞
𝑘=2
(10. Förenkla svaret.) Enligt Maclaurinutveckling ser man att summan motsvarar en exponentialfunktion.
𝑒𝑥 = �𝑥𝑘
𝑘!
∞
𝑘=0
= 1 + 𝑥 + �𝑥𝑘
𝑘
∞
𝑘=2
⇔ �𝑥𝑘
𝑘
∞
𝑘=2
= 𝑒𝑥 − 1 − 𝑥
1 + �𝑥𝑘
𝑘
∞
𝑘=2
= 1 + (𝑒𝑥 − 1 − 𝑥) = 𝑒𝑥 − 𝑥
𝐒𝐯𝐚𝐫: 𝒚 = 𝑒𝑥 − 𝑥
Anmärkning 1: Definitionen av fakultet ger 0! = 1
Anmärkning 2: Det svåraste steget vid beräkningen är att komma fram till ett uttryck för 𝑐𝑘. Här finns det ingen exakt metod, utan man får testa sig fram med olika värden på 𝑘 för att kunna sluta sig till en lösning.
Sida 39 av 44
Exempel: Uppgift 6 från Tenta 2012-03-10 TATA42
Låt differentialekvationen 𝑦′′ − 2𝑥𝑦′ − 4𝑦 = 0, med begynnelsevillkoren 𝑦(0) och 𝑦′(0) = 1, med en potensserieansats. Uttryck även lösningen med hjälp av elementära funktioner.
1. Ansätter 𝑦 som en potensserie. 𝑦 = �𝑐𝑘
∞
𝑘=0
𝑥𝑘 = 𝑐0 + 𝑐1𝑥 + 𝑐2𝑥2 + 𝑐3𝑥3 + 𝑐4𝑥4. ..
2. Bestämmer derivatorna: 𝑦′ och 𝑦′′ 𝑦′ = �𝑘𝑐𝑘
∞
𝑘=1
𝑥𝑘−1 = 𝑐1 + 2𝑐2𝑥1 + 3𝑐3𝑥2 + 4𝑐4𝑥3. ..
𝑦′′ = �𝑘(𝑘 − 1)𝑐𝑘𝑥𝑘−2∞
𝑘=2
= 2𝑐2 + 6𝑐3𝑥 + 12𝑐4𝑥2+. ..
3. Sätter in uttrycken i differentialekvationen.
𝑦′′ − 2𝑥𝑦′ − 4𝑦 = 0 ⇔
�𝑘(𝑘 − 1)𝑐𝑘𝑥𝑘−2∞
𝑘=2
− 2𝑥�𝑘𝑐𝑘
∞
𝑘=1
𝑥𝑘−1 − 4�𝑐𝑘
∞
𝑘=0
𝑥𝑘 = 0 ⇔
�𝑘(𝑘 − 1)𝑐𝑘𝑥𝑘−2∞
𝑘=2
−� 2𝑘𝑐𝑘
∞
𝑘=1
𝑥𝑘 −� 4𝑐𝑘𝑥𝑘∞
𝑘=0
= 0 ⇔
4. Anpassar uttrycken. Vi måste ha: a) Samma exponent till 𝑥 b) Summor med samma intervall.
�𝑘(𝑘 − 1)𝑐𝑘𝑥𝑘−2∞
𝑘=2
= �(𝑘 + 2)(𝑘 + 1)𝑐𝑘+2𝑥𝑘∞
𝑘=0
� 2𝑘𝑐𝑘
∞
𝑘=1
𝑥𝑘 = � 2𝑘𝑐𝑘
∞
𝑘=0
𝑥𝑘
5. Förenklar uttrycket. �(𝑘 + 2)(𝑘 + 1)𝑐𝑘+2𝑥𝑘∞
𝑘=0
−� 2𝑘𝑐𝑘
∞
𝑘=0
𝑥𝑘 −� 4𝑐𝑘𝑥𝑘∞
𝑘=0
= 0 ⇔
��(𝑘 + 2)(𝑘 + 1)𝑐𝑘+2 − 2𝑘𝑐𝑘 − 4𝑐𝑘�𝑥𝑘∞
𝑘=0
= 0 ⇔
6. Bestämmer uttryck för 𝑐𝑘 Från uppgiften får vi begynnelsevillkor:
� 𝑦(0) = 0 ⇒ 𝑐0 = 0
𝑦′(0) = 1 ⇒ 𝑐1 = 1
Genom att jämföra med 𝐻𝐿 ser man att koefficienten till 𝑥𝑘 måste vara noll. (𝑘 + 2)(𝑘 + 1)𝑐𝑘+2 − 2𝑘𝑐𝑘 − 4𝑐𝑘 = 0 ⇔ (𝑘 + 2)(𝑘 + 1)𝑐𝑘+2 = 4𝑐𝑘 + 2𝑘𝑐𝑘 ⇔
𝑐𝑘+2 =4 + 2𝑘
(𝑘 + 2)(𝑘 + 1) ∙ 𝑐𝑘
Med hjälp av formeln ovan ska man nu komma fram till ett samband. Först ser man att 𝑐𝑘 = 0 för jämna 𝑘. Detta kan visas:
𝒌 𝑐𝑘 𝑐𝒌+2 =4 + 2𝒌
(𝒌 + 2)(𝒌 + 1) ∙ 𝑐𝑘
𝒌 = 𝟎 𝑐0 = 0 𝑐2 =
4 + 2𝒌(𝒌 + 2)(𝒌 + 1) ∙ 0 = 0
𝒌 = 𝟐 𝑐2 = 0 𝑐4 =4 + 2𝑘
(𝑘 + 2)(𝑘 + 1) ∙ 0 = 0
Sida 40 av 44
Osv… Återstår då att kontrollera udda 𝑘.
𝒌 𝑐𝑘 𝑐𝒌+2 =4 + 2𝒌
(𝒌 + 2)(𝒌 + 1) ∙ 𝑐𝑘
𝒌 = 𝟏 𝑐1 = 1 𝑐3 =
4 + 2 ∙ 𝟏(𝟏 + 2)(𝟏 + 1) ∙ 1 =
63 ∙ 2
= 1
𝒌 = 𝟑 𝑐3 = 1 𝑐5 =4 + 2 ∙ 𝟑
(𝟑 + 2)(𝟑 + 1) ∙ 1 =10
5 ∙ 4=
12
𝒌 = 𝟓 𝑐5 =12
𝑐7 =4 + 2 ∙ 𝟓
(𝟓 + 2)(𝟓 + 1) ∙12
=14
7 ∙ 6∙
12
=1
2 ∙ 3
𝒌 = 𝟕 𝑐7 =1
2 ∙ 3 𝑐9 =
4 + 2 ∙ 𝟕(𝟕 + 2)(𝟕 + 1) ∙
12 ∙ 3
=18
9 ∙ 8∙
12 ∙ 3
=1
2 ∙ 3 ∙ 4
Här ser man att vi verkar få fakultet i nämnaren för termerna. Nu gäller det bara att hitta ett samband mellan nämnaren och värdet på 𝑘. 12!