Моделирование систем и процессов нефтепереработки и нефтехимии.Модель для имитации динамических рядов многомерного объекта, страница 2

Подберем оптимальные настройки блока спектра сигнала и построим графики исследовательского сигнала:

Рисунок 16 – Окно настройки блоков спектра сигнала

Рисунок 17 – Спектры входных сигналов

Рисунок 18 – График вход/выходпервого канала без помех (1000c)

Рисунок 19– График вход/выход второго канала без помех (1000 с)

Для первого выхода размах примерно равен 0,9 и второго выхода размах примерно равен 2,4. Мощность шума при котором размах помехи будет занимать 50% от размаха амплитуды сигнала без помехи то есть в нашем случае 0,45 и 1,2 соответсвенно.

Рисунок 20 – мощность белого шума для 50%

Рисунок 20 – График вход/выход первогоканала с помехой(0.5%)

Рисунок 21 – График вход/выход второго канала с помехой(0.5%)

Рисунок 22 – График вход/выход первого канала с помехой (5%)

Рисунок 23 – График вход/выход второго канала с помехой (5%)

Рисунок 24 – График вход/выход первого канала с помехой(50%)

Рисунок 25 – График вход/выход второго канала с помехой (50%)

Рисунок 26 – График выходного сигнала без помех

Рисунок 27 – График выходного сигнала с помехой (0.5%)

Рисунок 28 – График выходного сигнала с помехой (5%)

Рисунок 29 – График выходного сигнала с помехой (50%)


Импортируем следующие динамические ряды в SystemIdentification, где x0 иy0, x4иy4, x8иy8 – без помех; x1иy1,x5иy5,x9иy9 – 0.5% помеха; x2иy2, x6иy6,x10иy10 – 5% помеха; x3иy3,x7иy7,x11иy11 – 50% помеха.

Рисунок 30 – Динамические ряды

Вывод к первой части курсовой работы:

В первой части курсовой работы мы научились формировать модель для имитации динамических рядов многомерного объекта, обеспечивающих максимальную точность идентификации в условиях действия случайных помех. Сформировали 12 динамических рядов размерностью 500, 5000, 50000 (с разрешения преподавателя можно уменьшить размерность на порядок) без случайного шума и с шумом величиной 0.5, 5, 50%. Полученные динамические ряды идентичны. Шум на входное воздействие не влияет.


Задание 2: Исследование точности идентификации динамического ряда без помех с помощью ProcessModels и функций N4SID и PEM.

Импортируем динамический ряд размерностью 100000 в SystemIdentification

Проведем идентификцию параметров модели известной структуры через ProcessModels.

a)                                                     b)

Рисунок 18 – точность идентификации в processmodels

a)  y1        b) y2

Далее проведем предварительную идентификацию всего динамического ряда с помощью функции N4SID и найдем наиболее приемлемый порядок модели.

Рисунок 19 – Выбор порядка модели

Наиболее подходящим порядком модели в данном случае является 20-й порядок

Рисунок 20 – Точность идентификации функции  при 4,8,12,16,20 порядке

Рисунок 21 – точности идентификации по обоим каналам

Идентификация параметров модели неизвестной структуры функциями и N4SID в пространстве состояний: выбираем структуру, для этого на примере базового динамического ряда исследуем точность идентификации для порядков модели , уточнение точности вблизи наилучшего порядка 6,7,8,9,10 - для подтверждения уточненного лучшего порядка

Рисунок 22 – уточнение УЛП

Рисунок 23 – график точностей для 6,7,8,9,10 порядков по обоим каналам

Рисунок 21 – Окно System Identification

Рисунок 23 – Идентификация для УЛП с помощью функции PEM

Таблица 1 -Результаты идентификации