Решим алгебраические нелинейные уравнения, исследуем полученные корни и, подставив, нужный корень, получили коэффициенты передач.
Таблица 4 – Результаты расчета коэффициентов передачи
10% | 50% | 90% | |
y | 3,6 | 6 | 8,4 |
k | 1,209 | 1,417 | 1,598 |
2 ПОСТРОЕНИЕ ДИНАМИЧЕСКОЙ МОДЕЛИ ОБЪЕКТА
Динамическая характеристика объекта нужна для построения его динамической математической модели, которая описывает поведение объекта во времени, начиная с момента подачи входного сигнала и до момента, когда все переходные процессы заканчиваются.
Динамические характеристики, в свою очередь, подразделяются на временные и частотные.
Временными характеристиками звена или системы называют изменение во времени значений выходной величины при поступлении на вход некоторого типового воздействия. Наиболее важной временной характеристикой является реакция системы на единичное, мгновенное, скачкообразное изменение значения входной величины, так как этот режим очень часто возникает в системах регулирования, как при включении, так и при изменении заданного значения регулируемой величины.
Таким образом, под временной характеристикой системы будем понимать процесс изменения выходной величины в функции времени при переходе системы из одного равновесного состояния в другое в результате поступления на вход системы единичного, ступенчатого воздействия.
Для получения динамической переходной характеристики объекта регулирования необходимо:
а) задаться рядом значений времени t;
б) зарегистрировать значение выходного сигнала Yi в заданные моменты t, в результате интенсивного экспериментирования. Эти данные сведены в таблицу 5.
Таблица 5 – Динамическая характеристика объекта регулирования
i | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
t | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 9 |
Y | 0 | 0,1 | 0,5 | 0,7 | 0,82 | 0,91 | 0,975 | 0,99 | 1 |
Для получения аналитической зависимости, заданную таблично динамическую характеристику необходимо аппроксимировать выражением первого порядка. Затем по наименьшему значению суммы квадратов отклонений для характеристик без запаздывания и с запаздыванием нужно выбрать наиболее приближенную к экспериментальным данным динамическую характеристику.
Для статических объектов первого порядка без запаздывания будем иметь:
- дифференциальное уравнение
Т - постоянная времени, это время в течение которого выходная координата объекта достигла установившегося состояния, если бы изменялась с максимальной скоростью;
k - коэффициент передачи.
- передаточную функцию
- решение дифференциального уравнения
Для статических объектов первого порядка с запаздыванием будем иметь:
- дифференциальное уравнение
- передаточную функцию
- решение дифференциального уравнения
Для статических объектов N-го порядка без запаздывания, имеющих кратные (одинаковые корни), будем иметь:
- передаточную функцию
- решение дифференциального уравнения (переходный процесс)
Для статических объектов N-го порядка с запаздыванием будем иметь:
- передаточную функцию
переходный процесс
Общий процесс решения поставленной задачи будет выглядеть следующим образом:
- составляем систему уравнений для всей совокупности экспериментальных точек. Для объектов высоких порядков систему предварительно линеаризуем посредством замены переменных и решения в каждой точке нелинейного уравнения для этой введенной дополнительной переменной;
- составляем систему уравнений для всей совокупности экспериментальных точек. Для объектов высоких порядков систему предварительно линеаризуем посредством замены переменных и решения в каждой точке нелинейного уравнения для этой введенной дополнительной переменной;
- полученную систему уравнений решаем матричным методом наименьших квадратов и находим неизвестные и T.
Расчёт вручную
Решением дифференциального уравнения будет являться сумма общего и частного решений:
Найдем Yобщ(t):
Найдем Yчаст(t):
Подставим:
Найдем постоянную С:
Y(0)=Y0
тогда
;Подставим С:
По таблице при t = 0, Y(0) = 0, тогда:
, откуда и получим:Где
- установившееся значение, в нашем случае Yуст = Ymax.Найдем постоянную времени Т методом наименьших квадратов. Преобразуем выражение:
Прологарифмируем выражение:
Обозначим
Рассчитаем
для каждого момента времени ti и занесем в таблицу 6.Таблица 6 - Значения
i | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
1 | 0.9 | 0.5 | 0.3 | 0.18 | 0.09 | 0.025 | 0.01 | 0 |
Далее составляем систему алгебраических уравнений:
В систему не вошло уравнение для момента времени t9, так как ln0 не определен, и для момента времени t=0, так как в них
. Составляем матричное уравнение для:Составим матрицы:
Находим произведение
: