Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Методы обработки данных в научных исследованиях. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
07.09.2026
Размер:
2 Мб
Скачать
☆
ВЧ-фильтр будет иметь частоту среза 0,5*fs − 0,2*fs = 0,3*fs. Это должно приниматься во внимание при проектировании исходного НЧ-фильтра.
Проектирование ВЧ-фильтра с использованием
импульсной характеристики НЧ-фильтра
Метод инверсии спектра:
проектируем ФНЧ (линейная фаза, нечетное N); изменяем знак каждого коэффициента в импульсной ха-
рактеристике h(n);
 прибавляем 1 к коэффициенту в центре симметрии h(n).
Метод реверсирования спектра:
проектируем ФНЧ; изменяем знак каждого второго коэффициента в импульс-
ной характеристике h(n);
это переворачивает характеристику в частотной области
слева направо: значение по частоте, соответствующее 0, становится соответствующим 0,5, а 0,5 — соответствую­щим 0; т. е. если частота среза ФНЧ равна 0,2, то частота среза результирующего ФВЧ равна 0,3.
Рис. 12. Проектирование полосовых и режекторных фильтров
с использованием ФНЧ и ФВЧ
31
Полосовой и режекторный фильтры можно спроектировать, комбинируя надлежащим образом соответствующие НЧ- и ВЧ­фильтры (рис. 12). Полосовые фильтры проектируются посредством каскадного соединения НЧ- и ВЧ-фильтров. Вычисляя свертку двух индивидуальнх импульсных характеристик, получают эквивалент­ную импульсную характеристику каскадных фильтров.
Режекторный фильтр проектируется посредством параллельного подключения НЧ- и ВЧ-фильтров и суммирования сигналов с их выходов. Складывая индивидуальные импульсные характеристики, получают эквивалентную импульсную характеристику.
4.6. Фильтры с бесконечной импульсной характеристикой
БИХ-фильтры получили такое название, потому что их импульс­ные характеристики растянуты на бесконечном временном интерва­ле. Это объясняется тем, что данные фильтры являются рекурсив­ными, т. е. используют обратную связь. Хотя БИХ-фильтры могут реализовываться с меньшим, чем КИХ-фильтры, количеством вы­числений, у них нет таких характеристик, которыми обладают КИХ-фильтры. Более того, БИХ-фильтр не имеет линейной фазовой характеристики. Но вычислительные преимущества БИХ-фильтра теряются, когда выходной сигнал фильтра подвергается децимации, поскольку в этом случае всякий раз приходится вычислять заново значение выходной величины. БИХ-фильтры обычно реализуются с помощью звеньев второго порядка, которые называются биквадрат­ными фильтрами, потому что описываются биквадратными уравне­ниями в z-области. Фильтры высокого порядка проектируют, ис­пользуя каскадирование биквадратных звеньев. Например, фильтр шестого порядка требует трех биквадратных звеньев.
Фильтры с бесконечной импульсной характеристикой
Есть обратная связь (рекурсия).
Импульсная характеристика имеет бесконечную длитель-
ность.
Потенциально нестабильны.
Нелинейная фазочастотная характеристика.
Более эффективны, чем КИХ-фильтры.
Нет вычислительных преимуществ при децимации по выходу.
32
Обычно проектируются по характеристике аналогового филь-
тра и каскадным соединением звеньев второго порядка (би­квадратные фильтры).
Структура биквадратного БИХ-фильтра представлена на рис. 13. Нули формируются коэффициентами прямой связи b люса (порядок) — коэффициентами обратной связи a
, b1 и b2; а по-
0
и a2.
1
Рис. 13. Аппаратная реализация БИХ-фильтра
второго порядка (биквадратного) — форма 1
Общее уравнение цифрового фильтра описывает обобщенную передаточную функцию H(z), которая содержит полиномы и в чис­лителе, и в знаменателе. Корни знаменателя определяют располо­жение полюсов фильтра, а корни числителя характеризуют распо­ложение нулей. Хотя существует возможность создания непосред­ственно по этому уравнению БИХ-фильтра более высокого порядка (так называемая прямая реализация), накапливающиеся ошибки квантования (из-за арифметики с фиксированной точкой и конечной длины слова) могут вызывать неустойчивость работы фильтра и большие ошибки. По этой причине правильнее расположить кас-
33
кадно несколько биквадратных звеньев с соответствующими коэф­фициентами, чем использовать прямую форму реализации. Данные при вычислении биквадратных фильтров могут масштабироваться раздельно, а затем биквадратные звенья каскадируются для мини­мизации ошибок квантования коэффициентов и накапливающихся ошибок рекурсивного накопления. Каскадные биквадратные филь­тры работают более медленно, чем их эквиваленты прямой формы реализации, но они более устойчивы и в них минимизируются эф­фекты, связанные с арифметическими ошибками конечной разряд­ности данных.
Для первой прямой формы биквадратного звена необходимо ис­пользование четырех регистров. Эта конфигурация может быть за­менена эквивалентной схемой, которая называется второй прямой формой реализации и требует применения только двух регистров. Можно показать, что уравнения, описывающие биквадратный БИХ­фильтр второй прямой формы реализации, такие же, как и уравне­ния первой прямой формы реализации (рис. 14).
Рис. 14. Биквадратный БИХ-фильтр (форма 2)
4.7. Методы проектирования БИХ-фильтров
Популярный метод проектирования БИХ-фильтра сводится к тому, что сначала проектируется эквивалентный аналоговый
34
фильтр, а затем функция передачи H(s) преобразуется математи­чески в z-область, H(z). Проектирование фильтров более высоких порядков выполняется каскадированием биквадратных звеньев. Наиболее популярными аналоговыми фильтрами являются филь­тры Баттерворта, Чебышева, эллиптические и Бесселя.
Фильтр Баттерворта, не имеющий нулей частотной характе­ристики (также называемый фильтром с максимально плоской характеристикой), не создает пульсаций (неравномерности) в по­лосе пропускания и полосе задержки, т. е. обладает монотонной характеристикой в обеих полосах.
Фильтр Чебышева 1-го рода имеет более быстрый спад ча­стотной характеристики, чем фильтр Баттерворта (при равном порядке), и создает пульсации (неравномерность) в полосе про­пускания. Реже используются фильтры Чебышева 2-го рода с пульсациями (неравномерность) в полосе задержки, а не в по­лосе пропускания.
Эллиптический фильтр (фильтр Кауэра) имеет полюса и нули частотной характеристики и создает пульсации (неравномер­ность) и в полосе пропускания, и в полосе задержки. У этого фильтра более быстрый спад частотной характеристики, чем фильтр Чебышева, при том же числе полюсов (порядке). Эллип­тический фильтр часто используется там, где допускается не­сколько худшая фазовая характеристика.
Фильтр Бесселя (Томпсона), который не имеет нулей частот- ной характеристики, обладает оптимальной импульсной характе­ристикой и линейной фазовой характеристикой, но у него худ­ший спад частотной характеристики из всех типов обсуждавших­ся фильтров при том же числе полюсов (порядке).
Обзор популярных фильтров
Фильтр Характеристика
Нет нулей частотной характеристики и пульса-
Баттерворта
ций в полосе пропускания и задержки Максимально плоская характеристика (быстрый
спад без пульсаций)
35
Нет нулей частотной характеристики, пульса­ции в полосе пропускания и полосе задержки
Чебышева
1-го рода
Эллиптический
(Кауэра)
Бесселя
(Томпсона)
Выбор между КИХ- и БИХ-фильтрами может быть своего рода состязанием в проектировании, но несколько основных руководя­щих принципов дать можно. Как правило, БИХ-фильтры более эф­фективны, потому что требуют меньшего количества памяти и опе­раций умножения с накоплением. БИХ-фильтры могут быть разра­ботаны, основываясь на предыдущем опыте проектирования анало­говых фильтров. БИХ-фильтры могут приносить проблемы не­устойчивости, но это происходит реже, если проектируемые филь­тры высокого порядка реализуются как системы, состоящие из кас­кадов второго порядка.
С другой стороны, КИХ-фильтры требуют большего количества звеньев и, соответственно, операций умножения с накоплением для реализации частотной характеристики с заданной частотой среза, но при этом имеют линейную фазовую характеристику. КИХ-фильтры работают на конечном временном интервале данных, поэтому, если
Более короткая область перехода, чем у фильтра Баттерворта для данного порядка
Фильтр 2-го рода имеет пульсации в полосе за­держки, нет пульсаций в полосе пропускания
Имеет полюса и нули, пульсации и в полосе пропускания, и в полосе задержки
Более короткая область перехода, чем у фильтра Чебышева для данного порядка
Фазовая характеристика хуже Нет нулей частотной характеристики, нет пуль-
саций в полосе пропускания и задержки Оптимизирован по линейной фазовой и им-
пульсной характеристикам Самая длинная переходная характеристика из
всех фильтров данного порядка
4.8. Сравнение КИХ- и БИХ-фильтров
36
часть данных испорчена (например, в результате сбоев в работе АЦП), КИХ-фильтр будет «звенеть» только на временном интерва­ле, соответствующем N-1 отсчетам. БИХ-фильтр из-за наличия об­ратной связи будет «звенеть» значительно более длительный пери­од времени.
Сравнение КИХ- и БИХ-фильтров
КИХ-фильтр БИХ-фильтров
Более эффективны Менее эффективны Есть аналоговый эквивалент Нет аналогового эквивалента Могут быть нестабильными Всегда стабильны Нелинейная фазовая характе-
ристика Больше «звон» при наличии
ложных сигналов Доступны средства САПР Доступны средства САПР Децимация не влияет на эффек-
тивность
Если необходимы фильтры с крутым спадом и испытывается де­фицит во времени, отведенном для обработки, хороший выбор — эллиптические БИХ-фильтры. Если число операций умножения с накоплением не является чрезмерным и требуется линейная фаза, то следует использовать КИХ-фильтр.
Линейная фазовая характери­стика
Меньше «звон» при наличии ложных сигналов
Децимация увеличивает эффек­тивность
37
5. ПЕРЕДИСКРЕТИЗАЦИЯ ЦИФРОВЫХ СИГНАЛОВ
f
f
Передискретизация в цифровой обработке сигналов — измене­ние частоты дискретизации цифрового сигнала. Алгоритмы пере­дискретизации широко применяются при обработке звуковых сиг­налов, радиосигналов и изображений (передискретизация растрово­го изображения — это изменение его разрешения в пикселах).
Отсчеты сигнала, соответствующие новой частоте дискретиза­ции, вычисляются по уже имеющимся отсчетам и не содержат но­вой информации. Повышение частоты дискретизации называется интерполяцией, понижение — децимацией.
Согласно теореме Котельникова, любой непрерывный сигнал с конечным спектром (в котором спектральные составляющие соот­ветствуют частотам, выше или равным некоторой частоты
) мо-
0
жет быть представлен в виде отсчетов дискретного сигнала с часто-
2
той дискретизации
f . При этом такое преобразование явля-
0
d
ется взаимно однозначным, т. е. при соблюдении условий теоремы Котельникова по дискретному сигналу можно восстановить исход­ный сигнал с финитным спектром без искажений.
При передискретизации отсчеты сигнала, соответствующие од­ной частоте дискретизации, вычисляются по имеющимся отсчетам этого же сигнала, соответствующим другой частоте дискретизации (при этом предполагается, что обе частоты дискретизации соответ­ствуют условиям теоремы Котельникова). Идеальная передискрети­зация эквивалентна восстановлению непрерывного сигнала по его отсчетам с последующей дискретизацией его на новой частоте.
Рис. 15. Пример вычисления промежуточного отчета
в точке −0,5 с помощью идеального фильтра нижних частот
38
На рис. 15 показано вычисление промежуточного отсчета (в точке
x
x
−0,5) дискретного сигнала с помощью идеального фильтра нижних частот. Синяя кривая — исходный непрерывный сигнал, зеленая — импульсная характеристика идеального ФНЧ. Для интерполяции суммируются значения импульсной характеристики в точках отсчетов, умноженные на величину соответствующих отсчетов.
Точное вычисление значения исходного непрерывного сигнала в определенной точке производится следующим образом:
w
d
sin ( )
() ( ) ,
xt xt
()
t — i-й отсчет сигнала; it — момент времени, соответству-
где ющий этому отсчету; 2
i
wf
dd
i
i
— циклическая частота дискрети-
2
w
d
()
tt
2
tt
i
i
зации; x(t) — интерполированное значение сигнала в момент вре­мени t.
Функция
w
d
sin ( )
2
w
d
()
tt
2
tt
i
i
не является финитной, поэтому для вычисления значения сигнала в определенный момент времени с помощью вышеприведенного выражения необходимо обработать бесконечное число его отсчетов (как в прошлом, так и в будущем), что нереализуемо на практике. В реальной жизни интерполяция осуществляется с помощью других фильтров, при этом выражение для нее принимает следующий вид:
() ( ) ( ),
txthtt
ii
i
где h(t) — импульсная характеристика соответствующего восста­навливающего фильтра. Вид этого фильтра выбирается в зависимо­сти от задачи.
39
Прямое вычисление новых отсчетов сигнала по вышеприведен­ным формулам требует значительных вычислительных ресурсов и нежелательно для приложений реального времени. Существуют важные частные случаи передискретизации, для которых вычис­ление новых отсчетов производится проще:
децимация с целым коэффициентом (уменьшение частоты
дискретизации в целое число раз);
интерполяция с целым коэффициентом (увеличение частоты
дискретизации в целое число раз);
изменение частоты дискретизации в рациональное M/N чис-
ло раз (этот случай можно рассматривать как комбинацию двух предыдущих).
При таких ограничениях становится удобным применение стандартных реализаций цифровых фильтров для передискретиза­ции.
5.1. Выбор фильтра для передискретизации
Выбор функции h(t) обуславливается компромиссом между качеством передискретизации (т. е. близости ее к идеальной) и вычислительной сложностью этого процесса. В принципе, для передискретизации может быть использован любой фильтр ниж­них частот с необходимой частотой среза. КИХ-фильтры приме­няются для этих задач чаще, чем БИХ-фильтры, из-за возможно­сти построения КИХ-фильтров с линейной фазочастотной харак­теристикой.
Чаще всего при передискретизации используются следующие классы цифровых фильтров.
1. Фильтры, построенные исходя из критерия близости частот­ной характеристики к частотной характеристике идеального филь­тра нижних частот:
1.1) оконные sinc-фильтры — их импульсная характеристика h(t) получается путем умножения импульсной характеристики идеального ФНЧ на оконную функцию;
1.2) равноволновые фильтры Чебышева.
2. Классические способы интерполяции функций (часто при­меняются для изображений):
2.1) линейные интерполяторы;
40
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]