[62fbef]: wavelets / toshow_wfiltComp.m  Maximize  Restore  History

Download this file

76 lines (55 with data), 1.3 kB

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
function toshow_wfiltComp
load vonkoch;
f=vonkoch';
clear vonkoch;
%f = traindoppler;
figure(1);
clf;
plot(f);
scalegramBounds = [-200 0];
figure(3);
clf;
waitS =3;
J = 5;
w = waveletfb('db',5);
c = fwt(f,w,J);
figure(2);clf;plotwavc(c,'clim',scalegramBounds);
figure(3);freqzfb(fftshift(multid(w.h,J)));
pause(waitS);
w = waveletfb('optfs',2);
c = fwt(f,w,J);
figure(2);clf;plotwavc(c,'clim',scalegramBounds);
figure(3);freqzfb(multid(w.h,J));
pause(waitS);
w = waveletfb('dden',6);
c = fwt(f,w,J);
figure(2);clf;plotwavc(c,'clim',scalegramBounds);
figure(3);freqzfb(multid(w.h,J));
pause(waitS);
w = waveletfb('hden',3);
c = fwt(f,w,J);
figure(2);clf;plotwavc(c,'clim',scalegramBounds);
figure(3);freqzfb(multid(w.h,J));
pause(waitS);
w = waveletfb('dgrid',3);
c = fwt(f,w,J);
figure(2);clf;plotwavc(c,'clim',scalegramBounds);
figure(3);freqzfb(multid(w.h,J));
pause(waitS);
w = waveletfb('symds',3);
c = fwt(f,w,J);
figure(2);clf;plotwavc(c,'clim',scalegramBounds);
figure(3);freqzfb(multid(w.h,J));
pause(waitS);
w = waveletfb('dtree');
c = fwt(f,w,J);
figure(2);clf;plotwavc(c,'clim',scalegramBounds);
hdt = cell(2,J);
hdt{1} = w.h{1};
hdt{2} = w.h{2};
for jj = 2:J
hdt{1,jj} = w.h{3};
hdt{2,jj} = w.h{4};
end
figure(3);freqzfb(multid(hdt,J));
pause(waitS);