Тема: ПИД-регулятор. Методы параметрической оптимизации. Приближен-ные методы настройки. Принцип действия ПИД-регулятора. Примеры использования ПИД-регуляторов. Методы настройки ПИД-регуляторов. Настройка ПИД-регулятора с помощью MATLAB.
Исполнитель - Паздникова Анастасия Алексеевна, гр. 241-об.
Научный руководитель - Русинов Владислав Леонидович, СКБ "Промышленная робототехника и автоматизация".
Актуальность темы: широкое распространение ПИД-регуляторов в различных областях автоматического управления.
Цель работы: изучение принципов работы ПИД-регулятора, методов его настройки и применение программного обеспечения MATLAB для оптимизации параметров.
Задачи:1) Изучить принцип действия ПИД-регулятора и его составляющих (пропорциональной, интегральной и дифференциальной).
2) Описать различные методы параметрической оптимизации ПИД-регуляторов.
3) Рассмотреть приближенные методы настройки ПИД-регуляторов (методы Циглера-Никольса, Коха и др.).
4) Проанализировать примеры использования ПИД-регуляторов в различных системах.
5) Изучить методы настройки ПИД-регуляторов, включая ручные и автоматические методы.
6) Реализовать настройку ПИД-регулятора для конкретной системы с использованием MATLAB и оценить результаты.
ВВЕДЕНИЕ
В этой исследовательской работе рассматриваются принципы работы и применения пропорционально-интегрально-дифференциального (ПИД) регулятора.
ПИД-регуляторы являются наиболее распространенным типом регуляторов в системах автоматического управления, благодаря своей относительной простоте реализации и эффективности в широком спектре задач.
Работа будет посвящена анализу трех основных компонентов ПИД-регулятора - пропорционального, интегрального и дифференциального - и их вклада в общую производительность системы.
Будет исследовано влияние параметров настройки ПИД-регулятора (Kp, Ki, Kd) на характеристики системы, такие как точность, скорость реагирования, перерегулирование и устойчивость.
Кроме того, будут рассмотрены различные методы настройки ПИД-регулятора регулятора с использованием среды MATLAB/Simulink, сравнены их преимущества и недостатки, и продемонстрирована их эффективность на примере конкретной системы.
Результаты исследования могут быть использованы для оптимизации работы различных систем управления и повышения их эффективности.
Содержание
1 Теоретическая часть.
1.1 Принцип действия ПИД-регулятора.
1.1.1 Математическая модель ПИД-регулятора.
1.1.2 Объяснение влияния каждого компонента (P, I, D) на поведение системы.
1.1.3 Анализ устойчивости замкнутой системы с ПИД-регулятором.
1.2 Методы параметрической оптимизации.
1.2.1Ручная настройка или метод настройки "по графикам".
1.2.2Циглера-Никольса.
1.2.3Метод Коха.
1.2.4 Сравнительный анализ точности и скорости настройки различных приближенных методов.
1.3 Примеры использования ПИД-регуляторов.
2 Практическая часть
2.1 Выбор моделируемой системы.
2.2 Разработка модели в MATLAB/Simulink.
2.3 Настройка ПИД-регулятора. Описание использования различных методов настройки.
2.3.1 Ручная настройка,
2.3.2 Приближенные методы
2.3.3 Методы оптимизации
2.4 Сравнение результатов настройки различными методами и их анализ.
3 Заключение
4 Список литературы
Этот план является ориентировочным и может быть изменен в зависимости от глубины исследования и доступных ресурсов!!!
1 ТЕОРЕТИЧЕСКАЯ ЧАСТЬ1.1 Принцип действия ПИД-регулятора.ПИД регулятор - это механизм, используемый в управляющих контурах с обратной связью для поддержания в заданном интервале одного или нескольких параметров. [1]
Он регулирует некоторую изменяющуюся во времени величину процесса путем вычисления управляющего сигнала, который представляет собой сумму трех составляющих: пропорциональной, интегральной и дифференциальной. Отсюда и название. [1]
Регулируемую величину, которая отклоняется в результате действия возмущения, ПИД-регулятор возвращает в пределы допустимых значений.[2]
Основной принцип работы ПИД-регулятора заключается в следующем. Он сравнивает текущую регулируемую величину с заданным значением (уставкой). На основе их разницы алгоритм вычисляет управляющий сигнал и отправляет его на исполнительный механизм. Тот, в свою очередь, приводит объект управления к желаемому значению процесса (заданному значению).
Они учитывают фактическую величину, заданное значение, разность значений и скорость изменения контролируемых характеристик.
Такие устройства универсальны, при помощи ПИД-регуляторов можно реализовать любые законы регулирования.[2]
1.1.1 Математическая модель ПИД-регулятора.Математическая модель ПИД-регулятора описывает его выходной сигнал как функцию пропорциональную отклонению регулируемого параметра, интегралу от этого отклонения и скорости этого отклонения.
Ошибка (e(t)) определяется как разность между заданным значением (SP) текущим значением (PV):
e(t) = SP - PV(t)
Выходной сигнал ПИД-регулятора (u(t)) вычисляется по формуле:
где Kp - коэффициент пропорциональности, определяющий пропорциональную составляющую выходного сигнала относительно текущей ошибки;
Ki - коэффициент интегральной составляющей. Определяет интегральную составляющую, которая накапливает ошибку во времени;
Kd - коэффициент дифференциальной составляющей. Определяет дифференциальную составляющую;
∫
₀ᵗ e(τ)dτ - интеграл от ошибки по времени от 0 до t;
de(t)/dt - производная ошибки по времени.
Альтернативная форма записи (с использованием оператора Лапласа):
В частотной области (используя преобразование Лапласа), передаточная функция ПИД-регулятора имеет вид:
G(s) = Kp + Ki/s + Kd s ,
где s - комплексная переменная Лапласа.
1.1.2 Объяснение влияния каждого компонента (P, I, D) на поведение системы.ПИД-регулятор вычисляет управляющий сигнал, формируемый ПИД регулятором, сумма трех составляющих:
- Пропорциональной.Т. е. разности входного сигнала и сигнала обратной связи (сигнал рассогласования);
- Интегральной.Т.е интегралу сигнала рассогласования;
- Дифференциальной.Т.е производной сигнала рассогласования.
Схема системы управления с обратной связью с участием ПИД-регулятора представлена на рисунке 1.
Рисунок 1 - Схема системы управления с обратной связью с участием ПИД-регулятора
Kp - коэффициент пропорциональности. Определяет пропорциональную составляющую выходного сигнала относительно текущей ошибки. Большее значение Kp приводит к более быстрому реагированию, но может вызвать колебания;
2) Ki - коэффициент интегральной составляющей. Определяет интегральную составляющую, которая накапливает ошибку во времени. Используется для устранения статической ошибки (остаточной разницы между SP и PV). Большие значения Ki могут привести к перерегулированию;
3) Kd - коэффициент дифференциальной составляющей. Определяет дифференциальную составляющую, которая реагирует на скорость изменения ошибки. Снижает колебания и ускоряет установление. Большие значения Kd могут привести к повышенной чувствительности к шуму.В математической модели, в отличие от реального устройства, шумов не бывает. При настройке этого коэффициента инженер должен найти компромисс между качеством переходного процесса и чувствительностью к шумам.
1.1.3 Анализ устойчивости замкнутой системы с ПИД-регулятором.Анализ устойчивости замкнутой системы с ПИД-регулятором -- важная задача автоматического управления. Устойчивость означает, что система вернется в заданное состояние после возмущения, а не будет расходиться.
Выбор метода зависит от сложности системы и требований к точности анализа. Для простых систем может быть достаточно анализа корней характеристического уравнения или критерия Рауса-Гурвица. Для более сложных систем необходимы частотные методы или численное моделирование. В любом случае, важно помнить, что анализ устойчивости - это критически важная часть проектирования и настройки систем автоматического управления.
Существует несколько методов анализа устойчивости, наиболее распространенные из которых:
1) Анализ корней характеристического уравнения:Замкнутая система с ПИД-регулятором описывается передаточной функцией. Характеристическое уравнение получается из знаменателя передаточной функции замкнутой системы после приведения её к общему знаменателю.
Характеристическое уравнение замкнутой системы находится в знаменателе передаточной функции и записывается в виде:
1 + Wo6(s)WP(s) = 0,
где Wo6(s) - перед передаточная функция объекта системы;
WP(s) - перед передаточная функция регулятора.
Корни этого уравнения равны полюсам si передаточной функции замкнутой системы, которые определяют поведение системы. Устойчивость системы гарантирована, если все корни имеют отрицательные действительные части. Нахождение корней может быть выполнено аналитически (для простых систем) или численно (для сложных систем) с помощью программного обеспечения (например, MATLAB, Mathematica).
2) Критерии устойчивости Рауса-Гурвица - метод, который позволяет
определить устойчивость системы, не вычисляя корни характеристического уравнения, а анализируя коэффициенты полинома характеристического уравнения. Он предоставляет необходимые и достаточные условия устойчивости.
3) Частотный метод анализа (критерий Найквиста):А) Амплитудно-фазовая частотная характеристика (АФЧХ): Этот метод основан на анализе АФЧХ разомкнутой системы. Критерий Найквиста устанавливает условия устойчивости на основе анализа поведения АФЧХ в комплексной плоскости. Он позволяет определить наличие неустойчивых полюсов в замкнутой системе, не вычисляя их непосредственно. График АФЧХ строится путём подстановки `s = jω` в передаточную функцию разомкнутой системы и изменения частоты ω.
Б) Геометрическая интерпретация: Критерий Найквиста использует геометрическую интерпретацию в комплексной плоскости. Условие устойчивости формулируется в терминах числа обходов АФЧХ вокруг точки (-1, 0).
4) Метод корневого годографа:Построение корневого годографа - это графический метод демонстрирующий, как изменяются полюса замкнутой системы при изменении коэффициента усиления разомкнутой системы. Построение корневого годографа позволяет визуально оценить устойчивость системы и определить оптимальное значение коэффициента усиления.
5) Методы моделирования:Численное моделирование: Используется для анализа устойчивости сложных систем, для которых аналитическое решение затруднено. Системы моделируются с помощью специализированного программного обеспечения (MATLAB/Simulink, Python с библиотеками SciPy/NumPy). Моделирование позволяет исследовать реакцию системы на различные входные сигналы и возмущения.
1.2 Методы параметрической оптимизации. Методы параметрической оптимизации - это алгоритмы, которые находят оптимальные значения целевой функции, зависящей от одного или нескольких параметров. Они отличаются от методов непараметрической оптимизации тем, что предполагают знание некоторой структуры целевой функции (например, аналитическое выражение или параметрическая модель).
Выбор конкретного метода зависит от свойств целевой функции (дифференцируемость, выпуклость, наличие шума и т.д.) и вычислительных ресурсов.
1.2.1 Ручная настройка или метод настройки "по графикам".Часто можно услышать о простейшем подходе к настройке, когда вручную настраивают сперва Kp, предварительно обнулив остальные каналы, потом Ki, а последним Kd, оценивая результат работы чисто по графикам процессе (так называемый эмпирический и эвристический методы).
Такая инструкция встречается практически для всех ПЧ и даже сервоприводов, возможно это оправдано тем, что процесс быстрый и графики можно получить без особых затрат просто перебирая коэффициенты по логарифмической шкале, но для температуры так будет настраивать неудобно.
Важное замечание: даже для одного и того же объекта процедуру настройки ПИД регулятора следует повторять, если изменились параметры объекта (например, загрузка печи) или значительно изменилась рабочая температура (остывание и соответственно отток мощности зависит от градиента температур).
1.2.2 Тангенциальный методВ тангенциальном методе для нахождения ПИД коэффициентов используется переходная характеристика объекта, в случае терморегулирования кривая разогрева объекта.
Алгоритм тангенциального метода:
1. Задать значение мощности на нагрузке в режиме ручного управления контроллера. Мощность должна быть такой, чтобы установившаяся температура не превысила допустимого значения.
2. Дождаться установления температуры в этом режиме и в некоторый момент времени скачком изменить значение мощности (например, на 5%).
3. Записать процесс изменения температуры (на компьютере, на самописце или по точкам, вручную) до установления нового устоявшегося значения температуры.
4. Провести касательную к полученной кривой в точке, где наклон касательной будет максимален.
5. Определить «мёртвое» время процесса как промежуток времени, прошедший с момента изменения мощности до момента времени, определяемого точкой пересечения касательной и средней линии температуры первого устойчивого состояния.
6. Определить величину относительного наклона касательной по формуле:
К = (∆T/∆t)/∆N.∆T - изменение температуры в °С,
∆t - соответствующее изменение времени в секундах,
∆N - изменение мощности в %.
7. Установить значения коэффициентов Kp, Ki, Kd в соответствии с таблицей.
Таблица 1 - Параметры типовых регуляторов
Тип регулятора Kp Ki Kd
П-регулятор 100*τм*К - -
ПИ-регулятор 111*τм*К 0,54 Kp*/ T* 3.3*τм
ПИД-регулятор 83*τм*К 2 Kp*/ T* 0.25∗τм*
1.2.3 Метод колебаний (или второй частотный метод Зиглера-Никольса)Он основан на наблюдении реакции системы на ступеньку. Метод не гарантирует оптимальной настройки, но обеспечивает рабочую точку, с которой можно начинать дальнейшую настройку.
Процесс настройки:
1. Вначале устанавливаются I=0 и D=0. Регулятор работает только в пропорциональном режиме.
2. На вход системы подается ступенчатый сигнал (скачок) заданного значения.
3. Наблюдается реакция системы на ступенчатый сигнал. Важно получить устойчивую, затухающую колебательную реакцию. Если система не колеблется, необходимо увеличить коэффициент усиления P. Если колебания слишком затухают, то также следует увеличить коэффициент усиления P. Цель - получить кривую с устойчивыми колебаниями, не раскачивающимися до бесконечности.
4. Из графика реакции системы определяются два параметра:
Ku - Предельное усиление: это значение коэффициента усиления P, при котором система начинает демонстрировать устойчивые колебания постоянной амплитуды.
Pu- Предельный период: это период колебаний, наблюдаемых при Ku. Он измеряется в секундах или других единицах времени, соответствующих системе. Это время между двумя соседними пиками или впадинами колебаний.
5. Используя найденные Ku и Pu, вычисляются параметры ПИД-регулятора по следующим формулам из таблицы 2
Таблица 2 - Параметры типовых регуляторов
Тип регулятора Kp Ki Kd Ti Td
П-регулятор 0,50 Kp* - - - -
ПИ-регулятор 0,45 Kp* 0,54 Kp*/ T* - 0,83 T* -
ПИД-регулятор 0,60 Kp* 1,2 Kp*/ T* 0,075 Kp*/ T* 0,5 T* 0,125 T*
