ОСНОВЫ МОДЕЛИРОВАНИЯ СИСТЕМ …bek.sibadi.org/fulltext/ED578.pdf · 8 – Simulink - открывает окно браузера библиотек Simulink;

Министерство образования Российской Федерации Сибирская государственная автомобильно-дорожная академия

( СибАДИ )

В.С. Щербаков, А.А. Руппель, В.А. Глушец

ОСНОВЫ МОДЕЛИРОВАНИЯ СИСТЕМ АВТОМАТИЧЕСКОГО

РЕГУЛИРОВАНИЯ И ЭЛЕКТРОТЕХНИЧЕСКИХ

СИСТЕМ В СРЕДЕ MATLAB И SIMULINK

Учебное пособие

Омск Издательство СибАДИ

2003

2

УДК 681.06:621.3 ББК 32.973.26-018.2 Щ 61

Рецензенты: доктор технических наук, профессор С.Ф. Абдулин (ОмГИС), доктор технических наук, профессор И.Г. Браилов (ОмГТУ)

Работа одобрена редакционно-издательским советом академии и учебно-методическим советом факультета ТТМ в качестве учебного пособия по дисциплинам “Теория автоматического управления”, “Автоматизация технологических процессов и производств”, “Проектирование систем управления” для специальности 210200 – Ав-томатизация технологических процессов и производств.

Щербаков В.С., Руппель А.А., Глушец В.А. ОСНОВЫ МОДЕЛИРОВАНИЯ

СИСТЕМ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ И ЭЛЕКТРОТЕХНИЧЕСКИХ СИСТЕМ В СРЕДЕ MATLAB И SIMULINK: Учебное пособие. – Омск: Изд-во СибАДИ, 2003. – 160 с.

Предназначено для студентов всех форм обучения по специальностям, изучаю-

щим основы теории автоматического управления, теоретические основы электротехни-ки и электропривода.

Учебное пособие в простой и доступной форме знакомит читателей с про-граммными продуктами MАTLAB и Simulink. Приводятся типовые структурные схемы систем автоматического управления и их характеристики. Рассмотрены задачи анализа и синтеза как линейных, так и нелинейных систем автоматического управления. При-ведены примеры моделирования электротехнических схем. Изложение сопровождается подробно разработанными примерами, иллюстрациями, пиктограммами.

Ил. 165. Библиогр.: 20 назв.

ISBN 5 – 93204 – 154 – 4 © В.С. Щербаков, А.А. Руппель, В.А. Глушец, 2003 © Издательство СибАДИ, 2003

3

Оглавление

1. Знакомство с программным продуктом MATLAB......................

1.1.Основные сведения о системе MATLAB...................................... 1.2. Пользовательский интерфейс MATLAB...................................... 1.2.1. Командное меню.............................….................................. 1.2.2. Панель инструментов.......................................……............ 1.2.3. Панель запуска компонентов системы и браузер рабочего пространства................................................…..... 1.2.4. Окно истории сессии и браузер текущей директории…..

2. Пакет моделирования динамических систем Simulink............... 2.1.Назначение и функции пакета Simulink........................................ 2.1.1. Командное меню браузера библиотек...........................….. 2.1.2. Панель инструментов браузера библиотек Simulink......... 2.2. Блоки библиотеки Simulink........................................................... 2.2.1. Continuous – раздел непрерывных блоков....................….. 2.2.2. Discrete – раздел дискретных блоков............................….. 2.2.3. Function & Tables – раздел блоков функций и таблиц..… 2.2.4. Math – раздел математических блоков...........................… 2.2.5. Nonlinear – раздел нелинейных блоков........................….. 2.2.6. Signal & Systems – раздел блоков сигналов и систем...… 2.2.7. Sinks – раздел блоков получателей сигналов................…. 2.2.8. Sources – раздел блоков источников сигналов.................. 2.2.9. Библиотека Simulink Extras..............................................… 2.2.9.1. Additional Discrete – дополнительные дискретые блоки..……………………………………………… 2.2.9.2. Additional Linear – дополнительные линейные блоки........………………………………………….. 2.2.9.3. Дополнительные блоки Additional Sinks....……… 2.2.9.4. Блоки триггеров Flip Flops............................…….. 2.2.9.5. Раздел преобразований Transformations..........…... 2.3. Создание модели устройства в программе Simulink................... 2.3.1. Создание нового окна модели.......................................….. 2.3.2. Панель инструментов.....................................................….. 2.3.3. Командное меню.............................................................….. 2.3.4. Размещение блоков в окне модели................................…. 2.3.5. Выделение блока или группы блоков в окне модели..…. 2.3.6. Сохранение модели.........................................................…. 2.3.7. Установка параметров блоков, входящих в модель.....…. 2.3.8. Установка параметров моделирования..........................…. 2.3.9. Добавление надписей и текстовых комментариев………

6 6 7 8 9 10 11 12 12 14 15 16 16 16 17 17 19 20 22 22 23 24 24 25 26 27 27 27 28 29 30 31 31 31 31 33

4

2.3.10. Соединение блоков между собой..........................……... 2.3.11. Создание отвода линий...........................................…….. 2.3.12. Удаление соединений..............................................…….. 2.3.13. Изменение размеров блоков...................................…….. 2.3.14. Перемещение блоков и вставка блоков в соединение... 2.3.15. Печать текущей модели...........................................……. 2.3.16. Вставка модели в текстовые редакторы………………. 2.3.17. Создание простейшей модели..........................................

3. Вид переходных процессов и частотных характеристик типо-вых динамических звеньев...............................................................

3.1. Переходные характеристики звеньев.......................................... 3.2. Частотные характеристики звеньев............................................. 3.3. Получение переходных и частотных характеристик.................

3.4. Идеальное усилительное (безынерционное) звено.................... 3.5. Апериодическое (инерционное) звено........................................

3.6. Колебательное звено.........................................................……… 3.7. Интегрирующее звено.....................................................………. 3.8. Дифференцирующее звено.............................................………. 3.9. Звено с постоянным запаздыванием..............................………. 4. Способы соединения звеньев в структурную схему.....................

4.1. Последовательное соединение звеньев........................................ 4.2. Параллельное соединение звеньев............................................... 4.3. Звено, охваченное обратной связью............................................. 5. Решение задач анализа линейных систем автоматического управления........................................................................................... 5.1. Понятие устойчивости и критерии устойчивости систем ав-

томатического управления............................................................... 5.1.1. Критерий устойчивости по корням характеристического уравнения.........................................................................…. 5.1.2. Частотный критерий устойчивости Найквиста...... 5.2. Определение запасов устойчивости систем автоматического

управления............................................................................................. 5.3. Показатели качества переходных процессов...............................

5.4. Общая методика анализа линейных систем автоматического управления..............................................................................................

6. Решение задач синтеза линейных систем автоматического управления, коррекция их динамических свойств......................

6.1. Назначение и виды коррекции динамических свойств систем автоматического управления............................................................... 6.2. Получение передаточных функций регуляторов в пакете Simulink................................................................................................... 6.3. Общая методика синтеза линейных систем автоматического

34 34 35 35 35 35 36 36 41 41 42 43 51 55 57 58 60 62 65 65 65 66 67 67 69 70 73 75 80 87 87 90

5

управления.............................................................................................. 6.4. Приложение для синтеза корректирующих звеньев SISO De-sign Tool...........................................................................................…... 6.4.1. Командное меню................................................................... 6.4.2. Панель инструментов........................................................... 6.4.3. Работа с приложением SISO Design Tool........................... 6.4.4. Синтез по управляющему воздействию............................. 6.4.5. Синтез по возмущающему воздействию............................

7. Анализ простейших структур при случайных возмущающих воздействиях.........................................................................................

7.1. Понятие случайной величины и случайного процесса. Законы распределения случайных величин......................................................

7.2. Прохождение случайного сигнала через линейную систему.... 8. Нелинейные блоки пакета Simulink................................................

8.1. Блок ограничения Saturation......................................................... 8.2. Блок с зоной нечувствительности Dead Zone.............................. 8.3. Релейный блок Relay...................................................................... 8.4. Блок с ограничением скорости Rate Limiter................................ 8.5. Блок следящего квантования Quantizer........................................ 8.6. Блок фрикционных эффектов Coulombic and Viscous Friction.. 8.7. Блок люфта Backlash...................................................................... 8.8. Блоки-переключатели Switch........................................................

9. Пакет Nonlinear Control Design (NCD) Blockset.............................. 9.1. Основные правила работы с NCD.................................................

9.2. Оптимизация коэффициента передачи И-регулятора................. 9.3. Меню окна блока NCD Output....................................................... 9.4. Особенности решаемых задач.......................................................

10. Моделирование электротехнических систем................................ 10.1. Состав библиотеки Power System Blockset................................. 10.1.1. Состав библиотеки Electrical Sources.............................. 10.1.2. Соединительные элементы Connction............................. 10.1.3. Библиотека компонентов Elements.................................. 10.1.4. Состав библиотеки Power Electronics………………….. 10.1.5. Блоки измерений Measurements....................................... 10.1.6. Библиотека электрических машин Machines................. 10.1.7. Расширенная библиотека Powerlib Extras...................... 10.2. Пример расчета электрических цепей постоянного тока......... 10.3. Пример расчета однофазных электрических цепей синусои-дального тока..........................................................................………... 10.4. Пример моделирования асинхронного электрического дви-гателя с короткозамкнутым ротором...........................................……

Библиографический список.......................................................................

92 97 98 105 106 107 112 117 117 119 125 126 127 128 130 130 131 132 133 135 135 137 140 141 143 144 144 145 146 147 148 149 151 153 155 158 160

6

1. ЗНАКОМСТВО С ПРОГРАММНЫМ ПРОДУКТОМ MATLAB

1.1. Основные сведения о системе MATLAB Использование персональных компьютеров в учебном процессе позво-

ляет облегчить труд преподавателя и дать практические навыки в решении тех или иных задач, возникающих в процессе обучения. Для этих целей используется целый набор интегрированных программных систем и паке-тов программ, охватывающих практически весь спектр решаемых задач. В частности, средства для автоматизации математических расчетов: Reduce, Derive, TK Solver, Mathcad, Mathematica, Maple V, MATLAB и многие другие.

Система автоматизации математических расчетов MATLAB является одной из старейших и мощнейших в этой области. Название системы про-изошло от сокращения – МАТричная ЛАБоратория (MATrix LABoratory), но эта ориентация программы практически не ощущается теми пользова-телями, которые не интересуются непосредственно матричными вычисле-ниями, вследствие проработанности и продуманности синтаксиса языка программирования системы.

Система MATLAB была разработана фирмой Math Works, Inc. как язык программирования высокого уровня для математических, инженерных, на-учных и технических вычислений. В процессе развития система “обраста-ла” различными дополнениями в виде пакетов расширения, с увеличением количества которых росли и возможности системы. Но также несколько увеличилась сложность ее освоения на первых этапах, что несколько ком-пенсируется развитой системой помощи (на английском языке).

Как и всякая компьютерная программа, система MATLAB обладает как достоинствами, так и недостатками. Достоинства системы очевидны и час-тично описаны выше, лишь необходимо отметить легкость написания про-грамм в данной системе и открытость ее архитектуры, что позволяет легко изменять систему “под себя” и создавать новые пакеты расширений. К не-достаткам системы можно отнести: линейную запись выражений и фор-мул, что по сравнению, например, с системой Mathcad выглядит несколько неудобочитаемым; некоторую громоздкость системы, что создает сложно-сти в ее изучении на первых этапах.

В данный момент существует довольно большое количество литерату-ры, посвященной системе MATLAB и ее приложениям, но даже вся эта литература не может охватить все возможности и способности этого про-граммного продукта. В частности, недостаточно освещен вопрос о приме-нении системы MATLAB (а точнее, пакета расширения Simulink) для ре-шения задач анализа и синтеза систем автоматического управления (САУ). Поскольку данное учебное пособие посвящено исключительно этому во-

7

просу, в нем будет использоваться только пакет расширения Simulink и ра-бота непосредственно с системой MATLAB описываться не будет. Но так как пакет расширения Simulink не может функционировать без системы MATLAB, в следующем подразделе даны краткие рекомендации по поль-зованию системой MATLAB.

Также необходимо отметить, что данное учебное пособие ориентиро-вано на использование версии MATLAB 6.0 Release 12 и версии Simulink 4.0 (далее в тексте - без номера версии).

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

Windows 95 OSR2, Windows 98 или Windows 98 OSR2. Запуск системы MATLAB 6.0 производится из главного меню операционной системы либо двойным щелчком левой кнопки мыши на иконке MATLAB на рабочем столе.

После запуска появляется основное окно системы MATLAB. Это обычное окно приложений Windows, его можно перемещать, изменять его размер, открывать на весь экран и т. д.

1.2. Пользовательский интерфейс MATLAB

Интерфейс MATLAB 6.0 использует все возможности операционной системы Windows, следовательно, все окна, появляющиеся во время рабо-ты системы, можно перемещать, изменять их размеры, открывать на весь экран, сворачивать и закрывать.

После запуска системы появляется основное окно MATLAB (рис.1.1), которое состоит из пяти частей:

1 – командное меню; 2 – панель инструментов; 3 – панель запуска компонентов системы (Launch Pad) и браузер (редак-

тор) рабочего пространства (Workspace); 4 – окно истории сессии1 (Command History) и браузер текущей директо-

рии (Current directory); 5 – командное окно (Command Window).

1 Сессией здесь называется сеанс работы в системе MATLAB в прямом режиме ввода команд

8

Рис.1.1. Основное окно системы MATLAB

1.2.1. Командное меню

Командное меню (рис.1.2) предназначено для управления и настройки системы и содержит следующие пункты: File – операции открытия нового файла, загрузки, сохранения и печати файлов, стандартные для операционной системы Windows; Edit – операции редактирования сессии и работы с буфером обмена (команды Cut, Copy и Paste), стандартные для операционной системы Windows; View – управление внешним видом окон MATLAB; Web – доступ к интернет-ресурсам; Window – переключение между открытыми окнами MATLAB; Help – доступ к справочной системе и демонстрационным примерам MATLAB.

Рис.1.2. Командное меню системы MATLAB

1

2

3

4

5

9

Стандартный способ открытия меню мышью - это выбор пункта меню и щелчок по нему левой кнопки мыши. Кроме того, для открытия какой-либо позиции главного меню можно воспользоваться “горячими клавиша-ми” – при нажатой клавише Alt нажать клавишу символа, подчеркнутого в названии позиции. Например, чтобы открыть позицию Window, надо, удерживая клавишу Alt, нажать клавишу W.

Рис.1.3. Подменю командного меню сис-

темы MATLAB

Открытая позиция главного меню выглядит как подменю (рис.1.3), со-

держащее различные операции и команды. Выделенная команда или опе-рация исполняется при нажатии клавиши Enter. Выполнение команды можно также осуществить двойным щелчком мыши.

1.2.2. Панель инструментов

Панель инструментов (рис.1.4) системы MATLAB предназначена для быстрого доступа к основным функциям интерфейса, что осуществляется с помощью кнопок, выполняющих следующие функции: 1 – New M-file - выводит пустое окно редактора m-файлов; 2 – Open file - открывает окно для загрузки m-файла; 3 – Cut - вырезает выделенный фрагмент и помещает его в буфер обмена; 4 – Сору - копирует выделенный фрагмент в буфер; 5 – Paste - вставляет фрагмент из буфера в текущую строку ввода; 6 – Undo - отменяет предшествующую операцию; 7 – Redo - восстанавливает последнюю отмененную операцию; 8 – Simulink - открывает окно браузера библиотек Simulink; 9 – Help - открывает окно справки; 10 – Current Directory - выбор рабочего каталога системы.

О назначении кнопок говорят и всплывающие подсказки, появляющие-ся, когда стрелка мыши устанавливается и некоторое время удерживается на соответствующей кнопке.

Рис.1.4. Панель инструментов системы MATLAB

10

Наиболее важной кнопкой (для целей данного учебного пособия) явля-ется кнопка открытия окна браузера библиотек Simulink (см. рис.1.4, поз.8).

1.2.3. Панель запуска компонентов системы и браузер рабочего пространства

Панель запуска (рис.1.5) предназначена для быстрого доступа к компо-

нентам системы MATLAB (в том числе и к Simulink) и к их описаниям, а также к ссылкам Internet.

Рис.1.5. Панель запуска компонентов системы MATLAB

Браузер рабочего пространства (рис.1.6) предназначен для удобного просмотра хранящихся в оперативной памяти переменных, используемых системой при своей работе (пакет расширения Simulink также создает пе-ременные в рабочей области).

Рис.1.6. Браузер рабочего пространства

системы MATLAB

11

1.2.4. Окно истории сессии и браузер текущей директории

В окне истории сессии (рис.1.7) отображаются все команды, выполнен-ные при предыдущих запусках системы, которые хранятся до перезагрузки компьютера. При необходимости сессию можно сохранить.

Рис.1.7. Окно истории сессии системы MATLAB

Рис.1.8. Браузер текущей директории системы MATLAB

Браузер текущей директории (рис.1.8) отображает содержимое текущей директории (рабочего каталога) системы (см. рис.1.4, поз.10), а также ука-зывает тип файла и время его последнего изменения.

12

2. ПАКЕТ МОДЕЛИРОВАНИЯ ДИНАМИЧЕСКИХ СИСТЕМ SIMULINK

2.1. Назначение и функции пакета Simulink

Версия системы MATLAB 6.0 поставляется вместе с пакетом расшире-

ния Simulink 4.0, предназначенным для моделирования динамических сис-тем, модели которых составляются из отдельных блоков (компонентов).

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

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

Пакет Simulink является ядром интерактивного программного ком-плекса, предназначенного для математического моделирования линейных и нелинейных динамических систем и устройств, представленных своей функциональной блок-схемой, именуемой S-моделью, или просто моде-лью. При этом возможны различные варианты моделирования: во времен-ной области, в частотной области, с событийным управлением, на основе спектральных преобразований Фурье, с использованием метода Монте-Карло (реакция на воздействия случайного характера) и т. д.

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

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

13

Для запуска пакета Simulink необходимо нажать кнопку Simulink в па-нели инструментов главного окна системы MATLAB (см. рис.1.4, поз.8). При нажатии этой кнопки открывается окно браузера библиотек (рис.2.1).

Рис.2.1. Окно браузера библиотек пакета Simulink

В окне браузера содержится дерево компонентов библиотек Simulink. Для просмотра того или иного раздела библиотеки достаточно выделить его нажатием левой кнопки мыши, и в правой части окна Simulink Browser Library появится набор пиктограмм компонентов активного раздела биб-лиотеки.

Окно браузера библиотек состоит из пяти основных частей: 1 – командное меню; 2 – панель инструментов; 3 – поле информационных сообщений; 4 – окно разделов библиотеки;

1

2 3

4 5

6

14

5 – окно с входящими в эти разделы подразделами или блоками; 6 – строка состояния.

2.1.1. Командное меню браузера библиотек Командное меню содержит следующие пункты и подпункты:

1) File - работа с файлами библиотек: New – создание новой модели или нового раздела библиотеки, куда мо-гут входить как новые, так и уже существующие блоки; Open – открытие файла (окна) модели или нового раздела библиотеки; Preferences – редактирование свойств пакета Simulink. 2) Edit - добавление компонентов и их поиск: Add to the current model – добавить выделенный блок в текущее окно модели; Find block – найти блок в библиотеке пакета Simulink; Find next block – продолжить поиск (по критерию из предыдущего пунк-та меню). 3) View - управление элементами браузера библиотек Simulink: Toolbar – показать/убрать панель инструментов; Status bar – показать/убрать строку состояния; Description - показать/убрать поле информационных сообщений; Stay on top – включить/выключить режим браузера "поверх всех окон"; Collapse entire Browser – свернуть все разделы библиотек в окне разде-лов браузера библиотек; Expand entire Browser – развернуть все разделы библиотек в окне разде-лов браузера библиотек; Large/Small Icon – устанавливает размер иконок (большой/маленький) в окне, показывающем содержание разделов библиотек; Show parameters for selected block – показать параметры для выбранного блока (эквивалентно двукратному нажатию на нем левой кнопкой мышки). 4) Help - вывод справки по Simulink: Help on the selected block – справка по выбранному блоку; Simulink Help – справка по системе Simulink; Tip of the Day – выводит следующую подсказку в поле информационных сообщений.

15

2.1.2. Панель инструментов браузера библиотек Simulink

Панель инструментов (рис.2.2) содержит следующие кнопки: 1 – Create a new model - открывает новое окно для создания модели; 2 – Open a model - открывает существующую модель; 3 – Stay on top - установка статуса окна браузера "поверх всех окон"; 4 – Find block - найти блок.

Рис.2.2. Панель инструментов браузера библиотек Simulink

В поле информационных сообщений возникает название и краткое

описание выделенного раздела или блока. Если в окне разделов библиотек выделен какой-то раздел, то щелчок

правой кнопкой мыши выводит контекстное меню с одной командой, по-зволяющей вывести данный раздел в окне, расположенном справа от окна разделов, такой же эффект достигается двукратным нажатием левой кноп-ки мыши на нужном разделе.

На рис. 2.1 выделена основная библиотека Simulink (в левом окне) и показаны ее разделы (в правом окне). Основная библиотека Simulink со-держит достаточно большое количество самых разнообразных блоков, в рамках данного учебного пособия нужны лишь некоторые из них, приме-няемые в теории автоматического управления (ТАУ).

Кроме того, при нажатии правой кнопки мыши на конкретном блоке выводится небольшое контекстное меню, содержащее следующие коман-ды: 1) Add to the current model – добавить выделенный блок в текущую модель (также можно использовать сочетание клавиш Ctrl+I); 2) Help for the ‘block name’ block – справка по выделенному блоку (в ка-вычках содержится имя выделенного блока); 3) Go up a level – перейти в библиотеке на уровень выше; 4) Block parameters – параметры выделенного блока (достигается так же, двукратным нажатием левой кнопки мыши на блоке).

16

2.2. Блоки библиотеки Simulink

2.2.1. Continuous – раздел непрерывных блоков

Рис. 2.3. Раздел Continuous

В состав раздела непрерывных блоков входят (рис 2.3):

Derivative – дифференцирующий блок; Integrator – интегрирующий блок; Transport Delay – блок постоянного запаздывания; Transfer Fcn – блок, реализующий передаточную функцию, заданную в виде отношения полиномов (апериодическое и колебательное звенья); State-Space – линейная аналоговая система, заданная в виде уравнений состояния, т. е. в виде системы уравнений, представленной в форме Коши; Zero-Pole – линейная аналоговая система, заданная своими нулями и по-люсами; Memore – блок памяти, выполняющий задержку на один шаг модельного времени; Variable Transport Delay – блок памяти с переменной задержкой.

2.2.2. Discrete – раздел дискретных блоков

Discrete Transfer Fcn - блок, реализующий передаточную функцию, за-данную в виде отношения полиномов (апериодическое и колебательное звенья), для дискретных САУ; Discrete Time Integrator – интегрирующий блок для дискретных САУ. Выполняет численное интегрирование входного сигнала; Zero-Order-Hold –экстраполятор нулевого порядка; Unit Delay – блок задержки сигнала. Обеспечивает задержку входного сигнала на заданное число шагов модельного времени; Discrete State-Spase – блок задания дискретного звена матрицами его со-стояния;

17

Рис. 2.4. Раздел Discrete

Discrete Filter – блок задания дискретного звена через дискретную пере-даточную дробно-рациональную функцию относительно 1/Z; Discrete Zero-Pole – блок задания дискретного звена через указание зна-чений нулей и полюсов дискретной передаточной функции относительно 1/Z; First-Order Hold – экстраполятор первого порядка.

2.2.3. Function & Tables – раздел блоков функций и таблиц

Fcn – блок, реализующий в пакете Simulink любую функцию системы MATLAB, где в качестве параметра настройки можно ввести любое вы-числяемое выражение, аргументом которого выступает входной сигнал; MATLAB Fcn – подобен предыдущему блоку, но также позволяет реали-зовывать m-функции (программы, написанные на языке программирования MATLAB); Look-Up Table, (2-D) и (n-D) – блоки данных, заданных в виде таблиц (позволяют задать нелинейные блоки любого вида), отличаются друг от друга размерностью задающих таблиц. Все эти блоки выполняют различ-ного вида интерполяции.

2.2.4. Math – раздел математических блоков

Gain – усилительный блок; Sum – сумматор сигналов с различными знаками и любым количеством входов;

18

Product – вычислитель, формирующий на выходе результат умножения или деления двух и более входных сигналов. В качестве параметров на-стройки указывается число входов и вид выполняемой операции; Dot Product – звено, осуществляющее перемножение двух входных ве-личин, если они являются скалярами. Это звено вычисляет также сумму поэлементных произведений двух входных векторов одинаковой длины; Slider Gain – аналоговый усилитель с интерактивной настройкой; Matrix Gain – усилитель, на вход которого подается вектор; Mat Function – блок, позволяющий выбрать одну из математических функций в поле настройки и включить ее в модель; Trigonometric Function – звено формирования тригонометрических функций от входного сигнала. Выбор функции обеспечивается в поле на-стройки;

Рис. 2.5. Раздел Math

MinMax – блок выбирает минимальное или максимальное значение век-тора в соответствии с заданием поля настройки. Входной сигнал на блок задается числовым вектором. В окне настройки определяется также коли-чество входов;

19

Abs – блок, формирующий на выходе абсолютное значение входного сигнала; Sign – блок-реле, реагирующий на знак входного сигнала; Rounding Function – округление входного сигнала; Combinatorial Logic – блок обеспечивает преобразование входного сиг-нала в соответствии с формированной в окне настройки таблицей истинно-сти; Logical Operation? Relaition Operator – блоки производят известные ло-гические операции «и» и «или», количество входов задается в поле на-стройки; Bitwise Logical Operator – универсальный блок, реализующий любую ло-гическую функцию; Complex to Magnitude-Angle – блок, позволяющий выделить модуль и фазу входной комплексной величины; Magnitude-Angle to Complex – блок, преобразующий входную величину, заданную модулем и фазой в комплексную выходную величину; Complex to Real-Image, Real-Image to Complex – блоки, преобразующие комплексные величины из показательной формы в алгебраическую и об-ратно; Algebraic Constrain – блок, позволяющий в структурную модель вклю-чать систему алгебраических уравнений.

2.2.5. Nonlinear – раздел нелинейных блоков

Рис. 2.6. Раздел Nonlinear

Saturation – блок ограничения; Dead Zone – блок с зоной нечувствительности; Relay – релейный блок;

20

Rate Limiter – блок с ограничением скорости; Coulomb and Viscous Friction – блок фрикционных эффектов; Backlash – блок люфта; Switch и Multiport Switch – переключатели одно- и многовходовые; Quantizer – блок, обеспечивающий квантование входного сигнала по уровню. В системах управления такие блоки являются частью аналого-цифрового преобразователя; Manual Switch – блок, который переключается вручную. В процессе мо-делирования при помощи этого ключа удобно менять параметры и струк-туру модели; Coulumb & Viscous Friction – блок, реализующий характеристику трения в механических системах.

2.2.6. Signal & Systems – раздел блоков сигналов и систем

Sub System – блок-подсистема (полезен в случае громоздких систем, части которых можно с помощью данного блока заменить на один или не-сколько блоков-подсистем); In, Out – блоки входа/выхода, предназначенные для создания блоков-подсистем; Mux – микшер сигналов, собирающий несколько сигналов в одну шину; Demux – блок, выполняющий действие, обратное действию предыдуще-го блока; Enable, Trigger – блоки, предназначенные для логического управления работой модели; Bus Selector – блок выделяет из присоединенной к его порту шины тре-буемые сигналы. В окне настройки блока имеется два списка – входной и выходной; Selector – блок выбирает из входного вектора элементы, которые указа-ны в параметрах настройки; Merge – блок, осуществляющий объединение входных сигналов; Matrix Concatenation – блок, позволяющий векторный сигнал, представ-ленный развернутой строкой или столбцом, преобразовать к «свернутому» векторному сигналу; From, Goto Tag Visibility, Goto – блоки ( «Принять», «Признак видимо-сти», «Передать» ) используются совместно и предназначены для обмена между различными данными модели с учетом их доступности; Data Store Memory, Data Store Read, Data Store Write – «Память», «Чте-ние» и «Запись» данных, также используются совместно и обеспечивают хранение и передачу данных; Ground, Terminator – блоки используются в качестве «заглушек» для не-

21

использованных входных и выходных портов соответственно;

Рис. 2.7. Раздел Signal & Systems

Reshape – блок, позволяющий изменить размерность входного сигнала; Data Type Conversion – блок, обеспечивающий приведение типа данных входного сигнала к требуемому; Function-Call Generator – блок, обеспечивающий запуск подключенных к нему подсистем с заданной периодичностью; Configurable Subsystem – блок, реализующий функцию любой подсисте-мы, которая может быть библиотечной или созданной пользователем; Model info – блок, позволяющий получить информацию о модели; IC – блок, позволяющий установить начальное значение входного сиг-нала, значение которого задается в окне настройки; Width – блок, вычисляющий размерность сигнала на входе; Probe – блок, позволяющий получить на выходе необходимую информа-

22

цию о входном сигнале. Блок имеет один вход, число выходов зависит от числа исследуемых параметров входного сигнала; Signal Specification – блок, управляющий процессом моделирования в зависимости от параметров входного сигнала.

2.2.7. Sinks – раздел блоков получателей сигналов

Display – блок, отображающий цифровую информацию (уровень сигнала в данный момент времени); Scope – виртуальный осциллограф, предназначенный для получения вре-менных зависимостей (переходных процессов); XY Graph – виртуальный графопостроитель (построение фазовых порт-ретов и различных двухмерных графиков);

Рис. 2.8. Раздел Sinks

To File и To Workspace – блоки, передающие информацию в файл или в рабочую область MATLAB; Stop Simulink – остановка симуляции.

2.2.8. Sources – раздел блоков источников сигналов

Band-Limited White Noice – генератор белого шума; Chrip Signal – генератор сигнала с нарастающей частотой; Constant - источник постоянного воздействия, задающий константу (t < 0, y = 0; t 0, y = сonst); Discrete Pulse Generator – источник дискретных импульсов; From File и From Workspace - источником сигнала для этих блоков слу-жит файл или рабочая область MATLAB; Ramp – источник нарастающего воздействия;

23

Рис. 2.9. Раздел Sources

Random Number – источник случайного сигнала с нормальным распре-делением; Repeating Sequence – источник пилообразного сигнала; Sine Wave – источник синусоидального воздействия; Step – источник перепада сигнала (t < Step time, y = Initial value; t > Step time, y = Final value); Signal Generator – универсальный сигнал-генератор; Clock – источник времени моделирования; Digital Clock – цифровой источник времени.

2.2.9. Библиотека Simulink Extras

Библиотека Simulink Extras является дополнительной библиотекой па-

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

Библиотека Simulink Extras представлена на рис. 2.10. Работа с этой библиотекой ничем не отличается от работы с основной библиотекой.

24

Рис.2.10. Библиотека Simulink Extras

2.2.9.1. Additional Discrete – дополнительные дискретные блоки

Рис. 2.11. Дополнительные дискретные блоки

Дополнительные блоки Additional Discrete представлены всего четырь-

мя блоками – по два варианта известных нам блоков Discrete Transfer Fcn и Discrete Zero-Pole (рис 2.11). Их единственным отличием от описанных ра-нее блоков является возможность инициализации входов и состояний.

2.2.9.2. Additional Linear – дополнительные линейные блоки

Состав дополнительных линейных блоков раздела Additional Linear

показан на рис. 2.12. Блоки этого раздела можно разделить на две катего-рии: PID-контроллеры и блоки типа State-Shfct, Transfer Fnc и Zero-Pole, дополненные возможностями выходных сигналов и состояний.

Для анализа и синтеза систем управления наибольший интерес пред-ставляют PID-контроллеры. Первый из них PID-controller – это довольно универсальный блок, выходной сигнал которого задается операторным выражением:

OUT = P + I / s + Ds, где Р – входной сигнал; I – его интеграл; D – его производная.

25

Рис. 2.12. Дополнительные линейные

блоки

Параметр Р (по умолчанию 1) фактически задает коэффициент усиле-ния безынерционного усилительного звена. Параметр I задает пропорцио-нальность интегралу входного сигнала и, наконец, параметр D задает про-порциональность производной входного сигнала. Параметры P, I и D за-даются в таблице настройки PID-контроллера, поэтому задавая различные значения указанных параметров, можно получить П -, ПИ -, ПД - и ПИД - регуляторы .

Второй PID-controller (with Approximate Derivative) c улучшенной опе-рацией дифференцирования вычисляет выходной сигнал как:

OUT = Р + I / s + Ds/(1/Ns + 1).

За счет применения дополнительного параметра N улучшается вычис-

ление производной.

2.2.9.3. Дополнительные блоки Additional Sinks

Раздел дополнительных блоков Additional Sinks содержит ряд новых виртуальных регистраторов ( рис. 2.13): Auto Correlator – автокоррелятор (используется с пакетом Signal Process-ing Toolbox); Averaging Power Spectral Density – анализатор спектральной плотности мощности с усреднением;

26

Рис. 2.13. Дополнительные блоки

Additional Sinks

Averaging Spectrum Analyzer – спектральный анализатор с усреднени-ем; Cross-Correlator – кросс-коррелятор; Floating Point Bar – утилита построения гистограммы; Power Spectral Density – анализатор спектральной плотности мощности; Spectrum Analyser – анализатор спектра.

Данные блоки относятся к двум важным разделам моделирования – статистическому анализу и анализу спектров сигналов ( по уровню и по мощности ).

2.2.9.4. Блоки триггеров Flip Flops

Рис. 2.14. Раздел библиотеки Flip Flops

Раздел библиотеки Simulink Extras Flip Flops содержит следующие ос-

новные блоки (рис. 2.14): Clock – генератор тактовых импульсов; D Flip Flops, D Latch, J-K Flip Flops, R-S Flip Flops – 4 триггерных уст-ройства типа D -, J-K, R-S – триггеров.

27

2.2.9.5. Раздел преобразований Transformations Раздел преобразований Transformations (рис. 2.15) содержит 8 блоков

для осуществления типичных преобразований – температуры из градусов Цельсия в градусы Фаренгейта и наоборот; углов, выраженных в градусах, в углы, выраженные в радианах, и наоборот; прямоугольных систем коор-динат в полярные системы координат и наоборот.

Рис. 2.15. Блок Transformations Основные настройки вышеописанных блоков и правила работы с ними

описаны ниже. Любая модель, собираемая в программном пакете Simulink, должна со-

стоять из трех основных частей: источник сигнала, модель, приемник сиг-нала. Причем модель может состоять из любых блоков, описанных выше (кроме разделов источников и приемников сигналов).

2.3. Создание модели устройства в программе Simulink

2.3.1. Создание нового окна модели

Для создания окна новой модели нужно нажать кнопку Create a new

model (см. рис.2.1, поз.1) в панели инструментов браузера библиотек Simulink либо выполнить соответствующую команду из командного меню. Пустое окно модели показано на рис. 2.16.

28

Рис.2.16. Пустое окно модели (окно для создания новой модели)

2.3.2. Панель инструментов

Окно модели имеет собственную панель инструментов, содержащую 15 кнопок (рис. 2.17).

Рис.2.17. Панель инструментов окна модели

Кнопки имеют следующее назначение: 1 – New model - создание новой модели (открытие нового окна модели); 2 – Open model – загрузка ранее сохраненной модели; 3 – Save - сохранение текущей модели; 4 – Print - печать текущей модели; 5 – Cut - перенос выделенного объекта в буфер; 6 – Copy - копирование выделенного объекта в буфер; 7 – Paste - вставка объекта из буфера; 8 – Undo - отмена последней операции; 9 – Redo - восстановление последней отмененной операции; 10 – Library Browser - открытие браузера библиотек; 11 – Toggle model browser - открытие браузера модели (показывает нали-чие и состав подмоделей) в левой части окна модели; 12 – Go to parent system - переход в основную (родительскую систему), актив-но в случае работы в окне подмодели; 13 – Debug - переход в режим отладки модели; 14 – Start Simulation - запуск моделирования; 15 – Stop Simulation - остановка моделирования.

29

2.3.3. Командное меню

Несмотря на наличие панели инструментов, важная роль принадлежит командному меню окна модели Simulink. Оно дает более полный набор средств по управлению процессом моделирования. Меню (рис.2.18) со-держит следующие пункты и подпункты (приведены лишь нужные для ра-боты, с учетом целей данного учебного пособия):

Рис.2.18. Командное меню окна модели

1) File - операции с файлами S-моделей: New, Open и Close – стандартные операции создания, открытия и закры-тия файла (окна) модели. Следует отдельно отметить возможность созда-ния не только модели, но и раздела библиотеки, в который можно помес-тить как новые, так и уже существующие блоки (Ctrl+N, Ctrl+O и Ctrl+W соответственно); Save – сохранить модель (при первом сохранении запрашивается имя и место сохранения) (Ctrl+S); Save as… – сохранить модель под определенным именем и в определен-ном месте (можно использовать для полного копирования модели задани-ем другого имени для сохранения); Print… – распечатать содержимое окна модели (структурную схему мо-дели) (Ctrl+P); Print setup… – установки для печати; Exit MATLAB – выход из MATLAB, закрывает все открытые окна моде-лей и приложений системы, а также сам MATLAB (Ctrl+Q). 2) Edit - операции редактирования текущей модели: Can’t undo – отменить операцию отмены последнего действия (вместо слова ‘undo’ будет наименование последней операции) (Ctrl+Z); Can’t redo – операция, обратная предыдущей (Ctrl+Y); Cut, Copy и Paste – стандартные операции редактирования – вырезать в буфер, скопировать в буфер и вставить из буфера (работают с выделенны-ми объектами в окне моделирования) (Ctrl+X, Ctrl+C и Ctrl+V соответст-венно); Clear – очистить (стереть) выделенный объект (Delete); Select All – выделить все объекты в окне (удобно для копирования или очистки окна модели или подмодели) (Ctrl+A); Copy model to clipboard - поместить модель в буфер обмена ОС Windows (удобно в случае, когда необходимо поместить модель в какой-либо электронный документ);

30

Create subsystem – создать из выделенных в окне модели блоков подсис-тему (Ctrl+G). 3) View - управление видом окна модели: Go to parent - переход в предыдущую модель/подмодель; Toolbar – показать/убрать панель инструментов; Status bar – показать/убрать строку состояния; Show Library Browser – показать браузер библиотек; Zoom In/Out - увеличение/уменьшение размера блоков в окне модели (изменение масштаба изображения блоков); Fit system to view – выбрать оптимальный масштаб изображения; Normal (100%) – вернуть “нормальный” (заданный по умолчанию) мас-штаб изображения. 4) Simulation - операции запуска моделирования и его настройки: Start – запуск процесса моделирования (Ctrl+T); Stop – остановка процесса моделирования; Simulation parameters – параметры процесса моделирования (Ctrl+E). 5) Format - операции форматирования текущей модели (работают с выде-ленными блоками и надписями): Font… – выбор шрифта текста; Text alignment – операции выравнивания текста в текстовом окне; Flip name – поместить название блока над ним или под ним; Hide name – скрыть название блока; Flip block – изменить направление блока (вход-выход) (Ctrl+I); Rotate block – повернуть блок на 90 градусов по часовой стрелке (Ctrl+R); Show/Hide drop shadow – показать/скрыть отбрасываемую блоком тень; Show/Hide port labels – показать/скрыть метки портов для подсистем; Foreground color – выбрать цвет рамки, опоясывающей блок; Background color – выбрать цвет поля блока; Screen color – выбрать цвет заднего фона окна модели; 6) Tools - доступ к инструментальным средствам (зависит от комплекта ус-тановки системы): Linear analysis… – пакет дополнения “линейный анализ”. 7) Help - доступ к средствам справочной системы.

2.3.4. Размещение блоков в окне модели

Размещение блоков в окне модели производится следующим образом: в

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

31

ши, перетаскиваем блок в окно модели и отпускаем. В окне модели должна появиться пиктограмма этого блока. Либо нажав правую клавишу мыши на нужном нам блоке в библиотеке, выбираем в контекстном меню команду Add to ‘имя модели’.

2.3.5. Выделение блока или группы блоков в окне модели

Для выделения блока достаточно навести на него стрелку мыши и на-

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

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

Для выделения всех блоков модели можно воспользоваться командой Edit > Select All.

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

2.3.6. Сохранение модели

Можно сохранить созданную модель для последующего применения,

показа или модернизации. Для этого используется команда Save или Save As… меню File окна редактора моделей. Модель записывается в виде фай-ла с расширением .mdl.

2.3.7. Установка параметров блоков, входящих в модель

Для того чтобы вызвать окно модификации параметров блока, нужно

навести курсор мыши на изображение компонента и дважды щелкнуть ле-вой кнопкой мыши либо, нажав правую кнопку мыши, выбрать в контек-стном меню команду Block parameters.

2.3.8. Установка параметров моделирования

Прежде чем запустить модель, стоит ознакомиться с установкой общих

параметров моделирования. Для этого необходимо выполнить команду Simulation Parameters... в меню Simulation окна Simulink. Появится окно ус-тановки параметров моделирования (рис.2.19).

32

Рис.2.19. Окно задания параметров моделирования Это окно имеет ряд вкладок с довольно большим числом параметров.

Но с учетом задач данного учебного пособия необходимо знание только одной вкладки, которая открывается по умолчанию Solver (решатель). Эта вкладка позволяет установить параметры решающего устройства системы моделирования Simulink.

К числу важнейших параметров решателя относится время моделиро-вания - Simulation time. Оно задается начальным временем Start time (обычно 0) и конечным временем Stop time. Равенство Stop time бесконеч-ности (inf) означает, что моделирование будет происходить бесконечно долго, пока мы не прервем его с помощью кнопки окна модели или коман-ды Stop из командного меню. Однако в этом случае трудно получить раз-личимые осциллограммы работы устройства, поэтому рекомендуется зада-вать конечные значения Stop time.

Далее следуют опции решателя (Solver options), задающие: Type – тип решателя и метод моделирования.

Тип решателя может быть: Variable step – с переменным шагом по вре-мени и Fixed step – с постоянным шагом по времени (для дискретных сис-тем).

Справа от типа решателя выбирается метод моделирования, возможен выбор следующих методов: discrete – дискретный, ode45 и ode5 – метод Дорманда-Принса, ode23 – три варианта, включая метод Розенброка и ode113 – метод Адамса, ode4 – метод Рунге-Кутта и др.

33

Следующая группа параметров изменяется от выбора типа решателя и для переменного шага моделирования содержит: Max step size – максимальное значение шага моделирования; Min step size – минимальное значение шага моделирования; Initial step size – начальный шаг моделирования.

Для трех вышеописанных параметров по умолчанию устанавливается значение Auto (выбрать автоматически), но это значение можно изменять и устанавливать требуемое для данной задачи. Relative tolerance – относительная погрешность моделирования; Absolute tolerance – абсолютная погрешность моделирования.

Для постоянного шага моделирования: Fixed step size – значение фиксированного шага моделирования (при значении auto устанавливается шаг, заданный источником сигнала, и если источников несколько – устанавливается наименьший шаг моделирова-ния); Mode – режим работы (auto – выбрать автоматически; Single tasking – однозадачный режим; Multi tasking – многозадачный режим).

При моделировании сложных систем необходимо правильно устанав-ливать значение всех вышеописанных параметров, но при моделировании линейных САУ непрерывного действия (не дискретных) важно только ус-тановить конечное время моделирования, остальные параметры можно ос-тавить установленными автоматически. Вышесказанное справедливо и для моделирования дискретных систем, необходимо лишь выбрать тип реша-теля – с постоянным шагом моделирования и указать шаг (при необходи-мости).

