Используйте вектор полиномиальных коэффициентов для генерации процесса AR (4) путем фильтрации 1024 выборок белого шума. Сбросьте генератор случайных чисел для воспроизводимых результатов. Для оценки коэффициентов используется метод Юле-Уокера.

rng default

A = [1 -2.7607 3.8106 -2.6535 0.9238];

y = filter(1,A,0.2*randn(1024,1));

arcoeffs = aryule(y,4)

arcoeffs = 1×5

    1.0000   -2.7262    3.7296   -2.5753    0.8927

Генерируют 50 реализаций процесса, каждый раз изменяя дисперсию входного шума. Сравните оцененные по Юле-Уокеру отклонения с фактическими значениями.

nrealiz = 50;

noisestdz = rand(1,nrealiz)+0.5;

randnoise = randn(1024,nrealiz);
noisevar = zeros(1,nrealiz);

for k = 1:nrealiz
    y = filter(1,A,noisestdz(k) * randnoise(:,k));
    [arcoeffs,noisevar(k)] = aryule(y,4);
end

plot(noisestdz.^2,noisevar,'*')
title('Noise Variance')
xlabel('Input')
ylabel('Estimated')

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

Y = filter(1,A,noisestdz.*randnoise);

[coeffs,variances] = aryule(Y,4);

hold on
plot(noisestdz.^2,variances,'o')
hold off
legend('Single channel loop','Multichannel','Location','best')

Используйте вектор полиномиальных коэффициентов для генерации процесса AR (2) путем фильтрации 1024 выборок белого шума. Сбросьте генератор случайных чисел для воспроизводимых результатов.

rng default

y = filter(1,[1 -0.75 0.5],0.2*randn(1024,1));

Используйте метод Юле-Уокера, чтобы подогнать модель AR (10) к процессу. Выведите и постройте график коэффициентов отражения. Только первые два коэффициента лежат за пределами 95% доверительных границ, что указывает на то, что модель AR (10) значительно переоценивает временную зависимость в данных. Для получения дополнительной информации см. Выбор заказа AR с частичной автокорреляционной последовательностью.

[ar,nvar,rc] = aryule(y,10);

stem(rc)
xlim([0 11])
conf95 = sqrt(2)*erfinv(0.95)/sqrt(1024);
[X,Y] = ndgrid(xlim,conf95*[-1 1]);
hold on
plot(X,Y,'--r')
hold off
title('Reflection Coefficients')