Exponenta Banner
exponenta event banner

rootmusic

Алгоритм Root MUSIC

Свернуть все на странице

Синтаксис

w = rootmusic(x,p)
[w,pow] = rootmusic(x,p)
[w,pow] = rootmusic(___,'corr')
[f,pow] = rootmusic(___,fs)

Описание

w = rootmusic(x,p) оценивает содержимое частоты в входном сигнале x и возвращает w, вектор частот в рад/выборке. Можно задать подпространство сигнала размерности используя входной параметр p.

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

[w,pow] = rootmusic(x,p) возвращает вектор частот w и соответствующую степень сигнала в векторе pow.

пример

[w,pow] = rootmusic(___,'corr') форсирует входной параметр x интерпретируется как корреляционная матрица, а не как матрица сигнальных данных. Для этого синтаксиса x должна быть квадратной матрицей, и все ее собственные значения должны быть неотрицательными. Этот синтаксис может включать входные параметры из предыдущего синтаксиса.

Примечание

Можно разместить 'corr' где угодно после p.

[f,pow] = rootmusic(___,fs) возвращает вектор частот f рассчитывается в Гц. Вы задаете частоту дискретизации fs в Гц.

Примеры

свернуть все

Открыть Live Script

Оцените амплитуды для двух синусоидов в шуме. Разделение между синусоидами меньше, чем разрешение периодограммы, 2π/N радианы/выборка. Используйте автокорреляционную матрицу как вход для rootmusic.

rng default
n = (0:99)';
frqs = [pi/4 pi/4+0.06];

s = 2*exp(1j*frqs(1)*n)+1.5*exp(1j*frqs(2)*n)+ ...
    0.5*randn(100,1)+1j*0.5*randn(100,1);

[~,R] = corrmtx(s,12,'mod');
[W,P] = rootmusic(R,2,'corr')
W = 2×1

    0.7946
    0.8917


P = 2×1

    4.1535
    0.7797

Входные параметры

свернуть все

Входной сигнал, заданный как вектор или матрица. Если x является вектором, затем рассматривается как одно наблюдение сигнала. Если x является матрицей, каждая строка x представляет собой отдельное наблюдение сигнала. Например, каждая строка является одним выходом массива датчиков, как при обработке массивов, таким что x'*x является оценкой корреляционной матрицы.

Для комплексных входных данных x, pow и w имеют одинаковую длину. Для реальных входных данных x, длину соответствующего вектора степени pow является 0.5*length(w).

Примечание

Можно использовать выход corrmtx чтобы сгенерировать такой массив x.

Поддержка комплексного числа: Да

Размерность Подпространства, заданная как действительное положительное целое число или двухэлементный вектор. Если p является действительным положительным целым числом, тогда это рассматривается как размерность подпространства. Если p является двухэлементным вектором, вторым элементом p представляет порог, который умножается на λ мин, наименьшее оцененное собственное значение матрицы корреляции сигнала. Собственные значения ниже порога λ min * p(2) назначаются подпространству шума. В этом случае p(1) задает максимальную размерность подпространства сигнала. Дополнительный пороговый параметр во второй записи в p обеспечивает вам больше гибкости и управления в назначении шумовых и сигнальных подпространств.

Частота дискретизации, заданная как положительная скалярная величина. Можно задать частоту дискретизации fsв Гц. Если вы задаете fs как пустой вектор [], частота дискретизации по умолчанию равна 1 Гц.

Выходные аргументы

свернуть все

Выходные частоты в рад/сэмпле, возвращенные как вектор. Длина вектора w - вычисленная размерность подпространства сигнала.

Степень сигнала, возвращенная как вектор.

Выходные частоты в Гц, возвращенные как вектор. Вы задаете частоту дискретизации fs в Гц. Если вы задаете fs с пустым вектором [] частота дискретизации по умолчанию равна 1 Гц.

Совет

Если входной сигнал x реально, и нечетное число синусоидов определяется pотобразится сообщение об ошибке:

Real signals require an even number p of complex sinusoids.

Алгоритмы

Алгоритм классификации нескольких сигналов (MUSIC), используемый rootmusic является тем же самым, что и используемый pmusic. Алгоритм выполняет собственный пространственный анализ матрицы корреляции сигнала в порядок для оценки частотного содержимого сигнала.

Переменное различие между pmusic и rootmusic является:

  • pmusic возвращает псевдоспектр на всех частотных выборках.

  • rootmusic возвращает оцененный дискретный частотный спектр вместе с соответствующими оценками степени сигнала.

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

Расширенные возможности

.

См. также

| | |

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

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

Поддержка