Документация

Оцените модель Используя nlarx

Используйте nlarx и создать и оценить нелинейную модель ARX. После каждой оценки подтвердите модель путем сравнения его с другими моделями и симуляции или предсказания ответа модели.

Основная оценка

Начните с самой простой оценки с помощью m = nlarx(data,[na nb nk]). Например:

load iddata1;
% na = nb = 2 and nk = 1 
m = nlarx(z1,[2 2 1])

m =
Nonlinear ARX model with 1 output and 1 input
 Inputs: u1
 Outputs: y1
 Standard regressors corresponding to the orders:
   na = 2, nb = 2, nk = 1
 No custom regressor
 Nonlinear regressors:
   y1(t-1)
   y1(t-2)
   u1(t-1)
   u1(t-2)
 Nonlinearity: wavenet with 1 unit

Sample time: 0.1 seconds

Status:                                          
Estimated using NLARX on time domain data "z1".  
Fit to estimation data: 68.83% (prediction focus)
FPE: 1.975, MSE: 1.885

Второй входной параметр [na nb nk] задайте порядки модели и задержки. По умолчанию средство оценки нелинейности является сетью вейвлета (см. wavenet страница с описанием), который берет все стандартные регрессоры в качестве входных параметров к ее линейным и нелинейным функциям. m idnlarx объект.

Для систем MIMO, nb, nf, и nk ny-by-nu матрицы. Смотрите nlarx страница с описанием для получения дополнительной информации об оценке MIMO.

Задайте различное средство оценки нелинейности (например, сигмоидальная сеть).

M = nlarx(z1,[2 2 1],'sigmoid');

Создайте nlarxOptions набор опции и конфигурирует Focus свойство минимизировать ошибку симуляции.

opt = nlarxOptions('Focus','simulation');
M = nlarx(z1,[2 2 1],'sigmoid',opt);

Сконфигурируйте регрессоры модели

Стандартные регрессоры

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

Пользовательские регрессоры

Исследуйте включая пользовательские регрессоры в нелинейной структуре модели ARX. Пользовательские регрессоры в дополнение к стандартным регрессорам модели (см. Нелинейные Порядки Модели ARX и Задержку).

Используйте polyreg или customreg создать пользовательские регрессоры в терминах переменных ввода - вывода модели. Можно задать пользовательские регрессоры с помощью CustomRegressors свойство idnlarx класс или addreg добавлять пользовательские регрессоры к существующей модели.

Например, сгенерируйте регрессоры как полиномиальные функции вводов и выводов:

load iddata1
m = nlarx(z1,[2 2 1],'sigmoidnet');
getreg(m) % displays all regressors

Regressors:
    y1(t-1)
    y1(t-2)
    u1(t-1)
    u1(t-2)

% Generate polynomial regressors up to order 2:
reg = polyreg(m)

4x1 array of Custom Regressors with fields: Function, Arguments, Delays, Vectorized

Добавьте полиномиальные регрессоры к CustomRegressors.

m = addreg(m,reg);
getreg(m)

Regressors:
    y1(t-1)
    y1(t-2)
    u1(t-1)
    u1(t-2)
    y1(t-1).^2
    y1(t-2).^2
    u1(t-1).^2
    u1(t-2).^2

m теперь включает полиномиальные регрессоры.

Можно также задать произвольные функции переменных ввода и вывода. Например:

load iddata1
m = nlarx(z1,[2 2 1],'sigmoidnet','CustomReg',{'y1(t-1)^2','y1(t-2)*u1(t-3)'});
getreg(m) % displays all regressors

Regressors:
    y1(t-1)
    y1(t-2)
    u1(t-1)
    u1(t-2)
    y1(t-1)^2
    y1(t-2)*u1(t-3)

Добавьте полиномиальные регрессоры к CustomRegressors.

m = addreg(m,reg);
getreg(m) % polynomial regressors

