Проектирование симметричных полосовых фильтров и развязанных резонансных звеньев

№ 12’2008
PDF версия
В статье приводится методика расчета полосовых фильтров с симметричной АЧХ полиномиального типа Баттерворта и Чебышева в виде каскадного соединения развязанных резонансных звеньев. Методика рассчитана на произвольные параметры задаваемой АЧХ. При реализации могут быть использованы звенья любой схемотехники. Приведены примеры наиболее популярных схем звеньев и даны методы их расчета.

Постановка задачи проектирования полосового фильтра

Задача проектирования полосового фильтра (ПФ) задается, как правило, параметрами амплитудно–частотной характеристики (АЧХ), представленной в обобщенном виде на рис. 1.

Рис. 1. АЧХ полосового фильтра в терминах обобщенных параметров
Рис. 1. АЧХ полосового фильтра в терминах
обобщенных параметров

Данная АЧХ характеризуется следующими параметрами:

  • ωНП, ωВП — соответственно нижняя и верхняя границы полосы пропускания;
  • λ — величина допустимой неравномерности АЧХ ПФ в пределах полосы пропускания;
  • ωНЗ, ωВЗ — соответственно нижняя и верхняя границы зон задерживания;
  • ω0 — средняя частота полосы пропускания;
  • КНЗ, КВЗ — соответственно нижнее и верхнее значение коэффициента передачи в зонах задерживания.

Если величина отношения fВП/fНП<5, то такой ПФ относится к классу узкополосных. Его оптимальная реализация с точки зрения минимально необходимого числа элементов строится в виде каскадного соединения развязанных резонансных звеньев (контуров). Очевидно, что это предопределяет полиномиальный характер реализуемой АЧХ. Такие ПФ называют симметричными полосовыми фильтрами (СПФ). Симметричность АЧХ в данном случае подразумевается в геометрическом смысле и определяется выполнением следующих условий:

 

omega_0 = sqrt(omega_НП x omega_ВП) = sqrt(omega_НЗ x omega_ВЗ),

(1)
  КНЗ = КВЗ. (2)

Простейшим СПФ является резонансное звено 2–го порядка, передаточная функция которого имеет вид:

 

T(j omega)=M(j omega/omega_0 x d)/(1+j omega/omega_0 x d-j omega^2/omega_0^2),

(3)

где М — масштабный множитель (коэффициент плоского усиления); d — затухание, величина, обратная добротности резонансного контура:

  1/d = Q = ω0/(ωВП−ωНП) = ω0/Δω. (4)

Здесь ширина полосы пропускания (fВП−fНП) определяется по уровню 0,707 (−3 дБ).

Решение задачи аппроксимации — определение параметров передаточной функции — производят через ФНЧ–прототип. Для этого передаточную функцию ПФ преобразуют с помощью реактансного (частотного) преобразования в передаточную функцию ФНЧ–прототипа. Покажем этот переход на примере передаточной функции одиночного резонансного контура. Для этого разделим числитель и знаменатель выражения (3) на числитель, в результате чего получим:

 

