ESPRIT

Пример формирования подмасивов в методе 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

Литература

Content Disclaimer

Informasi ini disarikan dari Wikipedia dan disajikan kembali untuk tujuan edukasi. Konten tersedia di bawah lisensi CC BY-SA 3.0. Kami tidak bertanggung jawab atas ketidakakuratan data yang bersumber dari kontribusi publik tersebut.

  1. The information displayed on this website is sourced in part or in whole from Wikipedia and has been adapted for the purpose of restating it. We strive to provide accurate and relevant information, however:
  2. There is no guarantee of absolute accuracy. Wikipedia is an open, collaborative project that can be edited by anyone, so information is subject to change.
  3. It is not intended to constitute professional advice. The content displayed is for informational and educational purposes only. For important decisions (e.g., medical, legal, or financial), please consult a professional.
  4. Content copyright. Wikipedia is licensed under the Creative Commons Attribution-ShareAlike License (CC BY-SA). This means that content may be reused with appropriate attribution and shared under a similar license.
  5. Responsible use. Any risk arising from the use of information from this website is entirely the responsibility of the user.