selstruc

Выберите порядок модели для моделей ARX одно выхода

Синтаксис

nn = selstruc(v)
[nn,vmod] = selstruc(v,c)

Описание

Примечание

Используйте selstruc для систем одно выхода только. selstruc поддержки и одно вход и несколько - входные системы.

selstruc функция должна помочь выбрать структуру модели (порядок) от информации, содержавшейся в матричном v полученный как выход из arxstruc или ivstruc.

Значение по умолчанию c 'plot'. График показывает процент выходного отклонения, которое не объяснено моделью в зависимости от количества используемых параметров. Каждое значение показывает лучшее пригодное для того количества параметров. Путем нажатия в графике можно исследовать, какие порядки представляют интерес. Когда вы нажимаете Select, переменную nn экспортируется в рабочую область как оптимальная структура модели для вашего выбора количества параметров. Могут быть сделаны несколько выбора.

c = 'aic' не дает графиков, но возвращается в nn структура, которая минимизирует

$\begin{matrix} V_{\mod} = \log (V (1 + \frac{2 d}{N})) \\ = \log (V) + \frac{2 d}{N}, N ≫ d \end{matrix}$

где V функция потерь, d является общим количеством параметров в рассматриваемой структуре, и N является количеством точек данных, используемых для оценки. $\log (V) + \frac{2 d}{N}$ Информационный критерий Акэйка (AIC). Смотрите aic для получения дополнительной информации.

c = 'mdl' возвращается в nn структура, которая минимизирует критерий Минимальной длины описания (MDL) Риссэнена.

$V_{mod} = V (1 + \frac{d \log (N)}{N})$

Когда c равняется численному значению, структура, которая минимизирует $V_{\mod} = V (1 + \frac{c d}{N})$

выбран.

Выходной аргумент vmod имеет тот же формат как v, но это содержит логарифмы соответственно модифицированных критериев.

Примеры

свернуть все

Сгенерируйте комбинации Порядка Модели и оцените модель ARX Используя метод IV

Скрипт Open Live Script

Создайте наборы данных оценки и валидации

load iddata1;
ze = z1(1:150);
zv = z1(151:300);

Сгенерируйте комбинации порядка модели для оценки, указав диапазоны для порядков модели и задержек.

NN = struc(1:3,1:2,2:4);

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

V = ivstruc(ze,zv,NN);

Выберите порядок модели с лучшей подгонкой к данным о валидации.

order = selstruc(V,0);

Оцените модель ARX выбранного порядка.

M = iv4(ze,order);

Сгенерируйте комбинации Порядка Модели и оцените мультивходную модель ARX

Скрипт Open Live Script

Создайте наборы данных оценки и валидации.

load co2data;
Ts = 0.5; % Sample time is 0.5 min
ze = iddata(Output_exp1,Input_exp1,Ts);
zv = iddata(Output_exp2,Input_exp2,Ts);

Сгенерируйте комбинации порядка модели для:

na= 2:4
nb= 2:5 для первого входа и 1 или 4 для второго входа.
nk= 1:4 для первого входа и 0 для второго входа.

NN = struc(2:4,2:5,[1 4],1:4,0);

Оцените модель ARX для каждой комбинации порядка модели.

V = arxstruc(ze,zv,NN);

Выберите порядок модели с лучшей подгонкой к данным о валидации.

order = selstruc(V,0);

Оцените модель ARX выбранного порядка.

M = arx(ze,order);

Представлено до R2006a

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