Фундаментальная матрица Коши
Фундаментальная матрица системы линейных однородных дифференциальных уравнений - матрица, столбцы которой образуют фундаментальную систему решений этой системы.
Содержание
Определение
Фундаментальная матрица Коши $$X(t,\tau)$$ - решение задачи Коши \begin{equation} \left\{ \begin{aligned} & \frac{\partial X(t,\tau)}{\partial t} = A(t)X(t,\tau),\\ & X(\tau,\tau) = I. \end{aligned} \right. \end{equation}
$$X(t,\tau) = [x^1(t,\tau),\ldots,x^n(t,\tau)]$$, где $$x^j$$ - решение \[ \left\{ \begin{aligned} & \dot x(t) = A(t)x(t),\\ & x(t_0)=x^0. \end{aligned} \right. \]
Свойства
- Полугрупповое свойство: $$X(t,\tau) = X(t,s)X(s,\tau)$$.
$$\left\{ \begin{aligned} & \dot x(t)=A(t)x(t),\\ &x(\tau)=\xi. \end{aligned} \right.$$
$$\qquad \Downarrow$$
$$x(t)=X(t,\tau)\xi \Rightarrow \begin{aligned} & x(t) = X(t,s)x(s)\\ & x(s)=X(s,\tau)\xi \end{aligned}$$
|
$$X(t,\tau)\xi = X(t,s)X(s,\tau)\xi \quad \forall \xi.$$ |
- $$\tau = t$$
$$I=X(t,s)X(s,t)$$
|
$$X(s,t)=X^{-1}(t,s).$$ |
- Рассмотрим $$I=X(t,\tau)X(\tau,t)$$.
Продифференцируем по $$\tau$$: $$0=\frac{\partial X(t,\tau)}{\partial\tau}X(\tau,t)+X(t,\tau)\frac{\partial X(\tau,t)}{\partial\tau}$$. При этом $$\frac{\partial X(\tau,t)}{\partial\tau} = A(\tau)X(\tau,t)$$.
Отсюда следует, что
|
$$\frac{\partial X(t,\tau)}{\partial\tau}=-X(t,\tau)A(\tau)$$. |
- $$S(t,\tau)=X^T(\tau,t)$$
$$\frac{\partial S(t,\tau)}{\partial t}=-A^T(t)S(t,\tau)$$
$$\frac{\partial X(\tau,t)}{\partial t}=-X(\tau,t)A(t)$$
$$\Rightarrow $$
|
$$S(t)=-A^T(t)S(t)$$ - сопряжённая система. |
Система дифференциальных уравнений с постоянными коэффициентами
В случае $$A(t) = A = const$$ матрицант равен $$X(t) = e^{A(t-t_0)},$$ где $$e^{As}$$ - матричная экспонента (см. статью про матричный экспоненциал), следовательно, матрица Коши:
$$C(t,s) = e^{A(t-s)},$$
$$C(t,s) = X(t-s).$$
Таким образом, в этом случае для получения матрицы Коши достаточно подставить $$(t-s)$$ в качестве аргумента матрицанта.
Общее решение системы линейных неоднородных дифференциальных уравнений с постоянными коэффициентами имеет вид:
\[
x(t)=e^{A(t-t_0)}x(t_0)+\int\limits_{t_0}^{t}e^{A(t-s)}B(s)ds.
\]
Пример
Рассмотрим систему $$(1)$$ при $$A = \begin{bmatrix} -2 & 12 \\-1 & 5 \end{bmatrix}$$:
$$
\left\{
\begin{aligned}
& \dot x_1(t) = -2x_1(t) + 12 x_2(t),\\
& \dot x_2(t) = -x_1(t) + 5 x_2(t) + au,\\
& x_1(0) = x_1^0,\\
& x_2(0) = x_2^0,
\end{aligned}
\right.
$$
где $$u$$ - управление, $$a, x_1^0, x_2^0$$ - некоторые константы.
Найдём матричный экспоненциал.
1. Найдём собственные значения и собственные векторы матрицы $$A$$.
\[
\begin{vmatrix}
A - \lambda I
\end{vmatrix}
=
\begin{vmatrix}
-2 -\lambda & 12 \\ -1 & 5-\lambda
\end{vmatrix}
= \lambda^2-3\lambda+2 \Rightarrow \lambda_1=1, \quad \lambda_2=2.
\]
При $$\lambda = 1$$ собственный вектор $$(4,1)^T$$. При $$\lambda = 2$$ собственный вектор $$(3,1)^T$$.
2. Выпишем ФСР.
$$ x(t) = c_1e^{2t} \begin{bmatrix} 3 \\ 1 \end{bmatrix} + c_2e^t \begin{bmatrix} 4 \\ 1 \end{bmatrix}$$.
Найдём константы $$c_1$$ и $$c_2$$.
$$\left\{ \begin{aligned} 3c_1+4c_2=1, \\ c_1+c_2=0. \end{aligned} \right. \Rightarrow \left\{ \begin{aligned}& c_1=-1, \\& c_2=1. \end{aligned} \right.$$
$$\qquad$$
$$\left\{ \begin{aligned} 3c_1+4c_2=0, \\ c_1+c_2=1. \end{aligned} \right. \Rightarrow \left\{ \begin{aligned}& c_1=4, \\& c_2=-3. \end{aligned} \right.$$
\[ \text{Тогда } e^{At} = \begin{bmatrix} -3e^{2t}+4e^t & 12e^{2t}-12e^t \\ -e^{2t}+e^t & 4e^{2t}-3e^t \end{bmatrix}. \]
Далее по формуле Коши найдём $$x_1$$ и $$x_2$$:
$$x(t) = e^{At}x(0) + \int\limits_0^{t} e^{(t-\tau)}B(\tau)d\tau.$$ Здесь $$B = \begin{bmatrix} 0 \\ a \end{bmatrix}.$$
Вычислим $$e^{(t-\tau)}B(\tau):$$ $$\begin{bmatrix} a(12e^{2(t-\tau)}-12e^{(t-\tau)}) \\ a(4e^{2(t-\tau)}-3e^{(t-\tau)}) \end{bmatrix}.$$
Проинтегрируем $$\int\limits_0^{t} e^{(t-\tau)}B(\tau)d\tau = \begin{bmatrix} 6a \\ a \end{bmatrix}.$$
Итого получим
$$
\begin{aligned}
& x_1 = x_1^0(-3e^{2t}+4e^t) + x_2^0(12e^{2t}-12e^t)+6a,\\
& x_2 = x_2^0(-e^{2t}+e^t) + x_2^0(4e^{2t}-3e^t)+a.
\end{aligned}
$$
Нахождение
Один из возможных способов нахождения фундаментальной матрицы можно найти на вкладке "Приложение преобразования Лапласа"
Пример применения
Формула Коши для траектории линейной системы \[ \left\{ \begin{aligned} & \dot x(t) = A(t)x(t) + B(t)u(t)+c(t),\\ & x(t_0)=x^0. \end{aligned} \right. \] выглядит следующим образом: \[ x(t)=X(t,t_0)x^0+\int\limits_{t_0}^{t} X(t,\tau)[B(\tau)u(\tau)+c(\tau)]d\tau. \]
