Вейвлет Хаара

редактировать

Вейвлет Хаара Две итерации двумерного вейвлета Хаара на изображении Ленна. Исходное изображение подвергается высокочастотной фильтрации, что дает три субизображения с коэффициентами детализации (вверху справа: по горизонтали, внизу слева: по вертикали и внизу справа: по диагонали). Затем оно подвергается фильтрации нижних частот и масштабируется с уменьшением, в результате чего получается фрагмент изображения с коэффициентами аппроксимации (вверху слева); процесс фильтрации повторяется еще раз на этом приближенном изображении.

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

Последовательность Хаара была предложена в 1909 году Альфредом Хааром. Хаар использовал эти функции, чтобы дать пример ортонормированной системы для пространства квадратично интегрируемых функций на единичном интервале [0, 1]. Изучение вейвлетов и даже термина «вейвлет» появилось гораздо позже. Как частный случай вейвлета Добеши, вейвлет Хаара также известен как Db1 .

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

Материнская функция вейвлета вейвлета Хаара ψ (t) {\ displaystyle \ psi (t)}\ psi (t) можно описать как

ψ (t) = {1 0 ≤ t < 1 2, − 1 1 2 ≤ t < 1, 0 otherwise. {\displaystyle \psi (t)={\begin{cases}1\quad 0\leq t<{\frac {1}{2}},\\-1{\frac {1}{2}}\leq t<1,\\0{\t_dv{otherwise.}}\end{cases}}}\ psi (t) = {\ begin {cases} 1 \ quad 0 \ leq t <{\ frac {1} {2}}, \\ - 1 {\ frac {1 } {2}} \ leq t <1, \\ 0 {\ t_dv {в противном случае.}} \ End {cases}}

Его функция масштабирования φ (t) {\ displaystyle \ varphi (t)}\ varphi (t) можно описать как

φ (t) = {1 0 ≤ t < 1, 0 otherwise. {\displaystyle \varphi (t)={\begin{cases}1\quad 0\leq t<1,\\0{\t_dv{otherwise.}}\end{cases}}}{\ displaystyle \ varphi (t) = {\ begin {cases} 1 \ quad 0 \ leq t <1, \\ 0 {\ t_dv {в противном случае.}} \ End {case}}}
Содержание
  • 1 Функции Хаара и система Хаара
  • 2 Свойства вейвлетов Хаара
  • 3 Система Хаара на единичном интервале и родственные системы
    • 3.1 Система Фабера – Шаудера
    • 3.2 Система Франклина
  • 4 Матрица Хаара
  • 5 Преобразование Хаара
    • 5.1 Введение
    • 5.2 Свойство
    • 5.3 Преобразование Хаара и обратное преобразование Хаара
    • 5.4 Пример
    • 5.5 Приложение
  • 6 См. Также
  • 7 Примечания
  • 8 Ссылки
  • 9 Внешние ссылки
    • 9.1 Преобразование Хаара
Функции Хаара и система Хаара

Для каждой пары n, k целых чисел в Z, функция Хаара ψ n, k определяется на вещественной прямой Rформулой

ψ n, k (t) = 2 n / 2 ψ (2 n t - k), t ∈ R. {\ displaystyle \ psi _ {n, k} (t) = 2 ^ {n / 2} \ psi (2 ^ {n} tk), \ quad t \ in \ mathbf {R}.}\ psi _ {{n, k}} (t) = 2 ^ {{n / 2}} \ psi (2 ^ {n} tk), \ quad t \ in {\ mathbf {R}}.

Эта функция поддерживается на правом открытом интервале I n, k = [k 2, (k + 1) 2), то есть исчезает вне этого интервала. Он имеет интеграл 0 и норму 1 в гильбертовом пространстве L(R),

∫ R ψ n, k (t) dt = 0, ‖ ψ n, k ‖ L 2 (R) 2 = ∫ R ψ n, k ( t) 2 dt = 1. {\ displaystyle \ int _ {\ mathbf {R}} \ psi _ {n, k} (t) \, dt = 0, \ quad \ | \ psi _ {n, k} \ | _ {L ^ {2} (\ mathbf {R})} ^ {2} = \ int _ {\ mathbf {R}} \ psi _ {n, k} (t) ^ {2} \, dt = 1.}\ int _ {{{\ mathbf {R}}}} \ psi _ {{n, k}} (t) \, dt = 0, \ quad \ | \ psi _ {{n, k}} \ | _ {{L ^ {2} ({\ mathbf {R}})}} ^ {2} = \ int _ { {{\ mathbf {R}}}} \ psi _ {{n, k}} (t) ^ {2} \, dt = 1.

Функции Хаара попарно ортогональны,

∫ R ψ n 1, k 1 (t) ψ n 2, k 2 (t) dt = δ n 1, n 2 δ k 1, к 2, {\ displaystyle \ int _ {\ mathbf {R}} \ psi _ {n_ {1}, k_ {1}} (t) \ psi _ {n_ {2}, k_ {2}} (t) \, dt = \ delta _ {n_ {1}, n_ {2}} \ delta _ {k_ {1}, k_ {2}},}\ int _ {{{\ mathbf {R}}}} \ psi _ {{n_ {1}, k_ {1}}} (t) \ psi _ {{n_ {2}, k_ {2}} } (t) \, dt = \ delta _ {{n_ {1}, n_ {2}}} \ delta _ {{k_ {1}, k_ {2}}},

где δ i, j представляет Дельта Кронекера. Вот причина ортогональности: когда два поддерживающих интервала I n 1, k 1 {\ displaystyle I_ {n_ {1}, k_ {1}}}I_ {{n_ {1}, k_ {1}}} и I n 2, k 2 {\ displaystyle I_ {n_ {2}, k_ {2}}}I _ {{n_ {2}, k_ {2}}} не равны, то они либо не пересекаются, либо меньшая из двух опор, скажем I n 1, k 1 {\ displaystyle I_ {n_ {1}, k_ {1}}}I_ {{n_ {1}, k_ {1}}} , содержится в нижней или верхней половине другого интервала, на котором функция ψ N 2, k 2 {\ displaystyle \ psi _ {n_ {2}, k_ {2}}}\ psi _ {{n_ {2}, k_ {2}}} остается постоянным. В этом случае следует, что произведение этих двух функций Хаара кратно первой функции Хаара, следовательно, произведение имеет интеграл 0.

Система Хаара на действительной прямой является набор функций

{ψ n, k (t); n ∈ Z, k ∈ Z}. {\ displaystyle \ {\ psi _ {n, k} (t) \ ;; \; n \ in \ mathbf {Z}, \; k \ in \ mathbf {Z} \}.}\ {\ psi _ {{n, k}} (t) \ ;; \; n \ in {\ mathbf {Z}}, \; k \ in {\ mathbf {Z}} \}.

Это завершено в L (R ): система Хаара на линии является ортонормированным базисом в L (R ).

Свойства вейвлета Хаара

Вейвлет Хаара имеет несколько примечательных свойств:

  1. Любая непрерывная вещественная функция с компактной опорой может быть равномерно аппроксимирована линейными комбинациями из φ (T), φ (2 T), φ (4 T),…, φ (2 nt),… {\ displaystyle \ varphi (t), \ varphi (2t), \ varphi (4t), \ dots, \ varphi (2 ^ {n} t), \ dots}{\ displaystyle \ varphi (t), \ varphi (2t), \ varphi (4t), \ dots, \ varphi (2 ^ {n} t), \ dots} и их смещенные функции. Это распространяется на те функциональные пространства, где любую функцию в нем можно аппроксимировать непрерывными функциями.
  2. Любая непрерывная действительная функция на [0, 1] может быть аппроксимирована равномерно на [0, 1] линейными комбинациями постоянной функции 1, ψ (t), ψ (2 t), ψ (4 t),…, ψ (2 nt),… {\ displaystyle \ psi (t), \ psi (2t), \ psi (4t), \ dots, \ psi (2 ^ {n} t), \ dots}\ psi (t), \ psi (2t), \ psi (4t), \ dots, \ psi (2 ^ {n} t), \ dots и их смещенные функции.
  3. Ортогональность в виде

∫ - ∞ ∞ 2 (n + n 1) / 2 ψ (2 nt - k) ψ (2 n 1 t - k 1) dt = δ n, n 1 δ k, k 1. {\ displaystyle \ int _ {- \ infty} ^ {\ infty} 2 ^ {(n + n_ {1}) / 2} \ psi (2 ^ {n} tk) \ psi (2 ^ {n_ {1} } t-k_ {1}) \, dt = \ delta _ {n, n_ {1}} \ delta _ {k, k_ {1}}.}\ int _ {{- \ infty}} ^ {{\ infty}} 2 ^ {{(n + n_ {1}) / 2}} \ psi (2 ^ {n} tk) \ psi (2 ^ {{n_ {1}}} t-k_ {1}) \, dt = \ delta _ {{n, n_ {1}} } \ delta _ {{k, k_ {1}}}.

Здесь δ i, j представляет дельту Кронекера. Двойственная функция к ψ (t) - это сама функция ψ (t).

  1. Функции вейвлета / масштабирования с различным масштабом n имеют функциональную взаимосвязь: поскольку
φ (t) = φ (2 t) + φ (2 t - 1) ψ (t) = φ (2 t) - φ (2 T - 1), {\ Displaystyle {\ begin {align} \ varphi (t) = \ varphi (2t) + \ varphi (2t-1) \\ [. 2em] \ psi (t) = \ varphi (2t) - \ varphi (2t-1), \ end {align}}}{\ displaystyle {\ begin {align} \ varphi (t) = \ varphi (2t) + \ varphi (2t-1) \\ [.2em] \ psi (t) = \ varphi (2t) - \ varphi (2t-1), \ end {align}}}
следует, что коэффициенты шкалы n могут быть вычислены с помощью коэффициентов шкалы n + 1:
Если χ вес (К, N) знак равно 2 N / 2 ∫ - ∞ ∞ Икс (T) φ (2 NT - К) dt {\ Displaystyle \ чи _ {ш} (к, п) = 2 ^ {п / 2 } \ int _ {- \ infty} ^ {\ infty} x (t) \ varphi (2 ^ {n} tk) \, dt}{\ displaystyle \ chi _ {w} (k, n) = 2 ^ {n / 2} \ int _ {- \ infty} ^ {\ infty} x (t) \ varphi (2 ^ {n} tk) \, dt}
и X w (k, n) = 2 n / 2 ∫ - ∞ ∞ Икс (T) ψ (2 NT - К) dt {\ Displaystyle \ mathrm {X} _ {w} (к, п) = 2 ^ {п / 2} \ int _ {- \ infty} ^ {\ infty} x (t) \ psi (2 ^ {n} tk) \, dt}{\ displaystyle \ mathrm {X} _ {w} (k, n) = 2 ^ {n / 2} \ int _ {- \ infty} ^ {\ infty} x (t) \ psi (2 ^ {n} tk) \, dt}
, тогда
χ w (k, n) = 2 - 1/2 (χ w (2 k, п + 1) + χ вес (2 к + 1, п + 1)) {\ Displaystyle \ чи _ {ш} (к, п) = 2 ^ {- 1/2} {\ bigl (} \ чи _ { w} (2k, n + 1) + \ chi _ {w} (2k + 1, n + 1) {\ bigr)}}\ chi _ {w} (k, n) = 2 ^ { {-1/2}} {\ bigl (} \ chi _ {w} (2k, n + 1) + \ chi _ {w} (2k + 1, n + 1) {\ bigr)}
X w (k, n) = 2 - 1/2 (χ w (2 k, n + 1) - χ w (2 k + 1, n + 1)). {\ displaystyle \ mathrm {X} _ {w} (k, n) = 2 ^ {- 1/2} {\ bigl (} \ chi _ {w} (2k, n + 1) - \ chi _ {w } (2k + 1, n + 1) {\ bigr)}.}\ mathrm {X} _ {w} (k, n) = 2 ^ {{- 1/2}} {\ bigl (} \ chi _ {w} (2k, n + 1) - \ chi _ {w} (2k + 1, n + 1) {\ bigr)}.
Система Хаара на единичном интервале и связанные системы

В этом разделе обсуждение ограничено единичным интервалом [0, 1] и функциям Хаара, которые поддерживаются на [0, 1]. Система функций, рассмотренная Хааром в 1910 году и называемая в этой статье системой на [0, 1], состоит из подмножества вейвлетов Хаара, определенных как

{t ∈ [0, 1] ↦ ψ n, k (t); n ∈ N ∪ {0}, 0 ≤ k < 2 n }, {\displaystyle \{t\in [0,1]\mapsto \psi _{n,k}(t)\;;\;n\in \mathbb {N} \cup \{0\},\;0\leq k<2^{n}\},}\ {t \ in [0,1] \ mapsto \ psi _ {{n, k}} (t) \ ;; \; n \ in \ math bb {N} \ cup \ {0 \}, \; 0 \ leq k <2 ^ {n} \},

с добавлением постоянной функции 1 на [0, 1].

В терминах гильбертова пространства эта система Хаара на [0, 1] является полной ортонормированной системой, т. Е. ортонормированным базисом, для пространства L ([0, 1]) квадратично интегрируемых функций на единичном интервале.

Система Хаара на [0, 1] - с постоянной функцией 1 в качестве первого элемента, за которой следуют функции Хаара, упорядоченные в соответствии с лексикографическим упорядочением пар (n, k) - далее монотонный базис Шаудера для пространства L ([0, 1]), когда 1 ≤ p < ∞. This basis is безусловный при 1 < p < ∞.

Имеется связанная система Радемахера, состоящая из сумм функций Хаара,

rn (t) = 2 - n / 2 ∑ k = 0 2 n - 1 ψ n, К (T), T ∈ [0, 1], N ≥ 0. {\ Displaystyle r_ {n} (т) = 2 ^ {- n / 2} \ sum _ {k = 0} ^ {2 ^ {n } -1} \ psi _ {n, k} (t), \ quad t \ in [0,1], \ n \ geq 0.}r_ {n} (t) = 2 ^ {{- n / 2}} \ sum _ {{k = 0}} ^ {{2 ^ {n} -1}} \ psi _ {{n, k}} (t), \ quad t \ in [0,1], \ n \ geq 0.

Обратите внимание, что | r n (t) | = 1 на [0, 1). Это ортонормированная система, но не полная. На языке теории вероятностей последовательность Радемахера является экземпляром последовательности независимых Бернулли случайных величин со средним значением . 0. Неравенство Хинчина выражает тот факт, что во всех пространствах L ([0, 1]), 1 ≤ p < ∞, the Rademacher sequence is эквивалентно базису единичного вектора в ℓ. В частности, замкнутая линейная оболочка последовательности Радемахера в L ([0, 1]), 1 ≤ p < ∞, is , изоморфна.

Система Фабера – Шаудера

Система Фабера – Шаудера - это семейство непрерывных функций на [0, 1], состоящее из постоянной функции 1, и кратных неопределенных интегралов функций в системе Хаара на [0, 1], выбранных с нормой 1 в максимальной норме. Эта система начинается с s 0= 1, затем s 1 (t) = t - неопределенный интеграл, равный нулю в 0 функции 1, первого элемента системы Хаара на [0, 1]. Далее, для любого целого числа n ≥ 0 функции s n, k определяются формулой

sn, k (t) = 2 1 + n / 2 ∫ 0 t ψ n, k (u) du, t ∈ [0, 1], 0 ≤ k < 2 n. {\displaystyle s_{n,k}(t)=2^{1+n/2}\int _{0}^{t}\psi _{n,k}(u)\,du,\quad t\in [0,1],\ 0\leq k<2^{n}.}s _ {{n, k}} (t) = 2 ^ {{1 + n / 2}} \ int _ { 0} ^ {t} \ psi _ {{n, k}} (u) \, du, \ quad t \ in [0,1], \ 0 \ leq k <2 ^ {n}.

Эти функции s n, k непрерывны, кусочно-линейны, поддерживаются интервалом I n, k, который также поддерживает ψ n, k. Функция s n, k равна 1 в средней точке x n, k интервала I n, k, линейная на обеих половинах этого интервала.. Он везде принимает значения от 0 до 1.

Система Фабера – Шаудера - это базис Шаудера для пространства C ([0, 1]) непрерывных функций на [0, 1]. Для каждого f в C ([0, 1]) частичная сумма

fn + 1 = a 0 s 0 + a 1 s 1 + ∑ m = 0 n - 1 (∑ k = 0 2 m - 1 am, ksm, k) ∈ C ([0, 1]) {\ displaystyle f_ {n + 1} = a_ {0} s_ {0} + a_ {1} s_ {1} + \ sum _ {m = 0} ^ {n-1} {\ Bigl (} \ sum _ {k = 0} ^ {2 ^ {m} -1} a_ {m, k} s_ {m, k} {\ Bigr)} \ in C ( [0,1])}f _ {{n + 1}} = a_ {0} s_ {0} + a_ {1} s_ {1} + \ sum _ {{m = 0}} ^ {{n-1}} {\ Bigl (} \ sum _ {{k = 0}) } ^ {{2 ^ {m} -1}} a _ {{m, k}} s _ {{m, k}} {\ Bigr)} \ в C ([0,1])

разложения в ряд функции f в системе Фабера – Шаудера является непрерывной кусочно-линейной функцией, которая согласуется с f в 2 + 1 точках k 2, где 0 ≤ k ≤ 2. Далее формула

fn + 2 - fn + 1 = ∑ k = 0 2 n - 1 (f (xn, k) - fn + 1 (xn, k)) sn, k = ∑ k Знак равно 0 2 n - 1 an, ksn, k {\ displaystyle f_ {n + 2} -f_ {n + 1} = \ sum _ {k = 0} ^ {2 ^ {n} -1} {\ bigl ( } f (x_ {n, k}) - f_ {n + 1} (x_ {n, k}) {\ bigr)} s_ {n, k} = \ sum _ {k = 0} ^ {2 ^ { n} -1} a_ {n, k} s_ {n, k}}f _ {{n + 2}} - f _ {{n + 1}} = \ sum _ { {k = 0}} ^ {{2 ^ {n} -1}} {\ bigl (} f (x _ {{n, k}}) - f _ {{n + 1}} (x _ {{n, k }}) {\ bigr)} s _ {{n, k}} = \ sum _ {{k = 0}} ^ {{2 ^ {n} -1}} a _ {{n, k}} s _ {{ n, k}}

позволяет шаг за шагом вычислять расширение f. Поскольку f является равномерно непрерывным, последовательность {f n } равномерно сходится к f. Отсюда следует, что разложение f в ряд Фабера – Шаудера сходится в C ([0, 1]), и сумма этого ряда равна f.

Система Франклина

Система Франклина получается из системы Фабера – Шаудера с помощью процедуры ортонормировки Грама – Шмидта. Поскольку система Франклина имеет ту же линейную оболочку, что и система Фабера – Шаудера, эта оболочка плотна в C ([0, 1]), следовательно, в L ([0, 1]). Таким образом, система Франклина является ортонормированным базисом для L ([0, 1]), состоящим из непрерывных кусочно-линейных функций. П. Франклин доказал в 1928 г., что эта система является базисом Шаудера для C ([0, 1]). Система Франклина также является безусловным базисом Шаудера для пространства L ([0, 1]), когда 1 < p < ∞. The Franklin system provides a Schauder basis in the дисковая алгебра A (D). Это было доказано в 1974 году Бочкаревым, после того как существование базиса дисковой алгебры оставалось открытым более сорока лет.

Конструкция Бочкарева базиса Шаудера в A (D) выглядит следующим образом: пусть f комплексная липшицева функция на [0, π]; тогда f является суммой ряда косинусов с абсолютно суммируемыми коэффициентами. Пусть T (f) - элемент A (D), определенный комплексным степенным рядом с теми же коэффициентами,

{f: x ∈ [0, π] → ∑ n = 0 ∞ an cos ⁡ (nx)} ⟶ {T (f): z → ∑ n = 0 ∞ anzn, | z | ≤ 1}. {\ Displaystyle \ влево \ {е: х \ в [0, \ пи] \ вправо \ сумма _ {п = 0} ^ {\ infty} а_ {п} \ соз (пх) \ вправо \} \ longrightarrow \ влево \ {T (f): z \ rightarrow \ sum _ {n = 0} ^ {\ infty} a_ {n} z ^ {n}, \ quad | z | \ leq 1 \ right \}.}\ left \ {f: x \ in [ 0, \ pi] \ rightarrow \ sum _ {{n = 0}} ^ {\ infty} a_ {n} \ cos (nx) \ right \} \ longrightarrow \ left \ {T (f): z \ rightarrow \ sum _ {{n = 0}} ^ {\ infty} a_ {n} z ^ {n}, \ quad | z | \ leq 1 \ right \}.

Базис Бочкарева для A (D) образуют образы под T функций из системы Франклина на [0, π]. Эквивалентное описание Бочкарева для отображения T начинается с расширения f до четной липшицевой функции g 1 на [−π, π], отождествляемой с липшицевой функцией на единичной окружности T. Затем, пусть g 2 будет сопряженной функцией для g 1, и определим T (f) как функцию в A (D), значение которой на Граница T области D равна g 1 + ig 2.

При работе с 1-периодическими непрерывными функциями или, скорее, с непрерывными функциями f на [0, 1], такими что f (0) = f (1), удаляют функцию s 1 (t) = t из системы Фабера – Шаудера, чтобы получить периодическую систему Фабера – Шаудера . Периодическая система Франклина получается ортонормировкой из периодической системы Фабера – Шаудера. Можно доказать результат Бочкарева о A (D), доказав, что периодическая система Франклина на [0, 2π] является базисом для банахова пространства A r, изоморфного A (D). Пространство A r состоит из сложных непрерывных функций на единичной окружности T, чья сопряженная функция также непрерывна.

Матрица Хаара

Матрица Хаара 2 × 2, связанная с вейвлетом Хаара, равна

H 2 = [1 1 1 - 1]. {\ displaystyle H_ {2} = {\ begin {bmatrix} 1 1 \\ 1 -1 \ end {bmatrix}}.}H_ {2} = {\ begin {bmatrix} 1 1 \\ 1 -1 \ end {bmatrix}}.

Используя дискретное вейвлет-преобразование, можно преобразовать любую последовательность (a 0, a 1,…, a 2 n, a 2 n + 1) {\ displaystyle (a_ {0}, a_ {1}, \ dots, a_ {2n}, a_ {2n + 1})}(a_ {0}, a_ {1}, \ dots, a _ {{2n}}, a _ {{2n + 1}}) четной длины в последовательность двухкомпонентных векторов ((a 0, a 1),…, (a 2 n, a 2 n + 1)) {\ displaystyle \ left ( \ left (a_ {0}, a_ {1} \ right), \ dots, \ left (a_ {2n}, a_ {2n + 1} \ right) \ right)}\ left (\ left (a_ {0}, a_ {1} \ right), \ dots, \ left (a _ {{2n}}, a_ { {2n + 1}} \ right) \ right) . Если умножить каждый вектор вправо на матрицу H 2 {\ displaystyle H_ {2}}H_ {2} , получится результат ((s 0, d 0),…, (sn, dn)) {\ displaystyle \ left (\ left (s_ {0}, d_ {0} \ right), \ dots, \ left (s_ {n}, d_ {n} \ right) \ right)}\ left (\ left (s_ {0}, d_ {0} \ right), \ dots, \ left (s_ {n}, d_ {n} \ right) \ right) одной стадии быстрого вейвлет-преобразования Хаара. Обычно разделяют последовательности s и d и продолжают преобразовывать последовательность s. Последовательность s часто называют частью средних значений, тогда как последовательность d известна как часть деталей.

Если длина последовательности кратна четырем, можно построить блоки из 4 элементов и преобразовать их в аналогично с матрицей Хаара 4 × 4

H 4 = [1 1 1 1 1 1 - 1 - 1 1 - 1 0 0 0 0 1 - 1], {\ displaystyle H_ {4} = {\ begin { bmatrix} 1 1 1 1 \\ 1 1 -1 -1 \\ 1 -1 0 0 \\ 0 0 1 -1 \ end {bmatrix}},}H_ {4 } = {\ begin {bmatrix} 1 1 1 1 \\ 1 1 -1 -1 \\ 1 -1 0 0 \\ 0 0 1 -1 \ end {bmatrix}},

, которое объединяет два этапа быстрого вейвлет-преобразования Хаара.

Сравните с матрицей Уолша, которая является нелокализованной матрицей 1 / –1.

Как правило, матрицу Хаара 2N × 2N можно получить с помощью следующего уравнения.

H 2 N = [HN ⊗ [1, 1] IN ⊗ [1, - 1]] {\ displaystyle H_ {2N} = {\ begin {bmatrix} H_ {N} \ otimes [1,1] \ \ I_ {N} \ otimes [1, -1] \ end {bmatrix}}}H _ {{2N}} = {\ begin {bmatrix} H _ {{N}} \ otimes [1,1] \\ I _ {{N}} \ otimes [1, -1] \ end {bmatrix}}
где IN = [1 0… 0 0 1… 0 ⋮ ⋮ ⋱ ⋮ 0 0… 1] {\ displaystyle I_ {N} = {\ begin {bmatrix} 1 0 \ dots 0 \\ 0 1 \ dots 0 \\\ vdots \ vdots \ ddots \ vdots \\ 0 0 \ dots 1 \ end {bmatrix}}}I _ {{N}} = {\ begin {bmatrix} 1 0 \ dots 0 \\ 0 1 \ dots 0 \\\ vdots \ vdots \ ddots \ vdots \\ 0 0 \ dots 1 \ end {bmatrix}} и ⊗ {\ displaystyle \ otimes}\ otimes - это произведение Кронекера.

произведение Кронекера из A ⊗ B {\ displaystyle A \ otimes B}A \ otimes B , где A {\ displaystyle A}A - это матрица m × n, а B {\ displaystyle B}B - Матрица ap × q, выражается как

A ⊗ B = [a 11 B… a 1 n B ⋮ ⋱ ⋮ am 1 B… amn B]. {\ displaystyle A \ otimes B = {\ begin {bmatrix} a_ {11} B \ dots a_ {1n} B \\\ vdots \ ddots \ vdots \\ a_ {m1} B \ dots a_ {mn} B \ end {bmatrix}}.}A \ otimes B = {\ begin {bmatrix} a _ {{ 11}} B \ dots a _ {{1n}} B \\\ vdots \ ddots \ vdots \\ a _ {{m1}} B \ dots a _ {{mn}} B \ end {bmatrix}}.

Ненормированная матрица Хаара из 8 точек H 8 {\ displaystyle H_ {8}}H_ {8} показана ниже

H 8 = [1 1 1 1 1 1 1 1 1 1 1 1 - 1 - 1 - 1 - 1 1 1 - 1 - 1 0 0 0 0 0 0 0 0 0 1 1 - 1 - 1 1 - 1 0 0 0 0 0 0 0 0 1 - 1 0 0 0 0 0 0 0 0 1 - 1 0 0 0 0 0 0 0 0 1 - 1]. {\ Displaystyle Н_ {8} = {\ BEGIN {bmatrix} 1 1 1 1 1 1 1 1 \\ 1 1 1 1 -1 -1 -1 -1 \\ 1 1 -1 -1 0 0 0 0 \\ 0 0 0 0 1 1 -1 -1 \\ 1 -1 0 0 0 0 0 0 \\ 0 0 1 -1 0 0 0 0 \\ 0 0 0 0 1 -1 0 0 \\ 0 0 0 0 0 0 1 -1 \ end {bmatrix}}.}H _ {{8}} = {\ begin {bmatrix} 1 1 1 1 1 1 1 1 1 1 \ \ 1 1 1 1 -1 -1 -1 -1 \\ 1 1 -1 -1 0 0 0 0 \\ 0 0 0 0 0 1 1 -1 -1 \\ 1 -1 0 0 0 0 0 0 \\ 0 0 1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.

Обратите внимание, что приведенная выше матрица является ненормализованной матрицей Хаара. Матрица Хаара, требуемая преобразованием Хаара, должна быть нормализована.

Из определения матрицы Хаара H {\ displaystyle H}Hможно заметить, что, в отличие от преобразования Фурье, H {\ displaystyle H}Hимеет только действительные элементы (например, 1, -1 или 0) и не является симметричным.

В качестве примера возьмем 8-точечную матрицу Хаара H 8 {\ displaystyle H_ {8}}H_ {8} . Первая строка H 8 {\ displaystyle H_ {8}}H_ {8} измеряет среднее значение, а вторая строка H 8 {\ displaystyle H_ {8}}H_ {8} измеряет низкочастотную составляющую входного вектора. Следующие две строки чувствительны к первой и второй половине входного вектора соответственно, что соответствует умеренным частотным компонентам. Остальные четыре строки чувствительны к четырем частям входного вектора, которые соответствуют высокочастотным компонентам.

Преобразование Хаара

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

Введение

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

Преобразование Хаара выводится из матрицы Хаара. Пример матрицы преобразования Хаара 4x4 показан ниже.

Ч 4 = 1 2 [1 1 1 1 1 1 - 1 - 1 2 - 2 0 0 0 0 2 - 2] {\ displaystyle H_ {4} = {\ frac {1} {2}} {\ begin {bmatrix} 1 1 1 1 \\ 1 1 -1 -1 \\ {\ sqrt {2}} - {\ sqrt {2}} 0 0 \\ 0 0 {\ sqrt {2}} - {\ sqrt {2}} \ end {bmatrix}}}H_ {4} = {\ frac {1 } {2}} {\ begin {bmatrix} 1 1 1 1 \\ 1 1 -1 -1 \\ {\ sqrt {2}} - {\ sqrt {2}} 0 0 \\ 0 0 {\ sqrt {2}} - {\ sqrt {2}} \ end {bmatrix}}

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

Сравните с преобразованием Уолша, которое также равно 1 / –1, но не локализовано.

Свойство

Преобразование Хаара имеет следующие свойства

1. Нет необходимости в умножении. Для этого требуются только дополнения, а в матрице Хаара много элементов с нулевым значением, поэтому время вычислений невелико. Это быстрее, чем преобразование Уолша, матрица которого состоит из +1 и -1.
2. Длина входа и выхода одинакова. Однако длина должна быть степенью двойки, то есть N = 2 k, k ∈ N {\ displaystyle N = 2 ^ {k}, k \ in \ mathbb {N}}N = 2 ^ {k}, k \ in {\ mathbb {N}} .
3. Его можно использовать для анализа локализованных характеристик сигналов. Благодаря свойству ортогональности функции Хаара можно анализировать частотные компоненты входного сигнала.

Преобразование Хаара и обратное преобразование Хаара

Преобразование Хаара y n функции ввода n x n равно

yn = H nxn {\ displaystyle y_ {n} = H_ {n} x_ {n}}y_{n}=H_{n}x_{n}

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

H = H *, H - 1 = HT, т.е. HHT = I {\ displaystyle H = H ^ {*}, H ^ {- 1} = H ^ {T}, {\ text {ie}} HH ^ {T} = I}{\ displaystyle H = H ^ {*}, H ^ {- 1} = H ^ {T}, {\ text {т.е. }} HH ^ {T} = I}
где I {\ displaystyle I}I - единичная матрица. Например, когда n = 4
H 4 TH 4 = 1 2 [1 1 2 0 1 1 - 2 0 1 - 1 0 2 1 - 1 0 - 2] ⋅ 1 2 [1 1 1 1 1 1 - 1 - 1 2 - 2 0 0 0 0 2 - 2] = [1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1] {\ displaystyle H_ {4} ^ {T} H_ {4} = { \ frac {1} {2}} {\ begin {bmatrix} 1 1 {\ sqrt {2}} 0 \\ 1 1 - {\ sqrt {2}} 0 \\ 1 -1 0 {\ sqrt {2}} \\ 1 -1 0 - {\ sqrt {2}} \ end {bmatrix}} \ cdot \; {\ frac {1} {2}} {\ begin {bmatrix} 1 1 1 1 \\ 1 1 1 -1 -1 \\ {\ sqrt {2}} - {\ sqrt {2}} 0 0 \\ 0 0 {\ sqrt {2}} - {\ sqrt {2}} \ end {bmatrix}} = {\ begin {bmatrix} 1 0 0 0 \\ 0 1 0 0 \\ 0 0 1 0 \\ 0 0 0 1 \ end {bmatrix}}}H_ {4} ^ {{T}} H_ {4} = {\ frac {1} {2}} {\ begin {bmatrix} 1 1 {\ sqrt {2}} 0 \\ 1 1 - {\ sqrt {2}} 0 \\ 1 -1 0 {\ sqrt { 2}} \\ 1 -1 0 - {\ sqrt {2}} \ end {bmatrix}} \ cdot \; {\ frac {1} {2}} {\ begin {bmatrix} 1 1 1 1 \\ 1 1 -1 -1 \\ {\ sqrt {2}} - {\ sqrt {2}} 0 0 \\ 0 0 {\ sqrt {2}} - {\ sqrt {2}} \ end {bmatrix}} = {\ begin {bmatrix } 1 0 0 0 \\ 0 1 0 0 \\ 0 0 1 0 \\ 0 0 0 1 \ end {bmatrix}}

Таким образом, обратное преобразование Хаара равно

xn = HT yn {\ displaystyle x_ {n} = H ^ {T} y_ {n}}x _ {{n}} = H ^ {{T}} y_ { {n}}

Пример

Коэффициенты преобразования Хаара = 4-точечный сигнал x 4 = [1, 2, 3, 4] T {\ displaystyle x_ {4} = [1,2,3, 4] ^ {T}}x _ {{4}} = [1,2,3,4] ^ {{T}} можно найти как

y 4 = H 4 x 4 = 1 2 [1 1 1 1 1 1 - 1 - 1 2 - 2 0 0 0 0 2 - 2] [1 2 3 4] = [5–2–1 / 2–1 / 2] {\ displaystyle y_ {4} = H_ {4} x_ {4} = {\ frac {1} {2}} { \ begin {bmatrix} 1 1 1 1 \\ 1 1 -1 -1 \\ {\ sqrt {2}} - {\ sqrt {2 }} 0 0 \\ 0 0 {\ sqrt {2}} - {\ sqrt {2}} \ end {bmatrix}} {\ begin {bmatrix} 1 \\ 2 \\ 3 \\ 4 \ end {bmatrix}} = {\ begin {bmatrix} 5 \\ - 2 \\ - 1 / {\ sqrt {2}} \\ - 1 / {\ sqrt {2}} \ end {bmatrix}}}y _ {{4}} = H_ {4} x_ {4} = {\ frac {1} {2}} {\ begin {bmatrix} 1 1 1 1 \\ 1 1 -1 -1 \\ {\ sqrt {2}} - {\ sqrt {2}} 0 0 \\ 0 0 {\ sqrt {2}} - {\ sqrt {2}} \ end {bmatrix}} {\ begin {bmatrix} 1 \\ 2 \\ 3 \\ 4 \ end {bmatrix}} = {\ begin { bmatrix} 5 \\ - 2 \\ - 1 / {\ sqrt {2}} \\ - 1 / {\ sqrt {2}} \ end {bmatrix}}

Входной сигнал может затем можно полностью восстановить с помощью обратного преобразования Хаара

x 4 ^ = H 4 T y 4 = 1 2 [1 1 2 0 1 1 - 2 0 1 - 1 0 2 1 - 1 0 - 2] [5 - 2 - 1/2 - 1/2] = [1 2 3 4] {\ displaystyle {\ hat {x_ {4}}} = H_ {4} ^ {T} y_ {4} = {\ frac {1} { 2}} {\ begin {bmatrix} 1 1 {\ sqrt {2}} 0 \\ 1 1 - {\ sqrt {2}} 0 \\ 1 -1 0 {\ sqrt {2}} \\ 1 -1 0 - {\ sqrt {2}} \ end {bmatrix}} {\ begin {bmatrix} 5 \\ - 2 \\ - 1 / {\ sqrt {2}} \\ - 1 / {\ sqrt {2}} \ end {bmatrix }} = {\ begin {bmatrix} 1 \\ 2 \\ 3 \\ 4 \ end {bmatrix}}}{\ hat {x _ {{4}}}} = H _ {{4}} ^ {{T}} y _ {{4}} = {\ frac {1} {2}} {\ begin { bmatrix} 1 1 {\ sqrt {2}} 0 \\ 1 1 - {\ sqrt {2}} 0 \\ 1 -1 0 {\ sqrt {2}} \\ 1 -1 0 - {\ sqrt {2}} \ end {bmatrix}} {\ begin {bmatrix} 5 \\ - 2 \\ - 1 / {\ sqrt {2}} \\ - 1 / {\ sqrt {2}} \ end {bmatrix}} = {\ begin {bmatrix} 1 \\ 2 \\ 3 \\ 4 \ end {bmatrix}}

Приложение

Современные камеры способны создавать изображения с разрешением в диапазоне десятков мегапикселей. Эти изображения необходимо сжать перед сохранением и передачей. Преобразование Хаара можно использовать для сжатия изображений. Основная идея состоит в том, чтобы передать изображение в матрицу, в которой каждый элемент матрицы представляет пиксель изображения. Например, матрица 256 × 256 сохраняется для изображения 256 × 256. JPEG Сжатие изображений включает разрезание исходного изображения на фрагменты изображения 8 × 8. Каждое суб-изображение представляет собой матрицу 8 × 8.

Требуется двумерное преобразование Хаара. Уравнение преобразования Хаара: B n = H n A n H n T {\ displaystyle B_ {n} = H_ {n} A_ {n} H_ {n} ^ {T}}{\ displaystyle B_ {n} = H_ {n} A_ {n} H_ {n} ^ {T}} , где A n {\ displaystyle A_ {n}}A_ {n} - матрица × n, а H n {\ displaystyle H_ {n}}H_ {n} - n- точечное преобразование Хаара. Обратное преобразование Хаара: A n = H n TB n H n {\ displaystyle A_ {n} = H_ {n} ^ {T} B_ {n} H_ {n}}{\ displaystyle A_ {n} = H_ {n} ^ {T} B_ {n} H_ {n}}

См. Также
Примечания
Ссылки
Внешние ссылки
Викискладе есть средства массовой информации, связанные с вейвлетом Хаара.

Преобразование Хаара

Последняя правка сделана 2021-05-22 10:03:47
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте