Традиция

Фильтр Чебышёва

Фильтр Чебышёва — один из типов линейных аналоговых или цифровых фильтров, отличительной особенностью которого является более крутой спад амплитудно-частотной характеристики (АЧХ) и существенные пульсации амплитудно-частотной характеристики на частотах полос пропускания (фильтр Чебышёва I рода) и подавления (фильтр Чебышёва II рода), чем у фильтров других типов. Фильтр получил название в честь известного русского математика XIX века Пафнутия Львовича Чебышёва, так как характеристики этого фильтра основываются на многочленах Чебышёва.

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

Фильтры Чебышёва обычно используются там, где требуется с помощью фильтра небольшого порядка обеспечить требуемые характеристики АЧХ, в частности, хорошее подавление частот из полосы подавления, и при этом гладкость АЧХ на частотах полос пропускания и подавления не столь важна.

Различают фильтры Чебышёва I и II родов.

Фильтр Чебышёва I родаПравить

 
АЧХ фильтра Чебышёва I рода четвёртого порядка с ω 0 = 1 \omega_0=1 и ε = 1 \varepsilon=1

Это более часто встречающаяся модификация фильтров Чебышёва. Амплитудно-частотная характеристика такого фильтра n n -го порядка задаётся следующим выражением: G n ( ω ) = | H n ( j ω ) | = 1 1 + ε 2 T n 2 ( ω ω 0 ) G_n(\omega) = \left | H_n(j \omega) \right | = \frac{1}{\sqrt{1+\varepsilon^2 T_n^2\left(\frac{\omega}{\omega_0}\right)}}

где ε \varepsilon — показатель пульсаций, ω 0 \omega_0 — частота среза, а T n ( x ) T_n(x) — многочлен Чебышёва n n \! -го порядка.

В полосе пропускания такого фильтра видны пульсации, амплитуда которых определяется показателем пульсации (англ. ripple factor) ε \varepsilon . В полосе пропускания многочлены Чебышёва принимают значения от 0 до 1, поэтому коэффициент усиления фильтра принимает значения от максимального G = 1 \! G=1 до минимального G = 1 / 1 + ε 2 G=1/\sqrt{1+\varepsilon^2} . На частоте среза ω 0 \omega_0 коэффициент усиления имеет значение 1 / 1 + ε 2 1/\sqrt{1+\varepsilon^2} , а на частотах выше неё продолжает уменьшаться с увеличением частоты. (Примечание: обычное определение частоты среза как частоты, когда ЛАЧХ имеет значение −3 дБ в случае фильтра Чебышёва не работает).

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

Пульсации в полосе пропускания часто задаются в децибелах:

Пульсации в дБ = 20 log 10 ⁡ 1 1 + ε 2 20 \log_{10}\frac{1}{\sqrt{1+\varepsilon^2}} .

Например, пульсации амплитудой в 3 дБ соответствуют ε = 1 \varepsilon = 1 \! .

Более крутой спад характеристики может быть получен если допустить пульсации не только в полосе пропускания, но и в полосе подавления, добавив в передаточную функцию фильтра нулей на мнимой оси j ω j\omega в комплексной плоскости. Это однако приведёт к меньшему эффективному подавлению в полосе подавления. Полученный фильтр является эллиптическим фильтром, также известным как фильтр Кауэра.

Полюсы и нулиПравить

 
Логарифм модуля амплитудной характеристики фильтра Чебышёва I рода 8-го порядка на плоскости комплексной частоты ( s = σ + j ω s=\sigma+j\omega ) при ε = 0 , 1 \varepsilon=0,\!1 и ω 0 = 1 \omega_0=1 . Белые пятна — это полюсы фильтра. Они расположены на эллипсе с полуосью 0,3836… по действительной оси и 1,071… по мнимой оси. Полюсы передаточной функции фильтра расположены в левой полуплоскости. Чёрный цвет соответствует коэффициенту усиления менее 0,05, белый соответствует коэффициенту усиления более 20.

Для простоты примем частоту среза равной единице. Полюсы ( ω p m ) (\omega_{pm}) фильтра Чебышёва являются нулями его знаменателя. Используя комплексную частоту s s , получим: 1 + ε 2 T n 2 ( − j s ) = 0. 1+\varepsilon^2T_n^2(-js)=0.

Представив − j s = cos  Косинус  ( θ ) -js=\cos(\theta) и используя тригонометрическое определение многочленов Чебышёва, получим: 1 + ε 2 T n 2 ( cos  Косинус  ( θ ) ) = 1 + ε 2 cos  Косинус  2 ( n θ ) = 0. 1+\varepsilon^2T_n^2(\cos(\theta))=1+\varepsilon^2\cos^2(n\theta)=0.

Разрешим последнее выражение относительно θ \theta θ = 1 n arccos  Арккосинус  ( ± j ε ) + m π n . \theta=\frac{1}{n}\arccos\left(\frac{\pm j}{\varepsilon}\right)+\frac{m\pi}{n}.

Тогда полюсы фильтра Чебышёва определяются из следующего выражения: s p m = i cos  Косинус  ( θ ) = s_{pm}=i\cos(\theta)= = i cos  Косинус  ( 1 n arccos  Арккосинус  ( ± j ε ) + m π n ) . =i\cos\left(\frac{1}{n}\arccos\left(\frac{\pm j}{\varepsilon}\right)+\frac{m\pi}{n}\right).

Используя свойства тригонометрических и гиперболических функций, запишем последнее выражение в комплексной форме: s p m ± = ± sh ⁡ ( 1 n arsh ⁡ ( 1 ε ) ) sin  Синус  ( θ m ) + s_{pm}^\pm= \pm\,\mathop{\mathrm{sh}}\left(\frac{1}{n}\mathop{\mathrm{arsh}}\left(\frac{1}{\varepsilon}\right)\right)\sin(\theta_m)+ + j ch ⁡ ( 1 n arsh ⁡ ( 1 ε ) ) cos  Косинус  ( θ m ) , +j\mathop{\mathrm{ch}}\left(\frac{1}{n}\mathop{\mathrm{arsh}}\left(\frac{1}{\varepsilon}\right)\right)\cos(\theta_m) ,

где m = 1 , 2 , … , n m=1,\;2,\;\ldots,\;n и θ m = π 2 2 m − 1 n . \theta_m=\frac{\pi}{2}\,\frac{2m-1}{n}.

Это выражение можно рассматривать как параметрическое уравнение с параметром θ n \theta_n . Оно показывает, что полюсы лежат на эллипсе в s s -плоскости, причём центр эллипса находится в точке s = 0 s=0 , полуось действительной оси имеет длину sh ⁡ ( arsh ⁡ ( 1 / ε ) / n ) \mathop{\mathrm{sh}}(\mathop{\mathrm{arsh}}(1/\varepsilon)/n) , а полуось мнимой оси имеет длину ch ⁡ ( arsh ⁡ ( 1 / ε ) / n ) \mathop{\mathrm{ch}}(\mathop{\mathrm{arsh}}(1/\varepsilon)/n) .

Передаточная функцияПравить

Уравнение, выведенное выше, содержит полюсы, относящиеся к комплексному коэффициенту усиления фильтра G G . Для каждого полюса есть комплексно-сопряжённый, а для каждой комплексно-сопряжённой пары есть два полюса, отличающихся от них только знаком. Передаточная функция должна быть устойчивой, что означает, что её полюсы должны иметь отрицательную действительную часть, то есть лежать в левой полуплоскости комплексной плоскости. Передаточная функция в этом случае задаётся следующим выражением: H ( s ) = ∏ m = 0 n − 1 1 ( s − s p m − ) H(s)=\prod_{m=0}^{n-1}\frac{1}{(s-s_{pm}^-)}

где s p m − s_{pm}^- — только те полюсы, которые имеют отрицательную действительную часть.

Групповая задержкаПравить

 
Амплитуда и групповая задержка фильтра Чебышёва I рода пятого порядка с ε = 0 , 5 \varepsilon=0,\!5 . Видно, что в полосе пропускания и АЧХ и групповая задержка имеют пульсации, в полосе подавления этих пульсаций нет.

Групповая задержка определяется как производная фазы фильтра по частоте и является мерой искажения фазы сигнала на различных частотах. τ g = d d ω arg ⁡ ( H ( j ω ) ) \tau_g=\frac{d}{d\omega}\arg(H(j\omega))

Фазовые характеристикиПравить

 
Типовая ФЧХ и фазовая задержка фильтра Чебышёва I рода 10-го порядка.

