Идентифицируйте нелинейные модели черного ящика Используя приложение System Identification

Введение

Цели

Оцените и подтвердите нелинейные модели от single-input/single-output (SISO) данные, чтобы найти тот, который лучше всего представляет вашу системную динамику.

После завершения этого примера вы сможете выполнить следующие задачи с помощью приложения System Identification:

  • Импортируйте объекты данных из рабочей области MATLAB® в приложение.

  • Оцените и подтвердите нелинейные модели из данных.

  • Постройте и анализируйте поведение нелинейности.

Описание данных

Этот пример использует файл данных twotankdata.mat, который содержит данные временного интервала SISO для системы 2D бака, показанной в следующем рисунке.

Система 2D бака

В системе 2D бака вода льется через трубопровод в Бак 1, стекает в Бак 2 и оставляет систему через маленькое отверстие в нижней части Бака 2. Измеренный вход u (t) к системе является напряжением, применился к насосу, который подает воду в Бак 1 (в вольтах). Измеренный выход y (t) является высотой воды в более низком баке (в метрах).

На основе закона Бернулли, который утверждает, что вода, текущая через маленькое отверстие в нижней части бака, зависит нелинейно на уровне воды в баке, вы ожидаете, что отношение между входом и выходными данными будет нелинейно.

twotankdata.mat включает 3 000 выборок с шагом расчета 0,2 с.

Что такое нелинейные модели черного ящика?

Типы нелинейных моделей черного ящика

Можно оценить нелинейные модели черного ящика дискретного времени и для одно выхода и для нескольких - выходные данные временного интервала. Можно выбрать из двух типов нелинейных, структур модели черного ящика:

  • Нелинейные модели ARX

  • Хаммерстайн-Винер моделирует

Примечание

Можно оценить модели черного ящика Хаммерстайна-Винера из данных о вводе/выводе только. Эти модели не поддерживают данные timeseries, где нет никакого входа.

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

Что такое нелинейная модель ARX?

Нелинейная модель ARX состоит из регрессоров модели и средства оценки нелинейности. Средство оценки нелинейности включает и линейные и нелинейные функции, которые действуют на регрессоры модели, чтобы дать выход модели. Эта блок-схема представляет структуру нелинейной модели ARX в сценарии симуляции.

Программное обеспечение вычисляет нелинейную модель ARX выход y на двух этапах:

  1. Это вычисляет значения регрессора из текущих и прошлых входных значений и прошлых выходных данных.

    В самом простом случае регрессоры являются задержанными вводами и выводами, такими как u (t-1) и y (t-3). Подобные регрессоры называются стандартными регрессорами. Вы задаете стандартные регрессоры с помощью порядков модели и задержки. Для получения дополнительной информации смотрите Нелинейные Порядки Модели ARX и Задержку. Можно также задать пользовательские регрессоры, которые являются нелинейными функциями задержанных вводов и выводов. Например, u (t-1) *y (t-3). Чтобы создать набор полиномиальных регрессоров типа, используйте polyreg.

    По умолчанию все регрессоры являются входными параметрами и к линейному и к нелинейным функциональным блокам средства оценки нелинейности. Можно выбрать подмножество регрессоров как входные параметры к нелинейному функциональному блоку.

  2. Это сопоставляет регрессоры с выходом модели с помощью блока средства оценки нелинейности. Блок средства оценки нелинейности может включать линейные и нелинейные блоки параллельно. Например:

    F(x)=LT(x−r)+d+g(Q(x−r))

    Здесь, x является вектором регрессоров, и r является средним значением регрессоров x. LT(x)+d выход блока линейной функции и является аффинным когда d ≠ 0. d является скалярным смещением. g(Q(x−r)) представляет выход нелинейного функционального блока. Q является матрицей проекции, которая делает вычисления хорошо подготовленными. Точная форма F (x) зависит от вашего выбора средства оценки нелинейности. Можно выбрать из доступных средств оценки нелинейности, таких как сети древовидного раздела, сети вейвлета и многоуровневые нейронные сети. Можно также исключить или линейное или нелинейный функциональный блок от средства оценки нелинейности.

    При оценке нелинейной модели ARX программное обеспечение вычисляет значения параметра модели, такие как L, r, d, Q и другие параметры, задающие g.

Получившимися нелинейными моделями ARX является idnlarx объекты, которые хранят все данные модели, включая регрессоры модели и параметры средства оценки нелинейности. Для получения дополнительной информации об этих объектах, смотрите Нелинейные Структуры модели.

Что такое модель Хаммерстайна-Винера?

Эта блок-схема представляет структуру модели Хаммерстайна-Винера:

Где,

  • f является нелинейной функцией, которая преобразовывает входные данные u (t) как w (t) = f (u (t)).

    w (t), внутренняя переменная, является выходом блока Input Nonlinearity и имеет ту же размерность как u (t).

  • B/F является линейной передаточной функцией, которая преобразовывает w (t) как x (t) = (B/F) w (t).

    x (t), внутренняя переменная, является выходом блока Linear и имеет ту же размерность как y (t).

    B и F похожи на полиномы в линейной модели Output-Error. Для получения дополнительной информации о моделях Output-Error, смотрите то, Что Полиномиальные Модели?.

    Для ny выходные параметры и входные параметры nu, линейный блок является матрицей передаточной функции, содержащей записи:

    Bj,i(q)Fj,i(q)

    где j = 1,2,...,ny и i = 1,2,...,nu.

  • h является нелинейной функцией, которая сопоставляет выход линейного блока x (t) к системе выход y (t) как y (t) = h (x (t)).

Поскольку действия f на входном порте линейного блока, эта функция вызвана входная нелинейность. Точно так же, потому что действия h на выходном порте линейного блока, эта функция вызвана выходная нелинейность. Если ваша система содержит несколько вводов и выводов, необходимо задать функции f и h для каждого сигнала ввода и вывода. Вы не должны включать и вход и выходную нелинейность в структуре модели. Когда модель содержит только входную нелинейность f, это называется моделью Hammerstein. Точно так же, когда модель содержит только выходную нелинейность h, это называется Винеровской моделью.

Программное обеспечение вычисляет модель Хаммерстайна-Винера выход y на трех этапах:

  1. Вычислите w (t) = f (u (t)) от входных данных.

    w (t) является входом к линейной передаточной функции B/F.

    Входная нелинейность является статической функцией (без памяти), где значение выхода данное время t зависит только от входного значения во время t.

    Можно сконфигурировать входную нелинейность как сигмоидальную сеть, сеть вейвлета, насыщение, мертвую зону, кусочную линейную функцию, одномерный полином или пользовательскую сеть. Можно также удалить входную нелинейность.

  2. Вычислите выход линейного блока с помощью w (t) и начальные условия: x (t) = (B/F) w (t).

    Можно сконфигурировать линейный блок путем определения порядков числителя B и знаменатель F.

  3. Вычислите модель, выведенную путем преобразования выхода линейного блока x (t) с помощью нелинейного функционального h в качестве y (t) = h (x (t)).

    Подобно входной нелинейности выходная нелинейность является статической функцией. Можно сконфигурировать выходную нелинейность таким же образом как входную нелинейность. Можно также удалить выходную нелинейность, такую что y (t) = x (t).

Получившимися моделями является idnlhw объекты, которые хранят все данные модели, включая средства оценки нелинейности и параметры модели. Для получения дополнительной информации об этих объектах, смотрите Нелинейные Структуры модели.

Подготовка данных

Загрузка данных в рабочее пространство MATLAB

Загрузите выборочные данные в twotankdata.mat путем ввода следующей команды в Окне Команды MATLAB:

load twotankdata

Эта команда загружает следующие две переменные в браузер рабочего пространства MATLAB:

  • u входные данные, который является напряжением, применился к насосу, который подает воду в Бак 1 (в вольтах).

  • y выходные данные, который является высотой воды в Баке 2 (в метрах).

Создание iddata Объекты

Объекты данных System Identification Toolbox™ инкапсулируют и значения данных и свойства данных в одну сущность. Можно использовать команды System Identification Toolbox, чтобы удобно управлять этими объектами данных как одна сущности.

