кучка / Глазунова_диагностические_методы_1
.pdf
|
81 |
|
|
|
|
|
|
WaveletTransformerDiskr |
|
ColShag_a 30 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
ColShag_b 40 |
|
|
|
|
|
|
|
a0 1.2 |
|
|
|
|
|
|
|
dt 0.1 |
|
|
|
|
|
|
|
for j 0 ColShag_a |
1 |
|
|
|
|
|
|
for k 0 ColShag_b |
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
Cj k |
|
j |
j |
|
|
|
|
a0 |
|
a0 (k 20) |
t dt |
s(t dt ) dt |
|
|
|
t 1000 |
|
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
|
|
|
Вейвлетная плоскость сигнала |
|
|
|
|
|
WaveletTransformerDiskrT
Обратное вейвлет-преобразование
t 20 19.9 20 a0 1.2
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
82 |
|
|
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
i t |
|
|
|
|
|
|||||
|
|
|
|
0(t) e |
dt |
|
|
|
|
||||||||||||||
F |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
a029 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
||
C |
|
|
|
|
|
|
|
2 |
|
|
|
|
d |
|
|
|
|
||||||
0 |
|
F |
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
1 |
|
29 |
|
|
1 |
|
|
|
|
40 |
|
|
|
|
a0j |
a0j (b 20) t |
|
||
x(t) |
|
|
|
|
|
|
|
|
|
|
|
|
WaveletTransformerDiskr |
|
|
||||||||
|
|
|
|
C |
|
|
|
|
2 |
|
|
|
|
j b |
|
|
|
||||||
|
|
|
|
|
|
|
j |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
j 0 |
|
|
|
|
b 0 |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
(a0) |
|
|
|
|
|
|
|
||||||||||
Восстановленный сигнал |
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
s(t) |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
2 20 |
|
|
|
|
|
|
|
|
|
|
10 |
0 |
|
10 |
20 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
П 2.1.3. Вейвлет-преобразование MHAT ("мексиканская шляпа")
|
d |
2 |
|
|
t |
2 |
|
|
|
|
|
|
|
|
|
||
_mhat(t) |
dt2 |
exp |
2 |
|
|
0(t) _mhat(t)
|
1 |
|
|
|
t b |
|||
(a b t) |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|||
|
|
a |
|
|
|
a |
|
83 |
|
|
|
|
|
|
Прямое вейвлет-преобразование |
|
|
|
|
|
||
WaveletTransformerDiskr |
|
ColShag_a 30 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
ColShag_b 40 |
|
|
|
|
|
|
|
a0 1.2 |
|
|
|
|
|
|
|
dt 0.1 |
|
|
|
|
|
|
|
for j 0 ColShag_a |
1 |
|
|
|
|
|
|
for k 0 ColShag_b |
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
Cj k |
|
j |
j |
|
|
|
|
a0 |
|
a0 (k 20) |
t dt |
s(t dt ) dt |
|
|
|
t 1000 |
|
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
|
|
Вейвлетная плоскость сигнала
WaveletTransformerDiskrT
Обратное вейвлет-преобразование
t 20 19.9 20 a0 1.2
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
84 |
|
|
|
|
|
|
|
inf |
|
|
|
|
|
|
i t |
|
|
|
|
|
|
|||||
|
|
|
0(t) |
e |
dt |
|
|
|
|
||||||||||||
F |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
inf |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a029 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|||
C |
|
|
|
|
|
|
|
2 |
|
|
|
d |
|
|
|
|
|||||
0 |
|
F |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
1 |
|
29 |
|
|
1 |
|
|
|
|
|
40 |
|
|
|
a0j |
a0j (b 20) t |
|
||
x(t) |
|
|
|
|
|
|
|
|
|
|
WaveletTransformerDiskr |
|
|
||||||||
|
|
C |
|
|
|
|
2 |
|
|
|
|
j b |
|
|
|
||||||
|
|
|
|
|
j |
|
|
|
|
|
|||||||||||
|
|
|
|
j 0 |
|
|
|
|
|
|
b 0 |
|
|
|
|
|
|
||||
|
|
|
|
(a0) |
|
|
|
|
|
|
|
|
|
||||||||
Восстановленный сигнал |
|
|
|
|
|
||||||||||||||||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s(t) |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
2 20 |
|
|
|
|
|
|
|
|
|
|
10 |
|
0 |
10 |
20 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
П 2.1.4. Вейвлет-преобразование встроенными функциями |
Mathcad
Прямое вейвлет-преобразование
n 0 2 256 1 t 20 19.9 20
S s n 255
n 10
85
WaveletTransformerDiskr wave(S)
Вейвлетная плоскость сигнала
WaveletTransformerDiskr
Обратное вейвлет-преобразование
X iwave(WaveletTransformerDiskr )
Восстановленный сигнал
1
X n
0
s(t)
1
20 |
15 |
10 |
5 |
0 |
|
|
5 |
10 |
15 |
20 |
|
|
|
n 255 |
|
t |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
86
П 2.2. Вейвлет-преобразование модулированного сигнала
s(t) sin [0.2 (t)]2
t 20 19.9 20
|
1 |
|
|
|
|
|
|
|
|
s(t) |
0 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
20 |
15 |
10 |
5 |
0 |
5 |
10 |
15 |
20 |
|
|
|
|
|
t |
|
|
|
|
П 2.2.1. Вейвлет-преобразование Хаара
_Haar(t) |
|
|
1 |
|
if |
0 t |
1 |
|
|||||||
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|||
|
|
|
|
|
( 1) |
if |
1 |
t 1 |
|||||||
|
|
|
|
|
|
||||||||||
|
|
|
|
2 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
0 |
|
otherwise |
|||||||||
|
|
|
|
|
|||||||||||
0(t) _Haar(t) |
|
|
|
|
|
|
|
|
|||||||
|
|
1 |
|
|
|
|
t b |
||||||||
(a b t) |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
a |
|
|
|
a |
|
|
Прямое вейвлет-преобразование
|
87 |
|
|
|
|
|
|
WaveletTransformerDiskr |
|
ColShag_a 30 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
ColShag_b 40 |
|
|
|
|
|
|
|
a0 1.2 |
|
|
|
|
|
|
|
dt 0.1 |
|
|
|
|
|
|
|
for j 0 ColShag_a |
1 |
|
|
|
|
|
|
for k 0 ColShag_b |
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
Cj k |
|
j |
j |
|
|
|
|
a0 |
|
a0 (k 20) |
t dt |
s(t dt ) dt |
|
|
|
t 1000 |
|
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
|
|
|
Вейвлетная плоскость сигнала |
|
|
|
|
|
WaveletTransformerDiskrT
Обратное вейвлет-преобразование
t 20 19.9 20
a0 1.2
88
C 1
x(t) 1 |
|
29 |
|
|
|
1 |
|
|
40 |
|
|
|
|
|
|
a0j |
a0j (b 20) t |
|
|
|
|
|
|
|
|
|
|
|
|
WaveletTransformerDiskr |
|
||||||||
C |
|
|
|
|
2 |
|
|
|
|
|
j b |
|
|
|
|||||
|
|
|
j |
|
|
|
|
|
|
||||||||||
|
|
j 0 |
|
|
|
|
b 0 |
|
|
|
|
|
|
|
|
||||
|
|
(a0) |
|
|
|
|
|
|
|
|
|
||||||||
Восстановленный сигнал |
|
|
|
|
|
|
|||||||||||||
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s(t) |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 20 |
|
|
|
|
|
|
|
10 |
|
0 |
|
10 |
20 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
П 2.2.2. Вейвлет-преобразование WAVE |
|
|
|
|
|||||||||||||||
|
|
|
d |
|
|
|
|
|
t |
2 |
|
|
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
_wave(t) dt exp |
2 |
|
|
t exp |
2 t |
|
|
|
|
|
|||||||||
0(t) _wave(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
(a b t) |
1 |
|
|
|
|
t b |
|
|
|
|
|
|
|||||||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
a |
|
|
|
|
a |
|
|
|
|
|
|
|
|
Прямое вейвлет-преобразование
|
|
89 |
|
|
|
|
|
|
|
|
|
|
|
|
|
WaveletTransformerDiskr |
ColShag_a |
30 |
|
|
|
|
|
|
ColShag_b 40 |
|
|
|
|
|
|
|
a0 1.2 |
|
|
|
|
|
|
|
dt 0.1 |
|
|
|
|
|
|
|
for j 0 ColShag_a |
1 |
|
|
|
|
|
|
for k 0 ColShag_b |
|
|
|
|
||
|
|
1000 |
|
|
|
|
|
|
Cj k |
|
|
j |
j |
|
|
|
a0 |
|
a0 (k 20) |
t dt |
s(t dt ) dt |
||
|
|
t 1000 |
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вейвлетная плоскость сигнала |
|
|
|
|
|
WaveletTransformerDiskrT
Обратное вейвлет-преобразование
t 20 19.9 20
a0 1.2
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
90 |
|
|
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
|
i t |
|
|
|
|
|
|||||
|
|
|
|
|
0(t) e |
dt |
|
|
|
|
||||||||||||||
F |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
a029 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
||||
C |
|
|
|
|
|
|
|
|
|
2 |
|
|
d |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|||||||||||||||||
0 |
|
F |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
29 |
|
|
|
|
|
|
|
40 |
|
|
|
|
a0j |
a0j (b 20) t |
|
|||||
|
|
|
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|||||||||
x(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
WaveletTransformerDiskr |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
C |
|
a0j |
2 |
|
|
|
|
j b |
|
|
|
|||||||||
|
|
|
|
|
|
|
j 0 |
|
|
b 0 |
|
|
|
|
|
|
Восстановленный сигнал
2
x(t)
0
s(t)
2 |
|
|
|
|
|
|
|
|
|
20 |
10 |
0 |
10 |
20 |
t
П 2.2.3. Вейвлет-преобразование MHAT ("мексиканская шляпа")
|
d |
2 |
|
|
t |
2 |
|
1 |
|
|
|
|
|
|
|||
_mhat(t) |
dt2 |
exp |
2 |
|
exp 2 |
0(t) _mhat(t)
|
1 |
|
|
|
t b |
|||
(a b t) |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|||
|
|
a |
|
|
|
a |
t2 t2 exp 1 t2
2