×

Вы используете устаревший браузер Internet Explorer. Некоторые функции сайта им не поддерживаются.

Рекомендуем установить один из следующих браузеров: Firefox, Opera или Chrome.

Контактная информация

+7-863-218-40-00 доб.200-80
ivdon@ivdon.ru

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

Аннотация

Н. И. Гданский, А.В. Карпов, А.А. Бугаенко

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

Ключевые слова: Траектория перемещения, сплайны Лагранжа, сплайны Фергюссона, интерполяция

05.13.06 - Автоматизация и управление технологическими процессами и производствами (по отраслям)

Основным путем повышения эффективности оборудования является автоматизация основных и вспомогательных производственных операций. Выполнение последних, как правило, сопровождается недетерминированным изменением внешней нагрузки на приводах. В работе [1] предложено в таких случаях для управления перемещением по заданному закону в системе с одной степенью свободы φ использовать прогнозирование недетерминированной внешней нагрузки, представленной в виде сложной функции M(t,φ(t)) по времени t. Ее математическую модель представлена в виде  скалярного произведенияM (t,φ(t)) = (¯M k,¯φk(t)), в котором ¯M k - состоящий из частных производных по t и φ  силовой вектор порядка k (максимальный порядок производных), ¯φk(t) - соответствующий ему кинематический вектор.
В качестве промежуточной величины для расчета управляющего воздействия предложено использовать работу А, которую должен совершать привод на заданном временном отрезке управления [ti, ti+1]. Для принятой модели внешней нагрузки формула для А как функции времени t принимает общий вид:
.  (1)      
Практически малые угловые перемещения регистрируют при помощи инкрементных датчиков, выдающих дискретные импульсы при угловом смещении вала на фиксированный шаг, величину которого обозначим через h. Поэтому для сокращения формульных выражений при данном способе измерения угловых перемещений наряду с обычной работой рассмотрим приведенную работу, равную обычной, деленной на шаг h:    Безразмерные инкрементные угловые перемещения вала обозначим через y, а время (t -t0) относительно начала измерений - через t.
Для модели нагрузки первого порядка (k =1) выражение для приведенной работы  через вектор ¯M k имеет вид:
    (2)
При  k =2  выражение  (¯M 2,¯ Ψ2(τ))  имеет вид:
         
   (3)
Исходную информацию для расчета управляющего воздействия задают следующие данные обратной связи: 1) массив моментов времени {ti} с зарегистрированными сигналами инкрементного датчика; 2) массив значений работ {Аi} на отрезках [ti,ti+1]. Для получения из формул (¯Mk,¯Ψk(τ)) линейных уравнений относительно коэффициентов силового вектора¯M k на каждом отрезке управления [ti, ti+1] в эти зависимости должно быть подставлено уравнение некоторой кривой Ψi(τ), задающей перемещение управляемого звена внутри данного отрезка.
В работе [1] использована простейшая кусочная интерполяция траектории перемещения объекта управления при помощи сплайнов Лагранжа - отрезков прямых, проходящих через узловые точки. Однако численные расчеты показывают, что точность расчета интегралов в выражениях (2) и (3) при этом невысока и для получения лучшего приближения зависимости (¯Mk,¯Ψk(τ)) необходимо использовать другие виды степенных сплайнов Ψi(t) более высокого порядка по  t. Помимо численной близости к точным значениям интегралов в (2) и (3), второй существенной качественной характеристикой сплайнового приближения является трудоемкость алгоритмов расчета коэффициентов сплайнов. Она особенно актуальна при использовании в качестве вычислительных устройств в системах управления микроконтроллеров.
Интегралы в (2) и (3) представим их в обобщенном виде
   (4)
в котором степени p,q и r -неотрицательные целые числа.
При данном обозначении формулы (2) и (3) принимают вид:
1. k =1:
     (5)
2. k =2:
  (6)
