Integrate wavelet function psi (ψ)
[INTEG,XVAL] = intwave(
'wname'
,PREC)
[INTDEC,XVAL,INTREC] = intwave('wname'
,PREC)
[INTEG,XVAL] = intwave('wname'
,PREC)
[INTEG,XVAL] = intwave('wname'
,PREC,0)
[INTEG,XVAL] = intwave('wname'
)
[INTEG,XVAL] = intwave('wname'
,8)
intwave('wname
',IN2,IN3),
PREC = max(IN2,IN3)
intwave('wname'
,0)
intwave('wname'
,8,IN3)
intwave('wname'
)
intwave('wname'
,8)
[INTEG,XVAL] = intwave(
computes
the integral, 'wname'
,PREC)INTEG
, of the wavelet function ψ
(from −∞ to XVAL
values): for x in
XVAL
.
The function ψ is approximated on the 2PREC points grid XVAL
where PREC
is
a positive integer. 'wname'
is a character
vector containing the name of the wavelet ψ (see wfilters
for more information).
Output argument INTEG
is a real or
complex vector depending on the wavelet type.
For biorthogonal wavelets,
[INTDEC,XVAL,INTREC] = intwave(
computes
the integrals, 'wname'
,PREC)INTDEC
and INTREC
,
of the wavelet decomposition function ψdec and
the wavelet reconstruction function ψrec.
[INTEG,XVAL] = intwave(
is
equivalent to 'wname'
,PREC)[INTEG,XVAL] = intwave(
.'wname'
,PREC,0)
[INTEG,XVAL] = intwave(
is
equivalent to 'wname'
)[INTEG,XVAL] = intwave(
. 'wname'
,8)
When used with three arguments intwave('
and plots are given.wname
',IN2,IN3),
PREC = max(IN2,IN3)
When IN2
is equal to the special value 0, intwave(
is
equivalent to 'wname'
,0)intwave(
.'wname'
,8,IN3)
intwave(
is
equivalent to 'wname'
)intwave(
.'wname'
,8)
intwave
is used only for continuous analysis
(see cwt
for more information).
% Set wavelet name. wname = 'db4'; % Plot wavelet function. [phi,psi,xval] = wavefun(wname,7); subplot(211); plot(xval,psi); title('Wavelet'); % Compute and plot wavelet integrals approximations % on a dyadic grid. [integ,xval] = intwave(wname,7); subplot(212); plot(xval,integ); title(['Wavelet integrals over [-Inf x] ' ... 'for each value of xval']);
First, the wavelet function is approximated on a grid of 2PREC points
using wavefun
. A piecewise constant interpolation
is used to compute the integrals using cumsum
.