Формула Лиувилля

Толкование Перевод

Формула Лиувилля: Формула Лиуви́лля-Острогра́дского — формула, связывающая определитель Вронского (вронскиа́н) для решений дифференциального уравнения и коэффициенты в этом уравнении.

Пусть есть дифференциальное уравнение вида

$y^{(n)}+P_1(x)y^{(n-1)}+P_2(x)y^{(n-2)}+...+P_n(x)y=0,$

тогда $W(x)=W(x_0)e^{-\int_{x_0}^x P_1(\zeta)d\zeta}=Ce^{-\int P_1(x)dx},$ где $W(x)$ — определитель Вронского

Для линейной однородной системы дифференциальных уравнений

$y'(x) = A(x) y(x),$ где $A(x)$ — непрерывная квадратная матрица порядка $n$ , справедлива формула Лиувилля-Остроградского

$W(x)=W(x_0)e^{\int_{x_0}^x \mathop{\rm tr} A(\zeta)d\zeta},$ где $\mathop{\rm tr} A(x)$ — след матрицы $A(x)$

Содержание

1 Правило дифференцирования определителя размерности 2

2 Правило дифференцирования определителя размерности

3 Доказательство для уравнения второго порядка

4 Доказательство для линейной системы обыкновенных дифференциальных уравнений

5 Доказательство для линейного дифференциального уравнения произвольного порядка

6 Применение формулы Лиувилля-Остроградского

7 Пример

7.1 Используемая литература

Правило дифференцирования определителя размерности 2

Производная определителя $\Delta=\Delta(x)=\begin{vmatrix} a_{11}(x) & a_{12}(x) \\ a_{21}(x) & a_{22}(x) \end{vmatrix}=a_{11}a_{22}-a_{12}a_{21}$ по переменной х имеет вид $\frac{d \Delta}{d x}=(a_{11}a_{22}-a_{12}a_{21})'=a_{11}'a_{22}+a_{11}a_{22}'-a_{12}'a_{21}-a_{12}a_{21}'=\begin{vmatrix} a_{11}' & a_{12}' \\ a_{21} & a_{22} \end{vmatrix}+\begin{vmatrix} a_{11} & a_{12} \\ a_{21}' & a_{22}' \end{vmatrix}$

Правило дифференцирования определителя размерности $n$

Пусть $\Delta = \Delta(x) = \det \left( \begin{matrix} a_{11}(x) & a_{12}(x) & \dots & a_{1n}(x)\\ a_{21}(x) & a_{22}(x) & \dots & a_{2n}(x)\\ \vdots & \vdots & \ddots & \vdots\\ a_{n1}(x) & a_{n2}(x) & \dots & a_{nn}(x)\\ \end{matrix} \right)$

Тогда для производной $\Delta'(x)$ верно

$\Delta'(x) = \begin{vmatrix} a_{11}'(x) & a_{12}'(x) & \dots & a_{1n}'(x)\\ a_{21}(x) & a_{22}(x) & \dots & a_{2n}(x)\\ \vdots & \vdots & \ddots & \vdots\\ a_{n1}(x) & a_{n2}(x) & \dots & a_{nn}(x)\\ \end{vmatrix} + \begin{vmatrix} a_{11}(x) & a_{12}(x) & \dots & a_{1n}(x)\\ a_{21}'(x) & a_{22}'(x) & \dots & a_{2n}'(x)\\ \vdots & \vdots & \ddots & \vdots\\ a_{n1}(x) & a_{n2}(x) & \dots & a_{nn}(x)\\ \end{vmatrix} + \dots + \begin{vmatrix} a_{11}(x) & a_{12}(x) & \dots & a_{1n}(x)\\ a_{21}(x) & a_{22}(x) & \dots & a_{2n}(x)\\ \vdots & \vdots & \ddots & \vdots\\ a_{n1}'(x) & a_{n2}'(x) & \dots & a_{nn}'(x)\\ \end{vmatrix}$

(в $i$ -м слагаемом продифференцирована $i$ -я строка)

Доказательство

