Множество достижимости линейной управляемой системы, заданной при помощи ОДУ, без помехи. Внутренние оценки

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

Общий вид системы

Дана линейная система дифференциальных уравнений без помехи: \begin{equation} \label{1} \begin{cases} \dot x(t) = A(t)x(t) + B(t)u(t), \\ x(t) \in \mathcal{X}, \\ u(t) \in \mathcal{P}(t) \end{cases} \end{equation} Где \(\mathcal{P}(t)\) - непрерывное по Хаусдорфу многозначное отображение, \(\mathcal{P}(t) \subset conv\mathbb{R}^n; A(t) \in \mathbb{R}^{n \times n}, \ B(t) \in \mathbb{R}^{n \times m}, \ t \in [t_0, t_1], x \in \mathbb{R}^n, \ X \in \mathbb{R}^{n\times n}, \ q(t) \in \mathbb{R}^m, \ Q(t) \in \mathbb{R}^{m\times m} \). При этом многозначные отображения: \(\mathcal{X}\) и \(\mathcal{P}(t)\) - эллипсоиды: \[ \mathcal{X} = \mathcal{E}(x, X) \subset \mathbb{R}^n, \] \[ \mathcal{P}(t) = \mathcal{E}(q(t), Q(t)) \subset \mathbb{R}^m. \]

Эллипсоиды

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

Утверждение 1

  Выпуклое множество однозначно определяется своей  опорной функцией

Утверждение 2

  Эллипсоид с центром в точке q и матрицей Q в m-мерном вещественном пространстве можно определить как множество, построенное в направлении некоторого ненулевого вектора l:

\[\mathcal{E}(q, Q) = \{x \in \mathbb{R}^m \:|\: \langle x,\,l \rangle \leq \langle l,\,q \rangle + \langle l,\,Ql \rangle^\frac{1}{2}\},\]

Утверждение 3

 Для любого эллипсоида с центром в q, матрицей Q, и для любой матрицы A, такой что число ее столбцов совпадает с числом строк Q, выполнено тождество:

\[A\mathcal{E}(q,Q) = \mathcal{E}(Aq, AQA^{T}).\]

Теорема 1

 Для  суммы эллипсоидов можно получить внутреннюю эллипсоидальную оценку следующего вида:

\[ \sum\limits^n_{i=1} \mathcal{E}(q_i, Q_i) = \bigcup_{\parallel l \parallel = 1} \mathcal{E}_-(l), \] где \(\mathcal{E}_-(l) = \mathcal{E}(q_i, Q_-)\) — эллипсоид, построенный в направлении \(l\), с центром в \(q_i\) и матрицей \(Q_-(t) = Q'_*(t)Q_*(t)\), где \(Q_*(t) = \sum_{i=1}^n S_i(t) Q_i^\frac{1}{2}(t)\), а \(S_i(t)\) — некоторые ортогональные матрицы.

Внутренние эллипсоидальные оценки интеграла от многозначного отображения

Построим внутреннюю эллипсоидальную оценку интеграла от многозначного отображения, где в качестве образа выступает эллипсоид.

Получим эллипсоидальную оценку для суммы двух эллипсоидов.

Пусть даны эллипсоды: \(\mathcal{E_{1}}(0, Q_{1})\) и \(\mathcal{E_{2}}(0, Q_{2})\), где \(Q_{1}\) и \(Q_{2}\) - некоторые положительно определенные матрицы(случай невырожденных эллипсоидов).

Хотим найти \(\mathcal{E}\): \(\mathcal{E} \subseteq \mathcal{E_{1}}\) + \(\mathcal{E_{2}}\).

Введем обозначение: \( S[n] = (S_{1}, ... , S_{n}) \).

В качестве \(Q_-\) рассмотрим следующую матрицу: \(Q_-(S[2]) = (S_{1}Q^{1\2}_{1} + S_{2}Q^{1/2}_{2})^{'}(S_{1}Q^{1/2}+S_{2}Q^{1/2}_{2})\), где \(S_{1}\) и \(S_{2}\) - ортогональные матрицы: \(S'S\ = SS' = I\), \(I\) - единичная матрица. \(Q^{1/2}\) - квадратный корень из матрицы: у неотрицательно определенной матрицы \(Q\) \(\exists!\) квадратный корень, то есть симметричная матрица \(Q = Q^{1/2}Q^{1/2}\).

Покажем, что эллипсоид \(\mathcal{E_-}(0, Q_-) \) является внутренним по отношению к сумме \(\mathcal{E_{1}} + \mathcal{E_{2}}\), для этого оценим опорную функцию эллипсоида \(\mathcal{E_-}\).

\[ \mathcal{p^{2}}(l|\mathcal{E_-}(0, Q_-)) = \langle l, (S_{1}Q^{1/2}_{1} + S_{2}Q^{1/2}_{2})'(S_{1}Q^{1/2}_{1} + S_{2}Q^{1/2}_{2})l \rangle = \langle l, Q_{1}l \rangle + \langle l, Q_{2}l \rangle + \langle S_{1}Q^{1/2}_{1}l, S_{2}Q^{1/2}_{2}l \rangle \leq \] в силу неравенства Коши-Буняковского: \[ \leq \langle l, Q_{1}l \rangle + \langle l, Q_{2}l \rangle + 2 \langle l, Q_{1}l \rangle ^{1/2} \langle l,Q_{2}l \rangle ^{1/2} = (\langle l, Q_{1}l \rangle ^{1/2} + \langle l, Q_{2}l \rangle ^{1/2})^{2} = \mathcal{p^{2}}(l| \mathcal{E_{1}} + \mathcal{E_{2}}) \]

Заметим, что для \(l\) \(\neq\) \(0\) равенство достигается тогда и только тогда, когда аргументы скального произведения являются линейно зависимыми и сонаправленными, то есть: \(S_{1}Q^{1/2}_{1}l = \lambda S_{2}Q^{1/2}_{2}l \), где \(\lambda > 0\).

В силу того, что мы можем повернуть вектор \(Q^{1/2}l\) на любой угол, для любого \(l\) всегда можно выбрать такую матрицу \(S\), что в неравенстве для опорных функций будет выполняться равенство и, следовательно, существует точка, в которой будет происходить касание эллипсоидов в силу вложенности одного множества в другое.

В случае \(m\) эллипсоидов поступим по аналогии: \[ Q_-(s[m]) = (\sum_{k=1}^{m} S_{k}Q^{1/2}_{k})'(\sum_{k=1}^{m} S_{k}Q^{1/2}_{k}). \]

Внутренняя оценка для множества достижимости

Рассматривается следующая задача динамики: \begin{equation} \label{2} \begin{cases} \dot x = Ax + u, \\ x(t_0) \in \mathcal{E}_0(x_0, X_0), \\ u(t) \in \mathcal{E}(p(t), P(t)). \end{cases} \end{equation} Нам нужно найти внутреннюю оценку для множества достижимости, которое представимо в виде: \[ \mathcal{X[t]} = X(t,t_0)\mathcal{E_0}(x_0, X_0) + \int^{t}_{t_0}X(t, \tau)\mathcal{E}(p(\tau), \mathcal{P}(\tau)d{\tau} = \mathcal{E_0}(x_0, X(t, t_0)X_0X(t, t_0)^{T}) + \int^{t}_{t_0} \mathcal{E}(p(\tau), X(t, \tau)\mathcal{P}(\tau)X(t, \tau)^{T})d\tau \] Для системы \eqref{1} справедлива формула Коши: \[ x(t,t_0,x_0) = X(t,t_0)x_0 - \int^{t}_{t_0}X(t,\tau)B(\tau)u(\tau)d\tau, \] Где \(X(t,\tau)\) - фундаментальная матрица, удовлетворяющая системе: \[ \left\{\begin{aligned} & \frac{\partial X(t,\tau)}{\partial t} = A(t)X(t,\tau), \\ & X(\tau,\tau) = I. \end{aligned}\right. \] Тогда для матриц аппроксимирующего эллипсоида: \[ X\_ = (Q^{*})^{T}Q^{*}, где Q^{*} = (S_0(X_{0})^{1/2}(X(t, t_{0}))^{T} + \int^{t}_{t_0} \mathcal{S}(\tau)P^{1/2}(\tau)X(t, \tau)^{T}d\tau) \] Касание достигается в случае: \begin{equation} \label{3} \mathcal{S}(\tau)P^{1/2}(\tau)X(t, \tau)^{T}l(t) = \lambda(\tau)\mathcal{S_0}(X_{0})^{1/2}X(t, t_0)^{T}l(t), где \lambda(\tau) > 0, \forall \tau \in [t_0,t] \end{equation} В полученном равенстве у нас есть зависимость от t и \(\tau\): \[ \mathcal{S}(\tau) = \mathcal{S_t}(\tau), \lambda(\tau) = \lambda_{t}(\tau). \] Следовательно, мы не можем для фиксированного \(l\) построить хорошей оценки, так как для каждого отдельного \(t\) нам нужно будет делать пересчет, что влечет за собой большую вычислительную мощность. Сделаем особое преобразование для \(l(t)\): \[ l^{*}(t) = X(t_0,t)^{T}l_0 \] Подаставим это выражение в \eqref{3}: \[ \mathcal{S_t}(\tau)(P(\tau))^{1/2}X(t, \tau)^{T}X(T_0, t)^{1/2}l_{0} = \lambda_{t}(\tau)\mathcal{S_0}X_{0}^{1/2}(X(t, t_0))^{T}(X(t_0,t)^{T})l_0 \]

\[ \mathcal{S_t}(\tau)P^{1/2}(\tau)X^{T}(t_0, \tau)l_0 = \lambda_{t}(\tau)\mathcal{S_0}X^{1/2}_{0}l_{0} \]

Таким образом, в последнем равенстве у нас пропала зависимость от \(t\), следовательно, мы можем перейти к задаче Коши, которая допускает численное решение: \begin{equation} \begin{cases} \dot q\_(t) = Aq(t) \\ q\_(t_0) = q_{0} \end{cases} \end{equation}

\begin{equation} \begin{cases} \dot X \_ (t) = (\dot Q^{*}(t))^{T}(Q^{*}(t)) +(Q^{*}(t))^{T}(\dot Q^{*}(t)) \\ X\_(t_0) = X_0 \end{cases} \end{equation}

Матрица \(Q^{*}\) находится из следующей системы уравнений: \begin{equation} \begin{cases} \dot Q^{*}(t) = \mathcal{S}(t)P^{1/2}(t) + Q^{*}A^{T} \\ Q^{*}(t_0) = \mathcal{S_0}X^{1/2}_0 \end{cases} \end{equation}