Фазовые характеристики фильтра Чебышёва I рода — фазо-частотная характеристика (ФЧХ) и фазовая задержка — представлены на рисунке. Фазо-частотная характеристика показывает распределение по частоте смещения фазы выходного сигнала относительно входного. Фазовая задержка определяется как частное от деления фазо-частотной характеристики на частоту и характеризует распределение по частоте временного смещения выходного сигнала относительно входного. τ φ = arg ⁡ H ( j ω ) ω \tau_{\varphi}=\frac{\arg H(j\omega)}{\omega}

Временны́е характеристикиПравить

 
Типовые временные характеристики фильтра Чебышёва I рода 10-го порядка.

Временные характеристики фильтра Чебышёва I рода — импульсная переходная функция и переходная функция — представлены на рисунке. Импульсная переходная функция представляет собой реакцию фильтра на входной сигнал в виде дельта-функции Дирака, а переходная функция — реакцию на входное воздействие в виде единичной функции Хевисайда.

Фильтр Чебышёва II родаПравить

 
АЧХ фильтра Чебышёва II рода (фильтр низких частот) с ω 0 = 1 \omega_0=1 и ε = 0 , 01 \varepsilon=0,\!01

Фильтр Чебышёва II рода (инверсный фильтр Чебышёва) используется реже, чем фильтр Чебышёва I рода ввиду менее крутого спада амплитудной характеристики, что приводит к увеличению числа компонентов. У него отсутствуют пульсации в полосе пропускания, однако присутствуют в полосе подавления. Амплитудная характеристика такого фильтра задаётся следующим выражением: G n ( ω , ω 0 ) = 1 1 + 1 ε 2 T n 2 ( ω 0 / ω ) G_n(\omega,\;\omega_0) = \frac{1}{\sqrt{1+ \frac{1} {\varepsilon^2 T_n ^2 \left ( \omega_0 / \omega \right )}}}

В полосе подавления полиномы Чебышёва принимают значения от 0 до 1, из-за чего амплитудная характеристика такого фильтра принимает значения от нуля до 1 1 + 1 ε 2 \frac{1}{\sqrt{1+ \frac{1}{\varepsilon^2}}}

минимальной частотой, при которой достигается этот максимум является частота среза ω 0 \omega_0 . Параметр ε \varepsilon связан с затуханием в полосе подавления γ \gamma в децибелах следующим выражением: ε = 1 10 0 , 1 γ − 1 \varepsilon = \frac{1}{\sqrt{10^{0,\!1\gamma}-1}}

Для затухания на частотах полосы подавления в 5 дБ: ε = 0 , 6801 \varepsilon=0,\!6801 ; для затухания в 10 дБ: ε = 0 , 3333 \varepsilon=0,\!3333 . Частота f C = ω C / ( 2 π ) f_C=\omega_C/(2\pi) является частотой среза. Частота затухания в 3 дБ f H f_H связана с f C f_C следующим выражением: f H = f C ch ⁡ ( 1 n ch − 1 ⁡ 1 ε ) . f_H = f_C\,\mathop{\mathrm{ch}}\left(\frac{1}{n}\mathop{\mathrm{ch}}^{-1}\frac{1}{\varepsilon}\right).

Полюсы и нулиПравить

 
Логарифм модуля амплитудной характеристики фильтра Чебышёва II рода восьмого порядка на комплексной плоскости ( s = σ + j ω s=\sigma+j\omega ) с ε = 0 , 1 \varepsilon=0,\!1 и ω 0 = 1 \omega_0=1 . Белые пятна соответствуют полюсам, а чёрные — нулям. Показаны все 16 полюсов. 6 нулей (все нули второго порядка) показаны также, 2 находятся за пределами картинки (один на положительной мнимой оси, другой — на отрицательной мнимой оси). Полюсы передаточной функции фильтра — это полюсы, находящиеся в левой полуплоскости, нули передаточной функции — это нули модуля амплитудной характеристики фильтра Чебышёва, только не второго, а первого порядка. Чёрный цвет соответствует коэффициенту усиления менее 0,01, белый — коэффициенту усиления более 3.

Приняв частоту среза равной единице, получим выражение для полюсов ( ω p m ) (\omega_{pm}) фильтра Чебышёва: 1 + ε 2 T n 2 ( − 1 / j s p m ) = 0. 1+\varepsilon^2T_n^2(-1/js_{pm})=0.

