Математическое моделирование объектов и систем управления

Москва 2019

МИНИС ТЕРС ТВО НАУКИ И ВЫСШ ЕГО О Б РА З О ВА Н И Я РФ

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ 
ВЫСШЕГО ОБРАЗОВАНИЯ 
«НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ТЕХНОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ «МИСиС»

ИНСТИТУТ ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ  
И АВТОМАТИЗИРОВАННЫХ СИСТЕМ УПРАВЛЕНИЯ 
 
Кафедра автоматизации

Н.В. Осипова

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ 
ОБЪЕКТОВ И СИСТЕМ УПРАВЛЕНИЯ

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

Рекомендовано редакционно-издательским 
советом университета

№ 3312

УДК 004:001.573 
 
О-74

Р е ц е н з е н т 
д-р техн. наук, проф. Л.Д. Певзнер (МИРЭА)

Осипова Н.В.
О-74  
Математическое моделирование объектов и систем управления : учеб. пособие / Н.В. Осипова. – М. : Изд. Дом НИТУ  
«МИСиС», 2019. – 67 с.
ISBN 978-5-906953-66-7

Учебное пособие состоит из трех разделов. В нем рассмотрены основные 
задачи, возникающие при моделировании производственных процессов и систем управления: идентификация объекта управления, настройка регуляторов, 
фильтров Чебышева, Баттерворта, Бесселя, Кауэра и Калмана-Бьюси. 
Предназначено для студентов, обучающихся в магистратуре по направлению подготовки 27.04.04 «Управление в технических системах».

УДК 004:001.573

 Н.В. Осипова, 2019
ISBN 978-5-906953-66-7
 НИТУ «МИСиС», 2019

ОГЛАВЛЕНИЕ

Предисловие ..............................................................................................4
1. Идентификация объектов управления ................................................5
1.1. Идентификация динамических характеристик  
промышленного объекта ......................................................................... 5
1.2. Идентификация статических характеристик  
промышленного объекта ....................................................................... 14
2. Методы настройки регуляторов .........................................................30
2.1. Общие сведения о регуляторах ..................................................... 30
2.2. Метод разгонной характеристики ................................................. 32
2.3. Методы Циглера–Никольса и CHR ............................................... 39
3. Основные типы фильтров ..................................................................46
3.1. Синтез фильтров для систем автоматического управления ....... 46
Контрольные вопросы ........................................................................... 55
3.2. Настройка фильтра Калмана–Бьюси для оценки параметров 
движения объекта ................................................................................... 56
Заключение ..............................................................................................64
Библиографический список ...................................................................65

ПРЕДИСЛОВИЕ

Целью данного пособия является закрепление знаний, полученных 
при изучении курса лекций, на лабораторных работах и практических 
занятиях. Основная задача состоит в приобретении навыков по идентификации модели промышленного объекта, настройке регуляторов, 
фильтров и умении составлять математические модели в среде Matlab с использованием приложений Simulink, Statistics Toolbox и Signal 
Processing Toolbox.
Учебное пособие выступает в качестве основной литературы для 
освоения дисциплины «Математическое моделирование объектов и 
систем управления».
Разработка моделей необходима для проверки правильности работы объекта, описанного уравнениями, изучения поведения САУ с 
рассчитанными параметрами регуляторов и фильтров в условиях различного рода возмущений. Такой предварительный анализ систем 
обладает быстротой и позволяет экономить время и ресурсы, затраченные на проведение эксперимента на физическом объекте. В связи 
с этим издание учебного пособия, в котором изложена методика построения моделей объектов и систем, является актуальной задачей.
Новизна издания заключается в приведении примеров построения 
моделей идентификации систем, расчета новых типов регуляторов, 
фильтров, изучение которых не предусмотрено в программе бакалавриата.
В учебном пособии автор рассмотрел теоретические основы моделирования объектов и систем в приложениях Simulink, Statistics 
Toolbox, Signal Processing Toolbox, имеющих наиболее полный набор 
функций и доступных для самостоятельного изучения.
Более эффективному освоению учебного материала будут способствовать дополнительные источники, приведенные в списке литературы.
Учебное пособие включает теоретические сведения, схемы для 
моделирования, перечень заданий, контрольные вопросы для самопроверки. 

1. ИДЕНТИФИКАЦИЯ ОБЪЕКТОВ 
УПРАВЛЕНИЯ

1.1. Идентификация динамических 
характеристик промышленного объекта

