Диференцијалне једначине са константим коефицијентима

y"+py′+qy=f(x)Најпре решавамо хомогени део диференцијалне једначине : y"+py′+qy=0k2+pk+q=0 -карактеристична једначина k1/2 -корени карактеристичне једначине Ако је : k1≠k2   ⇒   yh=C1ek1x+C2ek2xk1=k2   ⇒   yh=C1ek1x+C2xek1xk1/2=α±βi   ⇒   yh=eαx(C1cos⁡βx+C2sin⁡βx)Ако је нехомогени део облика f(x)=eax(Pn(x)cos⁡bx+Qn(x)sin⁡bx)Тада је партикуларно решење облика : yp=xreax(SN(x)cos⁡bx+TN(x)sin⁡bx)где је r-вишеструкост a±bi као корена карактеристичне једначине, а N=max⁡{n,m}\begin{array}{l} \\ y" + py' + qy = f\left( x \right) \\ \\ {\text{Најпре решавамо хомогени део диференцијалне једначине : }} \\ y" + py' + qy = 0 \\ {k^2} + pk + q = 0 {\text{ -карактеристична једначина }} \\ {k_{1/2}} {\text{ -корени карактеристичне једначине }} \\ {\text{Ако је : }} \\ {k_1} \ne {k_2} {\text{ }}{\text{ }}{\text{ }} \Rightarrow {\text{ }}{\text{ }}{\text{ }} {y_h} = {C_1}{e^{{k_1}x}} + {C_2}{e^{{k_2}x}} \\ {k_1} = {k_2} {\text{ }}{\text{ }}{\text{ }} \Rightarrow {\text{ }}{\text{ }}{\text{ }} {y_h} = {C_1}{e^{{k_1}x}} + {C_2}x{e^{{k_1}x}} \\ {k_{1/2}} = \alpha \pm \beta i {\text{ }}{\text{ }}{\text{ }} \Rightarrow {\text{ }}{\text{ }}{\text{ }} {y_h} = {e^{\alpha x}}\left( {{C_1}\cos \beta x + {C_2}\sin \beta x} \right) \\ \\ {\text{Ако је нехомогени део облика }} \\ f\left( x \right) = {e^{ax}}\left( {{P_n}\left( x \right)\cos bx + {Q_n}\left( x \right)\sin bx} \right) \\ {\text{Тада је партикуларно решење облика : }} \\ {y_p} = {x^r}{e^{ax}}\left( {{S_N}\left( x \right)\cos bx + {T_N}\left( x \right)\sin bx} \right) \\ {\text{где је r-вишеструкост $a \pm bi$ као корена карактеристичне једначине, а $N = \max \left\{ {n,m} \right\}$}} \\ \end{array}

Задаци


1. Наћи опште решење диференцијалне једначине   y"+5y′+6y=xy" + 5y' + 6y = x.

Најпре решавамо хомогени део ове диференцијалне једначине

y"+5y′+6y=0y" + 5y' + 6y = 0

Решавајући њену карактеристичну једначину

k2+5k+6=0{k^2} + 5k + 6 = 0

добијамо да су њена решења

k1=−3  ∧  k2=−2{k_1} = - 3{\rm{ }} {\text{ }}{\text{ }} \wedge {\rm{ }}{\text{ }}{\text{ }} {k_2} = - 2

па је хомогено решење наше диференцијалне једначине

yh=C1e−3x+C2e−2x{y_h} = {C_1}{e^{ - 3x}} + {C_2}{e^{ - 2x}}

Да бисмо одредили партикуларно решење посматрамо нехомогени део

f(x)=xf\left( x \right) = x

Како нехомогени део треба да буде облика f(x)=eax(Pn(x)cos⁡bx+Qm(x)sin⁡bx)f\left( x \right) = {e^{ax}}\left( {{P_n}\left( x \right)\cos bx + {Q_m}\left( x \right)\sin bx} \right)

Закључујемо

a=0b=0}⇒  a±ib=0  ⇒  r=0\left. \begin{array}{l} a = 0\\ b = 0 \end{array} \right\} \Rightarrow {\text{ }}{\text{ }} a \pm ib = 0 {\text{ }}{\text{ }} \Rightarrow{\text{ }}{\text{ }} r = 0

n=1m=/}⇒  N=1\left. \begin{array}{l} n = 1\\ m = / \end{array} \right\} \Rightarrow{\text{ }}{\text{ }} N = 1

Па је партикуларно решење облика

yp=x0e0x((Ax+B)cos⁡0x+(Cx+D)sin⁡0x){y_p} = {x^0}{e^{0x}}\left( {\left( {Ax + B} \right)\cos 0x + \left( {Cx + D} \right)\sin 0x} \right)

односно yp=Ax+B{y_p} = Ax + B

Потребно је још одредити коефицијенте А и В.

То ћемо учинити убацујући партикуларно решење у почетну једначину.

yp′=A{y_p}' = A

yp"=0{y_p}" = 0

0+5A+6(Ax+B)=x0 + 5A + 6\left( {Ax + B} \right) = x

6Ax+5A+6B=x6Ax + 5A + 6B = x

Изједначавајући коефицијенте уз одговарајуће степене добијамо

6A=1  ⇒  A=166A = 1 {\text{ }}{\text{ }} \Rightarrow{\text{ }}{\text{ }} A = \frac{1}{6}

5A+6B=0  ⇒  B=−5365A + 6B = 0 {\text{ }}{\text{ }} \Rightarrow{\text{ }}{\text{ }} B = - \frac{5}{{36}}

па је партикуларно решење

yp=16x−536{y_p} = \frac{1}{6}x - \frac{5}{{36}}

