Реальная cepstrum и минимальная фазовая реконструкция
Действительный cepstrum является обратным преобразованием Фурье действительного логарифма величины преобразования Фурье последовательности.
Примечание
rceps
работает только на реальных данных.
rceps
является реализацией алгоритма 7.2 в [2], то есть
y = real(ifft(log(abs(fft(x)))));
Соответствующее оконцевание в области cepstral формирует восстановленный минимально-фазовый сигнал:
w = [1;2*ones(n/2-1,1);ones(1-rem(n,2),1);zeros(n/2-1,1)]; ym = real(ifft(exp(fft(w.*y))));
[1] Oppenheim, Alan V., and Ronald W. Schafer. Digital Signal Processing, Englewood Cliffs, NJ, Prentice Hall, 1975.
[2] Программы для цифровой обработки сигналов, IEEE Press, Нью-Йорк, 1979.