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

Вопросы синтеза цифровых законов управления морскими судами с учетом неопределенностей Текст научной статьи по специальности «Математика»

CC BY
48
29
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МОРСКОЕ СУДНО / ИНВАРИАНТНЫЙ ЭЛЛИПСОИД / ЦИФРОВОЕ УПРАВЛЕНИЕ / НЕОПРЕДЕЛЕННОСТИ

Аннотация научной статьи по математике, автор научной работы — Смирнов Михаил Николаевич, Смирнова Мария Александровна

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Вопросы синтеза цифровых законов управления морскими судами с учетом неопределенностей»

Смирнов М.Н.1, Смирнова М.А.2

1 Санкт-Петербургский государственный университет, г.Санкт-Петербург, ассистент кафедры

компьютерных технологий и систем, 51глгпоу-тп@таП . ги

2 Санкт-Петербургский государственный университет, г.Санкт-Петербург, ассистент кафедры

компьютерных технологий и систем, 51тлгпоуа-та@Ьк.ги

ВОПРОСЫ СИНТЕЗА ЦИФРОВЫХ ЗАКОНОВ УПРАВЛЕНИЯ МОРСКИМИ СУДАМИ С УЧЕТОМ НЕОПРЕДЕЛЕННОСТЕЙ

КЛЮЧЕВЫЕ СЛОВА

Морское судно, инвариантный эллипсоид, цифровое управление, неопределенности. АННОТАЦИЯ

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

Уровень развития современных компьютерных технологий, непрерывный рост вычислительных мощностей, появление новых инструментальных программных средств - все это способствует повсеместной автоматизации функционирования механических объектов различного рода с использованием автономных бортовых систем с цифровой реализацией.

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

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

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

Задача о подавлении внешних возмущений с известными характеристиками относится к основным проблемам теории управления и рассматривается в различных ее разделах и приложениях. В качестве примера можно привести такие работы, как [1 - 13]. Прикладные задачи управления морскими судами рассматриваются в работах [14-29]. В указанных работах известные методы синтеза адаптируются под специфические особенности, присущие динамике морских судов.

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

которой посвящена статья.

Рассмотрим цифровую линейную модель динамики судна

X [ к + 1 ]=Ax [ к ] + BS [ к ]+Dw [ к ], S [ к + 1 ]=u [ к ] + S [ к ], (1)

^ у [ к ] = Cx [к ],

где X £ E - вектор состояния (здесь он определяет отклонения от положения равновесия), Ö£ Em - вектор состояния исполнительных органов, w £ E - вектор внешних возмущающих

воздействий, u £E - вектор управляющих сигналов (управлений), A, B, D , C и M - матрицы соответствующих размерностей с постоянными компонентами.

Для простоты будем считать, что измеряемый выход у системы одновременно является и контролируемой переменной.

Будем формировать обратную связь для системы (1) в виде

1х[к]\

u[k] = Kxx[k]+KsS[k] = K ...... ,К=(КТ ! Kg), (2)

где матрицы KX и K& соответствующих размерностей имеют постоянные компоненты.

Обозначим через йк множество стабилизирующих регуляторов, для которых корни характеристического полинома

A3(z^) = det(En+mz-A0-B0K),^J^S\Bc^^ (3)

замкнутой системы (1), (2) лежат внутри единичного круга.

Рассмотрим также сужение множества стабилизирующих регуляторов, определяя его желаемыми модальными требованиями к замкнутой системе:

