coef_dcov2d.m

多小波matlab程序包

function [HH,HL,LH]=coef_dcov2D(pflt)
%[HH,HL,LH]=coef_dcov2D(pflt)
%
%  This function returns diagonal blocks of the covariance matrix 
%  corresponding to 5 levels of 2-dimensional multiwavelet transform 
%  with preprocessing. For details see [SW].
%
%
%  [SW]  V. Strela and A. T. Walden, "Signal and Image Denoising via Wavelet 
%        Thresholding: Orthogonal and Biorthogonal, Scalar and Multiple Wavelet 
%        Transforms", Imperial College, Statistics Section, 
%        Technical Report TR-98-01 (1998).
%
%
%  Input:                                               
%    pflt       string of characters, name of transform used to generate 
%               the covariance matrix; for admissible names see below 
%
%  Output:                                                  
%    HH         r^2 by 5*r^2 real array, diagonal blocks of the high-high part
%               (high-pass filtering in both directions) of the covariance 
%               matrix; r is the number of scaling functions;
%               blocks in HH are organized as follows:HH=[HH1 HH2 ... HH5], 
%               HHk corresponds to the k-th level of the transform 
%    HL         r^2 by 5*r^2 real array, diagonal blocks of the high-low part 
%               (low-pass horizontal and high-pass vertical filtering) 
%               of the covariance matrix; r is the number of scaling functions;
%               blocks in HL are organized as follows: HL=[HL1 HL2 ... HL5], 
%               HLk corresponds to the k-th level of the transform 
%    LH         r^2 by 5*r^2 real array, diagonal blocks of the high-low part 
%               (high-pass horizontal and low-pass vertical filtering) 
%               of the covariance matrix; r is the number of scaling functions;
%               blocks in LH are organized as follows: LH=[LH1 LH2 ... LH5], 
%               LHk corresponds to the k-th level of the transform 
%
%  Admissible Names of the Transforms:
%  (names of wavelet transforms are from coef.m, 
%   names of prefilters are from coef_prep.m) 
%   'ghmap'     'ghm' multiwavelet transform with 'ghmap' prefilter  
%   'ghmorap'   'ghm' multiwavelet transform with 'ghmorap' prefilter
%   'ghmrr'     'ghm' multiwavelet transform with oversampled preprocessing
%   'clap'      'cl' multiwavelet transform with 'clap' prefilter 
%   'clrr'      'cl' multiwavelet transform with oversampled  preprocessing 
%   'bih5ap'    'bih52s' multiwavelet transform with 'bih5ap' prefilter        
%   'bih52srr'  'bih52s' multiwavelet transform with oversampled preprocessing
%   'sa4ap'     'sa4' multiwavelet transform with 'sa4ap' prefilter 
%   'bi9'       'bi9' wavelet transform without preprocessing
%
%  Example of Usage:
%   [HH,HL,LH]=coef_dcov2D('bih52sap')                

% Author: Vasily Strela
% COPYRIGHT 1997,98 by Vasily Strela

if strcmp(pflt,'ghmap')
  HH=[   3.1651876e-01   0.0000000e+00   0.0000000e+00   0.0000000e+00   1.1109160e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   1.8273632e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   2.2061161e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   2.3292864e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00
   0.0000000e+00   3.8611238e-01   0.0000000e+00   0.0000000e+00   0.0000000e+00   1.3012684e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   1.9184746e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   2.2395353e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   2.3382910e+00   0.0000000e+00   0.0000000e+00
   0.0000000e+00   0.0000000e+00   3.8611238e-01   0.0000000e+00   0.0000000e+00   0.0000000e+00   1.3012684e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   1.9184746e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   2.2395353e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   2.3382910e+00   0.0000000e+00
   0.0000000e+00   0.0000000e+00   0.0000000e+00   4.7100769e-01   0.0000000e+00   0.0000000e+00   0.0000000e+00   1.5242372e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   2.0141286e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   2.2734608e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   2.3473304e+00];
  HL=[   7.6355369e-01   7.2446437e-02   0.0000000e+00   0.0000000e+00   1.5388861e+00   5.5783251e-02   0.0000000e+00   0.0000000e+00   2.0550199e+00   2.0980846e-02   0.0000000e+00   0.0000000e+00   2.2804442e+00   6.4230886e-03   0.0000000e+00   0.0000000e+00   2.3505417e+00   1.6919509e-03   0.0000000e+00   0.0000000e+00
   7.2446437e-02   6.9628343e-01   0.0000000e+00   0.0000000e+00   5.5783251e-02   1.5188212e+00   0.0000000e+00   0.0000000e+00   2.0980846e-02   2.0423369e+00   0.0000000e+00   0.0000000e+00   6.4230886e-03   2.2779038e+00   0.0000000e+00   0.0000000e+00   1.6919509e-03   2.3496387e+00   0.0000000e+00   0.0000000e+00
   0.0000000e+00   0.0000000e+00   9.3143778e-01   8.8375382e-02   0.0000000e+00   0.0000000e+00   1.8025700e+00   6.5341557e-02   0.0000000e+00   0.0000000e+00   2.1574820e+00   2.2026940e-02   0.0000000e+00   0.0000000e+00   2.3149894e+00   6.5203884e-03   0.0000000e+00   0.0000000e+00   2.3596285e+00   1.6984917e-03
   0.0000000e+00   0.0000000e+00   8.8375382e-02   8.4937668e-01   0.0000000e+00   0.0000000e+00   6.5341557e-02   1.7790670e+00   0.0000000e+00   0.0000000e+00   2.2026940e-02   2.1441666e+00   0.0000000e+00   0.0000000e+00   6.5203884e-03   2.3124105e+00   0.0000000e+00   0.0000000e+00   1.6984917e-03   2.3587220e+00];
  LH=[   7.6355369e-01   0.0000000e+00   7.2446437e-02   0.0000000e+00   1.5388861e+00   0.0000000e+00   5.5783251e-02   0.0000000e+00   2.0550199e+00   0.0000000e+00   2.0980846e-02   0.0000000e+00   2.2804442e+00   0.0000000e+00   6.4230886e-03   0.0000000e+00   2.3505417e+00   0.0000000e+00   1.6919509e-03   0.0000000e+00
   0.0000000e+00   9.3143778e-01   0.0000000e+00   8.8375382e-02   0.0000000e+00   1.8025700e+00   0.0000000e+00   6.5341557e-02   0.0000000e+00   2.1574820e+00   0.0000000e+00   2.2026940e-02   0.0000000e+00   2.3149894e+00   0.0000000e+00   6.5203884e-03   0.0000000e+00   2.3596285e+00   0.0000000e+00   1.6984917e-03
   7.2446437e-02   0.0000000e+00   6.9628343e-01   0.0000000e+00   5.5783251e-02   0.0000000e+00   1.5188212e+00   0.0000000e+00   2.0980846e-02   0.0000000e+00   2.0423369e+00   0.0000000e+00   6.4230886e-03   0.0000000e+00   2.2779038e+00   0.0000000e+00   1.6919509e-03   0.0000000e+00   2.3496387e+00   0.0000000e+00
   0.0000000e+00   8.8375382e-02   0.0000000e+00   8.4937668e-01   0.0000000e+00   6.5341557e-02   0.0000000e+00   1.7790670e+00   0.0000000e+00   2.2026940e-02   0.0000000e+00   2.1441666e+00   0.0000000e+00   6.5203884e-03   0.0000000e+00   2.3124105e+00   0.0000000e+00   1.6984917e-03   0.0000000e+00   2.3587220e+00];

elseif strcmp(pflt,'sa4ap')
 A=diag(ones(1,4));
 HH=[A A A A A];
 HL=HH;
 LH=HH;

elseif strcmp(pflt,'ghmrr')
  HH=[   1.8225000e-02   1.4316750e-02   1.4316750e-02   1.1246603e-02   4.0214622e-01  -6.7219900e-03  -6.7219900e-03   1.1236000e-04   1.5940325e+00   1.3256775e-03   1.3256775e-03   1.1025000e-06   2.0730240e+00  -1.4398000e-04  -1.4398000e-04   1.0000000e-08   2.2049280e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00
   1.4316750e-02   9.0450000e-02   1.1246603e-02   7.1053500e-02  -6.7219900e-03   6.1785235e-01   1.1236000e-04  -1.0327580e-02   1.3256775e-03   1.7118284e+00   1.1025000e-06   1.4236425e-03  -1.4398000e-04   2.1070753e+00   1.0000000e-08  -1.4634500e-04   0.0000000e+00   2.2137632e+00   0.0000000e+00   0.0000000e+00
   1.4316750e-02   1.1246603e-02   9.0450000e-02   7.1053500e-02  -6.7219900e-03   1.1236000e-04   6.1785235e-01  -1.0327580e-02   1.3256775e-03   1.1025000e-06   1.7118284e+00   1.4236425e-03  -1.4398000e-04   1.0000000e-08   2.1070753e+00  -1.4634500e-04   0.0000000e+00   0.0000000e+00   2.2137632e+00   0.0000000e+00
   1.1246603e-02   7.1053500e-02   7.1053500e-02   4.4890000e-01   1.1236000e-04  -1.0327580e-02  -1.0327580e-02   9.4926049e-01   1.1025000e-06   1.4236425e-03   1.4236425e-03   1.8383292e+00   1.0000000e-08  -1.4634500e-04  -1.4634500e-04   2.1416859e+00   0.0000000e+00   0.0000000e+00   0.0000000e+00   2.2226337e+00];