← Back to branch summary

~ubuntu-branches/ubuntu/vivid/octave-ltfat/vivid-proposed

« back to all changes in this revision

Viewing changes to inst/wavelets/wfilt_symds.m

  • Committer: Package Import Robot
  • Author(s): Rafael Laboissiere
  • Date: 2014-10-09 13:58:00 UTC
  • mfrom: (1.1.1)
  • Revision ID: package-import@ubuntu.com-20141009135800-5kw1gvovmqh34bmm
Tags: 2.0.1-1
* Imported Upstream version 2.0.1
* d/p/autoload-yes.patch: Refresh patch
* d/p/add-subdirs-to-loadpath.patch: Drop patch (fixed upstream)
* d/check.m: Use the upstream unit testing script
* d/rules: Fix permission of all *.m files

Show diffs side-by-side

added added

removed removed

Lines of Context:
inst/wavelets/wfilt_symds.m
1
 
function [h,g,a,info] = wfilt_symds(N)
 
1
function [h,g,a,info] = wfilt_symds(K)
2
2
%-*- texinfo -*-
3
3
%@deftypefn {Function} wfilt_symds
4
4
%@verbatim
5
5
%WFILT_SYMDS  Symmetric wavelets dyadic sibling
6
 
%   Usage: [h,g,a] = wfilt_symds(N);
7
 
%
8
 
%   [h,g,a]=WFILT_SYMDS(N) Returns symmetric dyadic sibling wavelet filters.
9
 
%   The redundancy is equal to 2.
 
6
%   Usage: [h,g,a] = wfilt_symds(K);
 
7
%
 
8
%   [h,g,a]=WFILT_SYMDS(K) with K in {1,2,3,4,5} returns symmetric 
 
9
%   dyadic sibling wavelet frame filters from the reference.
 
10
%
 
11
%   The returned filterbank has redundancy equal to 2 and it does not form
 
12
%   a tight frame.
10
13
%
11
14
%   Examples:
12
15
%   ---------
13
16
%   :
14
 
%     figure(1);
15
17
%     wfiltinfo('ana:symds3');
16
 
%     figure(2);
 
18
%
 
19
%   :
17
20
%     wfiltinfo('syn:symds3');
18
21
% 
19
22
%   References:
31
34
%@end deftypefn
32
35
 
33
36
% Copyright (C) 2005-2014 Peter L. Soendergaard <soender@users.sourceforge.net>.
34
 
% This file is part of LTFAT version 1.4.4
 
37
% This file is part of LTFAT version 2.0.1
35
38
%
36
39
% This program is free software: you can redistribute it and/or modify
37
40
% it under the terms of the GNU General Public License as published by
46
49
% You should have received a copy of the GNU General Public License
47
50
% along with this program.  If not, see <http://www.gnu.org/licenses/>.
48
51
 
 
52
% AUTHOR: Zdenek Prusa
 
53
 
49
54
info.istight = 0;
50
55
a = [2;2;2;2];
51
56
 
52
 
switch(N)
 
57
switch(K)
53
58
 case 1
54
59
    % Example 1. Not a tight frame!
55
60
    harr = [
230
235
end
231
236
 
232
237
g=mat2cell(garr,size(garr,1),ones(1,size(garr,2)));
 
238
g = cellfun(@(gEl) struct('h',gEl(:),'offset',-floor(numel(gEl)/2)+1),g,'UniformOutput',0);
 
239
 
233
240
h=mat2cell(flipud(harr),size(harr,1),ones(1,size(harr,2)));
 
241
h = cellfun(@(hEl) struct('h',hEl(:),'offset',-floor(numel(hEl)/2)+1),h,'UniformOutput',0);
234
242
 
235
243