односно коначно решење полазне диференцијалне једначине

Y=C1e−3x+C2e−2x+16x−536Y = {C_1}{e^{ - 3x}} + {C_2}{e^{ - 2x}} + \frac{1}{6}x - \frac{5}{{36}}


2. Наћи опште решење диференцијалне једначине   y"′−4y′=x2e2xy"' - 4y' = {x^2}{e^{2x}}

Најпре решавамо хомогени део ове диференцијалне једначине

y"′−4y′=0y"' - 4y' = 0

k3−4k=0{k^3} - 4k = 0

k(k2−4)=0k\left( {{k^2} - 4} \right) = 0

k1=0  ∧  k2=2  ∧  k3=−2{k_1} = 0{\rm{ }}{\text{ }}{\text{ }} \wedge {\rm{ }}{\text{ }}{\text{ }} {k_2} = 2{\rm{ }} {\text{ }}{\text{ }} \wedge {\rm{ }}{\text{ }}{\text{ }} {k_3} = - 2

Следи да је хомогено решење наше диференцијалне једначине

yh=C1e0x+C2e2x+C3e−2x{y_h} = {C_1}{e^{0x}} + {C_2}{e^{2x}} + {C_3}{e^{ - 2x}}

yh=C1+C2e2x+C3e−2x{y_h} = {C_1} + {C_2}{e^{2x}} + {C_3}{e^{ - 2x}}

Из нехомогеног дела f(x)=x2e2xf\left( x \right) = {x^2}{e^{2x}} закључујемо да је

a=2b=0}⇒  a±ib=2  ⇒  r=1\left. \begin{array}{l} a = 2\\ b = 0 \end{array} \right\} \Rightarrow {\rm{ }}{\text{ }}{\text{ }} a \pm ib = 2{\rm{ }}{\text{ }}{\text{ }} \Rightarrow {\rm{ }}{\text{ }}{\text{ }} r = 1

n=2m=0}⇒  N=2\left. \begin{array}{l} n = 2\\ m = 0 \end{array} \right\} \Rightarrow {\rm{ }}{\text{ }}{\text{ }} N = 2

Па је партикуларно решење облика

yp=x1e2x((Ax2+Bx+C)cos⁡0x+(Dx2+Ex+F)sin⁡0x){y_p} = {x^1}{e^{2x}}\left( {\left( {A{x^2} + Bx + C} \right)\cos 0x + \left( {D{x^2} + Ex + F} \right)\sin 0x} \right)

yp=xe2x(Ax2+Bx+C){y_p} = x{e^{2x}}\left( {A{x^2} + Bx + C} \right)

Одредимо још коефицијенте А,B,C и D

yp′=2e2x(Ax3+Bx2+Cx)+e2x(3Ax2+2Bx+C){y_p}' = 2{e^{2x}}\left( {A{x^3} + B{x^2} + Cx} \right) + {e^{2x}}\left( {3A{x^2} + 2Bx + C} \right)

yp′=e2x(2Ax3+(2B+3A)x2+(2C+2B)x+C){y_p}' = {e^{2x}}\left( {2A{x^3} + \left( {2B + 3A} \right){x^2} + \left( {2C + 2B} \right)x + C} \right)

yp"=2e2x(2Ax3+(2B+3A)x2+(2C+2B)x+C){y_p}" = 2{e^{2x}}\left( {2A{x^3} + \left( {2B + 3A} \right){x^2} + \left( {2C + 2B} \right)x + C} \right)

      +e2x(6Ax2+(4B+6A)x+2C+2B){\text{ }}{\text{ }}{\text{ }} {\text{ }} {\text{ }} {\text{ }} + {e^{2x}}\left( {6A{x^2} + \left( {4B + 6A} \right)x + 2C + 2B} \right)

yp"=e2x(4Ax3+(4B+12A)x2+(4C+8B+6A)x+4C+2B){y_p}" = {e^{2x}}\left( {4A{x^3} + \left( {4B + 12A} \right){x^2} + \left( {4C + 8B + 6A} \right)x + 4C + 2B} \right)

yp"′=2e2x(4Ax3+(4B+12A)x2+(4C+8B+6A)x+4C+2B){y_p}"' = 2{e^{2x}}\left( {4A{x^3} + \left( {4B + 12A} \right){x^2} + \left( {4C + 8B + 6A} \right)x + 4C + 2B} \right)

      +e2x(12Ax2+(8B+24A)x+4C+8B+6A){\text{ }}{\text{ }}{\text{ }}{\text{ }} {\text{ }}{\text{ }} + {e^{2x}}\left( {12A{x^2} + \left( {8B + 24A} \right)x + 4C + 8B + 6A} \right)

yp"′=e2x(8Ax3+(8B+36A)x2+(8C+24B+36A)x+12C+12B+6A){y_p}"' = {e^{2x}}\left( {8A{x^3} + \left( {8B + 36A} \right){x^2} + \left( {8C + 24B + 36A} \right)x + 12C + 12B + 6A} \right)

Уврштавајући партикуларно решење у полазну једначину

y"′−4y′=x2e2xy"' - 4y' = {x^2}{e^{2x}} добијамо

e2x(8Ax3+(8B+36A)x2+(8C+24B+36A)x+12C+12B+6A){e^{2x}}\left( {8A{x^3} + \left( {8B + 36A} \right){x^2} + \left( {8C + 24B + 36A} \right)x + 12C + 12B + 6A} \right)

−4e2x(2Ax3+(2B+3A)x2+(2C+2B)x+C)=x2e2x- 4{e^{2x}}\left( {2A{x^3} + \left( {2B + 3A} \right){x^2} + \left( {2C + 2B} \right)x + C} \right)= {x^2}{e^{2x}}