Вы, должно быть, уже загрузили выборочные данные в рабочее пространство MATLAB, как описано в Загрузке Данных в рабочее пространство MATLAB.

Используйте следующие команды, чтобы создать два iddata объекты данных, ze и zv, где ze содержит данные для оценки модели и zv содержит данные для проверки допустимости модели. Ts шаг расчета.

Ts = 0.2; % Sample time is 0.2 sec
z = iddata(y,u,Ts);
% First 1000 samples used for estimation
ze = z(1:1000);
% Remaining samples used for validation
zv = z(1001:3000);

Просмотреть свойства iddata объект, используйте get команда. Например:

get(ze)

Программное обеспечение MATLAB возвращает следующие свойства данных и значения:

               Domain: 'Time'
                Name: ''
          OutputData: [1000x1 double]
                   y: 'Same as OutputData'
          OutputName: {'y1'}
          OutputUnit: {''}
           InputData: [1000x1 double]
                   u: 'Same as InputData'
           InputName: {'u1'}
           InputUnit: {''}
              Period: Inf
         InterSample: 'zoh'
                  Ts: 0.2000
              Tstart: 0.2000
    SamplingInstants: [1000x0 double]
            TimeUnit: 'seconds'
      ExperimentName: 'Exp1'
               Notes: {}
            UserData: []

Чтобы изменить свойства данных, используйте запись через точку. Например, чтобы присвоить названия канала и модули, которые подписывают оси графика, введите следующий синтаксис в Окне Команды MATLAB:

% Set time units to minutes
ze.TimeUnit = 'sec';
% Set names of input channels
ze.InputName = 'Voltage';
% Set units for input variables
ze.InputUnit = 'V';
% Set name of output channel
ze.OutputName = 'Height';
% Set unit of output channel
ze.OutputUnit = 'm';

% Set validation data properties
zv.TimeUnit = 'sec';
zv.InputName = 'Voltage';
zv.InputUnit = 'V';
zv.OutputName = 'Height';
zv.OutputUnit = 'm';

Проверять что InputName свойство ze изменяется, введите следующую команду:

ze.inputname

Совет

Имена свойства, такие как InputName, не являются чувствительными к регистру. Можно также сократить имена свойства, которые запускаются с Input или Output путем замены u для Input и y для Output на имя свойства. Например, OutputUnit эквивалентно yunit.

Запуск приложения System Identification

Чтобы открыть приложение System Identification, введите следующую команду в Окне Команды MATLAB:

systemIdentification

Имя сеанса по умолчанию, Untitled, появляется в строке заголовка.

Импорт объектов данных в приложение System Identification

Можно импортировать объекты данных в приложение от рабочего пространства MATLAB.

Вы, должно быть, уже создали объекты данных, как описано в Создании iddata Объекты, и открыли приложение, как описано в Запуске Приложения System Identification.

Импортировать объекты данных:

  1. В приложении System Identification выберите Import data> Data object.

    Это действие открывает диалоговое окно Import Data.

  2. Введите ze в поле Object, чтобы импортировать данные об оценке. Нажмите Enter.

    Это действие вводит информацию об объекте в Поля данных Импорта.

    Нажмите More, чтобы просмотреть дополнительную информацию об этих данных, включая названия канала и модули.

  3. Нажмите Import, чтобы добавить значок под названием ze к приложению System Identification.

  4. В диалоговом окне Import Data введите zv в поле Object, чтобы импортировать данные о валидации. Нажмите Enter.

  5. Нажмите Import, чтобы добавить значок под названием zv к приложению System Identification.

  6. В диалоговом окне Import Data нажмите Close.

  7. В приложении System Identification перетащите данные о валидации значок zv к прямоугольнику Validation Data. Данные об оценке значок ze уже определяются в прямоугольнике Working Data.

    В качестве альтернативы щелкните правой кнопкой по zv значок, чтобы открыть диалоговое окно Data/model Info. Установите флажок Use as Validation Data. Нажмите Apply и затем Close, чтобы добавить zv к прямоугольнику Validation Data.

    Приложение System Identification теперь напоминает следующую фигуру.

Оценка нелинейных моделей ARX

Оценка нелинейной модели ARX с настройками по умолчанию

В этом фрагменте примера вы оцениваете нелинейную модель ARX с помощью структуры модели по умолчанию и опций оценки.

Вы, должно быть, уже подготовили данные, как описано в Подготовка данных. Для получения дополнительной информации о нелинейных моделях ARX, смотрите то, Что Нелинейная Модель ARX?

  1. В приложении System Identification выберите Estimate> Nonlinear models.

    Это действие открывает диалоговое окно Nonlinear Models.

    Вкладка Configure уже открыта, и Model type по умолчанию является Nonlinear ARX.

    Во вкладке Regressors Input Channels и Output Channels установили Delay на 1 и набор No. of Terms к 2. Модель вывела y (t), связан с входом u (t) через следующее нелинейное авторегрессивное уравнение:

    y(t)=f(y(t−1),y(t−2),u(t−1),u(t−2))

    f является средством оценки нелинейности, выбранным в Nonlinearity выпадающий список вкладки Model Properties, и является Wavelet Network по умолчанию. Номер модулей для средства оценки нелинейности определяется к Select automatically и управляет гибкостью нелинейности — больше модулей соответствует более гибкой нелинейности.

  2. Нажмите Estimate.

    Это действие добавляет модель nlarx1 к приложению System Identification, как показано в следующем рисунке.

    Диалоговое окно Nonlinear Models отображается сводные данные информации об оценке во вкладке Estimate. Fit (%) является среднеквадратичной погрешностью между результатами измерений и симулированным выходом модели: 100% соответствуют совершенной подгонке (никакая ошибка) и 0% к модели, которая не способна к объяснению ни одного изменения выхода и только среднего уровня.

    Примечание

    Fit (%) вычисляется с помощью набора данных оценки, а не набора данных валидации. Однако выходной график модели в следующем шаге сравнивает подгонку к набору данных валидации.

  3. В приложении System Identification установите флажок Model output.

    Это действие симулирует модель с помощью данных о контроле ввода, как введено для модели и строит симулированный выход сверху выходных данных о валидации.

    Область Best Fits Выходного графика Модели показывает что соглашение между выходом модели и выводом данных валидации.

Графический вывод поперечных сечений нелинейности для нелинейных моделей ARX

Выполните следующую процедуру, чтобы просмотреть форму нелинейности как функция регрессоров на Нелинейном графике Модели ARX.

  1. В приложении System Identification установите флажок Nonlinear ARX , чтобы просмотреть поперечные сечения нелинейности.

    По умолчанию график показывает отношение между выходными регрессорами Height(t-1) и Height(t-2). Этот график показывает регулярную плоскость в следующем рисунке. Таким образом отношение между регрессорами и выходом является приблизительно линейной плоскостью.

  2. В окне Nonlinear ARX Model Plot, набор Regressor 1 к Voltage(t-1). Установите Regressor 2 на Voltage(t-2). Нажмите Apply.

    Отношение между этими регрессорами и выходом нелинейно, как показано в следующем графике.

  3. Чтобы вращать поверхность нелинейности, выберите Style> Rotate 3D и перетащите график к новой ориентации.

  4. Чтобы отобразить 1D поперечное сечение для Регрессора 1, установите Регрессор 2 на none, и нажмите Apply. Следующий рисунок показывает получившуюся величину нелинейности для Регрессора 1, который представляет переключенный временем сигнал напряжения, Voltage(t-1).

Изменение нелинейной структуры модели ARX

В этом фрагменте примера вы оцениваете нелинейную модель ARX с определенной входной задержкой и настройками нелинейности. Как правило, вы выбираете порядки модели методом проб и ошибок, пока вы не получаете модель, которая производит точную подгонку к данным.

