Формула Лиувиля - Liouvilles formula - Wikipedia

В математика, Формула Лиувилля, также известное как тождество Абеля-Якоби-Лиувилля, представляет собой уравнение, которое выражает детерминант из квадратная матрица решение системы первого порядка однородных линейные дифференциальные уравнения через сумму диагональных коэффициентов системы. Формула названа в честь Французский математик Джозеф Лиувиль. Формула Якоби обеспечивает другое представление того же математического отношения.

Формула Лиувилля является обобщением Личность Авеля и может использоваться, чтобы это доказать. Поскольку формула Лиувилля связывает различные линейно независимый решения системы дифференциальных уравнений, это может помочь найти одно решение из другого (ых), см. пример приложения ниже.

Формулировка формулы Лиувилля

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

${\displaystyle y'=A(x)y}$

на интервал я из реальная линия, куда А(Икс) за Икс ∈ я обозначает квадратную матрицу размерности п с настоящий или же сложный записи. Позволять Φ обозначим матричнозначное решение на я, что означает, что каждый Φ (Икс) квадратная матрица размерности п с реальными или сложными записями и производная удовлетворяет

${\displaystyle \Phi '(x)=A(x)\Phi (x),\qquad x\in I.}$

Позволять

${\displaystyle \mathrm {tr} \,A(\xi )=\sum _{i=1}^{n}a_{i,i}(\xi ),\qquad \xi \in I,}$

обозначить след из А(ξ) = (а_я, j(ξ))_{я, j ∈ {1,...,п}}, сумма его диагональных элементов. Если след А это непрерывная функция, то определитель Φ удовлетворяет

${\displaystyle \det \Phi (x)=\det \Phi (x_{0})\,\exp {\biggl (}\int _{x_{0}}^{x}\mathrm {tr} \,A(\xi )\,{\textrm {d}}\xi {\biggr )}}$

для всех Икс и Икс₀ в я.

Пример приложения

Этот пример показывает, как формула Лиувилля может помочь найти общее решение системы однородных линейных дифференциальных уравнений первого порядка. Учитывать

${\displaystyle y'=\underbrace {\begin{pmatrix}1&-1/x\\1+x&-1\end{pmatrix}} _{=\,A(x)}y}$

на открытом интервале я = (0, ∞). Предположим, что простое решение

${\displaystyle y(x)={\begin{pmatrix}1\\x\end{pmatrix}},\qquad x\in I,}$

уже найден. Позволять

${\displaystyle y(x)={\begin{pmatrix}y_{1}(x)\\y_{2}(x)\end{pmatrix}}}$

обозначим другое решение, тогда

${\displaystyle \Phi (x)={\begin{pmatrix}y_{1}(x)&1\\y_{2}(x)&x\end{pmatrix}},\qquad x\in I,}$

является квадратным матричным решением указанного выше дифференциального уравнения. Поскольку след А(Икс) равен нулю для всех Икс ∈ я, Из формулы Лиувилля следует, что определитель

${\displaystyle c_{1}:=\det \Phi (x)=x\,y_{1}(x)-y_{2}(x),\qquad x\in I,}$

(1)

на самом деле постоянная, не зависящая от Икс. Записывая первую составляющую дифференциального уравнения для у, получаем, используя (1) который

${\displaystyle y'_{1}(x)=y_{1}(x)-{\frac {y_{2}(x)}{x}}={\frac {x\,y_{1}(x)-y_{2}(x)}{x}}={\frac {c_{1}}{x}},\qquad x\in I.}$

Следовательно, интегрированием мы видим, что

${\displaystyle y_{1}(x)=c_{1}\ln x+c_{2},\qquad x\in I,}$

с участием натуральный логарифм и постоянная интеграции c₂. Решение уравнения (1) за у₂(Икс) и заменяя у₁(Икс) дает

${\displaystyle y_{2}(x)=x\,y_{1}(x)-c_{1}=\,c_{1}x\ln x+c_{2}x-c_{1},\qquad x\in I,}$

что является общим решением для у. Со специальным выбором c₁ = 0 и c₂ = 1 мы восстанавливаем простое решение, с которого начали, выбор c₁ = 1 и c₂ = 0 дает линейно независимое решение. Следовательно,

${\displaystyle \Phi (x)={\begin{pmatrix}\ln x&1\\x\ln x-1&x\end{pmatrix}},\qquad x\in I,}$

это так называемое фундаментальное решение системы.

Доказательство формулы Лиувилля

Мы опускаем аргумент Икс для краткости. Посредством Формула Лейбница для определителей, производная определителя Φ = (Φ_я, j)_{я, j ∈ {0,...,п}} можно вычислить, дифференцируя одну строку за раз и взяв сумму, т.е.

${\displaystyle (\det \Phi )'=\sum _{i=1}^{n}\det {\begin{pmatrix}\Phi _{1,1}&\Phi _{1,2}&\cdots &\Phi _{1,n}\\\vdots &\vdots &&\vdots \\\Phi '_{i,1}&\Phi '_{i,2}&\cdots &\Phi '_{i,n}\\\vdots &\vdots &&\vdots \\\Phi _{n,1}&\Phi _{n,2}&\cdots &\Phi _{n,n}\end{pmatrix}}.}$

(2)

Поскольку матричнозначное решение Φ удовлетворяет уравнению Φ '= АΦ, для каждого элемента матрицы Φ '

${\displaystyle \Phi '_{i,k}=\sum _{j=1}^{n}a_{i,j}\Phi _{j,k}\,,\qquad i,k\in \{1,\ldots ,n\},}$

или для всей строки

${\displaystyle (\Phi '_{i,1},\dots ,\Phi '_{i,n})=\sum _{j=1}^{n}a_{i,j}(\Phi _{j,1},\ldots ,\Phi _{j,n}),\qquad i\in \{1,\ldots ,n\}.}$

Когда мы вычитаем из я ^th ряд линейной комбинации

${\displaystyle \sum _{\scriptstyle j=1 \atop \scriptstyle j\not =i}^{n}a_{i,j}(\Phi _{j,1},\ldots ,\Phi _{j,n}),}$

всех остальных строк, то значение определителя не изменится, следовательно,

${\displaystyle \det {\begin{pmatrix}\Phi _{1,1}&\Phi _{1,2}&\cdots &\Phi _{1,n}\\\vdots &\vdots &&\vdots \\\Phi '_{i,1}&\Phi '_{i,2}&\cdots &\Phi '_{i,n}\\\vdots &\vdots &&\vdots \\\Phi _{n,1}&\Phi _{n,2}&\cdots &\Phi _{n,n}\end{pmatrix}}=\det {\begin{pmatrix}\Phi _{1,1}&\Phi _{1,2}&\cdots &\Phi _{1,n}\\\vdots &\vdots &&\vdots \\a_{i,i}\Phi _{i,1}&a_{i,i}\Phi _{i,2}&\cdots &a_{i,i}\Phi _{i,n}\\\vdots &\vdots &&\vdots \\\Phi _{n,1}&\Phi _{n,2}&\cdots &\Phi _{n,n}\end{pmatrix}}=a_{i,i}\det \Phi }$

для каждого я ∈ {1, . . . , п} линейностью определителя по каждой строке. Следовательно

${\displaystyle (\det \Phi )'=\sum _{i=1}^{n}a_{i,i}\det \Phi =\mathrm {tr} \,A\,\det \Phi }$

(3)

к (2) и определение следа. Осталось показать, что из этого представления производной следует формула Лиувилля.

Исправить Икс₀ ∈ я. Поскольку след А предполагается непрерывной функцией на я, она ограничена на каждом замкнутом и ограниченном подынтервале я и поэтому интегрируема, следовательно

${\displaystyle g(x):=\det \Phi (x)\exp \left(-\int _{x_{0}}^{x}\mathrm {tr} \,A(\xi )\,{\textrm {d}}\xi \right),\qquad x\in I,}$

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

${\displaystyle g'(x)={\bigl (}(\det \Phi (x))'-\det \Phi (x)\,\mathrm {tr} \,A(x){\bigr )}\exp {\biggl (}-\int _{x_{0}}^{x}\mathrm {tr} \,A(\xi )\,{\textrm {d}}\xi {\biggr )}=0,\qquad x\in I,}$

из-за производной в (3). Следовательно, грамм должен быть постоянным на я, так как иначе мы получили бы противоречие с теорема о среднем значении (применяется отдельно к действительной и мнимой части в комплексном случае). С грамм(Икс₀) = det Φ (Икс₀), Формула Лиувилля следует из решения определения грамм за det Φ (Икс).

Рекомендации

• Чиконе, Кармен (2006), Обыкновенные дифференциальные уравнения с приложениями. (2-е изд.), Нью-Йорк: Springer-Verlag, стр. 152–153, ISBN  978-0-387-30769-5, МИСТЕР  2224508, Zbl  1120.34001
• Тешл, Джеральд (2012), Обыкновенные дифференциальные уравнения и динамические системы., Провиденс: Американское математическое общество, МИСТЕР  2961944, Zbl  1263.34002