2.3.9. Добавление надписей и текстовых комментариев

В программировании существует мнение, что хороший программист –

тот, кто пишет необходимое количество комментариев к своей программе, что позволяет впоследствии без труда разобраться, что и как работает. Это справедливо и для моделей, созданных в пакете Simulink. Модели, не со-держащие текстовых комментариев, не наглядны и трудно воспринимают-ся другими людьми. Но и перенасыщение комментариями модели нежела-тельно.

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

Для изменения подписи к блокам моделей необходимо установить мышь в область надписи и щелкнуть левой кнопкой мыши - в подписи появится курсор ввода и ее можно будет редактировать.

34

Чтобы убрать надпись, нужно выделить ее (кстати, как и любой другой объект) и выполнить команду Edit > Clear или нажать клавишу Delete на клавиатуре.

В связи с тем, что система MATLAB является англоязычной (это отно-сится и к пакету Simulink), она плохо воспринимает русский язык. Поэтому при наличии большого числа комментариев на русском языке возможны сбои в работе системы.

2.3.10. Соединение блоков между собой

Соединение блоков между собой производится с помощью мыши. Бло-

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

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

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

Особо стоит отметить возможность задания наклонных линий соедине-ний при нажатой клавише Shift.

2.3.11. Создание отвода линий

Часто возникает необходимость сделать отвод от уже созданной линии.

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

При нажатой клавише Shift отвод строится наклонными линиями.

35

2.3.12. Удаление соединений

Для удаления соединительной линии достаточно выделить ее щелчком мыши и выполнить команду Clear или Cut (или нажать клавишу Delete на клавиатуре).

2.3.13. Изменение размеров блоков

Simulink имеет расширенные возможности редактирования блок-схем.

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

Растягивается только графическое изображение (пиктограмма) блока, а размеры его названия в виде текстовой надписи не изменяются.

2.3.14. Перемещение блоков и вставка блоков в соединение

Блок, участвующий в соединении, можно перемещать в окне модели,

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

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

2.3.15. Печать текущей модели

Команда File > Print... выводит окно печати модели. Это окно содержит

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

Команда Print Setup... в меню File окна Simulink открывает окно на-

36

стройки принтера. Это окно операционной системы Windows, поэтому ес-ли она русифицирована, то не стоит удивляться появлению в окне русскоя-зычных надписей. Параметры этого окна очевидны. Кнопка «Свойства» открывает окно свойств выбранного принтера.

2.3.16. Вставка модели в текстовые редакторы

Для вставки модели в различные текстовые или графические редакторы

необходимо выполнить команду Edit > Copy model to clipboard – модель скопируется в виде картинки в буфер обмена ОС Windows и может быть вставлена в любой текстовый или графический редактор путем нажатия клавиш Shift+Insert в открытом файле текстового редактора.

2.3.17. Создание простейшей модели

Любая модель, создаваемая в пакете Simulink, должна состоять из трех

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

единичного ступенчатого воздействия (блок Constant из раздела Sources) и виртуального осциллографа (блок Scope из раздела Sinks). Модель будет представлена проводником, который является идеальным усилительным звеном с коэффициентом усиления, равным 1.

Для создания модели и ее моделирования необходимо выполнить сле-дующие действия:

1) создать новое окно модели; 2) поместить в это окно блок Constant из раздела Sources (рис.2.20); 3) поместить в окно модели блок Scope из раздела Sinks и расположить

его справа от блока источника сигнала Constant (рис.2.21); 4) соединить блоки Constant и Scope между собой (рис.2.22); 5) установить необходимые параметры моделирования (Simulation

parameters) Stop time = 2 с (остальные параметры оставляем без измене-ния);

6) запустить процесс моделирования; 7) посмотреть результаты моделирования, дважды щелкнув левой

кнопкой мыши на блоке Scope (рис.2.23); 8) при необходимости сохранить модель.

37

Рис.2.20. Источник сигнала, помещенный в окно модели

Рис.2.21. Источник сигнала и приемник сигнала, помещенные в окно модели

Рис.2.22. Простейшая модель, готовая к процессу моделирования

38

Рис.2.23. Результаты моделирования В окне виртуального осциллографа (см. рис. 2.23) можно выделить две

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

Панель инструментов состоит из следующих кнопок (слева направо): 1 − Zoom – изменение масштаба графика переходного процесса по осям X-Y, с сохранением пропорций; 2 − Zoom X-axis – изменение масштаба графика переходного процесса по оси X; 3 − Zoom Y-axis – изменение масштаба графика переходного процесса по оси Y; 4 − Autoscale – автоматический подбор масштаба; 5 − Save current axes settings – сохранить текущие настройки блока; 6 − Properties – свойства блока;

39

7 − Print – печать графика переходного процесса. Окно свойств блока Scope (рис.2.24) содержит две вкладки: General –

основные свойства и Data history – сохранение данных. В первой вкладке находятся настройки:

1) Axes – осей: Number of axes – количество осей, при изменении этого параметра изме-няется количество одновременно отображаемых графиков (на разных эк-ранах, друг под другом) и соответственно количество входов; Time range – пределы временного интервала; Tick labels – настройка подписей по осям: bottom axis only – только под нижним графиком, all – под всеми графиками и none – не подписывать; 2) Sampling – установка временных соотношений: Decimation – в десятичных долях времени, по умолчанию 1; Sample time – в тактах эталонного времени (сэмплах), по умолчанию 0 – непрерывное время.

Рис.2.24. Окно свойств блока Scope, вкладка General

Во второй вкладке располагаются настройки хранения данных

(рис.2.25): Limit data point to last – включение/выключение предела и количества точек на графике (если моделируемая система работает достаточно долго и пропадает часть графика, то рекомендуется перед началом процесса моде-лирования убрать флажок с данной опции); Save data to workspace – сохранять результаты в рабочую область под именем, задаваемым в окне Variable name, и типы данных, задаваемых в окне Format (Structure with time – структура с временем, Structure – струк-тура и Array – массив).

40

Рис.2.25. Окно свойств блока Scope, вкладка Data history

Представленный на рис.2.23 результат моделирования − график пере-

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

Вместо блока Constant для получения единичного ступенчатого воздей-ствия возможно применение блока Step, но в последнем блоке необходимо изменить значение Step time с 1 на 0.

41

3. ВИДЫ ПЕРЕХОДНЫХ ПРОЦЕССОВ И ЧАСТОТНЫХ ХАРАКТЕРИСТИК ТИПОВЫХ ДИНАМИЧЕСКИХ ЗВЕНЬЕВ

Системы автоматического регулирования различного назначения име-

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

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

В теории автоматического регулирования можно выделить основные типовые динамические звенья: усилительное (безынерционное, пропор-циональное), апериодическое (инерционное), дифференцирующее, интег-рирующее, колебательное.

3.1. Переходные характеристики звеньев

Переходная, или временная характеристика (функция) звена представ-

ляет собой реакцию на выходе звена, вызванную подачей на его вход еди-ничного ступенчатого воздействия. Переходная характеристика обознача-ется h(t).

Единичное ступенчатое воздействие (единичная ступенчатая функция) – это воздействие, которое мгновенно возрастает от нуля до единицы и да-лее остается неизменным. График переходного процесса, представленный на рис.2.10, иллюстрирует вышесказанное.

Единичное ступенчатое воздействие обозначается 1(t) и может быть описано следующим равенством:

.0tпри1;0tпри0

)t(1 (3.1)

Таким образом, если обозначить входной сигнал x(t), а выходной – y(t), то h(t) – это выражение для y(t) при x(t) = 1(t).

42

3.2. Частотные характеристики звеньев

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

Таким образом, если на вход звена направленного действия подано гармоническое воздействие

x = xmax sin t, (3.2) где xmax – амплитуда, а – угловая частота этого воздействия, то по окончании переходного процесса на выходе звена будут существо-вать гармонические колебания с той же частотой, что и входные колеба-ния, но отличающиеся в общем случае по амплитуде и фазе, т.е. в устано-вившемся режиме выходная величина звена

y = ymax sin (t + ), (3.3)

где ymax – амплитуда выходных установившихся колебаний; – фазовый сдвиг между входными и выходными колебаниями.

В связи с этим различают амплитудно-частотные и фазово-частотные характеристики звеньев.

Амплитудно-частотная характеристика (АЧХ) выражает отношение амплитуды колебаний на выходе звена к амплитуде колебаний на его входе в зависимости от частоты выходного сигнала:

Α(ω) = yмакс / xмакс = f(ω). (3.4) Фазово-частотная характеристика (ФЧХ) выражает зависимость разно-

сти фаз между входными и выходными колебаниями звена от частоты входного сигнала:

(ω) = f(ω). (3.5) Опережению фазы соответствует > 0, а отставанию < 0. Амплитудную и фазовую частотные характеристики можно объеди-

нить в одну характеристику – амплитудно-фазовую частотную характери-стику (АФЧХ), используя Α(ω) и (ω) в качестве полярных координат (рис.3.1). Каждая точка АФЧХ соответствует определенному значению частоты ω.

Так же возможно построение АФЧХ в прямоугольной системе коорди-нат, при этом координатами будут показанные на рис.3.1 проекции U и V вектора A на соответствующие оси. Зависимости U(ω) и V(ω) называются соответственно действительной (вещественной) и мнимой частотными ха-рактеристиками.

При этом АФЧХ можно рассматривать как одну из форм записи пере-даточных функций для случая синусоидального воздействия на вход звена. Для получения W(jω) из передаточной функции W(р) необходимо заме-нить оператор Лапласа p на jω. Затем необходимо выделить действитель-ную U(ω) и мнимую V(ω) составляющие.

43

Рис.3.1. Амплитудно-фазовая частотная характеристика

Также при исследовании САУ часто применяются логарифмические

частотные и фазовые характеристики. Логарифмической амплитудной характеристикой (ЛАХ) называют за-

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

Логарифмической фазовой характеристикой (ЛФХ) называют зависимость аргумента частотной характеристики от логарифма частоты и обозначают ().

Вместе ЛАХ и ЛФХ называют логарифмическими частотными харак-теристиками (ЛЧХ).

3.3. Получение переходных и частотных характеристик

Для получения переходных характеристик в пакете Simulink можно

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

Поэтому специально для построения переходных и частотных характе-ристик и определения их показателей служит дополнение к пакету Simulink – LTI Viewer.

В окне LTI Viewer можно выделить следующие части (рис.3.2): 1 − командное меню; 2 − основное окно для вывода графиков; 3 − строка состояния (здесь также выводятся подсказки).

j

=

Vi

Ui

i

Ai

i

= 0

1 2

44

Рис.3.2. Окно LTI Viewer

Командное меню содержит следующие пункты и подпункты (здесь, как и ранее, приводятся только нужные для дальнейшей работы с приложени-ем пункты):

1) File – операции с LTI-файлами: New Viewer – новое окно LTI Viewer (Ctrl+N); Import… – импорт LTI-файла из рабочей области; Export… – экспорт LTI-файла в рабочую область и на диск; Toolbox Preferences… – общая настройка LTI Viewer (для всех когда-либо запускаемых LTI Viewer); Page Setup… – настройка печати; Print… – печать (Ctrl+P); Print to Figure – преобразование графика в изображение, с открытием внутреннего графического редактора MATLAB. 2) Edit – операции редактирования: Plot Configurations… – управление выводом графиков (от одного до шес-ти графиков в одном окне); Systems – управление выводом графиков (если активно несколько LTI-файлов), содержит в себе две опции: Refresh – “освежить” список LTI-файлов и Delete… – стереть один или группу LTI-файлов (из списка);

1

2

3

45

Line Styles… – управление стилей линий графиков (выбор типов и цве-тов линий и маркеров); Viewer Preferences… – настройка текущего LTI Viewer. 3) Simulink – моделирование (связь LTI Viewer с пакетом Simulink): Get Linearized Model – запустить процесс линеаризации (моделирова-ния); Remove Input/Output Points – удалить точки входа/выхода. 4) Window – управление окнами LTI Viewer; 5) Help – переход к файлам справки по LTI Viewer.

Пункты подменю File > Toolbox Preferences… и подменю Edit > Viewer Preferences… идентичны, но если изменения, внесенные в пункт подменю Toolbox Preferences…, применяются ко всем окнам LTI Viewer, запускае-мым после внесения этих изменений, пункт Viewer Preferences… действует только на текущее окно LTI Viewer.

Окно настройки LTI Viewer, вызываемое из меню File > Toolbox Preferences… (рис.3.3), состоит из четырех вкладок, каждая из которых со-держит в себе ряд настроек: 1) Units – единицы измерения: Frequency in – частота в рад/с или Гц; а также масштаб частоты – лога-рифмический или линейный; Magnitude in – амплитуда в дБ или абсолютных единицах; Phase in - фаза в градусах или радианах. 2) Style – управление стилем (внешним видом) окна LTI Viewer: Grids – включить/выключить сетку на графиках; Fonts – установка высоты и стиля шрифтов, используемых для подписи графиков в LTI Viewer’е; Colors – выбор цветов LTI Viewer. 3) Characteristics – управление характеристиками графиков: Show settling time within – показывать время переходного процесса с рамками в n %, где n – требуемое число процентов (по умолчанию 2, в ТАУ принято использовать 5% рамки); Show rise time from n1 to n2, % – показывать время нарастания сигнала с n1 до n2 в процентах. 4) SISO Tool – управление настройками для пакета синтеза линейных сис-тем SISO Design Tool: Compensator Format – формат корректирующего звена: с постоянными времени или нулями/полюсами передаточной функции.

Окно настройки LTI Viewer, вызываемое из меню Edit > Viewer Prefer-ences…, отличается от меню File > Toolbox Preferences… наличием вклад-ки Parameters и отсутствием вкладки SISO Design Tool (рис.3.4).

46

Рис.3.3. Окно настроек Toolbox Preferences

Рис.3.4. Окно настроек Viewer Preferences

Вкладка parameters (параметры) содержит входные параметры для пе-реходной и частотных характеристик – соответственно для вектора време-ни :

● Generate automatically – автоматический выбор;

47

● Define stop time – определение конечного времени; ● Define vector – определение в виде вектора [начальное время: шаг:

конечное время]; и для вектора частоты:

● Generate automatically – автоматический выбор; ● Define range – определение начального и конечного значений часто-

ты; ● Define vector – определение в виде вектора logspace[10^начальное

значение: 10^конечное значение: количество точек]. Для получения переходных и частотных характеристик в пакете

Simulink с помощью LTI Viewer необходимо выполнить следующие дейст-вия:

1) выполнить команду меню Tools > Linear analysis, при этом откроют-ся два окна: LTI Viewer и Model_Inputs_and_Outputs;

2) из окна Model_Inputs_and_Outputs перенести в окно модели элемен-ты Input Point (точка входа) и Output Point (точка выхода). Элемент Input Point разместить на линии соединения между источником сигнала и сле-дующим блоком, а элемент Output Point – на линии соединения между приемником сигнала и предыдущим блоком. На рис.3.5 показан пример для модели, представленной апериодическим звеном с коэффициентом усиления k = 10 и постоянной времени T = 2 c (здесь необходимо отметить, что в пакете Simulink для обозначения оператора Лапласа принято приме-нять s = d/dt);

3) в окне LTI Viewer выполнить команду меню Simulink > Get Linearized Model. В окне появится график переходной характеристики мо-делируемого звена (рис.3.6);

4) для получения логарифмических частотных характеристик следует нажать правой кнопкой мыши на поле графика (не на линии графика). В открывшемся меню выбрать пункт Plot Type > Bode (рис.3.7);

5) для получения амплитудно-фазовых частотных характеристик сле-дует нажать правой кнопкой мыши на поле графика (не на линии графика). В открывшемся меню выбрать пункт Plot Type > Nyquist (рис.3.8).

После того как в основном поле LTI Viewer появится график, возни-кает дополнительная возможность для редактирования и просмотра графи-ков – контекстное меню (рис.3.9), вызываемое по нажатию правой кнопки мыши на поле графика (не на линии графика).

Контекстное меню LTI Viewer содержит следующие пункты: 1) Plot type – тип графика, отображаемого в основном окне LTI Viewer (из-меняется воздействие, подаваемое на вход исследуемой модели): Step – переходный процесс (единичное ступенчатое воздействие); Impulse – импульсная переходная характеристика (единичный импульс); • Bode – логарифмические частотные характеристики ЛАХ и ЛФХ (гармо-

48

Рис.3.5. Расположение элементов Input Point (точка входа) и Output Point (точка выхода) на модели

Рис.3.6. Переходный процесс для системы, изображенной на рис.3.5

49

Рис.3.7. Логарифмические частотные характеристики для системы, изображенной на рис.3.5

Рис.3.8. Амплитудно-фазовые частотные характеристики для системы, изображенной на рис.3.5

50

Рис.3.9. Контекстное меню LTI Viewer

ническое воздействие); Bode Mag. – отдельно логарифмическая амплитудная характеристика ЛАХ (гармоническое воздействие); Nyquist – амплитудно-фазовая частотная характеристика АФЧХ (гармо-ническое воздействие); Pole-Zero – диаграмма распределения полюсов-нулей передаточной функции модели. 2) Systems – список моделей, которые возможно отобразить, с возможно-стью включения/выключения отображения отдельных моделей. 3) Characteristics – параметры отображаемой переходной или частотной ха-рактеристики (численные значения).

Для переходных характеристик: Peak response – максимальное значение, перерегулирование и время, на котором отмечено максимальное значение выходной величины; Settling time – время переходного процесса; Rise time – время нарастания выходного сигнала; Steady state – установившееся значение переходного процесса.

Для частотных характеристик: Peak response – максимальное значение амплитуды и частота, на которой зафиксировано это значение; Stability margins min/all – точки устойчивости (минимум/все), отобража-ется в первом случае: gain margin – запас по амплитуде, phase margin – за-пас по фазе, частота среза и частота фазового сдвига (at frequency в соот-ветствующей точке), вывод об устойчивости системы в замкнутом виде. Во втором случае появляется характеристика delay margin - запас по за-паздыванию. 4) Grid – включение/выключение сетки графиков. 5) Zoom – масштабирование графиков: X-Y – по двум осям, с сохранением пропорций; In-X – по оси X; In-Y – по оси Y; Out – вернуть первоначальные (автоматически определенные). 6) Properties… – свойства текущего графика (рис.3.10).

51

Рис.3.10. Окно свойств текущего графика

Окно свойств текущего графика содержит пять вкладок, содержимое которых изменяется в зависимости от того, какой график является теку-щим. В общем случае: 1) Labels – подписи на поле графика: Title – наименование графика; X-Label – подпись оси абсцисс; Y-Label – подпись оси ординат. 2) Limits – пределы (присутствует опция автомасштаба – Auto-scale): X-Limits – предельные (начальные и конечные) значения для оси абс-цисс; Y-Limits – предельные (начальные и конечные) значения для оси орди-нат. 3) Units – единицы измерения (аналогичные в пункте меню File > Toolbox Preferences…). 4) Style – управление стилем (внешним видом) окна LTI Viewer (аналогич-ные в пункте меню File > Toolbox Preferences…). 5) Characteristics – показать/убрать характеристики текущего графика (ана-логичные пункту Characteristics контекстного меню LTI Viewer).

3.4. Идеальное усилительное (безынерционное) звено

Усилительным (пропорциональным) называют звено, у которого вы-

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

52

у = kx, (3.6) где k - коэффициент передачи звена (в общем случае размерная величина).

Передаточная функция усилительного звена:

.k)s(X)s(Y)s(W (3.7)

Модель для исследования динамических свойств звена состоит из сле-дующих блоков: Step, Gain, Scope.

Полученная модель показана на рис.3.11.

Рис.3.11. Модель усилительного звена (k=1)

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

Единственным значимым для задач данного учебного пособия пара-метром является значение Gain – коэффициент передачи усилительного звена. Остальные параметры рекомендуется оставить без изменения.

Изменение коэффициента передачи незамедлительно сказывается на внешнем виде звена – изменяется цифра, которая и является тем самым ко-эффициентом передачи. В случае, если значение коэффициента передачи не вмешается в значок блока, оно заменяется на -k-.

Переходные и частотные характеристики звена представлены на рис.3.13. Из приведенных характеристик видно, что усилительное звено не оказывает влияние на фазовые частотные характеристики и не изменяет наклон амплитудных частотных характеристик, вследствие чего ампли-тудно-фазовая частотная характеристика идеального усилительного звена – точка. Для того чтобы убедиться, что оно влияет на установившееся зна-чение переходного процесса, изменим его коэффициент передачи, устано-вив k = 2 (рис.3.14).

Из характеристик, приведенных на рис.3.15, видно, что установившееся значение переходного процесса теперь равно 2, то есть установившееся значение выходного сигнала равно коэффициенту передачи усилительного

53

Рис.3.12. Окно параметров усилительного звена Gain

Рис.3.13. Переходные и частотные характеристики усилительного звена (k=1)

54

Рис.3.14. Модель усилительного звена (k=2)

Рис.3.15. Переходные и частотные характеристики усилительного звена (k=2) сигнала, Действительно для усилительного звена выходной сигнал y = kx, а в случае переходного процесса входным сигналом является x = 1(t), поэто-му yуст = k1 = k = 2.

Для логарифмической амплитудной характеристики A(1) = 20 log(k) = 20 log (2) = 6,02.

55

3.5. Апериодическое (инерционное) звено

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

Уравнение апериодического звена:

kxydtdyT , (3.8)

где Т – постоянная времени звена, характеризующая инерционность звена, с; k – коэффициент усиления звена.

В операторной форме: (Τs + 1) y = kx. (3.9)

Передаточная функция апериодического звена:

.1Ts

k)s(X)s(Y)s(W

(3.10)

Модель для исследования динамических свойств звена состоит из сле-дующих блоков: Step, Transfer Fcn, Scope.

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

щелкнуть на нем левой кнопкой мыши, появится окно параметров звена Transfer Fcn (рис.3.17).

Параметры звена Transfer Fcn: Numerator – числитель и Denominator – знаменатель передаточной функции в виде полиномов a0sn + a1sn-1 +…+ an-1s0 и b0sm + a1sm-1 +…+ am-1s0 соответственно. В общем случае с помощью этого звена можно задать любую передаточную функ-цию, представленную в виде отношения полиномов, причем n m, где n – порядок полинома в числителе, m – порядок полинома в знаменателе.

Рис.3.16. Модель апериодического звена (k = 1, T = 1)

56

Рис.3.17. Окно параметров звена Transfer Fcn

Рис.3.18. Переходные и частотные характеристики апериодического звена (k = 1, T = 1)

У апериодического звена в числителе находится один параметр – ко-эффициент усиления k, поэтому в окне Numerator для апериодического звена вводим [k].

Знаменатель апериодического звена Ts +1, то есть полином вида a0s1 + +a1s0, где a0 = T, а1 = 1. Поэтому в окне Denominator для апериодического звена вводим [T 1] (через пробел).

Переходные и частотные характеристики апериодического звена при-ведены на рис.3.18.

57

3.6. Колебательное звено

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

Уравнение колебательного звена в операторной форме записи: ,kxy1sTsT 2

221 (3.11)

где T1 и T2 – постоянные времени колебательного звена; k – коэффициент усиления колебательного звена.

Передаточная функция колебательного звена:

.1sTsT

k)s(X)s(Y)s(W

22

1 (3.12)

Постоянные времени колебательного звена T1 и T2 связаны зависимо-стью

= Т2 / 2Т1. (3.13) Коэффициент называют коэффициентом колебательности и, судя по

названию, он характеризует колебательность переходного процесса коле-бательного звена.

Если < 1, то переходный процесс звена – колебательный, и чем меньше единицы, тем колебательнее процесс.

Если 1, то получаем апериодическое звено второго порядка. Модель для исследования динамических свойств звена состоит из сле-

дующих блоков: Step, Transfer Fcn, Scope. Полученная модель показана на рис.3.19.

Рис.3.19. Модель колебательного звена (k = 1, T12 = 1, T2 = 1)

Методика изменения коэффициентов звена Transfer Fcn приведена в

п.3.5 для апериодического звена, единственное отличие – для колебатель -

58

Рис.3.20. Переходные и частотные характеристики колебательного звена (k = 1, T1

2 = 1, T2 =1)

ного звена в окне Denominator вводим [T12 T2 1] (также через пробел), по-

скольку в знаменателе колебательного звена находится полином вида a0s2 + a1s1 + a2s0, где a0 = T1

2, а1 = T1, a2 = 1. Переходные и частотные характеристики колебательного звена пред-

ставлены на рис.3.20.

3.7. Интегрирующее звено Интегрирующим называют звено, в котором выходная величина про-

порциональна интегралу во времени от входной величины. Уравнение интегрирующего звена:

sy = kx, (3.14) где k – коэффициент усиления интегрирующего звена.

В интегральной форме:

t

00.xdtxkx

sky (3.15)

При ступенчатом входном сигнале выходная величина линейно зависит от времени:

у = kxt = Kt, (3.16) где K = kx - постоянная величина; t - время.

59

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

Уравнение звена в операторной форме:

.xsky (3.17)

Передаточная функция интегрирующего звена:

.sk

)s(X)s(Y)s(W (3.18)

Рис.3.21. Модель интегрирующего звена (k = 1)

Рис.3.22. Переходные и частотные характеристики интегрирующего звена (k=1)

60

Модель для исследования динамических свойств звена состоит из сле-дующих блоков: Step, Gain, Integrator, Scope. В связи с тем, что передаточ-ная функция интегрирующего звена в пакете Simulink всегда 1/s, то для то- го, чтобы задать коэффициент усиления интегрирующего звена, использу-ется последовательное соединение усилительного (Gain) и интегрирующе-го звеньев (Integrator).

Полученная модель показана на рис. 3.21. Учитывая все вышесказанное, для изменения коэффициента усиления

интегрирующего звена необходимо изменять коэффициент передачи по-следовательно соединенного с ним усилительного звена (Gain). Параметры интегрирующего звена (Integrator) изменять не рекомендуется.

Переходные и частотные характеристики интегрирующего звена пред-ставлены на рис.3.22.

3.8. Дифференцирующее звено

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

Уравнение идеального дифференцирующего звена в операторной фор-ме записи:

у = ksx. (3.19) При подаче на вход звена ступенчатого сигнала на выходе получается

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

Передаточная функция идеального дифференцирующего звена:

.ks)s(X)s(Y)s(W (3.20)

Модель для исследования динамических свойств звена состоит из сле-дующих блоков: Step, Gain, Switched derivative for linearization (из раздела Linearization библиотеки Simulink Extras), Scope. Здесь используется звено Switched derivative for linearization (из раздела Linearization библиотеки Simulink Extras), а не звено Derivative из раздела Continuous, поскольку последнее не оказывает нужного влияния на частотные характеристики системы.

В связи с тем, что передаточная функция дифференцирующего звена в пакете Simulink всегда s, то для того чтобы задать коэффициент усиления интегрирующего звена, используется последовательное соединение усили-тельного (Gain) и дифференцирующего звеньев (Switched derivative for

61

linearization). Полученная модель показана на рис. 3.23.

Рис.3.23. Модель дифференцирующего звена (k =1)

Учитывая вышесказанное, для изменения коэффициента усиления дифференцирующего звена необходимо изменять коэффициент передачи, последовательно соединенного с ним усилительного звена (Gain). В пара-метрах блока Switched derivative for linearization устанавливаем: Switch value = 1 и Derivative constant = 107 (чем больше последнее значение, тем точнее звено, но дольше время моделирования).

Рис.3.24. Переходные и частотные характеристики дифференцирующего звена (k=1)

62

Причем логарифмические характеристики для данного звена можно

считать достоверными на интервале от 100 до частоты, на два порядка меньшей, чем значение Derivative constant, заданное в окне параметров блока Switched derivative for linearization (рекомендуем установить это зна-чение, равным 107, тогда график достоверен на интервале до 105).

Характеристики звена представлены на рис. 3.24.

3.9. Звено с постоянным запаздыванием

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

Звено определяется как запаздывающее, если оно описывается уравне-нием

y(t) = x(t – ), (3.21) где - время запаздывания.

Передаточная функция звена запаздывания:

.e)p(W s (3.22)

Модель для исследования динамических свойств звена составляется из следующих блоков: Step, Transport Delay, Scope (рис. 3.25).

Рис.3.25. Модель звена с постоянным запаздыванием (t =1)

Для изменения параметров звена необходимо сделать двойной щелчок левой кнопкой мыши по элементу Transport Delay в окне модели. В от-

63

крывшемся окне в поле Time delay (время запаздывания) вводим нужное значение.

Рис.3.26. Переходные и частотные характеристики звена с постоянным запаздыванием (t = 1, Pade Order = 2)

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

64

Переходные и частотные характеристики звена с постоянным запазды-ванием представлены на рис.3.26.

Из представленных характеристик видно, что данное звено не совсем точно отрабатывается в LTI Viewer – переходный процесс далек от иде-ального, но при просмотре переходного процесса при помощи блока Scope он выглядит близко к идеальному (рис.3.27). В версии MATLAB 6.5 для устранения этого недостатка в библиотеке Simulink Extras имеется допол-нительная группа блоков Linearisation (рис. 3.28), один из которых – Switch transport delay for linearisation предназначен для использования в моделях с LTI Viewer.

Рис. 3.28. Блоки Linearisation

Для увеличения точности данного звена необходимо изменять его на-

стройку Pade Order (порядок функции-аппроксимации Паде), чем больше это число, тем ближе переходный процесс к идеалу и тем дольше ведется расчет. Наиболее оптимальным числом является Pade Order = 10.

Логарифмическая фазовая частотная характеристика не уходит в бес-конечность и достоверна на протяжении 1,5 декады (до точки перегиба).

65

4. СПОСОБЫ СОЕДИНЕНИЯ ЗВЕНЬЕВ В СТРУКТУРНУЮ СХЕМУ

4.1. Последовательное соединение звеньев

При последовательном соединении звеньев выходная величина преды-

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

n

1ii sWsW . (4.1)

Рис.4.1. Последовательное соединение звеньев в пакете

Simulink

В пакете Simulink последовательное соединение получается простым соединением звеньев (рис.4.1).

4.2. Параллельное соединение звеньев

При параллельном соединении на вход всех звеньев подается общий

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

n

1ii sWsW . (4.2)

В пакете Simulink параллельное соединение получаем с помощью сум-матора (Sum), который находится в разделе Math основной библиотеки Simulink (рис.4.2).

Рис.4.2. Параллельное соединение звеньев в пакете Simulink

66

4.3. Звено, охваченное обратной связью

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

sWsW1

sWsW21

1

. (4.3)

Знак (−) в этой формуле ставится при положительной обратной связи, т.е. при сложении сигналов, а знак (+) при отрицательной обратной связи, т.е. при вычитании сигналов. Если второе звено отсутствует в цепи обрат-ной связи, то передаточная функция примет вид

sW1

sWpW1

1

. (4.4)

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

Рис.4.3. Звено, охваченное обратной связью,

в пакете Simulink В случае охвата обратной связью группы звеньев, что характерно для

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

,sW1

sWsWp

з

(4.5)

где W(s) – передаточная функция цепочки звеньев от входного до выход-ного воздействия; Wр(s) – передаточная функция разомкнутой системы, то есть цепочки из всех звеньев системы, получающейся после разрыва об-ратной связи.

Для системы, изображенной на рис. 4.3, W(s) = W1(s), а Wp(s) = W1(s)W2(s) и передаточная функция замкнутой системы превраща-ется в уравнение (4.3).

67

5. РЕШЕНИЕ ЗАДАЧ АНАЛИЗА ЛИНЕЙНЫХ СИСТЕМ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ

Задача анализа линейных систем автоматического управления включа-

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

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

5.1. Понятие устойчивости и критерии устойчивости систем

автоматического управления

Устойчивость – это свойство системы возвращаться в исходный или близкий к нему установившийся режим после всякого выхода из него в ре-зультате какого-либо воздействия.

Это понятие хорошо описывается тремя простыми примерами (рис.5.1): 1. Шар, находящийся на вершине горки (рис.5.1, а), является примером

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

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

3. Шар, находящийся на плоскости, является примером системы, нахо-дящейся на границе устойчивости, так как при приложении к нему конеч-ного возмущения он остановится, но уже в другом положении (рис.5.1, в).

Рис.5.1. Пример, иллюстрирующий понятие устойчивости Таким образом, если система неустойчива, то достаточно любого толч-

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

а) б) в)

68

системе, вызванный каким-либо воздействием, со временем затухает, и система вновь возвращается в устойчивое состояние.

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

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

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

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

(a0sn+a1sn-1+...+an-1)y=0. (5.1) При этом система автоматического управления будет устойчива, если

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

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

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

Критерии устойчивости делятся на две разновидности: алгебраические и частотные. Алгебраические критерии являются аналитическими, а час-

69

тотные – графоаналитическими. Критерии устойчивости позволяют также оценить влияние параметров системы на устойчивость в виде запасов ус-тойчивости.

В пакете Simulink возможно автоматическое применение трех критери-ев устойчивости:

1) по корням характеристического уравнения системы; 2) частотные критерии устойчивости Найквиста; 3) критерий устойчивости Никольса. В нашей стране применяются в основном первые два.

5.1.1. Критерий устойчивости по корням

характеристического уравнения Решение уравнения (5.1) представляет собой сумму слагаемых, вид ко-

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

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

Рис.5.2. Определение устойчивости САУ по корням характеристического уравнения

70

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

В пакете Simulink устойчивость системы по корням характеристиче-ского уравнения определяется с использованием пакета расширения LTI Viewer путем нажатия правой кнопки на поле графика и выбора пункта контекстного меню Plot Type > Pole/Zero (рис.5.2).

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

5.1.2. Частотный критерий устойчивости Найквиста

Частотный критерий устойчивости Найквиста позволяет оценить ус-

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

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

Для амплитудно-фазовой частотной характеристики это определение звучит следующим образом: для устойчивости замкнутой системы необ-ходимо и достаточно, чтобы АФЧХ разомкнутой системы, при измене-нии частоты от 0 до , не охватывала точку с координатами (−1, j0).

Если АФЧХ разомкнутой системы проходит через точку с координата-ми (−1, j0), то система находится на границе устойчивости.

Для логарифмических частотных характеристик это определение имеет другой вид: для устойчивости замкнутой системы необходимо и доста-точно, чтобы ЛАХ разомкнутой системы пересекла ось абсцисс раньше, чем ЛФХ, спадая окончательно, перешла за значение − (−180).

Точка, в которой ЛАХ пересекает ось абсцисс, называется частотой среза и обозначается с.

71

Точка, в которой ЛФХ пересекает значение − (−180) спадая оконча-тельно, называется частотой фазового сдвига и обозначается .

Для вышепринятых обозначений условия устойчивости по критерию Найквиста для ЛЧХ выглядят следующим образом:

1) если с < , то система является устойчивой; 2) если с > , то система является неустойчивой; 3) если с = , то система находится на границе устойчивости. В пакете Simulink для определения устойчивости системы с наличием

обратной связи (замкнутой системы) по критерию Найквиста необходимо разомкнуть цепь обратной связи, то есть превратить систему в разомкну-тую. Для этого нужно удалить одну из соединительных линий в цепи об-ратной связи и переставить точку выхода (для пакета LTI Viewer). На рис.5.3 в качестве примера приведена структурная схема замкнутой систе-мы стабилизации угловой скорости турбогенератора, а на рис.5.4 – струк-турная схема этой же системы, но в разомкнутом виде.

Рис.5.3. Пример замкнутой системы

Рис.5.4. Пример разомкнутой системы

В пакете Simulink устойчивость системы по критерию Найквиста опре-

деляется с использованием пакета расширения LTI Viewer с помощью на-жатия правой кнопки на поле графика и выбора из контекстного меню:

1) для АФЧХ: Plot Type > Nyquist (рис.5.5); 2) для ЛЧХ: Plot Type > Bode (рис.5.6).

72

Рис.5.5. Определение устойчивости САУ по АФЧХ

Рис.5.6. Определение устойчивости САУ по ЛЧХ

Из рис.5.5 видно, что АФЧХ разомкнутой системы не охватывает точку

с координатами (−1, j0) (обозначена на графике красным крестиком), сле-довательно, замкнутая система стабилизации угловой скорости турбогене-ратора, структурная схема которой приведена на рис.5.3 и 5.4, является ус-

73

тойчивой. Такой же вывод можно сделать из анализа ЛЧХ (см. рис.5.6), так как ЛАХ пересекает ось абсцисс раньше, чем ЛФХ, окончательно спадая, переходит за значение − (−180). Для данной системы c = 1,2 рад/с, а = 1,39 рад/с, то есть c < и, следовательно, система устойчива.

5.2. Определение запасов устойчивости систем

автоматического управления

Вывод об устойчивости системы автоматического управления не явля-ется достаточным для определения устойчивости, поскольку существует вопрос: «Насколько устойчива система?». Для ответа на данный вопрос необходимо определить так называемые показатели запасов устойчивости.

В случае применения частотных критериев устойчивости Найквиста та-ких показателей два (рис.5.7):

1) запас устойчивости по амплитуде L, измеряемый в децибелах (дБ) и показывающий, на сколько следует сдвинуть ЛАХ, чтобы система оказа-лась на границе устойчивости (см. рис.5.7, пунктирная линия);

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

Запасы устойчивости получаются следующим образом (рис.5.7): 1) из точки пересечения ЛФХ, при окончательном ее спаде, c линией,

проведенной на уровне −180 (частота фазового сдвига), поднимается пер-пендикуляр до пересечения с ЛАХ. Расстояние от этой точки до оси абс-цисс даст искомое значение запаса устойчивости по амплитуде (L, дБ);

2) из точки пересечения ЛАХ оси абсцисс (частота среза) опускается перпендикуляр до пересечения с ЛФХ, расстояние от этой точки до уровня −180 даст искомое значение запаса устойчивости по фазе (, ).

В пакете Simulink численные значения запасов устойчивости системы автоматического управления можно определить по ЛЧХ или АФЧХ. Для получения ЛЧХ нужно в пакете LTI Viewer нажать правой кнопкой мыши на график и выбрать пункт меню Plot Type > Bode (см. рис.5.6).

Для определения запасов устойчивости необходимо щелкнуть правой кнопкой мыши на поле графика, и выбрать пункт контекстного меню Char-acteristics > Stability Margins (Min).

Для получения численных значений запасов устойчивости необходимо подвести курсор мыши к выделенным на графике точкам (рис.5.8 и 5.9) для системы, представленной на рис.5.4.

При этом в появляющихся окошках указаны численные значения запа-сов и частота в соответствующих точках логарифмических частотных ха-рактеристик:

74

Рис.5.7. Определение запасов устойчивости

Рис.5.8. Определение запасов устойчивости по ЛЧХ

L, дБ

-,

L, дБ

,

lg , рад/с

lg , рад/с

0

0

-180

с

75

Рис.5.9. Определение запасов устойчивости по АФЧХ Gain margin (dB) – запас устойчивости по амплитуде, дБ; Phase margin (deg) – запас устойчивости по фазе, градусы; At frequency (rad/sec) – “на частоте”, рад/с.

Кроме того, делается вывод об устойчивости замкнутой системы Closed-loop stable? (Замкнутая система устойчива?), который может при-нимать два значения: Yes (Да) и No (Нет).

5.3. Показатели качества переходных процессов

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

Ранее, при построении характеристик систем автоматического управ-ления аналитически, с целью сокращения затрат времени пользовались

76

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

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

● Время регулирования (длительность переходного процесса) – интер-вал времени, по истечении которого отклонение переходной характеристи-ки от установившегося значения не превышает некоторой заданной вели-чины. Обычно отклонение управляемой величины от установившегося значения должно быть не более 5%.

В пакете Simulink время регулирования определяется в пакете расши-рения LTI Viewer с помощью нажатия правой кнопки на поле графика пе-реходного процесса (Step Response) и выбора пункта меню Characteristics > Settling Time (см. рис.5.10) для системы, представленной на рис.5.3.

При этом на графике появляется точка, и после того как подвести к ней указатель мыши, появляется окошко с числовым значением времени пере-ходного процесса (Settling Time).

Рис.5.10. Определение времени переходного процесса

В пакете расширения LTI Viewer по умолчанию установлено отклоне-

ние от заданной величины 2%. Для того чтобы изменить это значение, не-обходимо в командном меню LTI Viewer в пункте Edit выбрать Viewer

77

Preferences. На закладке Characteristics открывшегося окна установить зна-чение Show settling time within, равное 5%.

● Перерегулирование () – максимальное отклонение управляемой ве-личины от установившегося значения, выраженное в процентах, – характе-ризует колебательность переходного процесса.

%,100y

yy

уст

устmax

(5.2)

где ymax – максимальное значение выходного сигнала; yуст – установившее-ся значение выходного сигнала.

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

В пакете Simulink перерегулирование определяется в пакете расшире-ния LTI Viewer с помощью нажатия правой кнопки на графике переходно-го процесса (Step Response) и выбора пункта меню Characteristics > Peak Response (рис.5.11, для системы на рис.5.3).

При этом на графике появляется точка, после подвода к которой указа-теля мыши, появляется окошко с числовыми значениями: − Peak amplitude – максимальное значение выходного сигнала ymax; − Overshoot – перерегулирование , %; − At time – время, при котором зафиксировано максимальное значение вы-ходного сигнала.

Рис.5.11. Определение перерегулирования

78

● Установившееся значение выходного сигнала yуст, определяющее статическую точность системы. В системе MATLAB установившееся зна-чение определяется в пакете расширения LTI Viewer с помощью нажатия правой кнопки на графике переходного процесса (Step Response) и выбора пункта меню Characteristics > Steady State (рис.5.12, для системы на рис.5.3).

При этом на графике появляется точка, и после того как подвести к ней указатель мыши, появится окошко с числовым значением установившегося значения выходного сигнала (DC Gain), которое равно коэффициенту пе-редачи замкнутой системы.

● Статическая ошибка системы ст показывает на сколько отличается уровень выходного сигнала от заданного его значения:

ст = yуст – yзад , (5.3) где yзад – заданное значение выходного сигнала.

Поскольку для получения переходных процессов используется единич-ный ступенчатый сигнал, то для управляющего воздействия yзад = 1.

● Частота резонанса системы – частота, на которой в системе наблю-даются колебания р с наибольшей амплитудой Аз(р).

Рис.5.12. Определение установившегося значения выходного сигнала

79

Рис.5.13. Определение резонансной частоты и амплитуды системы

В пакете Simulink частоту резонанса можно найти с использованием

пакета расширения LTI Viewer с помощью нажатия правой кнопки на поле графика и выбора пункта меню Plot Type > Bode Mag., повторного нажатия

Рис.5.14. Определение показателя колебательности

80

правой кнопки на поле графика и выбора Characteristics > Peak Response (рис.5.13). Затем нужно подвести курсор мыши к появившейся точке. В ре-зультате появится окошко со значениями амплитуды (Peak gain) и резо-нансной частоты (At frequency) этой точки.

● Колебательность показывает вид переходного процесса и оценивает его колебательность, которая оценивается отношением соседних максиму-мов:

%,100yyK

2max

1max (5.4)

где ymax1 – значение первого пика переходной характеристики; ymax2 – значение второго пика переходной характеристики.

Для пакета Simulink удобно пользоваться отношением максимальной амплитуды выходного сигнала ymax и амплитуды соседнего с ним пика пе-реходной характеристики, который необходимо снять с графика переход-ного процесса (5.14, для системы на рис.5.3).

В случае, если показатель колебательности К = 1 – переходный процесс апериодический; если К > 1 – колебательный сходящийся, если К < 1, то имеет место расходящийся колебательный процесс. Чем показатель коле-бательности больше единицы, тем переходный процесс более колебатель-ный.

5.4. Общая методика анализа линейных систем автоматического управления

Для проведения анализа произвольной системы автоматического

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

1. Построить принципиальную или кинематическую схему САУ (рис.5.15) и выделить на ней отдельные функциональные элементы.

В схему входят следующие элементы: Т – турбина (объект управле-ния); МП – механическая передача; ЦБМ – центробежный маятник (вы-полняет функции датчика, определяющего действительную угловую ско-рость, задатчика угловой скорости и сумматора, вычисляющего ошибку регулирования); ГУ – гидроусилитель; ГСП – гидравлический сервопри-вод; ПВК – паровпускной клапан с регулируемым расходом; Г – генератор, приводимый во вращение турбиной.

Некоторой установившейся нагрузке генератора соответствует опреде-ленная угловая скорость турбины, положение грузов ЦБМ и поршня ГСП, а также величина открытия ПВК. Если нагрузка уменьшится, угловая ско-рость турбины возрастет, грузы ЦМБ, а вместе с ними и муфта перемес-тятся вверх (z). Переместится и золотник ГУ (s1), поршень ГСП пока оста-

81

ется неподвижным (s2 = 0), так как закрыты окна ГУ. Он приоткроет дос-

Рис.5.15. Принципиальная схема системы стабилизации угловой

скорости турбогенератора туп рабочего тела (жидкости) в верхнюю полость ГСП, поршень которого начнет перемещаться вниз и прикроет ПВК, который, в свою очередь, уменьшит подачу пара в турбину, вследствие чего угловая скорость турби-ны уменьшится. Одновременно, из-за наличия жесткой обратной связи между ГСП и ГУ (рычаг, связывающий ГСП, ГУ и ЦБМ), золотник ГУ пе-реместится вниз и прикроет окна ГУ и движение поршня ГСП прекратится. Система вновь будет в некотором установившемся режиме.

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

2. На основе принципиальной или кинематической схемы составить функциональную схему системы (рис.5.16).

Рис.5.16. Функциональная схема системы стабилизации угловой скорости турбогенератора

зад e z х1

х2

х m p ЦБМ ГУ ГСП ПВК Т

ОС

mн

ЦБМ

Т

Г

МП

z

s

ГСП

m ГУ

ПВК

Пар

Мн

82

На функциональной схеме показаны следующие сигналы: e – ошибка регулирования (отклонение); z – перемещение муфты ЦБМ; х1 – переме-щение золотника ГУ от ЦБМ; х2 – перемещение золотника ГУ от ГСП; х – суммарное перемещение золотника ГУ; m – перемещение ПВК, пропор-циональное перемещению поршня ГСП; p – давление пара на входе в тур-бину; – угловая скорость вращения вала турбины; зад – заданная угло-вая скорость вращения вала турбины; mн – момент нагрузки на валу тур-бины.

В любой момент времени с помощью сумматора выполняется сравне-ние заданного и фактического значений угловой скорости вращения вала турбины (вычисляется ошибка регулирования):

e = зад − . (5.5)

Если ошибка регулирования не равна нулю, то происходят следующие

перемещения: муфты ЦБМ (z), золотника ГУ (s), поршня ГСП (m), штока ПВК, вследствие чего изменяется давление на входе в турбину (p). Знаки перемещений зависят от знака ошибки регулирования.

В результате чего изменяется фактическая угловая скорость вращения вала турбины и ошибка регулирования стремится к нулю.

Таким образом, в данной схеме: ● ЦБМ исполняет функции: датчика угловой скорости, задатчика и

сумматора (сравнивающего элемента); ● ГУ и ГСП являются ПИ-регулятором; ● ПВК является исполнительным устройством, изменяющим подачу

пара в турбину; ● турбина Т является объектом управления. 3. Произвести математическое описание элементов функциональной

схемы в виде передаточных функций. Уравнения отдельных элементов в отклонениях имеют следующий вид:

;ekzdtdzT

dtzdT 122

22

1

;zkx 21 ;xxx 21

;xkdtdm

3 (5.6)

;mkpdtdpT 43

;mkpkdtdT н554

83

.mkx 62 Далее преобразуем уравнения в передаточные функции:

;

1sTsTk

sEsZsW

222

1

11

;ksZsSsW 2

12

;

sk

sSsMsW 3

3

;

1sTk

sMsPsW

3

44

(5.7)

;

1sTk

sPssW

4

55

;

1sTk

sMssW

4

5

н5

.ksMsSsW 6

26

Числовые значения параметров: k1 = 10; k2 = 0,6; k3 = 2; k4 = 0,8; k5 = 1;

k`5 = 0,2; k6 = 0,4; Т12 = 0,001 с2; Т2 = 0,2 с; Т3 = 5 с; Т4 = 0,3 с.

4. На основании функциональной схемы строим структурную схему сис-темы (рис.5.17).

Рис.5.17. Структурная схема системы стабилизации угловой скорости турбогенератора

зад e z х1

х2

х m p ЦБМ ГУ ГСП ПВК

Т

ОС

mн

W1 W2 W3 W4 W5

W6

W`5

84

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

5. Собрать структурную схему в пакете Simulink, используя типовые блоки (процесс сборки схемы и типовые блоки подробно описаны в треть-ей главе данного учебного пособия).

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

Для анализа системы по управляющему воздействию точки входа и выхода (Input and Output point) должны быть расположены так как, показа-но на рис.5.3 и 5.4, то есть точка входа располагается в месте приложения управляющего воздействия, а выхода – в месте расположения выходного сигнала.

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

6. Используя разомкнутую структурную схему системы, без возму-щающего воздействия (рис.5.4) с помощью частотного критерия Найкви-ста (см. рис.5.5 и 5.6) сделать вывод об устойчивости системы по управ-ляющему воздействию и получить значения запасов устойчивости (см. рис.5.8 и 5.9).

Для анализируемой системы можно сделать вывод об ее устойчивости, так как АФЧХ разомкнутой системы не охватывает точку (−1, j0) и частота среза (с = 1,2 рад/с) меньше частоты фазового сдвига ( = 1,39 рад/с).

Запасы устойчивости: - по амплитуде L = 2,37 дБ; - по фазе = 9,75 . 7. Используя замкнутую структурную схему системы без возмущаю-

щего воздействия (см. рис.5.3), получить показатели качества. Для анализируемой системы показатели качества: - установившееся значение выходной величины yуст = 0,923; - ошибка регулирования ст = 0,923 − 1 = − 0,077 (так как на входе

единичное ступенчатое воздействие, то yзад = 1); - перерегулирование = 76,3 %; - время регулирования tр = 36,9 с.

85

8. Провести анализ влияния возмущающего воздействия на переход-ный процесс. Для этого необходимо добавить в структурную схему звено с передаточной функцией W`5 (отрицательный знак в передаточной функции W`5 учитывается с помощью сумматора), источник единичного ступенча-того воздействия и разместить еще одну точку входа (Input Point1) (рис.5.18).

Рис.5.18. Структурная схема системы стабилизации угловой скорости

турбогенератора с возмущающим воздействием

При этом система анализируется также с использованием LTI Viewer, но в его окне теперь содержится два графика переходных процессов (рис.5.19) – от управляющего (левый) и от возмущающего (правый).

Рис.5.19. Анализ влияния возмущающего воздействия на

переходный процесс

86

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

- максимальное значение ymax = − 0,174; - длительность переходного процесса tпп = 36,3 с; - установившееся значение yуст = − 0,0154; - статическая ошибка ст = yуст - yзад = − 0,0154 - 0 = − 0,0154. При отработке системой автоматического управления возмущающего

воздействия эта система должна стремиться к тому, чтобы возмущающее воздействие не влияло на статический режим, то есть ст должна быть как можно меньше или, в идеале, равна нулю (поэтому для возмущающего воз-действия yзад = 0). Для этого используют различные способы коррекции (управление по отклонению и т.д.).

В рассматриваемом примере в результате приложения к системе одно-временно управляющего и возмущающего воздействий установившееся значение выходного сигнала будет равно сумме значений yуст для этих воз-действий: yуст = 0,923 – 0,0154 = 0,9076. Таким же образом можно подсчи-тать и суммарную статическую ошибку ст = − 0,077 − 0,0154 = − 0,0924.

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

Приведенная в примере система стабилизации угловой скорости турбо-генератора является системой судовой автоматики, к которой предъявля-ются следующие требования:

- запас устойчивости по амплитуде: не менее 20 дБ; - запас устойчивости по фазе: от 20 до 80; - перерегулирование: не более 10%; - статическая ошибка: не более ±5% (±0,05). Сравниваем полученные значения с требуемыми: - запас устойчивости по амплитуде: менее 20 дБ и 2,37 дБ; - запас устойчивости по фазе: от 20 до 80 и 9,75; - перерегулирование: 10% и 76,3%; - статическая ошибка: 0,05 и −0,077 – не учитывая возмущающее воз-

действие и −0,0924, учитывая возмущающее воздействие. Таким образом, анализируемая в примере система с указанными на-

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

87

6. РЕШЕНИЕ ЗАДАЧ СИНТЕЗА ЛИНЕЙНЫХ СИСТЕМ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ, КОРРЕКЦИЯ

ИХ ДИНАМИЧЕСКИХ СВОЙСТВ

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

Задача синтеза систем автоматического управления (коррекция их ди-

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

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

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

Корректирующие звенья по способу включения можно разделить (рис.6.1):

1) на последовательные (рис.6.1, а); 2) параллельные (рис.6.1, б); 3) стоящие в обратной связи (местные или главная обратные связи)

(рис.6.1, в).

Рис.6.1. Способы включения корректирующих звеньев в систему: Wо(s) – передаточная функция основных звеньев системы; Wк(s) – передаточная функция корректирующего звена

Действие корректирующих звеньев сводится к следующему: 1) введение в контур САУ пропорциональных воздействий либо

а) Wк(s) Wo(s)

б)

Wо(s)

Wк(s)

в)

Wо(s)

Wк(s)

88

воздействий по производным и/или интегралам (дифференцирующие и/или интегрирующие звенья);

2) введение корректирующих воздействий параллельно отдельным звеньям;

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

4) введение корректирующих воздействий в функции внешних воздей-ствий и их производных (управление по возмущению).

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

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

Пропорциональные (П), имеющие передаточную функцию:

,ksW пп (6.1)

где kп – коэффициент усиления пропорционального регулятора. Пропорционально-дифференцирующие (ПД):

,skksW дппд (6.2)

где kп – коэффициент усиления пропорциональной составляющей регуля-тора; kд – коэффициент усиления дифференцирующей составляющей регу-лятора.

Пропорционально-интегрирующие (ПИ):

,s

kksW иппи (6.3)

где kп – коэффициент усиления пропорциональной составляющей регуля-тора; kи – коэффициент усиления интегрирующей составляющей регулято-ра.

Пропорционально-интегродифференцирующие (ПИД):

,sks

kksW ди

ппид (6.4)

где kп – коэффициент усиления пропорциональной составляющей регуля-тора; kи – коэффициент усиления интегрирующей составляющей регулято-

89

ра; kд – коэффициент усиления дифференцирующей составляющей регуля-тора.

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

С учетом инерционности регуляторов (реальные регуляторы): Пропорциональные (П):

,1sT

ksWп

пп

(6.5)

где Tп – постоянная времени пропорционального регулятора.

Пропорционально-дифференцирующие (ПД):

,1sTskk

sWпд

дппд

(6.6)

где Tпд – постоянная времени пропорционально-дифференцирующего ре-гулятора.

Пропорционально-интегрирующие (ПИ):

,1sT

skk

sWпи

ип

пи

(6.7)

где Tпи – постоянная времени пропорционально-интегрирующего регуля-тора.

Пропорционально-интегродифференцирующие (ПИД):

,1sT

sks

kksW

пид

ди

ппид

(6.8)

где Tпид – постоянная времени пропорционально-интегро-дифференцирующего регулятора.

Применение идеальных регуляторов в реальных САУ возможно, если постоянная времени применяемого регулятора достаточно мала и можно принять, что она равна нулю, в этом случае все реальные регуляторы пре-вращаются в идеальные.

Корректирующие обратные связи делятся: − По знаку обратной связи:

90

а) на отрицательные; б) положительные.

Главные обратные связи, как правило, являются отрицательными. Ме-стные обратные связи могут иметь любой знак. − По составу: а) на жесткие; б) гибкие.

Жесткая обратная связь осуществляется статическим звеном, то есть когда Wос(0) 0. Гибкая обратная связь осуществляется дифференцирую-щим звеном и Wос(0) = 0, следовательно, в статике сигнал обратной связи отсутствует, то есть гибкая обратная связь действует только в динамике, с чем и связано ее название.

6.2. Получение передаточных функций регуляторов

в пакете Simulink

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

1. П-регулятор, представляет собой идеальное усилительное звено (блок Gain), коэффициентом усиления которого является kп (рис.6.2).

Рис.6.2. Структурная схема идеального П-регулятора в пакете Simulink

2. ПД-регулятор, представляет собой параллельное соединение идеаль-

ного усилительного звена (блок Gain) и последовательно соединенных идеального усилительного и дифференцирующего (блок Switched deriva-tive for linearization) звеньев (рис.6.3).

Рис.6.3. Структурная схема идеального ПД-регулятора в пакете Simulink

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

ального усилительного звена (блок Gain) и последовательно соединенных идеального усилительного и интегрирующего (блок Integrator) звеньев (рис.6.4).

91

Рис.6.4. Структурная схема идеального ПИ-регулятора в пакете Simulink

Знак в сумматоре зависит от знака интегрирующей составляющей. 4. ПИД-регулятор, представляет собой параллельное соединение иде-

ального усилительного звена (блок Gain), последовательно соединенных идеального усилительного и дифференцирующего (блок Switched deriva-tive for linearization) звеньев и последовательно соединенных идеального усилительного и интегрирующего (блок Integrator) звеньев (рис.6.5).

Рис.6.5. Структурная схема идеального ПИД-регулятора в пакете Simulink

Знаки в сумматоре зависят от знаков дифференцирующей и интегри-

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

димо добавить к соответствующей структурной схеме блок Тransfer Fcn. (последовательное соединение) для учета инерционности (рис.6.6 − 6.9).

Рис.6.6. Структурная схема реального

П-регулятора в пакете Simulink

92

Рис.6.7. Структурная схема реального ПД-регулятора

в пакете Simulink

Рис.6.8. Структурная схема реального ПИ-регулятора

в пакете Simulink

Рис.6.9. Структурная схема реального ПИД-регулятора

в пакете Simulink

6.3. Общая методика синтеза линейных систем автоматического управления

Как уже говорилось выше, задача синтеза САУ заключается в опреде-

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

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

93

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

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

Рассмотрим общую методику синтеза линейных САУ. 1. Определение основной (неварьируемой) части системы. Часть звеньев системы обычно однозначно определяется непосредст-

венно заданием на разработку системы. К ним относятся: объект управле-ния, смежные с объектом звенья, и определяемые ими исполнительные и чувствительные элементы управляющего устройства.

В некоторых случаях возможно варьировать исполнительные и чувст-вительные элементы (например, замена датчика одного типа на датчик другого типа), но параметры этих элементов останутся неизменными, то есть неварьируемыми. Разумеется, при выборе этих звеньев необходимо исходить из предъявляемых к САУ требований (по точности, быстродей-ствию).

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

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

2. Определение требуемого коэффициента передачи и порядка астатиз-ма системы.

Этот параметр САУ находится исходя из требований к точности в ус-тановившихся режимах при детерминированных воздействиях. В общем случае, при известном значении статической ошибки в установившемся режиме (s = 0), определяется требуемый коэффициент передачи системы. При этом, если коэффициент передачи оказывается настолько большим, что существенно затрудняет осуществление требований по устойчивости, целесообразно повысить порядок астатизма и этим свести до нуля статиче-скую ошибку системы (сделать ее не зависящей от коэффициента передачи системы). В результате чего становится возможным выбирать величину этого коэффициента исходя только из соображений устойчивости и каче-ства при нулевой статической ошибке.

Однако существуют системы, в которых статическая ошибка не должна

94

быть равна нулю (в частности, системы автоматического регулирования частоты вращения судовых дизелей – для предотвращения перегрузки), и более того, существуют САУ, в которых статическая ошибка должна быть величиной регулируемой. Такие системы не могут быть астатическими, и в случае, если в этих системах невозможно обеспечить устойчивость при требуемой статической ошибке, то необходимо изменение структуры кор-ректирующего звена либо введение дополнительных корректирующих звеньев.

На этом же этапе решается вопрос о применении управления по не-скольким или по одному выделенному (основному) возмущению.

Управление по возмущению целесообразно, если есть возможность достаточно просто измерить возмущение, по которому осуществляется управление, и если такое управление требует достаточно простых коррек-тирующих звеньев.

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

3. Выбор корректирующих звеньев и их параметров. а) Синтез по управляющему (основному) воздействию. В первую очередь осуществляется подборка корректирующего звена и

его параметров с целью обеспечения устойчивости системы и необходи-мых запасов качества. Наиболее удобным способом для этого является синтез с помощью логарифмических частотных характеристик, при кото-ром строятся требуемые частотные характеристики, из которых определя-ется передаточная функция требуемого корректирующего звена (вид кор-ректирующего звена и его параметры). В пакете Simulink такая коррекция выполняется с помощью расширения SISO Design Tool (инструмент для синтеза корректирующих звеньев САУ), работа с которым будет описана ниже.

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

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

95

При этом необходимо обеспечивать требуемый коэффициент передачи системы, полученный на предыдущем этапе синтеза, для поддержания требуемой статической ошибки, то есть в случае синтеза по частотным ха-рактеристикам изменять частоты сопряжений (на ЛАХ), вид и наклон ЛФХ, не изменяя начальной точки ЛАХ (L(1) = 20 lg k = const).

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

б) Синтез по возмущающим воздействиям. Корректирующие воздействия в функции внешних воздействий повы-

шают качество переходных процессов, воздействуя на объект управления либо в функции внешних воздействий, либо в функции их производных. Порядок производной зависит от звеньев, входящих в контур управления по возмущению, то есть как от объекта управления, так и от датчика и т.д. В результате введения в САУ управления по возмущающему воздействию получаем комбинированную САУ (рис.6.10): ЧЭ и ЧЭ’ – чувствительные элементы по управляющему и возмущающему воздействиям соответст-венно; УУ и УУ’ – управляющие устройства (корректирующие звенья) по управляющему и возмущающему воздействиям соответственно: ОУ и ОУ’ – объект управления по управляющему и возмущающему воздействиям соответственно; ИМ – исполнительный механизм.

Рис.6.10. Функциональная схема комбинированной САУ

g

УУ ОУ

ЧЭ

x

ЧЭ’ ОУ’

ИМ

f

УУ’

96

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

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

Для случая САУ, изображенного на рис.6.10:

Wxf(s) = Wоу’(s) + Wчэ’(s)Wуу’(s)Wим(s)Wоу(s). (6.9)

Идеальным является независимость выходной величины от возму-щающего воздействия, то есть

Wxf(s) = 0. (6.10)

Таким образом, приравняв нулю выражение (6.9), необходимо выра-

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

Wуу’(s) = – Wоу’(s) / (Wчэ’(s)Wим(s)Wоу(s)). (6.11)

Чаще всего полученное корректирующее звено оказывается очень

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

Wxf(0) = 0. (6.12)

При этом получаются относительно несложные корректирующие зве-

нья, которые легко реализуются на практике (чаще всего П- или ПД- регу-ляторы).

В общем случае в САУ может быть любое количество возмущающих воздействий. Коррекцию следует производить по тем из них, которые наи-более сильно влияют на объект управления, которые можно измерить и по которым возможно реализовать коррекцию.

97

6.4. Приложение для синтеза корректирующих звеньев SISO Design Tool

SISO Design Tool является приложением системы MATLAB пакета

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

Это приложение основано на применении LTI Viewer, входит в пакет Control System и запускается из окна Launch Pad (Панель запуска) основно-го окна MATLAB. Для запуска необходимо двукратное нажатие на опции Control System Toolbox и выбор (двукратным нажатием) SISO Design Tool, после чего должно появиться основное окно этого приложения.

Основное окно содержит в себе ( рис.6.11): 1 − командное меню; 2 − панель инструментов; 3 − панель корректирующего звена; 4 − панель структуры системы; 5 − корневой годограф; 6 − логарифмическую амплитудную характеристику; 7 − логарифмическую фазовую характеристику; 8 − панель подсказок.

Рис.6.11. Основное окно приложения SISO Design Tool

1

2 4

3

5

6

7

8

98

Рис.6.12. Панель структуры системы

После запуска SISO Design Tool следует сначала определить, как вы-

глядит модель, для которой составляется корректирующее звено. Она представлена на панели структуры системы (рис.6.12).

Здесь введены следующие обозначения: F – Prefilter (Префильтр) – передаточная функция совокупности звеньев, стоящих между точкой приложения входного сигнала и сумматором; C – Compensator (Корректирующее звено) – передаточная функция синте-зируемого звена; G – Plant (Объект) – передаточная функция совокупности звеньев, стоящих в прямой ветви; H – Sensor (Датчик) – передаточная функция совокупности звеньев, стоя-щих в обратной связи системы.

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

На этой же панели расположены две кнопки: − кнопка с обозначением “+/−” изменяет знак обратной связи системы; − кнопка с обозначением “FS” изменяет положение корректирующего зве-на в прямой ветви либо в обратной связи.

Таким образом, с помощью приложения SISO Design Tool возможно производить как синтез САУ, так и синтез местных отрицательных или по-ложительных обратных связей с корректирующими звеньями в них.

6.4.1. Командное меню

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

1) File − работа с файлами моделей. Import… – импорт моделей из рабочей области или с диска;

Диалоговое окно для импорта моделей (рис.6.13) содержит три поля: − System Name – имя синтезируемой системы. При изменении имени в этом поле изменяется название синтезируемой системы, что отражается в заголовке основного окна SISO Design Tool; − Import from – здесь указывается, откуда производится импорт модели: Workspace – из рабочей области; MAT-file – из МАТ-файла (сохранение части или всей рабочей области MATLAB и Simulink – из файлов моделей

99

Рис.6.13. Окно импорта моделей в SISO Design tool

Simulink, собранных из LTI-блоков. Если импорт производится с диска, то становится активным окно, в которое вводится путь и имя импортируемого файла, а также кнопка Browse, которая открывает стандартное окно для за-грузки файлов; − System Data – здесь указана структурная схема системы (рассмотрена выше) и указывается, какой именно блок импортируется. Кнопка Other… изменяет положение корректирующего звена (в прямой ветви или обрат-ной связи). Export… – экспорт моделей в рабочую область или на диск. Окно для экспорта моделей (рис.6.14) включает в себя список экспорти-руемых моделей, выполненный в виде таблицы, где: − Component – компоненты системы, которые возможно экспортировать, сюда входят части структурной схемы и модели разомкнутой (Open Loop) и замкнутой (Closed Loop) синтезируемых систем; − Model – формат экспортируемой модели (Current – текущий формат (за-дается импортом), State Space – в переменных состояния и др.); − Export As – имя, с которым будет экспортирована система. Также в этом окне справа расположены кнопки: − Export to Workspace – экспорт в рабочую область; − Export to Disk – экспорт на диск; − Help – вызов справки по экспорту и форматам моделей; − Close – закрыть окно экспорта.

100

Рис.6.14. Окно экспорта моделей из SISO Design tool Для осуществления экспорта в рабочую область или на диск необходи-

мо выбрать в таблице нужные модели (при нажатой кнопке Ctrl левой кнопкой мыши выделяются отдельные строки таблицы) и затем нажать со-ответствующую кнопку. Toolbox Preferences… – настройка приложения (полностью эквивалентна настройкам LTI Viewer). Print… – распечатка корневого годографа и графиков логарифмических частотных характеристик (на одном листе) (Ctrl+P). Print to Figure – передать графики корневого годографа и логарифмиче-ских частотных характеристик в редактор рисунков системы MATLAB. Close – закрыть приложение SISO Design Tool (Ctrl+W).

2) Edit – операции редактирования передаточной функции корректи-рующего звена (с помощью изменения графиков корневого годографа и логарифмических частотных характеристик). Undo – отменить последнюю операцию (Ctrl+Z). Redo – отменить последнюю команду Undo (Ctrl+Y). Root Locus… – редактирование (построение желаемого) корневого годо-графа. Bode… – редактирование (построение желаемых) логарифмических час-тотных характеристик.

Подменю для редактирования корневого годографа и логарифмических частотных характеристик содержат практически одинаковый набор пара-

101

метров: 1. Add… – добавить на соответствующем графике:

Real Pole – полюс, лежащий на вещественной оси (мнимая часть равна нулю); Complex Pole – комплексно-сопряженные полюса; Integrator – интегрирующее звено (W(s) = 1/s); Real Zero − нуль, лежащий на вещественной оси; Complex Pole – комплексно-сопряженные нули; Differentiator – дифференцирующее звено (W(s) = s); Lead – добавить пару полюс/ноль (W(s) = (1+T1s)/(1+T2s)), причем полюс лежит левее нуля; Lag – добавить пару ноль/полюс (W(s) = (1+T1s)/(1+T2s)), причем полюс лежит правее нуля; Notch – добавить комплексно-сопряженные пары пар нолей и полюсов (W(s) = (1+T1s+T3s2)/(1+T2s+T3s2)).

2. Delete Pole/Zero – стереть с корневого годографа или с логарифмиче-ских частотных характеристик ранее добавленные элементы.

3. Design Constraints… – (эта опция активна только для Root Locus) за-дание показателей для построения корневого годографа.

4. Show… – (эта опция активна только для Bode) настройка характери-стик, показываемых для логарифмических частотных характеристик (Am-plitude – амплитуда, Phase – фаза, Stability Margins – запасы устойчивости).

5. Grid – показать/убрать сетку. 6. Zoom – опции масштабирования (X-Y – по обоим осям, ln-X – по го-

ризонтальной оси, ln-Y – по вертикальной оси, Out – вернуть первоначаль-ный масштаб).

7. Properties – свойства соответствующего графика (адекватны на-стройкам LTI Viewer).

3) View – настройка внешнего вида SISO Design Tool. Root Locus – показать/убрать корневой годограф. Bode Diagram - показать/убрать логарифмические частотные характери-стики. System Data – просмотреть сведения по синтезируемой системе (рис.6.15).

В этом окне указывается имя синтезируемой системы (System Name) и перечисляются ноли (Zeros) и полюса (Poles) передаточных функций неиз-менных частей системы. Кнопка Show Transfer Function показывает пере-даточную функцию соответствующей части системы. Closed-Loop Poles – показывает полюса замкнутой системы. Design History – показывает все операции, произведенные над синтези-руемой системой, кнопка Save to Text File позволяет сохранить их в тек-стовом файле.

102

Рис.6.15. Окно сведений по синтезируемой системе (System Data)

4) Compensator – операции над корректирующим звеном (рис.6.16). Format… – вызывает окно установки формата корректирующего звена.

Рис.6.16. Окно установки формата корректирующего звена

103

Возможны два формата корректирующего звена (вид передаточной функции): Time-constant – с постоянными времени и Zero/pole/gain – ну-ли/полюса/общий коэффициент усиления.

Кроме того, опция Show plant/sensor poles and zeros включает или вы-ключает показ полюсов и нолей прямой ветви и корректирующего звена на графике логарифмических частотных характеристик. Edit… – редактирование (изменение) корректирующего звена (рис.6.17).

Рис.6.17. Окно редактирования корректирующего звена В этом окне вводится общий коэффициент усиления (Gain), выбирается

формат (Format) корректирующего звена и добавляются нули (Zeros) и/или полюса (Poles), лежащие на вещественной (Real) или комплексной оси (Complex). В этом же окне выводятся численные значения для уже добав-ленных в корректирующее звено компонентов (см. рис.6.17), которые можно изменять в соответствующих окнах, также напротив каждого ком-понента имеется пустой квадратик (под надписью Delete – стереть), поста-

Рис.6.18. Окно восстановления ранее сохраненных корректирующих звеньев

104

вив в котором галочку (левой кнопкой мыши) и нажав кнопку “Apple”, можно удалить этот компонент. Store – сохранить в памяти текущее корректирующее звено (запрашива-ется имя сохраняемого корректирующего звена). • Retrieve… – восстановить из памяти ранее сохраненное корректирующее звено (рис.6.18).

В появившемся окне в таблице указывается: Name – имя; Order – выс-ший порядок, Sample Time – время дискретизации для дискретного кор-ректирующего звена.

Кнопка Retrieve восстанавливает выбранное из списка корректирую-щее звено (выбрать можно только одно корректирующее звено, кнопка De-lete – удаляет из памяти, кнопка Help – вызывает справку, кнопка Cancel – отменяет восстановление). Clear – очистить текущее корректирующее звено (устанавливается иде-альное усилительное звено с коэффициентом усиления, равным единице).

3) Tools – набор инструментов для работы с SISO Design Tool Loop Responses… – вывод различных графиков синтезируемой системы с помощью LTI Viewer.

Этот пункт включает в себя следующие подпункты (рис.6.19): − Plant Output (Step) – график переходного процесса для замкнутой систе-мы от единичного ступенчатого воздействия;

Рис.6.19. Окно построения произвольных графиков

105

− Control Signal (Step) – график переходного процесса для замкнутой сис-темы от дельта-функции (единичного импульса); − Compensator Bode – логарифмические частотные характеристики коррек-тирующего звена; − Closed-Loop Bode – логарифмические частотные характеристики замкну-той системы; − Open-Loop Nyquist – амлитудно-фазовая частотная характеристика ра-зомкнутой системы; − Open-Loop Nichols – характеристика Николса для разомкнутой системы; − Custom – произвольные графики.

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

В разделе Plots (графики) показаны шесть возможных графиков, кото-рые будут расположены на одном окне. С помощью пункта Change to вы-бирается тип соответствующего графика: step – переходный процесс от единичного ступенчатого сигнала; impulse – переходный процесс от дель-та-функции (единичного импульса); bode – логарифмические частотные характеристики; bodemag – логарифмическая амплитудная характеристика; nyquist – амплитудно-фазовая частотная характеристика; nichols – диа-грамма Николса; pzmap – корневой годограф; none – нет графика.

В разделе Contents of Plot выбираются входная и выходная точки для замкнутой системы (также приводятся звенья, попадающие в контур) и для разомкнутой системы – входящие в нее звенья структурной схемы. Continuous/Discrete Conversions… – переход от непрерывной системы к дискретной и наоборот, выбирается время дискретизации (Sample time) и метод перехода для каждого звена синтезируемой схемы. Draw Simulink Diagram – составить в Simulink модель, адекватную син-тезируемой. Создает модель со структурной схемой, используемой для синтеза, каждый блок которой представлен подсистемой в переменных со-стояния.

4) Window – переходы по активным окнам системы MATLAB. 5) Help – справка по различным компонентам SISO Design Tool.

6.4.2. Панель инструментов На панель инструментов вынесены минимально необходимые функции

командного меню, они разделены на три секции вертикальными чертами. Первые шесть кнопок позволяют изменять графики (получать желаемые характеристики), причем одновременно может быть активна только одна из этих кнопок: − Default mode – перемещение ранее установленных на графиках точек (активна по умолчанию);

106

− Add real pole – добавить полюс, лежащий на вещественной оси; − Add real zero – добавить ноль, лежащий на вещественной оси; − Add complex pole – добавить комплексно-сопряженные полюса; − Add complex zero – добавить комплексно-сопряженные нули; − Delete pole/zero – удалить с графиков ненужный полюс/нуль.

Далее идут четыре кнопки для изменения масштаба графиков (активна может быть только одна из них): − Mouse zoom – увеличение выделенного с помощью мыши участка гра-фика (одновременно по горизонтальной и вертикальной осям); − Zoom in X – увеличение по горизонтальной оси; − Zoom in Y – увеличение по вертикальной оси; − Zoom out (full view) – вернуться к первоначальному масштабу.

Последняя кнопка (Context-sensitive help) предназначена для получения контекстной справки, то есть при активности этой опции, при подведении указателя мыши к какому-либо объекту появляется контекстная подсказка.

6.4.3. Работа с приложением SISO Design Tool

Если в результате анализа САУ был получен вывод о ее неработоспо-

собности с данным регулятором или с заданными его настройками, то есть сделан вывод о необходимости синтеза, то необходимо воспользоваться приложением SISO Design Tool.

Рассмотрим последовательность действий для синтеза системы автома-тического регулирования частоты (САРЧ) вращения вала судового дизеля, структурная схема которой приведена на рис.6.20.

Рис.6.20. Структурная схема САРЧ

САРЧ представляет собой САУ по ошибке, в которую входят:

1) Р – регулятор, представленный фазоимпульсным усилителем (ФИУ) и терристорным преобразователем (ТП); 2) ИМ– исполнительный механизм, представленный двигателем постоянного тока (ДПТ) и редуктором (РЕД); 3) ОУ – объект управления, который представляет собой двигатель

107

внутреннего сгорания (ДВС); 4) Д – датчик, измеряющий текущую частоту вращения, который является тахогенератором (ТГ).

Требуется получить следующие характеристики САРЧ: 1) запас устойчивости по амплитуде: не менее 20 дБ; 2) запас устойчивости по фазе: от 20 до 80 градусов; 3) статическая ошибка: от 0 до ±5%; 4) перерегулирование: не более 20%; 5) время переходного процесса: от 0,5 до 3 с.

Начальные параметры настройки регулятора: − коэффициент усиления ФИУ: всегда равен единице (не изменяется); − коэффициент усиления терристорного преобразователя kтп = 1; − постоянная времени терристорного преобразователя Tтп = 0,035.

Настройки регулятора (кроме коэффициента усиления ФИУ) могут из-меняться: коэффициент усиления терристорного преобразователя – от 0,1 до 10, а постоянная времени терристорного преобразователя – от 0,001 до 0,01.

Анализ САРЧ позволил получить следующие характеристики (при на-чальных настройках регулятора): − запас устойчивости по амплитуде 15,3 дБ; − запас устойчивости по фазе 61,7 градусов; − статическая ошибка 0%; − перерегулирование 5,77%; − время переходного процесса 1,58 с.

Таким образом, САРЧ не проходит только по запасу устойчивости по амплитуде.

6.4.4. Синтез по управляющему воздействию

Для синтеза по управляющему воздействию необходимо произвести

следующие действия: 1. Представить систему в виде, пригодном для ее синтеза в SISO Design

Tool и перенести части структурной схемы в рабочую область или на диск. Для этого необходимо выделить: 1) Префильтр F – у данной системы нет звеньев, стоящих между точ-

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

2) Прямая ветвь G – для данной системы в прямой ветви находятся ИМ и ОУ. Их следует скопировать со структурной схемы САРЧ в новое окно модели и установить точки входа и выхода, как это показано на рис.6.21.

Затем необходимо вызвать окно LTI Viewer (в окне модели – Tools > Linear Analysis) и прорешать модель прямой ветви в этом приложении (в

108

окне LTI Viewer’а – Simulink > Get Linearized Model). При этом в памяти

Рис.6.21. Структурная схема прямой ветви

образуется LTI-модель, адекватная структурной схеме прямой ветви. Те-перь необходимо перенести эту модель в рабочую область MATLAB или на диск, для этого необходимо выполнить в окне LTI Viewer’а: File > Ex-port…, в появившемся окне (рис.6.22) выделить в разделе Export List нуж-ную нам систему (для рис.6.22 – reg_upr_1) и сохранить модель прямой ветви в рабочей области (Export to Workspace) или на диске (Export to Disk). В последнем случае потребуется введение имени для сохраняемой модели.

Рис.6.22. Окно экспорта модели LTI Viewer 3) Обратная связь H – у данной системы в обратной связи стоит датчик

(ТГ), имеющий передаточную функцию, равную единице; 4) Регулятор С – для данной системы регулятор состоит из ФИУ и ТП.

Порядок действий для переноса его в рабочую область или на диск такой же, как и для прямой ветви (рис.6.23).

2. Запустить приложение SISO Design Tool и импортировать в него по-лученные в предыдущем пункте модели.

Для запуска SISO Design Tool необходимо выполнить в панели Launch

109

Рис.6.23. Структурная схема регулятора

Pad основного окна MATLAB: Control System Toolbox > SISO Design Tool.

В окне SISO Design Tool выполнить: File > Import…, в появившемся окне можно задать имя системы (System Name) и следует выбрать, откуда ведется импорт – из рабочей области (Workspace) или с диска (MAT-File). В списке выбрать нужную нам систему и нажать кнопку со стрелкой, стоящую напротив нужного нам элемента системы (для прямой ветви – выбрать reg_upr_1 и нажать кнопку со стрелкой напротив G (прямая ветвь). После этого в полях элементов системы появится имя соответствующей модели. Так же импортируется и регулятор.

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

После того как экспорт моделей всех элементов завершен, необходимо нажать кнопку “ОК”.

После выполнения вышеописанных действий можно убедиться, что внешний вид графиков изменился и передаточная функция регулятора именно такая, какая нам нужна (в данном примере Wр(s) = kр (1/(1+Ts)). Передаточная функция регулятора (корректирующего звена) приводится в панели корректирующего звена (см. рис.6.11, поз.3).

3. Добиться с помощью SISO Design Tool необходимых характеристик системы.

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

Если нажать любой кнопкой мыши на неизменяемую часть системы (блоки, обозначенные желтым цветом), то можно увидеть окно сведений по синтезируемой системе (см. рис.6.15), вызываемое активизацией в окне SISO Design Tool позиции View > System Data. При нажатии на блоке кор-ректирующего звена появляется окно редактирования корректирующего звена (см. рис.6.17), вызываемое активизацией в окне SISO Design Tool по-зиции Compensator > Edit.

На логарифмических частотных характеристиках выводятся значения: G.M. – запас устойчивости по амплитуде; P.M. – запас устойчивости по фа-

110

зе; Freq. – значение частоты в соответствующей точке; Stable/Unstable loop – устойчивая/неустойчивая система.

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

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

Нажимая левой кнопкой мыши на характерные точки графиков, можно вызвать окна сведений по синтезируемой системе (см. рис.6.15) и редакти-рования корректирующего звена (см. рис.6.17). Нажатие правой кнопки мыши на поле графиков (не на линии) вызывает появление контекстного меню, полученного сведением пунктов командного меню Edit > Root Lo-cus… или Bode… и Compensator > Edit…, описанных выше. Вид контекст-ного меню зависит от того, на каком из графиков нажать правую кнопку мыши.

Изменяя положение логарифмической амплитудной характеристики вверх/вниз, можно добиться требуемого запаса устойчивости по амплиту-де. Для синтезируемой в примере системы, чтобы добиться запаса устой-чивости по амплитуде более 20 дБ, необходимо переместить логарифмиче-скую амплитудную характеристику вниз. Для примера, уменьшим коэф-фициент усиления корректирующего звена до значения 0,5 (рис.6.24) в па-нели или окне корректирующего звена. В этом случае запас устойчивости по амплитуде становится равным L = 21,3 дБ, а по фазе = 75,2 градуса, то есть показатели устойчивости удовлетворяют их требуемым значениям.

Для того, чтобы оценить показатели качества системы при этих на-стройках, необходимо получить график переходного процесса, выполнив в командном меню SISO Design Tool: Tools > Loop Responses… > Plant Out-put (Step). При этом появится окно с графиком переходного процесса (рис.6.25).

Полученные показатели качества: - статическая ошибка ст = 0%; - перерегулирование = 0%; - время переходного процесса tp = 42 с.

111

Рис.6.24. Пример настройки корректирующего звена

Рис.6.25. Переходный процесс замкнутой синтезируемой системы

Таким образом, с настройками регулятора kтп = 0,5 и Tтп = 0,035 с

112

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

Для примера, при настройках регулятора kтп = 0,85 и Tтп = 0,00145: - запас устойчивости по амплитуде L = 20,1 дБ; - запас устойчивости по фазе = 68,6 градусов; - статическая ошибка ст = 0%; - перерегулирование = 0,944%; - время переходного процесса tp = 1,21 с.

Можно попробовать еще улучшить показатели качества, но следует не забывать, что настройки регулятора лежат в пределах: коэффициент уси-ления терристорного преобразователя – от 0,1 до 10, а постоянная времени терристорного преобразователя – от 0,001 до 0,01.

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

6.4.5. Синтез по возмущающему воздействию

Синтез по возмущающему воздействию будет заключаться в добавке в

структурную схему синтезируемой системы управления по возмущению. Структурная схема САРЧ, с учетом возмущающего воздействия (измене-ния момента нагрузки на валу ДВС), но без управления по возмущению представлена на рис.6.26, где ОУ’ – передаточная функция объекта управ-ления по возмущающему воздействию.

По графику переходного процесса по возмущающему воздействию (рис.6.27) определены следующие показатели влияния возмущения на час-тоту вращения ДВС: - максимальный заброс частоты вращения ymax = -0,0377; - статическая ошибка по возмущению в

ст = 0%; - время переходного процесса по возмущению tв

пп = 1,37 с. Для определения изменения момента нагрузки можно воспользоваться

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

Wдм(s) = kдм = 1. (6.13)

113

Рис.6.26. Структурная схема САРЧ без управления по возмущению

Согласно (6.11), передаточная функция управляющего устройства (ре-

гулятора) по возмущению

Wуу’(s) = 0,002951s2 + 0,0227s. (6.14)

Рис.6.27. Переходный процесс по возмущающему воздействию

без управления по возмущению

114

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

Рис.6.28. Структурная схема САРЧ с управлением по возмущению Причем регулятор по возмущению РВ для удобства собран в блоке-

подмодели (рис.6.29).

Рис.6.29. Структурная схема регулятора по возмущению

Переходный процесс в системе по возмущающему воздействию с

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

115

Рис.6.30. Переходный процесс по возмущающему

воздействию с управлением по возмущению Небольшой скачок (ymax = −8,1110-5) в начале переходного процесса

обусловлен неточностью воспроизведения пакетом Simulink дифференци-рующего звена.

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

Wуу’(s) = 0,0227s. (6.15) Структурная схема регулятора по возмущению в этом случае будет вы-

глядеть так, как это показано на рис.6.31, и переходный процесс по возму-щающему воздействию, как это показано на рис.6.32.

Рис.6.31. Структурная схема регулятора по возмущению

116

Рис.6.32. Переходный процесс от возмущающего воздействия, с управлением по возмущению

В этом случае возмущающее воздействие влияет на частоту вращения

ДВС следующим образом: - максимальный заброс частоты вращения ymax = −0,0177; - статическая ошибка по возмущению в

ст = 0%; - время переходного процесса по возмущению tв

пп = 1,64 с. Таким образом, несмотря на то, что увеличилось время переходного

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

117

7. АНАЛИЗ ПРОСТЕЙШИХ СТРУКТУР ПРИ СЛУЧАЙНЫХ ВОЗМУЩАЮЩИХ ВОЗДЕЙСТВИЯХ

7.1. Понятие случайной величины и случайного процесса.

Законы распределения случайных величин

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

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

В системе MATLAB случайные процессы представлены двумя блока-ми: Uniform Random Number - источник случайного сигнала с равномер-ным распределением и Random Number - источник случайного сигнала с нормальным распределением. Оба блока находятся в разделе Sources биб-лиотеки Simulink.

Блок Uniform Random Number описывается следующим выражением:

если М < x < S;

,0

,MS

1)x(f

в противном случае, где М − минимальное значение уровня сигнала; S – максимальное значение уровня сигнала.

(7.1)

118

В параметрах блока Uniform Random Number (рис. 7.1) указываются значения минимального и максимального уровней сигнала (Minimum и Maximum), время дискретизации (Sample time), стартовое значение для ге-нератора случайных чисел (Initial seed).

Рис.7.1. Блок Uniform Random Number

Рис.7.2. Блок Random Number

119

Блок Random Number описывается аналитическим выражением:

2S22)Mx(

eS2

1)x(f

, (7.2)

где М − математическое ожидание; S2 − дисперсия.

В параметрах блока Random Number (рис.7.2) указывается математиче-ское ожидание (Mean), дисперсия (Variance), время дискретизации (Sample time), стартовое значение для генератора случайных чисел (Initial seed). Стартовое значение определяет вариант реализации выборки случайных чисел, для одного и того же числа в этой опции получается одна и та же выборка.

Генератор белого шума - Band Limited White Noise - служит для созда-ния шумового сигнала с заданной мощностью, равномерно распределен-ной по частоте.

7.2. Прохождение случайного сигнала через линейную систему

Для примера возьмем систему автоматического регулирования частоты

(САРЧ) вращения вала дизеля без управления по возмущению (см. рис.6.26) и в качестве возмущающего воздействия используем нормально распределенный сигнал (блок Random Number) с параметрами: математи-ческое ожидание Mean = 0,5; дисперсия Variance = 0,01; стартовое значе-ние Initial seed = 0 и время дискретизации Sample time = 0 (непрерывный сигнал). Кроме того, необходимо добавить регистрирующие блоки, пере-дающие информацию в рабочую область (блок To Workspace). Структур-ная схема показана на рис. 7.3.

Рис.7.3. Структурная схема системы стабилизации угловой скорости турбогенератора

со случайным возмущающим воздействием

120

Блок To Workspace служит для регистрации случайного сигнала, пода-ваемого на вход системы, а блок To Workspace1 – для регистрации выход-ного сигнала. Блок To Workspace имеет следующие параметры (рис.7.4):

1) Variable name – имя переменной, в которой будут сохраняться ре-зультаты;

2) Limit data points to last – предел количества точек запоминаемой ве-личины (inf – бесконечное число точек);