Таким образом, качество интерполирования величин работ (2) и (3) определяется точностью, с которой сплайны приближают следующие 6 интегралов, которые назовем траекторными:
I110(i), I120(i), I211(i), I210(i), I220(i), I230(i).  (7)
Поскольку линейные сплайны дают невысокую точность приближения,  для вычисления траекторных интегралов (7) рассмотрим возможности квадратичных и кубических локальных сплайнов, имеющих порядок 2 и 3, а также кубических интерполяционных сплайнов Фергюссона.
Отрезок [τii+1], на котором строится интерполирующий сплайн, для краткости назовем интерполируемым.  Наряду с ним рассмотрим также соседние с ним отрезки [τi-1i] и [τi+1i+2], которые назовем окрестностью интерполируемого отрезка.
Общее уравнение сплайна второго порядка для отрезка [τi, τi+1] для более удобного выполнения преобразований представим в виде: Ψi(t) = C0(i) + C1(i) τ+ C2(i) τ2/2, где C0(i), C1(i), C2(i)  - константы. К данному типу сплайнов отнесём:
- сплайн по усредненным первым производным в крайних точках отрезка интерполирования: C2(i) = (Ψi’(τi+1)- Ψi’(τi)) / li; C1(i)  = Ψi’(τi) - C2(i)τi; Ψi’ =(ΔΨi-1+ΔΨi)/(li-1+ li); Ψi’i+1 =(ΔΨi+ΔΨi+1)/(li+ li+1)),
- сплайн по взвешенным первым производным в крайних точках отрезка интерполирования: Ψ’ =(ΔΨi-1+2ΔΨi)/(li-1+ 2li); Ψ’(i+1)=(2ΔΨi+ΔΨi+1)/(2li+ li+1)).
Общее уравнение сплайна порядка 3 для отрезка [τi, τi+1] представим в виде:
Ψi(t) = C0(i) + C1(i) τ+ C2(i) τ2/2 + C3(i) τ3/3,     где C0(i), C1(i), C2(i), C3(i)   - константы.
К таким сплайнам относятся:
- сплайны Фергюссона[2],
-  сплайн по усредненным первым производным в крайних и средней точке отрезка интерполирования: C3(i)=(2Ψ’срi - Ψ’i+1 - Ψ’i)/(2τсрi2 - τi+12- τi2) = 2(Ψ’i+1 + Ψ’i - 2Ψ’срi)/l2i; C2(i)= 2[Ψ’срi -  Ψ’i  - C3(i)срi 2 - τi2 )]/li; C1(i) =Ψ’- C2(i)τi - C3(i)τi2; Ψ’i  =(ΔΨi-1+ΔΨi)/(li-1+ li); Ψ’i+1 =(ΔΨi+ΔΨi+1)/(li+ li+1); Ψ’срi= ΔΨi / li),
- сплайн по взвешенным первым производным в крайних точках отрезка интерполирования: Ψ’i   =(ΔΨi-1+2ΔΨi)/(li-1+ 2li); Ψ’i+1 =(2ΔΨi+DΨi+1)/(2li+ li+1); Ψ’срi = ΔΨi / li..
Введем следующие обозначения:
-¯mk-вектор частных производных, входящих в силовой вектор¯Mk порядка k,
- р – длина вектора¯mk,
-¯I (i) ,¯Iпр(i) - векторы точных интегралов при коэффициентах вектора¯mk для отрезка [τ ii+1] и приближенных интегралов, соответствующих выбранному методу интерполирования траектории.
Точное значение приведенной работы на отрезке [τ ii+1] равно:
   (8)
Приближенное значение приведенной работы на отрезке [τ ii+1] при выбранном методе интерполирования на заданной траектории выражается:
  (9)
Разность точного и приближенного значений работы на отрезке [τ ii+1]:

Используя неравенство Коши — Буняковского, оценим абсолютную величину разности работ, используя норму векторов в евклидовом пространстве:
      (10)
где
Оценим разность работ для совокупности n подряд идущих отрезков времени [τ0=0,τ1], [τ12],...,[τn-1n]. Так как при построении сплайнов наряду с интерполируемым отрезком необходимо использовать данные о функции на его окрестности, то суммирование применим по всем отрезкам кроме крайних с номерами 0 и (n -1):
    (11)
Полученная оценка точности расчета приведенной работы для всей траектории зависит от метода интерполирования (обозначим его Mth),  формы траектории (которую обозначим Tr) и модуля вектора точных значений интегралов I (i), от вектора частных производных¯mk и порядка модели k: Выделим в полученной формуле величину, зависящую преимущественно от метода ее интерполирования, формы траектории и порядка модели k:
    (12)
. Назовем величину относительной погрешностью расчета полной работы, а δA(Tr,Mth,k) - гарантированной точностью ее расчета. Примем последнюю в качестве оценочного критерия влияния метода интерполирования на относительную погрешность расчета работы  при выбранном методе интерполирования Mth на траектории Tr. Гарантированное качество выполнения интерполирования траектории заданным методом интерполирования возрастает с уменьшением величины δA(Tr,Mth,k).
Практически в методах управления с прогнозированием внешней нагрузки необходимо использовать модели нагрузки порядков от 0 до максимального значения k.=2 Поскольку в начало вектора коэффициентов¯Mk модели нагрузки максимального порядка k входят все коэффициенты векторов¯M1 всех моделей нагрузки более низких порядков q<k, то при минимальной величине δA(Tr,Mth,k), будут близки к минимальным значения δA(Tr,Mth,q). Поэтому для оценки качества интерполирования достаточно определить оптимальный метод интерполирования для максимального порядка k = 2.
В качестве тестовых рассмотрено 9 кривых. На каждой из них выбрано 16 точек (Ψi=ii),(i=0,...,n=15), которые моделируют инкрементный способ измерения перемещения с шагом 1 на 15, при котором выполняется условие: Ψi+1 - Ψi = 1,(i = 0, ... , n-1=14). Тестовыми являются следующие типовые виды перемещений:

  1. Линейное перемещение: Ψ1(τ)=τ1; τ(1)i=i.
  2. Квадратичное перемещение: Ψ2(τ)=τ2; τ(2)i=.
  3. Кубическое перемещение:  Ψ3(τ) =t3;  τ(3)i= .
  4. Парабола 4 степени:  Ψ4 (τ)=τ4;  τ(4)i=.
  5. Перемещение по графику квадратного корня: Ψ5(τ) = ; i= i2-1.
  6. Перемещение по графику кубического корня:  Ψ6(τ) = ;  τ(6)i= i3-1.
  7. Перемещение по графику корня 4 степени:  Ψ7 (τ) = ; τ(7)ii4-1.
  8. Перемещение по кубической параболе с точкой перегиба:  Ψ8 (τ) = τ2(8-τ/3)/45 .
  9. Перемещение по синусоиде:  Ψ9 (τ) = 16×sin(τ);  τ(9)i= arcsin(i/16).

  10. Вектор интегралов (точных и приближенных)¯I при коэффициентах вектора¯mk следующим образом выражается через траекторные интегралы:
    ¯I = (I110, I120+I211, I210, I220, I230). (13)
    Формула для гарантированной точности δA(Tr,Mth,2) расчета относительной приведенной работы в модели нагрузки второго порядка траектории Trs с номером s с использованием метода ее интерполированияMthr с номером r принимает вид:
    (14)
    где для выбранной траектории Trs:
    I(i)j – точное значение интеграла при коэффициенте вектора¯m2 с номером j для участка Trs с номеромi,
    I пр(i)j – аналогичное приближенное значение интеграла при коэффициенте вектора¯m2 с номером j для участка Trs с номером i, соответствующее методу интерполирования Mthr.Результаты расчетов приведены в таблице 1.

 

