ESPRIT

From Wikipedia (Ru) - Reading time: 3 min

Пример формирования подмасивов в методе 2D ESPRIT

ESPRIT (Estimation of Signal Parameters via Rotational Invariant Techniques) — численный метод, позволяющий оценить параметры суммы синусоид на фоне шумов по серии измерений. Является одним из спектральных алгоритмов сверхразрешения, использующих спектральное разложение корреляционной матрицы сигналов.

Метод ESPRIT предложили Полрэдж, Рой и Кейлейт в 1985 г. как развитие метода MUSIC[1][2].

Описание алгоритма

[править | править код]

Входным параметром алгоритма является корреляционная матрица сигналов антенной решетки R:

где К — количество временных отсчетов;

 — вектор значений в элементе антенной решетки;

 — эрмитово-сопряженная матрица к матрице ;

 — начальный момент времени.

Для работы алгоритма ESPRIT необходимо знать количество источников сигналов. Одним из возможных способов определения количества источников является использование спектрального разложения. За счет свойств корреляционной матрицы, ее спектральное и сингулярное разложения совпадают, поэтому в дальнейшем использование спектрального разложения подразумевает возможность использования вместо него сингулярного. Спектральное разложение корреляционной матрицы позволяет получить собственные числа и вектора сигнального и шумового подпространств:[3]

где и  — матрицы сигнального и шумового подпространств;

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

Матрицы и являются матрицами сигнальных векторов первой и второй подрешеток соответственно. Матрица является решением матричного уравнения :

Углы направлений прибытия сигналов определяются выражением:[3]

где d — расстояние между элементами антенной решетки.

Пример алгоритма

[править | править код]

Алгоритм реализации метода ESPRIT в MATLAB:

function esprit(y, model_order, number_of_sources):
    m = model_order
    n = number_of_sources
    create covariance matrix R, from the noisy measurements y. Size of R will be (m-by-m).
    compute the svd of R
    [U, E, V] = svd(R)
    
    obtain the orthonormal eigenvectors corresponding to the sources
    S = U(:, 1:n)                 
      
    split the orthonormal eigenvectors in two
    S1 = S(1:m-1, :) and S2 = S(2:m, :)
                                               
    compute P via LS (MATLAB's backslash operator)
    P = S1\S2 
       
    find the angles of the eigenvalues of P
    w = angle(eig(P))               
    return w

Применение

[править | править код]

Примечания

[править | править код]
  1. Paulraj, A.; Roy, R.; Kailath, T. (1985), "Estimation Of Signal Parameters Via Rotational Invariance Techniques - Esprit", Nineteenth Asilomar Conference on Circuits, Systems and Computers, pp. 83—89, doi:10.1109/ACSSC.1985.671426, ISBN 978-0-8186-0729-5
  2. Roy, R., & Kailath, T. (1989). Esprit-estimation of signal parameters via rotational invariance techniques. IEEE Transactions on Acoustics, Speech, and Signal Processing, 37(7), 984—995. https://doi.org/10.1109/29. 32276
  3. 1 2 Тагаев Т. И., Музафаров Д. А., Мохорт Д. Ю., Космачев В. М. Определение оптимальных тактико-технических характеристик антенной решетки для различения алгоритмом сверхразрешения ESPRIT целей, разнесенных на малое угловое расстояние // Антенны и распространение радиоволн. — 2021.
  4. Volodymyr Vasylyshyn. Direction of arrival estimation using ESPRIT with sparse arrays.// Proc. 2009 European Radar Conference (EuRAD). — 30 Sept.-2 Oct. 2009. — Pp. 246—249. — [1]
  5. Василишин В. И. Спектральный анализ методом ESPRIT при предварительной обработке данных методом SSA.// Системи обробки інформації.- 2015. -№ 15. — С. 12 −15. [2] Архивная копия от 25 марта 2022 на Wayback Machine

Литература

[править | править код]

Licensed under CC BY-SA 3.0 | Source: https://ru.wikipedia.org/wiki/ESPRIT
13 views |
↧ Download this article as ZWI file
Encyclosphere.org EncycloReader is supported by the EncyclosphereKSF