% function [cD cA] = getDWT(X,N,Name)
%
% finds the discrete wavelet transform at level N for signal X using the
% wavelet specified by Name.
%
%% Inputs:
% X: the original signal
% N: the decomposition level
% Name: the wavelet name to use
%
%% Outputs:
% cD is a N-row matrix containing the detail coefficients up to N levels
% cA is the same for the approximations
% This code was developed by David Springer for comparison purposes in the
% paper:
% D. Springer et al., "Logistic Regression-HSMM-based Heart Sound
% Segmentation," IEEE Trans. Biomed. Eng., In Press, 2015.
%
%% Copyright (C) 2016 David Springer
% dave.springer@gmail.com
%
% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% any later version.
%
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program. If not, see .
function [cD cA] = getDWT(X,N,Name)
%No DWT available for Morlet - therefore perform CWT:
if(strcmp(Name,'morl'))
c = cwt(X,1:N,'morl');
cD = c;
cA = c;
else
%Preform wavelet decomposition
[c,l] = wavedec(X,N,Name);
%Reorder the details based on the structure of the wavelet
%decomposition (see help in wavedec.m)
len = length(X);
cD = zeros(N,len);
for k = 1:N
d = detcoef(c,l,k);
d = d(:)';
d = d(ones(1,2^k),:);
cD(k,:) = wkeep1(d(:)',len);
end
cD = cD(:);
%Space cD according to spacing of floating point numbers:
I = find(abs(cD)