Abstract

Abstract

Рассматривается задача построения адаптивного фильтра (АФ) на основе вейвлета Добеши. АФ смоделирован на ЦВМ. Приводятся результаты моделирования, которые позволяют оценить качество работы АФ.

Full Text

Рассмотрим адаптивную фильтрацию наблюдаемого сигнала, получаемого на выходе измерительного устройства (ИУ). Предположим, что на вход АФ на основе вейвлета Добеши поступает полезный сигнал , на который накладывается погрешность ИУ , так что входной наблюдаемый сигнал имеет вид . (1) Детерминированный полезный сигнал описывается выражением , (2) где – неизвестная амплитуда гармонического сигнала; – известная частота, ; Гц; – период гармонических колебаний. Корреляционная функция погрешности описывается соотношением , (3) где – неизвестная дисперсия помехи ; – неизвестный коэффициент нерегулярности помехи . Таким образом, неизвестны параметр сигнала и параметры , корреляционной функции . Необходимо путем соответствующей обработки (фильтрации) уменьшить влияние погрешности ИУ . Из соотношения (1) имеем: , , (4) где ; ; ; ; – интервал дискретности, с которым измеряется сигнал . Из выражений (2), (3) получим , , (5) , , (6) где ; . Рассмотрим построение АФ на основе ортогонального вейвлета Добеши. Вейвлетное преобразование применяется, в частности, для фильтрации сигналов. В вейвлетном преобразовании в качестве весовых коэффициентов значений сигнала выступают вейвлетные функции [1]. Вейвлеты характеризуются своим временным и частотными образами. Временной образ определяется некоторой psi-функцией времени. А частотный образ определяется её Фурье-образом , который задает огибающую спектра вейвлета [2]. Фурье-образ определяется выражением , (7) где – круговая частота. Выбираем для адаптивной фильтрации сигнала , вейвлет Добеши. АФ на основе вейвлета Добеши будем строить на базе пакета расширения системы MatLab 7.0 Wavelet Toolbox [3]. Функция возвращает очищенный от шума сигнал , полученный ограничением вейвлет-коэффициентов преобразования входного сигнала . Строка задает правило выбора порога. Было выбрано правило , что означает выбор инверсного порога. Параметр принимал значение , что соответствует выбору гибкого порога для удаления шумов путем ограничения вейвлет-коэффициентов. Строка определяет мультипликативное пороговое перемасштабирование (если шум вне пределов [0,1] или не белый). Было выбрано значение , что соответствует перемасштабированию с использованием единственной оценки уровня шума, основанному на коэффициентах первого уровня. Параметр задает имя вейвлета. Было выбрано значение , что соответствует вейвлету Добеши. Параметр означает уровень декомпозиции или количество коэффициентов вейвлет-преобразования. Было выбрано . Таким образом, АФ на основе вейвлета Добеши в системе MatLab 7.0 реализуется программно в виде команды ; Сигнал на выходе АФ , представляет собой сигнал , на входе фильтра, очищенный от шума. Процедура удаления шума состоит в подавлении составляющей шума , в сигнале , и восстановлении составляющей , и включает в себя три шага [2]: 1) разложение. Выбор вейвлета и уровня декомпозиции . Вейвлет-разложение сигнала , ; 2) детализация. Для каждого уровня от до выбирается определенный порог и применяется гибкий порог для детализирующих коэффициентов; 3) восстановление. Вейвлет-восстановление, основанное на исходных коэффициентах аппроксимации на уровне , модификация детальных коэффициентов на уровнях от до . Определим ошибку фильтрации для АФ на основе вейвлета Добеши. Имеем , . м (8) Обозначим через оценку среднеквадратического значения ошибки фильтрации для АФ на основе вейвлета Добеши. Величина вычисляется по формуле . (9) Величина характеризует качество работы АФ на основе вейвлета Добеши. АФ на основе вейвлета Добеши был смоделирован на ЦВМ с помощью программного продукта MatLab 7.0 [3]. Были приняты следующие значения параметров: Параметр принимал значения 100,60. Параметр принимал значения 1, 2, 5. Результаты моделирования приведены в таблице и на рис. 1–4. Результаты моделирования 100 1 5 0,0586 17,06 2 5 0,2157 9,27 5 5 1,31 3,80 60 1 5 0,12 8,33 2 5 0,46 4,35 5 5 2,83 1,77 На рис. 1 для ; приведены графики . Рис. 1. Графики для ; На рис. 2 для ; показаны графики . На рис. 3 для ; приведены графики . На рис. 4 для ; показаны графики . Рис. 2. Графики для ; Рис. 3. Графики для ; Рис. 4. Графики для ; Таким образом, в статье выполнено построение АФ на основе вейвлета Добеши, произведена оценка качества работы этого фильтра.