
- •«Моделирование объектов и систем управления отрасли»
- •Введение
- •Методы моделирования объектов и систем управления
- •Физическое моделирование объектов и систем управления
- •Метод подобия
- •Метод аналогий
- •Математическое моделирование химико-технологических процессов
- •Сущность, свойства и этапы математического моделирования химико-технологического процесса
- •Основные виды математических моделей химико-технологических процессов
- •Блочный принцип разработки математических моделей химико-технологических процессов
- •Основные подходы получения математических моделей химико-технологических процессов
- •Классификация экспериментальных методов. Виды математических моделей.
- •Модели динамики (статистические) химико-технологического процесса
- •Определение динамических характеристик химико-технологических процессов при активном эксперименте
- •Проведение эксперимента по снятию переходных функций
- •Обработка результатов эксперимента по снятию переходных функций
- •Получение динамических характеристик объекта по переходным функциям
- •Определение динамических характеристик объектов при входных воздействиях в виде случайных сигналов
- •Свойства и характеристики стационарных случайных процессов
- •Постановка эксперимента по определению статистических характеристик процессов
- •Определение динамических характеристик объектов и систем управления методом моментов
- •Применение метода модулирующих функций для определения динамических характеристик объектов и су
- •Получение модели статики объектов и су
- •Определение моделей статики (уравнения регрессии) при пассивном эксперименте
- •Регрессионный анализ в матричной форме
- •Получение модели статики при активном эксперименте
- •Детерминированные (аналитические) модели химико-технологических процессов
- •Использование уравнений материального и теплового баланса для описания стационарных и нестационарных режимов протекания химико-технологических процессов
- •Модели структуры потоков, как основа построения математических моделей гидромеханики химико-технологического процесса
- •Модель идеального перемешивания
- •Модель идеального вытеснения
- •Диффузионные модели
- •Однопараметрическая диффузионная модель (одм)
- •Двухпараметрическая диффузионная модель
- •Ячеечная модель
- •Математическое описание тепловых процессов
- •Прямоточные кожухотрубчатые теплообменники
- •Противоточный кожухотрубчатый теплообменник
- •Математические модели теплообменников с учетом накопления тепла в стенке
- •Прямоточный кожухотрубчатый теплообменник
- •Противоточный кожухотрубчатый теплообменник
- •Математическая модель конденсатора (паровой теплообменник)
- •Передаточные функции теплообменных аппаратов
- •4.5.1 Конденсатор без учета накопления тепла в станке
- •4.5.2 Конденсатор с учётом накопления тепла в стенке
- •4.5.3 Кожухотрубчатый противоточный теплообменник
- •Математическое моделирование массообменных аппаратов
- •Математическая модель процесса газоабсорбции
- •Математическая модель процесса ректификации
- •Передаточная функция ректификационной колонны
- •Математическое моделирование математических процессов
- •Математическая модель процесса химического превращения (кинетическая модель)
- •Математическая модель химического реактора идеального перемешивания
- •Математическая модель реактора идеального вытеснения
- •Передаточные функции химических реакторов
- •Передаточная функция реактора идеального перемешивания простого типа
- •Реактор идеального вытеснения
- •Математическое моделирование динамики типовых контуров регулирования технологических параметров объекта
- •Свойства, характеристики и классификация объектов регулирования
- •Регулирование уровня в резервуарах
- •Регулирование расхода жидких сред
- •Регулирование расхода сыпучего материала
- •Регулирование концентрации вещества в смесителях. Регулирование концентрации вещества в смесителях
- •Регулирование температуры в смесителях идеального перемешивания
- •Регулирование состава смеси в реакторе с мешалкой
- •Методы и алгоритмы идентификации динамических систем
- •Общая классификация задач идентификации
- •Постановка задачи идентификации
- •Идентификация непрерывных объектов в классе моделей с дискретным временем
- •Задача непараметрической идентификации
- •Задачи параметрической идентификации
- •Безпоисковые алгоритмы идентификации с адаптивной моделью (баиам)
- •8.6.1 Общая структура баиам в пространстве сигналов
- •8.6.2. Общая структура баиам с дискретным временем и операторным описанием
- •Оптимизация химико-технологических с использованием математических моделей
- •Постановка задач оптимизации
- •Поисковые модели идентификации с адаптивной моделью (пмиам)
- •Аналитические методы оптимизации
- •Метод классического анализа функций
- •Метод неопределенных множителей Лагранжа
- •Линейное программирование
- •Нелинейное программирование
- •Численные методы оптимизации. Постановка задачи оптимизации
- •Реализация методов оптимизации в пакете MatLab (toolbox)
- •Оптимизация (минимизация) функции
- •9.5.5. Использование метода наименьших квадратов для решения задач оптимизации
Проведение эксперимента по снятию переходных функций
На
объекте устанавливается рабочий режим
и производится измерение всех параметров
технологического объекта. В момент
времени, когда
наносится выбранный вид и амплитуда
испытательного воздействия
и производится непрерывное измерение
,
которая через некоторый момент времени
будет иметь установившуюся величину.
Рис. 7 Реакция объекта на испытательное воздействие
Линейность
статической характеристики не означает
линейность динамической характеристики
в малом. Для подтверждения линейности
динамической характеристики в малом
необходимо убедиться в соблюдении
принципа суперпозиции для чего опыт
проводят при различных знаках и величинах
.
Далее вычисляют коэффициент усиления
.
Если
отличается от 20 до 40%, то необходимо
уменьшить
и провести повторное снятие функции. В
некоторых случаях гладкие переходные
функции при одном и том же режиме и при
одних и тех же фиксированных моментах
времени в одних и тех же опытах имеют
различные значения. Для проверки
стационарности динамических характеристик
необходимо эксперимент по снятию функции
повторить несколько раз через определенные
промежутки времени. Опыты необходимо
проводить при одних и тех же
и
.
Обработка результатов эксперимента по снятию переходных функций
Обработка гладких переходных функций.
Полученные
переходные функции
строятся в одном масштабе, но так как
величина
в опытах изменялась, то переходные
функции необходимо представить в виде
единичных переходных функций (ЕПФ):
Которые строятся в одном масштабе. Если разброс между ЕПФ 1÷2%, то для последующей обработки выбирается одна из ЕПФ. Если разброс больше 2÷3%, то выбирают усредненную ЕПФ:
По
данным ЕПФ определяется время запаздывания
(𝜏)
— отрезок времени для которого выполняется
следующее неравенство
,
где ∆ — значение ЕПФ равное ∆=(0,01÷0,02)
.
Рис. 8 График ЕПФ
∆ — определяется классом точности выбранных приборов.
Сглаживание экспериментальных переходных функций.
Во многих случаях функции получаются гладкими. То есть получается зависимость вида:
где
— помеха (нормальное распределение
величины с
).
генерируется либо в самом объекте, либо наводится в измерительной цепи.
Рис. 8 График случайного процесса
Цель
обработки 2.5.12 получение истинных
значений
в гладкой форме по которой можно выделить
динамические свойства объекта. Для
этого необходимо снять
экспериментальных переходных функций
с тем, чтобы
С
данной зависимости 2.5.13 необходимо снять
от 80÷100
.
В промышленности данный метод не
используется, а используется “сглаживание”
переходной функции
основанное на усреднении заданных
таблично отстоящих друг от друга
на
(таких ординат берется
).
При этом предполагается, что переходная
функция носит неколебательный характер,
то есть корни ХАУ вещественные и меньше
ноля.
Сглаживание переходных функций скользящим усреднением.
Сущность
метода заключается в последовательном
усреднении ординат переходной функции
на некотором интервале
,
где
— цело, чётное число.
Алгоритм согласования имеет следующий вид:
где
— оценка ординаты отнесенная к середине
интервала сглаживания;
— предыдущее значение величины;
— память
у линейного фильтра, АФХ которого имеет
вид:
то
есть данный фильтр не пропускает частоты
выше
откуда находят
.
Увеличение величины может привести к искажению переходной функции и потере части уже сглаженных координат (особенно в начале) при:
Так
как спектр частот
обычно не известен, а также не известна
спектральная плотность помехи
,
то величину
находят экспериментально. При этом,
чтобы получить значение всех координат
измерение
начинается до момента нанесения активного
воздействия
,
а также после завершения процесса.
Начальный участок определяет начальную
структуру, который будет аппроксимировать
данную переходную функцию, а последний
участок определяет коэффициент усиления
модели.
Сглаживание методом четвертых разностей.
Сущность
метода заключается в аппроксимации с
помощью метода наименьших квадратов
каждых пяти соседних ординат
параболой второго порядка, затем
вычисляется разность между средним
значением
из данных пяти ординат и параболой
второго порядка. Величина данной поправки
равна центральной четвертой разности
функции
:
Далее находим значение функции:
В данном случае мы теряем четыре ординаты:
Для того чтобы восполнить эти потери применяются следующие формулы:
Кроме рассмотренных методов сглаживания применяется сглаживание рядами Фурье и полиномами Чебышева, но они очень сложные.
Усреднение переходных функций методами четвертых разностей дает лучшие результаты, потому что этот метод основан на определении, что на интервале усреднения переходная функция имеет вид параболы.