双树复小波

这是一个双树复复小波变换的文件，里面有诸多相关示例，已在MATLAB2016a上仿真验证

function [h_0h_1] = daubcqf（NTYPE）
%    [h_0h_1] = daubcqf（NTYPE）; 
%
%    Function computes the Daubechies‘ scaling and wavelet filters
%    （normalized to sqrt（2））.
%
%    Input: 
%       N    : Length of filter （must be even）
%       TYPE : Optional parameter that distinguishes the minimum phase
%              maximum phase and mid-phase solutions （‘min‘ ‘max‘ or
%              ‘mid‘）. If no argument is specified the minimum phase
%              solution is used.
%
%    Output: 
%       h_0 : Minimal phase Daubechies‘ scaling filter 
%       h_1 : Minimal phase Daubechies‘ wavelet filter 
%
%    Example:
%       N = 4;
%       TYPE = ‘min‘;
%       [h_0h_1] = daubcqf（NTYPE）
%       h_0 = 0.4830 0.8365 0.2241 -0.1294
%       h_1 = 0.1294 0.2241 -0.8365 0.4830
%
%    Reference: “Orthonormal bases of Compactly Supported Wavelets“
%                CPAM Oct.89 
%

%File Name: daubcqf.m
%Last Modification Date: 1/2/96	15:12:57
%Current Version: daubcqf.m	1.15
%File Creation Date: 10/10/88
%Author: Ramesh Gopinath  
%
%Copyright: All software documentation and related files in this distribution
%           are Copyright （c） 1988  Rice University
%
%Permission is granted for use and non-profit distribution providing that this
%notice be clearly maintained. The right to distribute any portion for profit
%or as part of any commercial product is specifically reserved for the author.
%

if（nargin < 2）
  TYPE = ‘min‘;
end;
if（rem（N2） ~= 0）
  error（‘No Daubechies filter exists for ODD length‘）;
end;
K = N/2;
a = 1;
p = 1;
q = 1;
h_0 = [1 1];
for j  = 1:K-1
  a = -a * 0.25 * （j + K - 1）/j;
  h_0 = [0 h_0] + [h_0 0];
  p = [0 -p] + [p 0];
  p = [0 -p] + [p 0];
  q = [0 q 0] + a*p;
end;
q = sort（roots（q））;
qt = q（1:K-1）;
if TYPE==‘mid‘
  if rem（K2）==1  
    qt = q（[1:4:N-2 2:4:N-2]）;
  else
    qt = q（[1 4:4:K-1 5:4:K-1 N-3:-4:K N-4:-4:K]）;
  end;
end;
h_0 = conv（h_0real（poly（qt）））;
h_0 = sqrt（2）*h_0/sum（h_0）; 	%Normalize to sqrt（2）;
if（TYPE==‘max‘）
  h_0 = fliplr（h_0）;
end;
if（abs（sum（h_0 .^ 2））-1 > 1e-4） 
  error（‘Numerically unstable for this value of “N“.‘）;
end;
h_1 = rot90（h_02）;
h_1（1:2:N）=-h_1（1:2:N）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2178  1999-03-09 11:15  whmt1_0\daubcqf.m

     文件       5306  1999-03-09 11:15  whmt1_0\emhht.m

     文件       5296  1999-03-09 11:15  whmt1_0\emhlt.m

     文件       5336  1999-03-09 11:15  whmt1_0\emlht.m

     文件        437  1999-03-09 21:21  whmt1_0\example1.m

     文件        467  1999-03-09 11:15  whmt1_0\example2.m

     文件        492  1999-03-09 11:15  whmt1_0\example3.m

     文件         81  1999-03-09 11:15  whmt1_0\gauss.m

     文件       3619  2018-07-30 08:59  whmt1_0\hdenoise.m

     文件       1071  2001-08-18 00:20  whmt1_0\hmt1d.tar\denoise_tree.m

     文件        552  2001-08-18 00:20  whmt1_0\hmt1d.tar\dwtlevel1.m

     文件        501  2001-08-18 00:21  whmt1_0\hmt1d.tar\dwtmult1.m

     文件        440  2001-08-18 00:19  whmt1_0\hmt1d.tar\hmtdenoise_tree.m

     文件        501  2001-08-18 00:20  whmt1_0\hmt1d.tar\idwtlevel1.m

     文件        373  2001-08-18 00:20  whmt1_0\hmt1d.tar\idwtmult1.m

     文件       2184  2001-08-18 00:19  whmt1_0\hmt1d.tar\indmixmod.m

     文件       1667  2001-08-18 00:20  whmt1_0\hmt1d.tar\make_tree_data.m

     文件       2263  2001-08-18 00:20  whmt1_0\hmt1d.tar\README

     文件       3373  2001-08-18 00:19  whmt1_0\hmt1d.tar\trainEM_tree.m

     文件       3830  2001-08-18 00:19  whmt1_0\hmt1d.tar\updown_tree.m

     文件       3169  2001-08-18 00:19  whmt1_0\hmt1d.tar\viterbi_tree.m

     文件       1062  1999-03-09 11:15  whmt1_0\hmtdeno.m

     文件       1102  1999-03-09 11:15  whmt1_0\hmtdeno0.m

     文件       1632  1999-03-09 11:15  whmt1_0\hmtmodel.m

     文件    4325760  1999-03-09 11:15  whmt1_0\hmtmodel.mat

     文件       2258  1999-03-09 11:15  whmt1_0\hmttrain.m

     文件        802  1999-03-09 11:15  whmt1_0\Idwt.m

     文件       2737  1999-03-09 11:15  whmt1_0\Idwt2.m

     文件        936  1999-03-09 11:15  whmt1_0\Iidwt.m

     文件       2275  1999-03-09 11:15  whmt1_0\Iidwt2.m

............此处省略16个文件信息