Параметры авторегрессивной модели все-полюса — метод ковариации
Используйте вектор с полиномиальными коэффициентами, чтобы сгенерировать 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 = arcov(y,4)arcoeffs = 1×5
1.0000 -2.7746 3.8419 -2.6857 0.9367
Сгенерируйте 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)] = arcov(y,4); end plot(noisestdz.^2,noisevar,'*') title('Noise Variance') xlabel('Input') ylabel('Estimated')
Повторите процедуру с помощью многоканального синтаксиса функции.
Y = filter(1,A,noisestdz.*randnoise); [coeffs,variances] = arcov(Y,4); hold on plot(noisestdz.^2,variances,'o') hold off legend('Single channel loop','Multichannel','Location',"best")
