Среди адаптивных методов приближения сигнала наибольший практический интерес представляют апертурные методы, осуществляющие контроль абсолютно ошибки при определении избыточных отсчетов и выборе существенных, т. е. передаваемых ординат. Они нашли широкое применение в задачах оперативной передачи и обработки сигнала из-за высокого быстродействия и простоты реализации.
Принцип их действия заключается в последовательном продвижении по дискретным регулярным отсчетам, полученным после дискретизации непрерывного сигнала, до некоторого п-го отсчета, в котором отклонение аппроксимированной ординаты от исходной превышает некоторое значение, задаваемое апертурой d. Ордината, первой вышедшая за пределы коридора шириной d, принимается за условную существенную ординату. Кроме того, вводится понятие существенной ординаты, используемой для передачи, обработки или восстановления сигнала. Выбор существенной ординаты зависит от конкретной реализации алгоритма. Во всех алгоритмах используется апертура, фиксированная по величине (d=2e, где e-максимально допустимое отклонение) и центрированная относительно аппроксимирующей прямой (±e).
Наиболее простым в реализации является метод сравнения дискретных отсчетов сигнала с фиксированными уровнями (рис.1).
Рисунок 1- Метод сравнения с фиксированными уровнями.
Здесь сплошной линией показана реализация сигнала и ступенчатая аппроксимация сигнала, а кружками отмечены выборки, подлежащие передаче.
Более эффективны адаптивные процедуры апертурного сжатия, использующие плавающую апертуру d, которая на очередном участке аппроксимации устанавливается определенным образом относительно последней выбранной существенной ординаты сигнала, и таким образом отслеживается изменение амплитуды сигнала. Процедура поиска существенных ординат здесь определяется характером аппроксимации: интерполяцией или экстраполяцией. В качестве аппроксимирующих функций чаще применяют полиномы нулевой, реже — первой степени.
Для плавающей апертуры выбор существенной ординаты при поступлении каждого i-го отсчета в последовательностиосуществляется в результате сравнения разности определенных для конкретного алгоритма значений сигнала с апертурой d. Процедура экстраполяции нулевого порядка (ЭНП) иллюстрируется на рис.2, а интерполяции (ИНП)- на рис.3.
Рисунок 2- Экстраполяция нулевого порядка.
Апертурная аппроксимация сигнала может быть реализована также путем сравнения отсчетов сигнала с его представлением алгебраическими полиномами первой степени. Использование более сложных функций затрудняет обработку сигнала в реальном времени и не дает ощутимого выигрыша в коэффициенте сжатия.
В отличие от алгоритма ЭНП при экстраполяции первого порядка (ЭПП) ось апертуры, являющаяся экстраполирующей прямой, располагается по линии, соединяющей первую ординату нового участка аппроксимации с предсказанным значением предыдущей ординаты.
Рисунок 2- Интерполяция нулевого порядка.
Все отсчеты сигнала, попавшие в построенный таким образом коридор, считаются избыточными, а первый, вышедший за его пределы, начинает следующий участок аппроксимации. На выход алгоритма сжатия могут передаваться предсказанное значение последнего отсчета аппроксимирующей прямой и длительность соответствующего участка аппроксимации.
В алгоритме интерполяции первого порядка (ИПП), в отличие от ЭПП, ось апертуры с приходом каждого следующего отсчета меняет свое положение. Вначале она проходит через первую и третью ординату текущего участка аппроксимации. Если второй отсчет попал в апертуру, то он считается избыточным. С приходом следующего отсчета уравнение аппроксимирующей прямой будет задаваться уже первой и текущей ординатой, и так до тех пор, пока хотя бы одна из промежуточных ординат не выйдет за пределы коридора. Тогда текущая ордината начнет новый участок аппроксимации, а предыдущая вместе с параметром длительностибудет передана на выход алгоритма сжатия.
Графическая интерпретация алгоритмов сжатия ЭПП и ИПП приведена на рис.3 и 4.
Рисунок 3- Экстраполяция первого порядка.
Рисунок 4- Интерполяция первого порядка.
Плавающая апертура может задаваться постоянной, постоянной со сдвигом или переменной, что во многом определяет сложность и эффективность конкретного алгоритма сжатия.
Рассмотренные алгоритмы сжатия используют однопараметрическую адаптацию по интервалу аппроксимации, поскольку последний автоматически определяется при поиске существенных отсчетов. Дальнейшее повышение эффективности сжатия ЭКГ может быть достигнуто за счет использования двухпараметрической адаптации, позволяющей автоматически определять как длительность интервала аппроксимации, так и степень аппроксимирующего полинома. При разработке алгоритмов анализа ЭКС, основанных на структурных методах, использующих сегментацию сигнала, изменение порядка аппроксимации на отдельных его отрезках может стать информативным признаком распознавания элементов ЭКГ различной крутизны. Наиболее широкое применение среди двухпараметрических алгоритмов сжатия ЭКС получили алгоритм сжатия АZТЕС (AmplitudeZoneTimeEpochCoding) и CORTES.
В основе алгоритма AZ 'ТЕС лежит представление сигнала последовательностью горизонтальных и наклонных отрезков прямой, называемых линия L или склон S соответственно. На первом этапе сжатия используется алгоритм ИНП, на выходе которого каждая j-я линия представлена парой чисел , На втором этапе каждый интервал аппроксимации анализируетсяна возможность включения его в новый, задаваемый полиномом первой степени. Для этого вводятся следующие дополнительные условия. Линия L, называется плато, если ее длительность превышает некоторый порог а. Линия Lназывается экстремумом, если
.
Если линия L, — плато или экстремум, то она называется границей. Далее, если в последовательности линий L0, L1,...Lт линии L0 и Lт являются границами, а все промежуточные не относятся к таковым, то последовательность горизонтальных линий L1, L2,…Lm-1, объединяется вновый k-й участок аппроксимации,и заменяется слоном S (рис.5).
Таким образом достигается дополнительное сжатие тех участков ступенчатой кривой, которые
соответствуют большей крутизне сигнала. На выходе алгоритма каждый элемент сжатого описания задается тремя элементами:
-амплитудой (ординатой линии или размахом склона);
-длительностью соответствующего участка аппроксимации;
-признаком линия/склон L/S, в качестве которого можно, например, использовать знак длительности, принимая его отрицательным для склонов.
Для получения более точного описания ЭКС в области QRS-комплекса был разработан алгоритм СОRТЕS, представляющий собой комбинацию двух алгоритмов: алгоритма АZТЕС,
Рисунок 5- Алгоритм сжатия AZTEC.
используемого для кодирования изоэлектрических участков, и алгоритма выделения точек локальных экстремумов сигнала, применяемого для сжатия высокочастотных сегментов сигнала.Работа последнего основана на контроле знака производной исходного сигнала. На каждом шаге алгоритма рассматривается трехточечная выборка отсчетов (u0, u1, u2), в кото рой u0 принимается за исходную существенную и передается на выход алгоритма. Далее из u1 и u2 выбирается и1, если выполняется условие
и u2 в противном случае. Затем вновь переданная ордината (u1 или u2) принимается за исходную и процесс обработки повторяется. При таком способе кодирования в два раза снижается частота следования отсчетов, сохраняются все точки изменениязнака производной, однако теряется информация о реальных временных соотношениях существенных отсчетов, из-за чего возможны искажения сигнала.
Алгоритм СОRТЕS построен таким образом, что две описанные процедуры сжатия работают параллельно, но выходные данные алгоритма, использующего контроль производной исходного сигнала, передаются лишь в том случае, если длина горизонтальной линии, строящейся по алгоритму АZТЕС, не превышает некоторого порогового значения. Для восстановления сигнала в исходном виде необходимы дополнительные вычисления, связанные с переходом от одного способа кодирования к другому. Кроме того, следует проводить сглаживание восстановленного сигнала посредством низкочастотной фильтрации из-за возможныхискажений на изоэлектрических участках ЭКС.
Все рассмотренные методы и алгоритмы сжатия относятся в основном к наиболее распространенным в практике оперативной обработки ЭКС.