第5章_小波变换的matlab实现详解.ppt

小波分析示例 一维连续小波 1. coefs = cwt(s,scale,’wname’) 2. coefs = cwt(s,scale,’wname’,’plot’) c=cwt(noissin,1:48,db4,plot); 一维离散小波分解 系数重构 命令：upcoef 格式： 1. Y=upcoef(O,X,’wname’,N) 2. Y=upcoef(O,X,’wname’,N,L) 3. Y=upcoef(O,X,’Lo_R, Hi_R’,N) 4. Y=upcoef(O,X,’Lo_R,Hi_R’,N,L) 5. Y=upcoef(O,X,’wname’) 6. Y=upcoef(O,X,Lo_R,Hi_R) O=‘a’ 低频，O=‘d’ 高频 举例： A1=upcoef(a,cA1,db1,1,ls); D1=upcoef(d,cD1,db1,1,ls); subplot(1,2,1);plot(A1);title(Approximation A1) subplot(1,2,2);plot(D1);title(Detail D1) 逆变换恢复信号 命令：idwt 格式： 1. X=idwt(cA,cD,’wname’) 2. X=idwt(cA,cD,Lo_R, Hi_R) 3. X=idwt(cA,cD,’wname’,L) 4. X=idwt(cA,cD, Lo_R, Hi_R, L) 举例：A0=idwt(cA1,cD1,db1,ls); 多尺度一维分解 命令：wavedec 格式： [C, L]=wavedec(X,N,’wname’) [C, L]=wavedec(X,N,Lo_D,Hi_D) [C,L]=wavedec(s,3,db1); 低频系数提取 命令：appcoef 格式： 1. A=appcoef(C,L,’wname’,N) 2. A=appcoef(C,L,’wname’) 3. A=appcoef(C,L,Lo_R,Hi_R,N) 4. A=appcoef(C,L, Lo_R,Hi_R) 高频系数提取 命令：detcoef 格式： 1. A=detcoef(C,L,N) 2. A=detcoef(C,L) 举例 cA3=appcoef(C,L,db1,3); cD3=detcoef(C,L,3); cD2=detcoef(C,L,2); cD1=detcoef(C,L,1); 重构系数 命令：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) type=‘a’ 低频， type=‘d’ 高频 A3=wrcoef(a,C,L,db1,3); D1=wrcoef(d,C,L,db1,1); D2=wrcoef(d,C,L,db1,2); D3=wrcoef(d,C,L,db1,3); 重构原始信号 命令：waverec 格式： 1.X=waverec(C,L,’wname’) 2.X=waverec(C,L,Lo_R,Hi_R) 例子： A0=waverec(C,L,’db1’); 重构最大误差： Err=max(abs(s-A0)) 图形接口方式（GUI） 2. 二维离散小波 单尺度分解——dwt2 格式： 1. [cA1,cH1,cV1,cD1]=dwt2(X,’wname’) 2. [cA1,cH1,cV1,cD1]=dwt2(X,Lo_D,Hi_D) cA1, cH1水平;cV1垂直;cD1对角 应用： load wbarb; figure(1); image(X); colormap(map); colorbar; [cA1,cH1,cV1,cD1]=dwt2(X,’bior3.7’) 重构系数 命令：upcoef2 格式： 1. Y=upcoef2(O,X,’wname’,N,S) 2. Y=up