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

y′=f(a1x+b1y+c1a2x+b2y+c2)ако је δ=∣a1b1a2b2∣≠0уводимо смену: x=X+α⇒dx=dXy=Y+β⇒dy=dYгде су α и β решења система a1α+b1β+c1=0a2α+b2β+c2=0ако је δ=∣a1b1a2b2∣=0уводимо смену: a2x+b2y+c2=t, где је t=t(x)\begin{array}{l} y' = f\left( {\frac{{{a_1}x + {b_1}y + {c_1}}}{{{a_2}x + {b_2}y + {c_2}}}} \right) \\ {\text{ако је }} \delta = \left| {\begin{array}{cccccccccccccccccccc} {a_1}&{b_1}\\{a_2}&{b_2} \end{array}} \right| \ne 0 \\ {\text{уводимо смену: }} \\ x = X + \alpha {\rm{ }} \Rightarrow {\rm{ }}dx = dX \\ y = Y + \beta {\rm{ }} \Rightarrow {\rm{ }}dy = dY \\ {\text{где су }} \alpha {\text{ и }} \beta {\text{ решења система }} \\ {a_1}\alpha + {b_1}\beta + {c_1} = 0 \\ {a_2}\alpha + {b_2}\beta + {c_2} = 0 \\ \\ {\text{ако је }} \delta = \left| {\begin{array}{cccccccccccccccccccc} {a_1}&{b_1}\\{a_2}&{b_2} \end{array}} \right| = 0 \\ {\text{уводимо смену: }} {a_2}x + {b_2}y + {c_2} = t,{\rm{ }}{\text{ где је }}t = t\left( x \right) \end{array}

Задаци


1. Наћи опште решење диференцијалне једначине  y′=x−y+1x+y−3y' = \frac{{x - y + 1}}{{x + y - 3}}.

Како је δ=∣1−111∣=1−(−1)=2≠0\delta = \left| {\begin{array}{cccccccccccccccccccc} 1&{ - 1}\\ 1&1 \end{array}} \right| = 1 - \left( { - 1} \right) = 2 \ne 0

и решења система α−β+1=0α+β−3=0‾\begin{array}{cccccccccccccccccccc} {\alpha - \beta + 1 = 0}\\ {\underline {\alpha + \beta - 3 = 0} } \end{array}

α=1,β=2\alpha = 1,\beta = 2

уводимо смену:

x=X+1⇒dx=dXx = X + 1{\rm{ }} \Rightarrow {\rm{ }}dx = dX

y=Y+2⇒dy=dYy = Y + 2{\rm{ }} \Rightarrow {\rm{ }}dy = dY

при чему важи да је y′=dydx=dYdX=Y′y' = \frac{{dy}}{{dx}} = \frac{{dY}}{{dX}} = Y'

Y′=X+1−Y−2+1X+1+Y+2−3=X−YX+YY' = \frac{{X + 1 - Y - 2 + 1}}{{X + 1 + Y + 2 - 3}} = \frac{{X - Y}}{{X + Y}}

Y′=X(1−YX)X(1+YX)Y' = \frac{{X\left( {1 - \frac{Y}{X}} \right)}}{{X\left( {1 + \frac{Y}{X}} \right)}}

На овај начин смо добили хомогену диференцијалну једначину коју решавамо уобичајеном сменом,

YX=t⇒Y=tX⇒Y′=t′X+t\frac{Y}{X} = t \Rightarrow Y = tX \Rightarrow Y' = t'X + t

t′X+t=1−t1+tt'X + t = \frac{{1 - t}}{{1 + t}}

dtdXX=1−t−t−t21+t\frac{{dt}}{{dX}}X = \frac{{1 - t - t - {t^2}}}{{1 + t}}

(1+t)dt1−2t−t2=dXX  /∫\frac{{\left( {1 + t} \right)dt}}{{1 - 2t - {t^2}}} = \frac{{dX}}{X}{\rm{ }} {\text{ }} {\text{ }} /\int {}

∫(1+t)dt1−2t−t2=∫dXX\int {\frac{{\left( {1 + t} \right)dt}}{{1 - 2t - {t^2}}} = \int {\frac{{dX}}{X}} }

∣1−2t−t2=s(−2−2t)dt=ds(1+t)dt=ds−2∣\left| {\begin{array}{cccccccccccccccccccc} {1 - 2t - {t^2} = s}\\ {\left( { - 2 - 2t} \right)dt = ds}\\ {\left( {1 + t} \right)dt = \frac{{ds}}{{ - 2}}} \end{array}} \right|

∫ds−2s=∫dXX\int {\frac{{\frac{{ds}}{{ - 2}}}}{s}} = \int {\frac{{dX}}{X}}

−12ln⁡∣s∣=ln⁡X+ln⁡C- \frac{1}{2}\ln \left| s \right| = \ln X + \ln C

s=(XC)−2s = {\left( {XC} \right)^{ - 2}}

1−2t−t2=1X2C21 - 2t - {t^2} = \frac{1}{{{X^2}{C^2}}}

1−2YX−(YX)2=1X2C21 - 2\frac{Y}{X} - {\left( {\frac{Y}{X}} \right)^2} = \frac{1}{{{X^2}{C^2}}}

X2−2XY−Y2X2=1X2C2\frac{{{X^2} - 2XY - {Y^2}}}{{{X^2}}} = \frac{1}{{{X^2}{C^2}}}

Kako je x=X+1⇒X=x−1y=Y+2⇒Y=y−2\begin{array}{cccccccccccccccccccc} {x = X + 1{\rm{ }} \Rightarrow {\rm{ }}X = x - 1}\\ {y = Y + 2{\rm{ }} \Rightarrow {\rm{ }}Y = y - 2} \end{array}

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

(x−1)2−2(x−1)(y−2)−(y−2)2=C1{\left( {x - 1} \right)^2} - 2\left( {x - 1} \right)\left( {y - 2} \right) - {\left( {y - 2} \right)^2} = {C_1}


2. Наћи опште решење диференцијалне једначине  (2x+y+1)dx−(4x+2y−3)dy=0\left( {2x + y + 1} \right)dx - \left( {4x + 2y - 3} \right)dy = 0.

(2x+y+1)dx=(4x+2y−3)dy\left( {2x + y + 1} \right)dx = \left( {4x + 2y - 3} \right)dy

2x+y+14x+2y−3=dydx\frac{{2x + y + 1}}{{4x + 2y - 3}} = \frac{{dy}}{{dx}}

y′=2x+y+14x+2y−3y' = \frac{{2x + y + 1}}{{4x + 2y - 3}}

δ=∣2142∣=4−4=0\delta = \left| {\begin{array}{cccccccccccccccccccc} 2&1\\ 4&2 \end{array}} \right| = 4 - 4 = 0

4x+2y−3=t⇒y=t−4x+32⇒y′=12t′−24x + 2y - 3 = t \Rightarrow y = \frac{{t - 4x + 3}}{2} \Rightarrow y' = \frac{1}{2}t' - 2

2(2x+y)−3=t2\left( {2x + y} \right) - 3 = t

2x+y=t+322x + y = \frac{{t + 3}}{2}

12t′−2=t+32+1t\frac{1}{2}t' - 2 = \frac{{\frac{{t + 3}}{2} + 1}}{t}

12dtdx=t+3+22t+2\frac{1}{2}\frac{{dt}}{{dx}} = \frac{{\frac{{t + 3 + 2}}{2}}}{t} + 2

12dtdx=t+5+4t2t\frac{1}{2}\frac{{dt}}{{dx}} = \frac{{t + 5 + 4t}}{{2t}}

dtdx=5+5t2t⋅2\frac{{dt}}{{dx}} = \frac{{5 + 5t}}{{2t}} \cdot 2

tdt5+5t=dx/∫\frac{{tdt}}{{5 + 5t}} = dx{\rm{ }}/\int {}

∫t5+5tdt=∫dx\int {\frac{t}{{5 + 5t}}dt = \int {dx} }

15∫t+1−11+tdt=x+C\frac{1}{5}\int {\frac{{t + 1 - 1}}{{1 + t}}dt = x + C}

15(∫dt−∫dtt+1)=x+C\frac{1}{5}\left( {\int {dt} - \int {\frac{{dt}}{{t + 1}}} } \right) = x + C

15t−15ln⁡∣t+1∣=x+C\frac{1}{5}t - \frac{1}{5}\ln \left| {t + 1} \right| = x + C

4x+2y−35−15ln⁡∣4x+2y−3+1∣=x+C\frac{{4x + 2y - 3}}{5} - \frac{1}{5}\ln \left| {4x + 2y - 3 + 1} \right| = x + C

−x+2y−35−15ln⁡∣4x+2y−2∣=C\frac{{ - x + 2y - 3}}{5} - \frac{1}{5}\ln \left| {4x + 2y - 2} \right| = C

3. Наћи опште решење диференцијалне једначине  y′=2(y+2x+y−1)2y' = 2{\left( {\frac{{y + 2}}{{x + y - 1}}} \right)^2}.

Како је δ=∣0111∣=0−1=−1≠0\delta = \left| {\begin{array}{cccccccccccccccccccc} 0&1\\ 1&1 \end{array}} \right| = 0 - 1 = - 1 \ne 0

и решења система β+2=0α+β−1=0‾\begin{array}{cccccccccccccccccccc} {\beta + 2 = 0}\\ {\underline {\alpha + \beta - 1 = 0} } \end{array}

