Обработка электрического сигнала с помощью фильтрации (стр. 1 из 2)

ПЕНЗЕНСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

Кафедра «Медицинские приборы и оборудование»

Курсовой проект

на тему:

«Обработка электрического сигнала с помощью фильтрации»

по дисциплине «Методы обработки биомедицинских сигналов»

Пенза 2005

Введение

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

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

1. Основы анализа электрических сигналов

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

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

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

2. Исследование полосового фильтра

Полосовой (или полосно-пропускающий) фильтр представляет собой устройство, которое пропускает сигналы в диапазоне частот с шириной полосы BW, расположенной приблизительно вокруг центральной частоты ω₀ (рад/с). На рисунке 1 изображена идеальная и реальная амплитудно-частотные характеристики полосового фильтра.

Рисунок 1

В реальной характеристике частоты ω_L и ω_U представляют собой нижнюю и верхнюю частоты среза и определяют полосу пропускания ω_L≤ω₀≤ω_U и ее ширину BW= ω_U – ω_L.

В полосе пропускания амплитудно-частотная характеристика никогда не превышает некоторого определенного значения, например А₁. Существует также две полосы задерживания, где значение амплитудно-частотной характеристики никогда не превышает заранее выбранного значения, наример А₂. Диапазоны частот пежду полосами задерживания и полосой пропускания образуют нижнюю и верхнюю переходные области, в которых характеристика является монотонной.

Отношение Q= ω₀/BW характеризует качество самого фильтра и является мерой его избирательности. Высокому значению Q соответствует относительно узкая, а низкому значению Q – относительно широкая ширина полосы пропускания. Коэффициент усиления фильтра К определяется как значение его амплитудно-частотной характеристики на центральной частоте.

Передаточная функция полосно-пропускающего фильтра имеет следующий вид:

Для исследования полосового фильтра, ввели в командное окно пакета MATLAB оператор «rlcdemo». Открылось окно, предназначенное для построения аналогового фильтра. Выбрали полосовой фильтр с последовательным соединением соединением элементов RLC: R=4,5; L=2,5; C=0,5 (рисунок 2).

Рисунок 2

Исходя из значений системных параметров рассчитали коэффициенты полиномов числителя и знаменателя:

G(s)= 1.8s

S²+1.8s+0.8

Средствами MATLAB построли прередаточную функцию системы sys:

>> sys=tf([0 1.8 0], [1 1.8 0.8])

Transfer function:

1.8 s

–

s^2 + 1.8 s + 0.8

Затем мы исследовали характеристики фильтра.

Во-первых, это диаграмма Боде (ЛАЧХ И ФЧХ) (рисунок 3).

>> bode(sys)

Рисунок 3

Во-вторых, это реакция системы на единичное воздействие (рисунок 4).

>> step(sys)

Рисунок 4

В-третьих, это диаграмма Найквиста, то есть представление вещественной и мнимой частей характеристики в зависимости от частоты (рисунок 5).

>> nyquist(sys)

Рисунок 5

И наконец, это импульсная характеристика (рисунок 6).

>> impulse(sys)

Рисунок 6

Сравнивая полученные характеристики фильтра с представленными на рисунке 2, убедились, что они идентичны.

3. Реализация модели фильтра в Simulink

Перевели данные содержащие электрокардиограмму (приложение 1) из MSExcel в MATLAB и построили график кардиосигнала (рисунок 7).

>> plot(G)

Рисунок 7

Затем наложили на исходный сигнал помеху (рисунок 8):

n=length(G);

e=rand(n)*17225

for i=1:n

s(i)=e(i)+G(i)

end

>> plot(s)

Рисунок 8

r=0

for i=1:n;

r=r+1;

t(i)=r;

end

t1=t'

S=s’

В приложении Simulink собрали схему, реализующую модель данного фильтра (рисунок 9).

Рисунок 9

На вход фильтра подается сигнал с помехой. Значения коэффициентов полиномов передаточной функции взяты из рабочей области MATLAB. На рисунке 10 представленны результаты фильтрации и исходный сигнал с помехой.

Рисунок 10

Заметим, что результат фильтрации неудовлетворителен, появляется значительный по амплитуде отрицательный выброс, в то время как амплитуда R-зубца уменьшается почти в четыре раза.

4. Реализация полосового фильтра в цифровой форме

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

В MATLAB существует библиотека (toolbox) Filter Design, предназначенная для решения задач проектирования, включая процедуры коррекции полученных результатов с учётом эффектов квантования. Библиотека имеет интерактивную графическую программу (GUI) fdatool, с помощью которой можно рассчитывать цифровые фильтры.

Сразу после загрузки программы отображается главное окно библиотеки. В нем мы задали параметры проектируемого фильтра (рисунок 11) в соотфетствии с полученными ранее ЛАЧХ И ФЧХ.

Рисунок 11

Воспользовавшись кнопками Import filter, Filter Coefficient, вывели значения коэффициентов фильтра (рисунок 12).

Рисунок 12

Последовательно выбирая команды Magnitude Recponse, Phase Recponse, Impulse Recponse, Step Recponse, определили характеристики фильтра.

Амплитудно-частотная характеристика фильтра (рисунок 13)

Рисунок 13

Фазо-частотная характеристика (рисунок 14)

Рисунок 14

Импульсная характеристика (рисунок 15)

Рисунок 15

Отклик на единичное воздействие (рисунок 16)

Рисунок 16

5. Реализация модели фильтра в Simulink

Нажав на кнопку

мы получаем созданный фильтр в виде блока-элемента библиотеки Simulink. Дважды щелкнув мышью на изображении этого блока можно получить его внутреннюю структуру (рисунок 17)

Рисунок 17

На рисунке 18 представлена схема, реализующая модельданного фильтра.

Рисунок 19

На вход фильтра подавали тот же сигнал, что и при исследовании аналогового фильтра

На рисунке 20 представлены результаты фильтрации и исходный сигнал с помехой.