Проектирование фильтров Iir для высокоточного оборудования для измерения и калибровки

Введение в фильтры IIR для точного измерения

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

В отличие от своих аналогов с конечным импульсным откликом (FIR), фильтры IIR используют обратную связь для достижения заданного затухания полосы остановки или ширины перехода с гораздо меньшим количеством коэффициентов. Эта эффективность уменьшает использование памяти и умножает меньше раз на образец, что имеет решающее значение во встроенных системах, которые обрабатывают несколько каналов одновременно. Однако эта эффективность достигается за счет: повышенной чувствительности к квантованию коэффициентов, потенциальной нестабильности и нелинейной фазовой реакции. Для калибровочных применений, где фазовая линейность или постоянная задержка группы могут быть важны, инженер должен тщательно оценить, какая топология IIR и метод проектирования наилучшим образом соответствуют ограничениям измерения.

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

Основы фильтров IIR

Фильтр IIR определяется его дифференциальным уравнением, которое включает в себя как условия подачи-передачи, так и условия обратной связи:

y[n] = b₀ x₁ x[n-1] + ... + b_M x _{— a1 — a2 — ... — aN y}

]

Коэффициенты обратной связи (a_]k) вводят полюса в передаточной функции фильтра, позволяя более крутое выкатывание на полюс по сравнению с фильтрами FIR. Ноли (из терминов b_]k) обеспечивают дополнительную гибкость для формирования частотной реакции. Поскольку выход зависит от предыдущих выходов, импульсная реакция простирается бесконечно — отсюда «Реакция бесконечного импульса».

Преимущества перед фильтрами FIR

Для высокоточных измерений основными преимуществами IIR-фильтров являются:

• Низкий порядок эквивалентной производительности: Фильтр IIR может достигать заданной полосы пропускания и затухания полосы остановки с порядком фильтра в два-пять раз ниже, чем фильтр FIR. Это снижает вычислительную нагрузку и задержку, оба из которых являются ценными в калибровочных петлях реального времени.
• Меньший след памяти: Меньшее количество коэффициентов означает меньшее количество памяти в DSP с фиксированной точкой или FPGA. В многоканальных системах это сохранение быстро умножается.
• Более узкие полосы перехода: Используя эллиптические или чебышевские методы, фильтры IIR могут достигать чрезвычайно резких отсечек, что полезно при разделении близко расположенных частотных компонентов в сигналах датчика (например, отторжение линии электропередачи 50/60 Гц при сохранении измерения низкой частоты).

Торговые компромиссы для решения

Недостатками необходимо управлять:

• Нелинейная фаза: Большинство фильтров IIR вводят фазовое искажение, которое изменяется с частотой. Для задач измерения, требующих сохраненной формы формы волны (например, рефлектометрия временной области или анализ импульсов удара), может потребоваться фильтр FIR или IIR Бесселя.
• Чувствительность к стабильности: Поляки вблизи единичного круга делают фильтр чувствительным к квантованию коэффициента.Незначительные ошибки округления могут вытолкнуть полюс за пределы единичного круга, вызывая колебания.
• Круглый шум: В реализациях с фиксированной точкой петля обратной связи может усиливать шум квантования, если структура фильтра не выбрана тщательно (например, с использованием прямой формы I против прямой формы II, перенесённой).

Ключевые соображения дизайна для высокоточного оборудования

Разработка фильтра ИКД для калибровочного оборудования выходит за рамки выбора частот отсечки и затухания стоп-диапазона. Следующие соображения часто являются разницей между фильтром, который работает в моделировании, и фильтром, который надежно работает на скамейке.

Стабильность и размещение полюсов

Стабильность не подлежит обсуждению. Все полюса должны находиться внутри единичного круга в z-плане. Для высокоточной работы инженеры обычно добавляют запас прочности — например, ограничивая радиус полюса до 0,95 или менее — для учета дрейфа температуры и округления коэффициента. При реализации фильтра высокого порядка стандартной практикой является разбиение функции передачи на секции второго порядка (SOS) и обеспечение того, чтобы каждая секция индивидуально соответствовала критерию стабильности.

Отзывы и групповая задержка

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

Затухание стоп-диапазона и динамический диапазон

Высокоточное измерительное оборудование часто должно отклонять шум, который находится на десятки или сотни дБ ниже амплитуды сигнала. Например, 24-битный АЦП имеет динамический диапазон, превышающий 120 дБ, но фильтр не должен ограничивать этот динамический диапазон. Затухание стоп-диапазона цифрового фильтра должно быть не менее 6 дБ ниже уровня шума АЦП, чтобы избежать ухудшения эффективного количества бит системы (ENoB). Это требование часто подталкивает конструкторов к эллиптическим или чебышевским фильтрам типа II, которые обеспечивают высокое затухание стоп-диапазона с минимальным порядком.

Вычислительная эффективность и пропускная способность в реальном времени

В системе калибровки, которая обрабатывает несколько датчиков с высокой частотой выборки (например, 1 МГц или более), каждое умножение имеет значение. Структуры фильтра IIR, такие как транспонированная прямая форма II и каскад SOS, уменьшают количество состояний и операций на образец. Кроме того, используя коэффициент мощности двух, где это возможно, можно заменить умножение с сдвигами, хотя это приносит гибкость.

Квантирующие эффекты

Квантирование коэффициентов и арифметическое округление являются основными источниками ухудшения производительности. Для высокоточных фильтров предпочтительнее использовать арифметику с плавающей точкой двойной точности всякий раз, когда аппаратное обеспечение поддерживает ее. Когда требуется фиксированная точка - из-за ограничений ресурсов FPGA или недорогих микроконтроллеров - инженер должен масштабировать коэффициенты и внутренние состояния, чтобы избежать переполнения при сохранении точности сигнала. Такие методы, как округление коэффициента до ближайшего целого числа, с использованием “ мощность двух & #8221; масштабирование и добавление битов защиты в аккумуляторе являются стандартными практиками.

Методы проектирования высокоточных IIR-фильтров

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

Фильтры Butterworth

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

Чебышев Тип I и Тип II Фильтры

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

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

Дизайн рабочего процесса с современными инструментами

Инженеры обычно разрабатывают IIR-фильтры с использованием библиотек программного обеспечения, таких как MATLAB Signal Processing Toolbox, Python’s SciPy (] или специализированные инструменты проектирования DSP.

1. Укажите полосы пропускания и стоп-диапазоны, допустимую рябь (пассажирскую полосу) и затухание (стоп-диапазон).
2. Выберите прототип фильтра на основе компромиссов.
3. Преобразовать аналоговый прототип в цифровой ИДК с использованием билинейного преобразования или инвариантности импульсов (билайнар предпочтительнее из-за избегания псевдонимов).
4. Преобразовать полученную функцию передачи в формат SOS для реализации.
5. Квантировать коэффициенты к целевой арифметике (например, 16-бит, 24-бит или двойной) и имитировать частотную реакцию квантованного фильтра.
6. Проверяйте дизайн с реальными шумовыми и сигнальными условиями.

Практические проблемы реализации

Даже идеально спроектированный фильтр в двухточной плавающей точке может демонстрировать ухудшенную производительность при развертывании на аппаратных средствах с фиксированной точкой.

Ошибки квантирования коэффициентов

Когда коэффициенты округлены до ближайшего представления с фиксированной точкой, места полюсов могут смещаться. Для фильтра с низким пропуском высокого Q (например, узкополосной выемки) небольшая ошибка в радиусе полюса может привести к пересечению полюса единичным кругом, что приводит к колебаниям. Методы для смягчения этого:

• Используйте каскад секций второго порядка (SOS) для изоляции чувствительных полюсов.
• Конструкция фильтра с избыточной точностью, а затем круглые коэффициенты таким образом, чтобы гарантировать стабильность (например, с использованием ограничений полюс-радиус).
• Внедрить фильтр в одноточную плавающую точку на микроконтроллере с аппаратным FPU (например, ARM Cortex-M4/M7).

Перелив и масштабирование

В арифметике с фиксированной точкой внутренние состояния фильтра могут переполняться, если усиление фильтра не управляется должным образом. Следует вычислить пиковое значение импульсной реакции фильтра и применить коэффициенты масштабирования к входам или коэффициентам, чтобы гарантировать, что все внутренние узлы остаются в пределах диапазона. Многие процессоры DSP обеспечивают логику насыщения, но, опираясь на нее, могут вводить искажение вырезки. Лучший подход - предварительно масштабировать усиление фильтра до единства в полосе пропускания.

Ограниченные циклы и мертвые зоны

Фильтры ИИР, реализованные в арифметике с фиксированной точкой, могут демонстрировать предельные циклы — устойчивые низкоуровневые колебания, вызванные ошибками округления в пути обратной связи. Это особенно проблематично при точном измерении, потому что колебания появляются в виде ложных тонов на выходе. Чтобы минимизировать предельные циклы:

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

Проверка и испытание

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

• Частотная проверка: Используйте синус с прометанием, чтобы проверить, соответствует ли отклик на величину спецификации дизайна в допустимых пределах.
• Шаговые ответы: Оценка превышения и времени оседания, которые имеют решающее значение для измерений, которые быстро меняются.
• Измерение шумового пола: При входе, заземленном или подключенном к точной ссылке, измеряйте выходную шумовую спектральную плотность. Любые предельные циклы или отроги будут отображаться как пики.
• Долгосрочная стабильность: Запуск фильтра в течение нескольких часов или дней при мониторинге дрейфа или нестабильности из-за округления коэффициента с течением времени.

Пример: разработка фильтра Notch для отказа от линии электропередач

Общим требованием в высокоточных измерениях является отказ от фундаментальных и гармонических линий электропередачи (например, 50/60 Гц и кратных). Узкий полосовой (выемочный) фильтр должен иметь глубокий нуль для подавления помех без искажения близлежащих частот. Эллиптический IIR-фильтр порядка 6 (три SOS) может достигать затухания 60 дБ при 60 Гц с переходной полосой всего 2 Гц с каждой стороны. Процедура проектирования в Python:

import scipy.signal as signal
fs = 1000.0 # sampling frequency
f0 = 60.0 # notch frequency
bw = 2.0 # bandwidth
sos = signal.iirnotch(f0, bw/f0, fs)

Это дает секции SOS, которые могут быть реализованы непосредственно. Для систем с фиксированной точкой коэффициенты от билинейного преобразования должны быть тщательно масштабированы и округлены. Полученный фильтр удаляет тон 60 Гц, оставляя измерение постоянного тока практически незатронутым, демонстрируя мощность конструкций IIR в задачах калибровки.

Заключение

Проектирование IIR-фильтров для высокоточного измерительного и калибровочного оборудования требует глубокого понимания как теоретических компромиссов, так и практических ограничений реализации. Инженер должен сбалансировать затухание стоп-диапазона, линейность фазы, стабильность и вычислительную нагрузку, одновременно учитывая неизбежные эффекты квантования коэффициента и конечной длины слова. Выбрав соответствующий прототип - будь то Баттерворт для плоских полос пропускания, Чебышев для крутого скатывания, эллиптика для максимального отклонения или Бессель для точности фазы - и используя каскадные структуры SOS с тщательным масштабированием, можно построить фильтры, которые повышают точность критических систем измерения.

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