Режекторные фильтры

Режекторный фильтр постоянной составляющей сигнала. Сконст- руируем простейший РЦФ, добавив к оператору (6.2.1) один полюс вне единичной окружности на малом расстоянии от нуля:

H_п(z) = G(1-z)/(1-az), z_p= 1/a. (6.2.2)

Допустим, что полюс помещен в точке z_p1= 1.01, при этом, а=0,99. Масштабный коэффициент G получим нормировкой H(z) к 1 на частоте Найквиста. Для приведенных условий G=0.995. Отсюда, при Dt=1:

H_п(z) = 0,995(1-z)/(1-0.99z),

y_k = 0.995(x_k-x_k-1)+ 0.99y_k-1.

Рис. 6.2.3.

Отображение нуля n₁ и полюса р₁ на z-плоскости и АЧХ фильтра для исключения постоянной составляющей приведены на рис.6.2.1. Коэффициент передачи сигнала на произвольной частоте w_i равен отношению длин векторов V_n1(z) и V_p1(z) соответственно из нуля и полюса до точки z(w_i) на единичной окружности и близок к единице для всех частот, за исключением нулевой:

|H_п(z)| = G V_n1(z)/V_p1(z).

Фазочастотная характеристика фильтра приведена на рис. 6.2.3 и определяется разностью фазовых углов векторов V_n1(z) и V_p1(z):

j_п(w) = j_n1-j_p1.

Режекторный фильтр произвольной частоты. При проектировании на подавление любой другой частоты ω_v нули и полюсы располагаются на соответствующем радиусе z-плоскости. Радиальный угол направления на нуль и полюс определяются выражением:

j_v = p·w_v/w_N. (6.2.3)

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

H_v(z) = G(z-z_n)(z-z_n*)/[(z-z_p)(z-z_p*)]. (6.2.4)

Нули фильтра располагаются на единичной окружности:

z_n = cos j_v + j sin j_v = Re z_n + j Im z_n. (6.2.5)

Полюсы - на полярном радиусе R:

z_p = R·cos j_v + j R·sin j_v = Re z_p + j Im z_p. (6.2.6)

Пример положения нулей (n2 и n2*) и полюсов (р2 и р2*) приведен на рис.6.2.1. Подставляя (6.2.5-6.2.6) в (6.2.4), получаем:

H_v(z) =, (6.2.7)

G = [1+(1+2Re z_p)/R²] / (2+2Re z_n). (6.2.8)

При приведении уравнения (6.2.7) в типовую форму:

H_v(z) =, (6.2.7')

b₀ = 1, b₁ = -2·Re z_n, b₂ = 1. (6.2.9)

a₁ = - (2·Re z_p)/R², a₂ = 1/R².

Соответственно, алгоритм вычислений:

y_k = G·(x_k+b₁·x_k-1+x_k-2) – a₁·y_k-1 – a₂·y_k-2. (6.2.10)

В качестве примера проведем расчет режекторного фильтра на сетевую частоту питания приборов f_s = 50 Гц, которая очень часто попадает в измеренные данные. При шаге дискретизации данных Dt = 0.001 сек радиальный угол на нули и полюса фильтра в z-плоскости:

f_N = 1/2Dt = 500 Гц, j = p·f_s/f_N = 0.1π.

Радиус полюса фильтра примем равным R = 1.01. Значения нуля и полюса:

z_n = cos j + j sin j = 0.951 + 0.309 j,

z_p = R·cos j_v + j R·sin j_v = 0.961 + 0.312 j.

Рис. 6.2.4.

Значение масштабного множителя G по (6.2.8):

G = 0.99.

Значения коэффициентов передаточной функции:

b₁ = -2·Re z_n = -1.902,

a₁ = - (2·Re z_p)/R² = -1.883, a₂ = 1/R² = 0.98.

При подстановке коэффициентов в уравнение (6.2.7') и замене z = exp(-jω) может быть получена частотная передаточная функция фильтра, которая приведена на рис. 6.2.4:

Алгоритм фильтра:

y_k = 0.99·(x_k - 1.902·x_k-1 + x_k-2) + 1.883·y_k-1 – 0.98·y_k-2.

На рис. 6.2.5 приведен модельный входной сигнал фильтра, состоящий из суммы двух равных по амплитуде гармоник с частотой 50 и 53 Гц, и сигнал на выходе фильтра (смещен вверх). Справа на рисунке приведены спектры входного и выходного сигналов. Спектр выходного сигнала зарегистрирован после интервала установления реакции фильтра, который хорошо заметен на начальной части графика выходного сигнала. После установления сигнал на выходе фильтра практически полностью освобожден от гармоники 50 Гц.

Рис. 6.2.5.

Рис. 6.2.6.

При R → 1 ширина полосы подавления фильтра становится все более узкой, но при этом увеличивается длительность импульсной реакции фильтра и, соответственно, время установления фильтра при изменении спектра входного сигнала. В первом приближении значимая часть импульсной реакции режекторных фильтров равна (4÷5)/(R-1). Пример импульсной реакции для фильтра, вычисленного выше, приведен на рис. 6.2.6. Отклик фильтра получен при подаче на вход РЦФ импульса Кронекера. Для наглядности реакции на графике не показан начальный пик отклика (отсчет на нулевой точке), амплитуда которого равна значению G.