Матричные операции в вейвлетном базисе
Реферат - Математика и статистика
Другие рефераты по предмету Математика и статистика
(Rd), для любого j Z, f(x)Vj тогда и только тогда, когда
f(2x) Vj-1,
3. Для любого f L2(Rd), для любого k Zd, f(x)V0 тогда и только тогда, когда f(x-k)V0,
4. Существует масштабирующая (scaling) функция V0, что {(x-k)}kZd образует
базис Ритца в V0.
Для ортонормальных базисов можно переписать свойство 4 в виде:
4. Существует масштабирующая функция V0, что {(x-k)}kZd образует ортонормальный базис в V0.
Определим подпространство Wj как ортогональное дополнение к Vj в Vj-1,
, (1.2)
и представим пространство L2(Rd) в виде прямой суммы
(1.3)
Выбирая масштаб n, можем заменить последовательность (1.1) следующей последовательностью:
(1.4)
и получить
(1.5)
Если имеем конечное число масштабов, то, не нарушая общности, можно положить j=0 и рассматривать
, V0 L2(Rd) (1.6)
вместо (1.4). В числовой реализации подпространство V0 конечномерно.
Функция - так называемая масштабирующая (скейлинг-) функция. С ее помощью можно определить функцию - вейвлет - такую, что набор {(x-k)}kZ образует ортонормальный базис в W0. Тогда
, m=0..M-1. (1.7)
Из свойства 4 непосредственно следует, что, во-первых, функция может быть представлена в виде линейной комбинации базисных функций пространства V-1 . Так как функции {j,k(x)=2-j/2(2-jx-k)}kZ образуют ортонормальный базис в Vj, то имеем
. (1.8)
Вообще говоря, сумма в выражении (1.8) не обязана быть конечной. Можно переписать (1.8) в виде
, (1.9)
где
, (1.10)
а 2-периодическая функция m0 определяется следующим образом:
. (1.11)
Во-вторых, ортогональность {(x-k)}kZ подразумевает, что
(1.12)
и значит
(1.13)
и . (1.14)
Используя (1.9), получаем
(1.15)
и, рассматривая сумму в (1.15) по четным и нечетным индексам, имеем
. (1.16)
Используя 2-периодичность функции m0 и (1.14), после замены /2 на , получаем необходимое условие
(1.17)
для коэффициентов hk в (1.11). Заметив, что
(1.18)
и определив функцию следующим образом:
, (1.19)
где
, k=0,…,L-1 , (1.20)
или преобразование Фурье для
, (1.21)
где
, (1.22)
можно показать, что при каждом фиксированном масштабе jZ вейвлеты
{j,k(x)=2-j/2(2-jx-k)}kZ образуют ортонормальный базис пространства Wj.
Равенство (1.17) определяет пару квадратурных зеркальных фильтров (quadrature mirror filters, QMF) H и G, где и . Коэффициенты QMF H и G вычисляются с помощью решения системы алгебраических уравнений. Число L коэффициентов фильтра в (1.11) и (1.22) связано с числом исчезающих моментов М, и всегда четно.
Выбранный фильтр Н полностью определяет функции и и, таким образом, многомасштабный анализ. Кроме того, в правильно построенных алгоритмах значения функций и почти никогда не вычисляются. Благодаря рекурсивному определению вейвлетного базиса, все операции проводятся с квадратурными зеркальными фильтрами H и G, даже если в них используются величины, связаные с и .
2. БЫСТРОЕ ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЕ
После того, как вычислены коэффициенты hk и gk, т.е. выбран определенный вейвлет, можно проводить вейвлет-преобразование сигнала f(x), поскольку задан ортонормальный базис (j,k, j,k). Любая функция f(x)L2(R) полностью характеризуется ее вейвлет-коэффициентами разложения по этому базису и потому может быть представлена формулой
. (2.1)
Зададим все пределы суммирования в формуле (2.1). Функцию f(x) можно рассматривать на любом n-м уровне разрешения jn. Тогда разделение между ее усредненными значениями на этом уровне и флуктуациями вокруг них выглядят как
. (2.4)
На бесконечном интервале первая сумма может быть опущена, и в результате получается «