Documentation

This is machine translation

Translated by Microsoft
Mouse over text to see original. Click the button below to return to the English verison of the page.

symaux

Symlet wavelet filter computation

Syntax

W = SYMAUX(N,SUMW)
W = SYMAUX(N)
W = SYMAUX(N,1)
W = SYMAUX(N,0)
W = SYMAUX(N,1)

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

|

Introduced before R2006a

Was this topic helpful?