Генетическая настройка цифрового регулятора частоты вращения с двухрежимной фильтрацией

НПП - название подпункта/параметра НКП; ед. -единица измерения НПП; М - значение, соответствующее НПП.

Тогда текст описания будет иметь следующую структуру:

НКЩ: НПП1(ед.1)=М1;

НКП2: НПЩед.2)=М2;

...; НКПП: НППп(ед.п)=Мп.

Например, если вышел из строя блок аппаратуры, произошел отказ ПРД или АС, то в этом поле будут приведены перечень названий неработоспособных приборов и номера их комплектов.

Реализация когнитивного дополнения

Когнитивное дополнение интерфейса реализовано в среде разработки CodeGear Delphi 2007 for Win32. Необходимую скорость выполнения данного приложения обеспечивает использование пакета Graphics32 с открытыми исходными кодами. Пакет содержит набор функций, классов, компонентов и элементов управления, разработанных для высокоэффективного программирования графическими средствами. Дерево журнала событий создано при помощи компоненты Virtual Treeview с открытыми исходными кодами, являющейся одним из самых гибких и доступных средств создания контроля элементов управления в виде дерева и/или списка на сегодняшний день.

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

ностики и прогнозирования состоянии технических систем НС КИС. На основе интерфейса и его когнитивного приложения оператор визуально оценивает состояние НС КИС и окружающей среды. Целесообразность разработанных форм индикации состоит в обеспечении надежного и эффективного взаимодействия ЛПР с ЭВМ, в повышении оперативности за счет сокращения аналитической работы, компактности и наглядности геометрических форм представления информации. Экспериментальное программное обеспечение, реализующее когнитивное дополнение, встроено и испытано в составе НС КИС нового поколения [3].

Литература

1. Энциклопедия космонавтики. URL: http://kosmonavti-ca.ru/tag/kompleks-upravleniya/ (дата обращения: 23.03.2009).

2. Развитие наземного комплекса управления космическим аппаратом «KazSat» / Ю.М. Урличич [и др.] // Развитие космической деятельности в Республике Казахстан на 20052007 годы: матер. Междунар. конф., посвященной итогам выполнения Гос. прогр. Алматы. 2007. С. 19-20.

3. Смирнов С.В. Средства когнитивной графики для отображения текущего состояния наземных станций командно-измерительных систем // Авиакосмическое приборостроение. 2008. № 3. С. 47-57.

4. Вагин В.Н., Еремеев А.П. Базовые принципы конструирования интеллектуальных систем поддержки принятия решений реального времени для мониторинга и управления сложными техническими объектами // Использование методов искусственного интеллекта «...» (Переславль-Залесский, 26-27 ноября 2003 г.): тр. Третьего расшир. сем. М.: Физмат-лит, 2003. С. 79-97.

5. Хачумов В.М., Ксенофонтова Е.В. Образный анализ и диагностика сложных процессов // Математические методы распознавания образов (ММРО-11): докл. 11-й Всерос. конф. (Пущино, 23-29 ноября 2003 г.). М., 2003. С. 201-204.

ГЕНЕТИЧЕСКАЯ НАСТРОЙКА ЦИФРОВОГО РЕГУЛЯТОРА ЧАСТОТЫ ВРАЩЕНИЯ С ДВУХРЕЖИМНОЙ ФИЛЬТРАЦИЕЙ

Ф.В. Калугин

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

Ключевые слова: двухрежимный фильтр, настройка ПИД-регулятора, генетический алгоритм, коленчатый вал.

Необходимость обеспечения устойчивости холостого хода автомобильных дизелей, оснащенных электронными системами управления, обусловливает создание специальных систем управления (СУ). Данная проблема решалась и при доводке электронной системы ЭСУ-1 на дизелях серий КАМАЗ-740.60-360, КАМАЗ-740.61-320, КА-МАЗ-740.62-280, поскольку двухкаскадный пропорционально-интегрально дифференциальный (ПИД) регулятор автоматически не обеспечивал устойчивую работу дизеля при отключении на-

грузки [1]. Алгоритмы управления получаются более эффективными, если вместо специальной СУ холостого хода применять универсальные, то есть обеспечивающие устойчивую работу дизеля на всех нагрузочно-скоростных режимах, включая и режимы холостого хода, и торможения двигателем. В настоящей работе данная проблема решается введением модифицированного цифрового фильтра оборотов частоты [2] с переменной структурой и универсальной генетической настройкой ПИД-регулятора.

Система управления дизельным двигателем

Принципиальная схема регулирования частоты вращения коленчатого вала показана на рисунке 1. Текущая частота вращения О определяется величиной подачи топлива (положение рейки топливного насоса Ь), свойствами двигателя и нагрузкой (момент инерции нагрузки, момент сил внешнего сопротивления М2), факторами окружающей среды.

Измеряемое значение О подается на вход фильтра высокой частоты, затем отфильтрованная частота вращения ОР передается в качестве входного сигнала в цифровой регулятор частоты вращения двигателя. Фильтр высокой частоты может отсутствовать в контуре управления двигателем. Цифровой регулятор задает требуемое положение рейки Ьд для обеспечения стабильности скоростного режима, который определяется настройкой задаваемой частоты вращения О§. Исполнительный механизм (ИМ) устанавливает физическое положение рейки Ь, близкое к задаваемому значению Ьд. Текущая величина Ь отличается от Ьк вследствие инерции ИМ. Для стабилизации параметра Ь в ИМ используется обратная связь с цифровым регулятором положения рейки. Модификация цифрового двухрежимного фильтра оборотов частоты заключается в том, что не допускается большое отклонение отфильтрованного сигнала от реального: при разнице отфильтрованного сигнала от реального больше некоторого порога, характерного для двигателя, в качестве результата берется реальная частота.

Математическая модель

системы регулирования

Уравнение частоты вращения двигателя имеет

вид: I-=МЕ ,

где I - суммарный момент инерции двигателя и нагрузки, приведенной к коленчатому валу; Ь -время; МЕ - суммарный момент приложенных сил. МЕ можно представить в виде:

МЕ = М1 + М2 + М3,

где М1 - крутящий момент; М2 - момент сил внешнего сопротивления; М3 - момент сил внутреннего трения. Мгновенное значение крутящего момента можно представить в виде:

МД1)=ни-^с(1)пе + М'(е),

ПТа

где Ни - теплотворная способность топлива; ^ -число цилиндров в двигателе; Та - тактность двигателя; gc - цикловая подача топлива; пе - эффективный КПД; М'(0 - флуктуационная составляющая (шум двигателя), среднее по времени значение М'(0 равно нулю [3]. Цикловая подача топлива зависит от частоты вращения О и координаты рейки топливного насоса Ь. В первом приближении будем считать gc пропорциональной координате рейки и не зависящей от частоты:

^ = к8Ь,

где кё - коэффициент пропорциональности. С учетом

К=Н.^4,

и ПТ„ ,е'

где К1 - коэффициент пропорциональности крутящего момента к координате рейки; Т1 - запаздывание (время, необходимое на движение и сгорание топлива). Предположим, что момент сил внутреннего трения пропорционален частоте вращения М3 =-0О , где 0 - коэффициент сопротивления.

С учетом (2)-(5) уравнение момента импульса (1) принимает вид:

1й°

Уравнение регулятора частоты вращения

Для регулирования частоты оборотов О используем классический ПИД-регулятор

Ьк+1 = Кр¥+ К, X>+ КвФ,

где нижний индекс к соответствует характеристике в момент времени Ь=кЛе ; Ль - шаг квантования по времени; Кр, К1, Ко - постоянные времени регулятора (пропорциональная, интегральная и дифференциальная); -О, где О8 - зада-

ваемая угловая скорость (частота).

Прямое численное моделирование системы регулирования частоты вращения

На основе уравнений (6) и (7) строим систему уравнений, которая описывает замкнутую систему, состоящую из двигателя и регулятора (рис. 1). Эта система при единичном скачке на входе определяет кривую переходного процесса в системе.

Численное моделирование проводим следующим образом:

1) устанавливаем первые йь

•Qk-1 и

равные 0, где к т= round

/ \ х

T

Ч T

2) вычисляем частоту на к-м шаге =А№к_!+Лк-х;

3) вычисляем положение рейки на шаге к+1, Ьк+1=КрЧ+К, X>+КвЧ;

4) вычисляем приращение частоты на к-м ша-

„, К1Ьк_кт+ Мк + м2к -ео,

ге Чк =-^-А1, где в качестве Мк берем случайные величины в интервале [-2, 2];

5) так как задачей нашего исследования является имитация условий холостого хода, Мк=0 (внешняя нагрузка отключена);

6) увеличиваем шаг к, переходим к пункту 2.

Вследствие наложения на работу двигателя

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

