26.3. Вейвлетная очистка сигналов от шумов.

Рис. 26.3.1.

Вейвлетная очистка сигналов от шумов может использоваться для любых типов сигналов, но особенно эффективна для сигналов, имеющих в своем составе крутые перепады значений (скачки), местоположение которых является информацией, подлежащей сохранению. Для наглядности использования в этом процессе типовых функций преобразования wave(s) и iwave(S) в качестве примера выполним операцию очистки зашумленного меандра, представленного на рис. 26.3.1 на 500 точках координатной оси.

Подготовка преобразования. Для подготовки вейвлет-преобразования определяется количество уровней полной декомпозиции сигнала М, и массив дополняется (в данном случае нулями) до требуемой величины 2^М.

Рис.26.3.2.

Для оценки количества возможных уровней декомпозиции, в которых будет преобладать шум, можно выполнить БПФ массива. На рис. 26.3.2 представлен спектр массива в интервале 0- (256 отсчетов), т.е. на первой половине главного диапазона БПФ. При диадном делении спектра на каждом уровне декомпозиции, первый уровень детализирующих коэффициентов будет сформирован из высокочастотной части спектра сигнала от /2 до  (в односторонней физической шкале частот). Вторая часть спектра от 0 до /2 конвертируется в аппроксимирующие коэффициенты. На втором уровне декомпозиции аппроксимирующие коэффициенты диапазона 0-/2 также будут разделены пополам с преобразованием диапазона /4-/2 в детализирующие, а диапазона 0-/4 в аппроксимирующие коэффициенты 2-го уровня декомпозиции, и т.д. Это позволяет непосредственно по частотному спектру сигнала установить ориентировочную границу шумов и, соответственно уровни декомпозиции, в которых мощность шумов соизмерима и выше мощности сигнала.

Рис. 26.3.3.

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

Рис. 26.3.4.

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

Рис. 26.3.5.

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

После формирования новой строки qu данного уровня декомпозиции вейвлетного спектра можно заменить этой строкой соответствующий уровень разложения в полном массиве коэффициентов и визуально оценить результаты операции. Пример приведен на рис. 26.3.6.

Рис. 26.3.6.

После проведения аналогичной операции на втором декомпозиции получаем результат, приведенный на рис. 26.3.7. Он является конечным, т.к. попытка повторения операции на третьем уровне с данной моделью сигнала результатов не принесла.

Рис. 26.3.7.

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

Рис. 26.3.8.

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

ЛИТЕРАТУРА

2. Дремин И.Л. и др. Вейвлеты и их использование. / Успехи физических наук, 2001, т.171, № 5, стр. 465-501.

3. Дьяконов В., Абраменкова И. MATLAB. Обработка сигналов и изображений. Специальный справочник. – СПб.: Питер, 2002, 608 с.

5. Петухов А.П. Введение в теорию базисов всплесков. – СПб.: Изд. СПбГТУ, 1999, 132 с.

12. Смоленцев Н.К. Основы теории вейвлетов. Вейвлеты в Matlab. М.: LVR Пресс, 2005. – 304 с.

13. Переберин А.В. О систематизации вейвлет-преобразований. – Вычислительные методы и программирование, 2001, т.2.

14. Воробьев В.И., Грибунин В.Г. Теория и практика вейвлет-преобразования. – СПб, ВУС, 1999. 204 с.

39. Дьяконов В.П. Вейвлеты. От теории к практике. – М.: СОЛОН-Р, 2002. – 448 с.

Cайт автора Лекции Практикум

О замеченных ошибках и предложениях по дополнению: davpro@yandex.ru.

Copyright © 2009-2010 Davydov А.V.