CC BY
43
Раздел IV. Математическое моделирование и обработка данных
УДК 681.51.01
А.В. Семенов, А.Р. Гайдук
БИНОМИАЛЬНЫЕ УСЛОВИЯ КОМПЕНСАЦИИ ПОЛИНОМИАЛЬНЫХ
ВОЗДЕЙСТВИЙ
В работе сформулирована и доказана теорема об условиях астатизма произвольного порядка дискретных систем управления по отношению к задающему воздействию. Полученные условия относятся к коэффициентам числителя и знаменателя передаточной функции по этому воздействию. Условия получены методом коэффициентов ошибок и включают биномиальные коэффициенты, определяемые треугольником Паскаля. В работе представлена как аналитическая, так и графическая формы полученных условий астатиз-ма. В графической форме используется широко известный треугольник Паскаля. Приведены примеры условий астатизма различных порядков.
Полученные условия астатизма позволяют конструировать желаемые передаточные функции дискретных систем с высоким (произвольным) порядком астатизма при небольшом перерегулировании, колебательности и желаемой длительности переходных процессов.
Астатизм; условия астатизма; дискретная система; биномиальные коэффициенты; треугольник Паскаля.
А^. Semenov, А.& Gaiduk
BINOMINAL CONDITION OF POLYNOMIAL DISTURBANCES COMPENSATION
In the given paper arbitrary order astaticism conditions of discrete control systems with respect to the reference disturbance is stated and proven. Obtained conditions are the coefficients of numerator and denominator of the transfer function of this disturbance. Conditions are obtained by the error ratio method. They include the binomial coefficients which determined by the Pascal triangle. The paper provides both analytic and graphical presentation form of obtained astaticism conditions. Pascal triangle was applied in the graphic form. Examples astaticism conditions of different orders are described.
Obtained astaticism conditions allow to design desired transfer functions of discrete systems with high (arbitrary) order of astaticism with small overshoot, vibration and the desired duration of the transition.
Aastaticism; astaticism conditions; discrete system; binomial coefficients; Pascal triangle.
Введение. Одним из основных критериев качества процесса регулирования является точность систем автоматического управления (САУ) как непрерывных, так и дискретных. Астатические системы воспроизводят входные и подавляют возмущающие полиномиальные воздействия с ошибкой, равной нулю. Астатическими обычно создаются следящие системы, для которых важным является отработка без ошибок задающих полиномиальных воздействий заданной степени [1]. При проектировании мехатронных следящих систем актуальными являются задачи обеспечения нулевых ошибок по положению, скорости, ускорению и т.д., т.е. раз-
работка управления с астатизмом первого, второго, третьего порядка и более высоких порядков [1, 2, 3]. Особенно актуальными эти задачи являются в настоящее время в связи с широким распространением цифровых средств управления (микроЭВМ, микроконтроллеры, ПЛИС).
При синтезе цифровых регуляторов некоторым методам необходимы, в частности, желаемые дискретные передаточные функции (ПФ) замкнутых систем, обладающих заданными показателями качества [4]. Эти функции удобно формировать на основе стандартных ПФ непрерывных систем с помощью известного ZT -преобразования [5]. Однако в этом случае преобразованная функция сохраняет лишь показатели качества и свойство астатизма первого порядка. Свойства аста-тизма более высоких порядков - теряются [4].
Поэтому в данной работе ставится задача получения условий астатизма высоких порядков для дискретных систем, которые могли бы использоваться для конструирования соответствующих дискретных передаточных функций.
Постановка задачи. ПФ ДС в замкнутом состоянии Wyg(z) по задающему воздействию g(k) имеет следующий вид:
Wyg(Z) S0+S1z+S2z2^-----hSn-1zn~1 + Snzn . (1)
В этом выражении по условиям физической реализуемости т < п [1, 5]. Точность дискетных систем определяется величиной ошибок, с которой система воспроизводит задающие и подавляет возмущающие воздействия. Для определения ошибок ДС, вызванных полиномиальными воздействиями, удобнее всего применять метод коэффициентов ошибок [1, 5]. Вызванная некоторым задающим полиномиальным воздействием g(k) степени rg ошибка системы с ПФ (1), в соответствии с методом коэффициентов ошибок, определяется выражением
5g(rg) = |Coggk + ^ Д gj, (2)
где Д1 gk- i-я разность воздействия g Cjg - i-й коэффициент ошибки по этому воздействию; |*| - обозначение модуля [5].
Чтобы проектируемая система обладала астатизмом vg-ro порядка по отношению к задающему воздействию g(k), необходимо, чтобы ошибка, вызванная полиномиальным воздействием g(k) степени rg = (vg — 1) или меньше, была равна нулю, а ошибка, вызванная полиномиальным воздействием g(k) степени rg = vg не была равна нулю, но была постоянной. Выражение (2) для системы, которая имеет астатизм vg-ro порядка, можно представить в виде
5g = |C0ggte + S^QgДgk| = 0 rg — 1 (3)
Sg = |cV(^vggk| = const Ф 0, rg = Vg. (4)
Для случая, когда степень полиномиального воздействия g(k) равна rg = vg и vg < п, из (3) и (4) следует, что значения коэффициентов ошибок системы с астатизмом Vg-го порядка должны быть равны:
C(v-i)g = 0, при v = 1, Vg ; CVgg Ф 0. (5)
Следовательно, чтобы система обладала астатизмом Vg-го порядка, необходимо и достаточно, чтобы первые Vg — 1 коэффициентов ошибок этой системы были равны нулю, а коэффициент ошибки с индексом Vg не был равен нулю [1, 5].
В ряде работ приводятся алгебраические условия астатизма в виде условий на коэффициенты ПФ ДС в замкнутом состоянии либо только для астатизма первого порядка [5], либо в форме, неудобной для практического применения [1]. Позднее в работах [6, 7] были получены алгебраические условия астатизма произвольного порядка, которые для астатизма первого порядка и второго порядка сов-
падают с приведёнными в [1, 5]. Сформулируем теорему об алгебраических условиях астатизма произвольного порядка vg по задающему воздействию §(к), которым должны удовлетворять коэффициенты ПФ (1).
Теорема 1 (Об условиях астатизма дискретных систем). Дискретная система управления с ПФ (1) имеет астатизм vg-го порядка по отношению к задающему полиномиальному воздействию g(k) степени = vg, если коэффициенты этой ПФ удовлетворяют следующим условиям:
1^-1(1-1 )т = йи-хС-О^ , V = (6)
(7)
где - биномиальные коэффициенты, определяемые формулой = И/к! (1 - к)!.
Доказательство. Передаточная функция по ошибке линейной системы с ПФ (1), вызванной задающим воздействием, определяется по формуле [1,5]:
"^(г)=1 - "^(г). (8)
Признаком астатизма vg-го порядка ДС к задающему воздействию g(k) является возможность представления её ПФ (8) в виде
1)^, (9)
где ^N^(7) - некоторая рациональная дробь, все нули знаменателя которой по модулю строго меньше единицы.
Как известно [5], коэффициенты ошибки С^ ДС определяются путем разложения правой части выражения (9) в ряд по степеням двучлена (г — 1), т.е.
+ С^ —1) + С2g(z — 1)2 + •••+ С^ — 1)* + •••. (10)
Проведя в выражении (10) замену г на Я + 1, т.е. г = Я + 1, получим
+ С^Я + С^Я2 + ••• + С^ + •••. (11)
С другой стороны, из выражений (1) и (8) следует равенство
" (7)= ео+е12+-+ет_1гт~1 + етгт+-+еп_1гп~1+епгп
80+81г+82г2+-+8п-1гп~1+8пгп ' ( )
где еI - коэффициенты, определяемые по формуле:
("(5; — при I < т, при т < I < п.
Полагая в выражении (12) г = Я + 1, раскрывая скобки вида (Я + 1)1 по формуле бинома Ньютона, приведем подобные. В результате будем иметь:
(14)
где
ёк = , 4к = о^ . (15)
Приравнивая (11) и (14), последовательно найдем выражения для коэффициентов ошибок системы с ПФ (1)
п _ ё0 _ _ ё1-б1Со^ ^ _ £2-51С11-52СО1
~Т0' Clg - 80 ' C2g = Т0 ' - (16)
Подчеркнем, что коэффициенты ё^ и связаны с коэффициентами ^ и ПФ (1) формулами (13) и (15). Поэтому для определения коэффициентов ё1 воспользуемся условиями (6) и (7). Перенеся сумму левых частей выражений (6) и (7) вправо, соответственно получим
№ — Н^-Л-ОЪ - о, V = (17)
ТЗД— ТЗД)т * 0. (18)
Преобразуем суммы коэффициентов в левых частях выражений (17) и (18) следующим образом:
Н^-Л-ОФ — лд — ТГ=т+1а-1^ -0, V = (19)
в. -
''-1 Si:
(13)
Е^Д)№ - т) - Е^+Д* 0. (20)
С учетом (13) из выражений (19) и (20) получаем
Еи-1&-1 )е1 = 0, V = )е1 * 0. (21)
Сравнивая суммы в выражении (21) и первом выражении (15), видим, что они равны при к = V -1. Перепишем выражения (21) с учетом обозначений (15) для ёк. В результате получим
ёу_1 = 0, V = 1,^; ё^ * 0. (22)
Подставляя значения и ё^, полученные в (22), в формулы (16), найдем, что = -'- = С = С = 0 * 0. Отсюда, согласно (5), следует, что если
"г
коэффициенты ПФ ДС, удовлетворяют условиям (6), (7) теоремы, то она имеет астатизм vg-го порядка. Доказательство окончено.
Подчеркнем, что число равенств типа (6) равно порядку астатизма vg. Отметим также, что выражения (6) можно использовать и для случая т = п. При этом т,п > vg.
Примеры условий астатизма. Приведем ряд условий, при которых система с ПФ (1) имеет заданный порядок астатизма. Если vg = 1, из выражения (6) следует условие вида = которое совпадает с приведенным в [1,5]. Если vg = 2, то условия (6) соответствуют двум равенствам: = £Г=о Ф;
¿^ = ¿5;, которые совпадают с приведенными в [5]. Если vg = 3, то условия (6) соответствуют трем равенствам: = ЕГ=оФ; = £¿^=2(2)^1 = £7=2(2)^1. Таким образом, используя биномиальные соотношения (6), можно записать алгебраические условия астатизма произвольного порядка дискретных систем с ПФ (1).
Биномиальные условия астатизма на основе треугольника Паскаля. Заметим, что в выражениях (6) и (7) «весовые множители» перед коэффициентами и ПФ "^(г) могут быть достаточно просто найдены по треугольнику Паскаля, предлагаемый вид которого приведен на рис. 1.
По сравнению с классическим видом треугольника Паскаля, принятого в математике, здесь добавлены параметры Уг,т,п систем управления, определенные в соответствии с выражениями (1) и (6), (7). Вдоль «правого ребра» треугольника располагаются значения V от 1 до vg. Вдоль «левого ребра» треугольника располагается индекс I весовых множителей. Фактически весовые множители v-го условия астатизма образуют в треугольнике Паскаля у-ю «левую диагональ», которая параллельна «левому ребру» треугольника.
Каждый весовой множитель (5,-1) находится на пересечении v-й «левой диагонали» треугольника с ¿-й строкой. Например, весовой множитель (з_х) из 3-го равенства (6) находится на пересечении 3-й «левой диагонали» треугольника с 7-й строкой и равен 21; а весовой множитель (|_х) из 4-го равенства (6) находится на пересечении 4-й «левой диагонали» треугольника с 6-й строкой и равен 20 (см. пример 1 на рис. 1).
Рассмотрим пример получения алгебраических условий астатизма 6-го порядка для системы управления с ПФ вида (1) при т = 6, п = 7 по треугольнику Паскаля. В соответствии с заданием, поставим маркеры т = 6, п = 7, vg = 6 и отметим на треугольнике Паскаля границы поиска весовых множителей: «левая диагональ» с номером v = vg = 6, строка с номером I = п = 7. Поскольку vg = 6, то и условия будут состоять из шести равенств вида (6), которые примут вид:
80 + ¿1 + ¿2 + ¿3 + £4 + £5 + ¿6 + £7 = + VI + Л2 + Лз + Л4 + Лв + Лв' 8г + 252 + 3£3 + 4^4 + 5^5 + 686 + 787 = + 2?72 + 3г]3 + 4?74 + 5г]5 + 6г]6: 82 + 3£3 + 6^4 + 1055 + 15£6 + 21£7 = ц2 + 3?7з + 6?74 + + 15^6; £3 + 4^4 + 10^5 + 20£6 + 3587 = + 4г}4 + 10^5 + 20^6; 84 + 5^5 + 15£6 + 35^7 = + 5т75 + 15^6; 85 + 6£6 + 21^7 = + 6?76. Заметим, что коэффициенты пятого равенства полученных условий взяты из пятой «левой диагонали» треугольника (см. пример 2 на рис. 1).
Рис. 1. Диаграмма получения алгебраических условий астатизма дискретных
систем управления
Подчеркнём, что для выполнения требований к порядку астатизма коэффициенты ПФ Wyg (7) (1) должны удовлетворять одновременно всем полученным равенствам.
Заключение. Полученные в работе биномиальные условия компенсации полиномиальных воздействий позволяют: определить порядок астатизма ДС с заданной ПФ; сконструировать ПФ ДС, которая должна обладать заданным порядком астатизма. Последнее можно выполнить, например, путем изменения коэффициентов числителя таким образом, чтобы выполнялись условия приведенной выше теоремы об астатизме.
3.
БИБЛИОГРАФИЧЕСКИМ СПИСОК
Красовский А.А., Поспелов Г.С. Основы автоматики и технической кибернетики. - М.-Л.: Госэнергоиздат, 1962.
Семенов А.В., Геложе Ю.А. Моделирование системы слежения с астатизмом второго порядка // Материалы 4-й Международной научно-технической и научно-методической конференции «Проблемы современной системотехники». - Таганрог, 2010. - С. 58-62. Семенов А.В., Геложе Ю.А., Чуйков В.М., Семерников А.А. Автоматическая следящая система с астатизмом третьего порядка и цифровым управлением. Мехатроника, автоматизация, управление (МАУ-2009) // Материалы Международной научно-технической конференции. - Таганрог: Изд-во ТТИ ЮФУ, 2009. - С. 39-42.
4. Семенов А.В., Гайдук А.Р. Метод построения желаемых передаточных функций дискретных систем с высоким порядком астатизма // Известия ЮФУ. Технические науки. - 2013. - № 2 (139). - С. 14-20.
5. Гайдук А.Р. Теория автоматического управления: Учебник. - М.: Высшая школа, 2010.
6. Семенов А.В., Гайдук А.Р. Алгебраические условия астатизма дискретных систем. Радиоэлектронные средства передачи и приёма сигналов и визуализации информации // Материалы Второй Всероссийской конференции. - М.-Таганрог: Изд-во ТТИ ЮФУ, 2012. - С. 81-84.
7. Семенов А.В. Условия астатизма произвольного порядка дискретных систем управления на основе треугольника Паскаля // Математические проблемы современной теории управления системами и процессами: Материалы Международной молодежной конференции. - Воронеж: ИПЦ Научная книга, 2012. - С. 52-57.
Статью рекомендовал к опубликованию д.т.н., профессор Р. А. Нейдорф.
Гайдук Анатолий Романович - Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Южный федеральный университет»; e-mail: gaiduk_2003@mail.ru; 347904, г. Таганрог, ул. Слесарная 26, кв. 2; тел.: 88634626287; кафедра систем автоматического управления; д.т.н.; профессор.
Семенов Александр Валерьевич - e-mail: semenov-av@rambler.ru; 347922, г. Таганрог, ул. Шевченко, 2; тел.: 88634312350; ведущий конструктор.
Gaiduk Anatoly Romanovich - « Taganrog Institute of Technology - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education "Southern Federal University"; e-mail: gaiduk_2003@mail.ru; 26, Slesarnaya street, app. 2, Taganrog, 347904, Russia; phone: +78634626287; the department of automatic control systems; dr. of eng. sc.; professor.
Semenov Alexander Valerevich - e-mail: semenov-av@rambler.ru; 2, Shevchenko street, Taganrog, 347922, Russia; leading designer .
УДК 504.064.36
С.С. Алхасов, Л.П. Милешко, А.А. Целых
ОПРЕДЕЛЕНИЕ КОНЦЕНТРАЦИЙ ИОНОВ ТЯЖЁЛЫХ МЕТАЛЛОВ ПОСРЕДСТВОМ БЛОКА ОБРАБОТКИ ДАННЫХ МУЛЬТИСЕНСОРНОЙ СИСТЕМЫ ДЛЯ МОНИТОРИНГА ВОДНЫХ СРЕД
Рассматриваются мультисенсорные системы (МСС), которые уже начинают находить широкое применение в экологическом мониторинге. В МСС для мониторинга водных сред наиболее часто в качестве сенсоров используют ионоселективные электроды (ИСЭ). Блок обработки данных - программный компонент МСС. Он состоит из двух модулей. Первый разделяет вектор-столбцы разностей потенциалов на кластеры посредством сети Кохонена. Второй модуль выполняет количественную идентификацию с помощью радиальных базисных сетей. Обучающий массив концентраций формируется генератором псевдослучайных чисел. Для анализа точности вычислений созданы дополнительные скрипт-программы. Показано, что точность вычисления концентраций ионов зависит от величины коэффициента селективности ИСЭ. Даны рекомендации по подбору ИСЭ для МСС.
Мультисенсорная система; экологический мониторинг; тяжёлые металлы; ионоселективные электроды; искусственные нейронные сети; сети Кохонена; радиальные базисные сети.
