Loading [MathJax]/jax/element/mml/optable/BasicLatin.js
8 svar
162 visningar
Victordobado 24
Postad: 27 dec 2020 16:03

Differentialekvation

Jag ska lösa denna ekvation:

y''-3y+2y = e^x

yh kunde jag lösa, men yp blev en aning svårare. Jag skrev ansatsen yp = Ae^x, men får inget värde på A. Hur ska man tänka då?

Moffen 1877
Postad: 27 dec 2020 16:08

Hej!

Vad fick du för yh?

Victordobado 24
Postad: 27 dec 2020 16:11

yh fick jag C1e^x + C2e^2x, där C1 och C2 är godtyckliga konstanter. När jag försökte lösa yp så började med att derivera ansatsen och sätta in första och andraderivatan i ekvationen, men A värdet försvinner.

Skaft 2373 – F.d. Moderator
Postad: 27 dec 2020 16:14

Prova att göra koefficienten "mer generell" med ett förstagradspolynom istället för en konstant: yp=(ax+b)ex

Albiki 5096 – Fd. Medlem
Postad: 27 dec 2020 18:13

Hej,

Metoden Variation av koefficienter går ut på att skapa en ansats där konstanten A ersätts med en funktion A(x). Därefter identifierar du vilka funktioner A(x) som passar ihop med den givna differentialekvationen.

Här ansätter du alltså partikulärlösningen y(x)=A(x)ex som ger förstaderivatan

    y'

och andraderivatan

    y''(x)={A''(x)+2A'(x)+A(x)}exy^{''}(x) = \{A^{''}(x)+2A^\prime(x)+A(x)\}e^x.

Insättning i den givna differentialekvationen ger en motsvarande differentialekvation som funktionen A(x)A(x) måste uppfylla.

    A''x-A'x-Ax=0.\displaystyle A^{''}\left(x\right)-A^\prime\left(x\right)-A\left(x\right)=0.

Lösningar till denna ekvation är

    A(x)=C1eax+C2ebxA(x)=C_1e^{ax}+C_2e^{bx}

där a=1-52a=\frac{1-\sqrt{5}}{2} och b=1+52b=\frac{1+\sqrt{5}}{2}.

Victordobado 24
Postad: 27 dec 2020 18:33
Skaft skrev:

Prova att göra koefficienten "mer generell" med ett förstagradspolynom istället för en konstant: yp=(ax+b)exy_p = (ax+b)e^x

Jag testade göra det här, men tyckte att det blev bökigt i slutet för jag visste inte hur jag skulle avläsa mina a och b värden efter att jag bröt ut e^x från båda leden.

Skaft 2373 – F.d. Moderator
Postad: 27 dec 2020 19:01

Om du har yp=(ax+b)exy_p = (ax+b)e^x så kan derivatorna skrivas

yp'=aex+ypyp''=2aex+ypy_p^\prime = ae^x + y_p\\ y_p^{\prime\prime} = 2ae^x + y_p

Ekvationens vänsterled y''-3y'+2yy^{\prime\prime}-3y^\prime+2y blir då:

(2aex+yp)-3(aex+yp)+2yp(2ae^x + y_p)-3(ae^x + y_p)+2y_p

dvs bara -aex-ae^x, alla ypy_p tar ut varann. För att detta ska vara lika med högerledet exe^x måste alltså a=-1a=-1. bb kan däremot väljas godtyckligt, och när du lägger ihop yh+ypy_h+y_p kan bb "smälta ihop" med C1C_1 som också är en godtycklig koefficient till exe^x.

Victordobado 24
Postad: 27 dec 2020 19:46
Skaft skrev:

Om du har yp=(ax+b)exy_p = (ax+b)e^x så kan derivatorna skrivas

yp'=aex+ypyp''=2aex+ypy_p^\prime = ae^x + y_p\\ y_p^{\prime\prime} = 2ae^x + y_p

Ekvationens vänsterled y''-3y'+2yy^{\prime\prime}-3y^\prime+2y blir då:

(2aex+yp)-3(aex+yp)+2yp(2ae^x + y_p)-3(ae^x + y_p)+2y_p

dvs bara -aex-ae^x, alla ypy_p tar ut varann. För att detta ska vara lika med högerledet exe^x måste alltså a=-1a=-1. bb kan däremot väljas godtyckligt, och när du lägger ihop yh+ypy_h+y_p kan bb "smälta ihop" med C1C_1 som också är en godtycklig koefficient till exe^x.

Om yp = (ax+b)*e^x,  borde inte yp' = e^x(ax+b+a), eller är jag helt ute och cyklar?

Skaft 2373 – F.d. Moderator
Postad: 27 dec 2020 19:56

Jo, det stämmer. Men det kan du skriva om som:

yp'=(ax+b+a)ex=(ax+b)ex+aexy_p^\prime = (ax+b+a)e^x = (ax+b)e^x+ae^x

Och eftersom yp=(ax+b)exy_p = (ax+b)e^x, så får vi att

yp'=yp+aexy_p^\prime = y_p+ae^x

Att "återanvända" ypy_p för att beskriva derivatorna är inte nödvändigt, men det kan underlätta. Förenklar förenklandet, så att säga =)

Svara
Close