小波变换matlab程序,第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 格式：