Авторефераты по всем темам  >>  Авторефераты по разным специальностям  

На правах рукописи

Федеральное государственное бюджетное образовательное учреждение

высшего профессионального образования

Национальный исследовательский университет МЭИ

Кошоева Бибигуль Бейшенбековна

Синтез алгоритмов численного дифференцирования,

работающих в реальном времени, с увеличенным шагом дискретизации для целей управления

Специальность 05.13.01. - ФСистемный анализ, управление и обработка информацииФ (по отраслям: энергетика, приборостроение, информатика, производственные процессы)

АВТОРЕФЕРАТ

диссертации  на соискание ученой степени

кандидата технических наук

Москва - 2012

Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Национальный исследовательский университет МЭИ

Научный руководитель:       доктор технических наук, профессор кафедры управление и информатика НИУ МЭИ

Колосов Олег Сергеевич

Официальные оппоненты:         доктор технических наук, профессор кафедры Автоматика  НИЯУ МИФИ

                        Рыбин Виктор Михайлович

кандидат технических наук, профессор кафедры автоматизированных систем управления тепловыми процессами Института тепловой и атомной энергетики НИУ МЭИ

                Волгин Владимир Владимирович

                       

Ведущая организация:                 ФГБОУ ВПО Московский государственный технический

                                       университет радиотехники, электроники и автоматики

Защита состоится  25.10.2012г. в  16  часов 00 мин. на заседании диссертационного совета Д 212.157.08 в государственном бюджетном образовательном учреждении высшего профессионального образования Национальный исследовательский университет УМЭИФ по адресу: 111250, г. Москва, Красноказарменная ул., д. 14, в малом актовом зале НИУ УМЭИФ.

С диссертацией можно ознакомиться в библиотеке НИУ МЭИ.

Автореферат разослан У  Ф _________________ 2012 года

Ученый секретарь

диссертационного совета Д 212.157.08

кандидат технических наук

доцент  Д.Н.Анисимов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

Актуальность работы

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

С другой стороны погрешность оценки производной во временной области не столь важна для обеспечения требуемого качества и устойчивости работы замкнутой системы с применением алгоритмов численного дифференцирования в стандартных регуляторах, таких, например, как ПД и ПИД-регуляторы.

В работах по автоматическому управлению справедливо отмечается, что как уменьшение шага дискретизации, так и усложнение алгоритма численного дифференцирования повышают точность оценки производной. Современная микропроцессорная техника, с помощью которой реализуются, алгоритмы управления, постоянно совершенствуется в части повышения быстродействия и соответственно,  уменьшения шага дискретизации при реализации алгоритмов управления численными методами. Вследствие этого даже простейшие алгоритмы численного дифференцирования по своим свойствам приближаются к идеальному дифференцированию. Это обстоятельство способствуют усилению высокочастотных помех за регулятором, приводит к излишне частому и неоправданному срабатыванию регулирующей аппаратуры, снижает качество и надежность работы систем. Для ослабления влияния высокочастотных шумов обычно вставляют отдельные фильтры, которые включаются последовательно с регуляторами. Однако последовательное включение фильтров с регуляторами порой существенно снижает быстродействие всей системы и увеличивает длительность переходного процесса. Поэтому, например,  на тепловых объектах с большими постоянными времени такое решение может привести к неприемлемо большой длительности и без того длительного переходного процесса.

Таким образом, актуальным является анализ и разработка методики синтеза структуры и параметров алгоритмов численного дифференцирования для систем управления. При этом важным является исследование следующих задач:

• оценка влияния длительности шага дискретизации алгоритма численного дифференцирования на устойчивость, качество и помехозащищенность проектируемой системы;
• привязка понятия требуемой точности алгоритма численного дифференцирования по отношению к рассматриваемой полосе пропускания системы;
• выбор структуры алгоритма численного дифференцирования применительно к особенностям проектируемой системы.

Целью диссертационной работы является разработка принципов построения, синтез структуры и параметров дискретных алгоритмов численного дифференцирования для одноконтурных систем управления с  увеличенным шагом дискретизации.

       Задачи исследования:

1. Анализ существующих подходов к реализации дискретных алгоритмов численного дифференцирования в частотной области  в составе одноконтурных систем управления при наличии шумов измерений и возмущений.
2. Проведение комплексного анализа частотных свойств алгоритмов численного дифференцирования на базе производной интерполяционного полинома Ньютона в зависимости от структуры алгоритма и шага дискретизации.
3. Разработка методики синтеза алгоритма дифференцирования в составе одноконтурных систем управления при наличии шумов измерений и возмущений.
4. Исследование эффективности алгоритма регулирования с точки зрения частоты переключений исполнительных механизмов при сохранении требуемой точности системы и качества функционирования локальной системы на примере существующей системы горячего водоснабжения на центральном тепловом пункте (ЦТП) с алгоритмом ПД-регулирования, в котором реализуется алгоритм численного дифференцирования с увеличенным шагом дискретизации.

Научная новизна:

1. Определены частотные свойства алгоритмов численного дифференцирования на базе производной интерполяционного полинома Ньютона в зависимости от структуры алгоритма и шага дискретизации, обеспечивающие уменьшения влияния высокочастотных помех в выходном сигнале алгоритмов управления на работу регулирующей аппаратуры.
2. Разработана методика синтеза алгоритмов численного дифференцирования в составе дискретных алгоритмов ПД и ПИД-регулирования для одноконтурных систем управления, определяющая структуру алгоритма дифференцирования и максимально допустимый его шаг дискретизации.

На защиту выносятся:

1. вопросы анализа и синтеза параметров дискретных алгоритмов дифференцирования в составе алгоритмов ПИД и ПД-регулирования, которые ориентированы на обеспечение необходимой величины запаса по фазе на частоте среза для ЛАЧХ разомкнутой системы.
2. методика оценки параметров алгоритмов дискретных дифференциаторов в составе алгоритмов регулирования, базирующаяся на требовании обеспечении заданного запаса по фазе на частоте среза разомкнутой системы.

       Практическая значимость:

Разработанная методика синтеза алгоритмов численного дифференцирования в составе дискретных алгоритмов ПД и ПИД-регулирования для одноконтурных систем управления, определяющая структуру алгоритма дифференцирования и максимально допустимый шаг дискретизации, которая позволяет повысить помехозащищенность систем управления, уменьшить число срабатываний регулирующей аппаратуры в единицу времени при сохранении точности регулирования.

Публикации. По результатам работы имеется 6 публикации. В том числе одна в журнале, входящем в перечень ВАК.

Апробация работы. Результаты работы обсуждены на следующих конференциях и семинарах:6-я научная конференция Управление и информационные технологии (УИТ-2010). Октябрь, 2010 г., Санкт - Петербург; XIX международный научно-технический семинар УСовременные технологии в задачах управления, автоматики и обработки информацииФ. Сентябрь, 2010 г., Алушта; Международная научно-техническая конференция Прикладная математика и механика: проблемы и перспективы. Ноябрь, 2010г, г.Бишкек, Кыргызстан; XX международный научно-технический семинар УСовременные технологии в задачах управления, автоматики и обработки информацииФ. Сентябрь, 2011 г., Алушта; Международная конференция Информационные технологии и математическое моделирование в науке, технике и образовании. Октябрь, 2011г., г.Бишкек, Кыргызстан.

Структура и объем работы. Диссертация состоит из введения, 4 глав, заключения, списка литературы из 70 наименований и приложения, изложенных на 150 страницах машинописного текста, содержит 55  рисунков, 33 таблицы.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ:

Во введении приводится актуальность темы, рассматриваются предпосылки для проведения исследований, приводится научная новизна и практическая значимость работы.

В первой главе диссертации рассматриваются алгоритмы управления, в которых чаще всего используются алгоритмы численного дифференцирования - ПД и ПИД-алгоритмы регулирования: особенности их использования в локальных системах управления, реализация операции дифференцировании в непрерывных и дискретных алгоритмах регулирования, влияние возмущений на качество работы систем с алгоритмами регулирования.

В основном методики настройки параметров алгоритмов ПД и ПИД-регулирования базируются на анализе частотных характеристик замкнутой системы путем анализа показателя колебательности М системы, запаса устойчивости, функции чувствительности и т.д. Такой подход фактически разделяет синтез оптимальных настроек регулятора системы при отработке заданного вида типового входного воздействия (обычно это скачок) и вопросы минимизации влияния высокочастотных помех.

