rooteig

Состав частоты и мощности с помощью метода собственного вектора

Описание

пример

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

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

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

Примечание

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

[f,pow] = rooteig(___,fs) возвращает вектор из частот f вычисленный в Гц. Вы предоставляете частоту дискретизации fs в Гц. Если вы задаете fs как пустой вектор [], значения по умолчанию частоты дискретизации к 1 Гц.

Примеры

свернуть все

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

rng default
n = 0:99;   
s = exp(1i*pi/2*n)+2*exp(1i*pi/4*n)+exp(1i*pi/3*n)+randn(1,100);

X = corrmtx(s,12,'mod'); 
[W,P] = rooteig(X,3)
W = 3×1

    0.7883
    1.5674
    1.0429

P = 3×1

    4.1748
    1.0572
    1.2419

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

свернуть все

Входной сигнал в виде вектора или матрицы. Если 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 Гц.

Алгоритмы

Метод собственного вектора используется rooteig совпадает с используемым peig. Алгоритм выполняет eigenspace анализ корреляционной матрицы сигнала, чтобы оценить содержимое частоты сигнала.

Различие между peig и rooteig :

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

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

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

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

| | |

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