Так как считается КЕПСТР?

[Dims]

В Википедии написано, что так:

signal → FT → abs() → square → log → FT → abs() → square → power cepstrum

(то есть, два раза прямое Фурье)

В справке Матлаба написано, что так:

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

то есть, без модуля и квадратов, а Фурье сперва прямое, а потом обратное.

По логике, должно быть два раза прямое Фурье, так как идея в том, что бы гармоники на спектре распознать.

Но формулы прямого и обратного Фурье очень похожи, а я всё подзабыл и запутался -- не могу понять. То ли имеются в виду немного разные преобразования? То ли они одинаковы просто по-разному записаны?

Пасиб.

[venco]

ЕМНИП, прямое и обратное преобразования Фурье отличаются комплексной сопряжённостью и, возможно, множителем. Сопряжённость не важна, т.к. после каждого преобразования идёт abs(), ну а множитель может быть и одинаковым, зависит от реализации.

[ewert]

venco в сообщении #299536 писал(а):

, ну а множитель может быть и одинаковым, зависит от реализации.

В Матлабе, насколько помню, реализация стандартна, т.е. симметрична. Т.е. прямое и обратное преобразования различаются исключительно комплексным сопряжением. Т.е. и там, и там множитель -- это единица делить на корень из двух пи.

(вот только в каком из двух преобразований у них плюс, а в каком минус -- в упор не помню... Но это, в конце концов, можно и методом научного тыка определить)

[Pirit]

Подскажите пожалуйста.
Есть сигнал на 1 секунду, 11264 отсчета(значения).
Разбиение делаем по фреймам 512 отсчетов.
Модуль и логарифм от модуля при нахождении спектра от чего брать?
Модуль от каждого значения во фрейме, и логарифм от модуля каждого значения во фрейме?