Определение

Уравнение первого порядка, разрешённое относительно производной:

dydx=M(x,y)N(x,y),\frac{dy}{dx} = \frac{M(x, y)}{N(x, y)},

называется уравнением с разделяющимися переменными, если функции MM и NN можно представить в виде:

M(x,y)=M1(x)M2(y),N(x,y)=N1(x)N2(y).M(x, y) = M_1(x) M_2(y), \qquad N(x, y) = N_1(x) N_2(y).

Тогда уравнение принимает вид

dydx=M1(x)M2(y)N1(x)N2(y).\frac{dy}{dx} = \frac{M_1(x) M_2(y)}{N_1(x) N_2(y)}.

Приведение к разделённому виду и общий интеграл

При условии, что N1(x)0N_1(x) \neq 0 и M2(y)0M_2(y) \neq 0, можно разделить переменные:

M1(x)N1(x)dx=N2(y)M2(y)dy.\frac{M_1(x)}{N_1(x)}\,dx = \frac{N_2(y)}{M_2(y)}\,dy.

Интегрируя, получаем общий интеграл (решение в неявной форме):

M1(x)N1(x)dxN2(y)M2(y)dy=C,\int \frac{M_1(x)}{N_1(x)}\,dx - \int \frac{N_2(y)}{M_2(y)}\,dy = C,

где CC — произвольная постоянная.

Часто уравнение записывают в виде

y=f1(x)f2(y),y' = f_1(x) f_2(y),

где f1(x)=M1(x)N1(x)f_1(x) = \frac{M_1(x)}{N_1(x)},   f2(y)=M2(y)N2(y)\; f_2(y) = \frac{M_2(y)}{N_2(y)}.


Теорема о существовании и единственности

Пусть функции f1(x)f_1(x) и f2(y)f_2(y) непрерывны в области R=(a,b)×(c,d)R = (a, b) \times (c, d) и f2(y)0f_2(y) \neq 0 во всей области. Тогда через каждую точку (x0,y0)R(x_0, y_0) \in R проходит единственное решение уравнения y=f1(x)f2(y)y' = f_1(x) f_2(y).

Это решение задаётся формулой

x0xf1(s)ds=y0ydtf2(t).\int_{x_0}^x f_1(s)\,ds = \int_{y_0}^y \frac{dt}{f_2(t)}.

Критерий особого решения

При приведении к разделённому виду мы делим на M2(y)M_2(y). Поэтому могут быть потеряны постоянные решения yηy\equiv\eta, где M2(η)=0.M_2(\eta)=0.

Решение y=ηy=\eta рассматривается на тех промежутках, на которых N1(x)0N_1(x)\neq0, поскольку в точках N1(x)=0N_1(x)=0 правая часть исходного уравнения не определена.

Выясним, при каких условиях решение yηy\equiv\eta является особым. Возьмём произвольную опорную точку

(x0,y0)R,y0η,(x_0,y_0)\in R, \qquad y_0\neq\eta,

и запишем общий интеграл с переменными верхними пределами:

x0xM1(s)N1(s)dsy0yN2(t)M2(t)dt=C\int_{x_0}^{x} \frac{M_1(s)}{N_1(s)}ds - \int_{y_0}^{y} \frac{N_2(t)}{M_2(t)}dt = C

При этом интегральная кривая, проходящая через точку (x0,y0)(x_0,y_0), соответствует значению C=0C=0.

Обозначим левую часть равенства через

ψ(x,y)=x0xM1(s)N1(s)dsy0yN2(t)M2(t)dt\psi(x,y) = \int_{x_0}^{x} \frac{M_1(s)}{N_1(s)}ds - \int_{y_0}^{y} \frac{N_2(t)}{M_2(t)}dt

Рассмотрим поведение функции ψ(x,y)\psi(x,y) при приближении yy к значению η\eta.

Случай 1: интеграл расходится

Пусть

limyηy0yN2(t)M2(t)dt=+\lim_{y\to\eta} \int_{y_0}^{y} \frac{N_2(t)}{M_2(t)}dt = +\infty

или

limyηy0yN2(t)M2(t)dt=.\lim_{y\to\eta} \int_{y_0}^{y} \frac{N_2(t)}{M_2(t)}dt = -\infty.

Тогда при фиксированном конечном xx величина ψ(x,y)\psi(x,y) также стремится к бесконечности с соответствующим знаком.

Следовательно, решение yηy\equiv\eta входит в формулу общего интеграла как предельный случай при бесконечном значении постоянной: C+C\to+\infty или CC\to-\infty. Такое решение является частным.


Случай 2: интеграл сходится

Пусть

limyηy0yN2(t)M2(t)dt=A,\lim_{y\to\eta} \int_{y_0}^{y} \frac{N_2(t)}{M_2(t)}dt = A,

где AA — конечное число.

Возьмём произвольную точку x1,ηx_1,\eta постоянного решения y=ηy=\eta. При yη,xx1y\to\eta, x\to x_1 из формулы (2) получаем

ψ(x,y)x0x1M1(s)N1(s)dsA.\psi(x,y) \longrightarrow \int_{x_0}^{x_1} \frac{M_1(s)}{N_1(s)}ds-A.

Обозначим это значение через

C(x1)=x0x1M1(s)N1(s)dsA.C(x_1) = \int_{x_0}^{x_1} \frac{M_1(s)}{N_1(s)}ds-A.

Полученное значение в общем случае зависит от выбранной точки x1x_1. Следовательно, различным точкам прямой y=ηy=\eta соответствуют различные значения CC.

Но одна интегральная кривая общего семейства должна определяться одним и тем же постоянным значением CC. Поэтому невозможно подобрать одно фиксированное значение CC, конечное или бесконечное, которое давало бы всю прямую y=η.y=\eta. Значит, постоянное решение yηy\equiv\eta не входит в общий интеграл ни при каком значении произвольной постоянной и является особым решением.

Аналогичные рассуждения проводятся для вертикальной интегральной кривой

xξ,N1(ξ)=0.x\equiv\xi, \qquad N_1(\xi)=0.

При этом x=ξx=\xi рассматривается как решение уравнения, перевёрнутого к исходному:

dxdyN(x,y)M(x,y).\frac{dx}{dy} \frac{N(x,y)}{M(x,y)}.
Built with LogoFlowershow