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.
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 | wavenames |
© 1994-2005 The MathWorks, Inc.