armcov

Авторегрессивные параметры модели все-полюса — измененный метод ковариации

Синтаксис

a = armcov(x,p) [a,e] = armcov(x,p)

Описание

a = armcov(x,p) использует измененный метод ковариации, чтобы соответствовать th-порядку p авторегрессивная модель (AR) к входному сигналу, x, который принят, чтобы быть выводом системы AR, управляемой белым шумом. Этот метод минимизирует прямые и обратные ошибки прогноза в смысле наименьших квадратов. Выходной массив, a, содержит нормированные оценки системных параметров AR, A (z), в убывающих степенях z. a имеет p + 1 столбец. Если x является вектором, то a является вектором - строкой. Если a является матрицей, то коэффициенты вдоль n th строка модели a n th столбец x.

[a,e] = armcov(x,p) возвращает оценку отклонения, e, белого шумового входа к модели AR.

Примеры

свернуть все

Оценка параметра Используя измененный метод ковариации

Скрипт Open Live Script

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

rng default

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

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

arcoeffs = armcov(y,4)

arcoeffs = 1×5

    1.0000   -2.7741    3.8404   -2.6841    0.9360

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

nrealiz = 50;

noisestdz = rand(1,nrealiz)+0.5;

randnoise = randn(1024,nrealiz);

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

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

Повторите процедуру с помощью многоканального синтаксиса armcov.

realiz = bsxfun(@times,noisestdz,randnoise);

Y = filter(1,A,realiz);

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

hold on
plot(noisestdz.^2,variances,'o')

q = legend('Single channel loop','Multichannel');
q.Location = 'best';

Больше о

свернуть все

AR (p) модель

Позвольте y (n) быть широким смыслом стационарный вероятностный процесс, полученный путем фильтрации белого шумового входа с отклонением e с системной функцией A (z). Если _Py (^ejω) является степенью спектральная плотность y (n), то

$P_{y} (e^{j ω}) = \frac{e}{{| A (e^{j ω}) |}^{2}} = \frac{e}{{| 1 + \sum_{k = 1}^{p} a (k) e^{- j ω k} |}^{2}} .$

Поскольку метод характеризует входные данные с помощью модели все-полюса, верный выбор порядка модели, p, важен.

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

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

Поддержка

Памятка переводчика

1. Если смысл перевода понятен, то лучше оставьте как есть и не придирайтесь к словам, синонимам и тому подобному. О вкусах не спорим.

2. Не дополняйте перевод комментариями “от себя”. В исправлении не должно появляться дополнительных смыслов и комментариев, отсутствующих в оригинале. Такие правки не получится интегрировать в алгоритме автоматического перевода.

3. Сохраняйте структуру оригинального текста - например, не разбивайте одно предложение на два.

4. Не имеет смысла однотипное исправление перевода какого-то термина во всех предложениях. Исправляйте только в одном месте. Когда Вашу правку одобрят, это исправление будет алгоритмически распространено и на другие части документации.

5. По иным вопросам, например если надо исправить заблокированное для перевода слово, обратитесь к редакторам через форму технической поддержки.

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