Оцените модель AR и сравните ее ответ с измеренным выходом.

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

load iddata9 z9

Оцените четвертый порядок модель AR.

sys = ar(z9,4)

sys =
Discrete-time AR model: A(z)y(t) = e(t)                            
  A(z) = 1 - 0.8369 z^-1 - 0.4744 z^-2 - 0.06621 z^-3 + 0.4857 z^-4
                                                                   
Sample time: 0.0039062 seconds
  
Parameterization:
   Polynomial orders:   na=4
   Number of free coefficients: 4
   Use "polydata", "getpvec", "getcov" for parameters and their uncertainties.

Status:                                                
Estimated using AR ('fb/now') on time domain data "z9".
Fit to estimation data: 79.38%                         
FPE: 0.5189, MSE: 0.5108                               

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

Можно найти дополнительную информацию о результатах оценки путем исследования отчета оценки, sys.Report. Например, можно получить ковариацию параметра.

covar = sys.Report.Parameters.FreeParCovariance

covar = 4×4

    0.0015   -0.0015   -0.0005    0.0007
   -0.0015    0.0027   -0.0008   -0.0004
   -0.0005   -0.0008    0.0028   -0.0015
    0.0007   -0.0004   -0.0015    0.0014

Для получения дополнительной информации о просматривании отчета оценки см. Отчет Оценки.

Учитывая синусоидальный сигнал с шумом, сравните спектральные оценки метода Бурга с найденными использованием прямого обратного подхода.

Сгенерируйте выходной сигнал и преобразуйте его в iddata объект.

y = sin([1:300]') + 0.5*randn(300,1);
y = iddata(y);

Оцените четвертый порядок модели AR с помощью метода Бурга и с помощью прямого обратного подхода по умолчанию. Постройте спектры модели вместе.

sys_b = ar(y,4,'burg');
sys_fb = ar(y,4);
spectrum(sys_b,sys_fb)
legend('Burg','Forward-Backward')

Эти два ответа соответствуют тесно в большей части частотного диапазона.

Оцените модель ARI, которая включает интегратор в источник шума.

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

load iddata9 z9

Интегрируйте выходной сигнал.

y = cumsum(z9.y);

Оцените модель AR с 'IntegrateNoise' установите на true. Используйте метод наименьших квадратов 'ls'. Поскольку y вектор и не iddata возразите, задайте Ts.

Ts = z9.Ts;
sys = ar(y,4,'ls','Ts',Ts,'IntegrateNoise',true);

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

compare(y,sys,5)

Измените опции по умолчанию для AR функция.

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

load iddata9 z9

Измените опции по умолчанию так, чтобы функция использовала 'ls' приблизьтесь и не оценивайте ковариацию.

opt = arOptions('Approach','ls','EstimateCovariance',false)

opt = 
Option set for the ar command:

              Approach: 'ls'
                Window: 'now'
            DataOffset: 0
    EstimateCovariance: 0
               MaxSize: 250000

Description of options

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

sys = ar(z9,4,opt);

Получите отражательные коэффициенты и функции потерь при использовании метода Бурга.

Основанные на решетке подходы такой, как метод Бурга 'burg' и геометрическая решетка 'gl', вычислите отражательные коэффициенты и соответствующие значения функции потерь как часть процесса оценки. Используйте второй выходной аргумент, чтобы получить эти значения.

Сгенерируйте выходной сигнал и преобразуйте его в iddata объект.

y = sin([1:300]') + 0.5*randn(300,1);
y = iddata(y);

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

[sys,refl] = ar(y,4,'burg');
refl

refl = 2×5

         0   -0.3562    0.4430    0.5528    0.2385
    0.8494    0.7416    0.5960    0.4139    0.3904