Реферат: Операторы в вейвлетном базисе

🛑 👉🏻👉🏻👉🏻 ИНФОРМАЦИЯ ДОСТУПНА ЗДЕСЬ ЖМИТЕ 👈🏻👈🏻👈🏻

Белорусский государственный университет
Факультет прикладной математики и информатики

ВВЕДЕНИЕ………..………………………………………………………..3
1. МНОГОМАСШТАБНЫЙ АНАЛИЗ И ВЕЙВЛЕТЫ………………...5
2. БЫСТРОЕ ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЕ….……………………...9
3. ДВУМЕРНЫЕ ВЕЙВЛЕТЫ…………………………………………..12
4. МАТРИЧНЫЕ ОПЕРАЦИИ………………………………………….13
4.1. Матричное умножение………………………………………...13
4.2. Обращение матрицы…………………………………………...16
4.3. Вычисление экспоненты, синуса и косинуса от матрицы.….16
ЛИТЕРАТУРА……………………………………………………………18
Вейвлет-преобразование сигналов (wavelet transform)
, теория которого оформилась в начале 90-х годов, является не менее общим по областям своих применений, чем классическое преобразование Фурье. Принцип ортогонального разложения по компактным волнам состоит в возможности независимого анализа функции на разных масштабах ее изменения. Вейвлет-представление сигналов (функций времени) является промежуточным между полностью спектральным и полностью временным представлениями.
Компактные волны относительно независимо были предложены в квантовой физике, физике электромагнитных явлений, математике, электронике и сейсмогеологии. Междисциплинарные исследования привели к новым приложениям данных методов, в частности, в сжатии образов для архивов и телекоммуникаций, в исследованиях турбулентности, в физиологии зрительной системы, в анализе радарных сигналов и предсказании землетрясений. К сожалению, объем русскоязычной научной литературы по тематике вейвлет-преобразований (да и нейронных сетей) относительно невелик.
Базовая идея восходит к временам 200-летней давности и принадлежит Фурье: аппроксимировать сложную функцию взвешенной суммой простых функций, каждая из которых, в свою очередь, получается из одной функции-прототипа. Эта функция-прототип выполняет роль строительного блока, а искомая аппроксимация получается комбинированием одинаковых по структуре блоков. При этом, если "хорошая" аппроксимация получается при использовании небольшого числа блоков, то тем самым достигается значительное уплотнение информации. В качестве таких блоков Фурье использовал синусоиды с различными периодами.
Что прежде всего отличает вейвлет-анализ от анализа Фурье? Основным недостатком Фурье-преобразования является его "глобальная" чувствительность к "локальным" скачкам и пикам функции. При этом модификация коэффициентов Фурье (например, обрезание высоких гармоник с целью фильтрации шума) вносит одинаковые изменения в поведение сигнала на всей области определения. Это особенность оказывается полезной для стационарных сигналов, свойства которых в целом мало меняются со временем.
При исследовании же нестационарных сигналов требуется использование некоторых локализованных во времени компактных волн, коэффициенты разложения по которым сохраняют информацию о дрейфе параметров аппроксимируемой функции. Первые попытки построения таких систем функций сводились к сегментированию сигнала на фрагменты ("окна") с применением разложения Фурье для этих фрагментов. Соответствующее преобразование - оконное преобразование Фурье - было предложено в 1946-47 годах Jean Ville и, независимо, Dennis Gabor. В 1950-70-х годах разными авторами было опубликовано много модификаций времени-частотных представлений сигналов.
В конце 70-х инженер-геофизик Морли (Jean Morlet) столкнулся с проблемой анализа сигналов, которые характеризовались высокочастотной компонентой в течение короткого промежутка времени и низкочастотными колебаниями при рассмотрении больших временных масштабов. Оконные преобразования позволяли проанализировать либо высокие частоты в коротком окне времени, либо низкочастотную компоненту, но не оба колебания одновременно. В результате был предложен подход, в котором для различных диапазонов частот использовались временные окна различной длительности. Оконные функции получались в результате растяжения-сжатия и смещения по времени гаусиана. Морли назвал эти базисные функции вейвлетами (wavelets)
- компактными волнами. В дальнейшем благодаря работам Мейера (Yves Meyer), Добеши (Ingrid Daubechies), Койфмана (Ronald Coifman), Маллы (Stephane Mallat) и других теория вейвлетов приобрела свое современное состояние.
Среди российских ученых, работавших в области теории вейвлетов, необходимо отметить С.Б. Стечкина, И.Я. Новикова, В.И. Бердышева.
Определение 1. Многомасштабный анализ (
multiresolutional
analysis)
– разложение гильбертова пространства L 2
(
R

d
)
, d
³1
, в последовательность замкнутых подпространств
2. Для любого f
Î
L 2
(
R

d
)
, для любого j
Î
Z