(8A−8A)x3+(8B+36A−8B−12A)x2+(8C+24B+36A−8C−8B)x\left( {8A - 8A} \right){x^3} + \left( {8B + 36A - 8B - 12A} \right){x^2} + \left( {8C + 24B + 36A - 8C - 8B} \right)x

+(12C+12B+6A−4C)=x2+ \left( {12C + 12B + 6A - 4C} \right) = {x^2}

Изједначавајући коефицијенте уз одговарајуће степене добијамо

24A=1  ⇒  A=12424A = 1{\text{ }}{\text{ }} \Rightarrow{\text{ }}{\text{ }} A = \frac{1}{{24}}

36A+16B=0  ⇒  B=−33236A + 16B = 0 {\text{ }}{\text{ }} \Rightarrow {\text{ }}{\text{ }} B = - \frac{3}{{32}}

6A+12B+8C=0  ⇒  C=7646A + 12B + 8C = 0 {\text{ }}{\text{ }} \Rightarrow {\text{ }}{\text{ }} C = \frac{7}{{64}}

па је партикуларно решење

yp=xe2x(124x2−332x+764){y_p} = x{e^{2x}}\left( {\frac{1}{{24}}{x^2} - \frac{3}{{32}}x + \frac{7}{{64}}} \right)

односно коначно решење полазне диференцијалне једначине

Y=yh+yp=C1+C2e2x+C3e−2x+e2x(124x3−332x2+764x)Y = {y_h} + {y_p} = {C_1} + {C_2}{e^{2x}} + {C_3}{e^{ - 2x}} + {e^{2x}}\left( {\frac{1}{{24}}{x^3} - \frac{3}{{32}}{x^2} + \frac{7}{{64}}x} \right)

3. Наћи опште решење диференцијалне једначине   y"′−2y"=xsin⁡2xy"' - 2y" = x\sin 2x

Најпре решавамо хомогени део ове диференцијалне једначине

y"′−2y"=0y"' - 2y" = 0

k3−2k2=0{k^3} - 2{k^2} = 0

k2(k−2)=0{k^2}\left( {k - 2} \right) = 0

k1,2=0  ∧  k3=2{k_{1,2}} = 0{\rm{ }}{\text{ }}{\text{ }} \wedge{\text{ }}{\text{ }} {\rm{ }}{k_3} = 2

Следи да је хомогено решење наше диференцијалне једначине

yh=C1e0x+C2xe0x+C3e2x{y_h} = {C_1}{e^{0x}} + {C_2}x{e^{0x}} + {C_3}{e^{2x}}

yh=C1+C2x+C3e2x{y_h} = {C_1} + {C_2}x + {C_3}{e^{2x}}

Из нехомогеног дела f(x)=xsin⁡2xf\left( x \right) = x\sin 2x закључујемо да је

a=0b=2}⇒  a±bi=2i  ⇒  r=0\left. \begin{array}{l} a = 0\\ b = 2 \end{array} \right\} \Rightarrow {\text{ }}{\text{ }} a \pm bi = 2i {\text{ }}{\text{ }} \Rightarrow{\text{ }}{\text{ }} r = 0

n=0m=1}⇒  N=1\left. \begin{array}{l} n = 0\\ m = 1 \end{array} \right\} \Rightarrow {\text{ }}{\text{ }} N = 1

Па је партикуларно решење облика

yp=x0e0x((Ax+B)cos⁡2x+(Cx+D)sin⁡2x){y_p} = {x^0}{e^{0x}}\left( {\left( {Ax + B} \right)\cos 2x + \left( {Cx + D} \right)\sin 2x} \right)

yp=(Ax+B)cos⁡2x+(Cx+D)sin⁡2x{y_p} = \left( {Ax + B} \right)\cos 2x + \left( {Cx + D} \right)\sin 2x

Одредимо још коефицијенте А,B,C и D

yp′=Acos⁡2x+(Ax+B)(−2sin⁡2x)+Csin⁡2x+(Cx+D)2cos⁡2x{y_p}' = A\cos 2x + \left( {Ax + B} \right)\left( { - 2\sin 2x} \right) + C\sin 2x + \left( {Cx + D} \right)2\cos 2x

yp′=(A+2Cx+2D)cos⁡2x+(−2Ax−2B+C)sin⁡2x{y_p}' = \left( {A + 2Cx + 2D} \right)\cos 2x + \left( { - 2Ax - 2B + C} \right)\sin 2x

yp"=2Ccos⁡2x+(A+2Cx+2D)(−2sin⁡2x)−2Asin⁡2x+(−2Ax−2B+C)2cos⁡2x{y_p}" = 2C\cos 2x + \left( {A + 2Cx + 2D} \right)\left( { - 2\sin 2x} \right) - 2A\sin 2x + \left( { - 2Ax - 2B + C} \right)2\cos 2x

yp"=(4C−4Ax−4B)cos⁡2x+(−4A−4Cx−4D)sin⁡2x{y_p}" = \left( {4C - 4Ax - 4B} \right)\cos 2x + \left( { - 4A - 4Cx - 4D} \right)\sin 2x

yp"′=−4Acos⁡2x+(4C−4Ax−4B)(−2sin⁡2x)−4Csin⁡2x+(−4A−4Cx−4D)2cos⁡2x{y_p}"' = - 4A\cos 2x + \left( {4C - 4Ax - 4B} \right)\left( { - 2\sin 2x} \right) - 4C\sin 2x + \left( { - 4A - 4Cx - 4D} \right)2\cos 2x