T(j omega)=M 1/(1+j Q(omega/omega_0-omega_0/omega),

(5)

Введем в рассмотрение новую переменную, называемую относительной расстройкой:

  Y = (ω/ω0)−(ω0/ω). (6)

Данное выражение и описывает реактансное, или частотное преобразование.

В результате такого перехода к новой переменной передаточная функция резонансного конура второго порядка (5) преобразуется в передаточную функцию ФНЧ–прототипа первого порядка:

  T(jY) = M(1/(1+jQY)). (7)

Таким образом, формально решение задачи проектирования СПФ в виде каскадного соединения развязанных резонансных контуров сводится на первом этапе к определению порядка передаточной функции ФНЧ–прототипа, что с учетом билинейного характера частотного преобразования (6) определяет требуемое число резонансных контуров, а также добротностей и частот настроек контуров. На втором этапе реализации осуществляется выбор схемы резонансного звена и рассчитываются по соответствующим методикам величины его элементов.

Рассмотрим методику решения первой задачи, ограничиваясь полиномиальными АЧХ чебышевского и баттервортовского типов. Полиномиальные АЧХ характеризуются монотонным характером спада за пределами полосы пропускания. В пределах полосы пропускания чебышевская АЧХ имеет равноволновый колебательный характер, а баттервортовская — монотонный.

Порядок расчета чебышевского СПФ

Порядок расчета СПФ, заданного параметрами АЧХ (рис. 1), может быть представлен в виде следующего вычислительного алгоритма.

  1. Производится симметризация (в геометрическом смысле) заданных граничных частот АЧХ из условия выполнения соотношения (1).
    Примечание. Если строго задана центральная частота и ширина полосы пропускания, то одна из границ полосы уточняется в соответствии с равенством (1). Если же строго задана только ширина полосы пропускания, то точно вычисляется значение центральной частоты, и при выбранном значении ω или ωВЗ вычисляется другая частота границ зон задерживания. Как правило, за исходное значение принимается ωВЗ. Тогда:

      ω = ω02ВЗ.
    (8)
  2. Осуществляется переход к параметрам АЧХ ФНЧ–прототипа на основе частотного преобразования (1):
      YВП = (ωВП0)−(ω0ВП). (9)
      YВЗ = (ωВЗ0)−(ω0ВЗ). (10)

    Очевидно, что при этом преобразовании будет выполняться следующее условие:

      YВЗ = (ωВЗ0)−(ω0ВЗ) = −YНЗ =

    = (ω0)−(ω0).
    (11)
  3. Осуществляется переход к нормированной частоте АЧХ ФНЧ–прототипа:
      Ω = Y/YВП. (12)

    При этом получим:

      ΩВП = YВП/YВП = 1, ΩВЗ = YВЗ/YВП. (13)
  4. Определяется порядок ФНЧ–прототипа — количество необходимых резонансных контуров. Фактически здесь определяется порядок полинома Чебышева, использующегося для аппроксимации АЧХ ФНЧ–прототипа.
     

    (14)

    Знак неравенства в данном выражении предполагает округление получаемой величины до ближайшего большего целого числа. Очевидно, что такое округление предопределяет образование некоторого «запаса» по величине КВЗ или ΩВЗ по сравнению с исходно заданными значениями.
    Примечание. Если выполняются условия: КВЗ2 << 1; λ << 1, то можно вместо (14) воспользоваться упрощенной формулой:

     

    (15)

    Возникающая при этом ошибка не превышает 0,6% даже при КВЗ = 0,3 и λ = 0,3.

  5. Уточняется величина КВЗ вследствие округления n:
     

    K_ВЗ=1/sqrt(0,5lambda(Omega_ВЗ+sqrt(Omega_ВЗ^2-1))^(2n)+1).

    (16)

    Примечание. Уточнение величины КВЗ предполагает сохранение неизменными величин λ и ΩВЗ. Но если неизменными величинами считать λ и КВЗ, то уточнять следует величину ΩВЗ. Для этого необходимо решить уравнение:

      ΩВЗ = (α2+1)/2α;

    где α = (((1/KВЗ2)−1)/(0,5×λ))2n.
    (17)
  6. Определяется число звеньев ФНЧ–прототипа:
    Примечание. Звено первого порядка ФНЧ–прототипа соответствует одиночному резонансному контуру СПФ, настроенному на частоту ω0. Звено второго порядка ФНЧ–прототипа соответствует двум резонансным контурам СПФ с одинаковой добротностью, симметрично расстроенным влево и вправо относительно ω0.
  7. Определяются параметры звеньев второго порядка ФНЧ–прототипа:
      dq = 2Dsin(ψq), (18)
      Aq = B2cos2q)+D2sin2q), (19)
      B = 0,5(E+1/E), D = 0,5(E−1/E), (20)
     

    E=sqrt[2n]((1+lambda)/lambda+sqrt(((1+lambda)/lambda)^2-1)),

    (21)
      ψq = (π/2n)×[1+2(k−1)], k = 1, 2, … 2n. (22)

    Фактически значения углов ψq определяются только для первого квадранта.
    Примечание. Если выполняется условие λ << 1, тогда вместо (21) можно воспользоваться упрощенной формулой:

     

    E=sqrt[2n](2(1+lambda)/lambda).

    (23)

    Возникающая при таком упрощении ошибка может быть оценена следующими примерами:

    • при λ = 0,1−δE = 0,1%;
    • при λ = 0,3−δE = 0,5%.
  8. Определяются параметры (частоты настроек и добротность) двух симметрично расстроенных резонансных контуров СПФ, соответствующих звену второго порядка ФНЧ–прототипа.
    Значение добротности определяется следующим выражением:
     

    Q_q=sqrt((A_q+4/Y_ВП^2)/(2d_q^2)+sqrt(((A_q+4/Y_ВП^2)/(2d_q^2))^2-1/d_q^2xY_ВП^2)).

    (24)

    Частоты настроек определяются согласно

      ωpq′= Xqω0, ωpq′= ω0/Xq, (25)
     

    X_q=0,5xd_qQ_qY_ВП+sqrt((0,5xd_qQ_qY_ВП)^2-1

    (26)
  9. Определяются параметры резонансного контура, соответствующего звену первого порядка ФНЧ–прототипа. Частота настройки этого контура равна центральной частоте СПФ — ω0, а его добротность определяется соотношением:
      QS = 1/(D×YВП). (27)

    Примечание. Если проектируемый СПФ является узкополосным, то есть если выполняется условие

    ωВПНП ≤ 1,2…1,25,

    тогда при расчете добротности и частот настроек пары расстроенных контуров вместо формул (24) и (26) можно использовать упрощенные выражения:

      Qq = (Xq/YВП)×(1/(Dsin(ψq)); (28)
      Xq = 1+0,5YВПBcos(ψq). (29)

Порядок расчета баттервортовского СПФ

Расчет баттервортовского СПФ, АЧХ которого имеет монотонный характер в пределах полосы пропускания, в целом аналогичен рассмотренному выше. Отличие обусловлено лишь двумя обстоятельствами. Во–первых, решение задачи аппроксимации основано на применении полиномов Баттерворта. Поэтому порядок фильтра определяется иным выражением. Во–вторых, АЧХ баттервортовского ФНЧ–прототипа формально имеет неизменную величину неравномерности в полосе пропускания −3 дБ. Тем не менее, СПФ с баттервортовской АЧХ можно рассчитать для произвольного значения неравномерности. Все это отражено в приводимом далее порядке расчета, который основан на задании параметров АЧХ СПФ в виде рис. 1.

  1. В соответствии с (1) производится симметризация частотных параметров:
     

    omega_0 = sqrt(omega_ВПxomega_НП) = sqrt(omega_ВЗxomega_НЗ) = sqrt(omega_ВП*xomega_НП*).

    (30)

    Здесь введены две частоты, определяющие границы полосы пропускания для произвольного значения неравномерности λ (помечены значком *). Частоты ωНП и ωВП определены для неравномерности −3 дБ (λ = 0,293).

  2. Переход к параметрам АЧХ ФНЧ–прототипа осуществляется аналогично предыдущему, за тем исключением, что в данном случае вычисляем
      YВП* = (ωВП*/ω0)−(ω0ВП*). (31)

    Далее — по формулам (10), (13).

  3. Определение порядка ФНЧ–прототипа (числа резонансных контуров) производится в данном случае по следующей формуле:
     

    (32)

    Здесь также знак неравенства предполагает округление до ближайшего большего целого числа.

  4. Уточнение величины КВЗ вследствие округления n:
     

    К_ВЗ=1/sqrt(1+Omega_ВЗ^(2n)).

    (33)

Дальнейший порядок расчета полностью соответствует приведенному выше чебышевскому варианту СПФ: пункты 6, 7, 8, 9 сохраняются полностью с одним уточнением, связанным с необходимостью для данных пунктов вычисления значения YВП для λ = 0,293 (−3 дБ) в соответствии с выражением:

 

Y_ВП=Y_ВП*/sqrt[2n](2lambda).

(34)

Резонансные ARC–звенья и их расчет

В результате расчета по приведенной выше методике определено требуемое число резонансных контуров и их параметры: добротности и частоты настроек. Иными словами, полностью определена структура СПФ в виде рис. 2.

Рис. 2. Структура СПФ
Рис. 2. Структура СПФ

При каскадировании следует располагать звенья в порядке нарастания величины их добротности.

Вторым этапом проектирования СПФ является выбор схемы резонансного звена и расчет величин его элементов по заданным значениям добротности и частоты настройки.

Резонансные ARC–звенья подразделяются на низкодобротные (Q < 5), среднедобротные (Q < 20) и высокодобротные (Q > 20). Количество схемотехнических решений резонансных звеньев достаточно велико и насчитывает много десятков известных вариантов, построенных по различной идеологии [1–3]. В рамках данной работы ограничимся лишь несколькими примерами практических схем звеньев.

Низкодобротное резонансное звено на ОУ с отрицательной обратной связью [3]

Схема звена приведена на рис. 3.

Рис. 3. Схема низкодобротного резонансного звена
Рис. 3. Схема низкодобротного резонансного звена

Коэффициент передачи звена имеет вид, аналогичный выражению (3):

 

T(j omega)=-Mx(jomega/omega_0xd)/(1+jomega/omega_0xd-omega^2/omega_0^2).

(35)

Для случая равноемкостного варианта схемы С1 = С2 = С основные параметры звена определяются следующими соотношениями:

  • частота настройки:
     

    omega_0 = 1/sqrt(rxR_3xC^2);

    (36)
  • добротность:
     

    Q = 1/d = 0,5sqrt(R_3/r);

    (37)
  • масштабный множитель (коэффициент плоского усиления):
      M = R2/(R1+R2)×2Q2 = (r/R1)×2Q2. (38)

В этих выражениях

  r = (R1×R2)/(R1+R2). (39)

Отметим, что величина плоского усиления в данной схеме прямо пропорциональна квадрату добротности.

Расчет схемы сводится к определению величин элементов по заданным значениям Q, M и ω0 = 2πf0.

Методика расчета в предположении идеальности ОУ состоит в следующем.

  1. Из конструктивных и технологических соображений выбираем величину емкости С.
  2. На основе выражения (36) вычисляем вспомогательную величину α:
     

    alpha = sqrt(rxR_3 = 1/(omega_0xC).

    (40)

    Затем вычисляем:

      R3 = 2Qα, (41)
      r = α/Q. (42)
  3. Вычисляем:
      R1 = (r/M)×2Q, (43)
      R2 = (R1×r)/(R1−r). (44)

Замечание: при неудачном выборе величины С последнее выражение может дать отрицательную величину сопротивления R2. В этом случае следует вернуться на шаг № 1 и выбрать другое значение С.

Резонансное звено на основе обобщенного конвертора импеданса

Схема звена, приведенная на рис. 4, позволяет реализовать резонансную характеристику с добротностью от единиц до двух сотен в диапазоне частот до 1 МГц в зависимости от типа используемых ОУ[4].

Рис. 4. Резонансное звено на обобщенном конверторе импеданса
Рис. 4. Резонансное звено
на обобщенном конверторе импеданса

Передаточная функция звена имеет вид (3), где:

  • масштабный множитель (коэффициент плоского усиления):
      M = 1+(R1/RH); (45)
  • частота настройки (резонанса):
     

    omega_0=1/sqrt(CxC_3xR_2xR_4xR/R_1);

    (46)
  • добротность:
     

    Q=R sqrt(C/C_3xR_1/(R_2xR_4xR_Н)).

    (47)

Поскольку частота настройки не зависит от R, изменяя его, можно изменять добротность при сохранении неизменными значений частоты.

Порядок расчета звена в предположении идеальности ОУ состоит в следующем.

Задано: M, Q, ω0. Полагаем С3 = С — равноемкостный вариант.

  1. Вычисляем вспомогательный параметр β:
      β = M−1. (48)
  2. Вычисляем оценочную величину резистора R2:
      R2 = 1/(ω0×C). (49)

    Полученное значение округляем до ближайшего удобного номинала из ряда.

  3. Рассчитываем вспомогательную величину r:
      r = √β/(ω0×C). (50)
  4. Вычисляем:
      R4 = r2/R2, (51)
      R = Q(r/√β). (52)
  5. Резисторы R1 и Rн входят в выражения (45–47) своим отношением. Поэтому далее выбираем R1 из ряда номиналов величиной, близкой к величине R2, и затем определяем R:
      Rн = R1/β. (53)

Универсальное звено второго порядка

Другой схемой высокодобротного резонансного звена является универсальное звено второго порядка, реализованного на основании метода аналогового моделирования (рис. 5)[1].

Рис. 5. Схема универсального звена второго порядка
Рис. 5. Схема универсального звена второго порядка

Данная схема в зависимости от точки съема выходного сигнала реализует передаточную функцию резонансного контура (ВЫХ ПФ) или два варианта передаточной функции ФНЧ второго порядка (соответственно ВЫХ НЧ1 и ВЫХ НЧ2), которые различаются видом коэффициента числителя передаточной функции. Реализуемые звеном передаточные функции имеют соответственно следующий вид:

  TПФ (jω) = M1((j(ω/ω0)×d)/(P(jω))); (54)
  TНЧ1(jω) = M2 (1/(P(jω))); (55)
  TНЧ2(jω) = M3 (1/(P(jω))); (56)

где:

  • P(jω) = 1+j(ω/ω0)×d−(ω202) (57) — полином знаменателя второго порядка;
  • ω0 = 1/RC (58) — частота настройки звена (для ПФ — это частота резонанса);
  • d = R/R6 = 1/Q (59) — затухание звена, величина, обратная добротности.

Частота настройки и затухание определены при условии равнономинальности:

  • C1 = C2 = C;
  • R1 = R2 = R3 = R4 = R5 = R;
  • масштабный множитель ПФ — M1 = R6/R3; (60)
  • масштабный множитель ФНЧ–1 — M2 = (R1×R4)/(R3×R5); (61)
  • масштабный множитель ФНЧ–2 — M3 = R1/R3. (62)

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

  • задано М, Q, ω0 = 2πf0;
  • выбираем из конструктивных и/или технологических соображений величину емкости конденсатора С (желательно из ряда номиналов);
  • вычисляем на основе (58) величину R:
      R = 1/2πf0C; (63)
  • на основе (59) вычисляем величину сопротивления R6:
      R6 = R×Q; (64)
  • на основе (60) вычисляем R3:
      R3 = R6/M = R(Q/M). (65)

Примечание. Универсальные звенья рассмотренного типа в настоящее время выпускаются в виде микросхем рядом фирм, производящих микросхемы. Примером могут служит микросхемы фирмы MAXIM MAX274, MAX275. Микросхема МАХ274 содержит в одном корпусе два звена второго порядка, каждое из которых позволяет реализовать звено ФНЧ–2 или резонансное звено ПФ–2. Микросхема МАХ275 содержит в одном корпусе четыре таких универсальных звена.

Подробную информацию об этих и других микросхемах активных фильтров можно найти на сайте данной фирмы [5].

Активное LCR–звено резонансного контура второго порядка

При реализации СПФ, предназначенного для работы в области достаточно высоких частот — сотни килогерц, единицы мегагерц, требующих к тому же больших значений добротности составляющих резонансных контуров, вполне допустимо применение катушек индуктивности. Для таких частот современная промышленность выпускает широкий ассортимент готовых индуктивных компонентов с весьма высокими показателями, малыми габаритами (выпускаются даже катушки для поверхностного монтажа) и относительно низкой стоимостью.

Простейшим примером реализации резонансного контура в этом случае может служить схема, приведенная на рис. 6.

Рис. 6. Резонансное LCR-звено с усилением
Рис. 6. Резонансное LCR–звено с усилением

Передаточная функция данного звена имеет вид (3).

Частота настройки определяется параметрами LC–резонансного контура:

 

omega_0 = 2pi f_0 = 1/sqrt(LC).

(66)

Добротность резонансного контура есть отношение сопротивления R1 и характеристического сопротивления ρ:

 

Q = R/rho = R/sqrt(L/C).

(67)

Усилитель на ОУ выполняет роль буферного усилителя для обеспечения развязки при каскадном включении таких звеньев. Кроме того, он может использоваться как масштабный усилитель, причем результирующий коэффициент усиления фильтра в целом можно распределить между всеми входящими в фильтр усилительными фрагментами.

Коэффициент усиления отдельного усилителя в звене и, следовательно, масштабный множитель данного звена не влияет ни на частоту настройки, ни на добротность и определяется простым соотношением:

  K = 1+R3/R2. (68)

Заключение

Предложенная в статье методика, основанная на развитии идей, заложенных в работе[6], представляет собой хорошо алгоритмизированную последовательность простых вычислений. Вычислительная процедура легко реализуется в пакете «MathCAD» и позволяет определять с высокой точностью основные параметры резонансных контуров (частоты настроек и добротности) практически для любого порядка СПФ с полиномиальными АЧХ. Ограничения налагаются только на величину максимальной добротности требуемого резонансного контура, поскольку это определяется реализационными возможностями выбираемой схемотехники.

В приложении 1 приведен пример расчета чебышевского СПФ 6–го порядка, реализованный в пакете «MathCAD 2001», который можно найти на сайте журнала http://www.kite.ru/assets/rasschetChebyshevsky.mcd.

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

Литература

  1. Хейнлейн В. Е., Холмс В. Х. Активные фильтры для интегральных схем. М.: Связь, 1980.
  2. Справочник по расчету и проектированию ARCсхем / Под ред. проф. А. А. Ланнэ. М.: Радио и связь, 1984.
  3. Мошитц Г., Хорн П. Проектирование активных фильтров. М.: Мир, 1984.
  4. Коротков А. С., Михалев П. Г. А. с. 1385260 (СССР). Полосовой активный RC–фильтр // БИ. 1988. № 12.
  5. www.maxim-ic.com
  6. Славский Г. Н. Активные RC– и RCL–фильтры и избирательные усилители. М.: Связь, 1966.

Добавить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *