181121 (628925), страница 2
Текст из файла (страница 2)
Определяем сумму квадратов отклонений:
Рисунок 1 – График статической модели 1-го порядка
1.2 Статистическая модель объекта второго порядка
В целом ход действий аналогичен случаю для линейной модели. Модель объекта второго порядка описывается уравнением вида y=ax2+bx+c.
Для решения этой системы воспользуемся матричным методом наименьших квадратов.
Составим матрицы входных и выходных сигналов:
Таким образом, получили матричное уравнение:
,
где - матрица коэффициентов полинома второго порядка
Находим значение главного определителя:
Δ=314160
Подставляя матрицу поочередно в первый, второй и третий столбец матрицы
, находим вспомогательные определители:
Находим коэффициенты полинома:
Таким образом, получили полином второго порядка:
Для качественной оценки полученного полинома вычислим аналитические значения функции и сравним их с экспериментальными данными. Результаты сведем в таблице 3.
Таблица 3 – Результаты расчета
X | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 9 |
Yзад. | 3 | 4,1 | 5 | 6 | 7 | 7,5 | 7,8 | 8,2 | 9 |
Yаналит. | 3,155 | 4,265 | 5,261 | 6,143 | 6,991 | 7,565 | 8,105 | 8,531 | 9,041 |
ΔY | 0,155 | 0,165 | 0,261 | 0,143 | -0,089 | 0,065 | 0,305 | 0,331 | 0,041 |
ΔY2 | 0,024 | 0,027 | 0,068 | 0,020 | 0,008 | 0,004 | 0,093 | 0,110 | 0,002 |
Далее приведен проверочный расчет линейной аппроксимации на ЭВМ в программной среде MathCAD.
-
векторы данных;
-
длина вектора
- задание степени
- переход к созданию матрицы Вандермонда и подматрицы для решения системы уравнений;
- матрица коэффициентов системы уравнений;
- вектор правых частей системы уравнений;
- решение системы уравнений;
- коэффициент c;
- коэффициент b;
- коэффициент a;
- вычисление значений аппроксимирующей функции;
Определяем сумму квадратов отклонений:
Рисунок 2 – График статической модели 2-го порядка
1.3 Расчёт коэффициентов передачи объекта
Коэффициент передачи объекта показывает, в какую сторону и в какой степени происходит изменение сигнала при прохождении его через объект, то есть усилительные свойства объекта.
Коэффициент передачи определяется как производная от выходной величины:
Расчет коэффициента передачи производим при 10%, 50% и 90% номинального режима, из таблицы данных находим максимальное и минимальное значения сигнала на выходе объекта
Решим алгебраические нелинейные уравнения, исследуем полученные корни и, подставив, нужный корень, получили коэффициенты передач.
Таблица 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.
2.1 Динамическая модель объекта 1-го порядка без запаздывания
Расчёт вручную
Решением дифференциального уравнения будет являться сумма общего и частного решений:
Найдем Yобщ(t):
Найдем Yчаст(t):
Подставим: