Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

CALCULATION OF COMBUNATIONAL FREQUENCIES FOR FILTRATION OF BIOELECTRICAL POTENTIALS BY PHASE FILTER”

Avdeeva D.K. 1 Vylegzhanin O.N. 1 Penkov P.G. 1 Nebaba S.G. 1
1 National Research Tomsk Polytechnic University
The features of numerical realization of filter for registration of bioelectrical potentials were considered. Mentioned filter is based on the using of procedure of phase modulation of harmonic of special form impulse by harmonics of measurable signal with noise. It was shown that in order to receive the maximal value of modulation parameter, it is needed to ensure the hit of initial phase values of reference impulse as well as the impulse with measurable signal on the ends of quasy-linear part of relationship between phase value and impulse echelons ratio. The algorithm for calculation of reference impulse’s parameters and of multiplication factor for measurable signal implemented this condition was received. The methods for selection of combinational frequencies with amplitudes above preset threshold were suggested. The model calculation for assigned signal with random noise were carry out.
calculation of combinational frequencies.
phase modulation
filtration
measurement of bioelectrical potentials
При медицинских исследованиях биоэлектрических потенциалов возникает задача выделения изучаемых сигналов из регистрируемой смеси собственного потенциала и шумов электрода и регистрирующей системы. Учитывая случайный характер шумовой компоненты и отсутствие априорной информации о ее характеристиках, выделение информативной составляющей является достаточно сложной задачей.  Ранее [1] нами был предложен метод выделения когерентного сигнала на фоне случайной помехи, основанный на применении процедуры фазовой модуляции.

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

Как было показано в [2], подобный фильтр может быть построен с использованием опорного импульса специального вида, показанного на рис.1.

pic 

Рис. 1. Форма опорного импульса

Импульс имеет две прямоугольные ступени с амплитудами U1 и U2 соответственно, общая длительность ступеней немного меньше длительности интервала наблюдения ТИ.

Как показали результаты анализа, при соотношениях амплитуд ступеней сигнала и отношении длительности интервала измерения к общей длительности ступеней, близких к 1, увеличивается чувствительность фазы к изменению напряжения (амплитуды ступени U2). При этом зависимость фазы от D (отношение амплитуды второй ступени к амплитуде первой) M (отношение длительности интервала измерения к суммарной длительности ступеней ) имеет вид:

f

Если поместить на интервале от ТИ/2 до ТИ опорного импульса измеряемый сигнал с помехой, то любая его гармоническая составляющая wi будет модулировать по амплитуде и фазе гармонику опорного импульса в виде [3,4,5]:

f

 

, (1)

Где β- глубины модуляции ( βi0 - φi, , т.е. разность значений фазы опорного импульса и модулирующей гармоники), Ω­частоты, а φ­ начальные фазы модулирующих гармоник. Из выражения (1) следует, что модулируемое колебание содержит целый набор комбинационных частот, отстоящих от основной частоты выше и ниже на величину f.

Из приведенных рассуждений вытекает следующая схема применения фильтра. Формируется опорный импульс, на нижнюю ступень которого помещается фрагмент аддитивной суммы измеряемого сигнала и шума. Выполнив преобразование Фурье для этого фрагмента, можно представить его в виде линейной комбинации гармоник разложения Фурье, которые и будут модулировать опорный импульс. Важным параметром, определяющим эффективность описываемого фильтра, является глубина модуляции. Анализ показывает, что максимальная глубина фазовой модуляции может быть равна p (гармоники находятся в противофазе). Но в реальных расчетах рекомендуется обеспечивать попадание фазы опорного импульса и фазы импульса с измеряемым фрагментом сигнала на квазилинейный участок графика зависимости фазы от параметра φ(D) .  Для этого был предложен следующий алгоритм - выбираем значение M, первую ступень опорного импульса принимаем равной 1, а затем вычисляем значение D , которое обеспечивает попадание начальной фазы основной гармоники на нижний конец квазилинейного участка кривой φ(D). Для этого вычисляется максимальное значение второй производной φ(D) по D, которое соответствует значению D=0. Тогда квазилинейный участок кривой  соответствует значениям D≅0.10 от этого значения. Сформировав таким образом опорный импульс, вычислим среднее значение измеряемого фрагмента сигнала  и умножим этот фрагмент на коэффициент f, где f - среднее значение фрагмента сигнал плюс шум. Такой выбор размеров опорного импульса и измеряемого фрагмента сигнала обеспечивает значение разности начальных фаз основной гармоники равным около 2,5 - 2,9 радиана.

Из выражения (1) следует, что для k  модулирующих частот, имеющих по p порядков функций Бесселя, получится (2 p+1)k комбинационных частот.  Таким образом, даже для процессора с тактовой частотой 2 ГГц вычисление  амплитуд и частот комбинационных гармоник становится довольно сложной задачей. Хотя в выражении (1) оператор суммирования формально выполняется для индексов от минус бесконечности до плюс бесконечности, поскольку по определению β≤π, то в реальных вычислениях, значения функции Бесселя с увеличением порядка быстро убывают. Поэтому, учитывая только слагаемые по модулю больше 0.01, для суммирования достаточно пяти членов с порядками -4, -2,0,2,4.

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

Исходный опорный импульс имел следующие параметры: количество точек отсчета 64, количество точек отсчета с нулевой амплитудой 1, высота первой ступени 1, D= 0.8516. В качестве модельного сигнала использовался полином второго порядка с коэффициентами 1, 1, -1. К этому сигналу добавлялся случайный, нормально распределенный шум с нулевым математическим ожиданием и дисперсией 0,2. Спектр шума на интервале существования, то есть в течение длительности второго импульса, равной Ти/2, представлен на графике (рис.2). Полученный в результате импульс представлен на рис. 3.

pic 

Рис. 2. Спектр зашумленного сигнала

 pic   pic

 а)        б)

Рис. 3. Графики: а) опорного импульса с добавленным фрагментом сигнала с шумом в четвертом квадранте показана исходная вторая ступень опорного импульса, б) полезного сигнала (серая линия) и добавленного шума (красная линия).

 Вычисленная начальная фаза основной гармоники опорного импульса равнялась -1.22955727604695, а начальная фаза этой же гармоники для опорного импульса с фрагментом сигнала без шума равнялась 1.3131183875097. Таким образом, разность фаз составила -2.54267566355665.  При добавлении шума начальная фаза основной гармоники становится равной 1.35917893369518. На рис.4 показан график первой производной с нанесенными точками, соответствующими: слева опорному импульсу, справа - импульсу с добавленным сигналом. Как видно из рисунка, точки расположены максимально близко к концам квазилинейного участка зависимости фазы от параметра  D.  Рассчитанные по соотношению (1) комбинационные частоты и амплитуды приведены на рис. 4.

pic 

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

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

Таблица 1. Значения фаз и частот модулирующих гармоник

Фаза радиан

Частота w

1.22955727604695

6.28318530717959

2.34750653814966

12.5663706143592

1.5423478107961

25.1327412287183

0.293137066358017

50.2654824574367

2.80177884625832

62.8318530717959

0.27693940127947

75.398223686155

Таблица 2. Рассчитанные комбинационные частоты и амплитуды

Амплитуда усл. ед.

Частота w.

0.02045246

-124.66

0.02009389

-99.53

0.02045246

-74.40

-0.01007805

-61.83

0.02045246

-49.26

0.02009389

-24.13

-0.01007805

-11.56

-0.01004917

1.0

-0.01007805

13.57

0.020093895

26.13

0.020093895

51.26

-0.01007805

63.83

0.02045246

76.40

0.020093895

101.53

0.02045246

126.66

На рисунке 5 представлен результат применения описанного фильтра к зарегистрированному кардиоимпульсу. Как видно из рисунка, примененная фильтрация позволила заметно улучшить качество записи.

pic 

А)

pic

Б)

Рис.5. Кардиоимпульс а) зарегистрированный; б) после обработки описанным алгоритмом

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