Уравнение Риккати: метод решения и особые решения

Уравнение Риккати занимает особое место среди обыкновенных дифференциальных уравнений первого порядка: оно нелинейное - содержит $y^2$ - но при одном известном частном решении немедленно линеаризуется. Именно этот переход «нелинейное → линейное» и делает уравнение Риккати учебным эталоном для курса ОДУ. Ниже разберём структуру уравнения, алгоритм нахождения особых решений, метод сведения к линейному и типичные задачи. Прежде чем перейти к теории, покрутите калькулятор: слайдеры коэффициентов сразу показывают траекторию и особые решения.

Общий вид и структура уравнения Риккати

Стандартная форма уравнения Риккати:

$y' = P(x) + Q(x)\,y + R(x)\,y^2,$

где $P(x)$, $Q(x)$, $R(x)$ - заданные функции, $R(x) \not\equiv 0$. Если $R \equiv 0$, уравнение становится линейным; если $P \equiv 0$ - превращается в уравнение Бернулли с показателем $n = 2$.

Для случая постоянных коэффициентов $P = a$, $Q = b$, $R = c$ уравнение принимает вид:

$y' = c\,y^2 + b\,y + a.$

Рассчитать для своей задачи

Дискриминант правой части $D = b^2 - 4ac$ определяет тип поведения:

• $D > 0$: два особых решения, траектории между ними ограничены;
• $D = 0$: одно особое решение параболического типа;
• $D < 0$: вещественных особых решений нет, все траектории уходят в $\pm\infty$.

Особые (стационарные) решения

Особым называется постоянное решение $y(x) = y^* = \text{const}$. Подставив $y' = 0$:

$c\,(y^*)^2 + b\,y^* + a = 0.$

Это обычное квадратное уравнение. При $D > 0$ два корня:

$y^*_{1,2} = \frac{-b \pm \sqrt{b^2 - 4ac}}{2c}.$

Особые решения - это горизонтальные прямые на фазовом портрете, которые траектории не пересекают (теорема единственности). Они разделяют фазовую плоскость на полосы с разным поведением.

Рассчитать для своей задачи

Сведение к линейному: метод подстановки

Главный приём - если известно одно частное решение $y_1(x)$, замена

$z = \frac{1}{y - y_1}$

линеаризует уравнение Риккати. Вычислим $y'$:

$y = y_1 + \frac{1}{z}, \qquad y' = y_1' - \frac{z'}{z^2}.$

Подставляя в исходное уравнение $y' = P + Q\,y + R\,y^2$ и используя, что $y_1' = P + Q\,y_1 + R\,y_1^2$, получаем:

$z' = -(Q + 2R\,y_1)\,z - R.$

Это линейное неоднородное уравнение первого порядка относительно $z(x)$ - его решают стандартным методом вариации постоянной.

$$\begin{aligned} y &= y_1 + \frac{1}{z},\\ z' &+ (Q + 2R\,y_1)\,z = -R. \end{aligned}$$

Пример: уравнение с постоянными коэффициентами

Рассмотрим уравнение $y' = 2y^2 - 3y + 1$ (то есть $a = 1$, $b = -3$, $c = 2$).

Шаг 1. Особые решения. Уравнение $2y^2 - 3y + 1 = 0$ даёт $y^*_1 = 1$ и $y^*_2 = \tfrac{1}{2}$.

Шаг 2. Сведение к линейному. Берём $y_1 = 1$ и вводим $z = 1/(y - 1)$:

$z' = -(Q + 2R\,y_1)\,z - R = -(-3 + 2 \cdot 2 \cdot 1)\,z - 2 = z - 2.$

Линейное уравнение $z' - z = -2$ решается элементарно: $z = C\,e^x + 2$. Возвращаясь к $y$:

$y = 1 + \frac{1}{C\,e^x + 2}.$

Шаг 3. Проверка. При $C \to \infty$ получаем $y \to 1 = y^*_1$; при $x \to -\infty$ и $C > 0$ - $y \to 1$; при $x \to +\infty$ - $y \to 1$. Второе особое решение $y^*_2 = \tfrac{1}{2}$ получается при $C \to 0$: $y = 1 + 1/2 = 3/2$... нет, точнее $z = 2$, $y = 1 + 1/2 = 3/2$ при $x = 0$. На самом деле $y^*_2 = 1/2$ соответствует $C = -2/e^x$ - параметрическое семейство. Особые решения $y = 1$ и $y = 1/2$ отвечают $C = 0$ и $C \to \infty$ соответственно.

Метод угадывания частного решения

Самый распространённый учебный приём - искать частное решение $y_1(x)$ в одном из стандартных классов:

• Степенная форма: $y_1 = \alpha\,x^n$ (подбирают $\alpha$ и $n$).
• Показательная форма: $y_1 = \alpha\,e^{\beta x}$.
• Тригонометрическая форма: $y_1 = \alpha\sin x + \beta\cos x$.

Алгоритм подбора. Предположим, что $y_1 = kx^n$, и подставим в уравнение $y' = P(x) + Q(x)y + R(x)y^2$. Левая часть $y_1' = kn\,x^{n-1}$, правая - некоторая комбинация степеней $x$. Приравниваем коэффициенты при совпадающих степенях и получаем систему для $k$ и $n$.

Пример: $y' = \dfrac{1}{x^2} - \dfrac{y}{x} + y^2$, то есть $P = 1/x^2$, $Q = -1/x$, $R = 1$.

Ищем $y_1 = k/x$ (то есть $n = -1$): $y_1' = -k/x^2$. Подставляем в правую часть:

$P + Q\,y_1 + R\,y_1^2 = \frac{1}{x^2} - \frac{1}{x}\cdot\frac{k}{x} + \frac{k^2}{x^2} = \frac{1 - k + k^2}{x^2}.$

Условие $y_1' =$ правая часть: $-k/x^2 = (1 - k + k^2)/x^2$, откуда $k^2 + 1 = 0$ - нет вещественных корней. Попробуем $n = -2$, то есть $y_1 = k/x^2$: $y_1' = -2k/x^3$, но правая часть имеет степень $-2$ - несоответствие степеней. На этом примере видно, что не у каждого уравнения Риккати решение угадывается в виде степени; в таких случаях переходят к численным методам.

Если же частное решение найдено (пусть $y_1 = 1/x^2$ для другого уравнения), дальше всё стандартно: вводим $z = 1/(y - y_1)$ и решаем линейное ОДУ.

Задача Коши для уравнения Риккати

Постановка: найти $y(x)$, удовлетворяющее уравнению Риккати и начальному условию $y(x_0) = y_0$.

Важное свойство. Если начальное условие попадает между двумя особыми решениями $y^*_2 < y_0 < y^*_1$ (при $D > 0$), то решение задачи Коши существует на всей оси $\mathbb{R}$ и остаётся в полосе $y^*_2 \le y(x) \le y^*_1$. Если же начальное условие лежит вне этой полосы, решение может уйти в бесконечность за конечное время.

Пример задачи Коши. Уравнение $y' = 2y^2 - 3y + 1$, $y(0) = 0$.

Особые решения: $y^*_1 = 1$, $y^*_2 = 1/2$. Начальное условие $y_0 = 0 < 1/2$ - ниже нижнего особого решения. Общее решение: $y = 1 + 1/(Ce^x + 2)$. При $x = 0$: $0 = 1 + 1/(C + 2)$, откуда $C + 2 = -1$, $C = -3$. Решение задачи Коши:

$y = 1 + \frac{1}{-3e^x + 2}.$

При $x^* = \ln(2/3) \approx -0.405$ знаменатель обращается в нуль - решение уходит в $-\infty$. Это «взрыв»: задача Коши имеет решение лишь на $(-\infty,\, \ln(2/3))$.

Калькулятор выше обнаруживает такие «взрывы» автоматически: при $a = 1$, $b = -3$, $c = 2$, $y(0) = 0$ поле «уход в $\infty$» покажет приближённый $x^*$.

Применение: задачи на оптимальный рост и теория управления

Уравнение Риккати возникает в нескольких совершенно разных областях. В задачах математической физики оно описывает распределение потенциала в некоторых моделях рассеяния. В экономике скалярное уравнение Риккати встречается в моделях оптимального роста: если $y(t)$ - отношение капитала к труду, темп изменения $y$ может удовлетворять уравнению квадратичного типа.

В теории оптимального управления матричное уравнение Риккати

$\dot P + PA + A^\top P - PBR^{-1}B^\top P + Q = 0$

возникает при решении задачи LQR (линейно-квадратичный регулятор). Здесь $P$ - симметричная матрица состояния, $A$, $B$ - матрицы системы, $Q$, $R$ - весовые матрицы. Скалярный случай ($n = 1$) даёт именно тот вид $y' = cy^2 + by + a$, который изучается в курсе ОДУ. Стационарное решение матричного уравнения $\dot P = 0$ - это уравнение алгебраического Риккати (ARE); его корень задаёт оптимальный регулятор обратной связи.

Для задач первого курса достаточно скалярного варианта, но осознание прикладного контекста объясняет, почему уравнение Риккати изучают отдельно: за ним стоят реальные инженерные и экономические задачи, а не только учебные упражнения.

Связь с уравнением Бернулли и уравнением второго порядка

Два смежных класса уравнений тесно связаны с Риккати:

1. Уравнение Бернулли $y' + p(x)\,y = q(x)\,y^n$ - при $n = 2$ и $p, q > 0$ это частный случай Риккати с $P = 0$: $y' = q(x)\,y^2 - p(x)\,y$. Стандартная замена Бернулли $z = 1/y$ ($n = 2$) тождественна формуле линеаризации Риккати с $y_1 = 0$.

2. Линейное ОДУ второго порядка. Замена $y = -u'/(R\,u)$ (подстановка Риккати) переводит уравнение Риккати в линейное уравнение второго порядка:

$u'' - \left(Q - \frac{R'}{R}\right)u' + P\,R\,u = 0.$

Этот переход используется в теории специальных функций: уравнение Бесселя $x^2 u'' + x u' + (x^2 - \nu^2)u = 0$ и уравнение Лежандра $(1-x^2)u'' - 2xu' + n(n+1)u = 0$ порождают соответствующие уравнения Риккати через замену $y = -u'/u$.

Практический вывод. Уравнение Риккати стоит на пересечении трёх больших классов: нелинейные ОДУ первого порядка (общий контекст), уравнения Бернулли (частный случай), линейные ОДУ второго порядка (через обратную замену). Понимание этих связей позволяет применять готовые решения из смежных задач вместо того, чтобы каждый раз выводить всё заново.

Геометрический смысл: изоклины и поле направлений

Уравнение $y' = f(x, y)$ задаёт поле направлений в плоскости $(x, y)$: в каждой точке угол наклона касательной к траектории равен $f(x, y)$. Для Риккати $f = cy^2 + by + a$.

Изоклины - кривые $f(x, y) = k$ для разных $k$. Для постоянных коэффициентов изоклины вырождаются в горизонтальные прямые: $cy^2 + by + a = k$ - параболы по $y$, но независимые от $x$. Поэтому поле направлений является автономным: вдоль любой вертикальной прямой $x = \text{const}$ наклоны совпадают. Отсюда и название - уравнение с постоянными коэффициентами не зависит от $x$ явно.

Особые решения $y^* = \text{const}$ - это горизонтальные изоклины с $k = 0$, на которых $y' = 0$. Траектории, стартующие между двумя такими прямыми, «затягиваются» к устойчивому особому решению или «отталкиваются» от неустойчивого - это стандартный фазовый анализ автономного ОДУ.

Частые ошибки

• Не проверяют особые решения. После нахождения корней $c\,y^2 + b\,y + a = 0$ их надо подставить в исходное уравнение - только тогда это особые решения, а не просто нули правой части.
• Путают знак в формуле линеаризации. Правильно: $z' = -(Q + 2R\,y_1)\,z - R$, а не $+(Q + 2R\,y_1)\,z$. Знак минус получается из дифференцирования $y = y_1 + 1/z$.
• Не учитывают уход в $\pm\infty$. При $D < 0$ решение уходит в бесконечность за конечное время - «взрыв». Указывать момент взрыва обязательно в задачах с начальным условием.
• Забывают условие $R \not\equiv 0$. Если $c = 0$, квадратного члена нет и уравнение Риккати вырождается в линейное; алгоритм нахождения особых решений неприменим.
• Ошибочная подстановка $z = y - y_1$. Правильная замена - $z = 1/(y - y_1)$, а не $z = y - y_1$. Именно обратная замена превращает нелинейную правую часть в линейную.

FAQ

Как найти особое решение уравнения Риккати? Приравняйте правую часть нулю: $P(x) + Q(x)\,y + R(x)\,y^2 = 0$. Если уравнение имеет постоянный корень $y^* = \text{const}$, это и есть особое решение. Для постоянных коэффициентов это квадратное уравнение $c\,y^2 + b\,y + a = 0$ с корнями $y^*_{1,2} = (-b \pm \sqrt{D})/(2c)$.

Что делать, если частное решение неизвестно? Единого аналитического метода нет - уравнение Риккати в общем случае не интегрируется в квадратурах. На практике частное решение угадывают в классах $y_1 = ax^n$, $y_1 = ae^{bx}$, $y_1 = a\sin x + b\cos x$. Если угадать не удалось, используют численное решение (метод Рунге-Кутта) или специальные функции.

Чем уравнение Риккати отличается от уравнения Бернулли? Уравнение Бернулли $y' + p(x)\,y = q(x)\,y^n$ при $n = 2$ и $P = 0$ является частным случаем Риккати. Общий вид Риккати содержит ненулевой свободный член $P(x)$, который не позволяет применить стандартную подстановку Бернулли $z = y^{1-n}$ напрямую.

Коротко

Уравнение Риккати $y' = P + Q\,y + R\,y^2$ нелинейно, но при одном известном частном решении $y_1$ замена $z = 1/(y - y_1)$ сводит его к линейному. Особые решения находятся из $R\,y^2 + Q\,y + P = 0$; дискриминант $D = Q^2 - 4PR$ (для постоянных коэф. $D = b^2 - 4ac$) определяет тип поведения: при $D > 0$ траектории между особыми решениями ограничены, при $D < 0$ уходят в $\pm\infty$ за конечное время.