yp"′=(−12A−8Cx−8D)cos⁡2x+(−12C+8Ax+8B)sin⁡2x{y_p}"' = \left( { - 12A - 8Cx - 8D} \right)\cos 2x + \left( { - 12C + 8Ax + 8B} \right)\sin 2x

Уврштавајући партикуларно решење у полазну једначину

y"′−2y"=xsin⁡2xy"' - 2y" = x\sin 2x

(−12A−8Cx−8D)cos⁡2x+(−12C+8Ax+8B)sin⁡2x+\left( { - 12A - 8Cx - 8D} \right)\cos 2x + \left( { - 12C + 8Ax + 8B} \right)\sin 2x +

+(−8C+8Ax+8B)cos⁡2x+(8A+8Cx+8D)sin⁡2x=xsin⁡2x+ \left( { - 8C + 8Ax + 8B} \right)\cos 2x + \left( {8A + 8Cx + 8D} \right)\sin 2x = x\sin 2x

(−12A−8Cx−8D−8C+8Ax+8B)cos⁡2x+\left( { - 12A - 8Cx - 8D - 8C + 8Ax + 8B} \right)\cos 2x +

+(−12C+8Ax+8B+8A+8Cx+8D)sin⁡2x=xsin⁡x+ \left( { - 12C + 8Ax + 8B + 8A + 8Cx + 8D} \right)\sin 2x = x\sin x

−12A−8Cx−8D−8C+8Ax+8B=0- 12A - 8Cx - 8D - 8C + 8Ax + 8B = 0

−12C+8Ax+8B+8A+8Cx+8D=x- 12C + 8Ax + 8B + 8A + 8Cx + 8D = x

−12A−8D−8C+8B=0- 12A - 8D - 8C + 8B = 0

−8C+8A=0- 8C + 8A = 0

−12C+8B+8A+8D=0- 12C + 8B + 8A + 8D = 0

8A+8C=18A + 8C = 1

A=116A = \frac{1}{{16}}

C=116C = \frac{1}{{16}}

B=332B = \frac{3}{{32}}

D=−116D = - \frac{1}{{16}}

Па је партикуларно решење наше једначине

yp=(116x+332)cos⁡2x+(116x−116)sin⁡2x{y_p} = \left( {\frac{1}{{16}}x + \frac{3}{{32}}} \right)\cos 2x + \left( {\frac{1}{{16}}x - \frac{1}{{16}}} \right)\sin 2x

Односно коначно решење

Y=C1+C2x+C3e2x+(116x+332)cos⁡2x+(116x−116)sin⁡2xY = {C_1} + {C_2}x + {C_3}{e^{2x}} + \left( {\frac{1}{{16}}x + \frac{3}{{32}}} \right)\cos 2x + \left( {\frac{1}{{16}}x - \frac{1}{{16}}} \right)\sin 2x

4. Диференцијалну једначину   y"+y′tgx−ycos⁡2x=e2sin⁡xcos⁡2xy" + y'tgx - y{\cos ^2}x = {e^{2\sin x}}{\cos ^2}x   сменом   t=sin⁡xt = \sin x   свести на диференцијалну једначину са константним коефицијентима и решити је.

t=sin⁡xt = \sin x

cos⁡x=1−sin⁡2x=1−t2\cos x = \sqrt {1 - {{\sin }^2}x} = \sqrt {1 - {t^2}}

tgx=sin⁡xcos⁡x=t1−t2tgx = \frac{{\sin x}}{{\cos x}} = \frac{t}{{\sqrt {1 - {t^2}} }}

t=sin⁡xt = \sin x

dt=cos⁡xdxdt = \cos xdx

dtdx=cos⁡x=1−t2\frac{{dt}}{{dx}} = \cos x = \sqrt {1 - {t^2}}

y′=dydx=dydt⋅dtdx=y˙⋅cos⁡x=y˙⋅1−t2y' = \frac{{dy}}{{dx}} = \frac{{dy}}{{dt}} \cdot \frac{{dt}}{{dx}} = \dot y \cdot \cos x= \dot y \cdot \sqrt {1 - {t^2}}

y"=dy′dx=dy′dt⋅dtdx=ddt(y˙1−t2)⋅cos⁡x=y" = \frac{{dy'}}{{dx}} = \frac{{dy'}}{{dt}} \cdot \frac{{dt}}{{dx}} = \frac{d}{{dt}}\left( {\dot y\sqrt {1 - {t^2}} } \right) \cdot \cos x =

    =(y¨1−t2+y˙121−t2⋅(−2t))⋅1−t2={\text{ }}{\text{ }} {\text{ }}{\text{ }} = \left( {\ddot y\sqrt {1 - {t^2}} + \dot y\frac{1}{{2\sqrt {1 - {t^2}} }} \cdot \left( { - 2t} \right)} \right) \cdot \sqrt {1 - {t^2}} =

    =y¨(1−t2)−ty˙{\text{ }}{\text{ }} {\text{ }}{\text{ }} = \ddot y\left( {1 - {t^2}} \right) - t\dot y

Уврштавајући у почетну диференцијалну једначину добијамо:

y¨(1−t2)−ty˙+y˙1−t2⋅t1−t2−y(1−t2)=e2t(1−t2)\ddot y\left( {1 - {t^2}} \right) - t\dot y + \dot y\sqrt {1 - {t^2}} \cdot \frac{t}{{\sqrt {1 - {t^2}} }} - y\left( {1 - {t^2}} \right) = {e^{2t}}\left( {1 - {t^2}} \right)

