rooteig

Частота и содержимое степени использование метода собственного вектора

Синтаксис

[w,pow] = rooteig(x,p)
[f,pow] = rooteig(...,fs)
[w,pow] = rooteig(...,'corr')

Описание

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

  • Строка или вектор-столбец, представляющий одно наблюдение за сигналом

  • Прямоугольный массив, для которого каждая строка x представляет отдельное наблюдение за сигналом (например, каждой строкой является один вывод массива датчиков, как в обработке матриц), такой, что x'*x является оценкой корреляционной матрицы

Примечание

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

Можно задать второй входной параметр p как также:

  • Скалярное целое число. В этом случае размерностью подпространства сигнала является p.

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

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

Длина векторного w является вычисленной размерностью подпространства сигнала. Для входных данных с действительным знаком x длина соответствующего вектора степени pow дают

length(pow) = 0.5*length(w)

Для входных данных с комплексным знаком x, pow и w имеют ту же длину.

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

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

Примечание

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

Примеры

свернуть все

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

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

Алгоритмы

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

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

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

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

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

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

| | |

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