Существуют также методы, основанные на анализе амплитудных и фазовых частотных характеристиках разомкнутой системы, когда по виду частотных характеристик в областях низких, средних и высоких частот проще оценить влияние разных вариантов настроек. Последние более удобны и особенно для систем, настраиваемых на работу без перерегулирования или с небольшим (не более 20%) перерегулированием. Для таких систем основным показателем является необходимый запас по фазе на частоте среза разомкнутой системы. Подобная ориентация на обеспечение надлежащего запаса по фазе на частоте среза разомкнутой системы позволяет четче синтезировать требуемую точность системы и обеспечить практическую реализацию дифференцирующей части алгоритма управления. Любые нюансы в частотном представлении различных объектов управления несложно проследить и учесть в окрестности частоты среза разомкнутой системы, не прибегая к утомительной процедуре отыскания оптимальных соотношений в коэффициентах настроек регулятора.

В общем виде структура системы с объектом управления Р, регулятором R, сигналом уставки r, управляющим сигналом u, возмущающим воздействием d, высокочастотным шумом измерений n и выходным сигналом y показана на рис.1.

Обычно d относят к относительно низкочастотным возмущениям, которые влияют на точность системы и в выходном сигнале системы y проявляются тем меньше, чем выше общий коэффициент передачи системы по контуру RP на низких частотах. Однако выходной сигнал регулятора u имеет составляющую сигнала возмущения, который не ослабляется контуром регулирования во всем диапазоне частот. Это приводит к нежелательному износу регулирующей аппаратуры из-за ее частого срабатывания. Аналогичное влияние на выходной сигнал регулятора на высоких частотах оказывают высокочастотные шумы. Высокочастотная помеха (шум измерений) n в области высоких частот не ослабляется замкнутым контуром системы, если модуль . По этой причине в состав дифференцирующей части алгоритма регулирования вводят фильтр, ограничивающий усиление в области высоких частот.

Рис.1. Общий вид структуры исследуемой системы

Вместе с тем, известно, что частотные характеристики дискретных алгоритмов дифференцирования на базе первых слагаемых производной интерполяционного полинома Ньютона так же ограничивают коэффициент передачи дифференциатора в области высоких частот и кроме этого существенно уменьшают коэффициент передачи его в окрестностях частот кратных частоте дискретизации . То есть фильтрующая способность алгоритмов численного дифференцирования (АЧД) в области высоких частот (при рациональном выборе шага дискретизации) может оказаться выше, чем у алгоритма дифференцирования с бесконечно малым шагом в связке с фильтром низких частот.

Указанное свойство АЧД может позволить отказаться от постановки дополнительного фильтра низких частот в состав дискретного дифференциатора. Функции дифференциатора и одновременно фильтра низких частот с более высоким качеством может выполнять АЧД с большим интервалом дискретизации. Однако интервал этот не может быть слишком большим из-за риска потери устойчивости системой. Важным моментом при этом является также обоснование целесообразности усложнения структуры самого алгоритма численного дифференцирования при стремлении увеличить точность операции взятия производной входного сигнала в алгоритме управления.

Вторая глава диссертации посвящена алгоритмам численного дифференцирования, работающих в реальном времени, и их свойствам. В главе показывается, что известные алгоритмы численного дифференцирования базируются на определенных суммах первых слагаемых производной интерполяционных полиномов Бесселя, Ньютона или Стирлинга. Однако для целей управления в режиме Уon-lineФ подходят алгоритмы взятия производной на основе интерполяционного полинома Ньютона, так как в них используются только предыдущие отсчеты.

Остановимся подробнее на интерполяционном полиноме Ньютона по переменной для функции f(t):

                       (1)

При t=t0 (u=0)

первая производная интерполяционного полинома Ньютона (1) имеет вид:

(2)

Дальнейшие преобразования позволяют получить ряд АЧД, работающие в реальном времени на базе сумм различного числа первых слагаемых ряда (2) (Табл.1).

Таблица 1

Алгоритмы численного дифференцирования

на базе производной полинома Ньютона

Число слагаемых	Алгоритм
1
2
3
4
5
6

Сравнение свойств АЧД (Табл.1) в частотной области по их логарифмическим амплитудно-частотным (ЛАЧХ) и фазо - частотным (ФЧХ) характеристикам позволяет оценить точность дифференцирования в полосе пропускания по отношению к идеальному дифференцированию и установить влияние интервала дискретизации на точность оценки производной.

Под полосой пропускания условно будем понимать область низких частот разомкнутой ЛАЧХ системы, ограниченной частотой среза (). В полосе пропускания частотные свойства АЧД должны с определенной точностью совпадать частотными свойствами идеального дифференцирования. Для такого АЧД и определяются оптимальные настройки ПИД-регулятора по известным методикам. Например, с использованием показателя колебательности М замкнутой системы.

За полосой пропускания, когда замкнутый контур системы не оказывает ослабляющего влияния на уровень высокочастотных помех, представляют интерес сами частотные свойства АЧД в области высоких частот с точки зрения усиления или ослабления помех.

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

      (3)

где - коэффициент передачи дифференцирующего звена.

Не нарушая общности в данном разделе принят равным единице.

Частотные характеристики вычисляются путем последовательного применения преобразований Лапласа и Фурье к алгоритму.

Комплексный коэффициент усиления АЧД  на базе алгоритма с использованием первого слагаемого (Тал. 1) будет иметь вид:

  (4)

Далее анализируется модуль и фаза АЧД (или дискретного дифференциатора) на базе первого слагаемого от производной интерполяционного полинома Ньютона (далее для краткости будем писать: на базе одного слагаемого).

В работе получено и обосновано достаточно простое условие для оценки верхней частоты АЧД (4), для которой точность совпадения модуля частотной характеристики (ЧХ) АЧД по сравнению с идеальным будет не хуже - 0.04%, а по фазе - 1.3%:

                (5)

На рис.2 показаны частотные характеристики АЧД на базе одного слагаемого при интервалах дискретизации и . А на рис.3 Годографы дискретных алгоритмов дифференцирования с разным числом слагаемых.

Рис.2. ЛАЧХ и ФЧХ дискретного алгоритма дифференцирования на базе одного слагаемого для интервала дискретизации: =1с и =0.1с (кривые 1 и 2)

Анализ построенных частотных характеристик демонстрирует естественную их периодическую повторяемость в интервалах частот, задаваемых границами:

,  где n = 1, 2, 3,  . Эта же периодичность иллюстрируется годографами алгоритмов дифференцирования, показанными на рис.3, где отмечены      

                       а)                                                        б)

                                                       в)

Рис.3. Годографы дискретных алгоритмов дифференцирования при с: а) на базе 1, 2, 3-х слагаемых, б) на базе 1,2,3,4,5,6 слагаемых; метки на годографах соответствуют: (), в)увеличенный участок  для .

точки на годографах, соответствующие правым границам возможной полосы пропускания полосы пропускания определяемой (5) и экстремумам ЛАЧХ . Несложно показать, что эти экстремумы во всех случаях оказываются на положительной части действительной оси комплексной плоскости. За границей полос пропускания, в окрестностях частот , где n = 2, 3,    происходит заметное ослабление высокочастотных помех. Следует отметить, что с уменьшением шага дискретизации растет коэффициент передачи на высоких частотах. Аналогичный эффект наблюдается у дискретных дифференциаторов с ростом числа слагаемых в алгоритме.

На рис. 4 и 5 показаны относительные погрешности АЧД с разным числом слагаемых

по отношению к идеальному дифференцированию для граничной частоты (5) и для менее жесткого условия :

Рис.4. Оценка погрешности АЧД по модулю

Рис.5. Оценка погрешности АЧД по фазе

Для оценки свойств АЧД в области высоких частот предлагается сопоставить амплитудно-частотную характеристику АЧД с аналогичной характеристикой, которой обладает непрерывное реальное дифференцирующее звено с комплексным коэффициентом передачи (6), а так же АЧД с малым шагом дискретизации и фильтром низких частот на выходе:

, (6)

где эквивалентная постоянная времени эквивалентного непрерывного реального дифференцирующего звена. Коэффициент соответственно для АЧД с одним, двумя и тремя слагаемыми будет - 2, 4, 6,67. На рис.6 представлены соответствующие частотные характеристики.

Рис.6. Частотные характеристики АЧД на базе двух слагаемых для и эквивалентного непрерывного дифференцирующего звена с (соответственно кривые 1 и 2).

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

Рис.7 иллюстрирует более высокую эффективность АЧД с разным числом слагаемых по сравнению с непрерывным аналогом при пропускании сигнала типа Белый шум

Рис.7. СКО сигналов на входах и выходах эквивалентного аналогового и дискретного (с разным числом слагаемых) алгоритмов дифференцирования.