y¨(1−t2)−ty˙+y˙t−y(1−t2)=e2t(1−t2)\ddot y\left( {1 - {t^2}} \right) - t\dot y + \dot yt - y\left( {1 - {t^2}} \right) = {e^{2t}}\left( {1 - {t^2}} \right)

y¨(1−t2)−y(1−t2)=e2t(1−t2)  /:(1−t2)\ddot y\left( {1 - {t^2}} \right) - y\left( {1 - {t^2}} \right) = {e^{2t}}\left( {1 - {t^2}} \right){\text{ }}{\text{ }} /:\left( {1 - {t^2}} \right)

y¨−y=e2t\ddot y - y = {e^{2t}}

Овим смо нашу диференцијалну једначину свели на диференцијалну једначину са константним коефицијентима, коју сада лако решавамо.

y¨−y=0\ddot y - y = 0

k2−1=0{k^2} - 1 = 0

k2=1{k^2} = 1

k1/2=±1{k_{1/2}} = \pm 1

yh=C1e1t+C2e−1t=C1esin⁡x+C2e−sin⁡x{y_h} = {C_1}{e^{1t}} + {C_2}{e^{ - 1t}} = {C_1}{e^{\sin x}} + {C_2}{e^{ - \sin x}}

f(t)=e2tf\left( t \right) = {e^{2t}}

a=2b=0}⇒  a±ib=2  ⇒  r=0\left. \begin{array}{l} a = 2\\ b = 0 \end{array} \right\} \Rightarrow {\text{ }}{\text{ }} a \pm ib = 2 {\text{ }}{\text{ }} \Rightarrow{\text{ }}{\text{ }} r = 0

n=0m=0}⇒  N=0\left. \begin{array}{l} n = 0\\ m = 0 \end{array} \right\} \Rightarrow{\text{ }}{\text{ }} N = 0

yp=Ae2t{y_p} = A{e^{2t}}

yp′=2Ae2t{y_p}' = 2A{e^{2t}}

yp"=4Ae2t{y_p}" = 4A{e^{2t}}

4Ae2t−Ae2t=e2t4A{e^{2t}} - A{e^{2t}} = {e^{2t}}

3Ae2t=e2t3A{e^{2t}} = {e^{2t}}

A=13A = \frac{1}{3}

yp=13e2t=e2sin⁡x3{y_p} = \frac{1}{3}{e^{2t}} = \frac{{{e^{2\sin x}}}}{3}

Y=yh+yp=C1esin⁡x+C2e−sin⁡x+e2sin⁡x3Y = {y_h} + {y_p} = {C_1}{e^{\sin x}} + {C_2}{e^{ - \sin x}} + \frac{{{e^{2\sin x}}}}{3}

5. Диференцијалну једначину   (1+ex)y"−(1+3ex+e2x)y′=(1+ex)3(x+ex+exeex)\left( {1 + {e^x}} \right)y" - \left( {1 + 3{e^x} + {e^{2x}}} \right)y' = {\left( {1 + {e^x}} \right)^3}\left( {x + {e^x} + {e^x}{e^{{e^x}}}} \right)   сменом   t=x+ext = x + {e^x}   свести на диференцијалну једначину са константним коефицијентима и решити је.

t=x+ext = x + {e^x}

dt=(1+ex)dxdt = \left( {1 + {e^x}} \right)dx

dtdx=1+ex\frac{{dt}}{{dx}} = 1 + {e^x}

y′=dydx=dydt⋅dtdx=y˙(1+ex)y' = \frac{{dy}}{{dx}} = \frac{{dy}}{{dt}} \cdot \frac{{dt}}{{dx}} = \dot y\left( {1 + {e^x}} \right)

