|
ru.algorithms
- RU.ALGORITHMS ----------------------------------------------------------------
From : Evgenij Masherov 2:5020/175.2 03 Apr 2003 09:42:52
To : Alex Astafiev
Subject : Minimum-phase фильтры
--------------------------------------------------------------------------------
Mon Mar 31 2003 04:11, Alex Astafiev wrote to Aleksey Vaneev:
AV>> Кто-нибудь может объяснить конкретнее, что мне нужно сделать:
AV>>
AV>> 1. Create an equiripple filter, that will meet |H(z)|^2.
AV>> 2. Zero-pad the FIR. Four times the filter length works well.
AV>> 3. Calculate the DFT (or FFT) of the above.
AV>> 4. Calculate 0.25 * log(|H(k)|^2).
AV>> 5. Calculate IDFT (with scaling factor) of the results.
AV>> 6. Multiply pointwise by the homomorphic filter lmin[n] = 2u[n] -
AV>> d[n], where d[n] = dirac delta function. 7. Calculate DFT of results.
AV>> 8. Calculate complex exponential of results. 9. Calculate IDFT (with
AV>> scaling factor) of the results. 10. Keep half (or half + 1 if orig
AV>> length was odd) of the results as the minimum phase filter.
AV>>
AV>> Совершенно не могу понять пункт 6.
AV>> Остальное, вроде бы, очевидно. Тем
AV>> более, что Евгений Машеров объяснил что есть log() и exp() от
AV>> комплексного спектра.
AA> Интересно, можно как нибудь
AA> 1. Синтезируем АЧХ
AA> 2. Перевод Im+Re -> Mag + Phase
AA> 3. Phase = 0.
AA> 4. перевод Mag + Phase -> Im+Re
AA> 5. Сдвижка.
AA> 6. Окно.
AA> 7. IDFT
AA> По идее, зануление Phase должно сломать АЧХ при п.4 или я не прав?
AA> Вопрос как занулить фазу?
П.6 что-то мне напоминает получение минимальной фазы, как это определяют
применительно к кепстру. Там очень похожий набор множителей. Hо можно ссылку
на первоисточник, я что-то не все уловил. Возможно, по полному тексту пойиу...
Евгений Машеров АКА СанитарЖеня
--- ifmail v.2.15dev5
* Origin: FidoNet Online - http://www.fido-online.com (2:5020/175.2)
Вернуться к списку тем, сортированных по: возрастание даты уменьшение даты тема автор
|
