| Wavelet Toolbox | |
Syntax
wavemngr('add',FN,FSN,WT,NUMS,FILE)
wavemngr('add',FN,FSN,WT,NUMS,FILE,B)
wavemngr('del',N)
wavemngr('restore')
wavemngr('restore',IN2)
OUT1 = wavemngr('read')
OUT1 = wavemngr('read',IN2)
OUT1 = wavemngr('read_asc')
Description
wavemngr is a type of wavelets manager. It allows you to add, delete, restore or read wavelets.
wavemngr('add',FN,FSN,WT,NUMS,FILE) or wavemngr('add',FN,FSN,WT,NUMS,FILE,B) or wavemngr('add',FN,FSN,WT,{NUMS,TYPNUMS},FILE) or wavemngr('add',FN,FSN,WT,{NUMS,TYPNUMS},FILE,B), add a new wavelet family to the toolbox.
FSN = Family Short Name (string of length equal or less than four characters)
WT = 1, for orthogonal wavelets
WT = 2, for biorthogonal wavelets
WT = 3, for wavelet with scaling function
WT = 4, for wavelet without scaling function
WT = 5, for complex wavelet without scaling function
If the family contains a single wavelet, NUMS = ' '.
|
|
|
NUMS is a string containing a blank separated list of items representing wavelet parameters.
|
: NUMS = '1.1 1.3 ... 4.4 5.5 6.8' |
NUMS is a string containing a blank separated list of items representing wavelet parameters, terminated by the special sequence **.
|
: NUMS = '1 2 3 4 5 6 7 8 9 10 **' |
shan |
: NUMS = '1-1.5 1-1 1-0.5 1-0.1 2-3 **' |
TYPNUMS specifies the wavelet parameter input format: 'integer' or 'real' or 'string'; the default value is 'integer'.
db |
: TYPNUMS = 'integer' |
bior |
: TYPNUMS = 'real' |
shan |
: TYPNUMS = 'string' |
FILE = MAT-file or M-file name (string). See usage in the "Examples" section.
B = [lb ub] specifies lower and upper bounds of effective support for wavelets of
type = 3, 4 or 5.
This option is fully documented in Chapter 7 of the User's Guide, "Adding Your Own Wavelets".
wavemngr('del',N), deletes a wavelet or a wavelet family. N is the Family Short Name or the Wavelet Name (in the family). N is a string.
wavemngr('restore') or wavemngr('restore',IN2), restore previous or initial wavelets. If nargin = 1, the previous wavelets.asc ASCII-file is restored; otherwise the initial wavelets.asc ASCII-file is restored. Here IN2 is a dummy argument.
OUT1 = wavemngr('read') returns all wavelet family names.
OUT1 = wavemngr('read',IN2) returns all wavelet names, IN2 is a dummy argument.
OUT1 = wavemngr('read_asc') reads wavelets.asc ASCII-file and returns all wavelets information.
Examples
% List initial wavelets families.
wavemngr('read')
ans =
===================================
Haar haar
Daubechies db
Symlets sym
Coiflets coif
BiorSplines bior
ReverseBior rbio
Meyer meyr
DMeyer dmey
Gaussian gaus
Mexican_hat mexh
Morlet morl
Complex Gaussian cgau
Shannon shan
Frequency B-Spline fbsp
Complex Morlet cmor
===================================
% List all wavelets.
wavemngr('read',1)
ans =
===================================
Haar haar
===================================
Daubechies db
------------------------------
db1 db2 db3 db4
db5 db6 db7 db8
db9 db10 db**
===================================
Symlets sym
------------------------------
sym2 sym3 sym4 sym5
sym6 sym7 sym8 sym**
===================================
Coiflets coif
------------------------------
coif1 coif2 coif3 coif4
coif5
===================================
BiorSplines bior
------------------------------
bior1.1 bior1.3 bior1.5 bior2.2
bior2.4 bior2.6 bior2.8 bior3.1
bior3.3 bior3.5 bior3.7 bior3.9
bior4.4 bior5.5 bior6.8
===================================
ReverseBior rbio
------------------------------
rbio1.1 rbio1.3 rbio1.5 rbio2.2
rbio2.4 rbio2.6 rbio2.8 rbio3.1
rbio3.3 rbio3.5 rbio3.7 rbio3.9
rbio4.4 rbio5.5 rbio6.8
===================================
Meyer meyr
===================================
DMeyer dmey
===================================
Gaussian gaus
------------------------------
gaus1 gaus2 gaus3 gaus4
gaus5 gaus6 gaus7 gaus8
gaus**
===================================
Mexican_hat mexh
===================================
Morlet morl
===================================
Complex Gaussian cgau
------------------------------
cgau1 cgau2 cgau3 cgau4
cgau5 cgau**
===================================
Shannon shan
------------------------------
shan1-1.5 shan1-1 shan1-0.5 shan1-0.1
shan2-3 shan**
===================================
Frequency B-Spline fbsp
------------------------------
fbsp1-1-1.5 fbsp1-1-1 fbsp1-1-0.5 fbsp2-1-1
fbsp2-1-0.5 fbsp2-1-0.1 fbsp**
===================================
Complex Morlet cmor
------------------------------
cmor1-1.5 cmor1-1 cmor1-0.5 cmor1-1
cmor1-0.5 cmor1-0.1 cmor**
===================================
In the following example, new compactly supported orthogonal wavelets are added to the toolbox. These wavelets, which are a slight generalization of the Daubechies wavelets, are based on the use of Bernstein polynomials and are due to Kateb and Lemarié in an unpublished work.
% Add new family of orthogonal wavelets.
% You must define:
%
% Family Name: Lemarie
% Family Short Name: lem
% Type of wavelet: 1 (orth)
% Wavelets numbers: 1 2 3 4 5
% File driver: lemwavf
%
% The function lemwavf.m must be as follows:
% function w = lemwavf(wname)
% where the input argument wname is a string:
% wname = 'lem1' or 'lem2' ... i.e.,
% wname = sh.name + number
% and w the corresponding scaling filter.
% The addition is obtained using:
wavemngr('add','Lemarie','lem',1,'1 2 3 4 5','lemwavf');
% The ascii file 'wavelets.asc' is saved as
% 'wavelets.prv', then it is modified and
% the MAT file 'wavelets.inf' is generated.
% List wavelets families.
wavemngr('read')
ans =
===================================
Haar haar
Daubechies db
Symlets sym
Coiflets coif
BiorSplines bior
ReverseBior rbio
Meyer meyr
DMeyer dmey
Gaussian gaus
Mexican_hat mexh
Morlet morl
Complex Gaussian cgau
Shannon shan
Frequency B-Spline fbsp
Complex Morlet cmor
Lemarie lem
===================================
% Remove the added family.
wavemngr('del','Lemarie');
% List wavelets families.
wavemngr('read')
ans =
===================================
Haar haar
Daubechies db
Symlets sym
Coiflets coif
BiorSplines bior
ReverseBior rbio
Meyer meyr
DMeyer dmey
Gaussian gaus
Mexican_hat mexh
Morlet morl
Complex Gaussian cgau
Shannon shan
Frequency B-Spline fbsp
Complex Morlet cmor
===================================
% Restore the previous ascii file
% 'wavelets.prv', then build
% the MAT-file 'wavelets.inf'.
wavemngr('restore');
% List restored wavelets.
wavemngr('read',1)
ans =
===================================
Haar haar
===================================
Daubechies db
------------------------------
db1 db2 db3 db4
db5 db6 db7 db8
db9 db10 db**
===================================
Symlets sym
------------------------------
sym2 sym3 sym4 sym5
sym6 sym7 sym8 sym**
===================================
Coiflets coif
------------------------------
coif1 coif2 coif3 coif4
coif5
===================================
BiorSplines bior
------------------------------
bior1.1 bior1.3 bior1.5 bior2.2
bior2.4 bior2.6 bior2.8 bior3.1
bior3.3 bior3.5 bior3.7 bior3.9
bior4.4 bior5.5 bior6.8
===================================
ReverseBior rbio
------------------------------
rbio1.1 rbio1.3 rbio1.5 rbio2.2
rbio2.4 rbio2.6 rbio2.8 rbio3.1
rbio3.3 rbio3.5 rbio3.7 rbio3.9
rbio4.4 rbio5.5 rbio6.8
===================================
Meyer meyr
===================================
DMeyer dmey
===================================
Gaussian gaus
------------------------------
gaus1 gaus2 gaus3 gaus4
gaus5 gaus6 gaus7 gaus8
gaus**
===================================
Mexican_hat mexh
===================================
Morlet morl
===================================
Complex Gaussian cgau
------------------------------
cgau1 cgau2 cgau3 cgau4
cgau5 cgau**
===================================
Shannon shan
------------------------------
shan1-1.5 shan1-1 shan1-0.5 shan1-0.1
shan2-3 shan**
===================================
Frequency B-Spline fbsp
------------------------------
fbsp1-1-1.5 fbsp1-1-1 fbsp1-1-0.5 fbsp2-1-1
fbsp2-1-0.5 fbsp2-1-0.1 fbsp**
===================================
Complex Morlet cmor
------------------------------
cmor1-1.5 cmor1-1 cmor1-0.5 cmor1-1
cmor1-0.5 cmor1-0.1 cmor**
===================================
Lemarie lem
------------------------------
lem1 lem2 lem3 lem4 lem5
===================================
% Restore initial wavelets.
%
% Restore the initial ascii file
% 'wavelets.ini' and initial
% MAT-file 'wavelets.bin'.
wavemngr('restore',0);
% List wavelets families.
wavemngr('read')
ans =
===================================
Haar haar
Daubechies db
Symlets sym
Coiflets coif
BiorSplines bior
ReverseBior rbio
Meyer meyr
DMeyer dmey
Gaussian gaus
Mexican_hat mexh
Morlet morl
Complex Gaussian cgau
Shannon shan
Frequency B-Spline fbsp
Complex Morlet cmor
===================================
% Add new family of orthogonal wavelets.
wavemngr('add','Lemarie','lem',1,'1 2 3','lemwavf');
% All command line capabilities are available for
% the new wavelets.
%
% Example 1: compute the four associated filters.
[Lo_D,Hi_D,Lo_R,Hi_R] = wfilters('lem3');
% Example 2: compute scale and wavelet functions.
[phi,psi,xval] = wavefun('lem3');
% Add a new family of orthogonal wavelets: special form
% for the GUI mode.
%
% The M-file lemwavf allows you to compute the filter for
% any order. If you want to get a popup of the form
% 1 2 3 **, associated with the family, then wavelets are
% appended for GUI mode using:
wavemngr('restore',0);
wavemngr('add','Lemarie','lem',1,'1 2 3 **','lemwavf');
% After this sequence, all GUI capabilities are available for
% the new wavelets.
% Note that the last command allows a short cut in the
% order definition only if possible orders are integers.
Caution wavemngr works on the current directory. If you add a new wavelet family, it is available in this directory only. Refer to Chapter 7 of the User's Guide, "Adding Your Own Wavelets".
Limitations
wavemngr allows you to add a new wavelet. You must verify that it is truly a wavelet. No check is performed either about this point or about the type of the new wavelet.
| | wavemenu | waverec | |