Условие стабилизируемости и коррекция движения линейных систем с импульсными воздействиями и запаздыванием Б.И. Ананьев1 Н.И. Желонкина2 abi@imm.uran.ru 312115@mail.ru 1–Институт математики и механики им. Н.Н. Красовского УрО РАН (Екатеринбург) 2–Уральский федеральный университет (Екатеринбург) Аннотация Исследуются свойства стабилизируемости решений линейной систе- мы дифференциальных уравнений с запаздыванием и обобщённым воздействием в матрице системы, а также предлагаются импульс- ные методы коррекции систем без запаздывания. Получены доста- точные условия, которые обеспечивают устойчивость и асимптоти- ческую устойчивость решений системы уравнений. Для задач кор- рекции систем без запаздывания с целью минимизации конечного функционала указаны алгоритмы, согласно которым определяются моменты коррекции и возникает невозрастающая последователь- ность прогнозируемых значений функционала. Результаты приме- няются в задачах навигации и биофизики. Ключевые слова: устойчивость; стабилизируемость; системы с запаздыванием; коррекция движения; импульсные управления. 1 Введение Системы с импульсным воздействием, в настоящее время, являются современным направлением при изуче- нии теории дифференциальных уравнений, имеющей применение в различных областях знаний. Они явля- ются приложением к задачам математического моделирования в технике, электрофизике, биологии [1,2,16], экономике [4]. Кроме того, такие процессы могут сопровождаться запаздыванием, обусловленным различ- ными причинами. Например, такой процесс можно наблюдать в системах с автоматическим регулировани- ем, когда необходимо учитывать время реакции системы на управляющие процессы, время прохождения сигнала. Для таких систем важным вопросом становится проблема стабилизируемости решения. Данная работа исследует устойчивость и асимптотическую устойчивость для такого рода систем, т.е. систем с импульсным воздействием в матрице системы и запаздыванием. Будем рассматривать ситуацию, когда од- нородная система без импульсов является неустойчивой, а устойчивость и асимптотическая устойчивость достигается за счёт импульсного воздействия. Подход в данной работе основан на [1], но там рассматрива- ется система без запаздывания, а мы усложняем исходную систему, вводя запаздывание. Заметим, что в ра- ботах [5,6], импульсное воздействие играло роль возмущения,а асимптотическая устойчивость достигалась c by the paper’s authors. Copying permitted for private and academic purposes. Copyright ° In: G.A. Timofeeva, A.V. Martynenko (eds.): Proceedings of 3rd Russian Conference "Mathematical Modeling and Information Technologies" (MMIT 2016), Yekaterinburg, Russia, 16-Nov-2016, published at http://ceur-ws.org 1 за счёт асимптотической устойчивости линейной однородной системы без импульсов. Главной особенно- стью рассматриваемой системы является наличие некорректной операции умножения разрывной функции на δ-функцию,которая получается за счёт обобщённого дифференцирования ступенчатой функции, а это приводит к возникновению импульса и разрыву траектории. Под решением, так же как и в [1, 5, 7], пони- мается поточечный предел последовательности гладких решений, порождённых гладкой аппроксимацией обобщённого воздействия, если этот предел не зависит от выбора аппроксимирующей последовательности. Рассматриваются также импульсные управляемые системы, описывающие в линейном приближении от- клонение фазового вектора от номинальной траектории. Такие системы часто возникают в задачах управ- ления летательными объектами, [8]. Общие методы решения задач управления с неполной информацией предложены, например, в [9–11]. Тем не менее, специальные ограничения на возмущения и подходящая формулировка задачи дают возможность получить более продвинутые результаты. Используются методы управления в условиях неопределённости, а также информационные и совместимые множества, [11–13]. Задачи управления наблюдениями рассматривались в [14]. В данной работе расширяются и дополняют- ся результаты [15] в случае импульсных управлений. Другие задачи коррекции движения исследовались в [17, 18]. 2 Задача стабилизации Рассмотрим линейную систему дифференциальных уравнений с запаздыванием, следующего вида: ³ m X ´ ẋ(t) = A + Dj (t)v̇j (t) x(t) + Aτ x(t − τ ), (1) j=1 с начальным условием x(t) = ϕ(t). где A, Aτ – постоянные n × n матрицы, Dj (t) (j ∈ 1, m) – непрерывные ограниченные взаимно коммутатив- ные матрицы-функции размерности n × n, vi – компоненты вектор функции v(t) = (v1 (t), v2 (t), ..., vm (t))T , ϕ(t) – начальная функция, которая является функцией ограниченной вариации, определённой на [t0 −τ, t0 ]. Будем предполагать, что v(t) = (v1 (t), v2 (t), ..., vm (t))T – кусочно-постоянные функции, имеющие разрывы в точках t0 < t1 < t2 < ... < tn < ... и эта последовательность не имеет точек сгущения. Функции vi (t) в точках разрыва будем считать непрерывными слева. При сделанных выше предположениях систему (1) можно записать в следующем виде: m X ẋ(t) = Ax(t) + Aτ x(t − τ ) + Dj (t) · x(t) · 4vj (ti ) · δ(t − ti ), (2) j=1 где 4vj (ti ) = vj (ti + 0) − vj (ti ) – скачки кусочно-постоянных функций vj (t), δ(t − ti ) – δ-функция Дирака [1], сосредоточенная в момент времени ti . Решение уравнения (2) , согласно [1, 2, 7], будут удовлетворять следующему интегральному уравнению Z t Z t X x(t) = ϕ(t0 ) + Ax(ξ)dξ + Aτ x(ξ − τ )dξ + S(ti , x(ti ), 4v(ti + 0)), t0 t0 ti 0, c ≥ 1), где Y (t, −s) – фундаментальная матрица системы (5). Далее переходим к вычислению скачка S(t1 , x(t1 ), 4v(t1 )) с помощью дифференциального уравнения (4), где значение z(1) даёт начальное усло- вие для движения системы (1) на следующем промежутке (t1 , t2 ], которое строится с помощью решения уравнения (5) с начальным условием x(t1 ) + S(t1 , x(t1 ), 4v(t1 )). Поставим в соответствие разрывной траектории x(t) непрерывную траекторию x∗ (t), которая будет определяться решением дифференциального уравнения ẋ∗ (s) = Ax∗ (s) + Aτ x∗ (s − τ ), на промежутке (ti +i, ti+1 +i], а на промежутке (ti +i−1, ti +i] как решение дифференциального уравнения (4). Пусть справедливо следующее неравенство: e i (s)k ≤ beλi s (b ≥ 1, λi < 0), kD (6) e i (s)− нормированная фундаментальная матрица системы (4). здесь D В результате получили систему с переменной структурой (4ωi−s ) [19] , решением которой является непрерывная вектор-функция x∗ (s). В результате, исходной задаче поставлена в соответствие система с переключениями, где в зависимости от интервала активируется соответствующая подсистема. Наглядно это демонстрирует следующий рисунок. Рис. 1: Слева изображена траектория с разрывами, справа – непрерывная траектория x∗ (s) 2.2 Условия стабилизируемости Предположим, что x(t)− аппроксимируемое решение уравнения (1), порождаемое начальной функцией ϕ(t), которая задаётся на [t0 − τ, t0 ]. Тогда согласно формуле Коши из [20], получим: Z 0 x(t) = Y (t, −t0 ) · ϕ(t0 ) + Y (t, −t0 − τ − s)Aτ ϕ(s)ds. −τ До точки разрыва аппроксимируемое решение будет определяться следующим уравнением: Z 0 X x(t) = Y (t, −t0 ) · ϕ(t0 ) + Y (t, −t0 − τ − s)Aτ · ϕ(s)ds + Y (t, ti ) · S(ti , x(ti ), 4v(ti )). −τ ti 0. (12) 0 Rb Здесь и далее | · | – евклидова норма и a |dU (t)| = var[a,b] U – полная вариация функции U . Цель управления – минимизировать терминальный функционал |Dx(T )|, где D ∈ Rd×n – заданная матрица. Выбор неопределённых параметров {x0 , v(·), w(·)} может мешать минимизации. 3.1 Информационные и совместимые множества Изложим кратко теорию оценивания в нашем случае. В случае невырожденной матрицы P0 информаци- онное множество X(t, y, u) ⊂ Rn , состоящее из всех векторов x = x(t), которые могут реализоваться в системе (9), (10), (11), описывается уравнениями n o X(t, y, u) = x ∈ Rn : |x − x̂(t)|2P (t) + h(t) ≤ 1 , x̂(0) = 0, h(0) = 0, ˙ x̂(t) = A(t)x̂(t) + B(t)u(t) + P −1 (t)G0 (t)R(t)(y(t) − G(t)x̂(t)), ḣ(t) = |y(t) − G(t)x̂(t)|2R(t) , (13) Ṗ (t) = G0 (t)R(t)G(t) − P (t)C(t)Q−1 (t)C 0 (t)P (t) − A0 (t)P (t) − P (t)A(t), P (0) = P0 . 4 При вырожденной матрице P0 используются уравнения [15, формулы (8),(9)]. Далее рассматриваем только невырожденный случай. Введём функции f (t) = y(t) − G(t)x̂(t), t ∈ [0, T ], V (t, x) = |x − x̂(t)|2P (t) + h(t). Использование обновляющей функции f (t) и сигнала y(t) (за вычетом слагаемого с управлением) абсо- лютно эквивалентно согласно (13). Функция V (t, x) имеет смысл минимума критерия типа (11) на [0, t] при краевом условии x = x(t) и заданном сигнале и управлении. Используются далее и совместимые множе- ства V(t, y, u) ⊂ Rn × Lq2 [t, T ] × Lm 2 [t, T ], состоящие из троек {(x(t), vt (·), wt (·))}, для которых существуют функции (v(·), w(·)), удовлетворяющие (11) и такие, что выход (10) на [0, t] с условием x = x(t) совпадает почти всюду с заданным сигналом y t (·). Здесь и далее vt (·) – сужение функции v(·) на [t, T ], а y t (·) – сужение функции y(·) на [0, t]. Совместимое множество описывается формулой ( Z T³ ) ´ V(t, y, u) = (x, vt , wt ) : V (t, x) + |vt (s)|2Q(s) + |wt (s)|2R(s) ds ≤ 1 . (14) t Соотношение между совместимым и информационным множеством даётся равенством X(t, y, u) = projRn V(t, y, u). Подробнее см. [15, Лемма 2, ф-ла (13)]. 3.2 Коррекция с предписанными моментами изменения управления Пусть λ : 0 < t1 < · · · < tN +1 = T – разбиение на [0, T ]. Моменты ti называются моментами Rt коррекции управления. Введём ресурс µ(t) = µ − 0 |dU (s)| управления и рассмотрим тройку p(t) = (x̂(t), h(t), µ(t)) как позицию системы в момент t. Заметим, что совместимое множество V(t, y, u) за- висит только от пары (x̂(t), h(t)). Переход между двумя смежными позициями pi = (x̂(ti ), h(ti ), µ(ti )) и pi+1 = (x̂(ti+1 ), h(ti+1 ), µ(ti+1 )) зависит от управления ui (t) = dUi (t)/dt и обновляющей функции fi (t) на [ti , ti+1 ]. Функция ограниченной вариации может иметь скачок в момент ti . В этом случае мы различаем © Rt ª скачок слева на [ti −, ti ] и справа на [ti , ti +]. Введём множества Ui (ν) = U (·) : tii+1 |dU (t)| ≤ ν , где ν ≤ µ, i ∈ 1 : N . Пусть символ XT (ut |V(t, y, u)) означает множество достижимости системы (9) из совместимого множества V (t, y, u), где wt = 0. Это множество состоит из векторов x(T ) при переборе в (14) совокупности V(t, y, u) c wt = 0. Задача 1. Найти функции Ui∗ (t) ограниченной вариации (u∗i (t) = dUi∗ (t)/dt) на интервалах [ti , ti+1 ] и неотрицательные числа νi , i ∈ 1 : N , дающие величину J ∗ = min min max . . . min max max |Dx|, (15) ν1:N U1 ∈U1 (ν1 ) f1 (·) UN ∈UN (νN ) fN (·) x∈XT (uN |V(tN ,y,u)) PN R ti+1 где i=1 νi ≤ µ, ti |fi (s)|2R(s) ds ≤ 1 − hi , для всех i ∈ 1 : N . Полагаем U (t) ≡ 0 на [0, t1 ]. Замечание 1. Так как (9), (10) линейны, мы имеем X(t, y, u) = z(t) + X(t, ỹ, 0), где ỹ(t) = y(t) − G(t)z(t) и ż(t) = A(t)z(t) + B(t)u(t), z(0) = 0. (16) Аналогично, имеем V(t, y, u) = (z(t), 0, 0) + V(t, ỹ, 0). С этого момента пишем множества с ỹ(·) и u(·) = 0 как X(t, ỹ) и V(t, ỹ), соответственно. Поэтому XT (ut | V(t, y, u)) = z(T ) + XT (0 | V(t, ỹ)) и величина (15) может быть записана как J ∗ = min inf max . . . min max max |D(z(T ) + x)| ν1:N U1 ∈U1 (ν1 ) f1 (·) UN ∈UN (νN ) fN (·) x∈XT (0|V(tN ,ỹ)) с теми же ограничениями, как и в (15). Задача 2. В любой предписанный момент ti , i ∈ 1 : N следует найти управление uTi ∗ (·), которое даёт решение задачи: max max |Dx| → min = ji (p), (17) fi (·) x∈XT (ui |V(ti ,y,u)) Ui RT RT при ограничениях ti |fi (s)|2R(s) ds ≤ 1 − h(ti ), ti |dUi (s)| ≤ µ(ti ), и далее совершить одношаговый прогноз Ji (p, u∗i ) = max ji+1 (p), (18) fi (·) 5 Rt где tii+1 |fi (s)|2R(s) ds ≤ 1−h(ti ) и u∗i (s) = dUiT ∗ (s)/ds на интервале [ti , ti+1 ]. Если Ji (p, u∗i ) < ji (p), сохраняем управление uTi ∗ на [ti , ti+1 ]. Иначе переходим к управлению uii+1∗ , минимизирующему величину (18) при Rt ограничении tii+1 |dU (t)| ≤ µ(ti ). Управление может быть не единственным. Если так, выбираем любой минимизатор. В примерах имеет смысл принять специальное решение. В этой задаче µ(t1 ) = µ, как и выше. 3.3 Коррекция с синтезом моментов коррекции В любой момент t с позицией p(t) определяем величину jt (p) по формуле (17), где ti заменяется на t. После этого делаем прогноз на интервал [t, τ ] подобно (18): Jt (p, τ, u∗ ) = max jτ (p), (19) ft (·) Rτ где t |ft (s)|2R(s) ds ≤ 1 − h(t) и u∗ (s) = dUtT ∗ (s)/ds на интервале [t, τ ]. Введём величину J∗ (t, p, u∗ ) = inf Jt (p, τ, u∗ ), (20) τ ∈[t,T ] которая определяет минимум худшего прогноза (19). Процедура начинается с t1 > 0 с µ(t1 ) = µ. Тогда если J∗ (t, p, u∗ ) < jt (p) для всех t близких к t1 и таких что t > t1 , мы продолжаем использовать u∗t1 до первого момента t = t2 > t1 , когда J∗ (t, p, u∗ ) = jt (p). В момент t2 программное минимаксное управление модифицируется, и мы переходим к оптимальному управлению uTt2∗ (·), дающему решение задачи подобно (17). Процедура повторяется на [t2 , T ]. Поэтому приходим к последовательности моментов коррекции {tk }, t1 < t2 < . . . ≤ T. Как увидим ниже, все минимумы и максимумы в формулах (15) – (20) достигаются. 4 Минимаксные решения задач коррекции Рассмотрим задачи в подсекции 3.2. Для краткости обозначим x̂(ti ) = x̂i , h(ti ) = hi , и µ(ti ) = µi , p(ti ) = pi . Введём функцию будущих потерь Wi (p, νi:N ) = min max . . . min max max |Dx|, Ui ∈Ui (νi ) fi (·) UN ∈UN (νN ) fN (·) x∈XT (uN |V(tN ,y,u)) PN R tj+1 где j=i νj ≤ µi , tj |fj (s)|2R(s) ds ≤ 1 − hj , для всех j ∈ i : N . Через Wi (p) обозначим величину minνi:N Wi (p, νi:N ). Нетрудно заметить, что J ∗ = W1 (p) в (15) и функции Wi (p, νi:N ), Wi (p), удовлетво- ряют рекуррентным соотношениям Wi (p, νi:N ) = min max Wi+1 (p, νi+1:N ), Wi (p) = min max Wi+1 (p), (21) Ui ∈Ui (νi ) fi (·) Ui (·) fi (·) R ti+1 R ti+1 где ti |fi (s)|2R(s) ds ≤ 1 − hi , ti |dUi (t)| ≤ µi , i ∈ 1 : N . Соотношения (21) имеют граничное условие n o WN +1 (p) = max 2 |Dx| = max l0 Dx̂(T ) + (1 − h(T ))1/2 |D0 l|P −1 (T ) . x:|x−x̂(T )|P (T ) ≤1−h(T ) |l|≤1 На последнем шаге соотношений (21), когда i = N , используя граничное условие, получаем ½ ³ ³ ´´1/2 ¾ 2 0 2 WN (p, νN ) = max r(l; tN )x̂N − max |r(l; s)B(s)|νN + (1 − hN ) λ(tN )(1 − |l| ) + |D l|P (T,tN ) , |l|≤1 s∈[tN ,T ] WN (p) = WN (p, µN ), 6 где используем обозначения r(l; s) = l0 DX(T, s), ∂P (t, s)/∂t = A(t)P (t, s) + P (t, s)A0 (t) + C(t)Q−1 (t)C 0 (t), P (s, s) = P −1 (s), λ(s) = max |D0 l|2P (T,s) . |l|≤1 Формула для WN (p, νN ) получена с использованием элементарного неравенства n o max k 1/2 A + (1 − hN − k)1/2 B = (1 − hN )1/2 (A2 + B 2 )1/2 , k∈[0,1−hN ] где A ≥ 0, B ≥ 0, и максимум достигается здесь в r∗ = (1−hN )A2 (A2 +B 2 )−1/2 . Оптимизация по f (·) выпол- RT нена при ограничении tN |f (s)|2R(s) ds = k. Если λmax – максимальное собственное число для симметриче- ской и положительно определённой матрицы Q, мы используем тот факт, что conc |l|Q на шаре совпадает с ¡ ¢1/2 λmax (1 − |l|2 ) + |l|2Q , [10, 11]. Здесь символ conc ϕ(l) означает минимальную вогнутую функцию, мажо- рирующую ϕ(l) на единичном шаре. Наконец, используется теорема о минимаксе, [22]. Здесь оптимальное ∗ управление UN (·) необходимо удовлетворяет соотношению Z T r(l∗ ; s)B(s)dUN ∗ (s) = − max |r(l∗ ; s)B(s)|νN , tN s∈[tN ,T ] ∗ где вектор l максимизирует в формуле для WN (p, νN ). Продолжая вычисления на следующих шагах, приходим к условиям оптимальности для задачи 1. На шаге i получаем ½ ¾ Wi (p, νi:N ) = max r(l; ti )x̂i − max |r(l; s)B(s)|νi + ϕi (l) , где |l|≤1 s∈[ti ,ti+1 ] ½ ½ Z ti+1 ¾¾ −1 0 (22) ϕi (l) = conc − max |r(l; s)B(s)|νi+1 + max r(l; s)P (s)G (s)R(s)fi (s)ds + ϕi+1 (l) , s∈[ti+1 ,ti+2 ] fi (·) ti i ∈ 1 : N − 1. R ti+1 Здесь ti |fi (s)|2R(s) ds ≤ 1 − hi . Оптимальное управление необходимо удовлетворяет соотношениям Z ti+1 r(l∗ ; s)B(s)dUi∗ (s) = − max |r(l∗ ; s)B(s)|νi , (23) ti s∈[ti ,ti+1 ] где l максимизирует в задаче (22). Величина функционала в задаче 1 равна J ∗ = minν1:N W1 (p, ν1:N ) = ∗ PN W1 (p), где i=1 νi = µ. Чтобы решить задачу 2, нужно вычислить величины (17), (18). Мы получаем ½ ³ ³ ´´1/2 ¾ 2 0 2 ji (p) = max r(l; ti )x̂i − max |r(l; s)B(s)|µi + (1 − hi ) λ(ti )(1 − |l| ) + |D l|P (T,ti ) , (24) |l|≤1 s∈[ti ,T ] ½ Z ti+1 Ji (p, u∗i ) = max ji+1 (p) = max r(l; ti )x̂i + r(l; s)B(s)dUi∗ (s) − max |r(l; s)B(s)|µi+1 fi (·) |l|≤1 ti s∈[ti+1 ,T ] ³ ³ ´´1/2 ¾ + (1 − hi ) λ(ti+1 )(1 − |l|2 ) + |D0 l|2P (T,ti ) . Процедура управления в задаче 2 начинается с i = 1 и приводит к последовательности позиций pi , где j1 (p) ≥ j2 (p) ≥ · · · ≥ jN (p). Действительно, сравним величины ji (p) и ji+1 (p). Если Ji (p, uTi ∗ ) < ji (p), мы получаем ji (p) > ji+1 (p). Иначе используем управление ui+1∗ i , которое минимизирует величину Ji (p, ui ) в (18). Поэтому ½ Z ti+1 ½ min Ji (p, ui ) = min max r(l; ti )x̂i + r(l; s)B(s)dUi (s) + conc − max |r(l; s)B(s)| ui (·) Ui ∈U (µi ) |l|≤1 ti s∈[ti+1 ,T ] µ Z ti+1 ¶ ³ ³ ´´1/2 ¾¾ 2 0 2 × µi − |dUi (s)| + (1 − hi ) λ(ti+1 )(1 − |l| ) + |D l|P (T,ti ) ≤ ji (p). ti 7 Имеем λ(ti+1 ) ≤ λ(ti ). Это следует из соотношения ∂P (T, s)/∂s = −X(T, s)P −1 (s)G0 (s)R(s)G(s)P −1 (s)X 0 (T, s), откуда в свою очередь следует, что норма симметрической и положительно определённой матрицы P (T, s) убывает по s. Если заменим λ(ti+1 ) на λ(ti ), левая часть неравенства возрастает, и в conc{·} мы имеем вогнутую функцию по l. Тогда min и max можно переставить по теореме о минимаксе. Итак, минимизация по Ui (·) даёт ji (p). Замечание 2. Процедура вычисления оптимальных управлений в задаче 1 более сложна, чем в задаче 2. Но мы можем упростить её, если немного увеличим функцию будущих потерь. Докажем неравенство Wi (p) ≤ ji (p), i ∈ 1 : N . Из (22), (24), имеем WN (p) = jN (p). По индукции предположим. что Wi+1 (p) ≤ ji+1 (p). Тогда из (21): ½ µ Wi (p) = min max Wi+1 (p) ≤ min max ji+1 (p) = min max r(l; ti )x̂i − max |r(l; s)B(s)| µi Ui (·) fi (·) Ui (·) fi (·) Ui (·) |l|≤1 s∈[ti+1 ,T ] Z ti+1 ¶ Z ti+1 ³³ ´´1/2 ¾ − |dUi (s)| + r(l; s)B(s)dUi (s) + (1 − hi ) λ(ti+1 )(1 − |l|2 ) + |D0 l|2P (T,ti ) . ti ti Если заменим здесь λ(ti+1 ) на λ(ti ), правая сторона неравенства возрастёт. После этого переставим max и min и проведём минимизацию по Ui (·). В результате получим ji (p). Отметим, что Ui (s) ≡ 0, если maxs∈[ti+1 ,T ] |r(l; s)B(s)| ≥ maxs∈[ti ,ti+1 ] |r(l; s)B(s)|. Иначе имеем Z ti+1 Z ti+1 r(l∗ ; s)B(s)dUi (s) = − max |r(l∗ ; s)B(s)|µi , |dUi (s)| = µi , (25) ti s∈[ti ,ti+1 ] ti где l∗ максимизирует в формуле (24). В этом случае ресурс управления исчерпывается в момент ti+1 . Замечание 3. Из результатов [23] следует, что все элементы проекции на Rd множества достижимости системы (16) могут быть достигнуты линейной комбинацией не более чем (d + 1) импульсными воздей- ствиями (δ-функциями). Поэтому, оптимальное управление Ui , удовлетворяющее соотношениям вида (23) или (25), имеет не более чем (d + 1) скачков. Если точка лежит на границе множества достижимости, соответствующее управление имеет только d скачков. 4.1 Пример Рассмотрим уравнение ẋ(t) = B(t)(u(t)+v(t)), t ∈ [0, T ], где диагональная матрица B(t) = [b1 (t), 0; 0, b2 (t)] ∈ RT R2 . Наблюдаем сигнал y(t) = x(t) + w(t), где 0 (2|v(t)|2 + |w(t)|2 )dt ≤ 1. 8 7 6 jt functional 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 10 t sec Рис. 2: Изменение функционала будущих потерь 8 √ √ Вектор-функция U (·) стеснена ограничениями (12) с µ = 1. Пусть b1 (s) √ = s, b2 (s) = T − s, D = [1, −1], T = 10. На отрезке √ [1/2, T ] реализуются √ возмущения v(t) = [1; 1] 2 sin(t)/7, w(t) = 0, а сигнал на [0, 1/2] равен y(t) = [1/ 8 + 1; 1/ 8 − 1]. Пусть на отрезке нанесена достаточно частая сетка, где возможна корректировка√ управления. Оптимальное управление√в задаче равно u∗i (s) = −Dx̂i B(s)D0 δ(s − T )/T , если |Dx̂i | ≤ T , и u∗i (s) = −sign(Dx̂i )B(s)D0 δ(s − T )/ T в противном случае. Таким образом, управление состоит из одного импульса, прилагаемого в конце процесса. Корректирование производится в конце процесса в зависимости от оценки x̂N , вычисляемой согласно измеренному сигналу. Изменение функционала показано на рисунке 2. Благодарности Исследование выполнено за счёт гранта Российского научного фонда (проект № 16-11-10146). Список литературы [1] S. T. Zavalishchin, A. N. Sesekin. Dynamic Impulse Systems: Theory and Applications. Kluwer Academic Publishers, Dordrecht, 1997. [2] B. M. Miller, E. Ya. Rubinovich. Discontinuous solutions in the optimal control problems and their representation by singular space-time transformations. Automation and Remote Control, 74: 1969– 2006, 2013. [3] V.A. Dychta, O.N. Samsonyuk. Optimalnoye impulsnoye upravleniye s prilozheniyami [Optimal Impulse Control with Applications]. Moscow, Fizmatlit, 2003. (in Russian) = В. А. Дыхта, О. Н. Самсонюк. Оптимальное импульсное управление с приложениями . Москва. Физматлит, 2003. [4] V.P. Maksimov. Positional Parry of Impulse Disturbances in Control Problem for Linear Retarded System. Vestnik Permskogo universiteta, Seriya: Ekonomika, 2(21): 6–14, 2014. (in Russian) = В. П. Максимов. Позиционное парирование импульсных возмущений в задаче управления ли- нейной системой с последействием. Вестник Пермского университета, Серия: Экономика, 2(21): 6–14, 2014. [5] I.A. Kornilov, A.N. Sesekin. On Stability of Linear Systems with Matrix Containing Generalized Functions. Vestnik UGTU-UPI, 3(33):386–388, 2004. (in Russian) = И. А. Корнилов, А. Н. Сесе- кин. Об устойчивости линейных систем с матрицей, содержащей обобщенные функции. Вестник УГТУ-УПИ, 3(33):386–388, 2004. [6] N.I. Zhelonkina, A.N. Sesekin, S.P. Sorokin. On Stability of Linear Retarded Systems with Impulse Disturbance in the Matrix of the System. Vestnik Udmurtskogo universiteta. Matematika. Mechanika. Komputernuye nauki, 1:40–46, 2011. (in Russian) = Н. И. Желонкина, А. Н. Сесекин, С. П. Сорокин. Об устойчивости линейных систем с импульсным воздействием в матрице системы и запаздыва- нием. Вестник Удмуртского университета. Математика. Механика. Компьютерные науки, 1:40–46, 2011. [7] A.N. Sesekin. Dynamical Systems with Nonlinear Impulse Structure. Sbornik nauchnych trudov. IMM UrO RAN, 6(2):497–514, 2000. (in Russian) = А. Н. Сесекин. Динамические системы с нелинейной импульсной структурой Сборник научных трудов. Тр. ИММ УрО РАН, 6(2):497–514, 2000. [8] R. Pratt (ed). Flight Control Systems. The American Institute of Aeronautics and Astronautics, Reston, VA 20191-4344, USA, 2000. [9] N. Krasovskiĭ, A. Subbotin. Game-Theoretical Control Problems. Springer-Verlag, New York, 1988. [10] A. Kurzhanski, I. Valyi. Ellipsoidal Calculus for Estimation and Control. SCFA. Birkhäuser, Boston, 1997. [11] A. Kurzhanski, P. Varaiya. Dynamics and Control of Trajectory Tubes: Theory and Computation. SCFA. Birkhäuser, Boston, 2014. 9 [12] D. Bertsecas, I. Rhodes. Recursive state estimation for a set-membership description of uncertainty. IEEE Trans. on Auto. Control, AC-16:117–128, 1971. [13] F. Schweppe. Uncertain Dynamic Systems. Prentice Hall, New York, 1973. [14] B. Ananyev. Some problems of observations’ control. AIP Conference Proceedings, 1404:235–241, 2011. [15] B. Ananyev, N. Gredasova. Positional and impulse strategies for linear problems of motion correction. AIP Conference Proceedings, 1789:34–39, 2016. [16] V.A. Dychta, O.N. Samsonyuk. Optimalnoye impulsnoye upravleniye s prilozheniyami [Optimal Impulse Control with Applications]. Moscow, Fizmatlit, 2003. (in Russian) = В. А. Дыхта, О. Н. Самсонюк. Оптимальное импульсное управление с приложениями . Москва. Физматлит, 2003. [17] N. Parusnikov, V. Morozov, V. Borzov. Zadachi korrektcii v inertcialnoi navigatcii [Problems of Correction in Inertial Navigation]. MGU, Moscow, 1982. (in Russian) = Н. Парусников, В. Морозов, В. Борзов. Задачи коррекции в инерциальной навигации. МГУ, Москва, 1982. [18] F. Chernousko, I. Ananievski, S. Reshmin. Control of Nonlinear Dynamical Systems. Methods and Applications. Springer-Verlag, Berlin, 2008. [19] D. Liberzon, A. Morse. Basic problems in stability and design of switched systems Control Syst. Mag. 19:59–70, 1999. [20] A. N. Sesekin, N. I. Zhelonkina. On the stability of linear systems with generalized action and delay. IFAC-PapersOnLine, Proceedings of the 18th IFAC World Congress Milano, 13404–13407, 2011. [21] R. Liptser, A. Shiryaev. Statistics of Random Processes, V.1 General Theory, V.2 Applications. 2-nd ed. Springer-Verlag, Berlin, 2001. [22] K. Fan. Fixed-point and minimax theorems in locally convex topological linear spaces. Proc. Nat. Acad. Sci., 38(2):121–126, 1952. [23] L. Newstadt. Optimization, a moment problem, and nonlinear programming. J. SIAM Control, 2(1):33–53, 1964. 10 A condition of stabilizability and motion correction for linear systems with impulse actions and delay Boris I. Ananyev Krasovskii Institute of Mathematics and Mechanics (Yekaterinburg, Russia) Natal’ya I. Zhelonkina Ural Federal University (Yekaterinburg, Russia) Abstract. Stability and asymptotic stability properties are investigated for linear systems with delay and with impulse action in the system matrix. Impulse methods of motion correction for linear systems without delay are suggested as well. A theorem is proved that provide a sufficient condition of stabilizability for such systems. For control systems without delay, a motion correction problem is considered with the aim to minimize the final functional. An algorithm is given that defines instants of correction. According to the algorithm, the nonincreasing sequence of predictable values of functional is arisen. The results are used in problems of navigation and biophysics. Keywords: stability; stabilizability; systems with dekay; motion correction; impulse controls. 11