2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Разложение решения матричного дифура
Сообщение21.10.2015, 19:20 
Аватара пользователя
Здравствуйте!

Дано:
$A, Y, X, Z, B$ - матрицы $2\times 2$
$Y' = A(t)Y, Y(0)=Y_0$
$A(t)=\alpha A_1(t)+ \beta A_2(t)$
Дальше время опускаю.
Представим решение в виде произведения $Y=XZ$, где $X $- решение $X'=BX$, $B=\alpha A_1$.

Подставим в исходное уравнение:
$(XZ)' = A \cdot XZ$
$X'Z + XZ' = A\cdot XZ$
$BXZ + XZ' = AXZ$
$XZ' = AXZ - BXZ$
$Z' = X^{-1}(AX-BX)Z$

В итоге у меня уравнение для $Z'$:
$Z' = (X^{-1}AX - X^{-1}BX)Z$,
а в статье, где это выводится, результат
$Z' = (X^{-1}AX - B)Z$.

В чём я ошибся?

 
 
 
 Re: Разложение решения матричного дифура
Сообщение21.10.2015, 23:56 
Аватара пользователя
А нет ли у Вас начального условия для $X'=BX$? У Вас правильно, а в статье используется допущение что $BX=XB$

 
 
 
 Re: Разложение решения матричного дифура
Сообщение22.10.2015, 00:41 
Аватара пользователя
Начальное условие есть: $X(0)=I$
В частном случае, который рассматривается, $B=\alpha(t) A = \alpha(t) \cdot \operatorname{diag}(1,-1)$, поэтому решение для $X$ будет экспоненциалом с экспонентами на диагонали.
В итоге в $X^{-1}BX$ будет произведение трёх диагональных матриц, в котором экспоненты красиво дадут 1, оставив только $B$.

Спасибо за помощь!

 
 
 
 Re: Разложение решения матричного дифура
Сообщение22.10.2015, 02:57 
Аватара пользователя
Никакого частного случая: $X=e^{Bt}$ коммутирует с $B$

 
 
 [ Сообщений: 4 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group