CC BY
60
Раздел I. Математические методы синтеза систем
УДК 62.50
С.В. Василенко
СИНТЕЗ НЕЛИНЕЙНЫХ ДИСКРЕТНЫХ РЕГУЛЯТОРОВ
В статье показан пример синтеза нелинейных дискретных регуляторов на основе метода подынтервалов и квазимодального управления.
Нелинейный; система; управление.
S.V. Vasilenko
SYNTHESIS OF NONLINEAR DISCRETE CONTROLERS
In this article an example of synthesis of nonlinear discrete controllers is given. This example is based on method of subintervals and quasimodal control.
Nonlinear; system; control.
Введение. Энергетика обладает рядом специфических особенностей, которые отличают ее от других отраслей производства. Процессы производства, , -текают одновременно и являются непрерывными. Практическое совпадение времени производства и потребления энергии обуславливает органическую зависимость между режимами работы энергетических предприятий и потреблением энергии на промышленных предприятиях, транспорте, в сельском хозяйстве и в бытовых целях.
Большая часть электрической и значительная часть тепловой энергии вырабатывается на электростанциях, объединенных в энергетические системы, в которые помимо электростанций входят линии электропередачи, подстанции и
,
процесса производства и распределения энергии.
Качество системы управления напрямую зависит от точности модели управляемого объекта или процесса, однако, большинство моделей используемых в настоящее время для синтеза законов управления энергоблоками являются линейными и поэтому содержат ряд качественных недостатков. Большинство регуляторов сейчас реализуется на основе цифровой электроники, поэтому полученные модели должны быть дискретными. Данная статья посвящена процедуре построения дискретной модели нелинейного объекта и последующему синтезу нелинейного дискретного регулятора.
.
моделям является максимальная точность при как можно большем шаге дискретизации, другим важным моментом является простота расчетов. Данным требованиям удовлетворяет, например, метод подынтервалов [1], который заключается в разбиении большого шага дискретизации на некоторое число малых интервалов (рис.1).
Рис. 1. Разбиение шага дискретизации на подынтервалы
Для построения дискретной модели данным методом исходные непрерывные уравнения объекта представляются в виде
х = Ах + Е(х)х + Ь(х)и .
(1)
<-1
где матрица А обязательно должна иметь обратную А Идея метода подынтервалов состоит в следующем:
• каждый шаг дискретизации разбивается на т равных подынтервалов длительностью Т (рис.1);
• на протяжении подынтервала Т уравнение объекта интегрируется по формуле Коши, причем матрицы ¥(х) и Ь(х) считаются постоянными;
• управление постоянно на всем шаге дискретизации Т ;
• дискретная модель ищется в виде +1 = А(Х^ )Х^ + Ь(Х^ )ик.
Как показано в [1], на основании непрерывной модели (1) можно получить дискретную модель вида:
хкТ+ц,т = Ак,№хкТ + Ьк,№икТ № = 1 т Т = (2)
где
Ак,№ = Л(Е + 0Е(хкТ+(№-1)т))Ак,(№-1);
(3)
Ьк,№ = Л((Е + 0Е(хкТ+(№-1)т))Ьк,(№-1) + 0Ь(хкТ+(№-1)т)) , (4)
(5)
(6)
Ч ,0
к ,о
Л = еАт, 0 = (Е-Л-1)А-1.
Полагая в выражениях (2) (4) /1 = т, получим, с учетом того, что
Т = тТ, искомую дискретную модель вида
Хк+1 = АкХк + Ькик , (7)
Раздел I. Математические методы синтеза систем
Ук = скХк, к = 0,1А-, (8)
где матрица Ак = Ак т (Хк) и вектор Ьк = Ьк т (Хк) вычисляются рекур-рентно с помощью выражений (3) - (5) при Ц = 1, т и определяются значением Хк = ХкТ.
Как было сказано выше, в случае применения метода подынтервалов управление постоянно на всем шаге дискретизации Т . Поэтому закон управле-
(7), (8).
Преобразование к каноническим формам. При определе нных условиях,
(1) , А Ь -
падают по форме с матрицей и вектором канонической управляемой формы системы с постоянными параметрами. Однако имеются и существенные отличия. -,
преобразование переменных состояния, так как элементы матрицы второго преобразования определяются по результатам первого преобразования. Причем преобразование к той или иной форме осуществляется по сути одними и теми же матрицами, только в различной последовательности. Во-вторых, что наиболее
,
функциями переменных состояния и часто зависят от производных по времени переменных состояния системы. Это приводит к зависимости от производных найденного таким путем управления, и, следовательно, правых частей уравнений синтезируемой системы. Однако, условия существования и устойчивости решений таких дифференциальных систем до сих пор неизвестны, что сужает область применения подобного подхода. В [2] было показано, что аналогичные преобразования можно провести и в случае нелинейных дискретных систем.
Приведение системы уравнений (7), (8) к каноническим формам можно осуществить с помощью преобразования Крылова -Луенбергера по входу или по выходу. Предположим, существуют матрицы Бк = £(Хк) ограниченные и неособенные при всех к = 0,1,2,.... Положим
Хк = £к-1~к , Хк+1 = £кХк+1, (9)
Хк - .
В новых переменных состояния система (7), (8) описывается уравнениями
вида:
~к+1 = Ак~к + Ькик , (10)
Ук = ХТ ~к , (11)
где
Х —1 Х —1 ~Т Т
Ак = Бк АкБк—1, Ьк = Бк Ьк , ~к = сТБк—1.(12)
Здесь Б—1 - матрица обратная к Бк, т.е. такая, что £—^Бк = Е.
Как показано в [2], система уравнений (10), (11) имеет каноническую форму Крылова по входу, если
'0 0 . 0 1 1 т Г ~т 1 с1к
II 1 0 . 0 ^... 1 II і ^ 0 > ~к - ~т с2к
0 0. 1 -Х(п-1)к _ 0 ~Т _с(п-1)к _
(13)
~т
гДе %1к = Х(хк,хк-1,хк-2,хк-з), с1к - некоторые нелинейные функции.
Система уравнений (10), (11) имеет каноническую форму Крылова по выходу, если
Ак -
0 1 0 і к і
0 0 0 , -[1 0 . 0 к ( II Ь2к
0 0 1 Ь(п - 2)к
«■к її; і 1 -а1к ■■ - п 1 к Ь(п-1)к
(14)
где а,к = а, (х^ Хк+и Хк, хк-1), Ьк - некоторые нелинейные функции.
В общем случае матрица перехода к канонической форме по входу и функции Х,к определяются выражениями [5]:
п - 2
Рк = Ьк АкЬк-1 АкАк-1Ьк-2 ••• ( ПАк-])Ьк-п+1
_ 3 =0
п -1
х=-Р(,+1)к( ПАк-3)Ьк-п 1, (15)
3 =0
матрица перехода к канонической форме по выходу и функции ак определяются выражениями:
п-1
А (п
т- -
(сТ )Т (сГ+1Ак+1)Т (4+2 Ак+2 Ак+1)Т ••• (сТ+ и-1 П Ак+(и - І ))Т
І-1
, а сТ+(п-1)( П Ак-(п - 3))Тк-1.
3=1
,
к КНФ или КУФ необходимо выполнить двукратное последовательное преобразование переменных состояния. Для перехода к КНФ сначала производится преобразование к канонической форме по выходу (14), (15), а затем - по входу (14), (16). Для перехода к КУФ эти преобразования нужно выполнить в обратном по.
нелинейной дискретной системы к КНФ или КУФ, что расширяет возможности анализа и синтеза подобных систем. Например, делает возможным синтез квази-модального управления [2].
(16)
Т
Синтез квазимодального управления. Уравнения нел инейной дискретной системы приведенной к КУФ с помощью преобразований Крылова - Луен-бергера имеют вид:
O І O O
xk+І - O 0 І xk + O
1 1 kO -a1k 1 k 1 і
uk .
(І?)
Для полученной системы (17) синтезируем квазимодальное управление. Закон управления будет иметь вид:
uk = — Kkxk , Kk k + Kмoд , (18)
где Kk - вектор коэффициентов управления, k - вектор компенсирующий
нелинейности системы, Kw0(,)- вектор коэффициентов модального управления.
С учетом (18) получим:
k = [- ^ - - ^]. (19)
Вектор модального управления имеет вид
= IX) Х1 Х2 ], (20)
где X - коэффициенты желаемого характеристического полинома системы. Тогда уравнения замкнутой системы будут иметь вид:
O і O O
xk+! - O 0 і xk + O
-Aq -АІ — A2_ і
gk.
(2 і)
Подставив (19), (20) в (18) и совершив обратное преобразование переменных, получим закон для модального управления исходной нелинейной дискретной системой
uk
KkPk—-\Tk—-!xk
(22)
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Гайдук А.Р., Василенко С.Я^горитмическая нелинейная модель энергоблока// Известия ТРТУ. 2004. - №7.
2. Василенко САСинтез дискретного квазимодального управления для нелинейных систем// Известия ТРТУ. 2008. - №7.
Василенко Сергей Валерьевич
Технологический институт федерального государственного образовательного учреждения высшего профессионального образования «Южный федеральный университет» в г. Таганроге Email: Salouma1@mail.ru
347928, Таганрог, ГСП 17А, Некрасовский, 44. Тел: 88634-371-689
ІІ
Vasilenko Sergey Valerevich
Taganrog Institute of Technology - Federal State-Owned Educational Establishment of Higher Vocational Education “Southern Federal University”
Email: Salouma1@mail.ru
44, Nekrasovsky, Taganrog, 347928. Phone: 88634-371-689 УДК 518.5.001.57
T.B. Камышникова, A.A. Афонин, B.B. Дурягина
ВЫВОД ДВУМЕРНОЙ МОДЕЛИ РАСПРОСТРАНЕНИЯ ЗАГРЯЗНЯЮЩИХ ПРИМЕСЕЙ В ВОДОХРАНИЛИЩЕ
Изучение гидрохимических характеристик вод водохранилищ является весьма актуальным, так как их изменчивость, как правило, является отражением возможных отклонений в функционировании водоема. И здесь математическое моделирование может оказать неоценимую помощь.
Модель «мелкой воды»; адвективно-диффузионное уравнение.
T.V. Kamyshnikova, A.A. Afonin, V.V. Duryagina
DISTRIBUTION OF CONSUMERS BETWEEN SOURCES OF ELECTRIC
ENERGY
Water basin hydro chemical characteristics’ studying is actual problem because of variability as possible deviations ’ reflection in water basin functioning. The best method of this problem solving is mathematical modeling.
Shallow water model; the advection-diffusion equation.
В последние десятилетия научные организации все больше внимания стали уделять исследованиям современных и ожидаемых антропогенных изменений гидрологических и гидрохимических условий водоемов, включая химические загрязнения воды. Загрязнение нарушает гидрохимический режим
,
воде. Основной источник загрязнения водных бассейнов - сброс недостаточно очищенных промышленных стоков, а также возвратных вод с .
относятся нефтепродукты, фенольные соединения, пестициды (хлор и фосфорорганические ядохимикаты), детергенты (или СПАВ-синтетические поверхностно-активные вещества).
Процессы самоочищения в водотоках бассейна идут с недостаточной
,
органические токсиканты не успевают минерализоваться в пределах водотока и около половины их объема поступает в водоем. Прибрежные воды водоема, в которые непосредственно поступают сточные и речные воды, содержащие
,
отличаются более высокими концентрациями загрязняющих веществ по сравнению с водами открытой части водохранилища. Эти концентрации меняются от года к году, как в зависимости от величин сбросов сточных вод, так и под влиянием течений, весьма изменчивых по времени [1].
Развитие методов расчета поля примесей, позволяющих выполнить анализ