Воспользуемся формулой полного разложения определителя

$\Delta(x) = \sum_{i_1, i_2, \dots, i_n} (-1)^{P(i_1, i_2, \dots, i_n)} a_{1 i_1}(x) a_{2 i_2}(x) \cdots a_{n i_n}(x)$

Сумма взята по всевозможным перестановкам чисел $1,2,\dots,n$ , $P(\cdot)$ - четность перестановки.

Дифференцируя это выражение по $x$ , получим

$\begin{align}[l] \Delta'(x) &= \sum_{i_1, i_2, \dots, i_n} (-1)^{P(i_1, i_2, \dots, i_n)} \frac{d\left(a_{1 i_1}(x) a_{2 i_2}(x) \cdots a_{n i_n}(x)\right)}{dx} = \\ &=\sum_{i_1, i_2, \dots, i_n} (-1)^{P(i_1, i_2, \dots, i_n)} \left( a_{1 i_1}'(x) a_{2 i_2}(x) \cdots a_{n i_n}(x) + \dots + a_{1 i_1}(x) a_{2 i_2}(x) \cdots a_{n i_n}'(x) \right) = \\ &=\sum_{i_1, i_2, \dots, i_n} (-1)^{P(i_1, i_2, \dots, i_n)} a_{1 i_1}'(x) a_{2 i_2}(x) \cdots a_{n i_n}(x) + \\ &+\sum_{i_1, i_2, \dots, i_n} (-1)^{P(i_1, i_2, \dots, i_n)} a_{1 i_1}(x) a_{2 i_2}'(x) \cdots a_{n i_n}(x) + \\ &+\dots + \\ &+\sum_{i_1, i_2, \dots, i_n} (-1)^{P(i_1, i_2, \dots, i_n)} a_{1 i_1}(x) a_{2 i_2}(x) \cdots a_{n i_n}'(x) \end{align}$

В каждой сумме продифференцированы элементы $i$ -й строки и только они. Заменив суммы определителями, получим

$\Delta'(x) = \begin{vmatrix} a_{11}'(x) & a_{12}'(x) & \dots & a_{1n}'(x)\\ a_{21}(x) & a_{22}(x) & \dots & a_{2n}(x)\\ \vdots & \vdots & \ddots & \vdots\\ a_{n1}(x) & a_{n2}(x) & \dots & a_{nn}(x)\\ \end{vmatrix} + \begin{vmatrix} a_{11}(x) & a_{12}(x) & \dots & a_{1n}(x)\\ a_{21}'(x) & a_{22}'(x) & \dots & a_{2n}'(x)\\ \vdots & \vdots & \ddots & \vdots\\ a_{n1}(x) & a_{n2}(x) & \dots & a_{nn}(x)\\ \end{vmatrix} + \dots + \begin{vmatrix} a_{11}(x) & a_{12}(x) & \dots & a_{1n}(x)\\ a_{21}(x) & a_{22}(x) & \dots & a_{2n}(x)\\ \vdots & \vdots & \ddots & \vdots\\ a_{n1}'(x) & a_{n2}'(x) & \dots & a_{nn}'(x)\\ \end{vmatrix}$

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

Пусть в уравнении $y''+p(x)y'+q(x)y=0$ функции $p(x),q(x)$ непрерывны на $[a;b]$ , а

$y_1=y_1(x),y_2=y_2(x)$ — решения данного уравнения.

Продифференцировав определитель Вронского получим

$\frac{d W}{d x}=\frac{d}{dx}\begin{vmatrix} y_{1} & y_{2} \\ y_{1}' & y_{2}' \end{vmatrix}=\begin{vmatrix} y_{1}' & y_{2}' \\ y_{1}' & y_{2}' \end{vmatrix}+\begin{vmatrix} y_{1} & y_{2} \\ y_{1}'' & y_{2}'' \end{vmatrix}$

Первое слагаемое равно 0, так как этот определитель содержит 2 одинаковые строки. Подставив

$y_1''=-py_1'-qy_1$

$y_2''=-py_2'-qy_2$

во второе слагаемое получим

$\frac{d W}{d x}=\begin{vmatrix} y_{1} & y_{2} \\ -py_1'-qy_1 & -py_2'-qy_2 \end{vmatrix}$

Прибавив первую строку, домноженную на q, ко второй, получим

$\frac{d W}{d x}=\begin{vmatrix} y_{1} & y_{2} \\ -py_1' & -py_2'\end{vmatrix}=-pW$

решения линейно независимы, поэтому

$W \ne 0 \to \frac{d W}{W}=-p dx$ — дифференциальное уравнение с разделяющимися переменными.

Интегрируя, получим

$\ln |W|=-\int P(x)dx +\ln|C| \to \ln\left|\frac{W}{C}\right|=-\int P(x)dx \to W=C e^{-\int P(x)dx}$

Доказательство для линейной системы обыкновенных дифференциальных уравнений

Пусть вектор-функции ${\mathbf y}_1(x), {\mathbf y}_2(x), \dots, {\mathbf y}_n(x)$ — решения линейной системы ОДУ. Введем матрицу $\Phi$ следующим образом

$\Phi(x) = || \begin{matrix} {\mathbf y}_1(x)& {\mathbf y}_2(x)& \dots& {\mathbf y}_n(x) \end{matrix} ||$

Тогда $W(x) \equiv \det \Phi(x)$ . Воспользуемся тем, что $y_i(x)$ - решения системы ОДУ, то есть ${\mathbf y}_i'(x) = A(x){\mathbf y}_i(x)$ .

В матричном виде последнее представимо в виде $|| \begin{matrix} {\mathbf y}_1'(x)& {\mathbf y}_2'(x)& \dots& {\mathbf y}_n'(x) \end{matrix} || = || \begin{matrix} A(x){\mathbf y}_1(x)& A(x){\mathbf y}_2(x)& \dots& A(x){\mathbf y}_n(x) \end{matrix} || = A(x) \Phi(x)$

или вводя производную от матрицы как матрицу из производных каждого элемента

$\Phi'(x) = A(x) \Phi(x)\,$

Пусть $\!\varphi_i(x)$ - $i$ -я строка матрицы $\Phi(x)\!$ . Тогда

$\varphi_i'(x) = \sum_{j=1}^n a_{ij}(x) \varphi_j(x)\,$

Последнее означает, что производная от $i$ -й строки матрицы $\Phi(x)\!$ есть линейная комбинация всех строк этой матрицы с коэффициентами из $i$ -й строки матрицы $A(x)\!$ . Рассмотрим определитель матрицы $\Phi(x)\!$ , в которой $i$ -я строка продифференцирована. Определитель не изменится, если из $i$ -й строки этой матрицы вычесть линейную комбинацию всех остальных строк.

$\left| \begin{matrix} \varphi_1(x)\\ \varphi_2(x)\\ \vdots\\ \varphi_i'(x)\\ \vdots\\ \varphi_n(x)\\ \end{matrix} \right| = \left| \begin{matrix} \varphi_1(x)\\ \varphi_2(x)\\ \vdots\\ \sum_{j=1}^n a_{ij}(x)\varphi_j(x)\\ \vdots\\ \varphi_n(x)\\ \end{matrix} \right| = \left| \begin{matrix} \varphi_1(x)\\ \varphi_2(x)\\ \vdots\\ \sum_{j=1}^n a_{ij}(x)\varphi_j(x) -\sum_{j\neq i} a_{ij}(x)\varphi_j(x) \\ \vdots\\ \varphi_n(x)\\ \end{matrix} \right| = \left| \begin{matrix} \varphi_1(x)\\ \varphi_2(x)\\ \vdots\\ a_{ii}(x) \varphi_i(x)\\ \vdots\\ \varphi_n(x)\\ \end{matrix} \right| = a_{ii}(x) W(x)$

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

$W'(x) = a_{11}(x) W(x) + a_{22}(x) W(x) + \dots + a_{nn}(x)W(x) = \mathop{\rm{tr}}A(x)W(x)$