Разработка систем автоматического управления технологическими 
процессами зачастую связана с процедурой идентификации модели 
объекта, которая заключается в определении ее параметров на основе 
данных результатов наблюдения или эксперимента [4, 9].
Выделяют три метода получения математических моделей: аналитический, экспериментальный и экспериментально-аналитический.
Аналитический метод позволяет формально представить физический процесс технологического объекта в виде математических уравнений, отражающих связи между его параметрами и переменными.
Экспериментальные методы делятся на активные и пассивные.
В первом случае на вход объекта подаются планируемые воздействия (ПВ) в виде ступенчатого 1 или 2, импульсного 3, гармонического 4 сигнала. Они также могут быть случайными, которые 
создаются генераторами шума (ГШ). При пассивном эксперименте 
специально никаких сигналов на систему не подается. Она сама вырабатывает случайные воздействия (СВ) (рис. 1.1.1) [7].

Объект 

управления 

1 

2 

3 

4 

ГШ 

СВ 

ПВ 

вход объекта управления 

выход объекта управления 

Рис. 1.1.1. Схема идентификации модели экспериментальным методом

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

Пример. В качестве объекта управления (ОУ) рассматривают технологический процесс перемешивания суспензии (рис. 1.1.2).

Зумпф

Клапан 
1 

Насос 

Суспензия 

Перемешиватель  

6 

3 

2 
 5

4 

Рис. 1.1.2. Схема технологического процесса перемешивания суспензии

Суспензия, расход которой регулируется клапаном, поступает в 
промежуточную емкость – зумпф, откуда насосом откачивается в перемешиватель. В нем расположена вращающаяся крыльчатка 6 с приводом 3. Перечислим контролируемые параметры и датчики, используемые в данной системе:
 – расходомеры 1 и 2 применяются для измерений расхода суспензии на входе зумпфа и выходе насоса;
 – измеритель объема 5 предназначен для контроля заполнения суспензией перемешивателя;
 – измеритель мощности 4 контролирует ее активную составляющую на приводном валу двигателя перемешивателя.
Для построения динамической модели данной системы определим 
передаточные функции ее отдельных элементов.
Исходные данные для расчета:
 – рабочее напряжение электропривода клапана U = 30 В;
 – скорость перемещения клапана υ = 100 %/мин; 
 – «зона нечувствительности» клапана Δ = ± 3 %;
 – номинальная мощность двигателя перемешивателя Pном = 880 кВт;
 – максимальная производительность насоса Qн.max = 8,8 м3/мин;
 – расход при максимальном открытии клапана Qк.max = 11,7 м3/мин;
 – рабочий объем зумпфа Vз = 14,6 м3;
 – рабочий объем перемешивателя Vп = 44 м3;
 – среднеквадратичное отклонение (СКО) шума измерений расхода суспензии в перемешиватель σ1 = 0,1 м3/мин;
 – СКО шума измерений заполнения перемешивателя σ2 = 1 м3;

– СКО шума измерений мощности двигателя перемешивателя 
σ3 = 10 кВт;
Рассмотрим аналитический метод идентификации ОУ.

Модель клапана

Клапан имеет исполнительный механизм, перемещение x которого 
описывается уравнением [5]:

 
υ 
1 ,
d
d
x
k U
t
=
=
 
(1.1.1)

где k1 – коэффициент передачи клапана.

Применяя преобразование Лапласа к выражению (1.1.1) получаем 
передаточную функцию клапана по каналу «напряжение U на входе 
электродвигателя – степень открытия x»:

 
1
1
( )
( )
( )
( )
.
( )
к
k
x p
px p
k U p
W
p
U p
p
=
→
=
=
  
(1.1.2)

Зная скорость хода v и рабочее напряжение электропривода U, 
можно определить k1 из формулы (1.1.1):

 
1

100
3,33
30
v
k
U
=
=
=
 %/(мин ⋅ В). 
(1.1.3)

При составлении модели клапана обычно передаточную функцию 
замыкают отрицательной обратной связью по положению [6]:

 
1
1
1

1

1
1
( )
.
1
0,3
1
1

к

k
W
p
p
k
p
p
k

=
=
=
+
+
+

  
(1.1.4)

Также необходимо учитывать «зону нечувствительности» Δ при 
перемещении рабочего механизма xн:

 
0,

,
.
н
 | |
,

 | |>

x
x
x
x

≤ ∆

=  − ∆
∆


 
(1.1.5)

На выходе данного элемента присутствует еще одно нелинейное 
звено типа «ограничение» степени открытия клапана от 0 до 100 %:

 
 
,
100 %,

100,
100 %.

н
н
о
н

 0
 

 
 

x
x
x
x

≤
≤

= 
>


  
(1.1.6)