В третьей главе  рассматривается методика синтеза АЧД, работающих в реальном времени с увеличенным шагом дискретизации для уменьшения влияния высокочастотных помех в тракте передачи управляющего сигнала и уменьшения числа срабатываний исполнительного элемента при сохранении требуемой точности и качества системы управления.

Методика разработана и ориентирована на синтез в частотной области по характеристикам системы в разомкнутом состоянии и может быть использована для двух возможных условий работы систем:

- Для систем с обеспечением максимального быстродействия и отсутствия перерегулирования (или с перерегулированием не более 20%) при отработке скачка. В этом случае объект управления может описываться упрощенными передаточными функциями либо в виде двух апериодических звеньев, либо в виде одного апериодического звена и звена запаздывания.

- Для систем с перерегулированием более 20% при отработке скачка. При этом повышаются требования к адекватности математического описания объекта управления с использованием современных методов структурной и параметрической идентификации.

На рис.8 представлена обобщенная структурная схема системы промышленной автоматики. Как правило, объекты промышленной автоматики представляют собой объекты с самовыравниванием. В работе показывается, что, если накладывать ограничение на динамику замкнутой системы в виде требования максимального быстродействия и перерегулирования не более 20%, объекты могут описываться в виде упрощенных передаточных функций вида (7) с параметрами, определяемыми по кривой разгона объекта, известными приемами:

или ,  (7)

где - коэффициент усиления объекта; - время запаздывания; То - постоянная времени разгона,

Объясняется это тем, что для таких систем необходим запас по фазе на частоте среза порядка . Поэтому суммарный фазовый сдвиг на частоте среза от запаздывания и от инерционных звеньев находится в пределах и, следовательно, сдвиг от инерционного звена с постоянной времени и от такого же запаздывания будут практически одинаковыми.

Рис.8. Обобщенная структурная схема исследуемой системы,

где X - уставка, Y - регулируемая координата, Z - регулирующее воздействие, F - возмущающие воздействия.

Частотные характеристики разомкнутых систем с объектами описываемые передаточными функциями (7) совпадают в окрестности частоты среза и гарантируют адекватное представление замкнутой системы.

Исполнительный орган для таких объектов представляет собой обычно некоторую управляющую (дозирующую) аппаратуру, которая приводится в действие исполнительным двигателем. Таким образом, к передаточной функции самого объекта добавляется передаточная функция исполнительного двигателя, которая содержит обязательно интегрирующее звено и, как минимум, еще и инерционное звено.

- передаточная функция двигателя совместно с регулирующей аппаратурой;

Поскольку постоянные времени исполнительного двигателя и датчика обычно много меньше постоянных времени объекта, то в анализе ими пренебрегают. В работе анализируется ряд систем в общем виде, начиная со второго порядка, и систем, с описанием объектов упрощенными передаточными функциями (7).

Проведено сравнение частотных характеристик алгоритмов ПД-регулирования и разомкнутых систем для вариантов применения дискретных и эквивалентных аналоговых АЧД в составе алгоритмов регулирования в анализируемых системах и показано совпадение качественных показателей таких систем благодаря совпадению их характеристик в окрестности частоты среза при условии обеспечения заданного запаса по фазе.

Управляющий сигнал Z воздействует на регулирующую аппаратуру либо непосредственно через пропорциональный усилитель мощности (линейный вариант системы), либо через трехпозиционное реле, включающее и выключающее исполнительный двигатель. Эти два случая представления регулирующей аппаратуры показаны на рис.9:

 

               а)                                                б)

Рис.9. Два варианта представления регулирующей аппаратуры.

Важнейшим качественным показателем, характеризующим работу системы, показанной на рис.8 является частота срабатывания регулирующей аппаратуры (рис.9) как в релейном, так и в непрерывном варианте. Последнее утверждение подразумевает наличие в структуре непрерывной системы определенной зоны нечувствительности, выход за пределы которой приводит к замыканию всего контура управления. Такая зона в принципе может быть как угодно малой и влияние высокочастотных помех можно оценивать по числу переходов через ноль сигнала управления Z (для чисто линейной системы).

Передаточная функция исследуемой системы относительно управляющего сигнала Z для области высоких частот имеет вид:

(8)

Учитывая требуемый запас по фазе, который гарантирует отсутствие перерегулирования и максимальное быстродействие в непрерывной линеаризованной системе, на примере конкретной системы иллюстрируется методика выбора параметров алгоритма ПД-регулирования для трех различных величин шага дискретизации: . Анализируемая нами система с рассматриваемыми параметрами описывается в разомкнутом состоянии последовательным соединением функциональных блоков с передаточными функциями:

