Дифференциальные уравнения первого порядка

Общий вид дифференциального уравнения первого порядка:

F(x,y,y′)=0F\left( {x,y,y'} \right) = 0.

Его общее решение имеет вид:

y=f(x,c)y = f\left( {x,c} \right).

Простейшие дифференциальные уравнения первого порядка.

  1. Дифференциальные уравнения с разделёнными перемеными: f1(x)dx=f2(y)dy,{f_1}\left( x \right)dx = {f_2}\left( y \right)dy, где множителем при dxdx является функция, зависящая только от xx, а множителем при dydy является функция, зависящая только от yy. Решение находится методом интегрирования обеих частей.∫f1(x)dx=∫f2(y)dy+C\int {{f_1}\left( x \right)dx} = \int {{f_2}\left( y \right)dy} + C
  2. Дифференциальные уравнения вида y′=f1(x)f2(y)dy,y' = {f_1}\left( x \right){f_2}\left( y \right)dy, где правая часть представляет собой произведение двух функций, из которых одна не зависит от xx, а вторая не зависит от yy, называется уравнением с разделяющимися переменными. Метод решения:∫dyf2(x)=∫f1(x)dx+C\int {\frac{{dy}}{{{f_2}\left( x \right)}}} = \int {{f_1}\left( x \right)dx} + C
  3. Дифференциальные уравнения с разделяющимися переменными, записанные в форме дифференциалов: f1(x)⋅f2(y)dx+f3(x)⋅f4(y)dy=0{f_1}\left( x \right) \cdot {f_2}\left( y \right)dx + {f_3}\left( x \right) \cdot {f_4}\left( y \right)dy = 0 для решения таких дифференциальных уравнений их надо привести к уравнениям с разделёнными переменными.∫f1(x)f3(x)dx+∫f4(y)f2(y)dy=C\int {\frac{{{f_1}\left( x \right)}}{{{f_3}\left( x \right)}}} dx + \int {\frac{{{f_4}\left( y \right)}}{{{f_2}\left( y \right)}}dy} = C

Функция f(x,y)f\left( {x,y} \right) называется однородной функцией n-го измерения, если при замене в ней переменных xx и yy соответственно на txtx и tyty, где  tt — произвольная величина (параметр) получается та же функция, умноженная на tn{t^n}, т.е. если выполняется условие: f(tx,ty)=tn⋅f(x,y)f\left( {tx,ty} \right) = {t^n} \cdot f\left( {x,y} \right) n- степень однородности уравнения.

Однородная функция степени nn представима в виде f(x,y)=xn⋅φ(yx)f\left( {x,y} \right) = {x^n} \cdot \varphi \left( {\frac{y}{x}} \right)

Однородная функция нулевой степени может быть записана в виде f(x,y)=φ(yx)f\left( {x,y} \right) = \varphi \left( {\frac{y}{x}} \right)

Если функции M(x,y)M\left( {x,y} \right) и N(x,y)N\left( {x,y} \right) однородные одной и той же степени nn, то дифференциальное уравнение M(x,y)dx+N(x,y)dy=0M\left( {x,y} \right)dx + N\left( {x,y} \right)dy = 0 называется однородным.

Уравнение y′=f(x,y)y' = f\left( {x,y} \right) называется однородным, если оно имеет вид: y′=φ(yx)y' = \varphi \left( {\frac{y}{x}} \right)

Очевидно, что f(x,y)f\left( {x,y} \right) однородная функция нулевого измерения.

Однородные уравнения приводятся к уравнению с разделяющимися перемеными при помощи подстановки.

t=yxt = \frac{y}{x} т.е. y=txy = tx и y′=t′x+ty' = t'x + t или в дифференциалах dy=tdx+xdtdy = tdx + xdt.

Линейные уравнения первого порядка.

Линейным дифференциальным уравнением первого порядка называется такое дифференциальное уравнение, в которое неизвестные функции  yy и y′{y'} входят в первых степенях и не перемножаются между собой.

Общий вид линейного уравнения первого порядка: y′+P(x)⋅y=Q(x)y' + P\left( x \right) \cdot y = Q\left( x \right)

Если Q(x)=0Q\left( x \right) = 0, то это линейное однородное уравнение с разделяющимися переменными.

Методы решения:метод Бернулли и метод Лагранжа.

Метод Бернулли

  1. Будем искать решение в виде y=U⋅Vy = U \cdot V, тогда y′=U′V+V′Uy' = U'V + V'U или dy=Vdu+Udvdy = Vdu + Udv (это подстановка Бернулли, где vv — вспомогательная функция.) Пример. xy′−2y=2x4x(U′V+V′U)−2UV=2x4xU′V+xV′U−2UV=2x4\begin{aligned}&xy' - 2y = 2{x^4} \\ &x\left( {U'V + V'U} \right) - 2UV = 2{x^4} \\ &xU'V + xV'U - 2UV = 2{x^4}\end{aligned}
  2. xU′V+U(xV′−2V)=2x4xU'V + U\left( {xV' - 2V} \right) = 2{x^4} найдём функцию VV таким образом, чтобы выражение в скобках было равно нулю. xV′−2V=0xdVdx=2V∫dVdx=2∫dxx\begin{aligned}&xV' - 2V = 0 \\ &x\frac{{dV}}{{dx}} = 2V \\ &\int {\frac{{dV}}{{dx}} = 2\int {\frac{{dx}}{x}} }\end{aligned} Интегрируя уравнение, получим ln⁡V=ln⁡x2⇒V=x2.\ln V = \ln {x^2} \Rightarrow V = {x^2}. Поскольку функция VV выбрана, чтобы удовлетворять определённому условию мы опускаем постоянную СС. Полученное выражение подставляем в исходное уравнение x⋅dUdx⋅x2=2x4dUdx=x⇒U=∫xdx=x2+C\begin{aligned}&x \cdot \frac{{dU}}{{dx}} \cdot {x^2} = 2{x^4} \\ &\frac{{dU}}{{dx}} = x \Rightarrow U = \int {xdx = {x^2}} + C\end{aligned} Объединив полученные выражения для VV UU в подстановке Бернулли, получим окончательное общее решение уравнения y=x2(x2+C).y = {x^2}\left( {{x^2} + C} \right).

Метод вариации произвольной постоянной. (Метод Лагранжа)

Покажем применение метода на том же примере.

xy′−2y=2x4xy' - 2y = 2{x^4}  Сначала решаем данное уравнение без правой части:xy′−2y=2x4xdydx=2yxy′−2y=0dyy=2dxxln⁡y=2ln⁡x+ln⁡Cy=C⋅x2\begin{aligned}&xy' - 2y = 2{x^4} \\ &x\frac{{dy}}{{dx}} = 2y \\ &xy' - 2y = 0 \\ &\frac{{dy}}{y} = 2\frac{{dx}}{x} \\ &\ln y = 2\ln x + \ln C \\ &y = C \cdot {x^2}\end{aligned} Пусть C=C(x)C = C\left( x \right) — некоторая неизвестная функция в исходном уравнении, тогда y=x2⋅C(x)y′=2x⋅C(x)+x2⋅C′(x)\begin{aligned}&y = {x^2} \cdot C\left( x \right) \\ &y' = 2x \cdot C\left( x \right) + {x^2} \cdot C'\left( x \right)\end{aligned} Подставляем в исходное уравнениеx⋅2x⋅C(x)+x⋅x2⋅C′(x)−2x2⋅C(x)=2x4∣÷(x2)2C(x)+x⋅C′(x)−2C(x)=2x2C′(x)=2x⇒C(x)=∫2xdx=x2+C\begin{aligned}&x \cdot 2x \cdot C\left( x \right) + x \cdot {x^2} \cdot C'\left( x \right) - 2{x^2} \cdot C\left( x \right) = 2{x^4}\left| { \div \left( {{x^2}} \right)} \right. \\ &2C\left( x \right) + x \cdot C'\left( x \right) - 2C\left( x \right) = 2{x^2} \\ &C'\left( x \right) = 2x \Rightarrow C\left( x \right) = \int {2xdx = {x^2}} + C\end{aligned} Подставляем полученное выражение в y=x2⋅C(x)y = {x^2} \cdot C\left( x \right) и получает окончательное решение y=x2⋅(x2+C).y = {x^2} \cdot \left( {{x^2} + C} \right).

Уравнение Бернулли.

Общий вид уравнения:y′+P(x)y=Q(x)yny' + P\left( x \right)y = Q\left( x \right){y^n} слева линейное выражение, а справа присутсвует множитель yn(n=const){y^n}\left( {n = const} \right). Умножим обе части на  1yn\frac{1}{{{y^n}}} 1yndydx+P(x)1yn−1=Q(x)\frac{1}{{{y^n}}}\frac{{dy}}{{dx}} + P\left( x \right)\frac{1}{{{y^{n - 1}}}} = Q\left( x \right) Применим подстановку

z(x)=1yn−1z\left( x \right) = \frac{1}{{{y^{n - 1}}}} и dzdx=z′=(1−n)⋅1yn⋅y′\frac{{dz}}{{dx}} = z' = \left( {1 - n} \right) \cdot \frac{1}{{{y^n}}} \cdot y' 1yn⋅y′=z′(1−n)\frac{1}{{{y^n}}} \cdot y' = \frac{{z'}}{{\left( {1 - n} \right)}} получим дифференциальное уравнение вида:z′(1−n)+P(x)z(x)=Q(x)dzdx+(1−n)P(x)z(x)=Q(x)(1−n)\begin{aligned}&\frac{{z'}}{{\left( {1 - n} \right)}} + P\left( x \right)z\left( x \right) = Q\left( x \right) \\ &\frac{{dz}}{{dx}} + \left( {1 - n} \right)P\left( x \right)z\left( x \right) = Q\left( x \right)\left( {1 - n} \right)\end{aligned}Это линейное уравнение I-го порядка, для его решения применяем, например, подстановку Бернулли.

Дифференциальное уравнение первого порядка F(x,y,y′)=0F(x, y, y') = 0 имеет общее решение y=f(x,C)y = f(x, C) с одной постоянной. Главные типы: уравнение с разделяющимися переменными, которое приводится к виду f1(x) dx=f2(y) dyf_1(x)\,dx = f_2(y)\,dy и решается интегрированием обеих частей; однородное y′=φ(yx)y' = \varphi\left(\frac{y}{x}\right), которое подстановка t=yxt = \frac{y}{x} сводит к предыдущему типу; линейное y′+P(x)y=Q(x)y' + P(x)y = Q(x), которое решают подстановкой Бернулли y=UVy = UV или методом Лагранжа; и уравнение Бернулли y′+P(x)y=Q(x)yny' + P(x)y = Q(x)y^n, которое подстановка z=y1−nz = y^{1-n} сводит к линейному.

Сначала определите тип уравнения, а уже потом выбирайте метод. В методе Бернулли функцию VV выбирают так, чтобы V′+PV=0V' + PV = 0, и для UU остаётся уравнение с разделяющимися переменными. Постоянную интегрирования пишите сразу, в одной части равенства, а её значение находите из начального условия. Без постоянной решение не будет общим.

510152025303520406080100t (мин)T (°C)комната
90°C
0,067
5 мин
Температура T(t)
70,1°C
Скорость остывания T′(t)
−3,35°C/мин
Кривая — решение T=20+(T0−20)e−ktT = 20 + (T_0 - 20)e^{-kt}. Чем горячее чай, тем быстрее он остывает: скорость T′T' пропорциональна разности T−20T - 20, поэтому кривая приближается к комнатной температуре, но не достигает её. Начальное положение соответствует примеру: точка на 5-й минуте показывает около 70 ∘C70\,^\circ\mathrm{C}.

Остывание тела, рост числа бактерий и снижение концентрации лекарства в крови описываются уравнениями первого порядка, поэтому ими постоянно пользуются инженеры, биологи и фармацевты.

Чай остывает

Чай температурой 90 ∘C90\,^\circ\mathrm{C} оставили остывать в комнате с температурой воздуха 20 ∘C20\,^\circ\mathrm{C}. По закону охлаждения Ньютона температура T(t)T(t) удовлетворяет уравнению T′=−k(T−20)T' = -k(T - 20), где tt — время в минутах. Через 5 минут температура чая равна 70 ∘C70\,^\circ\mathrm{C}. Когда он остынет до 50 ∘C50\,^\circ\mathrm{C}?

Решение

  1. Разделим переменные: dTT−20=−k dt\frac{dT}{T - 20} = -k\,dt.
  2. Интегрируя, получаем ln⁡(T−20)=−kt+C\ln(T - 20) = -kt + C, то есть T=20+Ae−ktT = 20 + Ae^{-kt}. Из T(0)=90T(0) = 90 следует A=70A = 70.
  3. Из T(5)=70T(5) = 70 следует 50=70e−5k50 = 70e^{-5k}, поэтому e−5k=57e^{-5k} = \frac{5}{7} и k=15ln⁡75≈0,0673k = \frac{1}{5}\ln\frac{7}{5} \approx 0{,}0673.
  4. Ищем tt, при котором 20+70e−kt=5020 + 70e^{-kt} = 50, то есть e−kt=37e^{-kt} = \frac{3}{7}. Логарифмируя, получаем t=1kln⁡73=5ln⁡73ln⁡75≈12,6t = \frac{1}{k}\ln\frac{7}{3} = \frac{5\ln\frac{7}{3}}{\ln\frac{7}{5}} \approx 12{,}6 мин.

Ответ: Чай остынет до 50 ∘C50\,^\circ\mathrm{C} примерно через 12,6 минуты после начала.

Задача для тренировки

Суп температурой 80 ∘C80\,^\circ\mathrm{C} остывает в комнате при 20 ∘C20\,^\circ\mathrm{C}, и через 10 минут его температура равна 60 ∘C60\,^\circ\mathrm{C}. Когда она станет 40 ∘C40\,^\circ\mathrm{C}?

Показать решение

T=20+60e−ktT = 20 + 60e^{-kt}. Из 60=20+60e−10k60 = 20 + 60e^{-10k} следует e−10k=23e^{-10k} = \frac{2}{3}, а из 40=20+60e−kt40 = 20 + 60e^{-kt} следует e−kt=13e^{-kt} = \frac{1}{3}. Поэтому t=10ln⁡3ln⁡1,5≈27,1t = \frac{10\ln 3}{\ln 1{,}5} \approx 27{,}1 минуты.