【文章內(nèi)容簡介】 bplot(626) plot(errx1) title(39。小波重構(gòu)誤差 39。) axis([0,1000,2e11,2e11])。 ?[lo_d,hi_d,lo_r,hi_r]=wfilters(39。db239。)。 ?[ca,cd]=dwt(s,lo_d,hi_d)。 ?x2=idwt(ca,cd,lo_r,hi_r)。 ?subplot(6,2,9)。 ?plot(x2)。 ?title(39。濾波器重構(gòu) 39。) ?errx2max=max(abs(sx2)) ?errx2=sx2。 ?subplot(6,2,10)。plot(errx2) ?title(39。濾波器重構(gòu)誤差 39。)。 ?axis([0,1000,2e11,2e11])。 小波濾波器 wfilters 格式 : (1)[LoD,HiD,LoR,HiR]=wfilters(‘wname’) (2)[f1,f2]=wfilters(‘wname’,’type’) [LO_D,HI_D,LO_R,HI_R] = WFILTERS(39。wname39。) putes four filters associated with the orthogonal or biorthogonal wavelet named in the string 39。wname39。. LO_D, the deposition lowpass filter HI_D, the deposition highpass filter LO_R, the reconstruction lowpass filter HI_R, the reconstruction highpass filter [F1,F2] = WFILTERS(39。wname39。,39。type39。) returns the following filters: LO_D and HI_D if 39。type39。 = 39。d39。 (Deposition filters) LO_R and HI_R if 39。type39。 = 39。r39。 (Reconstruction filters) LO_D and LO_R if 39。type39。 = 39。l39。 (Lowpass filters) HI_D and HI_R if 39。type39。 = 39。h39。 (Highpass filters) ’type’=‘d’ 分解濾波器 ’ type’=‘R’ 重構(gòu)濾波器 ’ type’=‘l’ 低通濾波器 ’ type’=‘h’ 高通濾波器 舉例 [lo_d,hi_d,lo_r,hi_r]=wfilters(39。haar39。)。 figure(1)； subplot(221)； stem(lo_d)； title(39。lod of haar39。)； subplot(222) stem(hi_d)； title(39。hid of haar39。) subplot(223)； stem(lo_r)； title(39。lor of haar39。) subplot(224) stem(hi_r) title(39。hir of haar39。) dwtmode 功能：離散小波變換拓展模式 格式： （ 1） dwtmode （ 2） dwtmode(‘mode’) 說明：當(dāng)對信號或圖像的邊緣進(jìn)行處理時，需要信號的邊緣進(jìn)行拓展。拓展模式有三種。該指令在進(jìn)行離散小波變換或小波包變換時，進(jìn)行模式拓展設(shè)定。 模式類型 類型說明 zdp 補零模式，缺省設(shè)定 sym 對稱延拓模式，即把邊緣值進(jìn)行復(fù)制 spd 平滑模式，對信號邊緣進(jìn)行某種平滑處理 wavedec 功能：多尺度一維小波分解（一維多分辨分析函數(shù)） 格式：（ 1） [c,l]=wavedec(x,n,’wname’) （ 2） [c,l]=wavedec(x,n,LoD,HiD) 用小波或分解濾波器對信號 X進(jìn)行一維多尺度分解 ，n為尺度和正整數(shù) 。 輸出參數(shù) c是由 組成 ， L是由 組成 。 11[ , , , , ]jjjc a c d c d c d?11[ , , , , ]jjjc a c d c d c d?的長度 的長度 的長度 的長度圖： xca 1 cd 1ca 2 cd 2ca 3 cd 3ca 3 cd 3 cd 2 cd 1C :L :ca 3 的長度 cd 3 的長度 cd 2 的長度 cd 1 的長度 x 的長度舉例 %多尺度一維離散小波變換； load sumsin。 s=sumsin。 subplot(611) plot(s)。 title(39。原始信號 39。) [c,l]=wavedec(s,3,39。db139。)。 subplot(613) plot(c)。 title(39。信號 s3尺度分解 39。)。 L= 125 125 250 500 1000 appcoef 功能：提取一維小波變換低頻系數(shù) 格式：（ 1） A＝ appcoef(c,l,’wname’,N) （ 2） A＝ appcoef(c,l,’wname’) （ 3） A＝ appcoef(c,l,LoR,HiR ) （ 4） A＝ appcoef(c,l,LoR,HiR ,N) 說明：該函數(shù)是一個一維小波分解函數(shù)，用于從 小波分解結(jié)構(gòu) [C,L]中提取一維信號的低頻 系數(shù)。 格式（ 1）計算尺度 N時的低頻系數(shù)， 格式（ 2）用于提取最后一個尺度的低頻系數(shù)， 格式（ 3）和（ 4）用濾波器提取低頻系數(shù)。 舉例 %提取一維小波變換低頻系數(shù)； load leleccum。 s=leleccum(1:2021) subplot(421) plot(s)。 title(39。原始信號 39。) [c,l]=wavedec(s,3,39。db139。)。 ca1=appcoef(c,l,39。db139。,1)。 subplot(445) plot(ca1) ylabel(39。ca139。)。 ca2=appcoef(c,l,39。db139。,2)。 subplot(4,8,17) plot(ca2)。 ylabel(39。ca239。)。 Detcoef 功能：提取一維信號小波變換高頻系數(shù) 格式：（ 1） d=detcoef(c,l,N) 提取 N尺度的高頻系數(shù)。 （ 2） d=detcoef(c,l)，提取最后一尺度的高頻系數(shù)。 舉例 %提取一維小波變換高頻系數(shù)； load leleccum。 s=leleccum(1:2021) subplot(421) plot(s)。 title(39。原始信號 39。) [c,l]=wavedec(s,3,39。db139。)。 cd1=detcoef(c,l,1)。 subplot(445) plot(cd1) ylabel(39。cd139。)。 cd2=detcoef(c,l,2)。 subplot(4,8,17) plot(cd2)。 ylabel(39。cd239。)。 九、 Waverec 功能：多尺度一維小波重構(gòu) 格式：（ 1） x=waverec(c,l,’wname’) （ 2） x=waverec(c,l,LoR,HiR) （ 3） x＝ waverec（ wavedec（ x,N,’wavename’),’ wavename’) 說明：該函數(shù)用指定的小波函數(shù)或重構(gòu)濾波器對 小波分解結(jié)構(gòu) (C,L)進(jìn)行多尺度一維小波重構(gòu)。 舉例 %多尺度一維小波重構(gòu)； load leleccum。 s=leleccum(1:3920) subplot(311) plot(s)。 title(39。原始信號 39。) [c,l]=wavedec(s,3,39。db539。)。 a=waverec(c,l,39。db539。) subplot(312) plot(a) title(39。重構(gòu)信號 39。) err=sa。 subplot(313) plot(err) title(39。誤差 39。) 十、 upwlev 功能：單尺度一維小波分解的重構(gòu) 格式： (1)[nc,nl,ca]=upwlev(c,l,’wname’) (2) [nc,nl,ca]=upwlev(c,l,LoR, HiR) 說明：該函數(shù)用于對小波分解結(jié)構(gòu) [C,L]進(jìn)行單尺度重構(gòu) ， 返回上一尺度的分解結(jié)構(gòu)并提取最后一尺度的低頻分量 。 %單尺度一維小波分解的重構(gòu)； load sumsin。 s=sumsin。 subplot(611) plot(s)。 title(39。原始信號 39。) [c,l]=wavedec(s,3,39。db139。)。 subplot(613) plot(c) title(39。尺度 3的小波分解結(jié)構(gòu) 39。) xlabel(39。尺度 3的低頻系數(shù)和尺度 1的高頻系數(shù) 39。) [nc,nl]=upwlev(c,l,39。db139。)。 subplot(615)。 plot(nc)。 title(39。尺度 2的小波分解結(jié)構(gòu) 39。) xlabel(39。尺度 2的低頻系數(shù)和尺度 1的高頻系數(shù) 39。) 等效于 [c,l]=wavedec(s,2,39。db139。)。 plot(c) NL=250 250 500 1000 L= 125 125 250 500 1000 十一、 Wrcoef 功能：對一維小波系數(shù)進(jìn)行單支重構(gòu) 格式： (1)x=wrcoef(‘type’,c,l,’wname’,N) (2)x=wrcoef(‘type’,c,l,LoR,HiR,N) (3)x=wrcoef(‘type’,c,l,’wname’) (4)x=wrcoef(‘type’,c,l,LoR,HiR) 說明：對一維信號的分解結(jié)構(gòu) [C,L]用指定的小波函數(shù)或重構(gòu)濾波器進(jìn)行重構(gòu) 。 當(dāng) ‘ type=a’時 ， 對信號的低頻部分進(jìn)行重構(gòu) ， 此時 N可以為 0；當(dāng)‘ type=d’時 ， 對信號的高頻部分進(jìn)行重構(gòu) ， 此時N為正整數(shù) 。 %對一維小波系數(shù)進(jìn)行單支重構(gòu)； load sumsin。 s=sumsin。 subplot(611) plot(s)。 title(39。原始信號 39。) [c,l]=wavedec(s,