polyreg

Полномочия и продукты стандартных регрессоров

Синтаксис

R = polyreg(model)
R = polyreg(model,'MaxPower',n)
R = polyreg(model,'MaxPower',n,'CrossTerm',CrossTermVal)

Описание

R = polyreg(model) создает массив R из полиномиальных регрессоров до степени 2. Если порядок модели ввел u и выход yнет данных=nb=2, и задержка nk=1, полиномиальными регрессорами является y (t −1) 2, u (t −1) 2, y (t −2) 2, u (t −2) 2. model idnlarx объект. Необходимо добавить эти регрессоры в model путем присвоения CustomRegressors model свойство или при помощи addreg.

R = polyreg(model,'MaxPower',n) создает массив R из полиномиальных регрессоров до степени n. Исключает сроки полномочий 1 и перекрестные условия, такие как y (t −1) *u (t−1).

R = polyreg(model,'MaxPower',n,'CrossTerm',CrossTermVal) создает массив R из полиномиальных регрессоров до степени n и включает перекрестные условия (продукты регрессоров стандартов) когда CrossTermVal 'on'. По умолчанию, CrossTermVal 'off'.

Примеры

свернуть все

Оцените нелинейную модель ARX с na=2, nb=2, и nk=1, и средство оценки нелинейности wavenet.

load iddata1
m = nlarx(z1,[2 2 1]);

Создайте полиномиальные регрессоры.

R = polyreg(m);

Оцените модель.

m = nlarx(z1,[2 2 1],'wavenet','CustomReg',R);

Просмотрите все регрессоры модели (стандартный и пользовательский).

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

Оцените нелинейную модель ARX с na=2, nb=1, и nk=1, и средство оценки нелинейности wavenet.

load iddata1
m = nlarx(z1,[2 1 1]);

Создайте полиномиальные регрессоры.

R = polyreg(m,'MaxPower',3,'CrossTerm','on')
16x1 array of Custom Regressors with fields: Function, Arguments, Delays, Vectorized

Если модель m имеет три стандартных регрессора aB и c, затем R включает условия a2, b2, c2, ab, ac, bc, a2b, a2c, ab2, abc, ac2, b2c, bc2, a3, b3, и c3.

Оцените модель.

m = nlarx(z1,[2 1 1],'wavenet','CustomReg',R);
Представленный в R2007a