[an error occurred while processing this directive]
|
А что нужно-то?
Кратко. Импульсная х-ка преобр. Гильберта = 1/t, т.е. атисимметричная и бесконечная,частотная = -j для положительных частот, +j для отрицательных.
Для синтеза КИХ фильтра, выполняющего преобр. Гильберта можно использовать, например, алгоритм Ремеза. Найти его можно в любой порядочной книжке по ЦОС, я бы рекомендовал Рабинера, Гоулда. В ней,
кроме всего прочего, не только алгоритм, но и его реализация, на Фортране помнится.
Можно поступить проще - взять обратное Фурье от част. характеристики преобразователя Гильберта и наложить на нее окно.
В том и другом случае приличные характреристики преобразования получаются при довольно большом порядке фильтра.
Выше речь шла о преобр. Гильберта для всей частотной оси от -pi до pi.
Часто этого не требуется, достаточно иметь преобразователь Гильберта в какой-то полосе частот. В этом случае проще синтезировать комплексный полосовой фильтр на требуемую полосу частот, причем можно как КИХ, так и БИХ.
После прохождения действительного через такой фильтр на выходе будут обе квадратурные компоненты, да еще и отфильтровнные.
Все это о фильтрации во временной области. Если задача позволяет, можно обойтись только частотной характеристикой Гильберта. Т.е. берете спектр сигнала,умножаете его на ЧХ Гильберта и получаете спектр сопряженного сигнала. Умножаете его на j, затем прибавляете спектр исходного сигнала и получаете спектр комплексного сигнала.
То же самое можно сделать проще - берете опять же спектр сигнала, обнуляете компоненты для отрицательных частот. Для сохранения масштаба величины компонент для положительных частот удваиваете.
В этих двух способах после обратного Фурье будут завороты.
E-mail: info@telesys.ru