Regressors:
    y1(t-1)
    y1(t-2)
    u1(t-1)
    u1(t-2)
    y1(t-1)^2
    y1(t-2)*u1(t-3)
    y1(t-1).^2
    y1(t-2).^2
    u1(t-1).^2
    u1(t-2).^2

Управляйте пользовательскими регрессорами с помощью CustomRegressors свойство idnlarx класс. Например, чтобы получить указатель на функцию первого пользовательского регрессора в массиве:

CReg1 = m.CustomReg(1).Function;

Просмотрите выражение регрессора.

m.CustomReg(1).Display

ans = 
'y1(t-1)^2'

Можно исключить все стандартные регрессоры и использовать только пользовательские регрессоры в структуре модели установкой na=nb=nk=0:

m = nlarx(z1,[0 0 0],'linear','CustomReg',{'y1(t-1)^2','y1(t-2)*u1(t-3)'});

В усовершенствованных приложениях можно задать усовершенствованные опции оценки для средств оценки нелинейности. Например, wavenet и treepartition классы обеспечивают Options свойство для установки таких опций оценки.

Задайте линейные и нелинейные регрессоры

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

m = nlarx(z1,[2 2 1]); % all standard regressors are 
                        % inputs to the nonlinear function
getreg(m); % lists all standard regressors

Regressors:
    y1(t-1)
    y1(t-2)
    u1(t-1)
    u1(t-2)

m = nlarx(z1,[4 4 1],sigmoidnet,'nlreg',[5 6 7 8]);

Этот пример использует getreg определить индекс каждого регрессора из полного списка всех регрессоров модели. Только регрессор числа 5 - 8 является входными параметрами к нелинейной функции - getreg показывает, что эти регрессоры являются функциями входной переменной u1. nlreg сокращение от NonlinearRegressors свойство idnlarx класс. В качестве альтернативы включайте только входные регрессоры в нелинейное использование функционального блока:

m = nlarx(z1,[4 4 1],sigmoidnet,'nlreg','input');

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

m = nlarx(z1,[4 4 1],sigmoidnet,'nlreg','search');

Этот поиск обычно занимает много времени. Можно отобразить поисковое использование прогресса:

opt = nlarxOptions('Display','on');
m = nlarx(z1,[4 4 1],sigmoidnet,'nlreg','search',opt);

После оценки используйте m.NonlinearRegressors просмотреть, какие регрессоры были выбраны автоматическим поиском регрессора.

m.NonlinearRegressors

ans = 1×4

     3     5     6     7

Сконфигурируйте средство оценки нелинейности

Задайте средство оценки нелинейности непосредственно в команде оценки как:

m = nlarx(z1, [2 2 1],'sigmoidnet');

m = nlarx(z1,[2 2 1],wavenet('num',5));

Эта оценка использует нелинейную модель ARX с нелинейностью вейвлета, которая имеет 5 модулей. Создать объект нелинейности прежде если это как вход к средству оценки нелинейности:

w = wavenet('num', 5);
m = nlarx(z1,[2 2 1],w);
% or
w = wavenet;
w.NumberOfUnits = 5;
m = nlarx(z1,[2 2 1],w);

Для систем MIMO можно задать различную нелинейность для каждого выхода. Например, чтобы задать sigmoidnet для первого выхода и wavenet для второго выхода:

load iddata1 z1
load iddata2 z2
data = [z1, z2(1:300)];
M = nlarx(data,[[1 1;1 1] [2 1;1 1] [2 1;1 1]],[sigmoidnet wavenet]);

Если вы хотите ту же нелинейность для всех выходных каналов, задаете одну нелинейность.

M = nlarx(data,[[1 1;1 1] [2 1;1 1] [2 1;1 1]],sigmoidnet);

Следующая таблица обобщает значения, которые задают средства оценки нелинейности.

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

Задайте многослойное использование нейронной сети:

m = nlarx(data,[na nb nk],NNet)

где NNet объект нейронной сети, что вы создаете использование программного обеспечения Deep Learning Toolbox™. Смотрите neuralnet страница с описанием.

Задайте пользовательскую сеть путем определения функции под названием gaussunit.m, как описано в customnet страница с описанием. Задайте пользовательский сетевой объект CNetw и оценка модель:

CNetw = customnet(@gaussunit);
m = nlarx(data,[na nb nk],CNetw)

SNL = sigmoidnet('LinearTerm','off');
m = nlarx(z1,[2 2 1],SNL);

m = nlarx(z1,[2 2 1],'linear');

lin_m = arx(z1,[2 2 1]);

M = nlarx(z1,[0 0 0],'linear','CustomReg',{'y1(t-1)^2','y1(t-2)*u1(t-3)'});

Итеративно совершенствуйте модель

Если ваша структура модели включает нелинейность, которая поддерживает итеративный поиск (см., Задают Опции Оценки для Нелинейных Моделей ARX), можно использовать nlarx совершенствовать параметры модели:

m1 = nlarx(z1,[2 2 1],'sigmoidnet');
m2 = nlarx(z1,m1); % can repeatedly run this command

Можно также использовать pem совершенствовать исходную модель:

m2 = pem(z1,m1);

Проверяйте поисковый критерий завершения m.Report.Termination.WhyStop . Если WhyStop указывает, что оценка достигла максимального количества итераций, попытайтесь повторить оценку и возможно задать большее значение для nlarxOptions.SearchOptions.MaxIterations опция оценки:

opt = nlarxOptions;
opt.SearchOptions.MaxIterations = 30;
m2 = nlarx(z1,m1,opt); % runs 30 more iterations
                               % starting from m1

Когда m.Report.Termination.WhyStop значением является Near (local) minimum, (norm( g) < tol или No improvement along the search direction with line search , подтвердите свою модель, чтобы видеть, соответствует ли эта модель соответственно данным. В противном случае решение может застрять в локальном минимуме поверхности функции стоимости. Попытайтесь настроить SearchOptions.Tolerance значение или SearchMethod опция в nlarxOptions наборе опции и повторение оценка.

Можно также попытаться тревожить параметры последней модели с помощью init (названная рандомизация) и совершенствование модели с помощью nlarx:

M1 = nlarx(z1, [2 2 1], 'sigm'); % original model
M1p = init(M1);	% randomly perturbs parameters about nominal values
M2 = nlarx(z1, M1p);  % estimates parameters of perturbed model

Можно отобразить прогресс итеративного поиска в Окне Команды MATLAB с помощью nlarxOptions.Display опция оценки:

opt = nlarxOptions('Display','on');
M2= nlarx(z1,M1p,opt);

Диагностируйте оценку

Если вы не получаете удовлетворительную модель после многих испытаний с различными структурами модели и настроек алгоритма, возможно, что данные плохи. Например, ваши данные могут пропускать важные переменные ввода или вывода и не достаточно отвечают на все рабочие вопросы системы.

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

Используйте nlarx, чтобы Оценить Нелинейные Модели ARX

В этом примере показано, как использовать nlarx оценить нелинейную модель ARX для измеренных данных о вводе/выводе.

Подготовьте данные к оценке.

load twotankdata
z = iddata(y, u, 0.2);
ze = z(1:1000); zv = z(1001:3000);

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

m1 = nlarx(ze,[2 2 1]);
m2 = nlarx(ze,[2 2 3]);
m3 = nlarx(ze,[2 2 3],wavenet('num',8));
m4 = nlarx(ze,[2 2 3],wavenet('num',8),...
                      'nlr', [1 2]);

Альтернативный способ выполнить оценку состоит в том, чтобы сконфигурировать структуру модели сначала, и затем оценить эту модель.

m5 = idnlarx([2 2 3],sigmoidnet('num',14),'nlr',[1 2]);
m5 = nlarx(ze,m5);

Сравните получившиеся модели путем графического вывода выходных параметров модели с измеренным выходом.

compare(zv, m1,m2,m3,m4,m5)