y"=dy′dx=ddx(y˙(1+ex))=dy˙dx(1+ex)+y˙d(1+ex)dx=y" = \frac{{dy'}}{{dx}} = \frac{d}{{dx}}\left( {\dot y\left( {1 + {e^x}} \right)} \right) = \frac{{d\dot y}}{{dx}}\left( {1 + {e^x}} \right) + \dot y\frac{{d\left( {1 + {e^x}} \right)}}{{dx}} =

    =dy˙dtdtdx(1+ex)+y˙ex=y¨(1+ex)(1+ex)+y˙ex{\text{ }}{\text{ }} {\text{ }}{\text{ }} = \frac{{d\dot y}}{{dt}} \frac{{dt}}{{dx}}\left( {1 + {e^x}} \right) + \dot y{e^x} = \ddot y \left( {1 + {e^x}} \right)\left( {1 + {e^x}} \right) + \dot y{e^x}

(1+ex)(y¨(1+ex)2+y˙ex)−(1+3ex+e2x)y˙(1+ex)=(1+ex)3(x+ex+ex+ex)\left( {1 + {e^x}} \right)\left( {\ddot y {{\left( {1 + {e^x}} \right)}^2} + \dot y{e^x}} \right) - \left( {1 + 3{e^x} + {e^{2x}}} \right)\dot y\left( {1 + {e^x}} \right) = {\left( {1 + {e^x}} \right)^3}\left( {x + {e^x} + {e^{x + {e^x}}}} \right)

y¨(1+ex)3+y˙ex(1+ex)−(1+3ex+e2x)y˙(1+ex)=(1+ex)3(t+et)\ddot y {\left( {1 + {e^x}} \right)^3} + \dot y{e^x}\left( {1 + {e^x}} \right) - \left( {1 + 3{e^x} + {e^{2x}}} \right)\dot y\left( {1 + {e^x}} \right) = {\left( {1 + {e^x}} \right)^3}\left( {t + {e^t}} \right)

y¨(1+ex)3+y˙(1+ex)(ex−(1+3ex+e2x))=(1+ex)3(t+et)\ddot y {\left( {1 + {e^x}} \right)^3} + \dot y\left( {1 + {e^x}} \right)\left( {{e^x} - \left( {1 + 3{e^x} + {e^{2x}}} \right)} \right) = {\left( {1 + {e^x}} \right)^3}\left( {t + {e^t}} \right)

y¨(1+ex)3+y˙(1+ex)(−1−2ex−e2x)=(1+ex)3(t+et)\ddot y {\left( {1 + {e^x}} \right)^3} + \dot y\left( {1 + {e^x}} \right)\left( { - 1 - 2{e^x} - {e^{2x}}} \right) = {\left( {1 + {e^x}} \right)^3}\left( {t + {e^t}} \right)

y¨(1+ex)3+y˙(1+ex)(−(1+ex)2)=(1+ex)3(t+et)\ddot y {\left( {1 + {e^x}} \right)^3} + \dot y\left( {1 + {e^x}} \right)\left( { - {{\left( {1 + {e^x}} \right)}^2}} \right) = {\left( {1 + {e^x}} \right)^3}\left( {t + {e^t}} \right)

y¨(1+ex)3−y˙(1+ex)3=(1+ex)3(t+et)  /:(1+ex)3\ddot y {\left( {1 + {e^x}} \right)^3} - \dot y{\left( {1 + {e^x}} \right)^3} = {\left( {1 + {e^x}} \right)^3}\left( {t + {e^t}} \right){\rm{ }}{\text{ }}{\text{ }} /:{\left( {1 + {e^x}} \right)^3}

y¨−y˙=t+et\ddot y - \dot y = t + {e^t}

Најпре решавамо хомогени део ове диференцијалне једначине.

y¨−y˙=0\ddot y - \dot y = 0

k2−k=0{k^2} - k = 0

k(k−1)=0k\left( {k - 1} \right) = 0

k1=0  ∧  k2=1{k_1} = 0{\rm{ }}{\text{ }}{\text{ }} \wedge {\rm{ }} {\text{ }}{\text{ }} {k_2} = 1

yh=C1e0t+C2e1t{y_h} = {C_1}{e^{0t}} + {C_2}{e^{1t}}

yh=C1+C2et{y_h} = {C_1} + {C_2}{e^t}

yh=C1+C2ex+ex{y_h} = {C_1} + {C_2}{e^{x + {e^x}}}

Нехомогени део раздвојићемо на две функције и на тај начин добити два партикуларна решења.

f1(t)=t{f_1}\left( t \right) = t

a=0b=0}⇒  a±ib=0  ⇒  r=1\left. \begin{array}{l} a = 0\\ b = 0 \end{array} \right\} \Rightarrow {\rm{ }}{\text{ }}{\text{ }} a \pm ib = 0{\rm{ }}{\text{ }}{\text{ }} \Rightarrow {\rm{ }}{\text{ }}{\text{ }} r = 1

n=1m=/}⇒  N=1\left. \begin{array}{l} n = 1\\ m = / \end{array} \right\} \Rightarrow {\rm{ }}{\text{ }}{\text{ }} N = 1

yp1=t1e0t((At+B)cos⁡0t+(Ct+D)sin⁡0t){y_{{p_1}}} = {t^1}{e^{0t}}\left( {\left( {At + B} \right)\cos 0t + \left( {Ct + D} \right)\sin 0t} \right)

yp1=t(At+B){y_{{p_1}}} = t\left( {At + B} \right)

y˙p1=2At+B{{\dot y}_{{p_1}}} = 2At + B

y¨p1=2A{{\ddot y}_{{p_1}}} = 2A

2A−(2At+B)=t2A - \left( {2At + B} \right) = t

−2A=1  ⇒  A=−12- 2A = 1 {\text{ }}{\text{ }} \Rightarrow {\text{ }}{\text{ }} A = - \frac{1}{2}

2A+B=0  ⇒  B=12A + B = 0 {\text{ }}{\text{ }} \Rightarrow{\text{ }}{\text{ }} B = 1

yp1=−12t2+t{y_{{p_1}}} = - \frac{1}{2}{t^2} + t

yp1=−(x+ex)22+x+ex{y_{{p_1}}} = - \frac{{{{\left( {x + {e^x}} \right)}^2}}}{2} + x + {e^x}

f2(t)=et{f_2}\left( t \right) = {e^t}

a=1b=0}⇒  a±ib=1  ⇒  r=1\left. \begin{array}{l} a = 1\\ b = 0 \end{array} \right\} \Rightarrow {\rm{ }{\text{ }}{\text{ }} }a \pm ib = 1{\rm{ }}{\text{ }}{\text{ }} \Rightarrow {\rm{ }}{\text{ }}{\text{ }} r = 1

n=0m=0}⇒  N=0\left. \begin{array}{l} n = 0\\ m = 0 \end{array} \right\} \Rightarrow {\rm{ }} {\text{ }}{\text{ }} N = 0

yp2=t1e1t(Acos⁡0t+Bsin⁡0t){y_{{p_2}}} = {t^1}{e^{1t}}\left( {A\cos 0t + B\sin 0t} \right)

yp2=Atet{y_{{p_2}}} = At{e^t}

y˙p2=Aet+Atet=et(A+At){{\dot y}_{{p_2}}} = A{e^t} + At{e^t} = {e^t}\left( {A + At} \right)

y¨p2=et(A+At)+etA=et(2A+At){{\ddot y}_{{p_2}}} = {e^t}\left( {A + At} \right) + {e^t}A = {e^t}\left( {2A + At} \right)