Последнее обыкновенное дифференциальное уравнение имеет решение

$W(x) = W(x_0) e^{\int_{x_0}^x\mathop{\rm{tr}}A(\zeta) d\zeta}$

Доказательство для линейного дифференциального уравнения произвольного порядка

Линейное дифференциальное уравнение $n$ -го порядка

$\!y^{(n)}(x) + P_1(x)y^{(n-1)}(x) + \dots + P_{n-1}(x) y'(x) + P_n(x) y(x) = 0$

эквивалентно следующей системе

$\begin{align} y_{n-1}'(x) &= -P_1(x) y_{n-1}(x) - \dots - P_{n-1}(x) y_1(x) - P_n(x) y_0(x)\\ y_{n-2}'(x) &= y_{n-1}\\ \vdots\\ y_{1}'(x) &= y_{2}\\ y_{0}'(x) &= y_{1}\\ \end{align}$

с матрицей $A(x)\!$ следующего вида

$A(x)=\left( \begin{matrix} 0 & 1 & 0 & \dots & 0\\ 0 & 0 & 1 & \dots & 0\\ 0 & 0 & \ddots & \ddots & 0\\ 0 & 0 & \dots & 0 & 1\\ -P_n(x) & -P_{n-1}(x) & \dots & -P_2(x) & -P_1(x) \\ \end{matrix} \right)$

Вронскианы исходного уравнения и системы совпадают, а след матрицы $A(x)\!$ равен $-P_1(x)\!$ . Подстановкой в формулу для системы получаем

$W(x) = W(x_0) e^{-\int_{x_0}^x P_1(\zeta) d\zeta}$

Применение формулы Лиувилля-Остроградского

Пусть известно решение $y_1(x)$ линейного обыкновенного дифференциального уравнения второго порядка. Используя формулу Лиувилля-Остроградского возможно найти линейно независимое от него решение $y_2(x)$ той же системы.

Распишем вронскиан:

$Ce^{-\int P(x)dx}=y_1y_2'-y_1'y_2=W.$

$\frac{W}{y_1^2}=\frac{y_1y_2'-y_1'y_2}{y_1^2}= \left(\frac{y_2}{y_1}\right)',$ поэтому

$\frac{y_2}{y_1}=\int \frac{W}{y_1^2} dx +B$ $\to y_2=y_1\left(\int \frac{W}{y_1^2} dx +B\right)= y_1\int \frac{C e^{-\int P(x)dx}}{y_1^2}dx+B y_1$

Так как для линейной независимости $y_1(x)$ и $y_2(x)$ достаточно $W \neq 0$ , приняв $C=1,\, B=0,$ получим $y_2=y_1\int \frac{ e^{-\int P(x)dx}}{y_1^2}dx.$

Пример

Пусть в уравнении $y''-\mathop{ \rm{tg}} x\, y'+2y=0$ известно частное решение $\!y_1=\sin x$ . Воспользовавшись формулой Лиувилля-Остроградского, получим

$y_2=\sin x\int \frac{dx}{\sin^2 x e^{-\int \tan x dx}}=\sin x \, \ln \left|\tan x +\frac{1}{\cos x}\right| \, -1.$

Тогда общее решение однородного уравнения $y=C_1\left(\sin x \,\ln \left|\tan x + \frac{1}{\cos x}\right| \, -1\right)+C_2 \sin x$

Используемая литература

Агафонов С. А.,Герман А. Д.,Муратова Т. В. Дифференциальные уравнения. Учебник для вузов -М. Изд-во МГТУ им. Баумана, 1999.-336с (Серия Математика в техническом университете ;Вып. VIII),Глава 5 параграф 2 .

Романко В. К. Курс дифференциальных уравнений и вариационного исчисления. - 2-е изд. - М.: Лаборатория Базовых Знаний, 2001 - 344 с.

Для улучшения этой статьи желательно^?:

Исправить статью согласно стилистическим правилам Википедии.

Проставив сноски, внести более точные указания на источники.

Проверить статью на грамматические и орфографические ошибки.

Категория:
Дифференциальные уравнения