;  ; ;        , ,

где - передаточная функция эквивалентного непрерывного ПД- регулятора, - передаточная функция дискретного ПД-регулятора, - ПФ объекта, - ПФ регулирующей аппаратуры, - ПФ датчика. В таблице 2 представлены частоты среза и запасы по фазе на них разомкнутых систем как для дискретных алгоритмов и для их аналогов. На рис. 10 представлены частотные характеристики разомкнутых систем. Показанные результаты свидетельствуют о том, что чем сложнее структура алгоритма дифференцирования, тем больший запас по фазе на частоте среза имеет система.

Далее в работе исследовано влияние случайных возмущений, приложенных к выходу управляемого объекта для систем с ПД-регуляторами, содержащими как  дискретные, так и эквивалентные аналоговые дифференциаторы, и показана более высокая эффективность применения дискретных дифференциаторов с увеличенным шагом благодаря уменьшению частоты срабатывания управляющей аппаратуры. Результаты этих исследований представлены на рис.11.

Таблица 2

Запасы по фазе на частотах среза разомкнутых систем при =2с.

Число слагаемых дискретного дифференциатора	Система с эквивалентным непрерывным дифференциатором	Система с дискретным дифференциатором
1сл	=0,24 =76,3о	=0,24 =76,6о
2сл	=0,238 =82,2о	=0,24 =88о
3сл	=0,235 =84,5о	=0,235 =88,9о

Рис.10. Частотные характеристики разомкнутых систем: кривая 1-система, где алгоритм ПД-регулирования реализован в виде непрерывного эквивалентного алгоритма дифференцирования; кривая 2 - система, где алгоритм ПД-регулирования реализован в виде дискретного алгоритма дифференцирования.

На рис.12 показано число срабатывания реле в той же систем. Но с трехпозиционным реле на входе исполнительной аппаратуры.

Результаты моделирования показывают, что при прохождении шумовых составляющих с известными спектральными свойствами через алгоритм ПД-регулирования наибольшее подчеркивание высоких частот производят непрерывные аналоги дискретных дифференциаторов, что приводит к увеличению частоты срабатывания исполнительной аппаратуры.

                       а)                                                б)

Рис.11.  Число переходов через ноль управляющего сигналаZ при Тф=5с  (-w- Непрерывный аналог алгоритма  ПД-регулирования,  -б- Дискретный алгоритм ПД- регулирования): а) =2с; б) =1с.

Результаты моделирования, представленные на Рис.11а,б показывают, что при =2с и =1с в системе с эквивалентным непрерывным аналогом ПД-регулятора число переходов через ноль управляющего сигнала больше, чем в системе с дискретным дифференциатором.

                       а)                                                б)

Рис.12.  Число срабатывания реле при Тф=5с (-w- Непрерывный аналог алгоритма ПД- регулирования,  -б- Дискретный алгоритм ПД-регулирования): а); =2с б) =1с.

Сформулирована методика синтеза АЧД с увеличенным шагом дискретизации в составе алгоритмов ПД и ПИД -регулирования для объектов промышленной автоматики:

1. На основании результатов идентификации отдельных элементов структуры и анализа частотных свойств разомкнутой линеаризованной системы известными способами производится синтез структуры алгоритмов ПД и ПИД Црегулирования и параметров его настроек с учетом требований технического задания.

При требовании отсутствия перерегулирования (или перерегулирования не более 20%) и максимального быстродействия при отработке ступенчатого входного воздействия запас по фазе на частоте среза разомкнутой системы должен быть более . При синтезе параметров настроек описание объекта допускается упрощенными передаточными функциями в виде последовательного соединения двух инерционных звеньев или одного инерционного звена и звена запаздывания.

Если для проектируемой системы допускается относительно высокий показатель колебательности в замкнутом состоянии, то в этом случае не следует описывать идентифицируемый объект упрощенными передаточными функциями, а необходимо воспользоваться более точными известными методами структурной и параметрической идентификации. Для этого варианта известными способами необходимо оценить минимально допустимый запас по фазе на частоте среза разомкнутой системы с учетом полученного математического описания по результатам идентификации объекта управления.

2. В случае использования непрерывного или релейного исполнительного устройства в сочетании с приводным реверсивным двигателем в анализируемой структуре появляется после ПД и ПИД - регулятора дополнительный интегратор. Частотные свойства последнего должны быть учтены заранее, на первой стадии анализа частотных свойств разомкнутой системы, синтеза структуры и коэффициентов настроек алгоритма ПД или ПИД - регулирования.
3. Для уменьшения влияния высокочастотных помех выбирается постоянная времени фильтра низких частот на выходе алгоритма численного дифференцирования. В непрерывном варианте это определяет выбор параметров реального дифференцирующего звена в канале дифференцирования.
4. На входе усилительно-преобразовательного устройства вводится зона нечувствительности , определяемая требованиями ТЗ по точности регулирования в стационарном режиме работы  в случае использования как непрерывного, так и  релейного исполнительного элемента.
5. В зависимости от ширины зоны нечувствительности находятся нормированные уровни срабатывания релейного элемента с использованием соотношения , где а - модуль зоны нечувствительности. При таком выборе уровней выходного сигнала реле максимум эквивалентного комплексного коэффициента передачи данного трехпозиционного реле не превышает единицы и с его учетом общий коэффициент передачи разомкнутой системы не должен превышать величину разомкнутого коэффициента передачи линеаризованной системы, найденного на первом этапе.
6. Для повышения эффективности регулятора в части подавления высокочастотных помех и упрощения его структуры  дифференцирующая часть, реализует алгоритм численного дифференцирования на базе первых (от одного до трех) слагаемых  производной интерполяционного полинома Ньютона с увеличенным шагом, определяемым из соотношения: , где k=2, 4, 6.67, - соответственно для алгоритма дифференцирования на базе одного, двух и трех слагаемых. Фильтр низких частот при этом из регулятора убирается.

В четвертой главе проводится синтез АЧД с увеличенным шагом дискретизации  для системы управления горячим водоснабжением (ГВС) на центральном тепловом пункте (ЦТП) по разработанной методике (ГВС-1) и проводится сравнение эффективности его работы в составе системы с результатами работы алгоритма ПД -регулирования для такой же системы, описанной в  работах журнала МАУ за 2010г (ГВС-2). Прототипом системы ГВС рассматриваются системы ряда ЦТП в Юго-Восточном округе г. Москвы.

Определены для систем ГВС-1 и ГВС-2 частоты среза и запасы по фазе при различных и сведены в таблицу 3. Проведенные ранее исследования системы ГВС-2 показывают, что для минимизации частоты срабатывания исполнительных механизмов в системе регулирования температуры ГВС на ЦТП  необходимо использовать дифференциатор на базе 3 слагаемых при шаге дискретизации или .

       На основании соотношения (5)        для системы ГВС-1  с шаг дискретизации должен быть не больше 7с, а для системы  ГВС-2 с должен быть .

Таблица 3.

Частоты среза и запасы по фазе

для систем  ГВС-1 и ГВС-2

Число слагаемых дискретного дифференциатора	Система ГВС-1		Система ГВС-2
	=2.2с	=7с	=2.2с	=7с
1 сл	=0.014 =78о	=0.015 =76.68о
3 сл			=0.036 =58.6о	=0.037 =61.65о

Таблица 4

Число срабатываний реле систем ГВС в

варианте построения исполнительной

аппаратуры с трехпозиционным реле,

tн=250-500с, реле:

Число слагаемых дискретного дифференциатора	Система ГВС-1		Система ГВС-2
	=2.2с	=7с	=2.2с	=7с
1 сл	9	6
3 сл			52	25

Из представленных таблиц следует, что:

- при увеличенном шаге дискретизации число срабатываний регулирующей аппаратуры уменьшается;

- число срабатываний в системе ГВС-2 с АЧД на базе 3 слагаемых как при =2.2с, так и при =7с, намного больше, чем в системе ГВС-1 с АЧД на базе 1 слагаемого.

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

Заключение.

1. В диссертационной работе проведен комплексный анализ свойств дискретных алгоритмов численного дифференцирования в частотной области, реализуемых на базе суммы, содержащей до шести первых слагаемых производной интерполяционного полинома Ньютона.
2. Показано, что дискретные алгоритмы дифференцирования, работающие в реальном времени в составе алгоритмов управления и реализуемые на базе ограниченного числа первых слагаемых производной интерполяционного полинома Ньютона, представляют собой как простейшие, известные алгоритмы дифференцирования, так и более точные и сложные алгоритмы.
3. Предложено оценивать точность дифференцирования для дискретных алгоритмов численного дифференцирования в составе алгоритмов регулирования замкнутых систем управления по максимальному отклонению частотных характеристик дискретного алгоритма дифференцирования от идеального в полосе нижних частот, ограниченной частотой среза разомкнутой системы.
4. Показано, что как уменьшение шага дискретизации дискретного алгоритма численного дифференцирования, так и увеличение суммы числа первых слагаемых производной интерполяционного полинома Ньютона в дискретном алгоритме численного дифференцирования увеличивает точность дифференцирования в заданной полосе нижних частот, ограниченной частотой среза разомкнутой системы, и одновременно поднимает коэффициент усиления для высоких частот.
5. Найдены условия и предложена методика оценки параметров дискретных алгоритмов численного дифференцирования в составе алгоритмов регулирования, базирующаяся на требовании обеспечении заданного запаса по фазе на частоте среза разомкнутой системы. Реализуемый при этом максимально допустимый шаг дискретизации уменьшает влияние высокочастотных помех на сигнал управления по сравнению с непрерывным аналогом подобного алгоритма дифференцирования, а также с вариантом использования дискретного алгоритма численного дифференцирования с малым шагом дискретизации совместно с фильтром низких частот.
6. Предложенная методика синтеза параметров дискретных алгоритмов дифференцирования позволяет пользоваться упрощенными аппроксимирующими передаточными функциями объекта управления при требовании отсутствия в системе перерегулирования и быстродействия близкого к максимальному.
7. Предложенная методика синтеза параметров дискретных алгоритмов дифференцирования позволяет отказаться от дополнительного фильтра высоких частот в составе алгоритмов регулирования и понизить чувствительность системы к высокочастотным помехам за счет определенного увеличения шага дискретизации в канале дифференцирования.
8. В системе ГВС использование предложенной методики синтеза дискретного алгоритма дифференцирования на базе первого слагаемого производной интерполяционного полинома Ньютона в составе алгоритма ПД - регулирования с увеличенным шагом дискретизации позволяет сократить число срабатываний регулирующей аппаратуры под действием случайных возмущений, что повышает ресурс работы системы.

Основные положения диссертации достаточно полно представлены в следующих публикациях:

1. Кошоева Б. Б., Колосов О.С.. Алгоритмы численного дифференцирования реального времени для задач автоматизации и управления. Ц Новые технологии, УМехатроника, автоматизация, управлениеФ, 2012, №2. Ц С.10-15.
2. Кошоева Б.Б., Колосов О.С. Выбор шага дискретизации в цифровых алгоритмах управления с дифференцированием. // Материалы 6-й научной конференции Управление и информационные технологии (УИТ-2010). - СПб.: ОАО Концерн ЦНИИ Электроприбор, 2010. С. 210 - 214.
3. Кошоева Б.Б., Колосов О.С. Алгоритмы численного дифференцирования реального времени в алгоритмах управления. Труды XIX международного научно-технического семинара УСовременные технологии в задачах управления, автоматики и обработки информацииФ. Сентябрь, 2010 г., Алушта. - М.: Издательский дом МЭИ, 2010. - С. 285-286.
4. Кошоева Б.Б. Алгоритмы  численного дифференцирования реального времени на базе интерполяционного полинома Ньютона. // Известия КГТУ им. И. Раззакова №22. Материалы Международной научно-технической конференции Прикладная математика и механика: проблемы и перспективы.  2011, г.Бишкек, Кыргызстан.- С. 259-265.
5. Кошоева Б.Б., Колосов О.С. Особенности синтеза алгоритмов численного дифференцирования реального времени для систем управления. - Труды XХ международного научно-технического семинара УСовременные технологии в задачах управления, автоматики и обработки информацииФ. Сентябрь, 2011 г., Алушта. - Пенза: Изд-во ПГУ, 2011. - С.114-115.
6. Кошоева Б.Б. Структуры и свойства алгоритмов численного дифференцирования реального времени.//Известия КГТУ им. И. Раззакова №24. Материалы Международной конференции Информационные технологии и математическое моделирование в науке, технике и образовании. Кыргызстан, г.Бишкек, 5-9 октября 2011. -С.48-54.

Авторефераты по всем темам  >>  Авторефераты по разным специальностям