Фильтр Баттерворта

Материал из Википедии — свободной энциклопедии

Линейные электронные фильтры
Фильтр Баттерворта
Фильтр Чебышёва
Эллиптический фильтр
Фильтр Бесселя
Фильтр Гаусса
Фильтр Лежандра
Править

Фильтр Баттерво́рта — один из типов электронных фильтров. Фильтры этого класса отделяются от других методом проектирования. Фильтр Баттерворта проектируется так, чтобы его амплитудная частотная характеристика была максимально гладкой на частотах полосы пропускания.

Подобные фильтры были впервые описаны британским инженером Стефаном Баттервортом в статье «О теории фильтрующих усилителей» (англ. On the Theory of Filter Amplifiers), в журнале Wireless Engineer в 1930 году.

Содержание

1 Обзор
2 Проектирование фильтра
- 2.1 Топология Кауэра
- 2.2 Топология Саллена-Кея
3 Сравнение с другими линейными фильтрами
4 Пример
5 См. также
6 Библиография
7 Ссылки

[править] Обзор

АЧХ фильтра Баттерворта максимально гладкая на частотах полосы пропускания и снижается практически до нуля на частотах полосы подавления. При отображении частотного отклика фильтра Баттерворта на логарифмической АФЧХ, амплитуда снижается к минус бесконечности на частотах полосы подавления. В случае фильтра первого порядка АЧХ затухает со скоростью −6 децибел на октаву (-20 децибел на декаду) (на самом деле все фильтры первого порядка независимо от типа идентичны и имеют одинаковый частотный отклик). Для фильтра Баттерворта второго порядка АЧХ затухает на −12 дБ на октаву, для фильтра третьего порядка — на −18 дБ и так далее. АЧХ фильтра Баттерворта — монотонно убывающая функция частоты. Фильтр Баттерворта — единственный из фильтров, сохраняющий форму АЧХ для более высоких порядков (за исключением более крутого спада характеристики на полосе подавления) тогда как многие другие разновидности фильтров (фильтр Бесселя, фильтр Чебышева, эллиптический фильтр) имеют различные формы АЧХ при различных порядках.

В сравнении с фильтрами Чебышева I и II типов или эллиптическим фильтром, фильтр Баттерворта имеет более пологий спад характеристики и поэтому должен иметь больший порядок (что более трудно в реализации) для того, чтобы обеспечить нужные характеристики на частотах полосы подавления. Однако фильтр Баттерворта имеет более линейную фазо-частотную характеристику на частотах полосы пропускания.

АЧХ для фильтров Баттерворта низких частот порядка от 1 до 5. Наклон характерстики — 20n дБ/декаду, где n — порядок фильтра.

Как и для всех фильтров при рассмотрении частотных характеристик используют фильтр низких частот, из которого легко можно получить фильтр высоких частот, а, включив несколько таких фильтров последовательно, — полосовой фильтр или режекторный фильтр.

Амплитудно-частотная характеристика $G(\omega) \!$ фильтра Баттерворта $n \!$ -го порядка может быть получена из передаточной функции $H(s) \!$ :

$G^2(\omega)=\left |H(j\omega)\right|^2 = \frac {G_0^2}{1+\left(\frac{\omega}{\omega_c}\right)^{2n}}$

где

$n \!$ — порядок фильтра
$\omega_c \!$ — частота среза (частота на которой амплитуда равна примерно −3dB)
$G_0 \!$ — коэффициент усиления по постоянной составляющей (усиление на нулевой частоте)

Легко заметить, что для бесконечных значений $n \!$ АЧХ становится прямоугольной функцией, и частоты ниже частоты среза будут пропускаться с коэффициентом усиления $G_0 \!$ , а частоты выше частоты среза будут полностью подавляться. Для конечных значений $n \!$ спад характеристики будет пологим.

С помощью формальной замены $s=+j\omega \!$ представим выражение $H(s)H(-s) \!$ в виде $|H(\omega)|^2 \!$ :

$H(s)H(-s) = \frac {G_0^2}{1+\left (\frac{-s^2}{\omega_c^2}\right)^n}$

Полюсы передаточной функции расположены на круге радиуса $\omega_c \!$ равноудалённо друг от друга в левой полуплоскости. То есть передаточную функцию фильтра Баттерворта можно определить лишь определением полюсов его передаточной функции в левой полуплоскости s-плоскости. $k \!$ -й полюс определяется из следующего выражения:

$-\frac{s_k^2}{\omega_c^2} = (-1)^{\frac{1}{n}} = e^{\frac{j(2k-1)\pi}{n}} \qquad\mathrm{k = 1,2,3, \ldots, n}$

откуда

$s_k = \omega_c e^{\frac{j(2k+n-1)\pi}{2n}}\qquad\mathrm{k = 1,2,3, \ldots, n}$

Передаточную функцию можно записать в виде:

$H(s)=\frac{G_0}{\prod_{k=1}^n (s-s_k)/\omega_c}$

Аналогичные рассуждения применимы и к цифровым фильтрам Баттерворта, с той лишь разницей, что соотношения записываются не для s-плоскости, а для z-плоскости.

Знаменатель этой передаточной функции называется полиномом Баттерворта.

[править] Нормированные полиномы Баттерворта

Полиномы Баттерворта могут записываться в комплексной форме, как показано выше, однако обычно они записываются в виде соотношений с вещественными коэффициентами (комплексно-сопряжённые пары объединяются с помощью умножения). Нормируются полиномы по частоте среза: $\omega_c=1 \!$ . Нормированные полиномы Баттерворта, таким образом, имеют следующую каноническую форму:

$B_n(s)=\prod_{k=1}^{\frac{n}{2}} \left[s^2-2s\cos\left(\frac{2k+n-1}{2n}\,\pi\right)+1\right]$ , $n \!$ — чётно

$B_n(s)=(s+1)\prod_{k=1}^{\frac{n-1}{2}} \left[s^2-2s\cos\left(\frac{2k+n-1}{2n}\,\pi\right)+1\right]$ , $n \!$ — нечётно

Ниже представлены коэффициенты полиномов Баттерворта для первых восьми порядков:

$n \!$	Коэффициенты полиномов $B_n(s) \!$
1	$(s+1) \!$
2	$\! s^2+1.4142s+1$
3	$\! (s+1)(s^2+s+1)$
4	$\! (s^2+0.7654s+1)(s^2+1.8478s+1)$
5	$\! (s+1)(s^2+0.6180s+1)(s^2+1.6180s+1)$
6	$\! (s^2+0.5176s+1)(s^2+1.4142s+1)(s^2+1.9319s+1)$
7	$\! (s+1)(s^2+0.4450s+1)(s^2+1.2470s+1)(s^2+1.8019s+1)$
8	$\! (s^2+0.3902s+1)(s^2+1.1111s+1)(s^2+1.6629s+1)(s^2+1.9616s+1)$

[править] Максимальная гладкость

Приняв $\! \omega_c=1$ и $\! G_0=1$ , производная амплитудной характеристики по частоте будет выглядеть следующим образом:

$\frac{dG}{d\omega}=-nG^3\omega^{2n-1} \!$

Она монотонно убывает для всех $\omega \!$ так как коэффициент усиления всегда положителен. Таким образом, АЧХ фильтра Баттерворта не имеет пульсаций. При разложении амплитудной характеристи в ряд, получим:

$G(\omega)=1 - \frac{1}{2}\omega^{2n}+\frac{3}{8}\omega^{4n}+\ldots$

Другими словами, все производные амлитудно-частотной характерситики по частоте до 2n-й равны нулю, из чего следует «максимальная гладкость».

[править] Спад характеристики на высоких частотах

Приняв $\omega_c=1\!$ , найдём наклон логарифма АЧХ на высоких частотах:

$\lim_{\omega\rightarrow\infty}\frac{d\log(G)}{d\log(\omega)}=-n$

В децибелах высокочастотная асимптота имеет наклон −20n дБ/декаду.

[править] Проектирование фильтра

Существует ряд различных топологий фильтра, с помощью которых реализуются линейные аналоговые фильтры. Эти схемы отличаются только значениями элементов, структура же остаётся неизменной.

[править] Топология Кауэра

Топология Кауэра использует пассивные элементы (ёмкости и индуктивности). Фильтр Баттеворта с заданной передаточной функцией может быть построен с в форме Кауэра 1 типа. k-й элемент фильтра задаётся соотношением:

$C_k = 2 \sin \left [\frac {(2k-1)}{2n} \pi \right ]$ ; k чётно

$L_k = 2 \sin \left [\frac {(2k-1)}{2n} \pi \right ]$ ; k нечётно

[править] Топология Саллена-Кея

Топология Саллена-Кея использует помимо пассивных также и активные элементы (операционные усилители и ёмкости). Каждый каскад схемы Саллена-Кея представляет собой часть фильтра, математически описываемую парой комплексно-сопряжённых полюсов. Весь фильтр получается последовательным соединением всех каскадов. В случае, если попадается действительный полюс, он должен быть реализован отдельно, обычно в виде RC-цепочки, и включён в общую схему.

Передаточная функция каждого каскада в схеме Саллена-Кея имеет вид:

$H(s)=\frac{1}{1+C_2(R_1+R_2)s+C_1C_2R_1R_2s^2}$

Нужно, чтобы знаменатель представлял собой один из множителей полинома Баттерворта. Приняв $\omega_c=1 \!$ , получим:

$\! C_1C_2R_1R_2=1\,$

$C_2(R_1+R_2)=2\cos\left(\frac{2k+n-1}{2n}\right)$

Последнее соотношение даёт две неизвестных, которые могут быть выбраны произвольно.

[править] Сравнение с другими линейными фильтрами

Рисунок ниже показывает АЧХ фильтра Баттерворта в сравнении с другими популярными линейными фильтрами одинакового (пятого) порядка:

Из рисунка видно, что спад АЧХ фильтра Баттерворта самый медленный из четырёх, однако он имеет и самую гладкую АЧХ на частотах полосы пропускания.

[править] Пример

Аналоговый фильтр Баттерворта низких частот (топология Кауэра) с частотой среза со следующими номиналами элементов: фарад, ом, и генри.

Аналоговый фильтр Баттерворта низких частот (топология Кауэра) с частотой среза $\omega_c = 1 \!$ со следующими номиналами элементов: $C_2=4/3 \!$ фарад, $R_4=1 \!$ ом, $L_1=1/2 \!$ и $L_3=3/2 \!$ генри.

Логарифмический график плотности передаточной функции H(s) на плоскости комплексного аргумента для фильтра Баттерворта третьего порядка с частотой среза $\omega_c = 1 \!$ . Три полюса лежат на круге единичного радиуса в левой полуплоскости.

Рассмотрим аналоговый низкочастотный фильтр Баттерворта третьего порядка с $C_2=4/3 \!$ фарад, $R_4=1 \!$ ом, $L_1=1/2 \!$ и $L_3=3/2 \!$ генри. Обозначив полное сопротивление ёмкостей C как 1/Cs и полное сопротивление индуктивностей L как Ls, где $s=\sigma+j\omega \!$ — комплексная переменная, и используя уравнения для расчёта электрических схем, получим следующую передаточную функцию для такого фильтра:

$H(s)=\frac{V_o(s)}{V_i(s)}=\frac{1}{1+2s+2s^2+s^3}$

АЧХ $G(\omega) \!$ задаётся уравнением:

$G^2(\omega)=|H(j\omega)|^2=\frac{1}{1+\omega^6}\,$

а ФЧХ задаётся уравнением:

$\Phi(\omega)=\arg(H(j\omega))\,$

Групповая задержка определяется как производная фазы по круговой частоте и является мерой искажений сигнала по фазе на различных частотах. Логарифмическая АЧХ $\log_{10}(G) \!$ такого фильтра не имеет пульсаций ни в полосе пропускания, ни в полосе подавления.

График модуля передаточной функции на комплексной плоскости ясно указывает на три полюса в левой полуплоскости. Передаточная функция полностью определяется расположением этих полюсов на единичном круге симметрично относительно действительной оси.

Заменив каждую индуктивность ёмкостью, а ёмкости — индуктивностями, получим высокочастотный фильтр Баттерворта.

$\log_{10}(G) \!$ и групповая задержка фильтра Баттерворта третьего порядка с частотй среза $\omega_c = 1 \!$

[править] См. также

[править] Библиография

Лукас В. А. Теория автоматического управления. — М.: Недра, 1990. — 416 с.: ил.
Справочник по теоретическим основам радиоэлектроники. Ред. Кривицкий Б. Х. Том 2. М.,: Энергия, 1977, с. 305—308.
Steven W. Smith, The Scientist and Engineer’s Guide to Digital Signal Processing, Second Edition, 1999, California Technical Publishing
C. Britton Rorabaugh, DSP Primer, McGraw-Hill, 1999.
Richard J. Higgins, Digital Signal Processing in VLSI, Prentice-Hall, 1990.
A. V. Oppenheim and R. W. Schafer, Digital Signal Processing, Prentice-Hall, 1975.
L. R. Rabiner and B. Gold, Theory and Application of Digital Signal Processing, Prentice-Hall, 1975.
John G. Proakis and Dimitris G. Manolakis, Introduction to Digital Signal Processing, MacMillian, 1988.
Fredrick J. Harris, On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform, Proc. IEEE, Vol. 66, No. 1, 1978 pp. 51-83.
B. Widrow and S.D. Stearns, Adaptive Signal Processing, Prentice-Hall, 1985.
S. Haykin, Adaptive Filter Theory, 3rd Edition, Prentice-Hall, 1996.
Michael L. Honig and David G. Messerschmitt, Adaptive Filters — Structures, Algorithms, and Applications, Kluwer Academic Publishers, Hingham, MA 1984.
J.D. Markel and A.H. Gray, Jr., Linear Prediction of Speech, Springer-Verlag, New York, NY, 1976.
L.R. Rabiner and R.W. Schafer, Digital Processing of Speech Signals, Prentice-Hall, 1978.

[править] Ссылки

Категория: Обработка сигналов