Полюса фильтра Чебышёва II рода представляют собой «инверсию» полюсов фильтра Чебышёва I рода: 1 s p m ± = ± sh ⁡ ( 1 n arsh ⁡ ( 1 ε ) ) sin  Синус  ( θ m ) + \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)+ + j ch ⁡ ( 1 n arsh ⁡ ( 1 ε ) ) cos  Косинус  ( θ m ) , +j\mathop{\mathrm{ch}}\left(\frac{1}{n}\mathop{\mathrm{arsh}}\left(\frac{1}{\varepsilon}\right)\right)\cos(\theta_m) ,

где m = 1 , 2 , … , n m=1,\;2,\;\ldots,\;n .

Нули ( ω z m ) (\omega_{zm}) фильтра Чебышёва II рода определяются из следующего соотношения:: ε 2 T n 2 ( − 1 / j s z m ) = 0. \varepsilon^2T_n^2(-1/js_{zm})=0.

Нули фильтра Чебышёва II рода являются «инверсией» нулей многочленов Чебышёва: 1 / s z m = cos  Косинус  ( π 2 2 m − 1 n ) , 1/s_{zm} = \cos\left(\frac{\pi}{2}\,\frac{2m-1}{n}\right),

где m = 1 , 2 , … , n m=1,\;2,\;\ldots,\;n .

Передаточная функцияПравить

Передаточная функция задаётся при помощи полюсов в левой полуплоскости комлексной плоскости, её нули совпадают с нулями модуля амплитудной характеристики, с тем лишь отличием, что их порядок равен 1.

Групповая задержкаПравить

 
Амплитудная характеристика и групповая задержка фильтра Чебышёва II рода пятого порядка с ε = 0 , 1 \varepsilon=0,\!1 .

Амплитудная характеристика и групповая задержка показаны на графике. Можно видеть, что пульсации амплитуды приходятся на полосу подавления, а не на полосу пропускания.

Фазовые характеристикиПравить

 
Типовая ФЧХ и фазовая задержка фильтра Чебышёва II рода 10-го порядка.

Фазовые характеристики фильтра Чебышёва II рода — фазо-частотная характеристика и фазовая задержка — представлены на рисунке. Фазо-частотная характеристика показывает распределение по частоте смещения фазы выходного сигнала относительно входного. Фазовая задержка определяется как частное от деления фазо-частотной характеристики на частоту и характеризует распределение по частоте временного смещения выходного сигнала относительно входного.

Временные характеристикиПравить

 
Типовые временные характеристики фильтра Чебышёва II рода 5-го порядка.

Временные характеристики фильтра Чебышёва II рода — импульсная переходная функция и переходная функция — представлены на рисунке. Импульсная переходная функция представляет собой реакцию фильтра на входной сигнал в виде дельта-функции Дирака, а переходная функция — реакцию на входное воздействие в виде единичной функции Хевисайда.

Цифровые фильтры ЧебышёваПравить

Фильтры Чебышёва часто реализуются в цифровой форме. Для того, чтобы от аналогового фильтра перейти к цифровому, необходимо над каждым каскадом фильтра осуществить билинейное преобразование. Весь фильтр получается путём последовательного соединения каскадов. Простой пример фильтра Чебышёва низких частот I рода чётного порядка:

Z-преобразование каждого каскада: S ( Z ) = a ( Z ) b ( Z ) = α 0 + α 1 ⋅ Z − 1 + α 2 ⋅ Z − 2 1 + β 1 ⋅ Z − 1 + β 2 ⋅ Z − 2 . 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}}.

Во временной области преобразование записывается как: y [ n ] = α 0 ⋅ x [ 0 ] + α 1 ⋅ x [ − 1 ] + α 2 ⋅ x [ − 2 ] − β 1 ⋅ y [ − 1 ] − β 2 ⋅ y [ − 2 ] 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]

