Accelerating the pace of engineering and science

# bswfun

Biorthogonal scaling and wavelet functions

## Syntax

[PHIS,PSIS,PHIA,PSIA,XVAL] = bswfun(LoD,HiD,LoR,HiR)
bswfun(LoD,HiD,LoR,HiR,ITER)
bswfun(LoD,HiD,LoR,HiR,'plot')
bswfun(LoD,HiD,LoR,HiR,ITER,'plot')
bswfun(LoD,HiD,LoR,HiR,'plot',ITER)

## Description

[PHIS,PSIS,PHIA,PSIA,XVAL] = bswfun(LoD,HiD,LoR,HiR) returns approximations on the grid XVAL of the two pairs of scaling function and wavelet (PHIA,PSIA), (PHIS,PSIS) associated with the two pairs of filters (LoD,HiD), (LoR,HiR).

bswfun(LoD,HiD,LoR,HiR,ITER) computes the two pairs of scaling and wavelet functions using ITER iterations.

bswfun(LoD,HiD,LoR,HiR,'plot') or bswfun(LoD,HiD,LoR,HiR,ITER,'plot') or bswfun(LoD,HiD,LoR,HiR,'plot',ITER) computes and plots the functions.

## Examples

```% Start from the Cohen-Daubechies-Feauveau wavelet
% and get the corresponding lifting scheme.
lscdf = liftwave('cdf3.1');

% Visualize the obtained lifting scheme.
displs(lscdf);

lscdf = {...
'p'             [ -0.33333333]              [-1]
'd'             [ -0.37500000 -1.12500000]  [1]
'p'             [  0.44444444]              [0]
[  2.12132034]  [  0.47140452]              []
};

% Transform the lifting scheme to biorthogonal
[LoD,HiD,LoR,HiR] = ls2filt(lscdf);

% Visualize the two pairs of scaling and wavelet
% functions.
bswfun(LoD,HiD,LoR,HiR,'plot');

```