α=3,β=−2\alpha = 3,\beta = - 2

уводимо смену:

x=X+3⇒dx=dXx = X + 3{\rm{ }} \Rightarrow {\rm{ }}dx = dX

y=Y−2⇒dy=dYy = Y - 2{\rm{ }} \Rightarrow {\rm{ }}dy = dY

при чему важи да је y′=dydx=dYdX=Y′y' = \frac{{dy}}{{dx}} = \frac{{dY}}{{dX}} = Y'

Y′=2(Y−2+2X+3+Y−2−1)2=2(YX+Y)2Y' = 2{\left( {\frac{{Y - 2 + 2}}{{X + 3 + Y - 2 - 1}}} \right)^2} = 2{\left( {\frac{Y}{{X + Y}}} \right)^2}

Да би могли раздвојити променљиве овом задатку прећићемо на хомогену диференцијалну једначину по XX

X′=12(X+YY)2X' = \frac{1}{2}{\left( {\frac{{X + Y}}{Y}} \right)^2}

па уводимо смену XY=t⇒X=tY⇒X′=t′Y+t\frac{X}{Y} = t \Rightarrow X = tY \Rightarrow X' = t'Y + t

t′Y+t=12(t+1)2t'Y + t = \frac{1}{2}{\left( {t + 1} \right)^2}

dtdYY=t2+2t+12−t\frac{{dt}}{{dY}}Y = \frac{{{t^2} + 2t + 1}}{2} - t

dtdYY=t2+12\frac{{dt}}{{dY}}Y = \frac{{{t^2} + 1}}{2}

2dtt2+1=dYY/∫\frac{{2dt}}{{{t^2} + 1}} = \frac{{dY}}{Y}{\rm{ }}/\int {}

2∫dtt2+1=∫dYY2\int {\frac{{dt}}{{{t^2} + 1}} = \int {\frac{{dY}}{Y}} }

2arctgt=ln⁡∣Y∣+ln⁡C2arctgt = \ln \left| Y \right| + \ln C

2arctgXY=ln⁡YC2arctg\frac{X}{Y} = \ln YC

2arctgx−3y+2=ln⁡(y+2)C2arctg\frac{{x - 3}}{{y + 2}} = \ln \left( {y + 2} \right)C

4. Наћи опште решење диференцијалне једначине  x−y−1+(y−x+2)y′=0x - y - 1 + \left( {y - x + 2} \right)y' = 0.

(y−x+2)y′=−x+y+1\left( {y - x + 2} \right)y' = - x + y + 1

y′=−x+y+1y−x+2y' = \frac{{ - x + y + 1}}{{y - x + 2}}

Како је у овој диференцијалној једначини δ=∣−11−11∣=−1+1=0\delta = \left| {\begin{array}{cccccccccccccccccccc} { - 1}&1\\ { - 1}&1 \end{array}} \right| = - 1 + 1 = 0

уводимо смену

−x+y+2=t⇒y=t+x−2⇒y′=t′+1- x + y + 2 = t \Rightarrow y = t + x - 2 \Rightarrow y' = t' + 1

−x+y=t−2- x + y = t - 2

t′+1=t−2+1tt' + 1 = \frac{{t - 2 + 1}}{t}

t′=t−1−ttt' = \frac{{t - 1 - t}}{t}

dtdx=−1t\frac{{dt}}{{dx}} = - \frac{1}{t}

tdt=−dx/∫tdt = - dx{\rm{ }}/\int {}

t22=−x+C\frac{{{t^2}}}{2} = - x + C

(−x+y+2)2=−2x+2C{\left( { - x + y + 2} \right)^2} = - 2x + 2C

x2+y2+4−2xy−4x+4y=−2x+2C{x^2} + {y^2} + 4 - 2xy - 4x + 4y = - 2x + 2C

(x−y)2−2x+4y=C1{\left( {x - y} \right)^2} - 2x + 4y = {C_1}

Једначина y′=f(a1x+b1y+c1a2x+b2y+c2)y' = f\left(\frac{a_1x + b_1y + c_1}{a_2x + b_2y + c_2}\right) није хомогена због слободних чланова c1c_1 и c2c_2. Прво израчунајте δ=∣a1b1a2b2∣\delta = \begin{vmatrix} a_1 & b_1 \\ a_2 & b_2 \end{vmatrix}. Ако је δ≠0\delta \neq 0, праве a1x+b1y+c1=0a_1x + b_1y + c_1 = 0 и a2x+b2y+c2=0a_2x + b_2y + c_2 = 0 секу се у једној тачки (α,β)(\alpha, \beta), коју налазите из система a1α+b1β+c1=0a_1\alpha + b_1\beta + c_1 = 0, a2α+b2β+c2=0a_2\alpha + b_2\beta + c_2 = 0. Смена x=X+αx = X + \alpha, y=Y+βy = Y + \beta премешта координатни почетак у ту тачку, слободни чланови нестају, а Y′=y′Y' = y', па добијате хомогену једначину коју решавате сменом YX=t\frac{Y}{X} = t.

Ако је δ=0\delta = 0, бројилац и именилац зависе од истог израза, па смена t=a2x+b2y+c2t = a_2x + b_2y + c_2 даје једначину која раздваја променљиве. На крају вратите X=x−αX = x - \alpha и Y=y−βY = y - \beta: то се често заборави.

1020304050−551525x (m)y (m)5x + y − 160 = 0x − 5y + 20 = 0Mдрон
1
Растојање од споменика
5,69 m
Угао θ
−360°
Угао лета према споменику
78,7°
Праве бројиоца и имениоца секу се у споменику M(30,10)M(30, 10), па смена x=X+30x = X + 30, y=Y+10y = Y + 10 ставља почетак баш у MM. Око MM решење је спирала r=20eθ5r = 20e^{\frac{\theta}{5}}: дрон стално лети под истим углом према MM, а сваким кругом растојање се смањи на e−2π5≈0,28e^{-\frac{2\pi}{5}} \approx 0{,}28 претходног.

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

Дрон кружи око споменика

Дрон снима споменик на тргу. Програм га води тако да његов правац лета увек заклапа исти угао са правцем ка споменику, па његова путања у координатном систему трга (јединица је метар) задовољава y′=5x+y−160x−5y+20y' = \frac{5x + y - 160}{x - 5y + 20}. Дрон креће из тачке (50,10)(50, 10) и спирално се приближава споменику. Где је споменик и колико је дрон удаљен од њега кад око њега направи један пун круг?

Решење

  1. Овде је δ=∣511−5∣=−26≠0\delta = \begin{vmatrix} 5 & 1 \\ 1 & -5 \end{vmatrix} = -26 \neq 0, па решавамо систем 5α+β−160=05\alpha + \beta - 160 = 0, α−5β+20=0\alpha - 5\beta + 20 = 0: α=30\alpha = 30, β=10\beta = 10. Праве бројиоца и имениоца секу се у споменику (30,10)(30, 10).
  2. Смена x=X+30x = X + 30, y=Y+10y = Y + 10 ставља координатни почетак у споменик и даје хомогену једначину Y′=5X+YX−5YY' = \frac{5X + Y}{X - 5Y}.
  3. Смена YX=t\frac{Y}{X} = t, Y′=t′X+tY' = t'X + t: t′X=5+t1−5t−t=5(1+t2)1−5tt'X = \frac{5 + t}{1 - 5t} - t = \frac{5(1 + t^2)}{1 - 5t}, па је (1−5t) dt1+t2=5 dXX\frac{(1 - 5t)\,dt}{1 + t^2} = \frac{5\,dX}{X}.
  4. Интеграцијом: arctg⁡t−52ln⁡(1+t2)=5ln⁡∣X∣+C\operatorname{arctg} t - \frac{5}{2}\ln(1 + t^2) = 5\ln|X| + C, тј. arctg⁡YX=5ln⁡X2+Y2+C\operatorname{arctg}\frac{Y}{X} = 5\ln\sqrt{X^2 + Y^2} + C. Са углом θ\theta и растојањем rr од споменика то је θ=5ln⁡r+C\theta = 5\ln r + C, односно спирала r=Aeθ5r = Ae^{\frac{\theta}{5}}.
  5. На почетку је r=20r = 20 и θ=0\theta = 0, па је A=20A = 20. Дрон се приближава споменику, па после пуног круга θ=−2π\theta = -2\pi и r=20e−2π5≈5,69r = 20e^{-\frac{2\pi}{5}} \approx 5{,}69 m.

Одговор: Споменик је у тачки (30,10)(30, 10), а после једног круга дрон је од њега удаљен око 5,69 m: сваким кругом растојање се смањи на око 28% претходног.

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

Одредите опште решење једначине y′=x+y+1x+y−1y' = \frac{x + y + 1}{x + y - 1}.

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

δ=0\delta = 0, па уводимо t=x+y−1t = x + y - 1: t′=1+y′=1+t+2t=2(t+1)tt' = 1 + y' = 1 + \frac{t + 2}{t} = \frac{2(t + 1)}{t}, тј. t dtt+1=2 dx\frac{t\,dt}{t + 1} = 2\,dx. Интеграцијом t−ln⁡∣t+1∣=2x+Ct - \ln|t + 1| = 2x + C, па је y−x−ln⁡∣x+y∣=C1y - x - \ln|x + y| = C_1.