Сгенерируйте выходные данные на основе заданной модели ARX и используйте выходные данные, чтобы оценить модель.

Задайте полиномиальную модель sys0 со структурой ARX. Модель включает входную задержку одной выборки, описанной как начальный нуль в B полином.

A = [1  -1.5  0.7];
B = [0 1 0.5];
sys0 = idpoly(A,B);

Сгенерируйте измеренный входной сигнал u это содержит случайный бинарный шум и сигнал ошибки e это содержит нормально распределенный шум. С этими сигналами симулируйте измеренный выходной сигнал y из sys0.

u = iddata([],idinput(300,'rbs'));
e = iddata([],randn(300,1));
y = sim(sys0,[u e]);

Объедините y и u в один iddata объект z. Оцените новую модель ARX с помощью z и те же полиномиальные порядки и вход задерживаются как исходная модель.

z = [y,u];
sys = arx(z,[2 2 1])

sys =
Discrete-time ARX model: A(z)y(t) = B(z)u(t) + e(t)
  A(z) = 1 - 1.524 z^-1 + 0.7134 z^-2              
                                                   
  B(z) = z^-1 + 0.4748 z^-2                        
                                                   
Sample time: 1 seconds
  
Parameterization:
   Polynomial orders:   na=2   nb=2   nk=1
   Number of free coefficients: 4
   Use "polydata", "getpvec", "getcov" for parameters and their uncertainties.

Status:                                          
Estimated using ARX on time domain data "z".     
Fit to estimation data: 81.36% (prediction focus)
FPE: 1.025, MSE: 0.9846                          

Выход отображает полином, содержащий предполагаемые параметры вместе с другими деталями оценки. Под Status, Fit to estimation data показывает, что предполагаемая модель имеет 1 шаг вперед точность предсказания выше 80%.

Оцените модель AR timeseries с помощью arx функция. Модель AR не имеет никакого измеренного входа.

Загрузите данные, которые содержат временные ряды z9 с шумом.

load iddata9 z9

Оцените четвертый порядок модель AR путем определения только na закажите в [na nb nk].

sys = arx(z9,4);

Исследуйте предполагаемое полиномиальные параметры и припадок оценки к данным.

param = sys.Report.Parameters.ParVector

param = 4×1

   -0.7923
   -0.4780
   -0.0921
    0.4698

fit = sys.Report.Fit.FitPercent

fit = 79.4835

Оцените параметры модели ARIX. Модель ARIX является моделью ARX с интегрированным шумом.

Задайте полиномиальную модель sys0 со структурой ARX. Модель включает входную задержку одной выборки, описанной как начальный нуль в B.

A = [1  -1.5  0.7];
B = [0 1 0.5];
sys0 = idpoly(A,B);

Симулируйте выходной сигнал sys0 использование случайного двоичного входа сигнализирует о u и нормально распределенный сигнал ошибки e.

u = iddata([],idinput(300,'rbs'));
e = iddata([],randn(300,1));
y = sim(sys0,[u e]);

Интегрируйте выходной сигнал и сохраните результат yi в iddata объект zi.

yi = iddata(cumsum(y.y),[]);
zi = [yi,u];

Оцените модель ARIX от zi. Установите аргумент пары "имя-значение" 'IntegrateNoise' к true.

sys = arx(zi,[2 2 1],'IntegrateNoise',true);

Предскажите выход модели с помощью предсказания с 5 шагами и сравните результат с yi.

compare(zi,sys,5)

Используйте arxRegul определить константы регуляризации автоматически и использовать значения для оценки модели FIR с порядком 50.

Получите lambda и R значения.

load regularizationExampleData eData;
orders = [0 50 0];
[lambda,R] = arxRegul(eData,orders);

Используйте возвращенный lambda и R значения для упорядоченной оценки модели ARX.

opt = arxOptions;
opt.Regularization.Lambda = lambda;
opt.Regularization.R = R;
sys = arx(eData,orders,opt);

Загрузите данные.

load iddata1ic z1i

Оцените модель ARX второго порядка sys и возвратите начальные условия в ic.

na = 2;
nb = 2;
nk = 1;
[sys,ic] = arx(z1i,[na nb nk]);
ic

ic = 
  initialCondition with properties:

     A: [2x2 double]
    X0: [2x1 double]
     C: [0 2]
    Ts: 0.1000

ic initialCondition объект, который инкапсулирует свободный ответ sys, в форме пространства состояний, к вектору начального состояния в X0. Можно включить ic когда вы симулируете sys с z1i входной сигнал и сравнивает ответ с z1i выходной сигнал.