Re: wavelet

*To*: mathgroup at smc.vnet.net*Subject*: [mg114730] Re: wavelet*From*: Bill Rowe <readnews at sbcglobal.net>*Date*: Tue, 14 Dec 2010 06:56:49 -0500 (EST)

On 12/13/10 at 3:50 AM, dauphinester at gmail.com (clansa) wrote: >The order StandardDeviation[Flatten[dwd[{___,1},"Values"]] compute >the standard deviation of all the wavelets coefficients dwd >But I want to calculate standard deviations of the wavelets >coefficients for each level of the decomposition Do you know the >solution? Here are a couple of ways In[1]:= data = RandomReal[1, {16}]; In[2]:= dwd = DiscreteWaveletTransform[data]; In[3]:= Normal[dwd] /. HoldPattern[a_ -> b_] :> If[Length[b] > 1, {a, StandardDeviation[b]}, {a, {}}] Out[3]= {{{0}, 0.342116}, {{1}, 0.342662}, {{0, 0}, 0.429588}, {{0, 1}, 0.297391}, {{0, 0, 0}, 0.428882}, {{0, 0, 1}, 0.297418}, {{0, 0, 0, 0}, {}}, {{0, 0, 0, 1}, {}}} In[4]:= StandardDeviation /@ Most[dwd[{___, 1}, "Values"]] Out[4]= {0.342662,0.297391,0.297418}