Оценка качества переходного процесса

Для оценки используем 4 параметра [4] (вид переходной характеристики и ее параметров изображен на рисунке 2):

1) 1;н - время нарастания, то есть время, за которое переменная у(0 возрастает с 0,1 до 0,9 установившегося значения;

2) 8 - «выброс» (максимальное превышение сигналом у(0 единичного уровня);

3) 1;3 - время затухания переходного процесса (время между моментом первого достижения сигналом у(0 единичного уровня и моментом, начиная с которого значения у(0 остаются внутри интервала [1±е]; е - некоторая постоянная);

4) - среднеквадратичное отклонение от задаваемого сигнала, рассчитывается на отрезке затухания переходного процесса.

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

Оценка переходной кривой по упомянутым выше параметрам (1н, 8, 1;3, представляет собой многокритериальную задачу. Для ее решения воспользуемся предложенным Р. Беллманом и Л. Заде методом слияния целей и ограничений [5]. Введем лингвистические переменные «время нарастания», «время затухания», «выброс», «равномерность затухания» и определим на каждой из них терм «приемлемое значение» как нечеткое множество с функцией принадлежности р,(х) трапецеидального вида. Пример функции принадлежности для терма «приемлемое значение» лингвистической переменной «выброс» приведен на рисунке 3.

Конкретные значения параметров переходной кривой (1н, 8, 1;3, Б() могут быть охарактеризованы степенью принадлежности ц(х) каждого из параметров к терму «приемлемое значение» соответствующей лингвистической переменной. В рассматриваемом случае представляется естественным принять в качестве оценки переходного процесса

Р

Заданы следующие степени принадлежности:

Н )=

0,

tH < 200

1

tН -1.5, tН £

400 1, t Н > 600

(0, S<0 I

p.(S)=-!-0,01S-1, Se (0,100) > 11, S>100 I

R(t з )=

R(Df)=

0, t3<0

——13-1, t3 £

1350 3 3 1, t3 >1350

Df < 0

1 45

Df-1, Df e

1,

Df > 45

Исследуя работу реального фильтра частоты вращения, оценим переходный процесс при уста-

новлении величины частоты оборотов О, равной 600 об/сек.

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

Нагрузка Двигатель

Исполнительный ^ ' я механизм

Цифровой регулятор положения рейки

Цифровой регулятор

Настройка

Рис. 4

Необходимо автоматизировать поиск оптимальных настроек ПИД-регулятора.

Схема поиска оптимальной кривой переходного процесса

Для поиска оптимальных настроек, на которых оценочная функция (8) достигает максимума, необходимо уметь оценивать переходный процесс

для определенных численных настроек (Кр,К, ,

К°°) ПИД-регулятора. Схема поиска оптимальной кривой переходного процесса выглядит следующим образом [4]:

1) эксперт задает функции принадлежности терма «приемлемое значение» для параметров Ьн, 8, Ьз, Б,;

2) устанавливаются начальные величины коэффициентов регулятора (Кр, К?, К°);

3)вводятся случайные значения поправок к начальным величинам коэффициентов (ЛКР, ДК,, ЛК0);

4) вычисляются новые значения коэффициентов КР = Кр +ДКР , К, = К? +ЛК,, К°° = КО + +ЛК°° ;

5) производится численное моделирование системы регулирования частоты вращения с коэффициентами, полученными в п. 4, при этом получается переходная характеристика, отображенная на рисунке 2;

6) вычисляются значения параметров Ьн, 8, ЬЗ, Б,;

7) вычисляется степень принадлежности ^(х) каждого из параметров к нечеткому множеству «приемлемое значение» соответствующей лингвистической переменной;

8) вычисляется общая оценка переходной кривой по выражению (8);

9) если вычисленная оценка в п. 8 является максимальной, то алгоритм останавливается, иначе повторяется процедура для новых значений

(Кр,К?,К°0).

Генетический подход

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

1) на протяжении к-й итерации сохраняем популяцию потенциальных решений

Р(к)={хк ,...,хП}, в начальный момент времени

можно наполнить 1-ю популяцию любыми при-

к

годными значениями; отдельной хромосомой х, является набор коэффициентов ПИД-регулятора (Кр,К0,К°);

2) каждое решение хк оцениваем мерой пригодности Г(хк) по выражению (8);

3) формируем к+1 поколение методом селекции решений в соответствии с мерой пригодности

- каждое решение хк отбирается в следующее

