CC BY
138
электронное научно-техническое издание
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС 77 - 30569. Государственная регистрация №0420900025. ISSN 1994-040S
Нейросетевые идентификация и управление механизмом параллельной кинематики типа «трипод» # 09, сентябрь 2010
авторы: Гоменюк С. М., Карпенко А. П.
УДК 519.6
[email protected] МГТУ им. Н.Э. Баумана
Введение
Серьезным недостатком классических многозвенных рычажных роботов-манипуляторов является их недостаточная жесткость, усложняющая управление этими манипуляторами, затрудняющая использование высокоэнергетического обрабатывающего инструмента и достижение высокой точности обработки. В значительной мере преодолеть указанные недостатки могут манипуляторы типа «хобот», построенные на основе многосекционных механизмов с параллельной структурой [1].
Манипуляторы типа «хобот» могут быть построены на основе механизмов параллельной кинематики с двумя, тремя, четырьмя и шестью степенями свободы [2]. В работе рассматривается манипулятор, секции которого представляют собой триподы -параллельные механизмы, обладающие двумя вращательными степенями свободы относительно двух осей, перпендикулярных продольной оси предыдущей секции манипулятора.
Управление многосекционным роботом-манипулятором является сложной задачей, в которой естественным образом выделяются стратегический, тактический и исполнительный уровни [3]. На стратегическом уровне выполняется планирование целевой конфигурации манипулятора, на тактическом уровне - синтез траектории перехода из текущего состояния манипулятора в целевое состояние, на исполнительном уровне - управление приводами секций манипулятора.
В работе рассматривается исполнительный уровень управления роботом-манипулятором типа «хобот», точнее говоря, рассматривается управление триподом, как секцией этого манипулятора.
Управление сложной динамической системой, которую представляет собой трипод, сопряжено с решением задачи идентификации этой системы [4]. Классические методы идентификации разделяются на методы непараметрической и параметрической идентификации.
При непараметрической идентификации чаще всего используют переходную характеристики системы, построенную на основе кривой разгона при ступенчатом изменении управления на входе системы. Затем эту характеристику аппроксимируют подходящим аналитическим выражением и на этой основе строят модельное дифференциальное уравнение. Одним из вариантов непараметрической идентификации являются частотные методы идентификации [5]. Параметрическая идентификация динамических систем основана на определении параметров передаточной функции, имеющей заданную структуру. В работе используется непараметрическая идентификация в среде модуля System Identification Toolbox программной системы MatLab [6].
Среди неклассических методов идентификации динамических систем в настоящее время наиболее популярной является идентификация с помощью многослойных нейронных сетей. В качестве нейросетевых модельных структур обычно используют следующие модификации линейных регрессионных моделей [7]: NNARX - нейросетевая авторегрессионная модель; NNARMAX - нейросетевая авторегрессионная модель скользящего среднего; NNSSIF -нейросетевая модель типа «обновления пространства состояний». Мы используем нейросетевую идентификацию на основе NNARX-модели, полученной средствами модуля Nеural Nеtwоrk Tооlbоx программного комплекса Matlab [8].
При функционировании трипода в составе многосекционного робота-манипулятора типа «хобот», динамические характеристики платформы трипода, а также силы и моменты, действующие на нее, могут изменяться в очень широких пределах. Поэтому система автоматического управления (САУ) триподом должна быть адаптивной [5].
В работе сначала рассмотрены классические и нейросетевые неадаптивные системы управления. В качестве классической САУ рассмотрена двухконтурная система на основе ПИД-регуляторов, в качестве нейросетевой САУ - одноконтурная система на
основе инверсной модели объекта управления [9]. Затем рассмотрены адаптивные САУ, построенные также на основе классических и нейросетевых систем управления.
Во всех случаях нейросетевые компоненты систем автоматического управления реализованы с использованием модуля Neural Network Toolbox.
Эффективность рассматриваемых моделей идентификации и систем автоматического управления исследуется с помощью модуля MatLab/Simulink [10].
1. Постановка задачи Механизм с параллельной кинематикой типа «трипод» состоит из неподвижного основания, подвижной платформы, трех штанг, каждая из которых
состоит из двух полуштанг и исполнительного звена, а также из центрального неподвижного стержня AB - рисунок 1. Легко показать, что число подвижностей
данного механизма W равно двум [2].
Л1, Л2, Л3 - сферические шарниры; B, B1, B2, B3 - шарниры Гука;
1,2,3 - исполнительные звенья Рисунок 1 - Схема трипода с двумя степенями свободы
В качестве исполнительного звена (привода) рассмотрим двигатель постоянного тока в комбинации с винтовой передачей, осуществляющей преобразование вращательного движения вала электродвигателя в поступательное движение штанги.
Введем следующие обозначения: Зде - суммарный момент инерции
механической части электродвигателя, приведенный к валу электродвигателя; Мдв
- момент на валу электродвигателя; Мтр - момент трения; L - эквивалентная
индуктивность якорной цепи; R - эквивалентное сопротивление якорной цепи; Iя -ток якоря; Е - ЭДС электродвигателя; КЕ - коэффициент электродвижущей силы; сод - угловая скорость вращения вала электродвигателя; Ктр - коэффициент вязкого трения; Км - коэффициент момента. Во введенных обозначениях модель двигателя постоянного тока можно записать в виде системы обыкновенных дифференциальных уравнений [11]
dm Д
J —- =МДВ-МТР, дв at
L^*. +1 R = e, dt я
(1)
где E = KEn-, Мтр = КтрП- , Mm = KMI*.
Модель (1) реализована в системе MatLab/Simulink в виде модели DC Machine блока SimPowerSystems [10].
Введем также следующие обозначения: T - момент трения в резьбе; Q - осевая нагрузка на винтовую передачу; dZ - средний диаметр резьбы; а - угол подъема винтовой линии; / - приведенный угол трения в передаче. В этих обозначениях модель винтовой передачи имеет вид [12]
T = Q^tg(а + 3) > Траб.
Кинематика трипода описывается системой алгебраических уравнений [2]
l = h(V ,v) = VZWj- -BA)2 , i, j = 1,2,3
(2)
(3)
где А ^ , Аа2 , АА3 - координаты X,Y, 2 шарнира А{ в инерциальной системе координат AXYZ соответственно; БА1, В А2, ВА3 - те же координаты шарнира В1; I= I((, (2) -длины подвижных штанг; (, (р2 - углы Эйлера, определяющие положение платформы относительно основания:
Bt(Vl,V2)
B fAv^
= T (V1,V2)
B
V BB у
Здесь ВБ1, ВБ2, ВБ3 - координаты шарнира В1 в системе координат Вхух, Т((, (2) матрица перехода от системы координат Вху2 к системе координат АХУ2 .
Динамические уравнения трипода имеют вид системы обыкновенных дифференциальных уравнений [2]
(4)
Jm Р1 (t) = F(t)cosrh2(t)BA(t)-F3(t)cosуъл(t)BA,
<
Jm V2V) = -Fi(t)cos/u(t)BA3, где Jпл - момент инерции платформы,
™ , ™ ^ BAj (p(tр))
cos 7i,j(t) = cos 7i,j (P (t), P2(t)) = 7, -7—.
' ' li (P(t),^2(t))
Здесь силы F1 (t), F3 (t) - усилия, развиваемые приводами штанг (активными полагаются первая и третья штанги).
Необходимо найти такие законы изменения усилий F1(t), F3(t), которые на решениях системы дифференциальных уравнений (4) обеспечивают заданные законы изменения фазовых переменных системы р1 (t), (р2 (t) .
Поскольку непосредственное измерение указанных фазовых переменных технически сложно, будем пересчитывать по формуле (3) требуемые законы движения платформы (p1 (t), р2 (t) в соответствующие законы изменения длин штанг l1 (t) = l2 (t), l3 (t). Положим, что эти длины измеряются без статических и динамических погрешностей.
Таким образом, ставится следующая задача: построить систему автоматического управления, формирующую законы изменения управляющих напряжений для приводов штанг, которые обеспечивают заданные законы изменения их длин l1 (t) = l2 (t), l3 (t).
Положим, что материал изделия - сталь и введем следующие обозначения (рисунок 1):
r0 - радиус основания, r0 = 0,3 м;
Г - радиус платформы, r1 = 0,28 м;
h - высота механизма, h = 0,75 м;
50 - толщина основания, 50 = 0,05 м;
51 - толщина платформы, 51 = 0,0475 м;
d - внутренний радиус штанги, d = 0,03 м;
D - внешний радиус штанги, D = 0,05 м.
Отметим, что указанные параметры трипода соответствуют параметрам промышленного робота FANUC F-200ÍB, построенного на основе гексапода [13].
2. Идентификация 2.1. Классическая идентификация. Положим, что длины двух штанг в рассматриваемой электромеханической системе фиксированы и управляемой является только одна штанга. Выполним непараметрическую идентификацию этой системы с использованием модуля MatLab/System Identification Toolbox. Результатом идентификации являются значения коэффициентов TZ,TP1,TP2,TP3,Td передаточной
функции этой системы
G(s) =--e -v. (5)
W s(1 + Tpis)(1 + Tp 2 s)(1 + TP3s) w
Переходная характеристика Iй (t) динамической системы (5) представлена на рисунке 2, который показывает практическое совпадение этой характеристики с характеристикой lC (t) динамической системы, описываемой уравнениями (1) - (4). Здесь и далее l* (t) - ступенчатое единичное входное воздействие.
-I-*
Рисунок 2 - Результаты классической идентификации: ступенчатое входное
воздействие /*(/)
Хорошо известно, что переходная характеристика полностью описывает только линейную стационарную систему. Поскольку рассматриваемая система не является ни линейной, ни стационарной, при входных воздействиях, отличных от единичного ступенчатого, модель может быть неадекватной. Данное утверждение иллюстрирует
рисунок 3, на котором представлена реакция динамической системы (5) на гармоническое входное воздействии ) с единичной амплитудой и частой, равной 1 с-1.
2.2. Нейросетевая идентификация рассматриваемой электромеханической системы выполнена с помощью МЫАКХ модели. Для построения обучающей выборки использовался ступенчатый единичный входной сигнал I* (^). Некоторые
результаты идентификации представлены на рисунках 4 - 6. Здесь и далее приняты
следующие обозначения: п, т - числа задержек входного и выходного сигналов соответственно; М - число скрытых слоев сети; Ni - числа нейронов в i -ом скрытом
слое; I = 1,2,...
ж
Iе (!)
¿(0
•
Рисунок 3 - Результаты классической идентификации: гармоническое входное
воздействие 1*а([)
'Ш М
ю в Б 4
i О
-г
1
^г \
Vй т
\ /С(1 х
■ц- ,4т \
+ ; i ; ;
t,c
Рисунок 4 - Результаты МЫАКХ-идентификации: ступенчатое входное воздействие /*(/);
п = 1, т = 2, М = 1, N1 = 6
и 12
И
в $
4
2 0
-Е
- г—— -1- + .................,+, - +
......
........ л \ ■
+ .....
........ + ¿г ;
......■+— !
+
+ + * .ХХ+1 щ 1
1 1 •
70 30 ¿0 г,С
Рисунок 5 - Результаты МЫАКХ-идентификации: ступенчатое входное воздействие /*(:);
п = 3, т = 2, М = 1, N1 = 6
> Г
л?
/
у #0
10 Ж X Ф 50 Ш
Рисунок 6 - Результаты МЫАКХ-идентификации: ступенчатое входное воздействие /*(;);
п = 1, т = 2, М = 1, N = 12, N2 = 6
Рисунки 4 - 6 показывают, что NNARX-модель с параметрами п = 1, т = 2, М = 1, N1 = 6 обеспечивает хорошую аппроксимацию переходной характеристики
системы.
Выходной сигнал NNARX-модели при гармоническом входном воздействии /*а (;) иллюстрирует рисунок 7. Сравнение рисунков 3, 7 показывает, что по
сравнению с NNARX-идентификацией, классическая непараметрическая идентификация обеспечивает более высокую точность.
л 1СЩ
Ж V
Рисунок 7 - Результаты МЫАКХ-идентификации: гармоническое входное воздействие /*а(;); п = 1, т = 2, М = 1, N = 6
3. Неадаптивные системы управления
3.1. Система управления на основе ПИД-регуляторов. Структура двухконтурной системы автоматического управления рассматриваемой электромеханической системы приведена на рисунке 8.
ЭД - электродвигатель; ВП - винтовая передача Рисунок 8 - Структурная схема неадаптивной САУ на основе ПИД-регуляторов
Simulink-модель рассматриваемой САУ представлена на рисунке 9. Модель состоит из моделей электродвигателя, винтовой передачи, модели кинематики и динамики трипода, а также из моделей регуляторов. В качестве моделей ПИД-регуляторов использован специализированный блок PID controller, входящего в состав систеы MatLab/Simulink [10]. Поскольку система управления предполагается цифровой, моделируется насыщение цифроаналоговых преобразователей этой системы.
Полагается, что платформа трипода обладает центрально-осевой симметрией, тензор инерции которой определятся матрицей
f 1,8 0 0 ^
0 1,8 0 0 0 3,6
(6)
2
где компоненты имеют размерность кг м .
Переходные процессы в САУ, полученные с помощью рассматриваемой Simulink-модели, представлены на рисунке 10.
Замечание 1. Здесь и далее в качестве возмущающего воздействия рассматривается момент Tвн = 5 Нм, приложенный в момент времени I = 0,5 с к валу электродвигателя.
Замечание 2. Здесь и далее на аналогичных рисунках приведены также графики соответствующих входных ступенчатых воздействий /* ), /3* ) .
Рисунок 10 показывает, что переходные процессы имеют большое перерегулирование. Кроме того, имеют место высокочастотные колебания длин штанг вблизи установившегося значения. Указанные эффекты объясняются наличием
перекрестных связей в системе. Отметим трудность подбора коэффициентов ПИД-регуляторов рассматриваемой САУ.
Поведение САУ в ситуации, когда нарушена центрально-осевая симметрия платформы, иллюстрирует рисунок 11. Из рисунка видно, что данное обстоятельство приводит к возрастанию перерегулирования, переходные процессы становятся не регулярными.
Рисунок 9 - Simulink-модель неадаптивной САУ на основе ПИД-регуляторов
о?: 016 сге
оге
071 0 7
щ
.............................
1
1!
Рисунок 10 - Переходные процессы в неадаптивной САУ на основе ПИД-регуляторов: в момент времени I = 0,5 с приложено внешнее возмущающее воздействие
Замечание 3. Здесь и далее, когда речь идет о нарушении центрально-осевой симметрии платформы, имеется в виду, что тензор инерции платформы гексапода задается не матрицей (6), а аналогичной несимметричной матрицей
(1,8 0 0,8^ 0 1,8 0 0,6 0 3,6
Рисунок 11 - Переходные процессы в неадаптивной САУ на основе ПИД-регуляторов: имеет место нарушение центрально-осевой симметрии платформы
3.2. Нейросетевая система управления. Структура нейросетевой системы управления одной из штанг трипода на основе его инверсной модели приведена на
рисунке 12. Отметим, что здесь, в отличие от рассмотренной выше системы, САУ является одноконтурной.
Г(§ —0
НС1 НС2 ИЗ ОУ
Ж
l(t)
НС1, НС2 - нейронные сети ПИД-регулятора и инверсной модели соответственно;
ИЗ - исполнительное звено; ОУ - объект управления Рисунок 12 - Структура нейросетевой неадаптивной системы управления
Инверсная модель представляет собой обратную модель объекта управления совместно с исполнительным звеном [9]. Эта модель получена при ступенчатых входных воздействиях различной амплитуды и имеет вид, представленный на рисунке 13. Соответствующая нейросетевая инверсная модель, полученная с помощью модуля Nural Nеtwоrk Tооlbоx, представлена на том же рисунке. Рисунок демонстрирует высокую точность нейросетевой аппроксимации инверсной модели - ошибка аппроксимации не превышает 0,5%..
Рисунок 13 -Инверсная характеристика объекта управления и исполнительного звена yinv (t) (сплошная линия) и ее нейросетевая аппроксимация ynn (t) (крестики):
n = 1, m = 2, M = 2, N = 13, N2 = 4
Simulink-модель нейросетевой неадаптивной системы управления приведена на рисунке 14. Здесь, как и в рассмотренной выше Simulink-модели классической САУ, ПИД-регулятор реализован с помощью блока РЮ ^п^оШг. Отметим, что в данном случае настройка регулятора выполняется значительно проще, поскольку все нелинейности
исполнительного звена и объекта управления компенсируются инверсной моделью. Переходные процессов в системе управления иллюстрирует рисунок 15.
Рисунок 14 ^тиНпк-модель нейросетевой системы управления
Рисунок 15 - Переходные процессы в нейросетевой неадаптивной системе управления
Рисунок 15 показывает, что нейросетевая неадаптивная система управления имеет существенно более низкое быстродействие по сравнению с системой на основе ПИД-регуляторов. С другой стороны, видно, что нейросетевая САУ обеспечивает существенно меньшую величину перерегулирования, а также отсутствие флюктуаций вблизи установившихся значений. Подчеркнем, что такая картина имеет место лишь при входных сигналах, "близких" к тестовым.
Эффективность нейросетевой неадаптивной САУ в ситуации, когда имеет место нарушение осевой симметрии платформы (7), иллюстрирует рисунок 16.
Рисунок 16 - Переходные процессы в нейросетевой неадаптивной системе управления: имеет место нарушение центрально-осевой симметрии платформы
4. Адаптивные системы управления
4.1. Адаптивная система управления на базе классических регуляторов. Обозначим Ki = ( kiP, kil, ki D ) вектор коэффициентов ПИД-регулятора,
управляющего длиной штанги ^; i = 1,3 , где ki P, kiI, kiD - коэффициенты при
пропорциональной, интегральной и дифференциальной составляющих регулятора
соответственно.
Будем использовать в качестве ошибки управления величину e(t) - отклонение выхода объекта управления ^ ^) от выхода эталонной модели ¡М ^) :
е ^) = ^) - ¡М ^); г = 1,3.
Таким образом, задача адаптивного управления длиной штанги ¡i состоит в подборе в процессе функционирования системы управления таких значений коэффициентов К{ = К{ ^), которые обеспечивают минимальное значение ошибки ei ^) .
Структура адаптивной системы управления представлена на ЭД -электродвигатель; ВП - винтовая передача
Рис 17. Система имеет три контура - адаптивный контур, а также контуры перемещения и скорости, идентичные соответствующим контурам неадаптивной САУ (рисунок 8).
ЭД - электродвигатель; ВП - винтовая передача Рисунок 17 - Структура классической адаптивной системы управления
Адаптивный контур построен на основе линеаризованной математической модели трипода (эталонной модели) и ПИД-регулятора (рисунок 18).
Рисунок 18. Simulink-модель адаптивного контура
Выход адаптивного контура является эталонным, так что функционирование контуров перемещения и скорости имеет целью компенсировать ошибку отклонения
выхода САУ от выхода эталонной модели. Отметим, что такая САУ позволяет отказаться от идентификации объекта управления в режиме реального времени.
Smu¡ink-модель классической адаптивной системы управления на основе эталонной модели представлена на рисунке 19. Здесь блоки BaseL3, BaseL1 соответствуют эталонной модели трипода.
ri-ri-
T3
Step3
-► " eetLeg3Leng L3UO) L3(t) -+£)-► PIDCsJTPID(sJT —».
12 T 13 PID Controllers PID Controllers
03 0.75
Vbltage
Aigel
T
w
à/
T a Fa eps
it/
Engine
TorqueToAcoeleration
Stepl
^ j-GetLeglLeng LI U1J) L1() PID(sT —PID@T —►
-- T L1 PID Controlled PID Controlled
Constants 0.75
Voltage
Aigel I
T
w
Englnel
T a
Fa eps
rqueToAccelerationl
* AL3 Position
> AL3 V^locrtjr SL3 Position
> AL3 Acoeleration
1+ ALI Position S LI Position
> ALI Vfelooity
> ALI Acoeleration SL3 Force
AL3 Force
> ALI Force SL1 Force
SL3
SL1
t
Tripod
0.75 Const.
n
H >T3 CD
a
s j=i
CD l-k vO
Рисунок 19 - Simulink-модель классической адаптивной системы управления
Переходные процессы в адаптивной системе управления, полученные с помощью данной Smulink-модели, представлены на рисунке 20. Сравнение этих переходных процессов с соответствующими процессами в неадаптивной САУ (рисунок 10) показывает, что в адаптивной САУ отсутствуют высокочастотные колебания. При том же времени переходного процесса, величина перерегулирования в адаптивной САУ значительно меньше (около 1%).
/
V
I
\
\
05 1 V! ?,<:
Рисунок 20 - Переходные процессы в классической адаптивной САУ: в момент времени I = 0,5 с приложено внешнее возмущающее воздействие
Переходные процессы в классической адаптивной САУ в ситуации, когда платформа трипода не обладает центральной симметрией, иллюстрирует рисунок 21.
l,m
цп
5 71
07
f 1 ■ " f"
щК
m
i
----------- ...... —i......... ........i........
f,C
Рисунок 21 - Переходные процессы в адаптивной САУ: имеет место нарушение центрально-осевой симметрии платформы
Сравнение рисунка 21 с рисунком 11, на котором представлены аналогичные переходные процессы в неадаптивной САУ, показывает, что адаптивная САУ менее чувствительна к изменениям параметров объекта управления.
4.2. Нейросетевая адаптивная система управления. Структурная схема нейросетевой адаптивной САУ одной из штанг трипода представлена на рисунке 22.
Здесь регулятор построен на основе нейросетевой инверсной модели объекта управления, а в качестве эталонной модели используется еще одна многослойная нейронная сеть. Задачи обучения обеих нейронных сетей по своей сути представляют собой задачи нейросетевой идентификации соответствующих динамических систем.
НС1, НС2 - нейронные сети эталонной и инверсной моделей соответственно;
ИЗ - исполнительное звено; ОУ - объект управления Рисунок 22 - Структура нейросетевой адаптивной САУ
Simulink-модель указанной системы управления представлена на рисунке 23. Модель аналогична модели неадаптивной нейросетевой САУ (рисунок 12) с тем отличием, что эталонной модели трипода в данном случае соответствует нейронная сеть (блоки ModelL3, ModelLl). Нейронные сети представляют собой сети типа NARX со следующими параметрами: n = 1; m = 2; M = 2; N1 = 12; N2 = 6.
Переходные процессы в нейросетевой САУ, полученные с помощью рассматриваемой Smu¡ink-модели, представлены на рисунке 24. Рисунок показывает, что переходные процессы в нейросетевой адаптивной САУ не отличаются от таких же процессов в нейросетевой неадаптивной системе. Приложение внешней нагрузки приводит к появлению статической ошибки. Таким образом, нейросетевая адаптивная САУ не является робастной (в отличие от аналогичной неадаптивной системы).
Переходные процессы в нейросетевой адаптивной системе управления в ситуации, когда нарушена центрально-осевая симметрия платформы, приведены на рисунке 25. Рисунок показывает, что, по сравнению с аналогичной неадаптивной системой, статическая ошибка в этом случае не возникает.
Рисунок 23 - Simulink-модель нейросетевой адаптивной системы управления
п г — : : : т : : ■
«Л "1..... ■'"Г"' :
0 71 ..... ■4---!........{- : : : :.......„|........|„ : ..... -1..... ..... : ; : .....*....... : ; ........
У : !
л 1 : V : :: ..л. : ............. Щ: г ; : : : ; 1 : [...........| : : ! ! .....!........ .....!......... :
; ; ..... г ■1..... : 9
[
** ч ■ • \ | 1 ; 1 4-4.......
|
"П....;........¿> 1 •
0« : ; I : 2 ::::::::
Рисунок 24 - Переходные процессы в нейросетевой адаптивной САУ: в момент времени I = 0,5 с приложено внешнее возмущающее воздействие
!
/щ { : :
А \ .........:........7........!........ : 1 1 ........|......... |
\ \ А..1........ ! .........1........ ! 1 ! .........I........:........|......... | ........^........
\ ! .....Д.;......... ; .....__._i_.__..... ; ! ! ! ♦ \--------- 1 \ \ .........\........*........;......... ! } 1 ;
:
.........\........ 1 " '1""
г : : ! : !
? ; \ ! :
х -
. < т,с
Рисунок 25 - Переходные процессы в нейросетевой адаптивной САУ: имеет место нарушение центрально-осевой симметрии платформы
Заключение
Результаты работы показывают, что при построении системы автоматического управления многосекционным роботом-манипулятором типа хобот на исполнительном уровне, перспективным является использование нейросетевых адаптивных систем
управления.
В развитии работы предполагается рассмотреть возможность повышения качества управления за счет использования нейро-нечетких САУ, исследовать эффективность САУ при учете стохастических ошибок в канале измерения, а также при учете стохастических внешних воздействиях на платформу механизма.
Авторы благодарят В.М. Буянкина и Ю.Т. Каганова за помощью в выполнении
работы.
Работа выполнена в рамках аналитической ведомственной целевой программы
«Развитие потенциала высшей школы (2009 - 2010 годы)», проект 2.1.2/1509.
Список литературы
1. Глазунов В.А., Колискор А.Ш., Крайнев А.Ф. Пространственные механизмы параллельной структуры.- М.: Наука, 1991.- 94 с.
2. Каганов Ю.Т., Карпенко А.П. Моделирование и оптимизация некоторых параллельных механизмов // Информационные технологии, Приложение, 2010, №5, с. 1-32.
3. Интеллектуальные системы автоматического управления / Под ред. И.М. Макарова, В.М. Лохина.- М.: Физматлит, 2001.- 576 с.
4. Гроп Д. Методы идентификации систем.- М.: Изд-во «Мир», 1979.- 294 с.
5. Александров А. Г. Оптимальные и адаптивные системы. М.: Высшая школа 2002.263 с.
6. Никульчев Е.В. Идентификация систем в Matlab 6 [Электронный ресурс]. (http://www.tspu.tula.ru/ivt/old_site/lcopy/Matlab_RU/systemidentific/book1/index.asp).
7. Пупков К.А., Егупов Н.Д. Методы робастного, нейронечеткого и адаптивного управления.- М.: Изд-во: МГТУ им. Н.Э. Баумана, 2002.- 254 с.
8. Маслобоев Ю.П. Введение в Neural Network Toolbox [Электронный ресурс]. (http://www.nsu.ru/matlab/MatLab_RU/ neuralnetwork/book1/index. asp).
9. Терехов В.А., Ефимов Д.В., Тюкин И.Ю. Нейросетевые системы управления.- М.: Высшая школа, 2002.- 183 с.
10. Черных И.В. Simulink: Инструмент моделирования динамических систем [Электронный ресурс]. (http://matlab.exponenta.ru/simulink/book1/index.php).
11. Онищенко Г.Б. Электрический привод: Учебник для вузов - М.: РАСХН. 2003.-320.
12. Решетов Д.Н. Детали машин: Учебник для вузов.- М.: Машиностроение, 1989.- 496.
13. Merlet J.P.. Parallel Robots. Solid mechanics and its applications.- Kluwer Academic Publishers, V. 74, 2000.
