Учебник. Понятие дифференциального уравнения

Рассмотрим движение тела массы m в вязкой среде с коэффициентом сопротивления k. По второму закону Ньютона можно записать: ma = –kv.

Так как ускорение – первая производная скорости, $a=\stackrel{ċ}{v}\mathrm{,}$ $m\stackrel{ċ}{v}=\mathrm{-}kv\mathrm{.}$

В это уравнение входит неизвестная величина v и ее производная по времени $\dot{v}\mathrm{.}$

Движение тела в вязкой среде в отсутствие движущей силы

Уравнения, подобные этому, нередко встречаются в физике, химии, экономике и других дисциплинах. Уравнение вида $f\left(x\mathrm{, }y\mathrm{, }{y}^{\mathrm{\prime }}\mathrm{, }{y}^{\mathrm{\prime \prime }}\mathrm{, }\mathrm{...}\mathrm{, }{y}^{\left(n\right)}\right)=0$ называется обыкновенным дифференциальным уравнением n-го порядка. Решением этого уравнения является произвольная функция y = y (x), подстановка которой в уравнение превращает его в верное тождество. График решения дифференциального уравнения называется интегральной кривой.

Простейшим дифференциальным уравнением является уравнение вида $y^{\mathrm{\prime }}=f\left(x\right)\mathrm{.}$

Чтобы его решить, нужно представить производную $y^{\mathrm{\prime }}$ как $\frac{dy}{dx}\mathrm{,}$ домножить обе части уравнения на dx и проинтегрировать обе части получившегося уравнения: $y=\int f\left(x\right)dx+C$

Как видно, это уравнение имеет бесконечное количество решений, отличающихся друг от друга на постоянную C. Выбрать конкретное решение уравнения можно, если знать начальные условия, например, точку, через которую проходит график функции y = y (x). Так, если известно, что y (x₀) = y₀, то подставляя это значение в общее решение $y=\int f\left(x\right)dx+C=F\left(x\right)+C\mathrm{,}$ получаем $y_0=F\left(x_0\right)+C\mathrm{,}$  откуда $C=y_0-F\left(x_0\right)$ и $y=y_0+F\left(x\right)-F\left(x_0\right)\mathrm{.}$ Это решение можно записать в виде $y=y_0+\underset{x_0}{\overset{x}{\int }}f\left(t\right)dt$.

Общим решением дифференциального уравнения называется функция y = y (x, C₁, C₂,…, C_n), зависящая от n констант, если она является решением дифференциального уравнения при любых значениях постоянных C₁, C₂,…, C_n.

Неоднозначность общего решения многих уравнений имеет простой физический смысл. Так, дифференциальное уравнение движения материальной точки массы m под действием силы F (второй закон Ньютона) $m\stackrel{ċċ}{x}=F$ не определяет однозначно закон движения этой точки: для этого необходимо знать его начальные скорость и координату.

Частным решением дифференциального уравнения называется решение, получаемое из общего решения путем придания определенного значения постоянным C_i.

Решения уравнения (y')² + y² = 1.

Для исследования решений дифференциального уравнения применяют метод фазовых траекторий. Для этого по осям абсцисс и ординат откладывают, соответственно, саму величину и ее производную. Фазовые траектории могут выглядеть самым причудливым образом даже для линейных дифференциальных уравнений; траектории для уравнения $\frac{dy}{dx}=\frac{y}{x}$  показаны на графике.

Наряду с частными существуют особые решения, которые нельзя получить из общего решения никакой подстановкой постоянных. Так, общим решением уравнения ${\left(y^{\mathrm{\prime }}\right)}^2+y^2=1$ будет семейство функций y = cos (x + C). Кроме того, это уравнение имеет два особых решения y = 1 и y = –1, в чем можно убедиться подстановкой. Графики особых решений в каждой точке касаются проходящего через эту точку частного решения.

Если задачу об отыскании всех решений дифференциального уравнения удается свести к вычислению конечного числа интегралов и производных от известных функций и к алгебраическим операциям, то говорят, что уравнение интегрируется в квадратурах. К сожалению, класс таких уравнений крайне узок: так, второй закон Ньютона $m\stackrel{ċċ}{x}=F$ интегрируется только тогда, когда F является функцией только от одной из переменных t, x, $\dot{x}\mathrm{.}$ Поэтому для исследования дифференциальных уравнений широко применяются приближенные и численные методы.

Во многих случаях приходится составлять дифференциальное уравнение, чтобы решить задачу. Рассмотрим несколько примеров.

1. Реактивное движение.

   Реактивное движение. Ракета, начальная масса которой равна m₀, движется вдали от небесных тел за счет газовой струи, выбрасываемой в направлении, противоположном движению, со скоростью u относительно ракеты. Начальная скорость ракеты равна v₀. Необходимо получить зависимость скорости ракеты от ее массы.

   Запишем закон сохранения импульса. Пусть в некоторый момент времени t импульс ракеты равнялся mv. За бесконечно малый промежуток времени dt скорость ракеты увеличилась от v до v + dv, а масса уменьшилась от m до m – dm, и ракета увеличила свой импульс на $\left(m-dm\right)\left(v+dv\right)\mathrm{,}$ выбросив назад dm топлива, импульс которого равен dm (v – u) (здесь и далее уравнение записывается в проекции на ось движения, u > 0). Получаем $mv=\left(m-dm\right)\left(v+dv\right)+dm\left(v-u\right)\mathrm{.}$

   Раскрывая скобки и приводя подобные члены, получаем $u\mathrm{ d}m=m\mathrm{ d}v-dm\mathrm{ d}v\mathrm{.}$ Слагаемое dm · dv бесконечно мало по сравнению с остальными слагаемыми, так как оно является произведением двух бесконечно малых величин; этим слагаемым можно пренебречь. Разделим обе части получившегося уравнения $u\mathrm{ d}m=m\mathrm{ d}v$ на dt: $u\frac{dm}{dt}=ma$

   Это уравнение называется уравнением движения тела с переменной массой. Это дифференциальное уравнение.

2. Истечение идеальной жидкости из сосуда.

   Истечение идеальной жидкости из сосуда

   Пусть из сосуда с постоянной площадью поперечного сечения S через отверстие в дне площадью S₀ вытекает вода. Пренебрегая вязкостью, записать уравнение изменения уровня жидкости в сосуде h.

   Скорость v истечения жидкости из сосуда, уровень жидкости в котором равен h, равна $v=\sqrt{2gh}$ (формула Торричелли). За бесконечно малый промежуток времени dt высота жидкости в сосуде уменьшается на dh. Этот же объем S dh истекает через отверстие площадью S₀. Можно записать, что за время dt через отверстие протекает «столб» воды dx = v dt. Таким образом $\mathrm{-}Sdh=S_{0}vdt=S_0\sqrt{2gh}dt$ (знак «минус» поставлен из-за того, что уровень воды в сосуде уменьшается). Таким образом, $\frac{dh}{dt}=\mathrm{-}\frac{S}{S_0}\sqrt{2gh}\mathrm{.}$ Это тоже дифференциальное уравнение.

Движение в воздушной среде