Г,

поколение с учетом вероятности р, =

4) некоторые решения из следующего поколения подвергаем действию генетических операторов (кроссовера и мутации) для образования новых решений [7];

5) алгоритм завершает работу при достижении заданных значений пригодности либо по исчерпании числа поколений.

При исследованиях используем модификацию ГА, которая сохраняет наилучшее решение [8].

Результат применения генетического подхода

В результате эксперимента при генетической настройке регулятора без применения фильтра частоты внутри контура управления и без учета флуктуационной составляющей были получены следующие параметры пропорционально-дифференциального регулятора: КР=7,74, К1=0, Ко=0,14. Неравномерность сигнала, зависящая от цикловой подачи топлива (так называемые горбы с периодом 88 тактов), «срезана». Это отрицательно влияет на работу двигателя - происходит более быстрый его износ.

В результате эксперимента при генетической настройке регулятора с применением фильтра

'

о

о

п

частоты внутри контура управления и без учета флуктуационной составляющей были получены следующие параметры пропорционально-дифференциального регулятора: Кр=5,4, К^0, К0=0. Данные о неравномерности сигнала, зависящие от цикловой подачи топлива, «не срезаны». Это положительно влияет на работу двигателя - более оптимально регулируется рабочая частота оборотов коленчатого вала.

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

Получили ПИД-регулятор с параметрами Кр=10,01, К1=2,81, Ко=0,12. Недостатки те же, что и без флуктуационной составляющей.

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

Получен ПИД-регулятор с параметрами Кр=14,12, К1=1,33, К0=0,07.

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

Рис. 6

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

Литература

1. Хрящев Ю.Е., Кирик В.В., Третьяков А.А Использование аппарата fuzzy-логики в управлении дизелем // Математические методы в технике и технологиях: сб. тр. XX междунар. науч. конф. Ярославль: Изд-во ЯГТУ, 2007. С. 310-312.

2. Калугин Ф.В. Двухрежимный алгоритм фильтрации частоты оборотов коленчатого вала дизельного двигателя // Мехатроника, Автоматизация, Управление. 2008. № 5. С. 35-41.

3. Крутов В.И. Автоматическое регулирование двигателей внутреннего сгорания. М.: Машиностроение, 1979. С. 616.

4. Денисенко В. ПИД-регуляторы: принципы построения и модификации // Современные технологии автоматизации. 2006. № 4. С. 66-74.

5. Giuseppe Munda. Multicriteria evaluation in a fuzzy environment. Physica-Verlag, 1995.

6. Chuck Karr. Genetic Algorithms for Fuzzy Controllers, AI Expert, 1991, Febriary, pp. 25-32.

7. Kim Chwee Ng, Yun Li. Design of Sophisticated Fuzzy Logic Controller Using Genetic Algorithms Proc. 3rd IEEE Int. Conf. On Fuzzy Systems, Orlando, FL, June, 1994. Vol. 3, pp. 1708-1712.

8. Goldberg D.E. Genetic Algorithms in Search, Optimization, and Machine Learning. Addison-Wesley, 1989.

ПАРАЛЛЕЛЬНАЯ РЕАЛИЗАЦИЯ МОДЕЛИ MAP-REDUCE С ИСПОЛЬЗОВАНИЕМ ШАБЛОННЫХ КЛАССОВ С++

А.А. Московский, к.х.н.; А.Ю. Первин; Е.О. Тютляева

В работе рассматривается подход к реализации модели параллельного программирования Map-Reduce с помощью шаблонов параллельного программирования, реализованных как шаблонные классы С++. Шаблоны обеспечивают высокоуровневую конструкцию параллельных программ, используя в качестве низкоуровневого средства параллельного программирования библиотеку T-Sim. Программисту предоставляется возможность настройки поведения шаблонов в зависимости от специфики алгоритмов обработки данных и объема входных данных с помощью классов стратегий, что влияет на производительность параллельного приложения. Использование шаблона позволяет получать хорошо масштабируемые параллельные приложения для самых различных прикладных областей.

Ключевые слова: шаблоны параллельного программирования, масштабируемые параллельные приложения, Map-Reduce.

Задача создания эффективных параллельных программ не только не потеряла своей актуальности с ростом числа вычислительных узлов в высокопроизводительных установках, а, наоборот,

приобрела еще большую остроту. Простое увеличение числа вычислителей не гарантирует рост производительности системы. Согласно закону Амдала, максимальный прирост производительно-