小波函数及常用指令.ppt

十、 upwlev 功能：单尺度一维小波分解的重构 格式：(1)[nc,nl,ca]=upwlev(c,l,’wname’) (2) [nc,nl,ca]=upwlev(c,l,Lo-R, Hi-R) 说明：该函数用于对小波分解结构[C,L]进行单尺度重构，返回上一尺度的分解结构并提取最后一尺度的低频分量。 第六十一页，共一百一十七页，2022年，8月28日 %单尺度一维小波分解的重构； load sumsin; s=sumsin; subplot(611) plot(s); title(原始信号) [c,l]=wavedec(s,3,db1); subplot(613) plot(c) title(尺度3的小波分解结构) xlabel(尺度3的低频系数和尺度3、2、1的高频系数) [nc,nl]=upwlev(c,l,db1); subplot(615); plot(nc); title(尺度2的小波分解结构) xlabel(尺度2的低频系数和尺度2、1的高频系数) 等效于[c,l]=wavedec(s,2,db1); plot(c) 第六十二页，共一百一十七页，2022年，8月28日 NL=250 250 500 1000 L= 125 125 250 500 1000 第六十三页，共一百一十七页，2022年，8月28日 十一、Wrcoef 功能：对一维小波系数进行单支重构 格式：(1)x=wrcoef(‘type’,c,l,’wname’,N) (2)x=wrcoef(‘type’,c,l,Lo-R,Hi-R,,N) (3)x=wrcoef(‘type’,c,l,’wname’) (4)x=wrcoef(‘type’,c,l,Lo-R,Hi-R) 说明：对一维信号的分解结构[C,L]用指定的小波函数或重构滤波器进行重构。当‘type=a’时，对信号的低频部分进行重构，此时N可以为0；当‘type=d’时，对信号的高频部分进行重构，此时N为正整数。 第六十四页，共一百一十七页，2022年，8月28日 %对一维小波系数进行单支重构； load sumsin; s=sumsin; subplot(611) plot(s); title(原始信号) [c,l]=wavedec(s, 5,sym4); a5=wrcoef(a,c,l,sym4,5); subplot(613) plot(a5) title(低频部分重构信号) a51=wrcoef(d,c,l,sym4,5); subplot(615) plot(a51) title(高频部分重构信号) 第六十五页，共一百一十七页，2022年，8月28日 第六十六页，共一百一十七页，2022年，8月28日 十二、upcoef 功能：一维系数的直接小波重构 格式：（1）y=upcoef(‘0’,x,’wname’,N) （2） y=upcoef(‘0’,x,’wname’,N,L) （3） y=upcoef(‘0’,x,Lo-R,Hi-R,N) （4） y=upcoef(‘0’,x,Lo-R,Hi-R,N,L) （5） y=upcoef(‘0’,x,’wname’) （6） y=upcoef(‘0’,x,Lo-R,Hi-R) 第六十七页，共一百一十七页，2022年，8月28日 说明：该函数用于一维小波分析，它用来计算向量X（信号系数）向上N步的重构小波系数，N为正整数。如果0＝a，对低频系数进行重构；如果0＝d，对高频系数进行重构；对于（2）和（4），则是对向量X中间长度为L部分进行重构。 第六十八页，共一百一十七页，2022年，8月28日 Load leleccum; s= leleccum(1:2000); Plot(s) title(原始信号); [c,l]=wavedec(s,3,db6); ca1=appcoef(c,l,db6,1); sca1=upcoef(a,ca1,db6,1); subplot(622);plot(sca1); title(尺度1的低频系数ca1 向上一步重构信号); axis([0,2000,200,600]); sca1l=upcoef(a,ca1,db6,1,1000); subplot(625);plot(sca1l); title(‘ca1向上一步只取1000 点重构信号); axis([0,2000,200,600]); cd1=detcoef(c,l,1); scd