УДК 681.51.011
Г.В. Масютина, В.Ф. Лубенцов
СТРУКТУРНО-ПА^МЕТРИЧЕСКИЙ СИНТЕЗ АДАПТИВНОЙ СИСТЕМЫ УПРАВЛЕНИЯ НА ОСНОВЕ НЕЧЕТКОЙ ЛОГИКИ
В работе предложен подход к структурно-параметрическому синтезу адаптивной системы управления и стабилизации сложного объекта на основе методов нечеткой логики и переменной структуры. Полученные результаты исследований системы подтверждают повышение эффективности управления, что в условиях возрастающей сложности объектов и процессов является актуальным.
Система; нечеткая логика; адаптивное управление; корректирующее воздействие; -.
G.V. Masyutina, V.F. Lubenzov
STRUCTURED-PARAMETRIC SYNTHESIS OF THE ADAPTIVE CONTROL SYSTEM BASED ON THE FUZZY LOGIC
The approach to structured-parametric synthesis of the adaptive control system of the complex object based on the methods of the fuzzy logic and variable structure is offered. The received results of research of system confirm increase of a management efficiency, that in conditions of growing complexity of objects and processes is actual.
System; fuzzy logic; adaptive control; correcting influence; PID-controller.
Современные методы построения систем автоматического управления основаны на использовании строгих математических моделей объектов. Однако для подавляющего большинства динамических систем построение точных математических моделей практически невозможно ввиду их плохой формализуемости. К тому же, эти системы могут функционировать в среде, свойства которой изменяются или же вообще не могут быть определены заранее. Структурно-параме^зический синтез
таких систем возможен на основе аппарата нечеткой логики и принципов перемен,
системы, сохраняя общие принципы адаптивного и робастного управления. Согласно этим принципам в основном контуре синтезируемой системы управления использован регулятор с аппроксимирующей непрерывной нелинейной функцией (АНФ-регулятор), который обладает определенной эффективностью по сравнению с обычными релейными и многопозиционными регуляторами [1], а в качестве корректирующего устройства к АНФ-ре^лятору, как показали предварительные исследования [2], целесообразно использовать нечеткий регулятор (фаззи-). - -ванный нелинейный закон управления U(e) = F1 (е) + F2 (е) = M1 signe + M2sign(e-a) (e - рассогласование между заданной g(t) и текущей регулируемой переменной x(t), F i (s) - релейная характеристика элемента без зоны нечувствительности "а" (ЗН), F2 (s) - релейная характеристика элемента с ЗН) функцией вида
где М - величина регулирующего воздействия; X - коэффициент усиления нелинейности в процессе управления.
Использование и-мерной линейной комбинации функций Е(е) с различными значениями М , а и X. даёт возможность аппроксимировать сложные нелинейные зависимости, используемые в качестве функций управления (законов регулирования).
[1 + exp(Л ■(е + i ■ а))] [1 + exp(-Л ■ (е-i - а))]
exp(Л■ (е + i ■ а)) - exp(-Л ■ (е- i -а))
Нечеткие регуляторы с фиксированной базой правил известны как статические нечеткие регуляторы и были разработаны в 1970-е годы как результат внедрения фундаментальных работ Л. Заде в задачи автоматического управления [3,4]. Однако эффективность применения фиксированной базы правил снижается на различных режимах, характерных, например, для систем программного управления и стабилизации, когда в переходном режиме имеют место значительные отклонения выходной переменной объекта управления от требуемого значения и требуется высокая точность стабилизации переменной на заданном уровне в уста. , для синтеза которого предлагается использовать несколько параллельно включен-
( - ).
На рис. 1 приведена структурная схема системы с двумя фаззи-корректорами и логическими блоками, обеспечивающими выбор и переключение баз правил (БП) с разными видами терм.
Рис. 1. Структурная схема адаптивной системы управления с двумя фаззи-корректорами
- . -
ными являются отклонение еф выходного сигнала хот заданного значения g(t) и производная (ёе/Ж) этого отклонения, терм-множество которых является следующим: {«мапое» (М), «норма» (Н) и «большое» (В)}. Для выходного сигнала перво-- 1 : « » ( ), « -» ( ) « » ( ). -корректора 2 примем следующие выражения: «сильно уменьшить (СМ)», «умень-» ( ), « » ( ), « » ( ) « » ( ). Оперируя входными и выходными переменными и учитывая, что чем больше отклонение хф от заданно го значения, тем больше должно быть корректирующее воздействие на коэффициент X передачи (усиления) АНФ-регулятора, составляем базу правил (табл. 1,2).
1
Совокупность правил представления для трех терм
Отклонение регулируемой переменной (е) Скорость изменения регулируемой переменной (^ М^)
М н в
М М м н
н м н в
в н в в
Таблица 2
Совокупность правил представления пятью термами
Отклонение регулируемой переменной (е) Скорость изменения регулируемой переменной (de/dt)
М н в
М СМ м н
н м н в
в н в св
При синтезе системы использованы фаззификация треугольными функциями принадлежности как наиболее простыми и алгоритм логического вывода по Мам-дани [5] как наиболее простой и позволяющий выносить суждение о том или ином параметре в виде утверждений, в отличие от их представления линейными функциями, характерного для нечеткого вывода Сугено. Для перехода от нечетких выводов к точному значению корректируемого параметра АНФ-регулятора используем формулу дефаззификации по методу центра тяжести.
, -
,
САУ будет возникать статическая ошибка. Для сведения ее к нулю и придания системе свойства астатизма требуется введение новых правил в базу знаний нечет. , базе традиционного ПИД-регулятора, является этап фаззификации его настроечных параметров, и, в частности, такого ключевого параметра для обеспечения астатических свойств системы как времени интегрирования Ти регулятора, так как термы «уменьшить Ти» или «существенно увеличить Ти» и т.п. не всегда могут соответствовать оптимальным значениям настроек регулятора, обеспечивающих хотя и устойчивый, но и одновременно не качественный переходный процесс. При этом до сих пор не установлено, какое количество решающих правил необходимо для достижения астатических свойств системы. В большинстве случаев количество правил экспоненциально возрастает с повышением точности [6]. Такой неограниченный рост приводит к существенному возрастанию вычислительной сложности нечеткой модели и, в конечном итоге, к практической нереализуемости.
Исходя из вышесказанного, в данной работе предлагается ввести в алгоритм управления нечеткой системы непрерывный интегратор. Для обеспечения требуемых показателей качества наряду с астатическими свойствами при наличии переходных и установившихся режимов в системе установлена целесообразность изменения коэффициента интегрирования с помощью следующего нелинейного закона:
Г а при £(t) >S; Ки=1/Ти = „
[в при £(t) <5,
где а, Р - коэффициенты настройки интегратора; s(t) - сигнал рассогласования; д - заданное значение сигнала рассогласования 0 = const).
В работе были получены переходные процессы при различных значениях параметра д (рис. 2). Удовлетворительные настройки интегратора а = 0,04 мин-1, Р = 0,1 мин-1 и значение S = 2 м.е. (м.е. - машинные единицы) получены в ходе математического моделирования. Установлено, что при неизменном количестве решающих правил нечеткий регулятор с переменным коэффициентом интегрирования Ки = f(e) позволяет улучшить показатели качества переходного процесса по сравнению с нечетким регулятором с постоянным значением Ки = 0,154 мин-1.
Опишем кратко работу системы. В системе задаются пороговые максимальный верхний Е1 и нижний Е2 уровни отклонения (е) выходной переменной от зада-
ния, а также задается максимально допустимое время устранения отклонения (Т~ (3... 5)ТР, где Тр - длительность переходного процесса).
Рис. 2. Графики переходных процессов в нечеткой СА У при изменении уставки д: 1- д = 0, 2- д = 1 м.е., 3- д = 2 м.е., 4- д = 4 м.е., 5- д = 5 м.е.
С учетом этого в блоке 1 вычисляются величины S1 = max(/ε/-E1) и S2 = max(/E2/-/£/ ), на основании которых оценивается необходимость подключения фаззи-корректора с той или иной базой правил. При этом, чем продолжительнее отклонение от порогового значения, тем большим должно быть корректирующее - . -лизованы в следующем виде:
если (/e/-Ej) > 0 и Т > ЗТр и БП=БП2 и Хк = f1(e, de/dt), = 1, = 2; если (/E2/-/zi)<0 и Т > ЗТр и БП^П1 и Хк = f2(e, ^e/dt),
то выбрать БП^П2, иначе БП^П1. (1)
Пусть в установившемся режиме на вход переключающего реле ПР поступает корректирующее воздействие Хк с выхода фаззи-корректора 1. При превышении ошибкой регулирования первого порогового значения Ei начинается отсчет времени блоком 2 и вычисление величин S1 и S2. В логическом блоке З проверяются (1). (1) З -
мандный сигнал на управляющий вход ПР, на второй вход которого подключен - 2.
фаззи-корректора 2 на вход АНФ-регулятора, изменяя настроечный параметр Хк в соответствии с новой базой правил. По истечении времени Т, равного ЗТР, вели-S1 S2 .
Модель нечеткого регулятора составлена средствами пакета «Fuzzy Logic Toobox». Нечеткие понятия, взятые выше в кавычки, представлены в модели термами лингвистических переменных, с которыми сопоставлены треугольные функ-.
на вход системы скачкообразного изменяющегося сигнала задания на трех вре-. ,
,
(t) d (t)/dt, -цессы без существенного перерегулирования с приемлемым временем регулирования. Такой характер переходных процессов объясняет несколько большее время переходных процессов по сравнению с системой управления с одной БП. В табл. З,4 приведены оценки показателей качества переходного процесса (максимальное ди-
max, , ,
).
ПИД-ре^лятора были рассчитаны с использованием критерия максимальной степе. -способ, предложенный в [7].
Анализ полученных результатов показывает, что качественные характеристики
переходного процесса в системе с нечетким АНФ-ре^лятором с двумя базами пра-, - -
регулятором. Следует также отметить, что у нечеткого АНФ-ре^лятора корректируется лишь всего один настроечный параметр, тогда как у ПИД-ре^лятора их три.
Таблица 3
Оценки показателей качества переходного процесса в системе с нечетким АНФ-ре^лятором с двумя базами правил для объекта с моделью
W(p) = Ко6 exp(-Tp) /(r2V + ТХp + 1)
Значения параметров модели Показатели качества переходного процесса в системе
объекта управления нечеткий нечеткий
АНФ-регулятор АНФ-регулятор
с двумя БП с одной БП
№ К05, ед. г, Ть ТА струк- А, а,% v, ед т„, А, ед о, % v, ед
вых/ед.в мин мин мин2 тура ед мин мин
X. БП
1 0,1411 4.0 25,50 233,48 2-1-1 0** 0 0 400 0,24 4,8 1.0 200
2 0,2298 4,1 27,96 281,57 2-1-1 0 0 0 300 0,64 12,8 0,44 425
3 0,0766 6.0 41.56 73,10 2-1-1 0 0 0 400 0,25 5.0 1.0 400
4 0,3830 6.0 41.56 73,10 1-2-2 1.0 12,5 1,0* 250 1.55 31,0 0,77 550
5 0,0766 6.0 14.35 625,50 2-1-1 0 0 0 400 0 0 0 400
6 0,3830 6.0 14.35 625,50 1-2-2 2.2 27,5 0,86 600 0,09 1.8 1,0 62,8
7 0.0766 6.0 41.56 625,50 2-1-1 0 0 0 600 0,27 5.4 1.0 325
8 0,3830 6,0 41,56 625,50 2-1-1 1.1 22,0 1,0 250 1,55 31,0 0,77 500
Примечание: * степень затухания \|/ = 1 соответствует апериодическому переходному процессу;
** нули в строке соответствуют монотонному переходному
( );
аналогично для табл.4.
Последовательность цифр 2-1-1, 1-2-2 и т.д. соответствует реализации первой, второй и первой или первой, второй и второй базы правил и т.д.
Совокупность всех правил для нечеткой переменной с функцией принадлежности с тремя термами обозначена базой правил БП1, а совокупность всех правил для нечеткой переменной с функцией принадлежности с пятью термами обозна-2.
4
Оценки показателей качества переходного процесса в системе с нечеткими регуляторами для объекта с моделью Ж(р) = Коб ехр(-тр) 1(Т1р2 + Т1 р + 1)
Значения параметров модели объекта ¡управления Максимальная степень устойчивости Показатели качества
тип регулятора
нечеткий АНФ-регулятор нечеткий ПИД-регулятор
№ К», ед. вых/ед.вх. т, мин Т|, мин Т;-. мин2 J А, ед о,% V, ед Тип, мин А, ед а, % V, ед Тал, МИН
1 0,1411 4,0 25,50 233,48 -0,141 0.24 4,8 1,0* 200 0,48 9,6 1,0 250
2 0.2298 4,1 27,96 281,57 -0,135 0.64 12,8 0.44 425 1,55 31.2 0,88 300
3 0,0766 6,0 41,56 73,10 -0,171 0,25 5,0 1,0 400 0,31 26.2 1,0 350
4 0.3830 6,0 41.56 73,10 -0,171 1,55 31,0 0,77 550 2.0 20.0 0,83 375
5 0,0766 6,0 14,35 625,50 -0,079 0" 0 0 400 0,41 8,2 0,51 700
6 0.3830 6,0 14.35 625.50 -0,079 0.09 1.8 1.0 62.8 1.9 38.0 0,74 574
7 0,0766 6,0 41.56 625,50 -0,092 0.27 5,4 1,0 325 0.4 8,0 1.0 300
8 0.3830 6,0 41.56 625,50 -0,092 1.55 31,0 0,77 500 4,23 84.6 0,21 700
Для обеспечения необходимой точности и астатических свойств системы достаточно скорректировать настройку интегратора, а структуру нечеткого регулятора
оставить без изменений, ограничив базу правил минимальным количеством термов.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Лубенцов В.Ф. Исследование динамики систем с непрерывными аппроксимирующими функциями управления // Наука и технологии: Труды XXV Российской школы по проблемам науки и технологий. 4.2. - М.: РАН, 2005. - С. 469-476.
2. Лубенцов В.Ф., Масютина Г.В. Коррекция регуляторов нелинейных систем с использо-
/ // -ная научная конференция 29.09-02.10.2009 г. Пятигорск. Сб. докладов. - Пятигорск: РИА на КМВ, 2009. - С. 220-223.
3. Заде Л. Понятие лингвистической переменной и его применение к принятию приближенных решений: Пер. с англ. - М.: Мир, 1976. - 165 с.
4. Zadeh L. The role of fuzzy logic in the management of uncertainty in expert systems // Fuzzy Sets a Systems. - 1983. - Vol. 11, № 3. - P. 553-557.
5. Леоненков А.В. Нечеткое моделирование в среде MATLAB и fuzzy TECH. - СПб. : БХВ - Петербург, 2005. - 736 с.
6. Кукса ПЛ. Обеспечение точности в нечетких системах [Электронный ресурс]. URL: http://paul.rutgers.edu/~phuksa/publications/fz - accuracy - iu-04.pdf.
7. . - : //
. - 2007. - 1. - . 78-88.
Масютина Галина Владимировна
Невинномысский технологический институт ГОУ ВПО «Северо-Кавказский государственный технический университет».
E-mail: [email protected].
357108, . , . , 1.
Тел.: 88655471343; 8865570401.
Лубенцов Валерий Федорович
Masyutina Calina Vladimirovna
Nevinnomyssk institute of technology "North Caucasian state technical university".
E-mail: [email protected].
1, Gagarina street, Nevinnomissk, 357108, Russia.
Lubenzov Valery Fedorovich
УДК 681.51
А.А. Колесников
СИНТЕЗ И МОДЕЛИРОВАНИЕ СИСТЕМ УПРАВЛЕНИЯ ОРБИТАЛЬНЫМ ДВИЖЕНИЕМ КОСМИЧЕСКИХ ЛЕТАТЕЛЬНЫХ
АППАРАТОВ
Решена задача управления полетом космического летательного аппарата с малой тягой с использованием системного закона гравитационного взаимодействия. Приведены процедура синтеза и результаты моделирования.
Моделирование; космический аппарат; нелинейная динамика; нелинейные системы; .