Система в дискретном времени имеет вид:

dt=24 c.

Зададим
,
, получим значения на выходах дискретной системы, которые совпадают с расчетом задания в п.4.
Таблица 6 Значение выходов дискретной системы
| Возмущение | Реакция выхода системы y(t) | |||||||||||
| u1=0.01
u2=0 | y1 y2 10-3 | 0 | 0 | 3.874 | 6.247 | 7.701 | 8.591 | 9.137 | 9.471 | 9.676 | 9.802 | 9.878 |
| 0 | 0 | -2.548 | -3.523 | -3.896 | -4.038 | -4.093 | -4.114 | -4.122 | -4.125 | -4.126 | ||
| такт | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | |

Рисунок 10 – Реакция выходов системы на возмущения u (t)
Построение графиков кривой разгона нелинейной системы
Данные для построения графиков получены в пункте 1.1.2
Для первого выхода пользуемся таблицей 1. Получившиеся графики можем сопоставить с графиками полученным в пункте 1.3.1, введя поправку на начальное значение параметра

Рисунок 11 – Реакция первого выхода на возмущения u1(t) в пункте 1.3.1

Рисунок 12 – Реакция первого выхода на возмущение для линеаризованной системы
Легко видеть, что эти график совпадают, что говорит о том, что линеаризация по первому выходу проведена на приемлемом уровне

Рисунок 14 – Реакция второго выхода на возмущения u1(t) полученного в пункте 1.3.1

Рисунок 13 – Реакция второго выхода на возмущения для линеаризованной системы
В данном случае имеет место погрешность которую можно связать с ошибкой вносимой кусочно – линейной аппроксимации.
Установившиеся состояния системы
Вычислить постоянное значение состояния системы в условиях

Т.к. установившееся значение предполагает отсутствие динамики, то систему можно записать в следующем виде


Идентификация многомерной математической модели по данным эксперимента
Активная идентификация
Для дискретной формы системы ( F , G , C ) из пункта 3. 1. провести реализацию системы.
Запишем систему в виде:


Подавая импульс по первому входу, рассчитаем:




Теперь имея экспериментальные данные, сгруппировав их в матрицы H и H 1 можем приступить к их обработки.




Из собственных векторов от (
) и (
) построим:






Для проверки идентификации найдем коэффициент передачи системы

Коэффициент передачи, вычисленный по исходным матрицам


Можно сделать вывод о том, что система идентифицирована, верно
Пассивная идентификация
Для дискретной формы системы (F, G, C) из пункта 3. 1. провести пассивную идентификацию системы, предполагая, что вектор входа изменяется соответственно таблице:
Таблица 7 Значение вектора входа для пассивной идентификации.
| Такт, n | 0 | 1 | 2 | 3 | 4 | 5 |
| U(n) | 0.01 | 0 | 0 | 0.04 | 0 | 0 |
| 0 | 0.01 | 0.02 | 0 | 0.03 | 0 |

Используя матрицы системы в дискретной форме для заданных значений вектора входа, рассчитаем значения вектора выхода

Результаты расчета сведем в таблицу:
| Такт, n | 1 | 2 | 3 | 4 | 5 | 6 |
| y(n) | 0.003935 | 0.006321 | 0.012 | 0.023 | 0.026 | 0.016 |
| -0.0026 | 0.022 | 0.053 | 0.0091 | 0.071 | 0.026 |
Используя данные эксперимента (Таблица 8) можем приступить непосредственно к определению параметров идентифицированной системы
Тогда


Для проверки идентификации найдем коэффициент передачи системы

Система идентифицирована, верно
Дата: 2019-07-30, просмотров: 263.