rootmusic

Корневой алгоритм MUSIC

Описание

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 в Гц.

Примеры

свернуть все

Оцените амплитуды для двух синусоид в шуме. Разделение между синусоидами меньше разрешения периодограммы, 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 λ, самое маленькое предполагаемое собственное значение корреляционной матрицы сигнала. Собственные значения ниже порога λ 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. Алгоритм выполняет eigenspace анализ корреляционной матрицы сигнала для того, чтобы оценить содержимое частоты сигнала.

Различие между pmusic и rootmusic :

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

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

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

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

| | |

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

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