Табл.1. Точность относительной погрешностью расчета полной работы


Номер траектории,
Метод интерполирования

1

2

3

4

5

6

7

8

9

δΣ

Квадратичные сплайны

1. Сплайн по усредненным производным в крайних точках

0

0,006988

0,007808

0,008549

0,002265

0,003655

0,003820

0,03561

0,02073

0,08942

2. Сплайн по взвешенным производным в крайних точках

0

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

0,0077665

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.
Анализ суммарной погрешности δΣ показывает, что наилучшие результаты дает использование кубических сплайнов:
1) интерполяционные сплайны Фергюссона (№3),  δΣ  = 0,02744,
2) локальные сплайны по усредненным производным в крайних и средней точке (№4),  δΣ  = 0,0315237,
3) локальные сплайны по взвешенным производным в крайних точках отрезка (№5),  δΣ  = 0,07014868.
Рассмотрим трудоемкость вычисления коэффициентов данных сплайнов. Поскольку вычисление коэффициентов у кубических локальных сплайнов по взвешенным производным в крайних точках отрезка (№5) более трудоемко, чем у аналогичных сплайнов (№4), а точность хуже, то  данные сплайны из анализа исключаем.
Отличительной особенностью используемых сплайнов является то, что их коэффициенты должны определяться для глобальной переменной.
Расчет n сплайнов Фергюссона для локальных переменных требует выполнения:  9n-3 сложения, 8n-3 умножений, 4n-2 деления. Однако при переходе к каноническому виду по глобальной переменной суммарные вычислительные затраты резко вырастают до: 28n-3 сложения, 36n-3 умножений, 6n-2 деления. Прямой расчет данных коэффициентов без использования полиномов Эрмита позволяет снизить вычислительные затраты до: 16n-18 сложений, 12n-5 умножений, 4n-2 деления.
Рассмотрим расчет n кубических локальных сплайнов по усредненным производным в крайних и средней точке (№4). Поскольку в расчетных формулах для их коэффициентов повторяется ряд величин, их целесообразно рассчитать до начала основного расчета:
1) τ(i)квi2,  (i = 0,1,...,n+1);  2) ΔΨi = Ψi - Ψi-1, (i = 0,1,...,n+1); 
3)  l i = τi - τi-1, (i = 0,1,...,n+1);  4)  Ψi’  =(ΔΨi-1+ΔΨi)/(li-1+ li) (i = 1,...,n+1) .
Расчетные формулы для каждого сплайна с номером i (i - ,1,...,n):
Ψ’срi = ΔΨi / li ;  τ срi = 0,5(τi + τi-1);    C3(i) = 2(Ψ’i+1 + Ψ’i - 2Ψ’срi)/l2i;
 C2(i)= 2[Ψ’срi -  Ψi’ - C3(i)срi 2 - τ(i)кв)]/li;      C1(i)i’ - C2(i)τi - C3(iτ(i)кв .
Суммарные вычислительные затраты составляют: 10n+4 сложения, 12n+3 умножения, 4n +1 деление.  Существенное снижение по сравнению с прямым расчетом коэффициентов сплайнов Фергюссона  получено только по сложению (6n). Число остальных операции с точностью до постоянных слагаемых одинаково. Поэтому оптимальным вариантом в задачах управления с прогнозированием является применение интерполяционных  сплайнов Фергюссона с использованием прямого расчета их коэффициентов по глобальной переменной.


Список литературы: 

  1.Марченко Ю.А. Адаптивный цифровой алгоритм программного управления в условиях переменной внешней нагрузки. // Химическое и нефтегазовое машиностроение. №12, 2010. – с. 34-36.
2. Гданский Н.И. Геометрическое моделирование и машинная графика. – М.: МГУИЭ, 2003 г. – 236 с.