Игры ⚽ Поможем сделать НИР

Полезное

Смотреть что такое "Формула Лиувилля" в других словарях:

Формула Лиувилля-Остроградского — Формула Лиувилля Остроградского формула, связывающая определитель Вронского (вронскиан) для решений дифференциального уравнения и коэффициенты в этом уравнении. Пусть есть дифференциальное уравнение вида y(n) + P1(x)y(n − 1) + P2(x)y(n − 2) … Википедия
Формула Родригеса — представляет собой: В геометрии формула поворота Родригеса Формула для получения серии выражений, повторяя дифференцирование какой то другой функции. Типичное применение: составление серии из ортогональных многочленов. Конкретней, для функции ω,… … Википедия
Формула Родрига — представляет собой: В геометрии формула поворота Родрига Формула для получения серии выражений, повторяя дифференцирование какой то другой функции. Типичное применение: составление серии из ортогональных многочленов. Конкретней, для функции , , и … Википедия
ЛИУВИЛЛЯ - ОСТРОГРАДСКОГО ФОРМУЛА — Л и у в и л л я формула, соотношение, связывающее вронскиан системы решений и коэффициенты линейного обыкновенного дифференциального уравнения. Пусть x1(t), . . ., xn(t) произвольная система прешений линейной однородной системы п го порядка с… … Математическая энциклопедия
Теорема Лиувилля — Теоремы, названные в честь Жозефа Лиувилля: Теорема Лиувилля об ограниченных целых аналитических функциях Теорема Лиувилля о конформных отображениях Теорема Лиувилля о приближении алгебраических чисел Теорема Лиувилля о сохранении фазового объёма … Википедия
ШТУРМА - ЛИУВИЛЛЯ ЗАДАЧА — задача, порожденная на конечном или бесконечном интервале ( а, b) изменения переменной хуравнением и нек рыми граничными условиями, где р(х) и r(х) положительны, l(х)действительна, а комплексный параметр. Начало глубокому изучению этой задачи… … Математическая энциклопедия
ОСТРОГРАДСКОГО - ЛИУВИЛЛЯ ФОРМУЛА — см. Лиувилля Остроградского формула … Математическая энциклопедия
Теорема Лиувилля об ограниченных целых аналитических функциях — У этого термина существуют и другие значения, см. Теорема Лиувилля. Теорема Лиувилля об ограниченных целых аналитических функциях: если целая функция комплексных переменных ограничена, то есть то есть константа. Обобщения Если … Википедия
Интегральная формула Коши — В этой статье отсутствует вступление. Пожалуйста, допишите вводную секцию, кратко раскрывающую тему статьи. Содержание … Википедия
ВРОНСКИАН — определитель Вроньского, определитель системы пвектор функций размерности п имеющий вид: В. системы n скалярных функций имеющих производные до ( п 1) го порядка включительно, есть определитель Это понятие было введено Ю. Вроньским [1]. Если… … Математическая энциклопедия

Словари и энциклопедии на Академике

Формула Лиувилля

Содержание

Правило дифференцирования определителя размерности 2

Правило дифференцирования определителя размерности $n$

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

Доказательство для линейной системы обыкновенных дифференциальных уравнений

Доказательство для линейного дифференциального уравнения произвольного порядка

Применение формулы Лиувилля-Остроградского

Пример

Используемая литература

Полезное

Смотреть что такое "Формула Лиувилля" в других словарях:

Поделиться ссылкой на выделенное

Словари и энциклопедии на Академике

Википедия

Формула Лиувилля

Содержание

Правило дифференцирования определителя размерности 2

Правило дифференцирования определителя размерности

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

Доказательство для линейной системы обыкновенных дифференциальных уравнений

Доказательство для линейного дифференциального уравнения произвольного порядка

Применение формулы Лиувилля-Остроградского

Пример

Используемая литература

Полезное

Смотреть что такое "Формула Лиувилля" в других словарях:

Поделиться ссылкой на выделенное

Прямая ссылка:

Правило дифференцирования определителя размерности $n$