Фильтр Чебышёва
Шаблон:Линейные электронные фильтры Фильтр Чебышёва[К 1] — один из типов линейных аналоговых или цифровых фильтров, отличительной особенностью которого является более крутой спад амплитудно-частотной характеристики (АЧХ) и существенные пульсации амплитудно-частотной характеристики на частотах полос пропускания (фильтр Чебышёва I рода) и подавления (фильтр Чебышёва II рода), чем у фильтров других типов. Фильтр получил название в честь известного русского математика XIX века Пафнутия Львовича Чебышёва, так как характеристики этого фильтра основываются на многочленах Чебышёва.
Фильтры Чебышёва обычно используются там, где требуется с помощью фильтра небольшого порядка обеспечить требуемые характеристики АЧХ, в частности, хорошее подавление частот из полосы подавления, и при этом гладкость АЧХ на частотах полос пропускания и подавления не столь важна.
Различают фильтры Чебышёва I и II родов.
Фильтр Чебышёва I рода
Это более часто встречающаяся модификация фильтров Чебышёва. Амплитудно-частотная характеристика такого фильтра [math]\displaystyle{ n }[/math]-го порядка задаётся следующим выражением:
- [math]\displaystyle{ G_n(\omega) = \left | H_n(j \omega) \right | = \frac{1}{\sqrt{1+\varepsilon^2 T_n^2\left(\frac{\omega}{\omega_0}\right)}} }[/math]
где [math]\displaystyle{ \varepsilon }[/math] — показатель пульсаций, [math]\displaystyle{ \omega_0 }[/math] — частота среза, а [math]\displaystyle{ T_n(x) }[/math] — многочлен Чебышёва [math]\displaystyle{ n }[/math]-го порядка.
В полосе пропускания такого фильтра видны пульсации, амплитуда которых определяется показателем пульсации (англ. ripple factor) [math]\displaystyle{ \varepsilon }[/math]. В полосе пропускания многочлены Чебышёва принимают значения от 0 до 1, поэтому коэффициент усиления фильтра принимает значения от максимального [math]\displaystyle{ G=1 }[/math] до минимального [math]\displaystyle{ G=1/\sqrt{1+\varepsilon^2} }[/math]. На частоте среза [math]\displaystyle{ \omega_0 }[/math] коэффициент усиления имеет значение [math]\displaystyle{ 1/\sqrt{1+\varepsilon^2} }[/math], а на частотах выше неё продолжает уменьшаться с увеличением частоты. (Примечание: обычное определение частоты среза как частоты, когда ЛАФЧХ имеет значение −3 дБ в случае фильтра Чебышёва не работает).
В случае аналогового электронного фильтра Чебышёва его порядок равен числу реактивных компонентов (например, индуктивностей), использованных при его реализации.
Пульсации в полосе пропускания часто задаются в децибелах:
Пульсации в дБ = [math]\displaystyle{ 20 \log_{10}\sqrt{1+\varepsilon^2} }[/math].
Например, пульсации амплитудой в 3 дБ соответствуют [math]\displaystyle{ \varepsilon = 1 }[/math].
Более крутой спад характеристики может быть получен если допустить пульсации не только в полосе пропускания, но и в полосе подавления, добавив в передаточную функцию фильтра нулей на мнимой оси [math]\displaystyle{ j\omega }[/math] в комплексной плоскости. Это однако приведёт к меньшему эффективному подавлению в полосе подавления. Полученный фильтр является эллиптическим фильтром, также известным как фильтр Кауэра.
Полюса и нули
![](https://cdn.xn--h1ajim.xn--p1ai/thumb.php?f=Chebyshev8_gain.png&width=300)
Для простоты примем частоту среза равной единице. Полюса [math]\displaystyle{ (\omega_{pm}) }[/math] фильтра Чебышёва являются нулями его знаменателя. Используя комплексную частоту [math]\displaystyle{ s }[/math], получим:
- [math]\displaystyle{ 1+\varepsilon^2T_n^2(-js)=0 }[/math].
Представив [math]\displaystyle{ -js=\cos(\theta) }[/math] и используя тригонометрическое определение многочленов Чебышёва, получим:
- [math]\displaystyle{ 1+\varepsilon^2T_n^2(\cos(\theta))=1+\varepsilon^2\cos^2(n\theta)=0 }[/math].
Разрешим последнее выражение относительно [math]\displaystyle{ \theta }[/math]
- [math]\displaystyle{ \theta=\frac{1}{n}\arccos\left(\frac{\pm j}{\varepsilon}\right)+\frac{m\pi}{n} }[/math].
Тогда полюса фильтра Чебышёва определяются из следующего выражения:
- [math]\displaystyle{ s_{pm}=i\cos(\theta)= }[/math]
- [math]\displaystyle{ =i\cos\left(\frac{1}{n}\arccos\left(\frac{\pm j}{\varepsilon}\right)+\frac{m\pi}{n}\right) }[/math].
Используя свойства тригонометрических и гиперболических функций, запишем последнее выражение в комплексной форме:
- [math]\displaystyle{ s_{pm}^\pm=
\pm\,\mathop{\mathrm{sh}}\left(\frac{1}{n}\mathop{\mathrm{arsh}}\left(\frac{1}{\varepsilon}\right)\right)\sin(\theta_m)+ }[/math]
- [math]\displaystyle{ +j\mathop{\mathrm{ch}}\left(\frac{1}{n}\mathop{\mathrm{arsh}}\left(\frac{1}{\varepsilon}\right)\right)\cos(\theta_m) }[/math],
где [math]\displaystyle{ m=1,\;2,\;\ldots,\;n }[/math] и
- [math]\displaystyle{ \theta_m=\frac{\pi}{2}\,\frac{2m-1}{n} }[/math].
Это выражение можно рассматривать как параметрическое уравнение с параметром [math]\displaystyle{ \theta_n }[/math]. Оно показывает, что полюса лежат на эллипсе в [math]\displaystyle{ s }[/math]-плоскости, причём центр эллипса находится в точке [math]\displaystyle{ s=0 }[/math], полуось действительной оси имеет длину [math]\displaystyle{ \mathop{\mathrm{sh}}(\mathop{\mathrm{arsh}}(1/\varepsilon)/n) }[/math], а полуось мнимой оси имеет длину [math]\displaystyle{ \mathop{\mathrm{ch}}(\mathop{\mathrm{arsh}}(1/\varepsilon)/n) }[/math].
Передаточная функция
Уравнение, выведенное выше, содержит полюса, относящиеся к комплексному коэффициенту усиления фильтра [math]\displaystyle{ G }[/math]. Для каждого полюса есть комплексно-сопряжённый, а для каждой комплексно-сопряжённой пары есть два полюса, отличающихся от них только знаком действительной части полюса. Передаточная функция должна быть устойчивой, что означает, что её полюса должны иметь отрицательную действительную часть, то есть лежать в левой полуплоскости комплексной плоскости. Передаточная функция в этом случае задаётся следующим выражением:
- [math]\displaystyle{ H(s)=\prod_{m=0}^{n-1}\frac{1}{(s-s_{pm}^-)} }[/math]
где [math]\displaystyle{ s_{pm}^- }[/math] — только те полюса, которые имеют отрицательную действительную часть.
Групповая задержка
Групповая задержка определяется как минус производная фазы фильтра по частоте и является мерой искажения фазы сигнала на различных частотах.
- [math]\displaystyle{ \tau_g=-\frac{d}{d\omega}\arg(H(j\omega)) }[/math]
Фазовые характеристики
![](https://cdn.xn--h1ajim.xn--p1ai/thumb.php?f=Chebyshev10_Phase_rus.png&width=250)
Фазовые характеристики фильтра Чебышёва I рода — фазо-частотная характеристика (ФЧХ) и фазовая задержка — представлены на рисунке. Фазо-частотная характеристика показывает распределение по частоте смещения фазы выходного сигнала относительно входного. Фазовая задержка определяется как частное от деления фазо-частотной характеристики на частоту и характеризует распределение по частоте временного смещения выходного сигнала относительно входного.
- [math]\displaystyle{ \tau_{\varphi}=\frac{\arg H(j\omega)}{\omega} }[/math]
Временны́е характеристики
![](https://cdn.xn--h1ajim.xn--p1ai/thumb.php?f=Chebyshev10_Time_Domain_rus.png&width=250)
Временные характеристики фильтра Чебышёва I рода — импульсная переходная функция и переходная функция — представлены на рисунке. Импульсная переходная функция представляет собой реакцию фильтра на входной сигнал в виде дельта-функции Дирака, а переходная функция — реакцию на входное воздействие в виде единичной функции Хевисайда.
Фильтр Чебышёва II рода
Фильтр Чебышёва II рода (инверсный фильтр Чебышёва) используется реже, чем фильтр Чебышёва I рода ввиду менее крутого спада амплитудной характеристики, что приводит к увеличению числа компонентов. У него отсутствуют пульсации в полосе пропускания, однако присутствуют в полосе подавления. Амплитудная характеристика такого фильтра задаётся следующим выражением:
- [math]\displaystyle{ G_n(\omega,\;\omega_0) = \frac{1}{\sqrt{1+ \frac{1} {\varepsilon^2 T_n ^2 \left ( \omega_0 / \omega \right )}}} }[/math]
В полосе подавления полиномы Чебышёва принимают значения от 0 до 1, из-за чего амплитудная характеристика такого фильтра принимает значения от нуля до
- [math]\displaystyle{ \frac{1}{\sqrt{1+ \frac{1}{\varepsilon^2}}} }[/math]
минимальной частотой, при которой достигается этот максимум является частота среза [math]\displaystyle{ \omega_0 }[/math]. Параметр [math]\displaystyle{ \varepsilon }[/math] связан с затуханием в полосе подавления [math]\displaystyle{ \gamma }[/math] в децибелах следующим выражением:
- [math]\displaystyle{ \varepsilon = \frac{1}{\sqrt{10^{0,\!1\gamma}-1}} }[/math]
Для затухания на частотах полосы подавления в 5 дБ: [math]\displaystyle{ \varepsilon=0,\!6801 }[/math]; для затухания в 10 дБ: [math]\displaystyle{ \varepsilon=0,\!3333 }[/math]. Частота [math]\displaystyle{ f_C=\omega_C/(2\pi) }[/math] является частотой среза. Частота затухания в 3 дБ [math]\displaystyle{ f_H }[/math] связана с [math]\displaystyle{ f_C }[/math] следующим выражением:
- [math]\displaystyle{ f_H = f_C\,\mathop{\mathrm{ch}}\left(\frac{1}{n}\mathop{\mathrm{ch}}^{-1}\frac{1}{\varepsilon}\right) }[/math].
Полюса и нули
![](https://cdn.xn--h1ajim.xn--p1ai/thumb.php?f=ChebyshevII8_gain.png&width=300)
Приняв частоту среза равной единице, получим выражение для полюсов [math]\displaystyle{ (\omega_{pm}) }[/math] фильтра Чебышёва:
- [math]\displaystyle{ 1+\varepsilon^2T_n^2(-1/js_{pm})=0 }[/math].
Полюса фильтра Чебышёва II рода представляют собой «инверсию» полюсов фильтра Чебышёва I рода:
- [math]\displaystyle{ \frac{1}{s_{pm}^\pm}=
\pm\,\mathop{\mathrm{sh}}\left(\frac{1}{n}\mathop{\mathrm{arsh}}\left(\frac{1}{\varepsilon}\right)\right)\sin(\theta_m)+ }[/math]
- [math]\displaystyle{ +j\mathop{\mathrm{ch}}\left(\frac{1}{n}\mathop{\mathrm{arsh}}\left(\frac{1}{\varepsilon}\right)\right)\cos(\theta_m) }[/math],
где [math]\displaystyle{ m=1,\;2,\;\ldots,\;n }[/math].
Нули [math]\displaystyle{ (\omega_{zm}) }[/math] фильтра Чебышёва II рода определяются из следующего соотношения::
- [math]\displaystyle{ \varepsilon^2T_n^2(-1/js_{zm})=0 }[/math].
Нули фильтра Чебышёва II рода являются «инверсией» нулей многочленов Чебышёва:
- [math]\displaystyle{ 1/s_{zm} = \cos\left(\frac{\pi}{2}\,\frac{2m-1}{n}\right) }[/math],
где [math]\displaystyle{ m=1,\;2,\;\ldots,\;n }[/math].
Передаточная функция
Передаточная функция задаётся при помощи полюсов в левой полуплоскости комплексной плоскости, её нули совпадают с нулями модуля амплитудной характеристики, с тем лишь отличием, что их порядок равен 1.
Групповая задержка
Амплитудная характеристика и групповая задержка показаны на графике. Можно видеть, что пульсации амплитуды приходятся на полосу подавления, а не на полосу пропускания.
Фазовые характеристики
![](https://cdn.xn--h1ajim.xn--p1ai/thumb.php?f=ChebyshevII10_Phase_rus.png&width=250)
Фазовые характеристики фильтра Чебышёва II рода — фазо-частотная характеристика и фазовая задержка — представлены на рисунке. Фазо-частотная характеристика показывает распределение по частоте смещения фазы выходного сигнала относительно входного. Фазовая задержка определяется как частное от деления фазо-частотной характеристики на частоту и характеризует распределение по частоте временного смещения выходного сигнала относительно входного.
Временные характеристики
![](https://cdn.xn--h1ajim.xn--p1ai/thumb.php?f=ChebyshevII10_Time_Domain_rus.png&width=250)
Временные характеристики фильтра Чебышёва II рода — импульсная переходная функция и переходная функция — представлены на рисунке. Импульсная переходная функция представляет собой реакцию фильтра на входной сигнал в виде дельта-функции Дирака, а переходная функция — реакцию на входное воздействие в виде единичной функции Хевисайда.
Цифровые фильтры Чебышёва
Фильтры Чебышёва часто реализуются в цифровой форме. Для того, чтобы от аналогового фильтра перейти к цифровому, необходимо над каждым каскадом фильтра осуществить билинейное преобразование. Весь фильтр получается путём последовательного соединения каскадов. Простой пример фильтра Чебышёва низких частот I рода чётного порядка[источник не указан 2804 дня]:
Z-преобразование каждого каскада:
- [math]\displaystyle{ S(Z) =\frac{a(Z)}{b(Z)}=\frac{\alpha_0 + \alpha_1 \cdot Z^{-1}+ \alpha_2 \cdot Z^{-2}}{1 + \beta_1 \cdot Z^{-1} + \beta_2 \cdot Z^{-2}} }[/math].
Во временно́й области преобразование записывается как:
- [math]\displaystyle{ y[n]=\alpha_0 \cdot x[0] + \alpha_1 \cdot x[-1] + \alpha_2 \cdot x[-2] - \beta_1 \cdot y[-1] - \beta_2 \cdot y[-2] }[/math]
Коэффициенты [math]\displaystyle{ \alpha_i }[/math] и [math]\displaystyle{ \beta_i }[/math] подсчитываются из коэффициентов [math]\displaystyle{ a_i }[/math] и [math]\displaystyle{ b_i }[/math]:
- [math]\displaystyle{ K = \mathop{\mathrm{tg}}\left( \pi \frac{\mbox{Frequency}}{\mbox{SampleRate}}\right) }[/math][источник не указан 2804 дня]
- [math]\displaystyle{ \mbox{temp}_i =\cos\frac{(2i+1)\pi}{n} }[/math][источник не указан 2804 дня]
- [math]\displaystyle{ b_i = \frac{1}{\mathop{\mathrm{ch}}^2\gamma-\mbox{temp}_i ^2} }[/math][источник не указан 2804 дня]
- [math]\displaystyle{ a_i = K \cdot b_i \cdot \mathop{\mathrm{sh}}\,\gamma \cdot 2\,\mbox{temp}_i }[/math][источник не указан 2804 дня]
- [math]\displaystyle{ \alpha_0 = K \cdot K }[/math]
- [math]\displaystyle{ \alpha_1 = 2 \cdot K^2 }[/math]
- [math]\displaystyle{ \alpha_2 = K \cdot K }[/math]
- [math]\displaystyle{ \beta_0^\prime = (a_i + K^2 + b_i) }[/math]
- [math]\displaystyle{ \beta_1^\prime = 2 \cdot (b_i - K^2) }[/math]
- [math]\displaystyle{ \beta_2^\prime = (a_i - K^2 - b_i) }[/math]
- [math]\displaystyle{ \beta_1 = \beta_1^\prime / \beta_0^\prime }[/math]
- [math]\displaystyle{ \beta_2 = \beta_2^\prime / \beta_0^\prime }[/math]
Для получения фильтра Чебышёва более высокого порядка, необходимо соединить последовательно несколько каскадов.
Сравнение с другими линейными фильтрами
Ниже представлены графики АЧХ фильтра Чебышёва I и II родов в сравнении с некоторыми другими фильтрами с тем же числом коэффициентов:
По графикам видно, что амплитудная характеристики фильтров Чебышёва имеет более крутой спад, чем у фильтров Баттерворта, но не такой крутой, как у эллиптического фильтра.
См. также
- Гребенчатый фильтр
- Фильтр (электроника)
- Фильтр с бесконечной импульсной характеристикой
- Цифровая обработка изображений
- Цифровая обработка сигналов
Комментарии
- ↑ Вопреки распространённому произношению старинной дворянской фамилии учёного — Чебышёв[1][2][3] — с ударением на первый слог (Чéбышев), обусловленному характерной для XX века тенденцией к обособлению фамилий на -ов/-ёв от исходных притяжательных прилагательных[2] и традиционным неразличением е/ё на письме, 4-е издание академического «Русского орфографического словаря» (2013), словарь ударений «Собственные имена в русском языке» (2001) и профильные академические издания, последовательно использующие ё при передаче имён и названий, фиксируют в качестве орфографической и орфоэпической нормы написание и произношение Чебышёв[4][5][6][7].
Примечания
- ↑ Чебышев Пафнутий Львович / Б. В. Гнеденко // Чаган — Экс-ле-Бен. — М. : Советская энциклопедия, 1978. — (Большая советская энциклопедия : [в 30 т.] / гл. ред. А. М. Прохоров ; 1969—1978, т. 29).. — В заголовке статьи: «Чебышев (произносится Чебышёв) Пафнутий Львович…»
- ↑ 2,0 2,1 Унбегаун, Б. О. Русские фамилии / пер. с англ. Л. В. Куркиной, В. П. Нерознака, Е. Р. Сквайрс; ред. Н. Н. Попов. — М. : Прогресс, 1989. — С. 349. — ISBN 5-01-001045-3.
- ↑ Калиткин, Н. Н. Численные методы : учебное пособие. — 2-е изд., испр. — СПб. : БХВ-Петербург, 2011. — С. 33 [чебышёвская система функций], 465 [чебышёвский набор шагов], 552 [критерий Чебышёва], 574 [многочлены Чебышёва]. — (Учебная литература для вузов). — ISBN 978-5-9775-0500-0.
- ↑ Чебышёв [многочлены Чебышёва, формула Чебышёва] ; чебышёвский // Русский орфографический словарь / Российская академия наук. Институт русского языка им. В. В. Виноградова; под ред. В. В. Лопатина, О. Е. Ивановой. — Изд. 4-е, испр. и доп. — М. : АСТ-ПРЕСС КНИГА, 2013. — С. 819. — (Фундаментальные словари русского языка). — ISBN 978-5-462-01272-3.
- ↑ Агеенко, Ф. Л. Чебышёв Пафнýтий // Собственные имена в русском языке : словарь ударений. — М. : Издательство НЦ ЭНАС, 2001. — С. 349. — ISBN 5-93196-107-0.
- ↑ Журнал вычислительной математики и математической физики. — М. : Издательство АН СССР, 1982. — Т. 22, № 1. — С. 142 [чебышёвский центр множества].
- ↑ Математический сборник. — М. : Наука, 2004. — Т. 195. — С. 29 [чебышёвский альтернанс], 56—57 [чебышёвский метод].
Библиография
- Кривицкий, Б. Х. Справочник по теоретическим основам радиоэлектроники. — М. : Энергия, 1977.
- Лукас, В. А. Теория автоматического управления. — M. : Недра, 1990.
- Daniels, Richard W. Approximation Methods for Electronic Filter Design. — New York : McGraw-Hill, 1974. — ISBN 0-07-015308-6.
- Higgins, Richard J. Digital Signal Processing in VLSI. — Paramus, NJ : Prentice-Hall, 1990. — ISBN 0-13-212887-X.
- Haykin, S. Adaptive Filter Theory. — 4th ed. — Paramus, NJ : Prentice-Hall, 2001. — ISBN 0-13-090126-1.
- Honig, Michael L.; Messerschmitt, David G. Adaptive Filters – Structures, Algorithms, and Applications. — Hingham, MA : Kluwer Academic Publishers, 1984. — ISBN 0-89838-163-0.
- Markel, J. D.; Gray Jr., A. H. Linear Prediction of Speech. — New York : Springer-Verlag, 1982. — ISBN 0-387-07563-1.
- Oppenheim, A. V.; Schafer, R. W. Digital Signal Processing. — Paramus, NJ : Prentice-Hall, 1975. — ISBN 0-13-214635-5.
- Proakis, John G.; Manolakis, Dimitris G. Introduction to Digital Signal Processing. — Paramus, NJ : Prentice-Hall, 1988. — ISBN 0-02-396815-X.
- Rabiner, L. R.; Schafer, R. W. Digital Processing of Speech Signals. — Paramus, NJ : Prentice-Hall, 1978. — ISBN 0-13-213603-1.
- Rabiner, L. R.; Gold, B. Theory and Application of Digital Signal Processing. — Paramus, NJ : Prentice-Hall, 1986. — ISBN 0-13-914101-4.
- Rorabaugh, Britton C. Approximation Methods for Electronic Filter Design. — New York : McGraw-Hill, 1999. — ISBN 0-07-054004-7.
- Smith, Steven W. The Scientist and Engineer’s Guide to Digital Signal Processing. — 2nd ed. — San-Diego : California Technical Publishing, 1999. — ISBN 0-9660176-4-1.
- Widrow, B.; Stearns, S. D. Adaptive Signal Processing. — Paramus, NJ : Prentice-Hall, 1985. — ISBN 0-13-004029-0.
Ссылки
- Как рассчитывать рекурсивные фильтры? (недоступная ссылка). НГУ. Дата обращения: 14 февраля 2014. Архивировано 30 ноября 2013 года.
- Классификация фильтров . Аналоговые измерительные устройства. Дата обращения: 14 февраля 2014.
- Лекция по цифровой фильтрации . DSP.sut.ru. Дата обращения: 14 февраля 2014. Архивировано 12 марта 2007 года.
- Расчёт фильтра Чебышёва первого рода с примерами . Dsplib.ru. Дата обращения: 14 февраля 2014.
- Расчёт фильтра Чебышёва второго рода с примерами . Dsplib.ru. Дата обращения: 14 февраля 2014.
- Фильтры нижних частот . Gaw.ru (Рынок Микроэлектроники). Дата обращения: 14 февраля 2014.
- Chebyshev Filters . Texas Instruments. Дата обращения: 14 февраля 2014. Архивировано 5 августа 2012 года. (англ.)
В статье есть список источников, но не хватает сносок. |