rceps

Вещественный кепстр и минимально-фазовое восстановление.

Библиотека

EngeeDSP

Синтаксис

y, ym = rceps(x) — возвращает вещественный кепстр y и восстановленную минимально-фазовую версию ym входной последовательности.

# x — входной сигнал
вещественный вектор

Details

Входной сигнал, заданный как вещественный вектор.

# y — вещественный кепстр
вектор

Details

Вещественный кепстр, возвращаемый в виде вектора.

# ym — минимально-фазовый вещественный кепстр
вектор

Details

Минимально-фазовый вещественный кепстр, возвращаемый в виде вектора.

Вещественный кепстр — это обратное преобразование Фурье от вещественного логарифма модуля преобразования Фурье последовательности.

Функция rceps работает только с вещественными данными.

Функция rceps — это реализация алгоритма 7.2 из [2]:

y = real(ifft(log(abs(fft(x)))))

Соответствующее оконное преобразование в кепстральной области формирует восстановленный минимально-фазовый сигнал:

w = [1.0; fill(2.0, n÷2 - 1); fill(1.0, 1 - n % 2); zeros(n÷2 - 1)]
ym = real(ifft(exp.(fft(w .* y))))

Oppenheim, Alan V., and Ronald W. Schafer. Digital Signal Processing, Englewood Cliffs, NJ, Prentice-Hall, 1975.
Programs for Digital Signal Processing, IEEE Press, New York, 1979.