Рівняння Ріккаті

Рівняння Ріккаті — нелінійне звичайне диференціальне рівняння виду:

Рівняння Ріккаті
Названо на честь	Якопо Ріккаті
Формула	$y'(x)=q_{0}(x)+q_{1}(x)\,y(x)+q_{2}(x)\,y^{2}(x)$
Підтримується Вікіпроєктом	Вікіпедія:Проєкт:Математика

$y'=q_{0}(x)+q_{1}(x)\,y+q_{2}(x)\,y^{2},\quad x\in I,\quad (*)$

де $q_{0}(x),q_{1}(x),q_{2}(x)$ — неперервні на деякому інтервалі $I$ функції. У випадку $q_{2}(x)\equiv 0$ рівняння Ріккаті є лінійним неоднорідним рівнянням першого порядку, а у випадку $q_{1}(x)\equiv 0$ — диференціальним рівнянням Бернуллі. В обох цих часткових випадках рівняння Ріккаті легко інтегрується. У загальному, рівняння Ріккаті не розв'язується у квадратурах.

Частковий випадок рівняння $(*)$ вигляду

$y'=ay^{2}+bx^{\alpha },\quad (**)$

де $a,b,\alpha \,$ — сталі, досліджував Якопо Ріккаті і яке називають спеціальним рівнянням Ріккаті.

Рівняння $(**)$ ще цікаве насамперед з огляду на такий факт. Даніель Бернуллі близько 1725 року встановив, що спеціальне рівняння Ріккаті допускає відшукання загального розв'язку в елементарних функціях, якщо $\alpha =-2\,$ або $\alpha ={\frac {-4n}{(2n-1)}}$, де $n$ — ціле число. У 1841 році Жозеф Ліувілль з'ясував, що при всіх інших значеннях $\alpha$ це рівняння вже не можна зінтегрувати в елементарних функціях. Його загальний розв'язок виражається через циліндичні функції.

Рівняння та його узагальнення на випадок систем диференціальних рівнянь мають важливі застосування в багатьох математичних дисциплінах.

Мотивація

Розглянемо таку задачу: описати закон руху тіла маси $m$, яке вільно падає під дією сили тяжіння, враховуючи опір повітря.

Нехай в деякий початковий момент часу $t_{0}$ тіло знаходиться на висоті $H$.

Шлях, який пройде тіло падаючи $\tau$ секунд (від моменту часу $t_{0}$, якому відповідає $\tau =0$) позначимо $S(\tau )$. Відповідно, в момент часу $t>t_{0}$ тіло буде знаходитиметься на висоті $H-S(t-t_{0})=H-S(\tau )$, де $\tau =t-t_{0}$.

Швидкість та прискорення тіла в момент $\tau$ позначимо через $V(\tau )$ та $a(\tau )$ відповідно.

На тіло, яке вільно падає, діють дві сили: сила земного тяжіння $F_{t}$ та сила опору повітря $F_{op}$. Як відомо, $F_{t}=mg$, де $m$ — маса тіла, а $g$ — прискорення вільного падіння. Сила опору повітря $F_{op}$ прямопропорційна квадрату швидкості тіла, тобто $F_{op}=k\cdot V^{2}(\tau )$, де $k$ — коефіцієнт пропорційності, який залежить від форми тіла, площі поперечного перерізу тіла та густини повітря.

Згідно другого закону Ньютона отримуємо, що

$ma(\tau )=F_{t}-F_{op}.$

Враховуючи, що $a(\tau )=V^{\prime }(\tau )$ та вирази для $F_{t},F_{op}$ одержуємо

$mV^{\prime }(\tau )=mg-kV^{2}(\tau ).$

Отримане рівняння поділимо на $m$ та перепишемо так:

$V^{\prime }(\tau )=-{\frac {k}{m}}V^{2}(\tau )+g.$

Це спеціальне рівняння Ріккаті $(**)$ у якому $a=-{\frac {k}{m}}$, $\alpha =0$, $b=g$. Розв'язавши отримане рівняння знаходимо $V(\tau )$. Тоді пройдений тілом шлях знаходимо за формулою $S(\tau )=\int _{0}^{\tau }V(\tau )d\tau$.

Властивості

Рівняння Ріккаті та його розв'язки володіють такими властивостями.

• Рівняння Ріккаті завжди можна зінтегрувати в квадратурах, якщо вдалося знайти хоча б один його частинний розв'язок.

Справді, якщо $y_{1}$ — частинний розв'язок рівняння, то заміна змінних $y=y_{1}+u$, де $u$ — нова невідома функція незалежної змінної $x$, зводить це рівняння до рівняння Бернуллі. Підставивши в $(*)$ $y=y_{1}+u$, дістанемо

$y_{1}'+u'=q_{0}+q_{1}\cdot (y_{1}+u)+q_{2}\cdot (y_{1}+u)^{2},$

Розкривши дужки та врахувавши, що $y_{1}'=q_{0}+q_{1}\,y_{1}+q_{2}\,y_{1}^{2}$ одержимо таке рівняння для нової функції $u$ :

$u'-(q_{1}+2\,q_{2}\,y_{1})\,u=q_{2}\,u^{2},$

яке є диференціальним рівнянням Бернуллі.

• Якщо $y_{1}(x),y_{2}(x)$ деякі часткові розв'язки рівняння Ріккаті, то загальний розв'язок визначається за формулою:

$y={\frac {Cy_{1}+U(x)y_{2}}{C+U(x)}},\quad U(x)=\exp \left[\int q_{2}(y_{1}-y_{2})dx\right].$

• Загальний розв'язок рівняння Ріккаті $(*)$ є раціональною функцією від сталої інтегрування, і навпаки, будь-яке диференціальне рівняння першого порядку, що володіє цією властивістю, є рівнянням Ріккаті.
• Якщо $y_{1}(x),\ldots ,y_{4}(x)$ — часткові розв'язки рівняння Ріккаті, що відповідають значенням $c_{1},\ldots ,c_{4}$ сталої інтегрування, то має місце така тотожність:

${\frac {y_{3}(x)-y_{1}(x)}{y_{3}(x)-y_{2}(x)}}:{\frac {y_{4}(x)-y_{1}(x)}{y_{4}(x)-y_{2}(x)}}\equiv {\frac {c_{3}-c_{1}}{c_{3}-c_{2}}}:{\frac {c_{4}-c_{1}}{c_{4}-c_{2}}}.$

З даної формули випливає, що можна побудувати загальний розв'язок рівняння Ріккаті, якщо відомо три його часткові розв'язки.

Зведення до лінійного рівняння другого порядку

Нелінійне рівняння Ріккаті завжди може бути зведене до лінійного диференціального рівняння другого порядку.

Розглянемо рівняння Ріккаті

$y'=q_{0}(x)+q_{1}(x)y+q_{2}(x)y^{2}\!$

за умови $q_{2}(x)\neq 0$ (Якщо $q_{2}(x)=0\,$, то рівняння Ріккаті стає лінійним рівнянням першого порядку і може бути легко розв'язане). Функція, визначена формулою $v=yq_{2}\,$, де $y$ — розв'язок вихідного рівняння Ріккаті, задовольняє рівняння Ріккаті виду:

$v'=v^{2}+R(x)v+S(x),\!$

де $S=q_{2}q_{0}$ і $R=q_{1}+\left({\frac {q_{2}'}{q_{2}}}\right)$. Справді

$v'=(yq_{2})'=y'q_{2}+yq_{2}'=(q_{0}+q_{1}y+q_{2}y^{2})q_{2}+v{\frac {q_{2}'}{q_{2}}}=q_{0}q_{2}+\left(q_{1}+{\frac {q_{2}'}{q_{2}}}\right)v+v^{2}.\!$

Виконавши заміну $v=-u'/u\,$, одержуємо що функція $u$ задовольняє лінійне рівняння другого порядку

$u''-R(x)u'+S(x)u=0\!$

оскільки

$v'=-(u'/u)'=-(u''/u)+(u'/u)^{2}=-(u''/u)+v^{2},$

а тому

$u''/u=v^{2}-v'=-S-Rv=-S+Ru'/u\!$

і остаточно

$u''-Ru'+Su=0.\!$

Розв'язавши отримане рівняння другого порядку відносно $u$ за формулою $y=-u'/(q_{2}u)\,$ одержуємо розв'язок вихідного рівняння Ріккаті.

Приклади

Сім'я розв'язків рівняння Ріккаті $y'=3e^{x}-y-e^{-x}y^{2}$

Приклад 1. Розв'язати рівняння Ріккаті

$y'=3e^{x}-y-e^{-x}y^{2}$

У даному рівнянні в нас $q_{0}(x)=3e^{x}$, $q_{1}(x)=-1$, $q_{2}(x)=-e^{-x}$.

Розв'язати рівняння Ріккаті можна спираючись на відомі властивості його розв'язку. Якщо нам вдасться підібрати один чи два часткові розв'язки, то зможемо потім записати і його загальний розв'язок, зокрема звівши до рівняння Бернуллі. З вигляду коефіцієнтів рівняння спробуємо це зробити шукаючи розв'язок у формі $y=Ae^{mx}$, де $A,m$ -- сталі, які потрібно знайти. Підставимо у рівняння:

$Ame^{mx}=3e^{x}-Ae^{mx}-e^{-x}A^{2}e^{2mx}\quad \Rightarrow$

$e^{mx}(Am+A^{2}e^{(m-1)x}+A)=3e^{x}.$

Прирівнявши показники експонент (поклавши $m=1$) та скоротивши дві частини рівності на $e^{x}$ отримуємо

$2A+A^{2}=3.$

Розв'язавши отримане квадратне рівняння щодо $A$ одержуємо два можливі значення числа $A$: $A=1$ та $A=-3$.

Отже, нам вдалося знайти два часткові розв'язки розглядуваного рівняння:

$y_{1}=-3e^{x},\quad y_{2}=e^{x}.$

Маючи два розв'язки, можемо записати загальний. Спочатку знаходимо

$U(x)=\exp \left[\int q_{2}(y_{1}-y_{2})dx\right]=\exp \left[\int -e^{-x}(-4e^{x})dx\right]=\exp \left[\int 4dx\right]=e^{4x}$

Тоді шуканий загальний розв'язок записується у вигляді

$y={\frac {Cy_{1}+U(x)y_{2}}{C+U(x)}}={\frac {-3Ce^{x}+e^{5x}}{C+e^{4x}}}.$

Матричне рівняння Ріккаті

Матричним рівнянням Ріккаті називається диференціальне рівняння:

${\frac {dX}{dt}}=XA(t)X+B(t)X+XC(t)+D(t).$

де $X$ — невідома матриця розмірів $n\times m$, а розміри матриць $A(t),B(t),C(t),D(t)$ відповідно $m\times n$, $n\times n$, $m\times m$, $n\times m$.

Матричне рівняння Ріккаті відіграє важливу роль в теорії лінійних гамільтонових систем, варіаційному численні, задачах оптимального управління, фільтрації, стабілізації керованих лінійних систем. Наприклад, оптимальне управління $u_{0}$ в задачі мінімізації функціонала:

$x^{T}(t_{1})\Phi x(t_{1})=\int _{t_{0}}^{t_{1}}\left[x^{T}(t)M(t)x(t)+u^{T}(t)N(t)x(t)\right]dt$

на розв'язках системи:

$x'=A(t)x+B(t)u,\quad x(t_{0})=x_{0},$

($n\times n$-матриці $\Phi ,M(t)$ симетричні і невід'ємноозначені, а $m\times m$ -матриця $N(t)$, додатноозначена при $t\in [t_{0},t_{1}]$), визначається формулою:

$u_{0}(t)=-N^{-1}(t)B^{T}(t)Z(t)x\,$

де $Z(t)$— розв'язок матричного рівняння Ріккаті:

$Z'=-ZA(t)-A^{T}(t)Z+ZB(t)N^{-1}(t)B^{T}(t)Z-M(t)\,$

з граничною умовою $Z(t_{1})=\Phi .\,$

У задачах управління на нескінченному інтервалі часу важливими є питання про існування у матричного рівняння Ріккаті невід'ємноозначеного обмеженого на $[t_{0},\infty ]$ розв'язку, про існування періодичного або майже періодичного розв'язку (у випадку періодичних або майже періодичних коефіцієнтів рівняння) і про способи наближеної побудови таких рішень.

Див. також

• Диференціальне рівняння Бернуллі
• Лінійне диференційне рівняння

Посилання

• Wolfram Math. World: Riccati Differential Equation

Література

• В. Ф. Зайцев, А. Д. Полянин. Справочник по обыкновенным дифференциальным уравнениям. М.: Физматлит, 2001.
• Лауфер М. Я. О решении уравнений Риккати // Лауфер М. Я. Избранные задачи математической физики. Сб. статей.— Северодвинск: НТО кораблестроителей им. акад. А. Н. Крылова, Севмашвтуз, Северодв. отд-ние Ломоносов. фонда, 2005.— стр. 137—140.— ISBN 5-7723-0605-9.
• Лионс Ж.-Л., Оптимальное управление Системами, описываемыми уравнениями с частными производными, пер. с франц., М., 1972;
• Самойленко А. М., Перестюк М. О., Парасюк I.О. Диференціальні рівняння: Підручник. — К.: Либідь, 2003 р. — 600 с.