symaux (Wavelet Toolbox)

Wavelet Toolbox

symaux

Symlet wavelet filter computation

Syntax

```
W = SYMAUX(N,SUMW)
```

W = SYMAUX(N)

Description

Symlets are the "least asymmetric" Daubechies' wavelets.

W = SYMAUX(N,SUMW) is the order N Symlet scaling filter such that SUM(W) = SUMW. Possible values for N are 1, 2, 3, ...

Note Instability may occur when N is too large.

W = SYMAUX(N) is equivalent to W = SYMAUX(N,1)

W = SYMAUX(N,0) is equivalent to W = SYMAUX(N,1)

Examples

% Generate wdb4 the order 4 Daubechies scaling filter.
wdb4 = dbaux(4)

wdb4 =

  Columns 1 through 7 

    0.1629    0.5055    0.4461   -0.0198   -0.1323    0.0218    0.0233

  Column 8 

   -0.0075


% wdb4 is a solution of the equation: P = conv(wrev(w),w)*2,
% where P is the "Lagrange à trous" filter for N=4.
% wdb4 is a minimum phase solution of the previous equation,
% based on the roots of P (see dbaux).
P = conv(wrev(wdb4),wdb4)*2;

% Generate wsym4 the order 4 symlet scaling filter.
% The Symlets are the "least asymmetric" Daubechies' 
% wavelets obtained from another choice between the roots of P.
wsym4 = symaux(4)

wsym4 =

  Columns 1 through 7 

    0.0228   -0.0089   -0.0702    0.2106    0.5683    0.3519   -0.0210

  Column 8 

   -0.0536


% Compute conv(wrev(wsym4),wsym4) * 2 and check that wsym4
% is another solution of the equation P = conv(wrev(w),w)*2.
Psym = conv(wrev(wsym4),wsym4)*2;
err = norm(P-Psym)

err =

  7.4988e-016

See Also
symwavf, wfilters

swt2 symwavf