Внедрение фильтров Iir в портативные медицинские устройства для точной фильтрации сигналов
Внедрение фильтров IIR в портативные медицинские устройства для точной фильтрации сигналов
В быстро развивающейся области портативных медицинских устройств точная фильтрация сигналов необходима для надежной диагностики и непрерывного мониторинга. Портативные устройства, такие как портативные мониторы ЭКГ, носимые гарнитуры ЭЭГ и беспроводные датчики ЭМГ, должны извлекать чистые биомедицинские сигналы из шумной среды при работе в условиях строгих ограничений мощности и размера. Фильтры с бесконечным импульсным откликом (IIR) стали краеугольным камнем в этих системах из-за их вычислительной эффективности и способности обеспечивать резкую частотную селективность с относительно небольшими коэффициентами. Эта статья предоставляет всеобъемлющий технический обзор того, как фильтры IIR внедряются в портативные медицинские устройства, охватывающие основы проектирования, соображения стабильности, проблемы обработки в реальном времени, оптимизацию мощности и будущие направления.
Основы фильтров IIR
Фильтры IIR (Infinite Impulse Response) — это цифровые фильтры, которые используют обратную связь от прошлых выходов, а также прошлых входов для формирования частотного ответа. Их функция передачи является рациональным полиномом в Z-домене:
H(z) = (b0 + b1z-1 + b2z-2 + ... + bmz-m) / (1 + a1z-1 + a2z-2 + ... + anz-n)
Рекурсивный характер знаменателя дает IIR-фильтрам бесконечную импульсную реакцию. Эта обратная связь позволяет фильтру достичь заданной частотной спецификации с гораздо меньшим по сравнению с фильтрами Finite Impulse Response (FIR) порядком. Для портативных медицинских устройств это напрямую приводит к меньшему умножению и добавлению на образец, что означает более низкое энергопотребление и снижение использования памяти.
Разница между фильтрами IIR и FIR
Выбор между IIR и FIR фильтрами является классическим компромиссом в цифровой обработке сигналов. FIR фильтры по своей сути стабильны, имеют линейную фазу (важно для сохранения морфологии формы волны в ЭКГ) и просты в проектировании с помощью простых методов окон. Однако они требуют высокого порядка фильтрации для достижения резких отсечек, потребляя больше энергии и памяти. IIR фильтры, с другой стороны, достигают аналогичной частотной производительности с 5-10-кратным меньшим коэффициентом, но они вводят фазовое искажение и могут стать нестабильными, если не тщательно разработаны. Для многих портативных медицинских приложений, где мощность в реальном времени является критической и небольшое фазовое искажение приемлемо (например, мониторинг частоты сердечных сокращений, удаление блуждания по базовой линии), IIR фильтры являются предпочтительным выбором.
Типы фильтров IIR, используемые в медицинских устройствах
Несколько классических аналоговых фильтров-прототипов преобразуются в цифровые IIR-фильтры с использованием билинейного преобразования или импульсной инвариантности. Наиболее распространенными типами, используемыми в портативных медицинских устройствах, являются:
Фильтры Butterworth
Фильтры Баттерворта обеспечивают максимально плоскую полосу пропускания с плавным откатом. Они идеально подходят для приложений, где важно сохранение амплитуды сигнала по полосе пропускания, например, в кондиционировании сигнала ЭКГ. Фильтр Баттерворта низкого прохода порядка 4 или 6 может эффективно удалять высокочастотный мышечный шум и электромагнитные помехи при сохранении сложной морфологии QRS.
Чебышев Фильтры
Фильтры типа I Чебышева имеют более крутой откат, чем Butterworth за счет полосы пропускания. Они полезны, когда требуется резкое отсечение, и небольшое количество полосы пропускания допустимо. В приложениях ЭЭГ требования к затуханию полосы остановки высоки, чтобы отклонить шум линий электропередач на 50/60 Гц, не затрагивая близлежащие альфа- или бета-диапазоны. Фильтры Чебышева могут соответствовать таким спецификациям с более низким порядком, чем Butterworth.
Эллиптические (Cauer) фильтры
Эллиптические фильтры предлагают самую острую полосу перехода для данного порядка, позволяя рябь как в полосе пропускания, так и в полосе остановки. Они очень эффективны в вычислительном отношении, но имеют наибольшее фазовое искажение. Они иногда используются в имплантируемых устройствах, где бюджет мощности чрезвычайно ограничен, а интересующий сигнал (например, ускоренный сердечный ритм) не чувствителен к нелинейности фазы.
Разработка фильтров IIR для биомедицинских сигналов
Процесс проектирования IIR-фильтров в портативных медицинских устройствах должен учитывать специфические характеристики целевого биосигнала и ограничения встроенного оборудования.
Спецификации фильтров для общих сигналов
- ECG (электрокардиограмма): Частотный диапазон 0,05—100 Гц. Типичные фильтры включают в себя высокий проход на 0,5 Гц для удаления исходного блуждания и низкий проход на 40—100 Гц для отбрасывания шума мышц и помех силовой линии. Фильтр выемки на 50/60 Гц часто реализуется в виде биквадратичной секции IIR.
- ЭЭГ (электроэнцефалограмма): Частоты от 0,5 до 60 Гц (дельта, тета, альфа, бета, гамма). Для полосовой фильтрации используются фильтры IIR для выделения конкретной полосы, например, альфа-полосного фильтра (8-12 Гц) с резким скатыванием, чтобы избежать перекрестных стыковок с другими полосами.
- ЭМГ (электромиография): Сигналы от 10 до 500 Гц. Фильтр IIR полосового прохода удаляет низкочастотные артефакты движения и высокочастотный шум. Поскольку сигналы EMG нестационарны, иногда используются адаптивные фильтры IIR.
- Оксиметрия импульса (PPG): Частотный диапазон 0,5—10 Гц. Фильтры ИИР низкого пропуска используются для извлечения фотоплетизмографической формы волны и удаления окружающего светового шума.
Билайнарная трансформация и преварпинг
Наиболее распространенным методом преобразования аналоговых фильтров-прототипов в цифровые IIR-фильтры является билинейное преобразование. Поскольку билинейное преобразование сжимает частотную ось, происходит искажение частоты. Для компенсации критические частоты (вырез, выемка) должны быть предварительно деформированы. Формула предварительного деформирования:
ωa = (2/T) * tan(ωd / 2)
где ω d — желаемая цифровая частота, а T — период отбора проб. Неспособность к пре-варп приводит к фильтрам с неверными частотами отсечения, что может ухудшить диагностическую точность.
Стабильность и численная точность в фильтрах IIR
Стабильность является первостепенной проблемой для IIR-фильтров в любом медицинском устройстве в реальном времени. Фильтр, который становится неустойчивым из-за квантования коэффициента или переполнения, может производить опасные выходные артефакты, которые могут быть неправильно истолкованы как клинические события.
Позиции полюсов на z-плане
Дискретный IIR-фильтр стабилен тогда и только тогда, когда все полюса лежат внутри единичного круга в z-плане. Полюсы являются корнями многочлена знаменателя. По мере увеличения порядка фильтра полюса группируются вблизи единичного круга для резких отсечек, делая стабильность более чувствительной к квантованию коэффициента. Конструкторы должны обеспечить достаточный запас (например, полюса внутри радиуса 0,95) для учета численных ошибок.
Проблемы реализации в фиксированной точке
Портативные медицинские устройства часто используют микроконтроллеры с фиксированной точкой без блоков с плавающей точкой для экономии затрат и мощности. Преобразование коэффициентов фильтра IIR и выполнение арифметики с конечной длиной слова может привести к:
- Ошибка квантования коэффициентов: Коэффициенты округлены до ближайшего репрезентабельного значения, что может отклонить фактический частотный отклик от цели проектирования.Использование каскадных секций второго порядка (SOS) вместо реализации прямой формы снижает чувствительность к квантованию коэффициентов.
- Квантовый шум: Округление или усечение промежуточных результатов вносит шум. Это можно смягчить с помощью более высокоточных аккумуляторов (например, 32-битных или 64-битных) и затем округлением обратно к входной длине слова.
- Переливные колебания: Если уровни сигнала не масштабируются должным образом, перелив может вызвать предельные циклы или большие колебания.
Каскадные секции второго порядка (SOS)
Стандартная методика внедрения IIR-фильтров во встроенные системы заключается в разложении фильтра высокого порядка в каскад биквадратических секций второго порядка. Каждая секция биквада реализует двухполюсный, двухнулевой фильтр. Преимущества включают:
- Улучшенная стабильность чисел: Каждая биквада имеет полюса, которые легче держать внутри единичного круга.
- Легче масштабировать : Промежуточные уровни могут быть скорректированы для предотвращения переполнения.
- Меньшая чувствительность к квантованию коэффициентов по сравнению с структурами прямой формы.
Для портативных медицинских устройств, секции биквада обычно реализуются с использованием транспонированной структуры прямой формы II, которая является более эффективной и менее подвержена переполнению, чем прямая форма I.
Оптимизация мощности для устройств с батарейным питанием
Потребление энергии является ключевым ограничением в портативных медицинских устройствах. Фильтры IIR способствуют экономии электроэнергии как за счет алгоритмической эффективности, так и за счет методов реализации.
Алгоритмическая эффективность
Поскольку фильтры IIR требуют меньше операций многократного накопления на образец, чем фильтры FIR для той же спецификации частоты, они уменьшают количество циклов процессора на образец. Для типичной выборки монитора ЭКГ на 500 Гц замена фильтра нижних частот FIR с пропуском 100 касаний на фильтр IIR 6-го порядка снижает арифметические операции более чем на 90%. Это позволяет микроконтроллеру проводить больше времени в режимах сна с низким энергопотреблением.
Низкоэнергетические стратегии аппаратного обеспечения
- Прямой доступ к памяти (DMA) на основе выборки: ADC может заполнить буфер через DMA без вмешательства ЦП. ЦП просыпается только тогда, когда буфер заполнен для запуска процедуры фильтра IIR, а затем возвращается в сон.
- Настраиваемое оборудование с фиксированной точкой : Некоторые микроконтроллеры включают аппаратный мультипликатор-аккумулятор (MAC), который может выполнять фильтрацию биквада в одном цикле команд, резко снижая энергию за операцию.
- Адаптивная частота отбора проб : Для сигналов, таких как частота сердечных сокращений, после получения хорошей оценки частота отбора проб может быть уменьшена. Коэффициенты фильтра IIR должны быть пересчитаны для новой скорости или используется конструкция фильтра с переменной скоростью.
Ограничения обработки в реальном времени
Портативные медицинские устройства должны обрабатывать сигналы в режиме реального времени для обеспечения немедленной обратной связи. Например, амбулаторный монитор ЭКГ должен обнаруживать аритмии в течение нескольких секунд. Фильтры IIR с их рекурсивной природой могут вводить групповую задержку. Конструкторы должны балансировать порядок фильтра (что увеличивает задержку) против требований времени отклика. Фильтр низкого пропуска 4-го порядка Баттерворта на частоте 40 Гц со скоростью выборки 250 Гц вводит групповую задержку около 8 образцов (32 мс), что приемлемо для большинства приложений мониторинга. Для приложений, требующих минимальной задержки (например, обнаружения захвата в режиме реального времени в ЭЭГ), могут использоваться фильтры более низкого порядка или предварительно вычисленные реакции шага.
Тематические исследования: фильтры IIR в конкретных портативных медицинских устройствах
Портативные мониторы ЭКГ
Современные портативные устройства ЭКГ (например, AliveCor, Biomonitor) в значительной степени полагаются на IIR-фильтры для кондиционирования сигнала.
- Фильтр высокого разрешения 0.5 Гц по методу Баттерворта порядка 2 для удаления базового дрейфа без удаления клинически значимой информации ST-сегмента.
- 40-Гц фильтр Чебышева с низким пропуском (заказ 4, 0,5 дБ ряби) для подавления мышечного шума при сохранении амплитуды QRS.
- Фильтр выемки 50/60 Гц реализован в виде биквада (типа Twin-T) для отбрасывания помех линии электропередач. Фильтр выемки часто имеет полосу пропускания 1-2 Гц, чтобы избежать искажения волн P и T.
Все три фильтра каскадируются как SOS-разделы на маломощном процессоре ARM Cortex-M0+, работающем на частоте 16 МГц. Общее время выполнения фильтра на образец составляет менее 20 циклов ЦП, оставляя достаточно времени для обнаружения QRS и передачи данных.
Носимые гарнитуры ЭЭГ
Носимые ЭЭГ-устройства (например, Emotiv, NeuroSky) используют IIR-фильтры для выделения конкретных полос мозговых волн в реальном времени. Например, конструкция с двумя фильтрами: сначала фильтр полосы частот 0,5-50 Гц Butterworth (порядок 4) удаляет смещение постоянного тока и высокочастотный шум. Затем отдельные фильтры полосы IIR извлекают дельта (0,5-4 Гц), тета (4-8 Гц), альфа (8-13 Гц), бета (13-30 Гц) и гамма (30-50 Гц). Каждый полосовой проход использует 4-й фильтр Чебышева типа II для остановки утечки из полосы. Фильтры работают на ядре DSP на 100 МГц, используя менее 5% нагрузки ЦП, что позволяет в реальном времени использовать интерфейс мозг-компьютер (BCI) приложения.
Беспроводные датчики EMG
Беспроводные поверхностные датчики EMG для протезирования и реабилитации часто используют адаптивные IIR-фильтры для удаления артефактов движения. Адаптивный IIR-фильтр наименьшего квадрата (LMS) настраивает выемку на доминирующей частоте артефакта движения, которая может варьироваться при движении конечностей. Адаптивная структура IIR использует форму решетки второго порядка для стабильности и низкой сложности. Реализация на маломощном Bluetooth SoC позволяет работать от батареи 8 до 12 часов при сохранении высокого качества сигнала.
Сравнение с фильтрами FIR: когда фильтры IIR лучше выбрать?
Хотя фильтры IIR популярны, они не всегда являются оптимальным выбором. В таблице ниже приведены ключевые компромиссы, имеющие отношение к портативным устройствам:
| Feature | IIR | FIR |
|---|---|---|
| Computational cost (for same spec) | Low | High |
| Memory usage | Low | High |
| Power consumption | Low | Higher |
| Phase linearity | Non‑linear | Linear possible |
| Stability | Conditional | Guaranteed |
| Design flexibility | Analog prototypes | Window / Parks‑McClellan |
Для приложений, где фазовая линейность имеет решающее значение (например, сохранение времени R-пика в ЭКГ для точного обнаружения QRS), FIR-фильтры могут быть предпочтительными, несмотря на более высокую мощность. Однако многие портативные устройства используют комбинацию: фильтр IIR низкого порядка для грубого снижения шума, за которым следует более короткий фильтр FIR для компенсации, если это необходимо. Современные системы с низким энергопотреблением часто реализуют линейно-фазовый IIR-фильтр с использованием метода фильтрации нулевой фазы (фильтрация вперед-назад), но который вводит задержку, которая может быть неприемлемой в реальном времени. В конечном счете, выбор зависит от конкретного сигнала, бюджета мощности и клинического требования.
Тестирование и валидация фильтров IIR в медицинских устройствах
Стандарты регулирования (IEC 60601, ISO 13485) требуют, чтобы цифровые фильтры в медицинских устройствах были проверены на безопасность и производительность.
- Измерение частотного отклика: Сравните фактический отклик фильтра (включая эффекты квантования коэффициентов) с техническими характеристиками конструкции.
- Стабильность по температуре и полный диапазон входных данных : моделирование экстремальных условий входных данных (например, больших всплесков артефактов) для обеспечения отсутствия нестабильности или предельных циклов.
- Анализ шума: Измерение шумового пола и соотношения сигнал/шум после фильтрации.
- Входные данные крайних случаев : Тестирование синусоидальными входами на критических частотах (вырез, выемка) для проверки производительности отторжения.
Дополнительное тестирование для реализации с фиксированной точкой включает в себя тесты переполнения, проверку масштабирования и моделирование бит-правды, сравнивающее ссылку с плавающей точкой с переносным кодом с фиксированной точкой.
Будущие направления
Эволюция портативных медицинских устройств продолжает расширять границы проектирования фильтров ИКД. К ключевым тенденциям относятся:
Адаптивные фильтры IIR
Адаптивные IIR-фильтры могут корректировать свои коэффициенты в реальном времени для отслеживания изменяющихся характеристик сигнала. Например, адаптивный выемочный фильтр может динамически блокировать изменения частоты линий электропередач (50±1 Гц) и поддерживать высокий отказ. Для быстрой конвергенции используется алгоритм рекурсивных наименьших квадратов (RLS), но его вычислительная сложность высока. Упрощенные варианты, такие как утечка LMS IIR-фильтра, реализуются на маломощных ускорителях нейронной сети.
Интеграция машинного обучения
Гибридные системы объединяют предфильтры IIR с классификаторами машинного обучения (свёрточные нейронные сети, SVM), работающими на одном устройстве. Фильтр IIR снижает пространственную размерность функции, позволяя классификатору работать с меньшим количеством ресурсов. Например, портативный монитор ЭКГ использует фильтр высокого пропуска IIR 4-го порядка на 0,5 Гц, затем CNN обнаруживает аритмии. Общее энергопотребление ниже, чем чистая обработка CNN необработанных данных.
Программно-определяемая фильтрация с помощью реконфигурируемого оборудования
Полевые программируемые затворные массивы (FPGA) и грубозерные реконфигурируемые массивы (CGRA) становятся достаточно энергоэффективными для портативных устройств. Они позволяют перенастраивать фильтры IIR на лету, позволяя одному устройству переключаться между режимом ЭКГ (0,5-40 Гц) и режимом ЭЭГ (0,5-60 Гц) путем загрузки различных наборов коэффициентов. Эта гибкость ценна для носимых патчей с несколькими датчиками.
Заключение
Фильтры IIR остаются критическим фактором для портативных медицинских устройств, обеспечивая эффективное кондиционирование сигналов в реальном времени, которое отвечает ограничениям мощности, размера и производительности современных носимых медицинских устройств. Понимая компромиссы между типами фильтров, тщательно проектируя для цифровой стабильности и используя оборудование с низким энергопотреблением, инженеры могут создавать надежные решения для фильтрации, которые повышают диагностику и безопасность пациентов. По мере развития технологий адаптивные и управляемые ИИ реализации фильтров IIR будут еще больше расширять возможности портативных медицинских устройств, принося сложный мониторинг в точки ухода и домашнюю среду.
Для дальнейшего чтения обратитесь к следующим авторитетным источникам:
- Kenny, P. A. & Zhu, S. (2019). «Digital Filter Design for Biomedical Signal Processing» IEEE Signal Processing Magazine, 36(4), 88-96.IEEE Xplore»
- Томпкинс, У. Дж. (2000). Биомедицинская обработка цифровых сигналов . Прентис-Холл. (глава 7: Фильтры IIR) Пирсон
- Texas Instruments (2021). «Доклад о применении: цифровая фильтрация с микроконтроллерами MSP430». TI SLAA512
- Натараджан Р. (2020). «Проектирование и внедрение фильтра IIR Notch для удаления помех ЭКГ в силовых линиях». Журнал медицинских систем , 44(5), 92. Спрингер