Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Met2009

.pdf
Скачиваний:
11
Добавлен:
14.05.2015
Размер:
479.23 Кб
Скачать

81

ееестьматематическогосреднее квадратичноеожиданияотклонение случайной величины t от личины. Чем меньше дисперсия,t итемназываетсякучнее к дисперсиейсреднемузначениюэтой ве-

располагаютсядисперсии значения случайной величины t. Поэтому, величина

функции t может служить характеристикой величины носителя

ψ.

Пример. Пусть

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

0,

 

t 6 [a, b]

 

 

 

 

 

ψ(t) =

 

b

a , t [a, b]

Тогда

 

= a+b

2

=

 

(b−a)2

 

 

 

 

 

 

t

 

 

 

 

 

 

2 ,

t

 

12 .

 

 

 

 

 

Теорема .3 (Соотношение неопределенности Гейзенберга).

Справедливо неравенство Гейзенберга

1

x ω

Доказательствоженных оператора. Определим в пространстве L2(R) два самосопря-

(Aψ)(t) = t

·

ψ(t), B = F−1AF.

 

 

Заметим, что Так как

(Bψ)(t) = 21πi dtd ψ(t).

ω

+∞ e−2πitωψ(t)dt =

−1

+∞

d

e−2πitωψ(t)dt =

2πi

 

 

 

Z−∞

+∞

Z−∞ dt

 

 

 

 

1

1

 

 

 

 

 

 

Z−∞ e−2πitω

 

dt =

 

F

 

 

.

 

 

2πi

dt

2πi

dt

Òî

tkψk2 = htψ, ψi = hAψ, ψi,

k ˆk2 h ˆ ˆi h i

ω ψ = ωψ, ψ = Bψ, ψ .

хоть произведение
a по частоте. Таким образом
ω

82 ГЛАВА . ФУРЬЕ И ВЕЙВЛЕТ ПРЕОБРАЗОВАНИЕ Справедливо равенство:

 

 

 

 

 

 

 

 

 

1 d

 

 

 

 

1 d

 

1

 

 

 

 

 

(AB − BA)ψ = t

 

 

 

 

 

 

ψ −

 

 

 

 

(tψ) = −

 

· ψ.

 

 

 

2πi

dt

2πi

dt

2πi

 

женнымиОператоры.ПоложимA è B в этой ситуации называют канонически сопря-

ственный оператор.

A1 =

 

 

 

 

B1 = ωE −1B

 

 

 

E

 

tE − A

 

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

, ãäå

 

тожде-

 

 

 

 

 

 

 

 

 

(A1B1 − B1A1)ψ = −

 

· ψ. Заметим, что

 

 

 

 

 

 

ˆ 2

2πi

tkψk

2

2

,

 

 

= kB1ψk

2. Например второе равенство

 

= kA1ψk

ωkψk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−B1ψk2 = h(

ωE − B) ψ, (

ωE − B) ψi =

 

 

 

 

= hF −1 (

ωE − A) Fψ, F −1 (

ωE − A) Fψi =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= h(ωE − A) ψ, (ωE − A) ψi =

 

 

 

 

 

Z +∞

|− |2| ˆ |2

=ω ω ψ(ω) dω

−∞

Следовательно

21π kψk2 = |h(A1B1 − B1A1)ψ, ψi| = |h(A1B1ψ, ψi − h(B1A1)ψ, ψi| = = |hB1ψ, A1ψi − hA1ψ, B1ψi| = |2 ImhB1ψ, A1ψi| ≤

p

≤ 2kB1ψkkA1ψk = 2kψk2 t ω.

Замечание. Соотношение неопределенности Гейзенберга означа- ет, что нельзя одновременно локализовать сигнал на временной оси и на оси частот. Основная идея вейвлет-анализа состоит в том, чтобы вести частотно-временной анализ сигнала с переменным разрешением ta по времени и

t ω будет постоянно, тем не менее получится более полная картина сигнала.

4. Непрерывное вейвлет преобразование

Определение .6. Пусть ψ L1(R) ∩ L2(R), функция ψ называется базовым вейвлетом если выполнено условие

2

Zψb(ω)

cψ = −∞ |ω| dω < ∞

Найдем преобразование Фурье от

83

ãäå ψb преобразование фурье функции ψ. Положим

p

ψθ,b(t) = |θ|ψ (θ (t − b)) , θ, b R.

ψθ,b(t)

 

Z

 

 

 

 

 

 

 

 

 

ψθ,b(ω) =

−∞ e−2πiωtp

|θ|

ψ (θ (t − b)) dt =

d

 

 

 

2πiω

u +b

 

1

 

 

 

 

 

 

=

 

|θ| Z−∞ e

(

θ

)

 

ψ (u) du =

 

θ

=

p

ω

1|θ|

 

 

 

e−2πiωbψ

θ

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

Определение .7. Непрерывным

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вейвлет преобразованием назы-

вается интегральное преобразование

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wf(θ, b) = Z−∞ f(t)

 

 

 

dt = p

 

Z−∞ f(t)

 

 

 

 

dt

 

 

 

 

 

 

ψθ,b(t)

|θ|

ψ (θ (t − b))

 

или из равенства Планшереля

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wf(θ, b) =

 

Z−∞ f(ω)ψθ,b (ω)dω =

 

 

 

 

 

 

 

 

 

1

b d

 

 

 

 

 

 

 

 

 

 

2πiωb

 

 

ω

 

 

=

 

 

 

 

Z−∞ f(ω)e

 

 

ψ

 

 

 

 

 

 

 

 

 

 

θ

 

 

 

 

|θ|

 

 

 

Теорема .4. Пусть

ψ

базовый

 

b

 

 

b

 

 

 

 

 

 

f, g

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вейвлет. Тогда для любых

 

 

 

 

 

 

 

 

 

L2(R)

Z Z

Wf(θ, b)Wg(θ, b)dθdb = cψ hf, gi

RR

Для любой точки t, в которой функция f непрерывна

 

f(t) = cψ ZR ZR Wf(θ, b)ψθ,b(t)dθdb

 

1

 

 

 

 

 

 

 

 

 

В случае вещественной функции ψ

 

 

 

 

 

 

 

 

 

cψ

 

 

 

 

 

Z0

ZR Wf(θ, b)Wg(θ, b)db =

 

hf, gi

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(t) =

 

Z0

ZR Wf(θ, b)ψθ,b(t)dθdb

cψ

84 ГЛАВА . ФУРЬЕ И ВЕЙВЛЕТ ПРЕОБРАЗОВАНИЕ Доказательство. В силу тождества Планшерелÿ

 

 

 

1

 

 

ω

 

W f(θ, b) =

 

 

 

 

 

Z−∞ f(ω)e2πiωbψ

 

 

dω = Fθ(b)

 

 

 

 

 

θ

 

 

 

 

|θ|

 

 

p

b

b

 

 

c

 

 

 

 

 

 

 

ãäå Fθ(ω) = f(ω)ψ

ω 1

 

 

θ

|θ| поэтому

 

 

 

 

b b

 

 

 

 

 

 

 

 

 

 

 

Z Z

Wf(θ, b)Wg(θ, b)dθdb

RR

ZZ

= dθ Fcθ(b)Gcθ(b)db =

RR

ZZ

 

 

=

 

 

 

 

dθ Fθ(ω)Gθ(ω)dω =

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

R

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

g(ω)ψ θ dω =

 

 

ZR

 

|θ|ZR f(ω)ψ θ

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

b

b

 

 

ω

 

 

 

 

 

 

 

b

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ

 

b

 

 

 

 

 

 

=

 

 

 

 

f(ω)

g(ω)

 

 

 

 

 

 

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

R

R

 

 

 

 

 

 

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

b

 

 

 

 

 

Z

 

 

 

b

|

 

 

|

 

 

 

 

 

 

 

 

 

 

 

òàê êàê

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ

ωθ

 

 

2

 

 

 

 

 

 

 

ψ (θ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

dθ = cψ

 

 

 

 

 

R

 

 

 

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b|

θ

|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

b

|

 

|

 

 

 

 

 

 

Z−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если функция ψ принимает действительные значения, то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ (ω) =

Z−∞ e−2πitωψ(t)dt,

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

2πitω

ψ(t)dt = ψ (−ω) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ (ω) =

Z−∞ e

 

 

 

 

 

 

 

 

 

 

 

 

 

b(ω)

 

 

 

=

 

 

 

ψ (−ω)

,

 

 

2

 

 

 

 

 

 

b

 

 

 

 

 

 

ψ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ

θ

 

 

 

 

 

 

 

 

 

b

 

ψ (θ)

 

 

 

 

 

cψ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

=

,

 

ω = 0

0

 

 

 

θ

 

 

 

 

 

 

 

 

 

 

0

 

b

θ

 

 

 

 

 

 

Z

 

 

b

|

 

|

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда следует первое утверждение. Положим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

t / [t0

− ε, t0

+ ε]

 

 

 

 

 

 

gε (t) =

 

 

 

t

[t0

 

ε, t0

+ ε]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

85

Тогда

 

0 h

 

 

εi

 

 

 

 

 

 

 

 

 

 

 

 

ε

 

f, g

=

0

),

 

 

 

 

 

 

 

 

 

lim

 

f(t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim W g

 

(θ, b)

=

lim

g

, ψ

θ,bi

=

ψ

 

(t

)

 

ε→0

 

ε

 

 

 

ε→0 h

ε

 

 

 

θ,b

0

 

 

подставляя в первое равенство получим второе утверждение.

Замечание. Можно показать, что функция F (a, b) есть преобра-

зование вейвлета от некоторой функции

 

f(t), в том и только в

том случае когда выполнено тождество

 

 

 

 

 

 

 

 

F (a0, b0) = ZR Z0

F (a, b)K (a, b; a0, b0)

a

 

 

 

 

 

 

 

 

 

 

 

dadb

 

ãäå

 

 

ZR a

a

a0

a0

0

0

 

K (a, b; a

, b

) =

 

1

ψ

t − b

 

1 ψ

b0

− t

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определение .8. Спектральной плотностью вейвлета назовем функцию

 

 

Eω(θ) = ψθ,b(ω)

2

= ψ

ω

 

2

1

 

 

 

 

 

θ

 

 

θ

 

 

Частотным коэффициентом

вейвлета назовем

 

число kψ определя-

емое из равенства

 

d

 

 

b

 

 

 

 

 

 

 

 

 

 

 

θ = kψω = 0θEω(θ)dθ = ω 0ψ (θ)2

2

 

θ2 .

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

Eω(θ)dθ

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

ψb(θ)

 

θ

 

 

 

R

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

Следовательно частотный коэффициент вейвлета равен

 

 

 

2

kψ = R0

ψ (θ)

 

2

.

 

 

 

 

 

θ2

 

 

ψ (θ)

 

θ

 

 

0

b

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

Коэффициентом частотной дисперсии вейвлета назовем число dψ определяемое из равенства

R

 

 

0

 

 

 

2

 

 

 

dψω2 =

0

 

θ

− θ Eω(θ)dθ

 

R

Eω(θ)dθ

 

 

 

 

 

 

 

 

86

ГЛАВА . ФУРЬЕ И ВЕЙВЛЕТ ПРЕОБРАЗОВАНИЕ

 

dψω2 = ω2

0ψ (θ)

2

θ3

 

0

ψ (θ)

2

θ

0ψ (θ)

2

θ2

2

 

 

.

 

 

 

R

b

 

 

 

R

b

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 R2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

ψ (θ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно коэффициент

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

частотной дисперсии вейвлета равен

 

dψ =

0

ψ (θ)

θ3

0ψ (θ)

2

θ

0

ψ (θ)

θ2

2

 

 

.

 

 

 

R

b

 

2

 

R

 

b

 

 

 

 

R2

b

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

ψ (θ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6. Дискретное Фурье и вейвлет преобразование

Определение .9. Пусть x = (x1, x2, . . . , xN ) вектор длины N. Дискретным преобразованием вектора x, называют вектор:

N

X

−i2π(k−1)(n−1)

, 1

 

k

 

N.

X(k) = x(n)e

N

 

 

 

 

 

n=1

Обратное преобразование определено формулой

N

x(n) =

1 X

N

k=1

X(k)e

−i2π(k−1)(n−1)

, 1 ≤ n ≤ N.

N

 

 

Прямое вычисление дискретного преобразования Фурье требует выполнения O(N2) арифметических операций. Если N =

n1n2 . .,. nm

 

 

 

 

требующий, то можно применить алгоритм быстрого преобразования

имеет порядок

1

2

+

· · ·m. Тогда сложность вычислений

öèéФурье. Он особенно удобен,O(N(nåñëè+ n

 

+ nm)) арифметических опера-

N = 2

O(N ln N).

Литература

[1]Канторович Л.В., Акилов Г.П. Функциональный анализ. М.: Наука, 1984. 752 с.

[2]Кириллов А.А., Гвишиани А.Д. Теоремы и задачи функционального анализ. М.: Наука, 1979. 384 с.

[3]Лаврентьев М.М., Савельев Л.Я. Теория операторов и некорректные задачи. Новосибирск: Изд-во Ин-та математики, 1999. - 702 с.

[4]Колмогоров А.М., Фомин В.Ф. Элементы функционального анализа и теории функций. М.: Наука, 1976. 544 с.

[5]Люстерник Л.А., Соболев В.В. Элементы функционального анализа. М.: Наука, 1982. 264 с.

[6]Треногин В.А., Писаревский Б.М., Соболева Т.С. Задачи и упражнения по функциональному анализу. М.: Наука, 1984.256 с.

[7]Треногин В.А. Функциональный анализ. М.: Наука, 1980. 496 с.

[8]Васильева А.Б., Тихонов М.А. Интегральные уравнения. М.: Èçä-âî ÌÃÓ, 1989. 156 ñ.

87

Саженков Александр Николаевич Вайгант Владимир Андреевич Матукевич Ольга Юрьевна Саженкова Татьяна Владимировна Славский Виктор Владимирович

ФУНКЦИОНАЛЬНЫЙ АНАЛИЗ практикум

ËÐ 020261 Í/Ê

Компьютерный набор О.А. Жданова Компьютерная верстка Н.С. Морозова Оригинал-макет представлен в авторской редакции

ПодписаноОфсетная печатьв печать.

УчФормат.-изд. 60ë. ×4.847./16.

Заказ

Тираж 100 экз.

Издательство Алтайского университета: 656099, Барнаул-99, ул. Димитрова, 66.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]