Вы, должно быть, уже оценили нелинейную модель ARX с настройками по умолчанию, как описано в Оценке Нелинейной Модели ARX с Настройками по умолчанию.

  1. В диалоговом окне Nonlinear Models кликните по вкладке Configure и кликните по вкладке Regressors.

  2. Для Voltage введите канал, дважды кликните соответствующую ячейку Delay, введите 3, и нажмите Enter.

    Это действие обновляет список Resulting Regressors, чтобы показать Voltage(t-3) и Voltage(t-4) — условия с минимальной входной задержкой трех выборок.

  3. Нажмите Estimate.

    Это действие добавляет модель nlarx2 к приложению System Identification и обновлениям Окно вывода Модели, чтобы включать эту модель. Диалоговое окно Nonlinear Models отображается новую информацию об оценке во вкладке Estimate.

  4. В диалоговом окне Nonlinear Models кликните по вкладке Configure и выберите вкладку Model Properties.

  5. В области Number of units in nonlinear block выберите Enter и введите 6. Этот номер управляет гибкостью нелинейности.

  6. Нажмите Estimate.

    Это действие добавляет модель nlarx3 к приложению System Identification. Это также обновляет Окно вывода Модели как показано в следующем рисунке.

Выбор подмножества регрессоров в нелинейном блоке

Можно оценить нелинейную модель ARX, которая включает только подмножество стандартных регрессоров, которые входят как входные параметры в нелинейный блок. По умолчанию все стандартные и пользовательские регрессоры используются в нелинейном блоке. В этом фрагменте примера вы только включаете стандартные регрессоры.

Вы, должно быть, уже задали структуру модели, как описано в Изменении Нелинейной Структуры модели ARX.

  1. В диалоговом окне Nonlinear Models кликните по вкладке Configure и выберите вкладку Regressors.

  2. Нажмите кнопку Edit Regressors, чтобы открыть диалоговое окно Model Regressors.

  3. Снимите следующие флажки:

    • Height(t-2)

    • Voltage(t-3)

    Нажмите OK.

    Это действие исключает переключенный временем Height(t-2) и Voltage(t-3) из списка входных параметров с нелинейным блоком.

  4. Нажмите Estimate.

    Это действие добавляет модель nlarx4 к приложению System Identification. Это также обновляет Окно вывода Модели.

Определение ранее оцененной модели с различной нелинейностью

Можно оценить ряд нелинейных моделей ARX путем создания систематических изменений к структуре модели и основывать каждую новую модель на настройке ранее предполагаемой модели. В этом фрагменте примера вы оцениваете нелинейную модель ARX, которая похожа на существующую модель (nlarx3), но с различной нелинейностью.

  1. В диалоговом окне Nonlinear Models выберите вкладку Configure. Нажмите Initialize. Это действие открывает диалоговое окно Initial Model Specification.

  2. В списке Initial Model выберите nlarx3. Нажмите OK.

  3. Кликните по вкладке Model Properties.

  4. В списке Nonlinearity выберите Sigmoid Network.

  5. В поле Number of units in nonlinear block введите 6.

  6. Нажмите Estimate.

    Это действие добавляет модель nlarx5 к приложению System Identification. Это также обновляет Выходной график Модели как показано в следующем рисунке.

Выбор лучшей модели

Лучшая модель является самой простой моделью, которая точно описывает динамику.

Чтобы просмотреть информацию о лучшей модели, включая порядок модели, нелинейность и список регрессоров, щелкают правой кнопкой по значку модели по приложению System Identification.

Оценка моделей Хаммерстайна-Винера

Оценка моделей Хаммерстайна-Винера с настройками по умолчанию

В этом фрагменте примера вы оцениваете нелинейные модели Хаммерстайна-Винера с помощью структуры модели по умолчанию и опций оценки.

Вы, должно быть, уже подготовили данные, как описано в Подготовка данных. Для получения дополнительной информации о нелинейных моделях ARX, смотрите то, Что Модель Хаммерстайна-Винера?

  1. В приложении System Identification выберите Estimate> Nonlinear models, чтобы открыть диалоговое окно Nonlinear Models.

  2. Во вкладке Configure выберите Hammerstein-Wiener в списке Model type.

    Вкладка I/O Nonlinearity открыта. Средством оценки нелинейности по умолчанию является Piecewise Linear с 10 модули для Input Channels и Output Channels, который соответствует 10 точки останова для кусочной линейной функции.

  3. Выберите вкладку Linear Block, чтобы просмотреть порядки модели и ввести задержку.

    По умолчанию порядки модели и задержка линейной модели ошибки на выходе (OE) являются nb=2, nf=3, и nk=1.

  4. Нажмите Estimate.

    Это действие добавляет модель nlhw1 к приложению System Identification.

  5. В приложении System Identification установите флажок Model output.

    Это действие симулирует модель с помощью данных о контроле ввода, как введено для модели и строит симулированный выход сверху выходных данных о валидации.

    Область Best Fits Окна вывода Модели показывает соглашение между выходом модели и выводом данных валидации.

Графический вывод нелинейности и линейной передаточной функции

Можно построить нелинейность ввода/вывода и линейную передаточную функцию модели на графике Хаммерстайна-Винера.

  1. В приложении System Identification установите флажок Hamm-Wiener, чтобы просмотреть график модели Хаммерстайна-Винера.

    График отображает входную нелинейность как показано в следующем рисунке.

  2. Кликните по yNL прямоугольнику в главном фрагменте окна Model Plot Хаммерстайна-Винера.

    График обновляется, чтобы отобразить выходную нелинейность.

  3. Кликните по прямоугольнику Linear Block в главном фрагменте окна Model Plot Хаммерстайна-Винера.

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

  4. В списке Choose plot type выберите Bode. Это действие отображает Диаграмму Боде линейной передаточной функции.

Изменение входной задержки модели Хаммерстайна-Винера

В этом фрагменте примера вы оцениваете модель Хаммерстайна-Винера с определенными настройками порядка и нелинейности модели. Как правило, вы выбираете порядки модели и задержки методом проб и ошибок, пока вы не получаете модель, которая производит удовлетворительную подгонку к данным.

Вы, должно быть, уже оценили модель Хаммерстайна-Винера с настройками по умолчанию, как описано в Оценке Моделей Хаммерстайна-Винера с Настройками по умолчанию.

  1. В диалоговом окне Nonlinear Models кликните по вкладке Configure и выберите вкладку Linear Block.

  2. Для Voltage введите канал, дважды кликните соответствующую ячейку Input Delay (nk), измените значение в 3, и нажмите Enter.

  3. Нажмите Estimate.

    Это действие добавляет модель nlhw2 к приложению System Identification и Окну вывода Модели обновляется, чтобы включать эту модель, как показано в следующем рисунке.

    Область Best Fits Окна вывода Модели показывает качество nlhw2 подгонка.

Изменение средства оценки нелинейности в модели Хаммерстайна-Винера

В этом фрагменте примера вы изменяете значение по умолчанию структура модели Хаммерстайна-Винера путем изменения ее средства оценки нелинейности.

Совет

Если вы знаете, что ваша система включает насыщение или мертво-зональную нелинейность, можно задать эти специализированные средства оценки нелинейности в модели. Piecewise Linear и Sigmoid Network средства оценки нелинейности для общего приближения нелинейности.

  1. В диалоговом окне Nonlinear Models кликните по вкладке Configure.

  2. Во вкладке I/O Nonlinearity, для Voltage введите, кликните по ячейке Nonlinearity и выберите Sigmoid Network из списка. Кликните по соответствующей ячейке No. of Units и установите значение к 20.

  3. Нажмите Estimate.

    Это действие добавляет модель nlhw3 к приложению System Identification. Это также обновляет Окно вывода Модели как показано в следующем рисунке.

  4. В диалоговом окне Nonlinear Models кликните по вкладке Configure.

  5. Во вкладке I/O Nonlinearity, набор Voltage вход Nonlinearity к Wavelet Network. Это действие устанавливает No. of Units быть определенным автоматически по умолчанию.

  6. Установите Height выход Nonlinearity к One-dimensional Polynomial.

  7. Нажмите Estimate.

    Это действие добавляет модель nlhw4 к приложению System Identification. Это также обновляет Окно вывода Модели как показано в следующем рисунке.

Выбор лучшей модели

Лучшая модель является самой простой моделью, которая точно описывает динамику.

В этом примере лучшая подгонка модели была произведена в Изменении Средства оценки Нелинейности в Модели Хаммерстайна-Винера.

Документация System Identification Toolbox
Поддержка
Для просмотра документации необходимо авторизоваться на сайте