CC BY
39
Оптимальное интерполирование типовых динамик в задаче управления с прогнозированием
Н. И. Гданский, А.В. Карпов, А.А. Бугаенко
Основным путем повышения эффективности оборудования является автоматизация основных и вспомогательных производственных операций. Выполнение последних, как правило, сопровождается недетерминированным изменением внешней нагрузки на приводах. В работе [1] предложено в таких случаях для управления перемещением по заданному закону в системе с одной степенью свободы ф использовать прогнозирование недетерминированной внешней нагрузки, представленной в виде сложной функции М(^,ф(0) по времени t. Ее математическую модель представлена в виде скалярного произведения М(^ф(^) = ( Мк, ф(0Х в котором Мк - состоящий из частных производных по t и ф силовой вектор порядка к (максимальный порядок производных), ф(^) - соответствующий ему кинематический вектор.
В качестве промежуточной величины для расчета управляющего воздействия предложено использовать работу А, которую должен совершать привод на заданном временном отрезке управления ^г, tг+1\. Для принятой модели внешней нагрузки формула для
А как функции времени I принимает общий вид:
Ч+1 __
А = ](М *, ф4 (0)Ф'(* )йг. (1)
Практически малые угловые перемещения регистрируют при помощи инкрементных датчиков, выдающих дискретные импульсы при угловом смещении вала на фиксированный шаг, величину которого обозначим через к. Поэтому для сокращения формульных выражений при данном способе измерения угловых перемещений наряду с обычной работой рассмотрим приведенную работу, равную обычной, деленной на шаг к: Д = Д /к .
Безразмерные инкрементные угловые перемещения вала обозначим через у, а время (^ - ^) относительно начала измерений - через т.
Для модели нагрузки первого порядка (к =1) выражение для приведенной работы
Д через вектор Мк имеет вид:
А = М(¥г+1 ~¥г ) + ~ | ¥ (т) • т • Лт+— I (У (т)) ■ Т ■ йт. (2)
т1 г Тг
При к =2 выражение Д( М2, у2(т)) имеет вид:
- дМ г дМ г
А = М(у+1 -¥г) + — I ¥(т) • т • йт + — I [(у'(т)) • т + у'(т) • у"(т) •т ] • йт +
— J ду •!
т т т
+1-2М ]¥'(т) •т2 • йт + — т|1(у'(т))2 •т2 • йт+1-2М т|+1(¥'(т))3 •т2 • йт. (3)
2 дт2 I дтду у У ” 2 ду2 }
Исходную информацию для расчета управляющего воздействия задают следующие данные обратной связи: 1) массив моментов времени {tг} с зарегистрированными сигналами инкрементного датчика; 2) массив значений работ {Аг} на отрезках [^, ^]. Для получения из формул Д ( М, ук(т)) линейных уравнений относительно коэффициентов силового вектора Мк на каждом отрезке управления [тг, тг+1] в эти зависимости должно быть подставлено уравнение некоторой кривой уг(т), задающей перемещение управляемого звена внутри данного отрезка.
В работе [1] использована простейшая кусочная интерполяция траектории перемещения объекта управления при помощи сплайнов Лагранжа - отрезков прямых, проходя-
щих через узловые точки. Однако численные расчеты показывают, что точность расчета интегралов в выражениях (2) и (3) при этом невысока и для получения лучшего приближения зависимости Д. ( Мк, ук(т)) необходимо использовать другие виды степенных
сплайнов уг(т) более высокого порядка по т. Помимо численной близости к точным значениям интегралов в (2) и (3), второй существенной качественной характеристикой сплай-нового приближения является трудоемкость алгоритмов расчета коэффициентов сплайнов. Она особенно актуальна при использовании в качестве вычислительных устройств в системах управления микроконтроллеров.
Интегралы в (2) и (3) представим их в обобщенном виде т+1
\тр • (¥(т)У • (¥'(т)Уйт = 1мК1). (4)
т
в котором степенир,ц и г -неотрицательные целые числа.
При данном обозначении формулы (2) и (3) принимают вид:
1. к =1:
\ дМ дМ
Д = М(Щ+! — Щ ) + • 1и((1) + • 1120(1). (5)
2. к =2:
Ч дМ дМ 1 д 2М д 2М 1 д 2М
Д = М (Щ+1 Щ1 ) + ^ • !ис(1) + (1120(1) + 1211(1) ) + 2 т ■1210(1) + ту 1220.1) 2 Щ2 ' 1230(1). (6)
Таким образом, качество интерполирования величин работ (2) и (3) определяется
точностью, с которой сплайны приближают следующие 6 интегралов, которые назовем
траекторными:
1110(г), 1120(г), 1211(г), 1210(г), 1220(г), 1230(г). (7)
Поскольку линейные сплайны дают невысокую точность приближения, для вычисления траекторных интегралов (7) рассмотрим возможности квадратичных и кубических локальных сплайнов, имеющих порядок 2 и 3, а также кубических интерполяционных сплайнов Фергюссона.
Отрезок [тг, тг+1], на котором строится интерполирующий сплайн, для краткости назовем интерполируемым. Наряду с ним рассмотрим также соседние с ним отрезки [тг.1, тг] и [тг+1, тг+2], которые назовем окрестностью интерполируемого отрезка.
Общее уравнение сплайна второго порядка для отрезка [тг, тг+1] для более удобного выполнения преобразований представим в виде: у^) = Сод + Сод т + Сод т2/2, где Сод, Сод, Сод - константы. К данному типу сплайнов отнесём:
- сплайн по усредненным первым производным в крайних точках отрезка интерполирования: С2(г) = (уг'(тг+1)- ^г'(тг)) / к Сод = ^г'(тг) - Сод^; V'г = (Дуг-1 + Ауг)/(/г-1 + /г) V'г+1 = (Ауг+Ауг+1)/(/г+ /г+1>),
- сплайн по взвешенным первым производным в крайних точках отрезка интерполирования: у ' гв = (Ауг-1+2Ауг)/(/г-1+ 21); у ' (г+1)в = (2Ауг+Ауг+1)/(2/г+ /г+1)).
Общее уравнение сплайна порядка 3 для отрезка [тг, тг+1] представим в виде:
^/(0 = С0(г) + С1(г) т + Сод т2/2 + С3(г) т3/3, где С0(г), Сх(г), С2(г), С3(г) - константы.
К таким сплайнам относятся:
- сплайны Фергюссона[2],
- сплайн по усредненным первым производным в крайних и средней точке отрезка интерполирования: С3(г)=(2у'срг - у 'г+1 - у 'г)/(2тсрг2 - тг+12- т,2) = 2(у'г+1 + у'г - 2у 'срг)//2г/ С2(г) = 2[у'срг
- у'г - С3(г) (т срг 2 - т? )]//; Сод =у'г - Сод^ - С3(гУГ2; у 'г = (Ауг-1+Ауг)/(/г-1+ /г); у'г+1 =
(Ауг+Ауг+1)/(/г+ /г+1); у'срг= А^г / /г),
- сплайн по взвешенным первым производным в крайних точках отрезка интерполирования: у 'г = (Ауг-1+2Ауг)/(/г-1+ 2/г); у 'г+1 = (2Ауг+Ауг+1)/(2/г+ /г+1); у 'срг = А^г / /г..
Введем следующие обозначения:
- т-вектор частных производных, входящих в силовой вектор М порядка к,
~к
- р - длина вектора т ,
- I (г) , їпр(і) - векторы точных интегралов при коэффициентах вектора гт для отрезка [т;,тг+і] и приближенных интегралов, соответствующих выбранному методу интерполирования траектории.
Точное значение приведенной работы на отрезке [тг,тг+1] равно:
Д = М (У,+1 - V,) + (тк, (8)
Приближенное значение приведенной работы на отрезке [тг-,тг+1] при выбранном методе интерполирования на заданной траектории выражается:
Арі,) = М(У,+і - V,) + (шк, 1«Р,)). (9)
Разность точного и приближенного значений работы на отрезке [тг-,тг+1]:
АД = Д - Дпр(і) = (тк,(І(І) - /Ир(0)).
Используя неравенство Коши — Буняковского, оценим абсолютную величину разности работ, используя норму векторов в евклидовом пространстве:
|ад 1=|(тк - 4,(1))| < тл • 14 - и, (10)
І -І =
где I ( 0 пр(°1 -у
р
X (І(І)Ї - Іпр(,)] )2. ]=1
Оценим разность работ для совокупности п подряд идущих отрезков времени [хо=0,Т1], [т1,т2],...,[тп_1,тп]. Так как при построении сплайнов наряду с интерполируемым отрезком необходимо использовать данные о функции на его окрестности, то суммирование применим по всем отрезкам кроме крайних с номерами 0 и (п -1):
п-2
АД =
п-2 п-2
<£ аД < |тк| • X
X (І(г)] - Іпр(г)] )2. (11)
]=1
Полученная оценка точности расчета приведенной работы для всей траектории зависит от метода интерполирования (обозначим его МлЩ, формы траектории (которую обозначим Тг) и модуля вектора точных значений интегралов /(*), от вектора частных производных шк и порядка модели к: Выделим в полученной формуле величину, зависящую преимущественно от метода ее интерполирования, формы траектории и порядка модели к:
п-2
АД
т1 •К(,)1
,=1 < --------
Х (І(-)] 1пр(і)] )
]=1 -----= ЬД(Тт, МґИ, к).. (12)
,=1
X (1(г)] )2 ]=1
м
Назовем величину -— !_ . относительной погрешностью расчета полной рабо. И • к «I
ты, а ЬЛ(ТгМ1Ь,к) - гарантированной точностью ее расчета. Примем последнюю в качестве оценочного критерия влияния метода интерполирования на относительную погрешность расчета работы при выбранном методе интерполирования М№ на траектории Тг. Г арантированное качество выполнения интерполирования траектории заданным методом интерполирования возрастает с уменьшением величины 5Л(ТгМ^,к).
Практически в методах управления с прогнозированием внешней нагрузки необходимо использовать модели нагрузки порядков от 0 до максимального значения к.=2 Поскольку в начало вектора коэффициентов М модели нагрузки максимального порядка к входят все коэффициенты векторов М1 всех моделей нагрузки более низких порядков
,=1
г=1
г=1
q<k, то при минимальной величине 5Л(Тг,Мй,к), будут близки к минимальным значения 5Л(Тг,Мй^). Поэтому для оценки качества интерполирования достаточно определить оптимальный метод интерполирования для максимального порядка к = 2.
В качестве тестовых рассмотрено 9 кривых. На каждой из них выбрано 16 точек (уг=/,тг),(/=0,...,и=15), которые моделируют инкрементный способ измерения перемещения с шагом 1 на 15, при котором выполняется условие: уг+1 - уг- = 1,(/ = 0, ... , п-1=14). Тестовыми являются следующие типовые виды перемещений:
1. Линейное перемещение: у1(т)=т ; Т(1);=/.
2 /7
2. Квадратичное перемещение: у2(т)=т ; ^(2)г г .
з П
3. Кубическое перемещение: у3(т) =т3; Т(3)г=л'г .
4 П
4. Парабола 4 степени: у4 (т)=т4; Т(4)г=л/г.
5. Перемещение по графику квадратного корня: у5(т) = +1; 1= /2-1.
6. Перемещение по графику кубического корня: у6(т) = 3т +1; Т(6)/= /3-1.
7. Перемещение по графику корня 4 степени: у7 (т) = 4т +1; т(7)/ = /4-1.
8. Перемещение по кубической параболе с точкой перегиба: у8 (т) = т2(8-т/3)/45 .
9. Перемещение по синусоиде: у9 (т) = 16-вт(т); Т(9)/= агс8т(//16).
Вектор интегралов (точных и приближенных) I при коэффициентах вектора т следующим образом выражается через траекторные интегралы:
I = (1110, 1120+1211, 1210, 1220, 1230). (13)
Формула для гарантированной точности 5Л(Тг,М^,2) расчета относительной приведенной работы в модели нагрузки второго порядка траектории Тг5 с номером ^ с использованием метода ее интерполирования Ы№г с номером г принимает вид:
13 Г5
(1 (г)] ~ 1пр(,); )2
М(Гг,,ЫАГ ,2) = —У ^
1
(14)
і =
где для выбранной траектории Тгі.
/(/■у - точное значение интеграла при коэффициенте вектора Ш2 с номером у для участка Тгц с номером І,
/ пр(і)] - аналогичное приближенное значение интеграла при коэффициенте вектора ш с номером у для участка Тг5 с номером /, соответствующее методу интерполирования Мйг.Результаты расчетов приведены в таблице 1.
Табл.1. Точность относительной погрешностью расчета полной работы
Номер траектории, Метод интерполирования 1 2 3 4 5 6 7 8 9 8е
Квадратичные сплайны
1. Сплайн по усредненным производным в крайних точках о 0,006988 0,007808 0,008549 0,002265 0,003655 0,003820 0,03561 0,02073 0,08942
2. Сплайн по взвешенным производным в крайних точках о 0,01874 0,02531 0,02849 0,001804 0,003191 0,003067 0,03068 0,02856 0,1398
Кубические сплайны
3. Сплайны Фергюссона 0 0,0000583 0,001715 0,0005241 0,01411 0,000677 0,007384 0,002843 0,00013 0,02744
4. Сплайн по усредненным производным в крайних и средней точке 0 0,001840 0,00263621 0,0027363 0,00070610 0,00120691 0,00146669 0,013165 5 6 6 7 Г'' 0 0, 0, 0,03152373
5. Сплайн по взвешенным производным в крайних точках отрезка 0 0,017588 0,00263336 0,0026575 0,00130052 0,00261488 0,00306716 0,0161371 0,0241498 0,07014868
Сравнение качества приближения по отдельным кривым показывает, что:
1) кубические интерполяционные сплайны Фергюссона (№3) дают лучший результат для 6 кривых из 8 рассмотренных,
2) кубические локальные сплайны по усредненным производным в крайних и средней точке (№4) дают лучший результат для 2 кривых из 8.
Анализ суммарной погрешности 5е показывает, что наилучшие результаты дает использование кубических сплайнов:
1) интерполяционные сплайны Фергюссона (№3), 5е = 0,02744,
2) локальные сплайны по усредненным производным в крайних и средней точке (№4), 5е = 0,0315237,
3) локальные сплайны по взвешенным производным в крайних точках отрезка (№5), 5е =
0,07014868.
Рассмотрим трудоемкость вычисления коэффициентов данных сплайнов. Поскольку вычисление коэффициентов у кубических локальных сплайнов по взвешенным производным в крайних точках отрезка (№5) более трудоемко, чем у аналогичных сплайнов (№4), а точность хуже, то данные сплайны из анализа исключаем.
Отличительной особенностью используемых сплайнов является то, что их коэффициенты должны определяться для глобальной переменной.
Расчет п сплайнов Фергюссона для локальных переменных требует выполнения: 9п-3 сложения, 8п-3 умножений, 4п-2 деления. Однако при переходе к каноническому виду по глобальной переменной суммарные вычислительные затраты резко вырастают до: 28п-3 сложения, 36п-3 умножений, 6п-2 деления. Прямой расчет данных коэффициентов без использования полиномов Эрмита позволяет снизить вычислительные затраты до: 16п-18 сложений, 12п-5 умножений, 4п-2 деления.
Рассмотрим расчет п кубических локальных сплайнов по усредненным производным в крайних и средней точке (№4). Поскольку в расчетных формулах для их коэффициентов повторяется ряд величин, их целесообразно рассчитать до начала основного расчета:
1) Т(г)кв=Тг2, (г = 0,1,...,п+1); 2) Дщ г = щ - Щ-1, (г = 0,1,...,п+1);
3) Iг = Тг- - Тг-1, (/ = 0,1,...,п+1); 4) щ’г = (Дуг--1+ДУг)/(1г-1+ I) (*' = 1,...,п+1) .
Расчетные формулы для каждого сплайна с номером г (г - ,1,...,п):
у срг Дуг / 1г ; ^срг 0,5(^г + ^г-1); С3(г) 2(у г+1 + у г - 2у срг)/1 г;
С2(г) 2[у срг - у г - С3(г) (^ срг - ^(г)кв)]/1г; С1(г) у г - С2(г)^г - С3(г^(г)кв .
Суммарные вычислительные затраты составляют: 10п+4 сложения, 12п+3 умножения, 4п +1 деление. Существенное снижение по сравнению с прямым расчетом коэффициентов сплайнов Фергюссона получено только по сложению (6п). Число остальных операции с точностью до постоянных слагаемых одинаково. Поэтому оптимальным вариантом в задачах управления с прогнозированием является применение интерполяционных сплайнов Фергюссона с использованием прямого расчета их коэффициентов по глобальной переменной.
Список литературы:
1.Марченко Ю.А. Адаптивный цифровой алгоритм программного управления в условиях переменной внешней нагрузки. // Химическое и нефтегазовое машиностроение. №12, 2010. - с. 34-36.
2. Гданский Н.И. Геометрическое моделирование и машинная графика. - М.: МГУИЭ,
2003 г. - 236 с.
