Пример 2. Моделирование временного отклика динамических систем

Постановка задачи

Определить временные характеристики звена (переходную и весовую характеристики) на примере колебательного звена второго порядка

. (3.1)

Решение.

Для построения схемы моделирования воспользуемся методом понижения производной (методом Кельвина). В нем можно выделить следующие этапы.

• Разрешаем уравнение относительно старшей производной

. (3.2)

• Полагаем старшую производную известной и выполняем ее последовательное интегрирование, получая все низшие производные и саму переменную y. Для этого потребуется два последовательно включенных интегратора, на выходах которых получим сигналы и y.
• Формируем старшую производную, используя уравнение (3.2). Для этого потребуется сумматор, складывающий сигналы и y,умноженные, соответственно, на коэффициенты -0,167 и -0,333.
• Для получения переходной и весовой характеристик на один из входов сумматора подаем единичный сигнал .
• Запускаем систему Matlab, а затем систему Simulink (см. п. 3.3).

Создаем окно новой модели File→New→Model.

• С помощью имеющихся блоков системы Simulink собираем структурную схему (рис.3.4).

Рис.3.4. Схема моделирования для снятия временных характеристик

Для сбора модели используются следующие блоки библиотеки Simulink:

Ø ступенчатая функция – библиотека Sources → Step;

Ø линейные преобразователи – библиотека Commonly Used Bloks → Gain;

Ø интеграторы – библиотека Continuous → Integrator;

Ø сумматор – библиотека → Commonly Used Bloks → Sum;

Ø виртуальные осциллографы – библиотека Sinks → Skope;

Ø производная - библиотека Continuous → Derivative.

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

Рис.3.5.

Для расчета модели необходимо выбрать пункт меню Simulation/Start или инструмент на панели инструментов.

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

При необходимости можно увеличить время расчета Stop time, с 10 с. По умолчанию, до, например, 100 с. Этот параметр задается в окне Simulation/Configuration Parameters. Время расчета задается указанием начального (Start time) и конечного (Stop time) значений времени расчета. Начальное время, как правило, задается равным нулю. Величина конечного времени задается пользователем исходя из условий решаемой задачи.

После повтора расчета результат выглядит, как показано на рис.3.6.

Рис.3.6. Переходная и весовая характеристики соответственно

Пример 3. Моделирование частотного отклика динамических систем

Постановка задачи

Построить ЛАЧХ (линейную амплитудно-частотную характеристику) и ЛФЧХ (линейную фазово-частотную характеристику) звена, заданного с помощью передаточной функции

.

Решение (Вариант 1)

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

Рис.3.7. Первый вариант схемы моделирования для снятия частотных характеристик

Для сбора модели используются следующие блоки библиотеки Simulink:

Ø синусоидальный сигнал – библиотека Sources → Sine Wave;

Ø динамическое звено – библиотека Continuous → Transfer Fcn;

Ø мультиплексор – библиотека Commonly Used Bloks→Mux;

Ø виртуальные осциллографы – библиотека Sinks → Skope.

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

задаются в блоке Sine Wave. Для расчетов амплитуду входного сигнала задаем равной единице.

Блок Mux предназначен для одновременного наблюдения на осциллографе входного и выходного сигналов.

Результат моделирования при представлен на рис.3.8.

Рис. 3.8. Входной сигнал (розовая линия) и отклик системы (желтая линия)

На экране виртуального осциллографа в общем случае будет наблюдаться изображение (рис.3.9)

Рис.3.9

Решение (Вариант 2)

Для получения динамических характеристик системы (рис.3.8) можно воспользоваться функциями Linear Analysis пакета расширения Control System Toolbox. Для этого необходимо выставить метку Input Point щелкнув правой кнопкой мыши по стрелке перед блоком Transfer Fcn и набрать Linearization Points / Input Point. Аналогично выставляется метка Output Point (рис.3.10).

Рис.3.10. Второй вариант схемы моделирования для снятия частотных характеристик

Далее, набрав команду Tools / Control Design / Linear Analysis, открывается окно Control and Estimation Tools Manager

Рис. 3.11

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

Далее нажимаем на кнопку Linearize Model и получаем окно с графическим отображением динамических характеристик (рис.3.12).

 

Рис.3.12 Переходной процесс

 

Для переключения между графиками нужно щелкнуть правой кнопкой мыши по графику и выбрать из подменю Plot Types/«необходимый график». Здесь « необходимый график » - одина из диграмм:

• Step – реакция на единичный скачок;
• Impulse – импульсная характеристика;
• Bode, Bode Magnitude – диаграммы Боде;
• Nyquist – диаграмма Найквиста;
• Nichols - диаграмма Никольса (логарифмическая амплитудно-фазовая диаграмма);
• Pole/Zero – картанулей и полюсов.

На рис.3.13 - 3.14 в качестве примера изображены диаграммы Найквиста и Боде.

Диаграмма Найквиста строится в системе координат: по оси абсцисс откладывается действительная часть передаточной функции при , по ординат – мнимая. Для получения годографа осуществляется изменение частоты от 0 до ∞.

 

Рис.3.13. Фазовый портрет – диаграмма Найквиста

 

Рис.3.14. Частотная характеристика – диаграмма Боде

 

⇐ Предыдущая1234Следующая ⇒

Что делает отдел по эксплуатации и сопровождению ИС? Отвечает за сохранность данных (расписания копирования, копирование и пр.)...

Что будет с Землей, если ось ее сместится на 6666 км? Что будет с Землей? - задался я вопросом...

Живите по правилу: МАЛО ЛИ ЧТО НА СВЕТЕ СУЩЕСТВУЕТ? Я неслучайно подчеркиваю, что место в голове ограничено, а информации вокруг много, и что ваше право...

Что делать, если нет взаимности? А теперь спустимся с небес на землю. Приземлились? Продолжаем разговор...

---

Не нашли то, что искали? Воспользуйтесь поиском гугл на сайте: