Линейные многомерные модели — различия между версиями

Материал из synset
Перейти к: навигация, поиск
Строка 158: Строка 158:
 
Если в пределе <math>\textstyle t\to\infty</math> у системы существует стационарный режим, то в этом случае матрица дисперсий <math>\textstyle \mathbf{D}</math> становится постоянной, а ковариация зависит только от разности времён <math>\textstyle \tau</math>.
 
Если в пределе <math>\textstyle t\to\infty</math> у системы существует стационарный режим, то в этом случае матрица дисперсий <math>\textstyle \mathbf{D}</math> становится постоянной, а ковариация зависит только от разности времён <math>\textstyle \tau</math>.
  
<math>\textstyle \bullet</math> Таким образом, алгоритм решения линейной задачи следующий:\\ <math>\textstyle \triangleright</math> Находим собственные значения и вектора матрицы <math>\textstyle \mathbf{A}</math>.\\ <math>\textstyle \triangleright</math> Записываем решение для средних () и выражаем <math>\textstyle \mu_k</math> через <math>\textstyle \mathbf{x}_0</math>.\\ <math>\textstyle \triangleright</math> При помощи соотношения <math>\textstyle \left[e^{\mathbf{A} t}\right]_{\alpha\beta}=\partial \overline{x}_\alpha/x_{0\beta}</math> находим <math>\textstyle e^{\mathbf{A} t}</math>.\\ <math>\textstyle \triangleright</math> Вычисляем матрицу дисперсий <math>\textstyle D_{\alpha\beta}</math>. \\
+
<math>\textstyle \bullet</math> Таким образом, алгоритм решения линейной задачи следующий:
 +
 
 +
* <math>\textstyle \triangleright</math> Находим собственные значения и вектора матрицы <math>\textstyle \mathbf{A}</math>.\\ <math>\textstyle \triangleright</math> Записываем решение для средних () и выражаем <math>\textstyle \mu_k</math> через <math>\textstyle \mathbf{x}_0</math>.
 +
 
 +
* <math>\textstyle \triangleright</math> При помощи соотношения <math>\textstyle \left[e^{\mathbf{A} t}\right]_{\alpha\beta}=\partial \overline{x}_\alpha/x_{0\beta}</math> находим <math>\textstyle e^{\mathbf{A} t}</math>.
 +
 
 +
* <math>\textstyle \triangleright</math> Вычисляем матрицу дисперсий <math>\textstyle D_{\alpha\beta}</math>. \\
  
 
----
 
----

Версия 19:45, 21 февраля 2010

Уравнение стохастического осциллятора << Оглавление >> Многомерие помогает одномерию


Найдём решение линейных стохастических уравнений (по — сумма):

Постоянный вектор можно убрать сдвигом . В решении делается обратный сдвиг. Поэтому будем изучать однородное уравнение, которое запишем в матричной форме:

где и — не зависящие от и времени матрицы.

Для определения среднего проще всего сразу воспользоваться соотношением ():

(EQN)

где — вектор начального значения. Если мы хотим "вернуть" , то потребуются две замены: и .

Монотонная зависимость от в матричной записи решения () обманчива. Рассмотрим стохастический осциллятор из предыдущего раздела:

(EQN)

В этом случае матрицу можно разбить на сумму двух матриц:

Невозможно разобрать выражение (синтаксическая ошибка): {\displaystyle \mathbf{A } = \begin{pmatrix} -\lambda & -\omega \\ \omega & -\lambda \\ \end{pmatrix} = \omega \cdot \mathbf{q} -\lambda \cdot \mathbf{1},\;\;\;\;\;\;где\;\;\;\; \mathbf{1}= \begin{pmatrix} 1 & 0 \\ 0 & 1 \\ \end{pmatrix}, \;\;\;\;\;\; \mathbf{q}= \begin{pmatrix} 0 & -1 \\ 1 & 0 \\ \end{pmatrix}.}

Несложно проверить, что:

Так как матрицы и коммутируют друг с другом (), экспонента суммы разбивается на произведение . Раскладывая второй множитель по степеням и учитывая аналогичное разложение для синуса и косинуса, решение можно представить в следующем виде:

(EQN)

Оно же выше было получено другим методом. Таким образом, "монотонная" зависимость от времени в матричных соотношениях вполне может превратиться в периодическую функцию.

Найдём более практичное, чем (), представление для решения линейного уравнения. Будем его искать в виде:

(EQN)

Постоянный вектор является собственным вектором матрицы , а параметр ""— её собственным значением. Перенося в левую часть, получаем систему однородных уравнений относительно , которая имеет ненулевое решение, только если её детерминант равен нулю:

Это уравнение называется характеристическим и является полиномом -той степени по . Обычно оно имеет различных решений . Часть из них может оказаться комплексными. Для каждого из них мы решаем уравнение () и находим собственные вектора . Внимание! Верхний индекс — это номер собственного вектора, а не его компонента.

Теперь общее решение для среднего значения вектора переменных состояния можно записать в следующем виде:

(EQN)

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

Если матрица симметрична, то собственные вектора можно выбрать ортогональными: (звёздочка — комплексное сопряжение). В этом случае .

Когда выражены через , можно найти явное представление экспоненты от матрицы. Действительно, из (), взяв производную по компонентам начального условия, имеем . В частности, если собственные вектора ортогональны (), то:

(EQN)

В качестве упражнения ( H) предлагается найти для матрицы 2x2, у которой . Необходимо это сделать прямым разложением экспоненты при помощи собственных значений.

Выразим теперь решение стохастической линейной системы через гауссовы переменные. Введём новый вектор , удовлетворяющий, по лемме Ито (), стр. \pageref{Ito_N_dim_mat}, следующему уравнению:

Матрица зависит только от времени, поэтому решение этого уравнения легко найти при помощи итерационного метода:

Сумма независимых гауссовых чисел снова пропорциональна гауссовому числу, которое удобно представить в виде суммы независимых величин (второе равенство). Найдём значения . Для этого вычислим среднее от :

Учитывая независимость случайных величин и , а также переходя к непрерывному пределу , получаем ():

или:

(EQN)

Напомню, что (см. стр. \pageref{math_mat_tensor}). Решение для запишем в матричном виде, учитывая, что при :

Поэтому, так как , окончательное решение системы линейных стохастических уравнений имеет вид:

(EQN)

где . Вектор представляет собой набор независимых случайных чисел с гауссовым распределением, имеющим нулевое среднее и единичную дисперсию, а — среднее значение (), (). В качестве упражнения ( H) предлагается найти матрицу для двухмерного осциллятора и проверить решение ().

Вычислим матрицу дисперсий:

Учитывая (), имеем:

(EQN)

Это соотношение можно ( H) сразу получить из уравнения для средних (), стр. \pageref{aver_F_n_dim_for_xx}, из которых следует матричное уравнение:

(EQN)

Если существует стационарный режим, то и уравнение () позволяет легко найти .

Распределение для имеет гауссовый вид, поэтому, зная матрицу дисперсий, можно записать марковскую плотность вероятности:

где — обратная матрица дисперсий и — средние значения динамических переменных. Они полностью определяют свойства процесса. В частности, характеристическая функция ( H):

позволяет легко находить моменты произвольных порядков.

При помощи (), () несложно ( H) найти ковариационную матрицу:

(EQN)

Если в пределе у системы существует стационарный режим, то в этом случае матрица дисперсий становится постоянной, а ковариация зависит только от разности времён .

Таким образом, алгоритм решения линейной задачи следующий:

  • Находим собственные значения и вектора матрицы .\\ Записываем решение для средних () и выражаем через .
  • При помощи соотношения находим .
  • Вычисляем матрицу дисперсий . \\

Уравнение стохастического осциллятора << Оглавление >> Многомерие помогает одномерию

Стохастический мир - простое введение в стохастические дифференциальные уравнения