, f(
x)
Î
V j

тогда и только тогда, когда
3. Для любого f
Î
L 2
(
R

d
)
, для любого k
Î
Z

d

, f(
x)
Î
V 0

тогда и только тогда, когда f(
x-
k)
Î
V 0

,
4. Существует масштабирующая (
scaling)
функция j
Î
V 0

, что {
j(
x-
k)}
k

Î

Z

d

образует
Для ортонормальных базисов можно переписать свойство 4 в виде:
4’. Существует масштабирующая функция j
Î
V 0

, что {
j(
x-
k)}
k

Î

Z

d

образует ортонормальный базис в V 0

.
Определим подпространство W j

как ортогональное дополнение к V j

в V j
-1

,
и представим пространство L 2
(
R

d
)
в виде прямой суммы
Выбирая масштаб n
, можем заменить последовательность (1.1) следующей последовательностью:
Если имеем конечное число масштабов, то, не нарушая общности, можно положить j=0
и рассматривать
вместо (1.4). В числовой реализации подпространство V 0

конечномерно.
Функция j
- так называемая масштабирующая (скейлинг-) функция. С ее помощью можно определить функцию y
- вейвлет - такую, что набор {
y(
x-
k)}
k

Î

Z

образует ортонормальный базис в W 0

. Тогда
Из свойства 4’ непосредственно следует, что, во-первых, функция j
может быть представлена в виде линейной комбинации базисных функций пространства V -1

. Так как функции {
j
j,

k
(
x)=2 -

j/2

j(2 -

j

x-
k)}
k

Î

Z

образуют ортонормальный базис в V j

, то имеем
Вообще говоря, сумма в выражении (1.8) не обязана быть конечной. Можно переписать (1.8) в виде
а 2p-периодическая функция m 0

определяется следующим образом:
Во-вторых, ортогональность {
j(
x-
k)}
k

Î

Z

подразумевает, что
и, рассматривая сумму в (1.15) по четным и нечетным индексам, имеем
Используя 2p-периодичность функции m 0

и (1.14), после замены x/2
на x
, получаем необходимое условие
для коэффициентов h k

в (1.11). Заметив, что
и определив функцию y
следующим образом:
можно показать, что при каждом фиксированном масштабе j
Î
Z

вейвлеты
{
y
j,

k
(
x)=2 -

j/2

y(2 -

j

x-
k)}
k

Î

Z

образуют ортонормальный базис пространства W j

.
Равенство (1.17) определяет пару квадратурных зеркальных фильтров (quadrature mirror filters,
QMF)
H
и G
, где и . Коэффициенты QMF H
и G
вычисляются с помощью решения системы алгебраических уравнений. Число L
коэффициентов фильтра в (1.11) и (1.22) связано с числом исчезающих моментов М
, и всегда четно.
Выбранный фильтр Н
полностью определяет функции j
и y
и, таким образом, многомасштабный анализ. Кроме того, в правильно построенных алгоритмах значения функций j
и y
почти никогда не вычисляются. Благодаря рекурсивному определению вейвлетного базиса, все операции проводятся с квадратурными зеркальными фильтрами H
и G
, даже если в них используются величины, связанные с j
и y
.
Сжатие операторов или, другими словами, представление их в разреженном виде в ортонормированном базисе непосредственно влияет на скорость вычислительных алгоритмов.
Нестандартная форма оператора Т
с ядром K(
x,
y)
достигается вычислением следующих выражений:
4.1 Оператор
d

/

dx

в вейвлетном базисе

Нестандартные формы некоторых часто используемых операторов могут быть вычислены явно. Построим нестандартную форму оператора d/
dx
. Матричные элементы , , матриц , , и матрицы , где i,
l,
j
Î
Z

для оператора d/dx легко вычисляются как
Кроме того, используя (1.8) и (1.19), имеем
Таким образом представление d/
dx
полностью определяется величинами или, другими словами, отображением d/
dx
на подпространство V 0


.
Предложение 4.1

.
1. Если существует интеграл (4.11), тогда коэффициенты , l
Î
Z

в (5.8) удлвлетворяют следующей системе линейных алгебраических уравнений:
2. Если , тогда система (4.15)-(4.16) имеет единственное решение с конечным числом ненулевых , а именно с и .

Замечание.

Если М
= 1
, тогда система (4.15)-(4.16) имеет единственное решение, но интеграл (4.11) может не быть абсолютно сходящимся. Для базиса Хаара ( ) , мы получаем простейший конечный дифференциальный оператор .
Замечание 2.

Заметим, что выражения (4.12) и (4.13) для и ( ) могут быть упрощены с помощью смены порядка суммирования в (5.10) и (5.11) и введения коэффициентов корреляции , и . Выражение для особенно просто: .
Для доказательства Предложения 4.1 можно обратиться к [2].
Для решения системы (4.15)-(4.16) можно также воспользоваться итерационным алгоритмом. Начать можно с и , а дальше итерировать, используя (4.15) для вычисления .
4.2 Оператор
d

n


/

dx

n


в вейвлетном базисе

Так же как и для оператора d/
dx
, нестандартная форма оператора d n
/
dx n

полностью определяется своим отображением на подпространство V 0


, т.е. коэффициентами
Предложение 4.2.

1. Если интеграл в выражении (4.18) существует, тогда коэффициенты , l
Î
Z

удовлетворяют следующей системе линейных алгебраических уравнений
2. Пусть M≥ (
n+1)/2
, где М – число исчезающих моментов. Если интеграл в (4.18) существует, тогда система (4.19)-(4.20) имеет единственное решение с конечным числом нулевых коэффициентов , а именно для . Также для четных n

Замечание 3.

Если M≥ (
n+1)/2
, тогда решение линейной системы в Предложении 2 может существовать, когда интеграл в (4.18) не является абсолютно сходящимся.
Интегральные уравнения второго рода

Линейное интегральное уравнение Фредгольма есть выражение вида
где ядро , а неизвестная функция f(
x)
и функция в правой части , . Для простоты будем рассматривать интервал и введём следующее обозначение для всех и :
Предположим, что {
φ 1
,
φ 1
,…}
– ортонормальный базис для ; ядро представимо в этом базисе в следующем виде:
где коэффициенты K ij

вычисляются по формуле
Аналогично функции f
и g
представимы в виде
где коэффициенты f i

и g i

вычисляются по формулам:
Интегральное уравнение в этом случае соответствует бесконечной системе уравнений
Представление ядра может быть урезано до конечного числа слагаемых, что приводит к представлению интегрального оператора R
:
который аппроксимирует K
. Тогда интегральное уравнение аппроксимируется системой n
уравнений с n
неизвестными:
[LO_D,HI_D,LO_R,HI_R] = wfilters(wname);
a(2*k-1)=a(2*k-1)+2*LO_D(i+1)*LO_D(i+2*k);
R(l,abs(2*l-n))=-a(n)+R(l,abs(2*l-n));
R(l,abs(2*l+n))=-a(n)+R(l,abs(2*l+n));
function [al, bet, gam]=difcoef(wname,N)
[LO_D,HI_D,LO_R,HI_R] = wfilters(wname);
% вычисление значений α l,

β l,

γ l


if(((2*i+k-k1+L)0))
al(i+L)=HI_D(k-L)*HI_D(k1-L)*R(2*i+k-k1+L)+al(i+L);
bet(i+L)=HI_D(k-L)*LO_D(k1-L)*R(2*i+k-k1+L)+bet(i+L);
gam(i+L)=LO_D(k-L)*HI_D(k1-L)*R(2*i+k-k1+L)+gam(i+L);
1.
Вейвлет Добеши с
M=2.
a 1
=1.1250
a 3
=-0.1250
r 1
=-0.6667
r 2
=0.0833

2.
Вейвлет
Добеши
с
M=3.
a 1
=1.1719 a 3
=-0.1953 a 5
=0.0234

r 1
=-0.7452 r 2
=0.1452 r 3
=-0.0146 r 4
=-0.0003

3.
Вейвлет
Добеши
с
M=4.
a 1
=1.19628906249870 a 3
=-0.23925781249914 a 5
=0.04785156250041 a 7
=-0.00488281249997

r 1
=-0.79300950497055
r 2
=0.19199897079726
r 3
=-0.03358020705113

r 4
= 0.00222404967066
r 5
=0.00017220619000
r 6
=-0.00000084085054

4.
Вейвлет Добеши с
M=5.
a 1
=1.21124267578280
a 3
=
-0.26916503906311 a 5
=0.06921386718738

a 7
=-0.01235961914130 a 9
=0.00106811523422

r 1
=-0.82590601185686 r 2
=0.22882018706986 r 3
=-0.05335257193327

r 4
=0.00746139636621 r 5
=-0.00023923581985 r 6
=-0.00005404730164

r 7
=-0.00000025241171 r 8
=-0.00000000026960

a 1
=1.22133636474683 a 3
=-0.29079437255810 a 5
=0.08723831176674

a 7
=-0.02077102661228 a 9
=0.00323104858448 a 11
=-0.00024032592766

r 1
=-0.85013666156022 r 2
=0.25855294414318 r 3
=-0.07244058999853

r 4
=0.01454551104340 r 5
=-0.00158856154379 r 6
=0.00000429689148

r 7
=0.00001202657519 r 8
=0.00000042069120 r 9
=-0.00000000289967

6.
Вейвлет Койфмана с
M=2.
a 1
=1.20035616471068
a 3
=-0.24753371156550
a 5
=0.05401594511476

a 7
=-0.00724698442340
a 9
=0.00043220193586
a 11
=-0.00002361577240

r
1
=-0.80177838961957
r 2
=0.20214744976459
r 3
=-0.03943577686925

r 4
=0.00404789045961
r 5
=-0.00008445623632
r 6
=0.00000255044096

r 7
=0.00000088836508
r 8
=0.00000000237860
r 9
=-0.00000000002099

a 1
=1.12499999999971 a 3
=-0.12499999999971

r 1
=-0.66666666666616 r 2
=0.08333333333308

a 1
=1.17187500000666 a 3
=-0.19531250000432 a 5
=0.02343749999766 r 1
=-0.74520547946903 r 2
=0.14520547945865 r 3
=-0.01461187214494 r 4
=-0.00034246575336

a 1
=1.19628906249990 a 3
=-0.23925781249985 a 5
=0.04785156249993

r 1
=-0.79300950497424 r 2
=0.19199897079876 r 3
=-0.03358020705098

r 4
=0.00222404967071 r 5
=0.00017220619000 r 6
=-0.00000084085054

1. Beylkin G. Wavelets and Fast Numerical Algorithms
2. Beylkin G. Wavelets, Multiresolution Analysis and Fast Numerical Algorithms
3. Beylkin G. In The Representation.of Operators in Bases of Compactly Supported Wavelets
4. Bradley K. Alpert A Class of Bases in L 2
for the Sparse Representation of Integral Operators
5. Дремин И.М., Иванов О.В., Нечитайло В.А. Вейвлеты и их использование // Успехи физических наук – 2001, №5. – С.465-500

Название: Операторы в вейвлетном базисе
Раздел: Рефераты по математике
Тип: реферат
Добавлен 04:49:57 08 октября 2005 Похожие работы
Просмотров: 44
Комментариев: 18
Оценило: 4 человек
Средний балл: 5
Оценка: неизвестно   Скачать

Срочная помощь учащимся в написании различных работ. Бесплатные корректировки! Круглосуточная поддержка! Узнай стоимость твоей работы на сайте 64362.ru
Если Вам нужна помощь с учебными работами, ну или будет нужна в будущем (курсовая, дипломная, отчет по практике, контрольная, РГР, решение задач, онлайн-помощь на экзамене или "любая другая" учебная работа...) - обращайтесь: https://clck.ru/P8YFs - (просто скопируйте этот адрес и вставьте в браузер) Сделаем все качественно и в самые короткие сроки + бесплатные доработки до самой сдачи/защиты! Предоставим все необходимые гарантии.
Привет студентам) если возникают трудности с любой работой (от реферата и контрольных до диплома), можете обратиться на FAST-REFERAT.RU , я там обычно заказываю, все качественно и в срок) в любом случае попробуйте, за спрос денег не берут)
Да, но только в случае крайней необходимости.

Реферат: Операторы в вейвлетном базисе
Реферат: Незаконне виготовлення підроблення використання або збут незаконно виготовлених одержаних чи
Реферат: Проблемы энергосбережения в системах электроснабжения
Рефераты На Тему Пушкина
Курсовая работа по теме Методы подбора персонала
Ескендір Эссе
Реферат: Внутрішній аудит витрат підприємства громадського харчування
Реферат: Essay About Criticism Of Shakespeare
Контрольная работа по теме Санитарно-эпидемиологические требования к отходам
Реферат: Отче по преддипломной практике в СЧ СУ при УВД по г. Ростову-на-Дону
Реферат по теме Взвод в наступлении с ходу
Реферат: Тюркское влияние на этническую динамику русского народа. Скачать бесплатно и без регистрации
Реферат: Электронные системы индивидуального управления отбора рабочих органов узорообразования
Дипломная работа по теме Эксперимент в социальных организациях
Дипломная работа по теме Организация методической работы и её роль в развитии учреждения социального обслуживания
Дневник Практики Помощник Врача Гигиениста Заполненный
Дипломная работа по теме Социальные проблемы современной сельской молодежи
Реферат по теме Неоклассический синтез Самуэльсона
Реферат: Подготовка населения в области защиты населения от чрезвычайных ситуаций. Скачать бесплатно и без регистрации
Заключительный Лист Реферата
Контрольная Работа Мордкович 10 Класс Ответы
Доклад: Медицина Древнего Египта
Реферат: Административные учреждения штатов и местное самоуправление в США
Реферат: Особенности права стран древнего востока