Пример расчета фильтра низких частот Баттеруорта

Рис. 10.1.2.

Начиная с этого параграфа, будем сопровождать рассмотрение теории последовательным расчетом фильтра низких частот с применением приводимых формул. Для расчета примем следующие исходные параметры фильтра:

- Шаг дискретизации данных Dt = 0.0005 сек. Частота Найквиста f_N = 1/2Dt = 1000 Гц, ω_N = 6.283·10³ рад.

- Граничная частота пропускания: f_p = 300 Гц, w_p = 1.885·10³ рад.

- Граничная частота подавления: f_s = 500 Гц, w_s = 3.142·10³ рад.

- Коэффициенты неравномерности: А_р = А_s = 0.1.

Расчет дополнительных параметров:

1. Значение d по формуле (10.1.3): d= 0.484.

2. Деформированные частоты по формуле (10.1.4): w_dp = 2.038·10³ рад. w_ds = 4·10³ рад.

3. Порядок фильтра по формуле (10.1.6): N = 4.483.

Для пояснения порядка расчетов при четном и нечетном порядке фильтра, принимаем N1=4, N2=5.

4. Частота среза по формуле (10.1.7): w_dc(N1) = 2.443·10³ рад (389 Гц), w_dc(N2) = 2.356·10³ рад (375 Гц).

5. По формуле H(w)= [1/(1+w^2N)]^1/2, w = ω/ω_dc, строим графики передаточных функций (рис.10.1.2).

Преобразование Лапласа. Переводим функцию |H(W)|² на координатную ось пространства преобразования Лапласа при p = jW, для чего достаточно подставить W = p/j:

|H(р)|² = 1/[1+(p/j)^2N]. (10.1.8)

Полюсы функции находятся в точках нулевых значений знаменателя:

1+(p/j)^2N = 0, p = j . (10.1.9)

Отсюда следует, что полюсы располагаются на единичной окружности в p-плоскости, а их местоположение определяется корнями уравнения (10.1.9). В полярных координатах:

p_n = j exp(jp(2n-1)/2N), n = 1,2,...,2N. (10.1.10)

p_n = j cos[p(2n-1)/2N] - sin[p(2k-1)/2N]. (10.1.10')

Рис. 10.1.2.