小波包频段划分供参习.doc

MATLAB小波包的频率分段我总觉得不是按照自然顺序排的，通过测得响应频谱所得的结果能明显看出来。可苦于找不到相关的说明，哎，踏破铁鞋无觅处得来全不费工夫，原来帮助里早都说了，只怪自己急功近利没好好学习，骑驴找驴，反倒浪费不少时间。 Building Wavelet Packets The computation scheme for wavelet packets generation is easy when using an orthogonal wavelet. We start with the two filters of length 2N, where h(n) and g(n), corresponding to the wavelet. Now by induction let us define the following sequence of functions: (Wn(x), n = 0, 1, 2, ...) by where W0(x) = φ(x) is the scaling function and W1(x) = ψ(x) is the wavelet function. For example for the Haar wavelet we have The equations become W0(x) = φ(x) is the Haar scaling function and W1(x) = ψ(x) is the Haar wavelet, both supported in [0, 1]. Then we can obtain W2n by adding two 1/2-scaled versions of Wn with distinct supports [0,1/2] and [1/2,1] and obtain W2n+1 by subtracting the same versions of Wn. For n = 0 to 7, we have the W-functions shown in Figure?6-37. Figure 6-37: Haar Wavelet Packets This can be obtained using the following command: [wfun,xgrid] = wpfun(db1,7,5); which returns in wfun the approximate values of Wn for n = 0 to 7, computed on a 1/25 grid of the support xgrid. Starting from more regular original wavelets and using a similar construction, we obtain smoothed versions of this system of W-functions, all with support in the interval [0, 2N-1]. Figure?6-38 presents the system of W-functions for the original db2 wavelet. Figure 6-38: db2 Wavelet Packets Wavelet Packet Atoms Starting from the functions and following the same line leading to orthogonal wavelets, we consider the three-indexed family of analyzing functions (the waveforms): . As in the wavelet framework, k can be interpreted as a time-localization parameter and j as a scale parameter. So what is the interpretation of n? The basic idea of the wavelet packets is that for fixed values of j and k, Wj,n,k analyzes the fluctuations of the signal roughly around the position and at various frequencies for the different admissible values of the