Научная статья на тему 'Алгоритм наблюдателя пониженного порядка системы управления курсом судна для оценки возмущений и шумов измерений'

Алгоритм наблюдателя пониженного порядка системы управления курсом судна для оценки возмущений и шумов измерений Текст научной статьи по специальности «Математика»

CC BY
634
55
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ДИНАМИЧЕСКАЯ МОДЕЛЬ СУДНА / НАБЛЮДАТЕЛЬ ПОНИЖЕННОЙ РАЗМЕРНОСТИ / ВЕКТОР ВЫХОДА / АЛГОРИТМ / МОДАЛЬНЫЙ МЕТОД СИНТЕЗА СИСТЕМ / DYNAMIC MODEL OF A SHIP / THE OBSERVER REDUCED DIMENSION VECTOR OUTPUT / ALGORITHM / MODAL SYNTHESIS METHOD OF SYSTEMS

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

Для повышения надежности и обеспечения требуемых качественных показателей системы автоматического управления курсом судна необходима достоверная информация о векторе состояния динамической модели, адекватной реальному объекту. При воздействии на судно внешних возмущений, приводящих к появлению помех в измеряемых сигналах, требуется обеспечить понижение уровня шумов путем фильтрации и применения математических датчиков информации. Датчики могут быть синтезированы с использованием динамических наблюдателей и оценивателей. Поскольку в условиях эксплуатации авторулевых вектор состояния не всегда доступен для измерения, необходимо управлять процессом по вектору выхода. Для оптимального управления требуется восстановление информации по вектору выхода. При движении судна в условиях волнения внешние возмущения приводят к возникновению дополнительных составляющих погрешности в оценке вектора состояния. При этом значительно ухудшаются качественные показатели судовых авторулевых комплексов, и для их повышения требуется использовать методы и вычислительные процедуры стохастической фильтрации.Рассматривается метод и алгоритм динамического наблюдателя пониженной размерности (наблюдателя Давида Луенбергера), отличающийся тем, что он позволяет уменьшить погрешность оценки вектора переменных состояния судового объекта, функционирующего в условиях изменяющейся внешней среды, путем построения модели состояния расширенной системы «объект среда». Оценка производится на основе алгоритмов линеаризации нелинейных моделей в пространстве состояния с использованием оптимальной фильтрации, инструментария матричной лаборатории и компьютерных технологий. Приведен пример расчета трехмерного динамического наблюдателя для оценки векторов возмущений и шумов измерений, выполненный на основе предложенного алгоритма. Результаты оценки переменных состояния и возмущений, полученные с использованием канонических форм уравнений состояния, соответствуют оптимальным (по быстродействию) переходным процессам в рассматриваемых динамических моделях при переводе из заданного начального состояния в установившееся состояние за минимальное время.

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

Похожие темы научных работ по математике , автор научной работы — Сахаров Владимир Васильевич, Чертков Александр Александрович, Сабуров Сергей Валерьевич

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

REDUCED-ORDER OBSERVER ALGORITHM COURSE MANAGEMENT SYSTEM OF THE SHIP TO ASSESS DISTURBANCE AND NOISE MEASUREMENTS

To increase reliability and ensure the required quality indicators of the automatic control system course ship need reliable information on the dynamic model of the State vector, adequate to the real object. When exposed to ship outside disturbances, resulting in noise in the measured signals are required to ensure that lowering noise by filtering and sensor application of mathematical information. Sensors can be synthesized using dynamic observers and evaluators. Since the conditions of use of autopilots vector state is always available for the dimension, you must manage the process on the vector output. For optimal control requires data recovery on vector output. When moving a ship in the midst of the excitement of external perturbations cause additional components of error in assessing the vector state. This significantly deteriorate the qualitative indices of ship autopilots complexes and their increasing need to use methods and computational procedures stochastic filtering.The method and algorithm of dynamic observers reduced dimensionality (Observer David Luenberger), notable for the fact that it reduces the error estimation of the vector of state variables of the channel object, functioning in the face of changing external Wednesday, by constructing a model of the extended status object system Wednesday. Score is based on algorithms of linearization of nonlinear state-space models using optimal filtration, laboratory matrix instrumentation and computer technology. An example of the calculation of the three-dimensional dynamic observer to evaluate the perturbation vectors and noise measurements made on the basis of the proposed algorithm. The results of the evaluation of state variables and perturbations produced using canonical forms of State equations correspond to the best (performance) the transitional processes in these dynamic models when translating from a given initial state in the settled condition in minimal time.

Текст научной работы на тему «Алгоритм наблюдателя пониженного порядка системы управления курсом судна для оценки возмущений и шумов измерений»

ВЕСТНИК*]

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ШЦ^

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА^,

ИНФОРМАЦИОННЫЕ ТЕХНОЛОГИИ И АВТОМАТИЗАЦИЯ НА ТРАНСПОРТЕ

DOI: 10.21821/2309-5180-2017-9-1-211-220

REDUCED-ORDER OBSERVER ALGORITHM COURSE MANAGEMENT SYSTEM OF THE SHIP TO ASSESS DISTURBANCE AND NOISE MEASUREMENTS

V. V. Saharov, A. A. Chertkov, S. V. Saburov

Admiral Makarov State University of Maritime and Inland Shipping, St. Petersburg, Russian Federation.

To increase reliability and ensure the required quality indicators of the automatic control system course ship need reliable information on the dynamic model of the State vector, adequate to the real object. When exposed to ship outside disturbances, resulting in noise in the measured signals are required to ensure that lowering noise by filtering and sensor application of mathematical information. Sensors can be synthesized using dynamic observers and evaluators. Since the conditions of use of autopilots vector state is always available for the dimension, you must manage the process on the vector output. For optimal control requires data recovery on vector output. When moving a ship in the midst of the excitement of external perturbations cause additional components of error in assessing the vector state. This significantly deteriorate the qualitative indices of ship autopilots complexes and their increasing need to use methods and computational procedures stochastic filtering.

The method and algorithm of dynamic observers reduced dimensionality (Observer David Luenberger), notable for the fact that it reduces the error estimation of the vector of state variables of the channel object, functioning in the face of changing external Wednesday, by constructing a model of the extended status object system Wednesday. Score is based on algorithms of linearization of nonlinear state-space models using optimal filtration, laboratory matrix instrumentation and computer technology. An example of the calculation of the three-dimensional dynamic observer to evaluate the perturbation vectors and noise measurements made on the basis of the proposed algorithm. The results of the evaluation of state variables and perturbations produced using canonical forms of State equations correspond to the best (performance) the transitional processes in these dynamic models when translating from a given initial state in the settled condition in minimal time.

Keywords: dynamic model of a ship, the observer reduced dimension vector output, algorithm, modal synthesis method of systems.

For citation:

Saharov, Vladimir V., Alexandr A. Chertkov, and Sergei V. Saburov. "Reduced-order observer algorithm course management system of the ship to assess disturbance and noise measurements." Vestnik Gosudarstvennogo universiteta morskogo i rechnogo flota imeni admirala S.O. Makarova 9.1 (2017): 211— 220. DOI: 10.21821/2309-5180-2017-9-1-211-220.

УДК. 621.398.04

2 О

7

АЛГОРИТМ НАБЛЮДАТЕЛЯ ПОНИЖЕННОГО ПОРЯДКА СИСТЕМЫ УПРАВЛЕНИЯ КУРСОМ СУДНА

ДЛЯ ОЦЕНКИ ВОЗМУЩЕНИЙ И ШУМОВ ИЗМЕРЕНИЙ О

ы у

В. В. Сахаров, А. А. Чертков, С. В. Сабуров К

ФГБОУ ВО «ГУМРФ имени адмирала С. О. Макарова», ^Г

Санкт-Петербург, Российская Федерация.

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

ЛВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

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

Рассматривается метод и алгоритм динамического наблюдателя пониженной размерности (наблюдателя Давида Луенбергера), отличающийся тем, что он позволяет уменьшить погрешность оценки вектора переменных состояния судового объекта, функционирующего в условиях изменяющейся внешней среды, путем построения модели состояния расширенной системы «объект - среда». Оценка производится на основе алгоритмов линеаризации нелинейных моделей в пространстве состояния с использованием оптимальной фильтрации, инструментария матричной лаборатории и компьютерных технологий. Приведен пример расчета трехмерного динамического наблюдателя для оценки векторов возмущений и шумов измерений, выполненный на основе предложенного алгоритма. Результаты оценки переменных состояния и возмущений, полученные с использованием канонических форм уравнений состояния, соответствуют оптимальным (по быстродействию) переходным процессам в рассматриваемых динамических моделях при переводе из заданного начального состояния в установившееся состояние за минимальное время.

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