Коэффициенты α i \alpha_i \! и β i \beta_i \! подсчитываются из коэффициентов a i a_i \! и b i \! b_i : K = tg ⁡ ( π Frequency SampleRate ) K = \mathop{\mathrm{tg}}\left( \pi \frac{\mbox{Frequency}}{\mbox{SampleRate}}\right) temp i = cos  Косинус  ( 2 i + 1 ) π n \mbox{temp}_i =\cos\frac{(2i+1)\pi}{n} b i = 1 ch 2 ⁡ γ − temp i 2 b_i = \frac{1}{\mathop{\mathrm{ch}}^2\gamma-\mbox{temp}_i ^2} a i = K ⋅ b i ⋅ sh γ ⋅ 2 temp i a_i = K \cdot b_i \cdot \mathop{\mathrm{sh}}\,\gamma \cdot 2\,\mbox{temp}_i α 0 = K ⋅ K \alpha_0 = K \cdot K α 1 = 2 ⋅ K 2 \alpha_1 = 2 \cdot K^2 α 2 = K ⋅ K \alpha_2 = K \cdot K
β 0 ′ = ( a i + K 2 + b i ) \beta_0^\prime = (a_i + K^2 + b_i) β 1 ′ = 2 ⋅ ( b i − K 2 ) \beta_1^\prime = 2 \cdot (b_i - K^2) β 2 ′ = ( a i − K 2 − b i ) \beta_2^\prime = (a_i - K^2 - b_i)
β 1 = β 1 ′ / β 0 ′ \beta_1 = \beta_1^\prime / \beta_0^\prime β 2 = β 2 ′ / β 0 ′ \beta_2 = \beta_2^\prime / \beta_0^\prime

Для получения фильтра Чебышёва более высокого порядка, необходимо соединить последовательно несколько каскадов.

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

Ниже представлены графики АЧХ фильтра Чебышёва I и II родов в сравнении с некоторыми другими фильтрами с тем же числом коэффициентов:

 

По графикам видно, что амплитудная характеристики фильтров Чебышёва имее более крутой спад, чем у фильтров Баттерворта, но не такой крутой, как у эллиптического фильтра.

См. такжеПравить

БиблиографияПравить

  • В.А. Лукас Теория автоматического управления. — M.: Недра, 1990.о книге
  • Б.Х. Кривицкий Справочник по теоретическим основам радиоэлектроники. — М.: Энергия, 1977.о книге
  • Richard W. Daniels Approximation Methods for Electronic Filter Design. — New York: McGraw-Hill, 1974. — ISBN 0070153086о книге
  • Steven W. Smith The Scientist and Engineer’s Guide to Digital Signal Processing. — Second Edition. — San-Diego: California Technical Publishing, 1999. — ISBN 0966017641о книге
  • Britton C. Rorabaugh Approximation Methods for Electronic Filter Design. — New York: McGraw-Hill, 1999. — ISBN 0070540047о книге
  • B. Widrow, S.D. Stearns Adaptive Signal Processing. — Paramus, NJ: Prentice-Hall, 1985. — ISBN 0130040290о книге
  • S. Haykin Adaptive Filter Theory. — 4rd Edition. — Paramus, NJ: Prentice-Hall, 2001. — ISBN 0130901261о книге
  • Michael L. Honig, David G. Messerschmitt Adaptive Filters — Structures, Algorithms, and Applications. — Hingham, MA: Kluwer Academic Publishers, 1984. — ISBN 0898381630о книге
  • J.D. Markel, A.H. Gray, Jr. Linear Prediction of Speech. — New York: Springer-Verlag, 1982. — ISBN 0387075631о книге
  • L.R. Rabiner, R.W. Schafer Digital Processing of Speech Signals. — Paramus, NJ: Prentice-Hall, 1978. — ISBN 0132136031о книге
  • Richard J. Higgins Digital Signal Processing in VLSI. — Paramus, NJ: Prentice-Hall, 1990. — ISBN 013212887Xо книге
  • A. V. Oppenheim, R. W. Schafer Digital Signal Processing. — Paramus, NJ: Prentice-Hall, 1975. — ISBN 0132146355о книге
  • L. R. Rabiner, B. Gold Theory and Application of Digital Signal Processing. — Paramus, NJ: Prentice-Hall, 1986. — ISBN 0139141014о книге
  • John G. Proakis, Dimitris G. Manolakis Introduction to Digital Signal Processing. — Paramus, NJ: Prentice-Hall, 1988. — ISBN 002396815xо книгеРегулярное выражение «ISBN» классифицировало значение «002396815x» как недопустимое.

СсылкиПравить


Первоисточник этой статьи был признан «хорошей статьёй» русского раздела Википедии.
Источник — https://traditio.wiki/w/index.php?title=Фильтр_Чебышёва&oldid=58870