3) Decimation – разрядность хранимой величины; 4) Sample time – шаг времени, через который происходит сохранение

значений; 5) Save format – формат, в котором происходит сохранение, может

быть: Structure with time – структура с временем; Structure – структура, Ar-ray – массив.

Рис.7.4. Настройка блока To Workspace В параметрах обоих блоков устанавливаем: Sample time = 0,1 и Save

format – Array, параметры Limit data points to last и Decimation оставляем без изменений. Имена переменных: для входного блока – In, а для выход-ного – Out, при этом имена переменных появляются на изображении блока на схеме (см. рис.7.3).

Время моделирования устанавливаем (в меню окна модели Simulation > Simulation Parameters… > Stop time), равным 50 с, этого достаточно, чтобы определить характеристики случайного сигнала.

Далее необходимо запустить процесс моделирования (в меню окна мо-дели Simulation > Start) и выполнить следующие действия:

121

1. В основном окне MATLAB выполнить следующие команды: t=0:0.1:50; – задаем вектор времени (от 0 до 50 с, с шагом 0,1); plot(t,In) – построить график случайного процесса, подаваемого

в качестве возмущающего воздействия. Оператор “;” стоящий после команды, позволяет отключить вывод ре-

зультатов этой команды на экран. 2. После выполнения этих команд появится окно Figure No.1 (рис.7.5),

в котором и отображается график случайного сигнала. Для того чтобы по-лучить характеристики этого сигнала, необходимо из меню окна Figure No.1 выполнить: Tools > Data Statistics, после чего появится новое окно - Data Statistics-1, в котором отображаются интересующие нас статистиче-ские параметры отдельно по осям графика (X – ось времени, Y – ось орди-нат случайного сигнала):

min – минимальное значение; max – максимальное значение; mean – математическое ожидание; std – среднеквадратическое отклонение (квадратный корень из дис-

персии, значение которой задано в параметрах блока Random Num-ber).

Рис.7.5. График случайного процесса и его параметры на входе в САРЧ При нажатии на квадратик напротив любой из величин появляются ли-

122

нии на графике, отображающие эту величину. Нас интересуют параметры случайного сигнала, поданного на вход

системы, то есть ось Y. Можно убедиться, что математическое ожидание близко к 0,5 и составляет 0,4973 (при увеличении количества точек оно стремится к заданному, равному 0,5) и среднеквадратическое отклонение, равное квадратному корню из дисперсии, составляет 0,1 ( 1,001,0 ).

Из-за особенностей пакета Data Statistics необходимо проверять значе-ние среднеквадратического отклонения, для этого в основном окне MAT-LAB выполняем команду std(In,1), при этом выдается значение 0,9999, что очень близко к заданной его величине.

3. Для анализа выходного сигнала необходимо выполнить следующие команды:

t1=3:0.1:50; – задаем вектор времени (от 3 до 50 с, с шагом 0,1), так как переходный процесс по возмущению составляет 1,37 с, то анализируем выход без влияния переходного процесса по управлению и берем не менее чем двукратный запас по времени;

Out1=Out(30:500); – вычленяем из выходного массива выходной сиг-нал, без влияния переходного процесса по управлению, то есть отбрасыва-ем 30 первых точек (первый операнд в скобках − это начальная точка, рав-ная начальному времени из предыдущей команды, деленному на шаг вре-мени, второй операнд – конечная точка, равная времени моделирования, деленному на шаг времени);

plot(t1,Out1) – построить график случайного процесса на выходе из САРЧ;

std(Out1,1) – уточняем значение среднеквадратического значения. 4. В окне Figure No.1 появляется график случайного процесса на выхо-

де из САРЧ, а в окне Data Statistics – характеристики этого процесса, без влияния переходного процесса по управлению (рис.7.6).

Получены следующие характеристики: минимальное значение равно 0,9978; максимальное значение равно 1,002; математическое ожидание равно 1; среднеквадратическое отклонение равно 6,755610-4 (уточненное

значение); дисперсия равна (6,755610-4)2 = 4,563810-7. Таким образом, в качестве возмущающего сигнала использовался слу-

чайный сигнал с параметрами: математическое ожидание, равное 0,5, то есть момент нагрузки на валу двигателя в среднем равен половине от мак-симально возможного; минимальный момент нагрузки равен 0,1833 от максимально возможного и максимальный момент нагрузки равен 0,7802 от максимально возможного.

123

Рис.7.6. График случайного процесса и его параметры на входе в САРЧ

При прохождении возмущающего сигнала через систему он изменяется и частота вращения от возмущающего воздействия в среднем не изменя-ется, так как математическое ожидание случайного процесса на выходе равно единице. Но в целом частота вращения изменяется от минимального значения, равного 0,9978, до максимального значения, равного 1,002, от номинальной частоты вращения.

Дисперсия случайного сигнала уменьшилась с 0,1 до 4,563810-7, то есть в 219115 раз. Таким образом, САРЧ служит своеобразным фильтром для случайного сигнала, проходящего через нее.

Для получения плотности вероятности распределения случайного сиг-нала служит команда histfit(Variable name, N), где операнд Variable name указывает на имя переменной, а операнд N устанавливает количество ин-тервалов гистограммы.

Для нашего примера значения этих операндов: 1) для входного сигнала – имя переменной In, количество интервалов

гистограммы, равное времени моделирования, деленному на шаг (500). Команда выглядит следующим образом: histfit(In,500) (рис.7.7);

2) для выходного сигнала – имя переменной Out1 (без переходного процесса по управлению), количество интервалов гистограммы, равное времени моделирования, деленному на шаг, за вычетом отбрасываемых

124

точек (470). Команда выглядит следующим образом: histfit(Out1,470) (рис. 7.8).

Рис.7.7. Плотность вероятности распределения случайного

сигнала на входе в САРЧ

Рис.7.8. Плотность вероятности распределения случайного

сигнала на выходе из САРЧ

125

8. НЕЛИНЕЙНЫЕ БЛОКИ ПАКЕТА Simulink

Пакет Simulink предназначен для моделирования любых систем, в том числе и нелинейных динамических систем. Раздел Nonlinear основной библиотеки Simulink, посвященный нелинейным компонентам, содержит наиболее распространенные нелинейные блоки (рис. 8.1), которые в теории автоматического управления имеют название «типовые нелинейные зве-нья».

Рис.8.1. Раздел нелинейных блоков библиотеки Simulink

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

126

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

8.1. Блок ограничения Saturation

Блок Saturation (звено с ограничением или насыщением) представляет

собой нелинейное устройство, сигнал на выходе которого равен входному сигналу до тех пор, пока не достигает порогов ограничения: верхнего Up-per limit или нижнего Lower limit, после этого сигнал перестает изменяться.

Статическая характеристика этого блока приведена на рис. 8.2.

Рис.8.2. Статическая характеристика блока с ограни-чением Saturation

Уравнение, описывающее этот блок, имеет вид:

,bXприB

,bXприB,bXприkX

Y (8.1)

где k = tg . Вышесказанное хорошо прослеживается при ограничении синусои-

дального сигнала (рис. 8.3). Для удобства сравнения сигналов – исходного и прошедшего через не-

линейный блок – используем блок Mux (мультиплексор) из раздела Signal & Systems библиотеки Simulink. На входы блока Mux (на рис. 8.3 он обо-значен черным прямоугольником) подаются соответствующие сигналы, а выход присоединяется к блоку Scope.

Как видно из рис. 8.3, окно параметров блока содержит лишь поля для установки верхнего (Upper limit) и нижнего (Lower limit) порогов ограни-чения.

B

-B

0 b -b X

Y

127

Рис.8.3. Блок ограничения Saturation и пример его применения

8.2. Блок с зоной нечувствительности Dead Zone

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

Статическая характеристика этого блока приведена на рис. 8.4.

Рис.8.4. Статическая характеристика блока c зоной нечувствительности Dead Zone

Уравнение, описывающее этот блок, имеет вид:

.bXпри)bX(k,bXпри0,bXпри)bX(k

Y (8.2)

Эта нелинейность моделируется блоком Dead Zone (рис. 8.5).

0 b

-b X

Y

128

Рис.8.5. Блок с зоной нечувствительности Dead Zone и пример его применения

Окно параметров этого блока содержит границы зоны нечувствитель-ности Start of dead zone и End of dead zone. Пo умолчанию они заданы рав-ными – 0,5 и 0,5. Флажки Saturate on integer owerflow (ограничение при пе-реполнении целых) и Treat as gain when linearizing (трактовать как ограни-чения при линеаризации) по умолчанию включены.

8.3. Релейный блок Relay

Релейный блок Relay имеет разрывную передаточную функцию с гис-

терезисом (или без него), подобную передаточной функции хорошо из-вестного триггера Шмитта или типовому нелинейному звену, изучаемому в курсе ТАУ – двухпозиционное реле с зоной нечувствительности.

Для этого звена, если сигнал на входе меньше некоторого порога, то на выходе получается сигнал одного уровня (обычно низкого), а если порог превышен, то сигнал на выходе становится другого уровня (обычно высо-кого). Если при спаде сигнала достигается другой порог, то сигнал на вы-ходе также скачком меняется.

Статическая характеристика этого блока приведена на рис. 8.6.

129

Рис.8.6. Статическая характеристика релейного блока Relay

Уравнение, описывающее этот блок, имеет вид:

.0Xесли,,bXприB,bXприB

;0Xесли,,bXприB,bXприB

Y (8.3)

Пример, приведенный на рис. 8.7, показывает работу релейного блока с

одинаковыми по абсолютной величине и очень малыми (eps) порогами при подаче синусоидального сигнала на вход.

Рис.8.7. Релейный блок Relay и пример его применения В окне параметров блока можно задать уровни сигнала на выходе при

включенном (Output when on) и выключенном (Output when on) состояни-ях, а также верхний и нижний пороговые уровни срабатывания (Switch on point и Switch off point соответственно). Их значения по умолчанию пред-ставлены на рис. 8.7.

0 b -b X

Y B

-B

130

8.4. Блок с ограничением скорости Rate Limiter

Блок с ограничением скорости Rate Limiter старается отслеживать

входной сигнал в условиях задания ограничений на скорость нарастания и спада сигнала (рис. 8.8).

Рис.8.8. Блок с ограничением скорости Rate Limiter и пример его применения

При работе блока вычисленное значение скорости изменения сигнала

сравнивается с установленным в окне параметров значением параметра Rising slew rate (скорость нарастания сигнала) и параметра Falling slew rate (скорость спада сигнала). По умолчанию заданы значения Rising slew rate и Falling slew rate, равные 0,5 и – 0,5 соответственно.

8.5. Блок следящего квантования Quantizer

Блок Quantizer служит для квантования меняющихся сигналов с одина-

ковым шагом по уровню (рис. 8.9). Блок имеет единственный параметр - шаг по уровню (по умолчанию

0,5). На рис. 8.9 показано квантование синусоидального сигнала. Можно отметить, что при большом шаге квантования его трудно на-

звать идеальным – точного слежения за уровнем входного сигнала нет.

131

Рис.8.9. Блок следящего квантования Quantizer и пример его применения

8.6. Блок фрикционных эффектов Coulombic and Viscous Friction

Блок фрикционных эффектов Coulombic and Viscous Friction служит

для моделирования фрикционных эффектов (рис. 8.10). Передаточная функция блока указана в окне установки его параметров.

В качестве параметра блока задается список смещений при фрикцион-ных эффектах и коэффициент передачи для приращений выходного сигнала.

Рис.8.10. Блок фрикционных эффектов Coulombic and Viscous Friction и пример его применения

132

8.7. Блок люфта Backlash

Блок Backlash имитирует эффект возникновения люфта (рис.8.12) и из-вестен в ТАУ как звено с мертвым ходом. Этот эффект создает передаточ-ную характеристику гистерезисного типа, которая представляется графи-чески в пиктограмме блока.

Статическая характеристика этого блока приведена на рис. 8.11.

Рис.8.11. Статическая характеристика блока люфта

Backlash

).XsignbX(kYили0Xпри)bX(k,0Xпри)bX(k

Y

(8.4)

Блок имеет два параметра: ширину диапазона Deadband width и на-чальный уровень сигнала на выходе Initial output (по умолчанию 1 и 0). Уровень Initial output является также срединным значением входного сиг-нала, a Deadband width определяет ширину петли гистерезиса передаточ-ной характеристики блока.

Рис.8.12. Блок люфта Backlash и пример его применения

0 b -b X

Y

133

8.8. Блоки-переключатели Switch

В основной библиотеке Simulink представлены три блока-переключателя:

1) Manual Switch (ручной переключатель) – осуществляет переключе-ние между двумя входами при двукратном нажатии на изображении блока (рис.8.13);

Рис.8.13. Блок Manual Switch и пример его применения

2) Switch (переключатель) – у этого блока переключение осуществля-ется также между двумя входами (первым и третим), причем сигнал про-ходит через первый вход, если значение на управляющем входе (второй вход) больше или равно значению, заданному в параметрах (Parameters) данного блока. В противном случае сигнал проходит через третий вход (рис.8.14);

3) Multiport switch (многовходовый переключатель) – осуществляет переключение между любым количеством входов (их количество задается в параметрах блока Number of inputs). Управляющим является первый вход, на который должен подаваться номер нужного входа (рис.8.15), при-чем значение на управляющем входе не должно быть меньше единицы; округление значения на управляющем входе не производится, то есть сиг-

134

нал идет через первый вход, если значение на управляющем входе лежит в пределах 1 xупр < 2, через второй вход, если 2 xупр < 3, и т.д.

Рис.8.14. Блок Switch и пример его применения

Рис.8.15. Блок Multiport Switch и пример его применения

135

9. ПАКЕТ NONLINEAR CONTROL DESIGN (NCD) BLOCKSET

Пакет прикладных программ для построения нелинейных систем управления Nonlinear Control Design (NCD) Blockset реализует метод ди-намической оптимизации. Это приложение представляет собой набор бло-ков, разработанных для использования в пакете Simulink и автоматически настраивающих параметры моделируемых систем, основываясь на опреде-ленных пользователем ограничениях на их временных характеристиках.

Пакет использует метод Click And Drag («щелкни и тяни») для измене-ния временных ограничений. Он реализует следующие возможности:

- легкую настройку переменных; - указание неопределенных параметров систем; - интерактивную оптимизацию; - моделирование методом Монте-Карло; - поддержка проектирования как одномерных, так и многомерных систем

управления; - моделирование подавления помех; - моделирование процессов слежения; - моделирование объектов с запаздыванием; - решение других задач управления.

Средства пакета являются весьма полезным дополнением к средствам моделирования систем управления.

9.1. Основные правила работы с NCD Пакет расширения NCD является частью пакета Simulink и наследует

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

На рис. 9.1 показано окно браузера библиотек, демонстрирующее дос-туп к библиотекам пакета NCD. Сокращение RMS означает сред-неквадратическую величину (root mean square).

Пакет NCD Blockset содержит следующие блоки (рис. 9.2): – блок CRMS (Continuous RMS); – блок DRMS (Discrete RMS); – блок NCD Outport (см. рис. 9.1).

136

Рис.9.1. Окно браузера библиотек

Рис.9.2. Раздел библиотеки NCD Blockset Рассматриваемые блоки могут применяться, в частности, в системах

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

137

Блок NCD Output является основным в рассматриваемом наборе бло-ков. Он имеет свое рабочее окно и меню и позволяет в интерактивном ре-жиме выполнять следующие операции:

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

- указывать параметры, подлежащие оптимизации; - указывать неопределенные параметры; - проводить параметрическую оптимизацию системы с учетом заданных

ограничений. Типовой сеанс работы в среде Simulink с использованием возможно-

стей и блоков NCD Blockset состоит из ряда стадий, описанных ниже. 1. В среде Simulink создается модель исследуемой динамической сис-

темы (в общем случае нелинейной). 2. Входы блоков NCD Outport соединяются с теми сигналами системы,

на которые накладываются ограничения. Этими сигналами могут быть, на-пример, выходы системы, их среднеквадратические отклонения и т. д.

3. В режиме командной строки MATLAB задаются начальные значения параметров, подлежащих оптимизации.

4. Двойным щелчком на пиктограмме NCD Output данные блоки «рас-крываются».

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

6. С помощью меню блока NCD Output задается интервал дискретиза-ции (один или два процента от длительности процесса моделирования) и указываются имена (идентификаторы) параметров системы, подлежащих оптимизации.

7. Задаются неопределенные параметры системы, указываются их но-минальные значения.

8. При необходимости сформированные ограничения сохраняются в виде файла с помощью команды меню Save (позднее они, естественно, мо-гут быть загружены с помощью команды Load).

9. Процесс оптимизации системы инициализируется нажатием кнопки Start.

9.2. Оптимизация коэффициента передачи И-регулятора

Рассмотрим пример параметрической оптимизации замкнутой системы автоматического регулирования (рис. 9.3).

Данная система представляет собой замкнутую структуру, состоящую: 1) из объекта регулирования, который в модели отображен звеном 2-го

порядка (блок Transfer Fcn) и звеном транспортного запаздывания (блок Transport Delay);

138

Рис.9.3. Пример использования системы NCD

2) интегрального регулятора (И-регулятора), отображаемого последо-вательно соединенными пропорциональным звеном с коэффициентом про-порциональности Kint и интегрирующим звеном с ограничением на вы-ходной сигнал (блок Limited Integrator);

3) контура обратной связи и звена сравнения Sum. В модель также введены источник входного сигнала в виде единичного

скачка и NCD-блок типа NCD Output, подключенный к выходу системы (здесь данный блок имеет имя NCD Outport 1).

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

- максимальное перерегулирование - не более 10 %; - время нарастания - не более 10 с; - длительность переходного процесса - не более 30 с. В командной строке MATLAB наберем: Kint=0.3 Далее двойным щелчком мыши откроем рабочее окно блока NCD

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

Заметим, что точную установку линий ограничения можно провести, выбирая требуемые линии с помощью щелчка левой кнопкой мыши (вы-бранная линия изменяет цвет - становится белой) и используя затем ко-манду меню Edit > Edit constraint. В появляющемся окне редактора ограни-чений (Constraint Editor) в строке снизу необходимо задать начальную и конечную точки прямой в формате [x1 y1 х2 у2], нажав затем кнопку Done.

Точно такие же действия можно обеспечить, не используя пункты ме-

139

ню, а просто щелкнув на соответствующей линии правой кнопкой мыши. Следующий этап - указание переменных, подлежащих оптимизации.

Выбор команды меню Optimization>Parameters приведет к открытию диа-логового окна задания настраиваемых параметров и интервала дискретиза-ции - Optimization Parameters (рис. 9.4).

Рис.9.4. Настройка блока NCD Outport

Рис.9.5. Результаты оптимизации

140

В верхнем поле данного окна (Tunable Variables) введем имя настраи-ваемого параметра Kint (если таких параметров несколько, то их имена разделяются пробелами или запятыми), а величину интервала дискретиза-ции (Discretization interval) установим равной 0,5 (ввод завершается нажа-тием кнопки Done).

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

Выведенные сообщения носят в основном служебный характер (глав-ное из них говорит о том, что процесс поиска решения завершился успеш-но - Optimization Converged Successfully). Что же касается формы переход-ного процесса для оптимизированной системы, то из рис. 9.5 видно, что она полностью соответствует заданным ограничениям.

Найденную оптимальную величину параметра Kint можно узнать, на-брав в командной строке имя данного параметра.

Удаление графиков процессов в окне блока NCD Output производится с помощью команды меню Edit > Delete plots (или одновременным нажатием клавиш Ctrl+X).

9.3. Меню окна блока NCD Output

Меню File (Файл) содержит стандартные команды Load (Загрузить),

Close (Закрыть), Save (Сохранить) и Print (Печатать). Действия, выполняе-мые при выборе любой из данных команд, относятся к графикам заданных временных ограничений, отображаемым в основном окне рассматриваемо-го блока. Меню Edit (Правка) содержит уже рассмотренные команды Edit constraint и Delete plots, а также команду Undo (Отменить), производящую отмену последнего действия.

Меню Options (Настройка) содержит следующие команды: - Initial response (Начальный отклик). Выбор данной команды приводит

к выводу (в основном окне блока NCD Output) отклика исследуемой сис-темы при начальных значениях ее параметров;

- Reference input (Задающий вход). Выбор данной команды приводит к открытию диалогового окна, в котором при желании - для вывода соответ-ствующего графика - можно указать параметры входного задающего сиг-нала системы. Никакого влияния на процессы моделирования и оптимиза-ции данные параметры не оказывают;

- Step response (Характеристики переходного процесса). Выбор этой команды приводит к открытию диалогового окна, в котором задаются па-раметры переходного процесса, такие как его длительность (Settling time), время нарастания (Rise time), максимальное перерегулирование (Percent overshoot), максимальное «недорегулирование» (Percent undershoot), соот-

141

ветствующие уровни определения данных характеристик (Percent settling и Percent rise, в процентах), а также начальное и конечное время моделиро-вания (Step time и Final time), желаемые начальное и конечное значения выхода (Initial output и Final output). Все эти параметры пользователь мо-жет изменять;

- Time range (Временной диапазон). Выбор данной команды приводит к открытию диалогового окна, в котором можно задать или изменить диапа-зон времени моделирования и метку оси времени, то есть параметры оси абсцисс;

- Y-Axis (Ось Y). То же, что и для предыдущей команды, но по отно-шению к оси ординат;

- Refresh (Обновить) - перерисовать все временные ограничения. Меню Optimization (Оптимизация) содержит рассмотренные ранее ко-

манды Parameters и Uncertainty, а также команду Start (Старт), выбор кото-рой запускает процесс моделирования и оптимизации системы (это анало-гично нажатию кнопки Start в окне Simulink или кнопки Start в панели кнопок в нижней части окна блока NCD Output). Команда Stop (Стоп) ос-танавливает процесс моделирования (аналогичное действие выполняется при нажатии кнопки Stop в панели кнопок блока NCD Output).

Наконец, последним является меню Style (Стиль). Здесь имеются сле-дующие команды:

- Grid (Сетка) - установка сетки на графике заданных ограничений; - Snap (Привязка) - при установке данного режима линии временных

ограничений можно проводить не под любым углом к оси абсцисс, а толь-ко под углом, кратным 22,5°;

- Hot-key help (Горячие клавиши) - вывод информации о «горячих» клавишах и их комбинациях;

- Readme.m (Файл справки) - вывод файла справки об окне блока NCD Output.

Из пяти кнопок панели, расположенной в нижней части окна рассмат-риваемого блока, четыре (Start, Stop, Help и Close) в дополнительных пояс-нениях не нуждаются. Что же касается пятой кнопки Split (Расщепить), то ее нажатие, при предварительном выборе с помощью левой кнопки мыши какой-либо ограничивающей линии, приводит к «расщеплению» этой ли-нии на две одинаковые по длине половинки с возможностью последующе-го редактирования отдельно каждой из них.

9.4. Особенности решаемых задач

При решении различных оптимизационных задач с помощью пакета

NCD Blockset следует иметь в виду следующие особенности этого пакета: - пакет может использоваться для решения задач оптимизации пере-

142

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

- оптимизация следящих систем проводится при входном сигнале типа единичного скачка;

- могут рассматриваться задачи оптимизации многомерных объектов с одновременным заданием временных ограничений на ряд сигналов систе-мы (что требует использования нескольких блоков NCD Output);

- пакет может также использоваться для оптимизации размещения ну-лей и полюсов передаточной функции на комплексной плоскости. Для это-го надо просто использовать блок Simulink Zero-pole или Transfer Fсn биб-лиотеки Continuous и объявить (с помощью диалогового окна настраивае-мых параметров) требуемые нули и полюса настраиваемыми, с заданием, разумеется, их предельных (наибольших и наименьших) значений.

143

10. МОДЕЛИРОВАНИЕ ЭЛЕКТРОТЕХНИЧЕСКИХ СИСТЕМ

В состав системы MATLAB входит пакет расширения Power System Blockset, позволяющий моделировать электротехнические и энергетиче-ские системы и устройства. Этот пакет ориентирован на моделирование технических устройств и систем вполне конкретного назначения и полно-стью совместим с пакетом Simulink.

Возможности пакета Power System Blockset прежде всего определяются компонентами входящих в него библиотек. Доступ к ним обычно осущест-вляется из среды Simulink (рис. 10.1). В окне браузера библиотек Simulink можно выбрать библиотеку Power System Blockset.

Рис. 10.1. Доступ к библиотекам Pоwer System Blockset

144

Можно также вызвать окно с разделами библиотеки, выполнив коман-ду powerlib в командной строке MATLAB.

10.1. Состав библиотеки Pоwer System Blockset

Рис. 10.2. Состав библиотеки Pоwer System Blockset

В состав библиотеки Pоwer System Blockset входят библиотеки сле-дующего направления: Electrical Sources – источники электрической энергии и сигналов; Elements – линейные и нелинейные компоненты электротехнических и электронных устройств; Power Electronics – устройства энергетической электроники; Machines – электрические машины; Connectors –подключающие устройства; Measurements – измерительные и контрольные устройства; Powerlib Extras – специальные энергетические устройства; Powergui – графический интерфейс пользователя пакета моделирования энергетических систем.

10.1.1. Состав библиотеки Electrical Sources

Пакет Pоwer System Blockset имеет модели источников, позволяющих

имитировать реальные источники электроэнергии. Дважды щелкнув мы-шью по пиктограмме библиотеки Electrical Sources, можно открыть окно этой библиотеки (рис. 10.3).

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

145

Рис. 10.3. Библиотека источников Electrical Sources

В составе библиотеки представлено пять типов источников электриче-

ской энергии: AC Current Source – источник переменного тока; AC Voltage Source – источник переменного напряжения; DC Voltage Source – источник постоянного напряжения; Controlled Current Source – регулируемый источник тока; Controlled Voltage Source – регулируемый источник напряжения.

10.1.2. Соединительные элементы Connectors

Состав соединительных элементов представлен на рис. 10.4.

Рис. 10.4. Библиотека соединительных элементов

Большинство соединительных элементов имеют вполне очевидное назначение. Например, заземляющие элементы служат для подключения

146

других элементов к общему проводу – земле (блок Ground). Эти элементы не имеют параметров.

10.1.3. Библиотека компонентов Elements

Основная библиотека компонентов (рис. 10.5) содержит ряд моделей,

имеющих достаточно универсальный характер. С помощью одной модели можно создать модели нескольких компонентов.

Рис. 10.5. Окно библиотеки компонентов Elements

Эта библиотека содержит несколько характерных компонентов: Series RLC Branch – последовательная RLC-цепь; Series RLC Load – последовательная RLC-цепь с нагрузкой; Parallel RLC Branch – параллельная RLC-цепь; Parallel RLC Load – параллельная RLC-цепь с нагрузкой; Linear Transformer – линейный трансформатор; Saturable Transformer – нелинейный трансформатор; Mutual Inductance – блок взаимной индуктивности; Surge Arrester – ограничитель пиковых напряжений; Breaker – выключатель управляемый;

147

PI Section Line – линия с сосредоточенными параметрами; Distributed Parameters Line – линия с распределенными параметрами.

В состав библиотеки входят две последовательные и две параллельные RLC-цепи. Эти цепи задаются тремя параметрами : сопротивлением R, идуктивностью L и емкостью C . У так называемых нагрузочных цепей до-полнительно задаются допустимые мощности рассеяния: активная для ре-зисторов и реактивные для индуктивности и конденсаторов.

Для ввода отдельных элементов (сопротивление R, идуктивность L и емкость C) можно использовать любую из RLC-цепей, задав параметрам значения, соответствующие отсутствию ненужных элементов. Например, если с помощью последовательной RLC-цепи нужно задать только рези-стор R, то надо задать L = 0 (индуктивность при этом исчезнет) и С = inf (inf означает бесконечное значение емкости, что превращает ее в провод-ник). Это правило модификации распространяется и на другие сложные компоненты.

10.1.4. Состав библиотеки Power Electronics

Библиотека рассматриваемого пакета включает достаточно представи-

тельный набор блоков коммутирующих устройств. Активизация пикто-граммы Power Electronics открывает окно с пиктограммами моделей управляемых ключей (рис. 10.6).

Рис. 10.6.Состав библиотеки компонентов

энергетической электроники

148

Представлены следующие типы ключей: Ideal Switсh – идеальный управляемый ключ; Mosfet – полевой транзистор с изолированными; Gto –запираемый транзистор (Gate turn off); Diode – полупроводниковый диод; Tryristor – упрощенная модель тиристора; Detailed thyristor – уточненная модель тиристора; IGBT – силовой биполярно-полевой модуль типа JGBT; Universal Bridge – универсальный модуль моста.

10.1.5. Блоки измерений Measurements

Эти блоки показаны на рис. 10.7. Блоки Voltage Measurements, Current

Measurements предназначены для соединений измерительных блоков биб-лиотеки Simulink c блоками пакета Pоwer System Blockset. Блок Impedance Measurement позволяет измерить частотную зависимость полного сопро-тивления между двумя точками исследуемой схемы.

Рис. 10.7. Библиотека Measurements

Особый интерес представляет блок Multimeter. Этот блок позволяет

измерить электрические переменные исследуемой схемы, в которой уста-новлены измерительные приборы Voltage Measurement, Current Measurement, Impedance Measurement. Окно настройки этого блока содер-жит два поля. В левом поле после нажатия кнопки Refresh появляются измеряемые пременные. Все или часть из них могут быть переведены во второе (правое) поле с помощью кнопки Select для измерения и регистра-ции результатов. Блок Multimeter своим выходом может быть подключен к внешним измерителям. Флажок Dishlay signals fn simulation stop позволяет

149

вывести измеряемые сигналы в виде переменных зависимостей.

10.1.6. Библиотека электрических машин Machines

Эта библиотека (рис. 10.8) содержит синхронные, асинхронные маши-ны и машины постоянного тока. Все машины могут быть представлены как в абсолютных, так и относительных единицах. Универсальный блок Measurement demultiplexer позволяет измерить требуемые переменные со-стояния машины.

Рис. 10.8. Библиотека Machines

На рис. 10.9 в качестве примера показана панель настройки асинхрон-ной машины.

В полях окна настройки задаются: в первом поле – тип ротора: короткозамкнутый или фазный;

150

во втором поле – координатная система отсчета;

Рис. 10.9. Окно настройки асинхронного электродвигателя

в третьем поле – мощность, действующее линейное напряжение и часто-та; в четвертом, пятом и шестом полях – параметры классической схемы замещения; в седьмом поле – момент инерции ротора, коэффициент трения и число

151

пар полюсов; в последнем поле – начальные условия.

Универсальный блок измерения Machines Measurements Demax позво-ляет измерять выбранные в окне настройки этого блока перменные состоя-ния электрической машины (рис. 10.10).

Рис. 10.10. Окно настройки блока Machines Measurements Demax

10.1.7. Расширенная библиотека Powerlib Extras

Библиотека Powerlib Extras представлена на рис. 10.11. Эта библиотека содержит шесть дополнительных библиотек: Measurement – библиотека дополнительных блоков измерения; Discrete Measurement – дискретные блоки измерений; Control Blocks – блоки управления; Discrete Control Blocks – дискретные блоки управления;

152

Three-Phase Library – библиотеки трехфазных цепей; Additional Mfchines – библиотека дополнительных машин.

Рис. 10.11. Библиотека Powerlib Extras

Многие дополнительные библиотеки содержат блоки, раширяющие

возможности основной библиотеки Power System Blockset. Особый интерес представляет блок Three-Phase Library, содержащий трехфазнае цепи раз-личного назначения (рис. 10.12).

Рис. 10.12. Библиотека Three-Phase Library

153

В этой библиотеке имеется набор последовательных и параллельных трехфазных нагрузок, заданных либо пассивными параметрами, либо зна-чениями активной и реактивной мощности; блок источника питания; трех-фазная индуктивно-связанная цепь; трехфазный трансформатор; блок, мо-делирующий повреждения в трехфазной сети ( блок 3-Phase Fault).

10.2. Пример расчета электрических цепей

постоянного тока Для примера рассмотрим многоконтурную электрическую схему, со-

держащую три контура, шесть сопротивлений и два источника постоянно-го напряжения (рис. 10.13).

Рис. 10.13. Трехконтурная электрическая схема

Для расчета подобных электрических цепей обычно применяют раз-

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

При использовании пакета расширения Power System Blockset расчет трехконтурной схемы значительно упрощается и сводится к составлению электрической схемы с помощью библиотек пакета Power System Blockset и задания параметров элементов схемы.

Для примера пусть заданы следующие значения элементов схемы: R1 = 20 Ом, R2 = 80 Ом, R3 = 100 Ом, Е1 = 100 В, R4 = 35 Ом, R5 = 150 Ом, R6 = 100 Ом, Е2 = 150 В.

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

Осуществляется запуск системы MATLAB 6.0 из главного меню опе-

DC

DC

R4 I4 R5

R6

R1

E1

E2

R3

R2

I5

I6I1

I3

I2

c

e

ba

154

рационной системы либо двойным щелчком левой кнопки мыши на иконке MATLAB на рабочем столе. После запуска появляется основное окно сис-темы MATLAB. Это обычное окно приложений Windows, его можно пере-мещать, изменять его размер, открывать на весь экран и т. д.

Для запуска пакета Simulink необходимо нажать кнопку Simulink в па-нели инструментов главного окна системы MATLAB (см. рис.1.4, поз.8). При нажатии этой кнопки открывается окно браузера библиотек (см. рис.2.1).

Создание модели осуществляется в соответствии с п. 2.2 и использова-нием библиотек Simulink – Sincks, Sources и библиотек Power System Blockset – Elements, Electrical Sources.

Составленная модель электрической схемы постоянного тока (см. рис. 10.13), имеет следующий вид – рис. 10.14:

Рис.10.14. Модель электрической схемы с результатами

расчета токов в ветвях схемы

При составлении схемы следует помнить, что все элементы библио-тек Simulink являются элементами направленного действия, т. е. у каждого элемента есть вход и выход. Поэтому необходимо следить за тем, чтобы выход предыдущего элемента был направлен ко входу последующего (рис. 10.15).

Рис. 10.15. Пример соединения элемнтов

Модель запускается на просчет нажатием на кнопку «►» на панели

155

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

10.3. Пример расчета однофазных электрических цепей

синусоидального тока

Порядок составления модели электрической цепи синусоидального то-ка аналогичен порядку составления модели цепи постоянного тока. Для примера рассмотрим электрическую схему синусоидального тока (рис 10.16).

Рис. 10.16. Однофазная разветвленная электрическая цепь синусоидального тока

Модель однофазной разветвленной цепи синусоидального тока , со-

ставленная с помощью пакетов расширения Power System Blockset – Elements, Electrical Sources, приведена на рис. 10.17.

На примере этой модели можно показать настройки элементов схемы. На рис. 10.18 показана иконка настройки измерителя напряжения. Основ-ными параметрами настройки являются: Cut – удалить; Copy – копировать в буфер; Clear – очистить; Mask parametrs..., Block parametr... и Block properties – при выборе этих команд активизируются окна настроек значений основных параметров и свойств устройства; Real-Time Workshop – позволяет создать подсистему или генерировать S – функцию; Look under mask – позволяет просмотреть внутреннюю модель устройт-ва, т.е. из каких подсистем состоит устройство; Format – открывает возможность работать с изображением элемента схемы: изменять шрифты надписей элемента (Font), перемещать надпись-название элемента относительно самого элемента (Flip name), убирать над-

R1

XL1

R3

XL3

R2

XC2

I

I3 I2

R4

U

156

Рис. 10.17. Модель однофазной разветвленной цепи синусоидального тока

Рис. 10.18. Иконка настройки измерителя напряжения

157

пись-название элемента (Hide name), разворачивать элемент на 1800 (Flip block), поворачивать элемент на 900 (Rotate block), оттенить изображение элемента (Show drop shadow); Foreground color – создание цветных линий, изображающих блок; Background color – закрашивание различными цветами тела элемента.

Результаты моделирования схемы можно увидеть на осциллографе, где показаны осциллограммы входного напряжения и токов в ветвях (рис. 10.19).

Рис. 10.19. Осциллограммы токов и напряжения

Определенный интерес представляет блок Powergui, представленный в

модели однофазной разветвленной цепи (см. рис. 10.17). Этот блок откры-вает все возможности графического интерфейса пользователя.

Открыть блок Powergui можно, щелкнув левой мышью по самому бло-ку и выбрав команду Open block. При этом открывается иконка Powergui (рис.10.20).

Интерфейс Powergui имеет несколько режимов работы: States –вывод данных о состоянии модели; Measurements –вывод данных измерений; Sources – вывод данных об источниках сигналов; Nonlinear – вывод данных о нелинейных параметрах.

Меню интерфейса имеет два пункта: Tools – доступ к средствам инициализации и контроля состояний; Report – генерация отчета о моделировании.

158

Рис. 10.20. Окно установки параметров Powergui

Вывод данных измерения становится возможным, если частота сину-

соидального тока источника напряжения или тока соответствует частоте в окне настройки Powergui. Сверху кнопки частоты находится кнопка выбо-ра максимального или действующего значения данных измерения.

Интерфейс Tools позволяет построить частотную характеристику об-щего сопротивления Zобщ ( режим – Impedance vs Frequency Measurement), а также вызвать блок Control System Toolbox (LTI Viewer), рассмотренный в разделе моделирования систем автоматического регулирования.

10.4. Пример моделирования асинхронного электрического

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

В составе библиотеки Power System Blockset имеются разделы, посвя-щенные моделированию синхронных и асинхронных машин, а также ма-шин постоянного тока.

Асинхронные машины получили широкое распространение в основном как двигатели для электрического привода. Это связано с их простотой и нетребовательностью. Пример моделирования асинхронного двигателя, питаемого от трехфазной сети переменного тока, показан на рис. 10.21. Окна установки параметров асинхронной машины и окно измерения ее па-раметров были приведены на рис.10.9 и 10.10.Последнее окно дает пред-ставление о системе параметров асинхронной машины, доступных для ре-гистрации с помощью виртуальных измерительных приборов.

159

Рис. 10.21. Модель асинхронной машины с короткозамкнутым

ротором

. Рис.10.22. Результаты моделирования асинхронной

машины при линейно нарастающей нагрузке . Приведенные на рис. 10.22 параметры имеет размерность: Speed – c-1, а

Moment – кгм.

0 60 120 180 240 300 Время, мин

160

Библиографический список

1. Гультяев А. Визуальное моделирование в среде MATLAB: Учебный курс.– СПб.: Питер, 2000.

2. Гультяев А. MATLAB 5.2. Имитационное моделирование в среде Windows.– СПб.: КОРОНА-принт, 1999.

3. Говорухин В., Цибулин В. Компьютер в математическом исследовании. Maple, MATLAB, Latex: Учебный курс.– СПб.: Питер, 2001.

4. Дьяконов В., Круглов В. MATLAB. Анализ, идентификация и моделирование систем: Специальный справочник. – СПб.: Питер, 2002. – 448 с.: ил.

5. Дьяконов В., Круглов В. Математические пакеты расширения MATLAB: Спе-циальный справочник. – СПб.: Питер, 2001. – 480 с.: ил.

6. Дьяконов В. Simulink: Специальный справочник. – СПб.: Питер, 2002. – 528 с.: ил.

7. Дьяконов В.П. MATHCAD 8/2000: Специальный справочник. –СПб.: Питер, 2000. – 590 с.

8. Дьяконов В. П. MATLAB 6: Учебный курс.– СПб.: Питер, 2001. 9. Дьяконов В. П. Компьютерная математика: Теория и практика.– М.: Нолидж,

2001. 10. Дьяконов В., Новиков Ю., Рычков В. Компьютер для студента: Самоучи-

тель.– СПб.: Питер, 2000. 11. Дьяконов В. П., Абраменкова И.В. MATLAB 5.0/5.3. Система символьной

математики.– М.: Нолидж, 1999. 12. Дьяконов В. П. Расширяемые системы для численных расчетов MATLAB//

Монитор- Аспект:– 1993:– №2. 13. Дьяконов В. П. Справочник по применению системы PC MATLAB.– М.:

Наука; Физматлит, 1993. 14. Данилов А.И. Компьютерный практикум по курсу “Теория управления”.

Simulink-моделирование в среде MATLAB: Учебное пособие /Под ред. А.Э. Софиева. – М.: МГУИЭ, 2002.– 128 с.

15. Лазарев Ю. MATLAB 5.0. Библиотека студента.– Киев: Ирина, BHV, 2000. 16. Медведев В.С., Потемкин В.Г. Control System Toolbox. MATLAB 5 для сту-

дентов.– М.: Диалог-МИФИ, 1999. 17. Мартынов Н.Н., Иванов А.П. MATLAB 5.х: Вычисления, визуализация, про-

граммирование.– М.: КУДИЦ-ОБРАЗ, 2000. 18. Потемкин В.Г. MATLAB 5 для студентов.– М.: Диалог-МИФИ, 1998. 19. Потемкин В.Г. Система инженерных и научных расчетов MATLAB 5.x.– М.:

Диалог-МИФИ, 2000. 20. Рудаков П.И., Сафонов В.И. Обработка сигналов и изображений. MATLAB

5.х. / Под общ. ред. В.Г. Потемкина.– М.: Диалог-МИФИ, 2000.

161

Учебное издание

Щербаков Виталий Сергеевич, Руппель Алексей Александрович,

Глушец Виталий Алексеевич

ОСНОВЫ МОДЕЛИРОВАНИЯ

СИСТЕМ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ

И ЭЛЕКТРОТЕХНИЧЕСКИХ СИСТЕМ

В СРЕДЕ MATLAB И SIMULINK

Учебное пособие

* * *

Подписано к печати 29.01.2004 г. Формат 60х90 1/16. Бумага ксероксная

Оперативный способ печати Гарнитура Times New Roman

Усл. п. л. 10,0, уч – изд. л. 10,0 Тираж 200 экз. Заказ

Цена договорная

Издательство СибАДИ 644099, Омск, ул. П.Некрасова, 10

Отпечатано в ПЦ издательства СибАДИ 644099, Омск, ул. П.Некрасова, 10