Оценка частоты методами подпространства

В этом примере показано, как разрешить близко расположенные синусоиды с помощью методов подпространства. Методы подпространства принимают гармоническую модель, состоящую из суммы синусоид, возможно объединяют в аддитивном шуме. В гармонической модели с комплексным знаком шум также с комплексным знаком.

Создайте сигнал 24 с комплексным знаком выборки в длине. Сигнал состоит из двух комплексных экпонент (синусоиды) с частотами 0,50 Гц и 0,52 Гц и аддитивным комплексным белым Гауссовым шумом. Шум имеет нулевое среднее значение и отклонение 0.22. В комплексном белом шуме и действительные и мнимые части имеют отклонение, равное 1/2 полное отклонение.

n = 0:23;
rng default
x = exp(1j*2*pi*0.5*n)+exp(1j*2*pi*0.52*n)+ ...
    0.2/sqrt(2)*(randn(size(n))+1j*randn(size(n)));

Используя periodogram, попытайтесь разрешить эти две синусоиды.

periodogram(x,rectwin(length(x)),128,1)

Периодограмма показывает размытый максимум около 1/2 Гц. Вы не можете разрешить две отдельных синусоиды, потому что разрешение частоты периодограммы равняется 1 / _ N _, где N является длиной сигнала. В этом случае 1 / _ N_ больше разделения этих двух синусоид. Нулевое дополнение не помогает разрешить два отдельных peaks.

Используйте метод подпространства, чтобы разрешить два близко расположенных peaks. В этом примере используйте метод корневой МУЗЫКИ. Оцените матрицу автокорреляции и введите матрицу автокорреляции в pmusic. Задайте модель с двумя синусоидальными компонентами. Постройте результат.

[X,R] = corrmtx(x,14,'mod');
[S,F] = pmusic(R,2,[],1,'corr');

plot(F,S,'linewidth',2)
xlim([0.46 0.60])
xlabel('Hz')
ylabel('Pseudospectrum')

Метод корневой МУЗЫКИ может разделить два peaks на уровне 0.5 и 0,52 Гц. Однако методы подпространства не производят оценки степени как степень спектральные оценки плотности. Методы подпространства являются самыми полезными для идентификации частоты и могут быть чувствительными к порядку модели misspecification.

Смотрите также

| |

Для просмотра документации необходимо авторизоваться на сайте