й5к = { K£йк: ^(K )£Cä, i = M+"m) , (4)

где (K) - корни характеристического полинома (3) замкнутой системы, CH - заданная область на комплексной плоскости. В частности, в качестве такой области можно принять круг C&={ z £ С1: |z |<ad) , где <Х d D (0, l) - заданное вещественное число, определяющее степень устойчивости замкнутой системы. Как и для аналоговых регуляторов [6, 11, 18], возможны и другие способы введения допустимой области [30].

Будем считать, что внешнее возмущающее воздействие ограничено по величине, т.е. справедливо соотношение

||w[к]||да<1, к=0,1,2,... . (5)

Определим функционал Jd= Jd (K) , характеризующий размер минимального инвариантного эллипсоида, включающего множество Iea реакций на внешние воздействия для замкнутой системы (1), (2):

Jd= Jd (K )= f (P (а о ))=tr (С з P (ao) C3) , (6)

где а0= arg min f (P(а)) , f (P(а )) = tr (С3P (а )C3) , P(а) - положительно определенное

0<а<1 1 '

решение линейного матричного уравнения

а 1— а

A3 P А3 — а P + -—- D3 D3=0 , (7)

A=A0 + B0K ,C3 = (C ! 0) , Dj =

Определение. Эллипсоид с центром в начале координат

£х= { хеЕп: хР-1х<11 , Р>0, будем называть инвариантным по переменной х с (по состоянию) для динамической системы

х = А х + D w,

с с с с '

е = С х .

с с

^ (t )||<1 при 0 < Г <с, если из условия х (0 следует х ^ )е £х для всех моментов времени 0 < ( <с .

Симметрическую матрицу Р при этом будем называть матрицей эллипсоида £ х .

В определении и далее по тексту факт положительной определенности матрицы отражается обозначением Р> 0 .

Иначе говоря, любая траектория системы х ^) , начинающаяся внутри эллипсоида £х , в каждый последующий момент времени будет оставаться внутри него, что качественно показано на рис. 1.

Рис. 1. Инвариантный эллипсоид

Рассмотрим задачу о выборе стабилизирующего регулятора (2), который минимизирует размер Jd инвариантного эллипсоида с учетом желаемых модальных свойств замкнутой системы

Jd=Jd(Kmin , (8)

K ^ k

где допустимое множество Qsk определяется формулой (4).

Если не принимать во внимание дополнительные модальные свойства, то мы приходим к постановке известной задачи

Jd= Jd(K И min , (9)

K

которая рассмотрена в работе [31 - 33], где приводится следующее ее решение:

Jd= Jd(K)= min Jd(K)=F(P(a,ß))=tr(C3P(a,ß)C3) , (10)

K £Qk

где (a,ß)= arg min F(P(a,ß)) , F(P(a ,ß)) = tr(C3P(a,ß)C3) ,P(a ,ß) -

a>0, ß>0

положительно определенное решение линейного матричного уравнения

A0 PAo-2 ß B0 Bo+ß2 B0 Bo( A0)"1 P-1 A-1 B0 Bo-a P + ^ D3 D=0 (11) K = argmin Jd (K )=-ß B0 (A0 )-1 P"1 (a,ß) . (12)

K £ö.

Введем обозначение для шуровского характеристического полинома замкнутой системы (1), (2), (12)

Лз(г,а,р (Еп+тх- А0-В0К(а,в)). (13)

Заметим, что решение задачи (8) удобно связать с параметризацией вещественными

_ г-,п+т

векторами у характеристических полиномов системы, замкнутой регуляторами из

допустимого множества Qsk , введенного соотношением (4).

Указанная параметризация множества характеристических полиномов определяет соответствующую параметризацию множества ^& регуляторов (2), которыми замыкается объект (1). Действительно, зададим произвольный вектор у е Еп6т и по указанным формулам построим полином А* (г, у). Для того, чтобы он был характеристическим полиномом для замкнутой системы необходимо и достаточно, чтобы матрица К коэффициентов регулятора (2) удовлетворяла тождеству

Л з( 2, К ^ (Еп+т* - Ао-Во К ) = л (г,у) .

_ | т-тХп+тХт

Если собрать все компоненты матрицы К в вектор к £ Е , то этому тождеству

будет соответствовать эквивалентная линейная система

Гк = 1л(у) ,А0

An -

■а ; вл ч«Го,

в

о

О

Е

где компоненты матрицы Г определяются только матрицами А0 и В0 , а вектор m(y) определяется коэффициентами полинома A*(s, у) и характеристического полинома матрицы А0 .

Заметим, что указанная система всегда совместна в силу условия полной управляемости. Она содержит П 6 m уравнений и имеет m X П 6 m X m неизвестных, т.е.

n

nc = mXn+mXm — n— m компонент вектора k , собранные в вектор hcGE c , могут быть

!1 a

Y,hcIGE , A = n + m+nc , и,

дополнительно задавая произвольный вектор hc , найдем соответствующее решение k (£ ) = k ( Y, hc) линейной системы. Таким образом, найдена матрица K = K (£ ) = K (Y, hc) коэффициентов регулятора (2) из множества Jsk , параметризованного векторами £ = { Y, hc| . На базе приведенной параметризации можно сформулировать следующее утверждение: Теорема. Для функционала Jd (K) в задаче (8) существует минимизирующая

последовательность {K(£i)} регуляторов (2), определяемая последовательностью векторов {вг-} такой, что

lim {K (£. )} = Ко = arg min Jä (K) , lim I Jd (K (£ ))|=Jd (K0 ) = J° .

i ^ да К GQsk i ^ да 1 >

Доказательство данного утверждения проводится аналогично доказательству теоремы для непрерывного времени, однако необходимо заметить, что для построения полинома с заданной степенью устойчивости используется следующая схема [29]:

Для любого вектора у d End степень устойчивости вспомогательного полинома, построенного по формулам

, I если nd - четное;

A',r.Vi= (14)

{(z - arf+1 (-у, а.Л)Л если nd - нечетное

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

не меньше наперед заданной величины а d D (0, l), и обратно, если степень устойчивости некоторого полинома А(s) не меньше величины аd D (0, l), то можно указать такой вектор у d ЕП , что А(z) к А*(z, у) , причем

Д(z,Y)=П(z2+a1(Y^d)Z +a0(уа)) , где (15)

i= 1

a1 (Y^d ) = —а-

exp

1 i-

2 4

Yi Yi 1 2 — 2" + — Y22

/ i-\

exp

Y21 Jy41 2 — 2- — ^ 4~" — Y2 2

. ^ ' ^ . и

0)'

(16)

a0(Yа)=а2ехР( —У^). i = 1,d> ad+1(yа) = adexP(—Yd

У = {Уц, У12 , У21, У 22 , •••, У А1, У А 2, Уа 01 •

Таким образом, приведенная теорема является основой для построения алгоритма решения задачи (8). Схему решения задачи (8) можно сформулировать в виде следующего алгоритма.

Алгоритм.

1. Взять любую точку уеЕп+т и построить вспомогательный полином Л (г,у) по формулам (15) - (16) либо по аналогичным формулам [30] для более сложного варианта модальных условий.

2. Сформировать систему линейных уравнений Гк = т( у) , где к - вектор, собранный из компонент матрицы К , обеспечивающую выполнение тождества

Лз( 2, К )=л (г,у) ,

которая всегда совместна и, если ее решение не является единственным, осуществить произвольный выбор вектора hc£ Е с свободных переменных (компонент матрицы К ) по отношению к этой системе.

3. После подстановки в систему Гк = т( у) принятого вектора £ = {у , ^ }е ЕХ, X = п +т +пс , найти ее решение к =к (£) , а, соответственно, и матрицу К = ~ ( £ ) .

4. Подставить найденную матрицу К = К (£) коэффициентов обратной связи (2) в соотношения (6), (7) и найти для замкнутой системы (1), (2) соответствующее значение J А= J а (К (£ ))= Jd ( £) размера минимального инвариантного эллипсоида.

5. С помощью любого допустимого численного метода решения задачи

Jd=Jd (^ ) = Jd (K (£)) + min

£ £EX

на безусловный экстремум, задать новую точку £ и, повторяя пункты 3,4, минимизировать функцию Jd(£) .

6. Процесс завершить после нахождения точки £o = arg min Jd (£) , определяя при этом К ££Ex

матрицу K0 = К ( £о) , которая принимается в качестве решения задачи

Jd=Jd (K)+ min

K к

обеспечивающего минимальное значение min Jd(K ) = Jd(K0)=Jd = Jd(£0) размера

K

минимального инвариантного эллипсоида.

Таким образом, в статье рассмотрены особенности синтеза цифровых законов управления морскими судами с учетом неопределенностей и предложен алгоритм решения задачи (8).

Литература

1. Веремей Е.И. Синтез законов многоцелевого управления движением морских объектов // Гироскопия и навигация. - 2009. - № 4. С. 3-14.

2. Лукомский Ю. А., Корчанов В. М. Управление морскими подвижными объектами.- СПб.: Элмор, 1996.- 320 с.

3. Лукомский Ю. А., Пешехонов В.Г., Скороходов Д.А. Навигация и управление движением судов. СПб.: «Элмор», 2002. 360с.

4. Мирошник И. В. Теория автоматического управления. Нелинейные и оптимальные системы. - СПб: Питер, 2005. -271 с.

+

5. Назин А.В., Назин С.А., Поляк Б. Т. О сходимости внешних эллипсоидальных аппроксимаций областей достижимости линейных дискретных динамических систем // АиТ. 2004. З8. С. 39-61.

6. Смирнов М.Н. Метод учета ограниченных внешних воздействий при синтезе обратных связей с многоцелевой структурой // Вестн. С.-Петерб. ун-та. Сер. 10: Прикладная математика, информатика, процессы управления. 2014. Вып. 2. С. 130-140.

7. Смирнова М.А. Обеспечение астатизма в системах управления движением морских судов // Вестн. С.-Петерб. унта. Сер. 10: Прикладная математика, информатика, процессы управления. 2014. Вып. 2. С. 141-153.

8. Смирнов М.Н. Информационная поддержка процесса обучения при моделировании системы управления шаром на наклонной направляющей // International journal of open information technologies. 2014.T.2, №3. С. 23-28.

9. Смирнова М.А. Вопросы информатизации обучения на примере программного комплекса цифрового управления роботом // International journal of open information technologies. 2014.T.2, №3. С. 29-34.

10. Смирнов М.Н. Использование современных информационных технологий для моделирования системы управления шаром на подвижной направляющей // Современные информационные технологии и ИТ-образование. М.: ИНТУИТ.РУ 2013. №9. С. 728-732.

11. Смирнов М.Н. Оптимизация управления подвижными объектами с ограниченными внешними возмущениями // Современные информационные технологии и ИТ-образование. М.: ИНТУИТ.РУ 2012. № 8. С. 1018-1024.

12. Смирнова М.А. Программный комплекс цифрового управления роботом, имитирующим башню танка // Современные информационные технологии и ИТ-образование.. М.: ИНТУИТ.РУ, 2013. №9. С. 733-737.

13. Смирнова М.А. Многоцелевое управление подвижными объектами в режиме реального времени // Современные информационные технологии и ИТ-образование. М.: ИНТУИТ.РУ 2012. №8. С. 1025-1032.

14. Смирнов М.Н. Динамическая компенсация ограниченных внешних возмущений в системе стабилизации курса судна / / Материалы XV конференции молодых ученых "Навигация и управление движением". 2013. С. 364-370.

15. Смирнова М.А. Синтез астатических законов управления с неполной обратной связью для морских подвижных объектов // Материалы XV конференции молодых ученых "Навигация и управление движением". 2013. С. 217223.

16. Поляк Б. Т., Щербаков П. С. Робастная устойчивость и управление. М.: Наука, 2002.

17. Смирнов М.Н., Федорова М.А. Компьютерное моделирование системы астатической стабилизации курса морского судна // Процессы управления и устойчивость: Труды 41-й международной научной конференции аспирантов и студентов. СПб, 2010. С. 495-500.

18. Смирнов М.Н. Алгоритм синтеза управлений, подавляющих ограниченные внешние воздействия на морское судно // Процессы управления и устойчивость: Труды 42-й международной научной конференции аспирантов и студентов. СПб, 2011. С. 362-368.

19. Федорова М.А. Синтез и компьютерное моделирование астатической системы управления курсом морского судна // Процессы управления и устойчивость: Труды 42-й международной научной конференции аспирантов и студентов. СПб, 2011. С. 368-374.

20. Смирнов М.Н., Смирнова М.А. Реализация программного комплекса для динамического управления нелинейным объектом // Процессы управления и устойчивость: Труды 44-й международной научной конференции аспирантов и студентов. СПб, 2013. C.297-301.

21. Смирнов М.Н., Смирнова М.А. Современные информационные технологии в процессе обучения технических специалистов // Процессы управления и устойчивость, 2014. T. 1. С. 397-400.

22. Арзуманян Н.К., Смирнова М.А., Смирнов М.Н. Моделирование системы управления тормозными усилиями автомобиля в Simulink // Процессы управления и устойчивость, 2015. T. 2. С. 357-363.

23. Какорин Н.С., Смирнов М.Н., Смирнова М.А. Моделирование распространения вирусного заболевания на основе данных о взаимодействии индивидов // Процессы управления и устойчивость, 2015. T. 2. С. 421-425.

24. Королев Е.А., Петрунин В.Н., Смирнов М.Н., Смирнова М.А. Разработка системы шифрования с применением стегано-криптографических методов / / Процессы управления и устойчивость, 2015. T. 2. С. 437-441.

25. Смирнов А.Н., Смирнов М.Н., Смирнова М.А. Управление аккомпанементом в процессе воспроизведения // Процессы управления и устойчивость, 2015. T. 2. С. 497-502.

26. Смирнов М.Н., Смирнова М.А. Вопросы синтеза стабилизирующих управлений при наличии неопределенных внешних возмущений // Процессы управления и устойчивость, 2015. T. 2. С. 503-508.

27. Арзуманян Н.К., Смирнова М.А., Смирнов М.Н. Синтез и моделирование обратных связей антиблокировочной системы автомобиля // Устойчивость и процессы управления: Материалы III международной конференции, посвященной 85-летию со дня рождения профессора, чл.-корр. РАН В.И. Зубова. Спб, 2015. С. 505-506.

28. Веремей Е.И., Смирнов М.Н., Смирнова М.А. Вопросы устойчивости и качества в управлении движением морских судов // Устойчивость и процессы управления: Материалы III международной конференции, посвященной 85-летию со дня рождения профессора, чл.-корр. РАН В.И. Зубова. Спб, 2015. С. 511-512.

29. Веремей Е.И. Алгоритмы решения одного класса задач HC -оптимизации систем управления // Известия РАН. Теория и системы управления. — 2011. — № 3. — С. 52 - 61.

30. Веремей Е. И. Обеспечение заданной степени устойчивости регуляторами с неполной информацией // Известия АН СССР. Техн. кибернетика. 1986. № 4. С. 123-130.

31. Назин С. А., Поляк Б. Т., Топунов М. В. Подавление ограниченных внешних возмущений с помощью метода инвариантных эллипсоидов. Автомат. и телемех., 2007, № 3, 106-125

32. Поляк Б. Т., Щербаков П.С. Трудные задачи линейной теории управления. Некоторые подходы к решению // АиТ. 2005. №5. С. 7-46.

33. Поляк Б.Т., Щербаков П. С. Техника D-разбиения при решении линейных матричных неравенств // АиТ. 2006. N11. С. 159-174.

i Надоели баннеры? Вы всегда можете отключить рекламу.