mwr                package:wavethresh                R Documentation

_M_u_l_t_i_p_l_e _d_i_s_c_r_e_t_e _w_a_v_e_l_e_t _t_r_a_n_s_f_o_r_m (_r_e_c_o_n_s_t_r_u_c_t_i_o_n).

_D_e_s_c_r_i_p_t_i_o_n:

     This function performs the reconstruction stage of Mallat's
     pyramid algorithm adapted for multiple wavelets (see Xia et
     al.(1996)), i.e. the discrete inverse _multiple_ wavelet
     transform.

_U_s_a_g_e:

     mwr(mwd, prefilter.type = mwd$prefilter, verbose = FALSE, start.level = 0,
             returnC = FALSE)

_A_r_g_u_m_e_n_t_s:

     mwd: A multiple wavelet decomposition object as returned by 'mwd'.

prefilter.type: Usually best not to change this (i.e. not to use a
          different prefilter on the reconstruction to the one used on
          decomposition).

 verbose: Controls the printing of "informative" messages whilst the
          computations progress. Such messages are generally annoying
          so it is turned off by default.

start.level: The level you wish to start reconstruction at. The is
          usually the first (level 0).

 returnC: If this is FALSE then a vector of the same length as the
          argument data supplied to the function 'mwd' that constructed
          the supplied 'mwd.object'. is returned, Ie. the reconstructed
          data. If true then the last level (highest resolution) C
          coefficients are returned in matrix form. This matrix has not
          been postprocessed. 

_D_e_t_a_i_l_s:

     The code implements Mallat's pyramid algorithm adapted for
     multiple wavelet decompositions (Xia et al. 1996). In the
     reconstruction the quadrature mirror filters G and H are supplied
     with C0 and D0, D1, ... D(J-1) (the wavelet coefficients) and
     rebuild C1,..., CJ. 

     The matrix CJ is postprocessed which returns the full
     reconstruction 

     If 'mwd.object' was obtained directly from 'mwd' then the original
     function can be reconstructued exactly. Usually, the 'mwd.object'
     has been modified in some way, for example, some coefficients set
     to zero by 'threshold'. Mwr then reconstructs the function with
     that set of wavelet coefficients. 

     See also Downie and Silverman, 1998

_V_a_l_u_e:

     Either a vector containing the final reconstruction or a matrix
     containing unpostprocessed coefficients.

_R_E_L_E_A_S_E:

     Version 3.9.6 (Although Copyright Tim Downie 1996)

_S_e_e _A_l_s_o:

     'accessC.mwd', 'accessD.mwd', 'draw.mwd', 'mfirst.last',
     'mfilter.select', 'mwd', 'mwd.object', 'plot.mwd', 'print.mwd',
     'putC.mwd', 'putD.mwd', 'summary.mwd', 'threshold.mwd', 'wd',
     'wr.mwd'.

_E_x_a_m_p_l_e_s:

     #
     # Decompose and then exactly reconstruct test.data
     #
     #tdecomp <- mwd(test.data)
     #trecons <- mwr(tdecomp)
     #
     # Look at accuracy of reconstruction
     # max(abs(trecons - test.data))
     #[1] 2.266631e-12
     #
     # See also the example of using this-is-escaped-codenormal-bracket68bracket-normal or mwr in
     # the this-is-escaped-codenormal-bracket70bracket-normal section of
     # the help for this-is-escaped-codenormal-bracket72bracket-normal.

