Фундаментальная матрица Коши

Материал из sawiki
Перейти к навигации Перейти к поиску

Фундаментальная матрица системы линейных однородных дифференциальных уравнений — матрица, столбцы которой образуют фундаментальную систему решений этой системы.

Определение

Фундаментальная матрица Коши $$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)$$ - сопряжённая система.

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

Рассмотрим неоднородную систему: \[ x'(t) = A(t)x(t)+B(t), \quad x(t)\in\mathbb{R}^n, \] где $$B(t)\in\mathbb{R}^n$$ — локально суммируемая функция.

В случае $$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} \]

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

Рассмотрим систему: \begin{equation} x'(t) = A(t)x(t), \quad x(t)\in\mathbb{R}^n. \end{equation}

Здесь $$A(t)$$ — периодическая с периодом $$T>0$$.
$$A(t+T)=A(t) \quad \forall t.$$
$$\tilde x[t]=x[t-kT]$$ — траектория (2). \[ \overset{\cdot}{\tilde x}[t]=\dot x[t-kT]=A(t-kT)x[t-kT]=A(t)\tilde x[t], \quad \forall k\in\mathbb{Z}. \] $$x[t]=X(t,t_0)x^0$$. Пусть $$t_0=0$$. Тогда $$x[t]=X(t,0)x^0$$. \[ X(t+kT,kT)x^0=\tilde x[t+kT]=x[t]=X(t,0)=x^0,\\ \tilde x[kT]=x[0]=x^0. \]

\[ X(t,0)=X(t+kT,kT)=X(t+kT,0)X(0,kT). \] Тогда

$$X(t+kT,0)=X(t,0)X(kT,0)=X(t,0)(X(T,0))^k.$$

Переобозначим $$\Phi(t)=X(t,0)$$. Тогда $$\Phi(T)=X(T,0)$$ — матрица монодромии.

$$\Phi(t+kT)=\Phi(t)(\Phi(T))^k.$$


Теорема Флоке

Если $$A(t+T)\equiv A(t),\quad \forall t.$$
Тогда существует $$\bar A=const\in\mathbb{C}^{n\times m}$$ и $$\psi(t)$$ такие, что $$\psi(t+T)=\psi(t), \quad \forall t,$$ $$|\psi(t)|\not=0$$, $$\psi(t)\in\mathbb{R}^{n\times n}$$ и $$\Phi(t)=\psi(t)e^{\bar At}.$$
Доказательство:
$$\Phi(0)=I=\psi(0).$$ Тогда $$\psi(T)=\psi(0)=I.$$
$$\Phi(T)=\psi(T)e^{\bar At}$$. Отсюда \begin{equation} e^{\bar At}=\Phi(T). \end{equation} Решением $$(3)$$ является $$\bar A=\frac1{T}\ln\Phi(T),\quad |\Phi(T)|\not=0$$. \begin{equation} \psi(t)=\Phi(t)\left(e^{\bar At}\right)^{-1}=\Phi(t)e^{-\bar At}. \end{equation} $$\psi(t+T)=\Phi(t+T)e^{-\bar A(t-T)}= \Phi(t)\underbrace{\Phi(T)e^{-\bar AT}}_{=I \text{ в силу } (3)}e^{-\bar At} = \Phi(t)e^{-\bar At}=\psi(t).$$ $$\blacksquare$$

Нахождение

Один из возможных способов нахождения фундаментальной матрицы можно найти на вкладке "Приложение преобразования Лапласа"

Пример применения

Формула Коши для траектории линейной системы \[ \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. \]