Для цитирования:

Сахаров В. В. Алгоритм наблюдателя пониженного порядка системы управления курсом судна для оценки возмущений и шумов измерений / В. В. Сахаров, А. А. Чертков, С. В. Сабуров // Вестник Государственного университета морского и речного флота имени адмирала С. О. Макарова. — 2017. — Т. 9. — № 1. — С. 11-220. DOI: 10.21821/2309-5180-2017-9-1-211-220.

Введение

Системы автоматизации и управления технологическими процессами на современных водоизмещающих судах характеризуются высокой размерностью и сложностью. Для обеспечения оптимального управления судовыми динамическими системами требуется получение достоверной информации о переменных состояния. Обеспечение заданной точности и надежности управления объектом при воздействии возмущений на измерения углов курса и крена, а также угловой скорости вращения судна требует компромиссных решений, связанных с выбором быстродействия системы и точности оценки переменных состояния. Вектор состояния не всегда доступен для измерений. Если измерение вектора состояния затруднено по каким-либо причинам, можно воспользоваться математическими датчиками, выполненными на основе наблюдателей неполной размерности. Измерения возмущений связаны с появлением дополнительных составляющих ошибки оценки переменных состояния. Наряду с инструментальными датчиками информации, применение математических датчиков позволяет идентифицировать параметры и структуру управляемых систем средствами и методами теории идентификации с оценками ^ на основе аналоговых и цифровых фильтров. В системах реального времени оценки должны производиться с использованием методов и вычислительных процедур стохастической и, в част-1 ности, биномиальной фильтрации.

« Развитие и совершенствование способов конструирования динамических наблюдателей,

г> предложенных в [1] - [7], а также их практическое применение на судах позволяют значительно

еЗ_ повысить эффективность и надежность эксплуатации объектов водного транспорта средствами

212] автоматизации, создать механизм алгоритмизации процессов энергосбережения. С учетом специфики воздействия на объект внешних сил и моментов в допустимой области вариации параметров нелинейные модели, как показано в [5], [6] и [10] - [12], приводятся к линейным моделям. Вместе с тем практическое использование наблюдателей для оценки векторов возмущений и шумов измерений на конкретных объектах также имеет свою специфику, которая должна учитываться в процессе синтеза [13] - [14].

Кроме наблюдателей полного порядка, или наблюдателей Калмана, у которых размерность вектора состояния равна размерности вектора состояния объекта, могут быть предложены также наблюдатели Давида Луенбергера пониженного порядка [4]. Такие наблюдатели используют информацию о состоянии объекта, непосредственно содержащуюся в измерениях вектора выхода. Это позволяет построить алгоритм оценки с помощью модели порядка (n - p), где p = rank (С). Обычно p = 1.

(i)

Основная часть

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

х (7) = Ах(7) + Ви (7); у ) = Сх(7),

где А — (п х п)-матрица; В — (п х £)-матрица и С — (р х п)-матрица постоянных коэффициентов; х() — (п х 1)-мерный вектор состояния; и© и у() — векторы управления и выхода соответствующей размерности. Измерения производятся для части переменных состояния, если р < п. Для системы управления с одним входом матрица В для системы с одним входом имеет размерность (п х 1).

Пусть ранг матрицы С равен р. Для упрощения уравнения выхода выполним преобразование базиса в уравнении (1). Выберем произвольную (п -р) х п-матрицу F так, чтобы матрица

\F

T =

C

была невырожденной. Это всегда возможно выполнить, так как rank С = p. Введем новый вектор состояния z (t) = T • х (t) и представим его в данном виде:

z (t) =

'w(t) 7 (t)

}n - p; } P,

где w (^) е Яп-р, у(^) е Яр, т. е. выходы объекта совпадают в выбранном базисе с последними р компонентами его вектора состояния. Выполнив преобразование базиса с матрицей Т, перейдем к уравнениям состояния:

W (t) У (t)

A

A12

A22

w(t) У (t)

+

Bi

u (t).

(2)

Из системы (2) можно выделить структуру (подсистему) порядка (n - p) с известными (доступными измерению) входами u(t), y(t). Заметим, для нее всегда можно выбрать требуемые коэффициенты характеристического уравнения.

Рассмотрим алгоритм синтеза наблюдателя пониженного порядка на примере наблюдателя (n - 1) порядка для системы со скалярным выходом (rank С = 1). Алгоритм синтеза наблюдателя пониженной размерности состоит из следующих шагов [4], [9].

1. Путем линейного преобразования уравнения состояния объекта (матрицы A, B, С) приведем преобразования к форме идентификационного канонического представления (ИКП), или наблюдаемой форме Давида Луенбергера [4], при которой матрица A является транспонированной матрицей Фробениуса:

Af =

0 0

0 1

0 0

(3) СИ

где ар а2 ап — коэффициенты характеристического многочлена объекта; матрица С равна

С = [0, ..., 0 11. (4)

rs. О

2. Выберем желаемые коэффициенты р. характеристического полинома наблюдателя:

(det /„-, - Ан) = s n-1 + piS n-2 + ... +

3. Строим матрицу преобразования P вида:

-р.-i

p=

К-1 -Р„-2

0

0

-Pi 1

P-1 =

I.-1

0

р.-1

в. - 2

р1 0 1

(5)

Нетрудно заметить, что если матрица приведена к виду ИКП, то в результате преобразования получим

A = PAP=

0 0 1 0

0 0 0 0

а также

B =

0 -p„_i К -Pi)PB_i - ап

0 -ри_2 (ai -Pi)Pn-2 - ап-i +РИ_1

1 -Pi (ai -Pi)Pr -a2 +P2 0 i -ai +pi

-вп -A

A -вп-2bn

Ьп - 1-P 1A b

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

(6)

(7)

2Í4J

Выделим из матрицы Ai верхнюю треугольную подматрицу А порядка (n - i) и расположим ее в верхнем левом углу матрицы А. Образуем первые (n - 1)-строк матрицы B, из которых сформируем матрицу (вектор-столбец) B. Введем обозначение верхних (n - 1)-строк последнего столбца матрицы А символом an. Тогда уравнение наблюдателя Луенбергера примет вид

x (t) = Ax (t) + any(t) + Bu(t).

Для оценки вектора состояния в исходном базисе, представленном в форме ИКП, образуем вектор x (t) = col {x (t ),y (t )} .Оценку состояния в базисе ИКП получим по формуле

x(t) = P-1 x (t).

Рассмотрим процедуру синтеза наблюдателя (n - 1)-порядка для оценки возмущений и шумов измерений. Сначала выберем модель управляемого объекта с возмущениями, которые интерпретируются решениями дифференциальных уравнений с изменяющимися начальными условиями, что позволяет косвенно оценить влияние внешних воздействий. С учетом полученных оценок будем иметь модель внешних воздействий, которую введем в состав расширенной динамической системы (объект, внешние возмущения) и получим наблюдатель. В результате оценки будут содержать как оценки вектора состояния объекта, так и внешних воздействий.

В тех случаях, когда внешние воздействия моделируются рядами, состоящими из многочленов вида

te'-P, (t),

i=1

процедура синтеза наблюдателя значительно упрощается. Здесь X. е Л — постоянные; P. (t) — многочлены с заданными коэффициентами. Упрощение обеспечивается за счет того, что такие модели представляют собой систему линейных дифференциальных уравнений с постоянными

коэффициентами, для которых, например, в системе MATLAB содержатся функции, предназначенные для исследования инвариантных во времени систем с помощью LTI (Linear Time Invariant) объектов. LTI-объекты можно исследовать во временной и частотной областях с помощью специальных функций практически в режимах прямых вычислений.

Рассмотрим алгоритм оценки внешних воздействий для модели объекта, заданной уравнениями состояния:

(8)

I х (I) = А(1) х(1) + В (I)ы (I) + f (I); 1 у (I) = С (I) х(г)+V (I),

при начальных условиях х(^0) = х0, ^ > Здесь х(1)еЯп — вектор состояния объекта; ы($ еЯт,у(0еЯ — векторы управления и выхода; А() Б(), С(1) — матрицы линейной нестационарной системы; /(¿), у(1) — воздействующие на объект, соответственно, внешние возмущения и шум (погрешность) измерений. В процессе синтеза наблюдателя выполняются измерения и(), у(/), неизмеряемыми переменными являются х(/), /(¿),

Представим внешние воздействия /(¿), у(1) в виде выходных процессов линейной системы, заданной в пространстве состояний уравнениями:

(9)

хс(*) = Ас(*) хс(*);

у Л*) = Сс (*) хс(хДО = хсо; * ^

где х() — вектор состояния среды размерности (п х 1); Ус(0 — вектор выхода модели внешних возмущений размерности ((пс + I) х 1).

Заметим, что у (Л = соНДЛ, А , С — известные матрицы. Составная матрица С

C =

'Cf C

содержит подматрицу С/ размерности (п х пс) и подматрицу Су размерности (I х пс), которые определяют связь хс(^) и /() с помехами в (1). В системе (9) вектор начальных состояний хс0 также неизвестен.

Введем расширенный («совокупный») вектор состояния системы «объект - среда»:

хо) = (х(),хс()) е япхп<. Объединив уравнения (8) и (9), получим уравнения расширенной системы в виде:

х (г) = А(г )х (г) + Ви (г); у (г) = С (г) х (г); х (г0) = х0; г > г0.

Матрицы состояния, управления и выхода этой системы имеют вид:

(10)

A =

A C,

0

B =

B

0

C = [C Cl.

Видно, что для динамической системы (10) размерности - = п + пс может быть предложен наблюдатель со следующей структурой:

x(t) = A(t)x(t) + B(t)u(t) + L(t) (y(0 - y(0); y(t) = C(t)x(t), x(L) = xn, t>L.

(11)

где х е^" — вектор состояния наблюдателя, служащий оценкой состояния модели «объект -среда»; уеЯ1 — вектор выхода наблюдателя, служащий оценкой выхода расширенной системы объекта и среды; Щ) — матрица (размерности (п х I)) коэффициентов обратной связи по невязке между выходами расширенной системы и наблюдателя.

Рассмотрим процедуру оценки влияния возмущений и шумов измерений на примере упрощенной модели судна:

2 О

7

Г215 I

«ВЕСТНИК

ЩШ ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

= и (') + М (),

(12)

где 1х — момент инерции судна относительно продольной оси; 8(?) — угол крена судна; и()) М(?) — управляющий и возмущающий моменты.

сС Э

Если вектор и(?) известен и угловая скорость крена юх ) = — (7) доступна для измерений,

с

то оценке подлежит возмущающий момент М(?).

Для переменной М(?), моделируемой линейной функцией М(?) = М0 + V? с неизвестными М0 и V, предлагается модель

| М () = V ();

Ы) = 0

с неизвестными начальными условиями М (0) и V (0), на основе которой в работе получено уравнение состояния вида (10) с матрицами А, В, С:

(13)

А =

0 1/ 1Х 0

0 0

в =

1/ 1х 0 0

; С = [1 0 0].

(14)

Построим модель наблюдателя (п - 1)-порядка для оценки возмущений и шумов измерений. Исходной моделью объекта будут служить уравнения расширенной системы (12), (13) с матрицами (14). На основе изложенного алгоритма уравнения состояния системы сначала приводятся к виду ИКП. С помощью формул (3) и (4) построим матрицы А, С канонического вида:

"0 0 0" А =10 0 0 1 0 _

Затем определяются матрицы наблюдаемости Q исходной системы и ( — системы вида ИКП. В результате находятся матрицы преобразования Т:

"10 0 1 Г0 0 1

; С = [0 0 1].

Q =

0 I

-1

0

0 0 I

; Q =

0 1 0 1 0 0

г> о

Т = Q^О =

0 0 I

-1

0 I;1

1 о

0 0

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

А У 5) = 5 2 +Т2О05 + 0о2,

где с помощью параметра задается быстродействие наблюдателя. Тогда Р: Матрицу преобразования Р приведем к каноническому виду (5):

Р =

"1 0 -О02

о 1 -л/2ос 0 0 1

Р-1 =

1

0 1 0 0

0 а2

л/2а

1

Теперь нетрудно с помощью (6) - (7) сформировать матрицы:

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

A = PAP=

B =

0

-П02 л/2П0

Матрица А содержит подматрицу А второго порядка и верхние две строки матрицы В, из которых образуем вектор-столбец ВВ. Тогда имеем:

0 -По2

1 -V2Q0 0 1

-"о С J-1

A =

0 -fi02 i -Яа

-V2n03 -n 2

в =

-n 0i-

—i '0

i

Приходим к уравнениям наблюдателя второго порядка:

Г (t ) = -n 02 X, (t ) - V2Q 03юх (t ) - Q021-lu (t ); 1X2 (t ) = 3cl (t ) - V2n0 x2 (t ) - n 02 юх (t ) - V2n0 rxlu (t ).

(15)

T _ t-ip-i _

0 0

0 4

1 0

I J2Q0 Ix "02

моделируется соотношениями:

|M (t ) = I^it ) + I^V2Q0 i2(t );

I Я0 = (t) + lA2^«-

(16)

Видно, что с помощью уравнений (15), (16) конструируется наблюдатель Луенбергера для рассматриваемой задачи.

Вычисления для рассматриваемого примера выполнены c помощью программы, составленной в кодах MATLAB с использованием инструментария Symbolic Toolbox.

Результаты моделирования системы (16) в среде MATLAB представлены на рис. 1 и 2, где M — возмущающий момент, изменяющийся во времени; M— момент, восстановленный с помощью динамического наблюдателя; v и v — соответственно исходное h восстановленное значения шумов измерений.

Mit)

2 о

7

Г217 I

150 200 250 300

Время t, с

Рис. 1. Переходный процесс оценивания возмущений в обобщенной системе

b

1

ЛВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

О 2 X

100 150 200

Время t, с

Рис. 2. Переходный процесс оценивания шумов измерений

300

Е)

Из графиков видно, что время переходного процесса в наблюдателе составляет не более 100 с. Параметры моделируемой системы: 1х = 1000 кг-м2; М(0) = 0,25 Н-м; у(0) = 5-10-3 Нм/с; = 0,05 1/с.

Выводы

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

Для линейной модели наблюдателя устойчивость гарантируется обоснованным выбором характеристического многочлена наблюдателя — многочлена Баттерворта второго порядка, корни которого мало отличаются от корней характеристического уравнения замкнутой системы управления судном. Выбор обобщенной модели наблюдателя путем введения векторов хД t), х2( t), х3( t) позволяет предложить способ мониторинга датчиков угловой скорости и угла дрейфа, являющихся переменными состояния объекта. Наблюдатели фактически являются математическими датчиками информации и могут использоваться как для восстановления недоступных измерению переменных состояния, так и переменных, измерение которых в реальных условиях затруднено.

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

СПИСОК ЛИТЕРАТУРЫ

1. Сазонов А. Е. Модальный метод синтеза наблюдателя для системы управления курсом судна / А. Е. Сазонов, В. В. Сахаров, А. А. Чертков // Вестник Государственного университета морского и речного флота имени адмирала С. О. Макарова. — 2016. — № 4 (38). — С. 211-223. DOI: 10.21821/2309-5180-2016-84-211-223.

ВЕСТНИК«)

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ^^

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

2. Сахаров В. В. Синтез оптимального оценивателя для системы управления судовым динамическим объектом / В. В. Сахаров, О. В. Шергина, А. А. Чертков // Вестник Государственного университета морского и речного флота имени адмирала С. О. Макарова. — 2013. — № 1 (20). — С. 26-31.

3. Ramdani N. Computing reachable sets for uncertain nonlinear monotone systems. / N. Ramdani, N. Meslem, Y. Candau // Nonlinear Analysis: Hybrid Systems. — 2010. — Vol. 4. — Is. 2. — Pp. 263-278. DOI: 10.1016/ j.nahs.2009.10.002.

4. Luenberger D. G. Introduction to Dynamic Systems, Theory, Models, and Applications / D. G. Luenberger. — New York: John Wiley & Sons, 1979. — 446 p. DOI: 10.1002/bs.3830260412.

