【正文】 title(39。小波包 [2,1]的系數(shù) 39。)。 subplot(324)。 plot(cfs22)。 title(39。小波包 [2,2]的系數(shù) 39。)。 subplot(325)。 plot(cfs31)。 title(39。小波包 [3,1]的系數(shù) 39。)。 subplot(326)。 plot(cfs32)。 title(39。小波包 [3,2]的系數(shù) 39。)。 0 200 400 600 800 1000 1200 1 050510原始信號0 50 100 150 200 250 300 2 0 1 001020小波包 [ 2 , 0 ] 的系數(shù)0 50 100 150 200 250 30050510小波包 [ 2 , 1 ] 的系數(shù)0 20 40 60 80 100 120 140 4 0 2 0020小波包 [ 3 , 0 ] 的系數(shù)0 20 40 60 80 100 120 140 1 050510小波包 [ 3 , 1 ] 的系數(shù)十六、 wprcoef 功能：小波包分解系數(shù)的重構(gòu)； 格式： x＝ wprcoef（ t， n） 說明： wprcoef是一個一維或二維的小波包分析函數(shù)，計算節(jié)點 n的小波包分解系數(shù)的重構(gòu)信號。 X＝ wprcoef（ t）＝ wprcoef（ t， 0） 該函數(shù)一次只能對一個節(jié)點進行重構(gòu)，不能同時對多點進行重構(gòu)，可以通過多次調(diào)用實現(xiàn)。 load noisdopp。 x=noisdopp(1:1000)。 figure(1) subplot(311) plot(x) title(39。原始信號 39。) t=wpdec(x,3,39。db139。,39。shannon39。)。 rcfs=wprcoef(t,[2,0])。 cfs21=wpcoef(t,[2,0])。 subplot(312) plot(cfs21) title(‘小波包節(jié)點（ 2， 0）系數(shù) 39。) subplot(313) plot(rcfs) title(‘重構(gòu)小波包節(jié)點（ 2， 0）信號 39。) 十七、 wpfun 功能：小波包函數(shù) 格式： [WPWS,X] = WPFUN(39。wname39。,NUM,PREC) putes the wavelets packets for a wavelet 39。wname39。 (see WFILTERS), on dyadic intervals of length 1/2^PREC. PREC must be a positive integer. Output matrix WPWS contains the W functions of index from 0 to NUM, stored rowwise as [W0。 W1。...。 Wnum]. Output vector X is the corresponding mon Xgrid vector. [WPWS,X] = WPFUN(39。wname39。,NUM) is equivalent to [WPWS,X] = WPFUN(39。wname39。,NUM,7). 十八、 wpsplt 功能：分解（分割）小波包 格式： t= wpsplt(t,n) returns the modified tree t corresponding to the deposition of the node n. [t,ca,cd] = wpsplt(t,n) with ca = approximation and cd = detail of node n for a 1D deposition. [t, ca,ch,cv,cd] = WPSPLT(T,N) with ca = approximation and ch, cv, cd = (Horiz., Vert. and Diag.) details of node n for a 2D deposition（二維小波變換） 舉例 load noisdopp。 x=noisdopp(1:1000)。 figure(1) subplot(311) plot(x) title(39。原始信號 39。) t=wpdec(x,3,39。db139。,39。shannon39。)。 plot(t) [wpt,wpd]=wpsplt(t,[3,0])。 plot(wpt) 十九、 wpjoin 功能：重新組合小波包 格式及說明： t = wpjoin(t,n) returns the modified tree t corresponding to a reposition of the node n. t = wpjoin(t) is equivalent to t = wpjoin(t,0). [t,x] = wpjoin(t,n) also returns the coefficients of the node n。 [t,x] = wpjoin (t) is equivalent to [t,x] = wpjoin(t,0) load noisdopp。 x=noisdopp(1:1000)。 figure(1) subplot(321) plot(x) title(39。原始信號 39。) t=wpdec(x,3,39。db139。,39。shannon39。)。 plot(t) [wpt,wpc]=wpjoin(t,[1,1])。 plot(wpt) figure(1) subplot(322) plot(wpc) title(39。節(jié)點 2的小波包分解系數(shù) 39。) 二十、 wpcutree 功能：剪切小波包分解樹 格式及說明： t = wpcutree(t,L) cuts the tree t at level L. In addition, [t,rn] = wpcutree(t,L) returns the vector rn which contains the indices of the reconstructed nodes. 舉例 load noisdopp。 x=noisdopp(1:1000)。 figure(1) subplot(211) plot(x) title(39。原始信號 39。) t=wpdec(x,3,39。db139。,39。shannon39。)。 plot(t) [wpt,rn]=wpcutree(t,2)。 plot(wpt) rn=3， 4， 5， 6 二十一、 besttree 功能：計算最佳樹 格式及說明： BESTTREE putes the optimal subtree of an initial tree with respect to an entropy type criterion. The resulting tree may be much smaller than the initial one. T = BESTTREE(T) putes the modified tree T corresponding to the best entropy value. [T,E] = BESTTREE(T) returns the best tree T and in addition, the best entropy value E. The optimal entropy of the node whose index is j1 is E(j). [T,E,N] = BESTTREE(T) returns the best tree T, entropy value E and in addition, the vector N containing the indices of the merged nodes. 舉例 load noisdopp。 x=noisdopp。 t=wpdec(x,3,39。db139。,39。shannon39。)。 x1=wprec(t)。 [wpt,wpca,wpcd]=wpsplt(t,[3,0])。 plot(wpt) [bt,e,n]=besttree(wpt)。 plot(bt) e39。 ans = +004 * Columns 1 through 8 Columns 9 through 11 n＝ 2 二十二、 bestlevt 計算完整最佳小波包樹 BESTLEVT putes the optimal plete subtree of an initial tree with respect to an entropy type criterion. The resulting plete tree may be of smaller depth than the initial one. T = BESTLEVT(T) putes the modified tree T corresponding to the best level tree deposition. [T,E] = BESTLEVT(T) returns best tree T and in addition, the best entropy value E. The optimal entropy of the node whose index is j1 is E(j). load noisdopp。 x=noisdopp。 t=wpdec(x,3,39。db139。,39。shannon39。)。 [wpt,wpca,wpcd]=wpsplt(t,[3,0])。 plot(wpt) [blt,e]= bestlevt(wpt)。 plot(blt) e39。 ans = +004 * Columns 1 through 8 Columns 9 through 15 二十三、 wp2wtree 從小波包樹中提取小波樹 WP2WTREE Extract wavelet tree from wavelet packet tree. t = wp2wtree(t) putes the modified tree T corresponding to the wavelet deposition tree. load noisdopp。 x=noisdopp。 t=wpdec(x,4,39。db139。,39。shannon39。)。 plot(t) wt=wp2wtree(t)。 plot(wt)