R{p) = R(—p). В этих обозначениях равенства (3.G) принимают вид системы из к уравнений: Я(1) = а1Я( 0) +ааЯ( 1) + --- + акЯ( * - 1 ) , Л(2) = ОіЯ( 1) + а2Я(0) + •■• + akR{k - 2), а последнее уравнение имеет вид R(k) —a.iR(k — 1) + a^R{k —2) + ■• ■+ akR(0). Полученную систему запишем в матричной форме. Введем обозначения: ( Я(0) Д( і ) • ■ R { k ~ 1) \ Ак = Я(1) Д(о) • R(k - 2) ^ R{k - 1) R(k - 2) • R{ 0) / ( Я(1) N / ax > fk = R{ 2) , a k = Я2 l m J \ ak J Теперь система (3.6) равносильна матричному уравнению Ak-cxk = /к• Решение можно записать в форме а к = АЦ1 Д . На практике используются лишь значения сигнала на некоторых интервалах, где сигнал можно считать стационарным. Значения автокорреляционной функции Rip) подсчитываются с помощью ДПФ. Что касается алгоритма решения системы уравнений. то существует более быстрый способ решения, не требующий обращения матрицы. Этот алгоритм носит название алгоритма Durbin'a. Он находит применение в случае, когда процедура реализуется на процессоре малой мощности. В настоящее время большинство пакетов использует обращение матрицы для нахождения коэффициентов aj.
RkJQdWJsaXNoZXIy MTExODQxMg==