4.2. Идеальные частотные фильтры.

Идеальным полосовым фильтром называется фильтр, имеющий единичную амплитудно-частотную характеристику в полосе от определенной нижней частоты _н до определенной верхней частоты _в, и нулевой коэффициент передачи за пределами этой полосы (для цифровых фильтров - в главном частотном диапазоне).

Импульсная реакция фильтра (коэффициенты оператора) находится преобразованием Фурье заданной передаточной функции H(). В общем случае:

h(nt) = (1/2) H() exp(jnt) d

Для получения вещественной функции импульсного отклика фильтра действительная часть передаточной функции должна быть четной, а мнимая - нечетной. Цифровые фильтры задаются в главном частотном диапазоне, границы которого (частота Найквиста _N) определяются интервалом дискретизации данных (_N = /t), подлежащих фильтрации, и соответственно определяют интервал дискретизации оператора фильтра (t = /_N). Для фильтров с нулевым фазовым сдвигом мнимая часть передаточной функции должна быть равна нулю, при этом оператор фильтра определяется косинусным преобразованием Фурье:

h(nt)= (1/) H() cos(n/_N) dn = 0,1,2,... (4.2.1)

Для идеального полосового фильтра H()=1 в полосе частот от _н до _в, и интеграл (4.2.1) вычисляется в этих пределах. Идеальные фильтры низких и высоких частот можно считать частными случаями идеальных полосовых фильтров с полосой пропускания от 0 до _в для низкочастотного и от _н до _N для высокочастотного фильтра.

При интервале дискретизации данных t, условно принимаемым за 1, главный частотный диапазон передаточных функций ограничивается значением частоты Найквиста от - до . Если на практике интервал дискретизации данных в физических единицах отличается от 1, то это сказывается только на изменении масштаба частотной шкалы передаточных функций.

Пример 1. t = 0.1 сек. f_N = 1/2t = 5 Гц. _N =/t = 10 .

Пример 2. x = 10 метров. f_N = 0.05 м^-1. _N= 0.1 .

Во всех дальнейших выражениях значение t, если это специально не оговорено, будем принимать равным 1.

При H()=A=1 в полосе пропускания _н -_в, и H()=0 за ее пределами, для идеальных симметричных полосовых НЦФ из (4.2.1) в общем виде получаем:

h(n) = (А/) [_в sinc(n_в) - _н sinc(n_н)], (4.2.2)

h_o = (_в - _н)/, h(n) = (sin n_в - sin n_н)/(n).

где sinc(n) = sin(n)/(n) - функция интегрального синуса (функция отсчетов), бесконечная по координате .

Рис. 4.2.1. Входные сигналы. Рис. 4.2.2. Спектр сигнала и границы фильтра.

На рис. 4.2.1 приведен пример сигнала однотональной балансной амплитудной модуляции (чистого – вверху, и с наложенными шумами внизу, мощность шумов равна мощности сигнала). Если информация заключена в частоте и амплитуде модулирующего сигнала, то полосовой фильтр выделения сигнала из шумов, спектр которого для одной модулирующей частоты приведен на рис. 4.2.2, в идеальном случае должен иметь плоскую частотную характеристику в границах возможных вариаций модулирующей частоты (от _н до _в).

Рис. 4.2.3. Оператор фильтра.

На рис. 4.2.3 приведен оператор полосового фильтра, вычисленный по (4.2.2) для приведенных выше условий, с ограничением по числу коэффициентов оператора до N=100. Как видно из рисунка, оператор затухает достаточно медленно и явно усечен, что должно сказаться на форме частотной характеристики фильтра. Все дальнейшие вычисления будут проводиться на продолжении данного примера.