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

УДК 656.256:581.51

А. Т. Когут, Н. Ю. Панфилова

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

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

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

Изучение нелинейных динамических систем относится к одной из сложных задач современной теории автоматического управления [1], где в настоящее время интенсивно развиваются два направления - траекторное управление [2] и топологические методы [3] анализа и синтеза многомерных систем. При проектировании систем управления сложными многосвязными нелинейными динамическими объектами [4], особенно в роботизированных и мехатронных комплексах, выделяют этап построения требуемых траекторий изменения физических координат, обеспечивающих согласованное движение объекта, и этап синтеза локальных регуляторов, поддерживающих заданные состояния объектов. Траекторное управление реализуется в локальных подсистемах, где необходимо сформировать такое управляющее воздействие, при котором состояние объекта или его выход изменялись бы в соответствии с требуемой траекторией или программой. Дополнительные условия, накладываемые на устойчивость и качество процессов в замкнутых системах, относят траекторное управление к классу обратных задач динамики [5; 6].

Топологический подход [1] или методы дифференциальной геометрии [3] позволяют теоретически найти такое преобразование, при котором в новых переменных состояния многомерный нелинейный объект описывается линейными уравнениями. Главная сложность заключается в существовании и определении диффеоморфизма между исходной и линеаризованной моделями. Практические результаты получены для аффинных объектов и траекторного управления, а задача синтеза локальных регуляторов решается методами обратных задач динамики [2; 5].

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

Рекуррентные соотношения широко используются в численных методах, поэтому и предлагается подход, основанный на применении в дискретном регуляторе алгоритмов в форме итерационных процедур. Исследуемый объект в локальной системе управления состоит из нелинейного элемента (НЭ) и линейной части [2; 5; 6]. В рабо-

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

Постановка задачи. Пусть дискретный стационарный объект описывается разностным уравнением

хк+1 = Ахк + Вф(ик),

Х0 = х (0) . (1) Здесь х - я-мерный вектор переменных состояния; А и В - матрица и вектор с известными постоянными параметрами размерами я х я и я х 1 соответственно; и - управляющее воздействие; ф(и) - непрерывная функция, являющаяся характеристикой нелинейного элемента, которая принадлежит классу ф(и) е С3, т. е. существуют ее первая, вторая и третья производные.

В системах траекторного управления для наблюдаемого вектора состояния хк+ 1 должно выполняться условие

хк +1 = &к + и (2)

где g - я-мерный вектор требуемого закона изменения состояния объекта (системы) или желаемая траектория движения.

В задачах обратной динамики [6; 7] подставляют равенство (2) в модель объекта (1) и по сути получают уравнение

gk +1 -Ахк -вф(ик) = 0 которое можно свести к скалярному

В+( gk +1 - Ахк )-ф(ик ) = °. (3)

В формуле (3) под В понимается псевдообратная матрица размером 1 х п, для которой необходимо, чтобы соблюдались условия

В+-В = 1 и В • В+ • Я = Я, где Я - я-мерный вектор.

Введем в рассмотрение следующую непрерывную функцию:

/ ^к+и хк, ик ) = В+^к+1 - Ахк )-ф(ик). (4)

Элементы векторов gk + 1 и хк известны, а численное значение ик необходимо определить, поэтому /(•) является функцией только одного аргумента и, и выражение (4) можно записать в виде

/(и) = В+(gk +1 -Ахк)-ф(и) . (5)

Несложно убедиться, что для частных производных выполняются равенства

/(|) (и) = -ф(|) (и) , I = 1,2, ... (6)

и функция/ (и) также принадлежит классу /(и) е С3.

Из сравнения выражения (5) для/(и) и уравнения (3) следует, что при нахождении управления и = ик должно выполняться равенство

/ (и ) = 0. (7)

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

Синтез рекуррентных алгоритмов формирования управляющих воздействий. Одним из наиболее известных численных процедур первого порядка, использующих только первую производную, является метод Ньютона [8] или Ньютона-Рафсона [9], когда формируется дискретная последовательность ик на основе следующего правила:

ик = ии -1 -[ Л'-1 ]-1 Л-1. (8)

Здесь как и в дальнейшем

/Л = /(0 (ик-,), I = 0,1, 2,3. (9)

Рекуррентные формулы методов второго порядка содержат как первую, так и вторую производную от функции/и). В общей теории численного решения нелинейного уравнения (7), приведенной и рассмотренной в работе [8], выделяются метод Чебышева, для которого строится итерационный процесс вида

[( г. -)’

/, -, (П -, )’ + 0.5/;., /I, и алгоритм Хэлли, вычисляющий приближение ик по формуле

ик = ик - 1 -

1-1 (11)

ии ----- ии 1

(10)

(п-1) -0,5/:_1 /к-1

Л -1Л -1.

Более общие методы второго порядка, использующие линеаризацию /(ик), предлагаются в [10], и рекуррентные процедуры на их основе имеют вид

ик = ик-1-[/к'-1]-1 [/к-1 + 0,5/к"-1 82]; (12)

= ик -1-[ Я -1 + 0,5 /к"-1 8к ]-1 /к -г

(13)

Конкретные алгоритмы определения управляющих воздействий ик в рекуррентных формах получаются, если в выражения (8), (10) и (11) подставить формулу (5) для функцииДи) и учитывать соотношение (6) для ее частных производных. Алгоритмы второго порядка, построенные на основе методов Чебышева и Хэлли, целесообразно реализовать в виде двухступенчатой процедуры. Первоначально, используя метод Ньютона-Рафсона, вычисляется значение у, по формуле (8), затем - 8к, а на второй ступени управляющее воздействие ик уточняется в соответствии с выражениями (12) или (13).

Свойства итерационных процедур в прикладной математике [9] оценивают по степени приближения дискретной последовательности ик к точному решению и* нелинейного уравнения (7) и для этого получают и исследуют аналитические зависимости е, = ик - и * от Ек -; = ик -1 - и *. Например, по показателю скорости сходимости метод Ньютона-Рафсона имеет квадратичную, а алгоритмы Чебышева и Хэлли - кубическую сходимость [8; 9]. Анализ малых отклонений е, -1 ® 0 позволяет также говорить и о более высокой точности методов второго порядка в определении и *, так как для них е, пропорциональна е, -1, в то время как у Ньютона-Рафсона - ек -1. При реализации рекуррентных процедур формирования ик в замкнутых системах управления точность оценивается по степени близости выходной переменной к входному или задающему воздействию [1; 2; 5].

Анализ точности процессов управления в замкнутый системах. Динамическая ошибка в замкнутой системе вк + 1 представляет собой я-мерный вектор и вычисляется по формуле

еь +1 = gk +1 - хк +1. (16)

Дискретное значение хк + 1 является реакцией объекта на управляющее воздействие ик. Допустим, что ик формируется на основе выражения (8) метода Ньютона-Раф-сона, тогда при условии, что 8к ® 0, модель (1) записывается в виде

хк +1 = Ахк +

(17)

В выражения (12) и (13) входит величина ^, которую рекомендуется определять следующим образом:

8к = - ик -1.

Здесь у- это к-е приближение к решению уравнения (7), вычисленное с помощью любого итерационного метода.

Применение процедуры первого порядка (например, Ньютона-Рафсона) к определению Ук позволяет записать, что

у, = ик-! -[ /к'-! ]-1 /к-.; (14)

8к =-[ /к'-1 ]-1 /к-1. (15)

Подстановка dk в форме (15) в выражение (12) приводит к рекуррентной процедуре (10), т. е. к методу Чебышева. Аналогично можно показать, что алгоритм Хэлли

(11) является частным случаем итерационного процесса вида (13).

+В [фк -1 +фк -1 8к + 0,5 фк -1 82 + °(83)].

Здесь и в дальнейшем приняты обозначения фк1 -1 =ф(1) (ик -1), I = 0, 1, 2, 3.

Разность управлений 8, определяется по формуле (15) и с учетом выражений (5) и (6) вычисляется следующим образом:

8 к = [фк - 1 ] 1 [ В+( gk +1 - Ахк )-фк - 1 ] . (18)

После подстановки 8к в уравнение (17) и несложных преобразований получается

х, +1 = gk +1 + 0,5 В фк -182 + О (83). Динамическая ошибка (16) при реализации в системе алгоритма управления (8) на основе метода первого порядка Ньютона-Рафсона определяется выражением

= -0,5 В ф, -1 82 + О (83 ) (19)

или приближенным равенством

4 +1 »-0,5 В ф, -1 8,, которое получено из формулы (19) при 8, ® 0 .

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

Duk = uk - uk -1-

+B

xt +1 = Axk + jk -1 + jk -1 Duk + 0,5 jk -1 Dui +

+1 jk3-1Du3 + 0 (Du4)

+B

xk+1 = gk+1 +

0,5 jk -1 Aul - 0,5 jk -1 52 +

+— j*^ Au3, + (

6

-о (дм;)

Представим Ди2 и Duk в виде

Д«2 = 52 -j -, ]-1 jk -, 52 + о (54) Дик =52 + о (54),

тогда

+B

Введем обозначение

ßii = 66 j*-1 - 0,5 (jk -1 )-1(jk -1)2-

,(22)

’- B ß22 53.

Auk =[jk -1 + 0,5 jk -1 5k ] X[B+ (gk +1 - Axk )-jk -1] ,

+B

По аналогии с записью (17) представим уравнение объекта (1) при Дик ® 0 в форме

0,5 jk Aul - 0,5 jk -1 52 + +0,25( jk - 1)-1(jk -1)2 5k +

+6 j*-1 Auk + 0 (Au 4)

(29)

(21)

Представим выражения Дм2 и Ди\, используя для Дик формулу (28), в следующем виде:

Ди2 = 52-[jk -1 ]-1 jk -152 + о (54);

Au, = 5, +

(22)

Оценим точность системы при формировании управления ик в соответствии с выражением (12) метода Чебышева. Разность управлений Ди, записывается в виде

Дик =[фк -1] *х

х[ В+ +1 - Ахк )-ф,-1 - 0,5фк -1 52 ],

где 8, определяется выражением (18).

Формула (22) с учетом равенства (18) переписывается как

Дик =8к- 0,5[фк -1]1 фк -152. (23)

После подстановки в уравнение (21) только Ди,, определяемого выражением (22), получается, что

O(54),

тогда состояние объекта хк + 1 записывается в форме

+B

6jk3-1 -0,25(jk-1 )-1 (jk-1)2]53 + 0(54). (30)

По аналогии с выражением (25) вводится

ß 21 = 6 j?-1 - 0,25 (jk -1 )-1(jk -1)2

(31)

и для вектора динамической ошибки при 8, ® 0 справедливо приближенное равенство

4+1 »-В в21 5,. (32)

Сравнительный анализ методов второго порядка предлагается осуществлять следующим образом. Допустим, что в системе реализован идеальный алгоритм управления второго порядка с использованием 82 (примеры таких итерационных методов приводятся в [9]), тогда ошибка такой квадратичной аппроксимации определяется приближенным равенством

(2)

1

’- 6 B jk3^1 53. 6

6 j^ - 0,5 (jk-1 )-1 (jk -1)2 ^53 + O (54 ). (24)

Точность системы при использовании в устройстве управления алгоритма Чебышева (при 8, ® 0 ) оценивается приближенным равенством

Допустим, что в системе реализуется алгоритм Хэл-ли, тогда

*-к +^ ^'Тк - 1”к- (33)

Введем для оценки степени близости приближенных алгоритмов к идеальной системе второго порядка разности

Де(21) = _(2) - Г21) и Де(22) = е(2) - Г22)

^к + 1 кк +1 Кк + 1й ^Кк + 1 Кк + 1 Кк +1-

Определим величину 83 из выражения (32) и подставим в формулы (31) и (26). После несложных преобразований получается

Д4+1 » 1,5[фк-1 фк3-1 ]-1 (фк-1 )2 е, (34)

Дек2+2)1»3 [фк -1 фк3-1 ]1 (фк -1 )2 4+1- (35)

Анализ выражений для дискретных значений приближенных ошибок (20), (26) и (32) при 8, ® 0 говорит о более высокой точности систем, в которых реализованы алгоритмы второго порядка. В методе Ньютона-Рафсо-на, в соответствии с формулой (20), динамическая ошиб-(27) ка ек +1 пропорциональна 82, а в методах второго поряд-

(25)

(26)

(22)

и е,

(21)

где 8, также определяется формулой (18).

Выражение (27) при 8, ® 0 представляется рядом Тейлора вида

Дик =8, - 0,5[фк -1 ] 1 фк -182 +

Г/ ч2Т'/ 42 , ч (28)

+0,25 (ф, -1) (ф, -1) 83 + О (84 ).

После подстановки в уравнение объекта (21) разности Дик в форме (27) получается

ка, согласно выражениям (26) и (32), ошибки г\ зависят от 52, т. е. на порядок меньше. Аналогичный вывод получен и для итерационных процедур нахождения корня нелинейного уравнения [9; 10].

Формулы для оценки норм векторов приведем в силу их компактности только для соотношений (34) и (35). Допустим, что

m1 = min |j'(u )|, m3 = min |j(3) (u )| тогда

1,5 М2

и M2 = max | j "(u )|

Ae(21) <-

‘2 He(2)

и De1

(22)

(2)

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

Пример. Рассмотрим объект управления первого порядка, описываемый разностным уравнением

хк+1 = ахк + ЬШик , х0 = 0. (36)

Выбраны следующие значения параметров: а = 0,6 и Ь = 0,9. Требуемая траектория движения - единичная функция, т. е. gk+1 = 1.

Характеристика НЭ ф(и) = Ш и принадлежит классу ф(и) О С и для нее выполняется свойство диффеоморфизма [2], поэтому существует аналитическое решение, и управляющее воздействие определяется по формуле

щ = агсШ [г1 (gk +1 - ах)], и0 = 0 .

В замкнутой системе переход из нулевого состояния х0 = 0 в конечное х1 = 1 будет происходить за один шаг.

Рассмотрим применение приближенных алгоритмов. Построим в соответствии с формулой (5) функцию /(и) для исследуемого объекта в виде

/ (и ) = Ь~‘ (gk+1- ахк)-Ли. (37)

Допустим, что управляющее воздействие ик будем формировать на основе итерационного метода Ньюто-на-Рафсона, поэтому, используя выражение (8), получим

ик = ик-1 + Ь~‘ с^ ик-1 [gk +1 -ахк -Ь*ик-1],

и0 = 0 . (38)

Алгоритмы второго порядка могут быть получены непосредственно из формул (10) и (11) либо из процедур

(12) и (13), в которых для исследуемого объекта величина 8, равна

5, = Г1 Л2 щ-1 +1 -ахк -ЬШщ-1].

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

вого состояния в конечное при использовании алгоритмов второго порядка происходит на один шаг быстрее.

Экспериментальные исследования для проверки адекватности полученных аналитических выражений проводились следующим образом. Затухание переходных процессов во всех системах происходит до пятого шага, поэтому при к = 5 происходит кратковременное возмущение состояния объекта на величину Дх. В системах в соответствии с алгоритмами (8), (12) или (13) вырабатывается управляющее воздействие ик, которое компенсирует появившееся отклонение Дх практически за два шага (см. рисунок). Теоретические значения ошибок рассчитывались для метода Ньютона-Рафсона по формуле (20), а алгоритмов Чебышева и Хэлли - соответственно по выражениям (26) и (31), в которых использовались 5к, полученные в результате моделирования. Экспериментальные ошибки вычислялись непосредственно по формуле (16).

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

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

Моделирование проводилось и для других значений отклонений состояния объекта и при Дх = 0,7 (см. табл. 2).

Произошло ухудшение показателя точности во всех системах и возросло расхождение между экспериментальными и теоретическими значениями, так как ошибки зависят от 8к и необходимо, чтобы 8к ® 0, а при больших Дх это условие не выполняется.

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

Переходные процессы в системах при различных методах управления

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

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

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

Библиографический список

1. Методы классической и современной теории автоматического управления : в 5 т. Т. 5. Методы современ-

ной теории автоматического управления / под ред. К. А. Пупкова, Н. Д. Егупова. М. : Изд-во Моск. гос. техн. ун-та им. Н. Э. Баумана, 2004.

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

3. Аграчев, А. А. Геометрическая теория управления / А. А. Аграчев, Ю. Л. Сачков. М. : Физматлит, 2005.

4. Нестационарные системы автоматического управления: анализ, синтез и оптимизация / под ред. К. А. Пупкова, Н. Д. Егупова. М. : Изд-во Моск. гос. техн. ун-та им. Н. Э. Баумана, 2007.

5. Ким, Д. П. Теория автоматического управления: в 2 т. Т. 2. Многомерные, нелинейные, оптимальные и адаптивные системы / Д. П. Ким. М. : Физматлит, 2004.

6. Крутько, П. Д. Обратные задачи динамики управляемых систем. Нелинейные системы / П. Д. Крутько. М. : Наука, 1988.

7. Рубан, А. И. Адаптивное управление с идентификацией / А. И. Рубан. Томск : Изд-во Том. ун-та, 1983.

8. Бахвалов, Н. С. Численные методы / Н. С. Бахвалов. М. : Наука, 1973

9. Трауб, Дж. Итерационные методы решения уравнений / Дж. Трауб. М. : Мир, 1985.

10. Когут, А. Т. Полиномиальная аппроксимация в некоторых задачах оптимизации и управления / А. Т. Когут. Омск : Ом. гос. ун-т путей сообщения, 2003.

Таблица 1

Теоретические и экспериментальные ошибки системы для различных алгоритмов управления при Ах = 0,1

k Метод Ньютона-Рафсона е® Метод Хэлли е(21) Метод Чебышева е/(22)

Теоретические Эксперимен- Теоретические Эксперимен- Теоретические Эксперимен-

тальные тальные тальные

6 2,3-10-3 2,2-10-3 1,4-10-4 1,3-10-4 2,2-10-4 2,2-10-4

7 2,0-10-3 2,0-10-3 1,1-Ю-4 1,2-10-4 1,6-10-4 1,7-10-4

8 2,3-10-6 2,3-10-6 1,3-10-11 1,4 -10-11 6,7-10-11 6,7-10-11

9 9,2-10-7 9,3-10-7 2,6-10-13 2,6-10-13 1,1-10-12 1,1-10-12

10 5,3-10-15 5,1-10-15 0 0 0 0

Таблица 2

Теоретические и экспериментальные ошибки системы для различных алгоритмов управления при Ак = 0,7

k Метод Ньютона-Рафсона е® Метод Хэлли е(21) Метод Чебышева е/(22)

Теоретические Эксперимен- тальные Теоретические Эксперимен- тальные Теоретические Эксперимен- тальные

6 0,109 2 0,073 5 0,047 1 0,034 1 0,075 2 0,061 8

7 0,033 8 0,028 5 0,019 6 0,019 0 0,012 9 0,013 3

8 1,38-10-6 1,38-Ю"6 7,13 -10-5 7,85-10-5 1,55-Ю-4 1,66-10-4

9 1,92-10-4 1,90-10-4 9,53-10-7 9,30-10-7 4,87-10-7 4,83-10-7

10 5,67-10-8 5,67-10-8 7,30-10-14 7,33-10-14 1,03-10-12 1,04-10-12

A. T. Kogut, N. Yu. Panfilova

FORMATION OF ALGORITHMS OF OPERATING INFLUENCES ON THE BASIS OF NUMERICAL METHODS

Possibility and technique of synthesis recurrent algorithms of control by nonlinear objects in the form of iterative procedures ofnumerical methods is considered. Accuracy of the closed systems is investigated at realisation in regulators offered methods. The comparative analysis of the received analytical expressions and results of imitating modelling, and also synthesised algorithms is carried out.

Keywords: object, algorithm, management, accuracy, method, equation.