et(2A+At)−et(A+At)=et{e^t}\left( {2A + At} \right) - {e^t}\left( {A + At} \right) = {e^t}

et(2A+At−A−At)=et{e^t}\left( {2A + At - A - At} \right) = {e^t}

A=1A = 1

yp2=tet{y_{{p_2}}} = t{e^t}

yp2=(x+ex)ex+ex{y_{{p_2}}} = \left( {x + {e^x}} \right){e^{x + {e^x}}}

Y=yh+yp1+yp2=C1+C2ex+ex−(x+ex)22+x+ex+(x+ex)ex+exY = {y_h} + {y_{{p_1}}} + {y_{{p_2}}} = {C_1} + {C_2}{e^{x + {e^x}}} - \frac{{{{\left( {x + {e^x}} \right)}^2}}}{2} + x + {e^x} + \left( {x + {e^x}} \right){e^{x + {e^x}}}

6. Сменом   t=arctgxt = arctgx   диференцијалну једначину   (1+x2)2y"+2x(1+x2)y′+y=4arctgx{\left( {1 + {x^2}} \right)^2}y" + 2x\left( {1 + {x^2}} \right)y' + y = 4arctgx   свести на диференцијалну једначину са константним коефицијентима и решити је.

t=arctgxt = arctgx

dt=11+x2dxdt = \frac{1}{{1 + {x^2}}}dx

dtdx=11+x2\frac{{dt}}{{dx}} = \frac{1}{{1 + {x^2}}}

y′=dydx=dydt⋅dtdx=y˙11+x2y' = \frac{{dy}}{{dx}} = \frac{{dy}}{{dt}} \cdot \frac{{dt}}{{dx}} = \dot y\frac{1}{{1 + {x^2}}}

y"=dy′dx=ddx(y˙11+x2)=dy˙dx11+x2+y˙d(11+x2)dx=y" = \frac{{dy'}}{{dx}} = \frac{d}{{dx}}\left( {\dot y\frac{1}{{1 + {x^2}}}} \right) = \frac{{d\dot y}}{{dx}}\frac{1}{{1 + {x^2}}} + \dot y\frac{{d\left( {\frac{1}{{1 + {x^2}}}} \right)}}{{dx}} =

    =dy˙dt⋅dtdx11+x2+y˙(−1(1+x2)22x)=y¨11+x211+x2−y˙2x(1+x2)2{\text{ }}{\text{ }} {\text{ }}{\text{ }} = \frac{{d\dot y}}{{dt}} \cdot \frac{{dt}}{{dx}}\frac{1}{{1 + {x^2}}} + \dot y\left( { - \frac{1}{{{{\left( {1 + {x^2}} \right)}^2}}}2x} \right) = \ddot y\frac{1}{{1 + {x^2}}}\frac{1}{{1 + {x^2}}} - \dot y\frac{{2x}}{{{{\left( {1 + {x^2}} \right)}^2}}}

(1+x2)2(y¨1(1+x2)2−y˙2x(1+x2)2)+2x(1+x2)y˙11+x2+y=4t{\left( {1 + {x^2}} \right)^2}\left( {\ddot y\frac{1}{{{{\left( {1 + {x^2}} \right)}^2}}} - \dot y\frac{{2x}}{{{{\left( {1 + {x^2}} \right)}^2}}}} \right) + 2x\left( {1 + {x^2}} \right)\dot y\frac{1}{{1 + {x^2}}} + y = 4t

y¨−2xy˙+2xy˙+y=4t\ddot y - 2x\dot y + 2x\dot y + y = 4t

y¨+y=4t\ddot y + y = 4t

y¨+y=0\ddot y + y = 0

k2+1=0{k^2} + 1 = 0

k1/2=±i{k_{1/2}} = \pm i

yh=e0t(C1cos⁡1t+C2sin⁡1t){y_h} = {e^{0t}}\left( {{C_1}\cos 1t + {C_2}\sin 1t} \right)

yh=C1cos⁡(arctgx)+C2sin⁡(arctgx){y_h} = {C_1}\cos \left( {arctgx} \right) + {C_2}\sin \left( {arctgx} \right)

yh=C111+x2+C2x1+x2{y_h} = {C_1}\frac{1}{{\sqrt {1 + {x^2}} }} + {C_2}\frac{x}{{\sqrt {1 + {x^2}} }}

f(t)=4tf\left( t \right) = 4t

a=0b=0}⇒  a±ib=0  ⇒  r=0\left. \begin{array}{l} a = 0\\ b = 0 \end{array} \right\} \Rightarrow {\rm{ }} {\text{ }}{\text{ }} a \pm ib = 0{\rm{ }} {\text{ }}{\text{ }} \Rightarrow {\rm{ }} {\text{ }}{\text{ }} r = 0

m=/n=1}⇒  N=1\left. \begin{array}{l} m = /\\ n = 1 \end{array} \right\} \Rightarrow {\rm{ }} {\text{ }}{\text{ }} N = 1

yp=t0e0t((At+B)cos⁡0t+(Ct+D)sin⁡0t){y_p} = {t^0}{e^{0t}}\left( {\left( {At + B} \right)\cos 0t + \left( {Ct + D} \right)\sin 0t} \right)

yp=At+B{y_p} = At + B

y˙p=A{{\dot y}_p} = A

y¨p=0{{\ddot y}_p} = 0

0+At+B=4t0 + At + B = 4t

A=4  ∧  B=0A = 4{\rm{ }} {\text{ }}{\text{ }} \wedge {\rm{ }} {\text{ }}{\text{ }} B = 0

yp=4t{y_p} = 4t

yp=4arctgx{y_p} = 4arctgx

Y=yh+yp=C111+x2+C2x1+x2+4arctgxY = {y_h} + {y_p} = {C_1}\frac{1}{{\sqrt {1 + {x^2}} }} + {C_2}\frac{x}{{\sqrt {1 + {x^2}} }} + 4arctgx

Једначина y′′+py′+qy=f(x)y'' + py' + qy = f(x) решава се у два корака: y=yh+ypy = y_h + y_p. Хомогени део решавате преко карактеристичне једначине k2+pk+q=0k^2 + pk + q = 0: за реалне различите корене yh=C1ek1x+C2ek2xy_h = C_1e^{k_1x} + C_2e^{k_2x}, за двоструки корен yh=(C1+C2x)ek1xy_h = (C_1 + C_2x)e^{k_1x}, а за k1/2=α±βik_{1/2} = \alpha \pm \beta i је yh=eαx(C1cos⁡βx+C2sin⁡βx)y_h = e^{\alpha x}(C_1\cos\beta x + C_2\sin\beta x).

Ако је f(x)=eax(Pn(x)cos⁡bx+Qm(x)sin⁡bx)f(x) = e^{ax}\left(P_n(x)\cos bx + Q_m(x)\sin bx\right), партикуларно решење тражите у облику yp=xreax(SN(x)cos⁡bx+TN(x)sin⁡bx)y_p = x^re^{ax}\left(S_N(x)\cos bx + T_N(x)\sin bx\right), где је N=max⁡{n,m}N = \max\{n, m\}, а rr вишеструкост броја a±bia \pm bi као корена карактеристичне једначине. Коефицијенте полинома SNS_N и TNT_N нађете уврштавањем ypy_p у једначину и изједначавањем коефицијената. Константе C1C_1 и C2C_2 из почетних услова одређујете тек у збиру yh+ypy_h + y_p. Честа грешка је r=0r = 0 кад a±bia \pm bi јесте корен: тада уврштавање даје немогућу једнакост.

0,511,52481216t (s)y (m)равнотежанајвише
4
20 m/s
Равнотежно истезање
2,5 m
Највеће истезање
12,81 m
Време до најниже тачке
0,91 s
Решење је збир yhy_h, осциловања, и yp=gω2y_p = \frac{g}{\omega^2}, равнотежног истезања (зелена линија), где је ω2=km\omega^2 = \frac{k}{m}. Тачка показује највеће истезање yp+yp2+v02ω2y_p + \sqrt{y_p^2 + \frac{v_0^2}{\omega^2}}; крива је нацртана док је уже затегнуто.

Овим једначинама инжењери описују осцилације опруга, ужади, амортизера и мостова, као и струју у колу са калемом и кондензатором; партикуларно решење је одговор система на спољашњу силу или напон.

Банџи скок

Скакач масе 75 kg скаче са платформе везан банџи ужетом крутости 300 N/m300\ \mathrm{N/m}. Кад уже почне да се затеже, скакач пада брзином 20 m/s. Нека је y(t)y(t) истезање ужета у метрима, а tt време у секундама од тог тренутка; отпор ваздуха занемарујемо, уже посматрамо као опругу која се истеже по Хуковом закону и узимамо g=10g = 10 m/s². Колико се уже највише истегне?

Решење

  1. По другом Њутновом закону 75y′′=75⋅10−300y75y'' = 75 \cdot 10 - 300y, тј. y′′+4y=10y'' + 4y = 10, уз почетне услове y(0)=0y(0) = 0 и y′(0)=20y'(0) = 20.
  2. Хомогени део: k2+4=0k^2 + 4 = 0, k1/2=±2ik_{1/2} = \pm 2i, па је yh=C1cos⁡2t+C2sin⁡2ty_h = C_1\cos 2t + C_2\sin 2t.
  3. Десна страна је 10=e0t⋅1010 = e^{0t} \cdot 10: број a±bi=0a \pm bi = 0 није корен, па је r=0r = 0, N=0N = 0 и yp=Ay_p = A. Уврштавањем 4A=104A = 10, тј. A=2,5A = 2{,}5.
  4. Опште решење је y=C1cos⁡2t+C2sin⁡2t+2,5y = C_1\cos 2t + C_2\sin 2t + 2{,}5. Из y(0)=0y(0) = 0 је C1=−2,5C_1 = -2{,}5, а из y′(0)=2C2=20y'(0) = 2C_2 = 20 је C2=10C_2 = 10.
  5. Дакле y=2,5−2,5cos⁡2t+10sin⁡2ty = 2{,}5 - 2{,}5\cos 2t + 10\sin 2t. Збир −2,5cos⁡2t+10sin⁡2t-2{,}5\cos 2t + 10\sin 2t највише је 2,52+102≈10,31\sqrt{2{,}5^2 + 10^2} \approx 10{,}31, па је највеће истезање 2,5+10,31≈12,812{,}5 + 10{,}31 \approx 12{,}81 m.

Одговор: Уже се највише истегне за око 12,8 m, а yp=2,5y_p = 2{,}5 m је равнотежни положај око ког скакач осцилује.

Задатак за вежбу

Одредите опште решење једначине y′′+y′−2y=6exy'' + y' - 2y = 6e^x.

Прикажи решење

k2+k−2=0k^2 + k - 2 = 0 даје k1=1k_1 = 1, k2=−2k_2 = -2. Пошто је a=1a = 1 корен, r=1r = 1 и yp=Axexy_p = Axe^x; уврштавањем 3Aex=6ex3Ae^x = 6e^x, па је A=2A = 2. Решење је y=C1ex+C2e−2x+2xexy = C_1e^x + C_2e^{-2x} + 2xe^x.