[МУЗЫКА]
[МУЗЫКА] Итак,
для нашей матричной экспоненты мы получили следующее выражение:
e в степени st = Q в -1,
здесь e в степени t √a² + d², 0, 0.
e в степени −t √a² + d².
Это наше собственное значение.
Здесь Q.
И для Q выражение
(1 √a² + d²
− d/a запятая.
Здесь нам все равно, где минус писать.
Просто можно умножить собственный вектор на произвольное число для удобства.
Просто чтобы минус перед длинным выражением не писать, я напишу у единицы.
Это будет одно и то же.
Всё, что я ни напишу в матрице Q,
какой множитель, он перекачается,
скомпенсируется соответственно матрицей Q в − 1,
которая находится стандартным образом,
который мы с вами знаем.
[БЕЗ_ЗВУКА]
[БЕЗ_ЗВУКА] Просто
элементарно по минорам можно идти к Q в −1.
Итак, в этом смысле оператор эволюции найден, чтобы решить эволюционную задачу,
нужно взять начальный вектор и помножить на него вот эту матрицу.
Это будет решение нашей задачи Коши, то есть задачи с начальными условиями.
Давайте все-таки получим некую качественную информацию.
Эта, как я уже говорил, система уравнений
не чистая какая-то фантазия формальная,
она описывает разбегания траекторий в гиперболическом течении.
Посмотрим, как ведет себя оператор эволюции в нашей матрице на больших
временах.
На больших временах эта экспонента пренебрежима.
Поэтому матрицы эволюции: Итак, t стремится к ∞.
Просто большое, не надо понимать, как бесконечная бесконечность.
Тогда [БЕЗ_ЗВУКА]
оператор эволюции имеет следующий простой вид.
Значит, вот экспонента убывает, ею пренебрегаем, то есть ставим здесь 0.
Эту экспоненту выносим на улицу.
[БЕЗ_ЗВУКА] И здесь
получаем, значит, наше Q в −1.
(1, 0, 0, 0) Q.
Я не буду явно теперь это произведение вычислять.
Единственное, что меня сейчас
будет интересовать, что, вообще говоря,
если посмотреть на наш результат, на эту диагональную матрицу.
Что можно сказать?
Что у нас есть одно направление, вдоль которого,
происходит растяжение и одно направление, вдоль которого, происходит сжатие.
Я тогда в предыдущей части не очень аккуратно нарисовал это
гиперболическое течение.
Давайте я его перерисую чуть более аккуратно.
Итак, у нас есть условное это направление-растяжение,
это направление-сжатие и течение устроено так.
Значит, вот
так вот.
То есть сюда у меня растягивается, а сюда сжимается.
Вот и поставим такой вопрос.
Если мы начальный вектор выбираем произвольно, случайно,
то с какой вероятностью мы получим экспоненциально
растущий на больших временах вектор?
Если совсем наивно рассуждать.
Есть два направления: одно растягивающее, другое — сжимающее.
Поэтому кажется, что половина направления растягивается, половина сжимается.
Но понятно, что это не так.
Для того чтобы ответить на этот вопрос, это то же самое,
что ответить на вопрос: с какой вероятностью вектор на
плоскости не имеет компоненты вдоль растягивающего направления?
Потому что даже если он имеют очень малую составляющую вдоль растягивающего
направления, эта составляющая будет экспоненциально
расти и на больших временах перебьет все остальное.
В этом смысле, если t велико, то найти вероятность
малого вектора, она экспоненциально мала.
Это нужно добиться того, чтобы компонента вектора вдоль
сжимающего направления, чтобы вектор практически весь лежал вдоль
сжимающего направления с экспоненциальной точностью.
Когда время t растет, то соответственно точность должна быть все больше.
И поэтому на больших временах с вероятностью 1 с точки зрения
выбора начального условия,
у нас будет экспоненциальный рост вектора x.
Если мы представляем себе этот вектор как начальное расстояние между двумя точками,
то с вероятностью 1 эти точки будут экспоненциально разбегаться.
Значит, такого рода поведение характерно
не только для такой гидродинамической постановки.
Можно говорить то, что экспоненциальное разбегание траектории
— это общее свойство так называемых стохастических систем.
В каком-то смысле с физической точки зрения это есть их определение.
Что если систему мы называем хаотической или стохастической,
это означает то, что возьмем два близких начальных условия,
выпустим фазовые траектории и эти фазовые траектории на больших временах будут
экспоненциально расходиться.
Хорошо, теперь рассмотрим совсем другой пример,
другую систему, в которой возникает некое новое качество.
Может быть даже более простая формальная система, чем эта.
Но качество в ней следующее, что матрица M,
которая была введена раньше, недиагонализуема.
Что значит недиагонализуема?
Я напоминаю элементарные сведения из линейной алгебры.
Матрица может быть приведена либо к диагональному виду то есть тогда после
преобразования подобия, после которого,
только на диагонали стоят любые элементы.
А вне диагонали стоят нули.
Либо матрица представляет собой блочную
диагональную и эти блоки являются жадановыми клетками.
Жаданова клетка возникает тогда,
когда есть некоторое вырождение, с точки зрения физики.
Конкретно вырождением мы называем совпадение собственных чисел.
Не то, чтобы всегда, но довольно часто.
И в физических системах это не такая большая экзотика.
Рассмотрим элементарный пример системы уравнений химической кинетики.
У нас есть три вещества.
Давайте я их обозначу: A, B и C.
Это концентрации реагентов.
Значит, вещество A распадается,
λ — это темп распада в единицу времени.
И распадается оно на B,
[БЕЗ_ЗВУКА] вещество B,
которое в свою очередь тоже распадается, сколько
получилось, с тем же самым темпом.
Темп — это λ, это количество событий в единицу времени.
Здесь количество распадов молекулы в единицу времени.
Здесь, соответственно, они распались, образовалось B.
Это приход.
Теперь B само распадается, пусть с тем же самым темпом.
И в конце концов,
продуктом распада B является вещество C.
Вот такая вот система уравнений.
Давайте запишем
ее в матричном виде.
[БЕЗ_ЗВУКА]
[БЕЗ_ЗВУКА]
[БЕЗ_ЗВУКА] Вот
такая вот матрица.
Все это умножить на вектор A, B, C.
Так, в следующей части
мы применим наш аппарат к анализу этой простой системы.
[МУЗЫКА]
[МУЗЫКА]
Фоминов Яков ВикторовичДоцент
Фейгельман Михаил ВикторовичГлавный научный сотрудник
Колоколов Игорь ВалентиновичПрофессор
Бурмистров Игорь СергеевичПрофессор
Скворцов Михаил АндреевичПрофессор
