dwt_filter_tcase.m

时间序列分析中很用的源码,书的原名为时间序列分析的小波方法.

function tc = dwt_filter_tcase
% dwt_filter_tcase -- munit test case to test dwt_filter.
%
%****f*  Tests.dwt/dwt_filter_tcase
%
% NAME
%   dwt_filter_tcase -- munit test case to test dwt_filter.
%
% USAGE
%   run_tcase('dwt_filter_tcase')
%
% INPUTS
%
%
% OUTPUTS
%   tc            = tcase structure for dwt_filter testcase.
%
% SIDE EFFECTS
%
%
% DESCRIPTION
%
%
% SEE ALSO
%   dwt_filter
%
  
% AUTHOR
%   Charlie Cornish
%
% CREATION DATE
%   2005-03-01
%
% COPYRIGHT
%   (c) 2004, 2005 Charles R. Cornish
%
% CREDITS
%
%
% REVISION
%   $Revision: 615 $
%
%***

%   $Id: dwt_filter_tcase.m 615 2006-04-27 21:45:17Z ccornish $

tc = MU_tcase_new(mfilename);
tc = MU_tcase_add_test(tc, @test_smoke_test);
tc = MU_tcase_add_test(tc, @test_no_arg);
tc = MU_tcase_add_test(tc, @test_invalid_wavelet);
tc = MU_tcase_add_test(tc, @test_arg_list);
tc = MU_tcase_add_test(tc, @test_invalid_wavelet_empty);
tc = MU_tcase_add_test(tc, @test_verify_haar_filter);
tc = MU_tcase_add_test(tc, @test_verify_d4_filter);
tc = MU_tcase_add_test(tc, @test_verify_la8_filter);

return

function test_smoke_test(mode)
  % Test Description:  
  %    Expected result: successful executiong
  [wtf] = dwt_filter('la8');
return

function test_no_arg(mode)
  % Test Description:  
  %    Expected error: WMTSA:InvalidNumArguments
  try
     dwt_filter;
  catch
    [errmsg, msg_id] = lasterr;
    MU_assert_error('MATLAB:nargchk:notEnoughInputs', msg_id);
  end
return

function test_invalid_wavelet(varargin)
  % Test Description:  
  %    Expected error:  Invalid wavelet
  try
    [wtf] = dwt_filter('not_a_wavelet');
  catch
    [errmsg, msg_id] = lasterr;
    MU_assert_error('WMTSA:invalidWaveletTransformFilter', msg_id);
  end
return

function test_arg_list(mode)
  % Test Description:  
  %    Expected result:  filter list
  dwt_filter('list');
return


function test_arg_verify_filter(mode)
  % Test Description:  
  %    Expected result: 
  exp_wtf.h = [0.7071 0.7071];
  exp_wtf.g = [0.7071 -0.7071];

  exp_wtf.L  = 2;
  exp_wtf.Name = 'Haar';
    
  [wtf] = dwt_filter('haar');
  
  MU_assert_isequal(exp_wtf.h, wtf.h);
  MU_assert_isequal(exp_wtf.g, wtf.g);
  MU_assert_isequal(exp_wtf.L, wtf.L);
  MU_assert_isequal(exp_wtf.Name, wtf.Name);
  
return

function test_invalid_wavelet_empty(varargin)
  % Test Description:  
  %    Expected error:  Invalid wavelet (empty string)
  try
    [wtf] = dwt_filter('');
  catch
    [errmsg, msg_id] = lasterr;
    MU_assert_error('WMTSA:invalidWaveletTransformFilter', msg_id);
  end
return


function test_verify_haar_filter(varargin)
  % Test Description:  
  %    Verify the results for the Haar filter.
  g_dwt = [0.7071067811865475 0.7071067811865475];
  exp_wtf.g = g_dwt;
  exp_wtf.h = flipdim(exp_wtf.g, 2) .* (-1).^([1:length(exp_wtf.g)]-1);
  exp_wtf.L = 2;
  exp_wtf.Name = 'Haar';
  [wtf] = dwt_filter('haar');
  MU_assert_isequal(exp_wtf.h, wtf.h);
  MU_assert_isequal(exp_wtf.g, wtf.g);
  MU_assert_isequal(exp_wtf.L, wtf.L);
  MU_assert_isequal(exp_wtf.Name, wtf.Name);
return

function test_verify_d4_filter(varargin)
  % Test Description:  
  %    Verify the results for the D4 filter.
  g_dwt = [0.4829629131445341  0.8365163037378077 ...
           0.2241438680420134 -0.1294095225512603];
  exp_wtf.g = g_dwt;
  exp_wtf.h = flipdim(exp_wtf.g, 2) .* (-1).^([1:length(exp_wtf.g)]-1);
  exp_wtf.L = 4;
  exp_wtf.Name = 'D4';
  [wtf] = dwt_filter('d4');
  MU_assert_isequal(exp_wtf.h, wtf.h);
  MU_assert_isequal(exp_wtf.g, wtf.g);
  MU_assert_isequal(exp_wtf.L, wtf.L);
  MU_assert_isequal(exp_wtf.Name, wtf.Name);
return

function test_verify_la8_filter(varargin)
  % Test Description:  
  %    Verify the results for the LA8 filter.
  g_dwt = [-0.0757657147893407 -0.0296355276459541, ...
           0.4976186676324578  0.8037387518052163, ...
           0.2978577956055422 -0.0992195435769354, ...
           -0.0126039672622612  0.0322231006040713];
  exp_wtf.g = g_dwt;
  exp_wtf.h = flipdim(exp_wtf.g, 2) .* (-1).^([1:length(exp_wtf.g)]-1);
  exp_wtf.L = 8;
  exp_wtf.Name = 'LA8';
  [wtf] = dwt_filter('la8');
  MU_assert_isequal(exp_wtf.h, wtf.h);
  MU_assert_isequal(exp_wtf.g, wtf.g);
  MU_assert_isequal(exp_wtf.L, wtf.L);
  MU_assert_isequal(exp_wtf.Name, wtf.Name);
  
  MU_assert_isequal(0,1);
return