5. Андриевский Б. Р. Избранные главы теории автоматического управления с примерами на языке MatLab / Б. Р. Андриевский, А. Л. Фрадков. — СПб.: Наука, 2000. — 475 с.

6. Yang J. A complete solution to a simple case of dynamic observer error linearization: New approach to observer error linearization / J. Yang, J. Back, J. H. Seo // IEICE transactions on fundamentals of electronics, communications and computer sciences. — 2011. —Vol. 94. — No. 1. — Pp. 424-429.

7. Митряшкин Ю. В. Линейные модели управляемых динамических систем. Ч. 1. Уравнения «вход -выход» и «вход - состояние - выход»: учеб. пособие / Ю. В. Митряшкин. — М.: Изд-во МГТУ им. Н. Э. Баумана, 2008. — 222 c.

8. Hui S. Observer design for systems with unknown inputs / S. Hui, S. H. Zak // International Journal of Applied Mathematics and Computer Science. — 2005. — Vol. 15. — Is. 4. — Pp. 431-446.

9. Chen X. State observer for a class of nonlinear systems and its application to machine vision / X. Chen and H. Kano // IEEE Transactions on Automatic Control. — 2004. — Vol. 49. — Is. 11. — Pp. 2085-2091. DOI: 10.1109/TAC.2004.837529.

10. Boutat D. On the transformation of nonlinear dynamical systems into the Extended Nonlinear Observable Canonical Form / D. Boutat, K. Busawon // International Journal of Control. — 2011. — Vol. 84. — Is. 1. — Pp. 94-106. DOI: 10.1080/00207179.2010.541285.

11. Zheng G. A single output dependent observability normal form / G. Zheng, D. Boutat, and J. P. Barbot // SIAM Journal on Control and Optimization. — 2007. — Vol. 46. — Is. 6. — Pp. 2242-2255. DOI:10.1137/050627137.

12. Tomera M. Nonlinear controller design of a ship autopilot / M. Tomera // International Journal of Applied Mathematics and Computer Science. — 2010. — Vol. 20. — Is. 2. — Pp. 271-280. DOI: 10.2478/v10006-010-0020-8.

13. Wang Y. Observer design using a generalized time-scaled block triangular observer form / Y. Wang, A. F. Lynch // Systems & Control Letters. — 2009. — Vol. 58. — Is. 5. — Pp. 346-352. DOI: 10.1016/j. sysconle.2008.12.005.

14. McGookin E. W. Ship steering control system optimisation using genetic algorithms / E. W. McGookin, D. J. Murray-Smith, Y. Li, T. I. Fossen // Control Engineering Practice. — 2000. — Vol. 8. — Is. 4. — Pp. 429-443. DOI: 10.1016/S0967-0661(99)00159-8.

REFERENCES

1. Sazonov, Anatoliy Efimovich, Vladimir Vasilevich Saharov, and Alexandr Alexandrovich Chertkov. "A modal method for the synthesis of an observer for heading control." Vestnik Gosudarstvennogo universiteta morskogo i rechnogo flota imeni admirala S. O. Makarova 4(38) (2016): 211-223. DOI: 10.21821/2309-5180-2016- ^ 8-4-211-223. i

N,

2. Saharov, V. V, O. V. Shergina, and A. A. Chertkov. "Optimal estimator synthesis for ship dynamic object o control." Vestnik gosudarstvennogo universiteta morskogo i rechnogo flota imeni admirala S. O. Makarova 1(20) T (2013): 26-31. i

3. Ramdani, Nacim, Nacim Meslem, and Yves Candau. "Computing reachable sets for uncertain nonlinear monotone systems." Nonlinear Analysis: Hybrid Systems 4.2 (2010): 263-278. DOI: 10.1016/j.nahs.2009.10.002. 4

4. Luenberger D. G. Introduction to Dynamic Systems, Theory, Models, and Applications. New York: John Wiley & Sons, 1979. DOI: 10.1002/bs.3830260412.

5. Andrievskij, B. R., and A. L. Fradkov. Izbrannye glavy teorii avtomaticheskogo upravlenija sprimerami na jazyke MatLab. SPb.: Nauka, 2000.

6. Yang, J., J. Back, and J. H. Seo. "A complete solution to a simple case of dynamic observer error linearization: New approach to observer error linearization." IEICE transactions on fundamentals of electronics, communications and computer sciences 94.1 (2011): 424-429.

ЛВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

7. Mitrjashkin, Ju. V. Linejnye modeli upravljaemyh dinamicheskih sistem. Chast 1. Uravnenija «vhod — vyhod» h «vhod — sostojanie — vyhod»: uchebnoeposobie. M.: izd-vo MGTU im. N. Je. Baumana, 2008.

8. Hui, Stefen, and Stanislaw H. Zak. "Observer design for systems with unknown inputs." International Journal of Applied Mathematics and Computer Science 15.4 (2005): 431-446.

9. Chen, Xinkai, and Hiroyuki Kano. "State observer for a class of nonlinear systems and its application to machine vision." IEEE Transactions on Automatic Control 49.11 (2004): 2085-2091. DOI: 10.1109/TAC.2004.837529.

10. Boutat, Driss, and Krishna Busawon. "On the transformation of nonlinear dynamical systems into the extended nonlinear observable canonical form." International Journal of Control 84.1 (2011): 94-106. DOI: 10.1080/00207179.2010.541285.

11. Zheng, Gang, Driss Boutat, and Jean-Pierre Barbot. "Single output-dependent observability normal form." SIAM Journal on Control and Optimization 46.6 (2007): 2242-2255. DOI:10.1137/050627137.

12. Tomera, M. "Nonlinear controller design of a ship autopilot." International Journal of Applied Mathematics and Computer Science 20.2 (2010): 271-280. DOI: 10.2478/v10006-010-0020-8.

13. Wang, Yebin, and Alan F. Lynch. "Observer design using a generalized time-scaled block triangular observer form." Systems & Control Letters 58.5 (2009): 346-352. DOI: 10.1016/j.sysconle.2008.12.005.

14. McGookin, Euan W., David J. Murray-Smith, Yun Li, and Thor I. Fossen. "Ship steering control system optimisation using genetic algorithms." Control Engineering Practice 8.4 (2000): 429-443. DOI: 10.1016/S0967-0661(99)00159-8.

ИНФОРМАЦИЯ ОБ АВТОРАХ

INFORMATION ABOUT THE AUTHORS

Сахаров Владимир Васильевич —

доктор технических наук, профессор ФГБОУ ВО «ГУМРФ имени адмирала С.О. Макарова»

198035, Российская Федерация, Санкт-Петербург,

ул. Двинская 5/7

e-mail: _saharov_@rambler.ru

Чертков Александр Александрович —

кандидат технических наук, доцент

ФГБОУ ВО «ГУМРФ имени

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

адмирала С.О. Макарова»

198035, Российская Федерация, Санкт-Петербург,

ул. Двинская 5/7

e-mail: chertkov51@mail.ru

Сабуров Сергей Валерьевич — аспирант

Научный руководитель:

Сахаров Владимир Васильевич

ФГБОУ ВО «ГУМРФ имени

адмирала С. О. Макарова»

198035, Российская Федерация, Санкт-Петербург,

ул. Двинская 5/7

e-mail: kaf_osnipr@gumrf.ru

Saharov, Vladimir V. —

Dr. of Technical Sciences, professor

Admiral Makarov State University

of Maritime and Inland Shipping

5/7 Dvinskaya Str., St. Petersburg 198035,

Russian Federation

e-mail: _saharov_@rambler.ru

Chertkov, Alexandr A. —

PhD, associate professor

Admiral Makarov State University

of Maritime and Inland Shipping

5/7 Dvinskaya Str., St. Petersburg 198035,

Russian Federation

e-mail: chertkov51@mail.ru

Saburov, Sergei V. — Postgraduate

Supervisor.

Saharov, Vladimir V.

Admiral Makarov State University

of Maritime and Inland Shipping

5/7 Dvinskaya Str., St. Petersburg 198035,

Russian Federation

e-mail: kaf_osnipr@gumrf.ru

Статья поступила в редакцию 9 ноября 2016 г.

Received: November 9, 2016.

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