How to extract the dyadic wavelet coefficeints from the DWT (decimated) 'wavedec' function output

Question

Terry 2012-10-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/49653-how-to-extract-the-dyadic-wavelet-coefficeints-from-the-dwt-decimated-wavedec-function-output

The coefficients output from the function 'wavedec' are the artifacts from the convolution of the input signal (vector) and the wavelet/scale functions. Thus the length of the coefficients from any scale is a function of the lengths of the wavelet/scale functions and the dyadic of the signal at the level being decomposed. As an example if the input signal is 64 samples long and a Daubechies wavelet with four vanishing moments is used (db4) then the lengths of the coefficients produced are 7,7,8,10,14,21, and 35 for ca1, cd1, cd2, ….. cd6. The actual lengths should be 1,1,2,4,8,16, and 32. This can be verified by performing the decomposition manually without convolution, ie using the’ sumproduct’ function with the proper wavelet/scale function in Excel for instance. How then do you select the dyadic coefficients from the output from ‘wavedec”? Sometimes is looks like taking the middle n coefficients for high numbered levels works and sometimes it looks like taking the last n coefficients works for low numbered levels. What’s up?s

added note: The dwtmode was 'sym' for the above discussion. When dwtmode = 'per' the dyadic coefficients result. Modify the question to mean for modes other than 'per'

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Wayne King 2012-10-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/49653-how-to-extract-the-dyadic-wavelet-coefficeints-from-the-dwt-decimated-wavedec-function-output#answer_60677

编辑：Wayne King 2012-10-2

在 MATLAB Online 中打开

You can set the dwtmode to 'per' and then use detcoef() or appcoef()

   dwtmode('per')
   x = randn(1024,1);
   [C,L] = wavedec(x,4,'db4');
   d1 = detcoef(C,L,1);

To get all the detail coefficients in a cell array

det = detcoef(C,L,'cells');

Or you can do something like:

[d1,d2,d3,d4] = detcoef(C,L,[1 2 3 4]);

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Terry 2012-10-4

Thanks for the reply but the question was really about the coefficients from decompositions from extended signals that do not have dyadic lengths. Sorry for the confusion. I added the note at the bottom after publishing to clarify.

After some thought I am wondering what the consequence is for applying a windowing function like a Blackman or other that will force stationarity on the signal and remove the edge effects from the wavelet transform. I suppose there will be some price to pay and it may be worse than doing nothing.

请先登录，再进行评论。