ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 23.01.2025
Просмотров: 2520
Скачиваний: 3
СОДЕРЖАНИЕ
В.Г. Новиков моделирование систем
Тема 3. Моделирование динамических процессов с помощью пакета MatLab “Simulink”29
1.2. Ввод и преобразование моделей
1.3.1.1. Примеры анализа во временной области
1.3.1.2. Примеры анализа в частотной области
1.3.2. Примеры анализа во временной области
1.3.3. Примеры анализа в частотной области
1.3.4. Примеры применения процедур для вычисления полюсов и нулей системы
1.4. Задание на самостоятельную работу
1.4.1. Структурная схема модели
1.4.2. Требования к исследованию модели
1.4.3. Исходные параметры модели
1.4.4. Отчетность по лабораторной работе
2.2. Общие средства фильтрации. Формирование случайных процессов
2.2.1. Общие основы линейной фильтрации
2.2.2. Формирование случайных процессов
2.3. Спектральный и статистический анализ
2.3.2. Примеры спектрального анализа
2.4. Задание на самостоятельную работу
2.4.1 Цифровая обработка сигналов
2.4.2. Формирование случайного процесса из белого шума
2.4.3. Формирование процесса как сумма гармоник
2.4.5. Отчетность по лабораторной работе
Тема 3. Моделирование динамических процессов с помощью пакета MatLab “Simulink”
3.1. Краткие сведения о подсистеме MatLab simulink
3.1.1. Запуск подсистемы simulink
3.1.3. Некоторые основные приемы подготовки и редактирования модели
3.5.1. Установка параметров моделирования
3.5.2. Выполнение моделирования
3.2.2.2. Результат составления модели
3.2.3. Результаты моделирования
3.3. Задание на самостоятельную работу
4.1.1. Особенности процессов в линейных нестационарных системах
4.1.2. Особенности процессов в нелинейных сау
4.1.3. Метод фазового пространства
4.2.1. Приведение линейного дифференциального уравнения к канонической форме
4.2.2. Пример исследования нестационарных и нелинейных процессов
4.2.3. Методика формирования модели в системе MatLab
4.2.2.3. Пример результатов исследования
4.3. Задание на самостоятельную работу
5. 2. Пример имитационного моделирования процессов с использованием подсистемы MatLab simulink
5.2.2. Формирование схемы моделирования
5.2.3. Подготовка к имитационному моделированию
5.2.3.1. Настройка блока Random Number - источника случайного сигнала с нормальным распределением.
5.2.3.3. Настройка блока To Workspace - блока записи в рабочую область Matlab
5.2.3.4. Установка параметров моделирования и сохранение модели
5.2.3.5. Представление результатов имитационного моделирования
5.2.4. Результаты моделирования
Примечание. Для обеспечения удобства чтения схем при размещении блоков в структурной схеме использован поворот отдельных блоков относительно исходного направления. По умолчанию SIMULINK задает направление вход – выход слева направо. (На рисунке изменено направление размещения блока Clock – источник времени и блока Look-Up Table - задание одномерной таблицы).
Данная операция выполняется следующим образом.
выделяется блок (левой кнопкой мыши) щелчком правой кнопки вызывается контекстное меню в подменю Format активизируется опция Rotate Block (поворот блока на 900), для изменения направления размещения (поворота на 1800) можно использовать опцию Flip Block.
Для уменьшения количества одновременно отображаемых на экране блоков, представленную схему можно свестикодному отдельномублоку – кподсистеме. Проделаем данную операцию с полученной схемой для формирования коэффициентаK0.
С этой целью:
л
евой
кнопкой мыши, не отпуская, выделяются
все блоки (см. рисунок 9)
правой
кнопкой вызывается контекстное меню
из этого меню выбирается опция Create
Subsystem
по этой опции все блоки выделяются в
одну подситем
для удобства можно изменить название
подсистемы.
Результат такого объединения представлен на рисунке 10.
Примечание. Состав подсистемы можно восстановить двумя щелчками левой кнопки мыши
Интегро-дифференцирующее звено W1 с переменной постоянной времени T1(t)
Приведенное на схеме (см. рисунок 8) операторное выражение W1представим в виде структурной схемы, которая будет аналогична представленной на рисунке 6.Схема, реализующая указанное звено, примет вид, представленный на рисунке 11.
Д
ля
удобства группу блоков на рисунке 11,
формирующих переменный коэффициентb1
можно объединить в один и присвоить
имяb1.
Результат представлен на рисунке 12.
Усилительное звено с переменным передаточным коэффициентом K1 (t)
Звено “набирается” таким же образом, как и звено K0 (t).
Результирующая схема представлена на рисунке 13.
П
осле
последующего объединения блоков схемы
в один блок получим соответствующую
подсистему, приведенную на рисунке 14.
3).Соединение всех звеньев в соответствии с исходной структурной схемой.
Помимо набранных
выше звеньев с переменными коэффициентами
дополним модель нелинейностью N(типа “насыщение”
),интегратором
,
а также сумматором
.
И
тоговая
схема для моделирования в нашем случае
может выглядеть так, как это представлено
на рисунке 15.
4). Исследование системы.
Для исследования процессов можно использовать любой из способов представления графиков:
- либо осциллографа Scope;
- либо графопостроительXY Graf;
- блок записи в рабочую областьMatLab To Workspase исследуемых переменных на выходе нужного блока споследующимпостроением графиков с использованием оператораplot.
4.2.2.3. Пример результатов исследования
Оценим отработку системой входных сигналов:
- единичная ступенька;
- синусоидальный сигнал с заданными параметрами.
Уровень ограничения заданной нелинейности типа “насыщение” для определенности ограничим уровнем 0.5.
Результаты моделирования в виде переходных процессов на выходе приведены на рисунках 16 и 17.
К
ак
следует из рассмотрения рисунка 17,
амплитуда выходного сигнала уменьшается
с течением времени, что вызвано
переменностью коэффициента передачиK0, K1
во времени.
Получим фазовую траекторию ошибки Δ при задании на входе синусоидального сигнала A sin(ω*t) при A=1 и ω=1/рад/c.
Для
этого, как сказано выше, необходимо
иметь помимо ошибки Δ
скорость dΔ/dt.
Ее легко получить, если воспользоваться
блоком вычисления производной Derivative
из
группы блоков Continuous.
На рисунке 18 приведена схема (аналогичная представленной на рисунке 15), снабженная блоком вычисления производной. А на рисунке 19 – фазовая траектория ошибки dΔ/dt=f(Δ).
4.3. Задание на самостоятельную работу
4.3.1. Исходная система
Дана нестационарная и нелинейная система, структурная схема которой представленная на рисунке 20.
г
деx, y– вход и выход системы;Wi
(s)–
операторные выражения передаточных
функций системы;Δ – ошибка
регулирования;N
– нелинейное звено.
Исходные параметры системы для различных вариантов заданы в таблицах 2 и 3.
4.3.2. Требования к работе
4.3.2.1.“Набрать ” модель с использованием пакетаSIMULINK.
Таблица 2. Характеристики объекта управления
|
t, c |
0 |
2 |
5 |
8 |
10 |
14 |
20 |
|
KОУ, м/с |
10 |
60 |
80 |
50 |
40 |
20 |
14 |
|
fОУ,Гц |
1.0 |
4.2 |
5.0 |
4.0 |
2.7 |
2.3 |
1.8 |
Характеристики динамических звеньев определяются следующими выражениями:
Корректирующее звено:
,
(28)
где
,(29)
а fОУ (t) – из таблицы 2, T2 и T3 - из таблицы 3.
Переменное усилительное звено:
,
(30)
где KОУ(t), fОУ (t) – известные зависимости из таблицы 2.
“Нестационарный” объект управления:
,
(31)
где
,
(32)
а KОУ(t) – из таблицы 2,ОУ – из таблицы 3.
Кинематическое звено
.
(33)
Нелинейное звено N – одним из двух видов
- первый, определяемый зависимостью
,
(34)
где
- нелинейность типа насыщение с порогами
ограничения на уровне±1;
kн – коэффициент передачи до ограничения сигнала,kн = 0.20;
- второй, определяемый зависимостью
,
(35)
где
- нелинейность типа ”насыщение” с
порогами ограничения на уровне±8,
ограничивающий входной сигнал Uвх.
Рекомендация:
При формировании нелинейности второго
типа воспользуйтесь блоком Trigonometric
Function
из раздела
блоков Math
Operations.
4.3.2.2. Провести исследование системы при задании на вход:
скачкообразного и
гармонического сигнала x=Aвх sin (2π fвх t)придвухзначениях частоты:fвх=0.5 Гц и fвх=1 Гц с заданной амплитудойAвх.
При этом провести анализ качества процессов на выходе системы:
- перерегулирование;
- время переходного процесса;
- построить фазовые траектории.
4.3.2.3. Результаты представить в виде графиков.