开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 波变换 > 查看源码
dcmddw1d.m
资源名称：wavelet_toolbox.rar　[点击查看]
上传用户：haiyisale
上传日期：2013-01-09
资源大小：3246k
文件大小：9k
源码类别：

波变换

开发平台：

Matlab

dcmddw1d.m：源码内容
  1. function varargout = dcmddw1d(varargin)
  2. %DCMDDW1D Demonstrates discrete 1-D wavelet tools in the Wavelet Toolbox.
  3. %
  4. % This is a slideshow file for use with wshowdrv.m
  5. % To see it run, type 'wshowdrv dcmddw1d', 
  6. %
  7. % See also APPCOEF, DETCOEF, DWT, IDWT, UPCOEF,
  8. %          WAVEDEC, WAVEREC, WRCOEF. 
  10. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
  11. %   Last Revision: 01-Jun-2001.
  12. %   Copyright 1995-2002 The MathWorks, Inc.
  13. %   $Revision: 1.14 $
  15. % Initialization and Local functions if necessary.
  16. if nargin>0
  17. action = varargin{1};
  18. switch action
  19.       case 'addHelp'
  20. % Add Help Item.
  21. %---------------
  22.         hdlFig = varargin{2};
  23. wfighelp('addHelpItem',hdlFig,'Wavelet Decomposition','DW1D_DECOMPOS');
  24.           
  25.   case 'auto'    , wshowdrv('#autoMode',mfilename,'close');
  26.   case 'gr_auto' , wshowdrv('#gr_autoMode',mfilename,'close');
  27.   case 'getFigParam'
  28. figName  = 'Orthogonal Wavelets 1-D';
  29. showType = 'mix6';
  30. varargout = {figName,showType};
  31. end
  32. return
  33. end
  35. if nargout<1,
  36.   wshowdrv(mfilename)
  37. else
  38.   idx = 0;
  39.   %========== Slide 1 ==========
  40.   idx = idx+1;
  41.   slide(idx).code = {
  42. 'figHandle = gcf;',
  43. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow; h = [];',   
  44. '' };
  46.   slide(idx).text = {
  47. '',
  48. ' Press the "Start" button to see a demonstration of',
  49. ' 1-D wavelet tools.',
  50. '',
  51. ' This demo uses Wavelet Toolbox functions.',
  52. ''};
  54.   %========== Slide 2 ==========
  55.   idx = idx+1;
  56.   slide(idx).code = {
  57. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',   
  58. 'load leleccum; s = leleccum(1:3920);',
  59. 'ls = length(s);',
  60. 'clear leleccum',
  61. 'h(1) = subplot(2,1,1); plot(1:ls,s,''r'')',
  62. 'title(''Original signal s.'');',
  63. 'set(h(1),''Xlim'',[1,ls]);',
  64. '' };
  66.   slide(idx).text = {
  67. ' % Load original 1D signal.',
  68. '',
  69. '        load leleccum; s = leleccum(1:3920);',
  70. '        ls = length(s);',
  71. ''};
  73.   %========== Slide 3 ==========
  74.   idx = idx+1;
  75.   slide(idx).code = {
  76. '[ca1,cd1] = dwt(s,''db1'');',
  78. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow; h = [];',   
  79. 'h(1) = subplot(2,1,1); plot(1:ls,s,''r'')',
  80. 'title(''Original signal s.'');',
  81. 'set(h(1),''Xlim'',[1,ls]);',
  83. 'h(3) = subplot(2,2,3); plot(1:length(ca1),ca1);',
  84. 'title(''Approx. coef. : ca1'');',
  85. 'set(h(3),''Xlim'',[1,length(ca1)]);',  
  86. 'h(4) = subplot(2,2,4); plot(1:length(cd1),cd1);',
  87. 'title(''Detail coef. : cd1'');',
  88. 'set(h(4),''Xlim'',[1,length(cd1)]);',   
  89. '' };
  91.   slide(idx).text = {
  92. ' % Perform one step decomposition of s using db1.',
  93. '',
  94.   '        [ca1,cd1] = dwt(s,''db1'');',
  95. ''};
  97.   slide(idx).info = 'dwt';
  99.   %========== Slide 4 ==========
  100.   idx = idx+1;
  101.   slide(idx).code = {
  102. 'a1 = upcoef(''a'',ca1,''db1'',1,ls);',
  103.     'd1 = upcoef(''d'',cd1,''db1'',1,ls);',
  105. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',   
  106. 'h(1) = subplot(3,1,1); plot(1:ls,s,''r'');',
  107. 'title(''Original signal s.'');',
  108. 'h(2) = subplot(3,1,2); plot(1:ls,a1);',
  109. 'title(''Approximation : a1.'');',
  110. 'h(3) = subplot(3,1,3); plot(1:ls,d1);',
  111. 'title(''Detail : d1.'');',
  112. 'set(h(1:2),''Xtick'',[],''Xticklabel'',[]);', 
  113. 'set(h,''Xlim'',[1,ls]);', 
  114. ''};
  116.   slide(idx).text = {
  117. ' % Perform one step reconstruction of  ca1 and cd1.',
  118. '',
  119. '        a1 = upcoef(''a'',ca1,''db1'',1,ls);',
  120. '        d1 = upcoef(''d'',cd1,''db1'',1,ls);',
  121. ''};
  123.   slide(idx).info = 'upcoef';
  125.   %========== Slide 5 ==========
  126.   idx = idx+1;
  127.   slide(idx).code = {
  128. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',   
  129. 'h(1) = subplot(3,1,1); plot(1:ls,s,''r'');',
  130. 'title(''Original signal s.'');',
  131. 'h(2) = subplot(3,1,2); plot(1:ls,a1+d1);',
  132. 'title(''Approximation + detail : a1  + d1.'');',
  133. 'set(h(1),''Xtick'',[],''Xticklabel'',[]);',
  134. 'set(h,''Xlim'',[1,ls]);', 
  135. ''};
  137.   slide(idx).text = {
  138. ' % Now plot a1 + d1.',
  139. ''};
  141.   slide(idx).info = 'upcoef';
  143.   %========== Slide 6 ==========
  144.   idx = idx+1;
  145.   slide(idx).code = {
  146. 'a0 = idwt(ca1,cd1,''db1'',ls);',
  147.   
  148. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',   
  149. 'h(1) = subplot(3,1,1); plot(1:ls,s,''r'');',
  150. 'title(''Original signal s.'');',
  151. 'h(2) = subplot(3,1,2); plot(1:ls,a1+d1);',
  152. 'title(''Approximation + detail : a1  + d1.'');',
  153. 'set(h(1:2),''Xtick'',[],''Xticklabel'',[]);',
  154.   
  155. 'h(3) = subplot(3,1,3); plot(1:ls,a0);',
  156. 'title(''Reconstructed signal : a0.'');',
  157. 'set(h,''Xlim'',[1,ls]);', 
  158. ''};
  160.   slide(idx).text = {
  161. ' % Invert directly decomposition of s using coefficients.',
  162. '',
  163. '        a0 = idwt(ca1,cd1,''db1'',ls);',
  164. ''};
  166.   slide(idx).info = 'idwt';
  168.   %========== Slide 7 ==========
  169.   idx = idx+1;
  170.   slide(idx).code = {
  171. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',   
  172. '[c,l] = wavedec(s,3,''db1'');',
  173. 'h(1) = subplot(5,1,1); plot(1:ls,s,''r'');',
  174. 'title(''Original signal s.'');',
  175. 'set(h,''Xlim'',[1,ls]);', 
  176. ''};
  178.   slide(idx).text = {
  179. ' % Perform decomposition at level 3 of s using db1.',
  180. '',
  181. '        [c,l] = wavedec(s,3,''db1'');',
  182. ''};
  184.   slide(idx).info = 'wavedec';
  186.   %========== Slide 8 ==========
  187.   idx = idx+1;
  188.   slide(idx).code = {
  189. 'ca3 = appcoef(c,l,''db1'',3);',      
  191. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',   
  192. 'h(1) = subplot(5,1,1); plot(1:ls,s,''r'');',
  193. 'set(h(1),''Xtick'',[],''Xticklabel'',[],''Xlim'',[1,ls]);', 
  194. 'h(2) = subplot(5,8,9); plot(1:length(ca3),ca3);',
  195. 'subplot(5,1,1); title(''Original signal s and ca3.'')',
  196. ''};
  198.   slide(idx).text = {
  199. ' % Extract approximation coefficients at level 3,',
  200. ' % from wavelet decomposition structure [c,l].',
  201. '',
  202. '        ca3 = appcoef(c,l,''db1'',3);',
  203. ''};
  205.   slide(idx).info = 'appcoef';
  207.   %========== Slide 9 ==========
  208.   idx = idx+1;
  209.   slide(idx).code = {
  210. 'cd3 = detcoef(c,l,3);',
  211. 'cd2 = detcoef(c,l,2);',      
  212. 'cd1 = detcoef(c,l,1);',
  214. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',   
  215. 'h(1) = subplot(5,1,1); plot(1:ls,s,''r'');',
  216. 'set(h(1),''Xtick'',[],''Xticklabel'',[],''Xlim'',[1,ls]);', 
  217. 'h(2) = subplot(5,8,9); plot(1:length(ca3),ca3);',
  219. 'subplot(5,1,1); title(''Original signal s , ca3 and cd3.'')',
  220. 'set(h(2),''Xtick'',[],''Xticklabel'',[]);', 
  221. 'h(3) = subplot(5,8,17); plot(1:length(cd3),cd3);',
  222. 'pause(1)',
  223. 'subplot(5,1,1); title(''Original signal s , ca3, cd3 and cd2.'')',
  224. 'set(h(3),''Xtick'',[],''Xticklabel'',[]);', 
  225. 'h(4) = subplot(5,4,13); plot(1:length(cd2),cd2);',
  226. 'pause(1)',
  227. 'subplot(5,1,1); title(''Original signal s , ca3, cd3, cd2 and cd1.'')',
  228. 'set(h(4),''Xtick'',[],''Xticklabel'',[]);', 
  229. 'h(5) = subplot(5,2,9); plot(1:length(cd1),cd1);',
  230. ''};
  232.   slide(idx).text = {
  233. ' % Extract detail coefficients at levels 1, 2 and 3,',
  234. ' % from wavelet decomposition structure [c,l].',
  235. '',
  236. '        cd3 = detcoef(c,l,3);',
  237. '        cd2 = detcoef(c,l,2);',
  238. '        cd1 = detcoef(c,l,1);',
  239. ''};
  241.   slide(idx).info = 'detcoef';
  243.   %========== Slide 10 ==========
  244.   idx = idx+1;
  245.   slide(idx).code = {
  246. 'a3 = wrcoef(''a'',c,l,''db1'',3);',
  248. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',   
  249. 'h(1) = subplot(5,1,1); plot(s,''r'');',
  250. 'title(''Original signal s and approximation of level 3.'');',
  251. 'h(2) = subplot(5,1,2); plot(1:ls,a3);',
  252. 'set(h(1),''Xtick'',[],''Xticklabel'',[]);',
  253. 'set(h,''Xlim'',[1,ls]);',
  254. ''};
  256.   slide(idx).text = {
  257. ' % Reconstruct approximation at level 3,',
  258.   ' % from wavelet decomposition structure [c,l].',
  259. '',
  260. '        a3 = wrcoef(''a'',c,l,''db1'',3);',
  261. ''};
  263.   slide(idx).info = 'wrcoef';
  265.   %========== Slide 11 ==========
  266.   idx = idx+1;
  267.   slide(idx).code = {
  268. 'd3 = wrcoef(''d'',c,l,''db1'',3);',
  269. 'd2 = wrcoef(''d'',c,l,''db1'',2);',
  270. 'd1 = wrcoef(''d'',c,l,''db1'',1);',
  272. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',   
  273. 'h(1) = subplot(5,1,1); plot(s,''r'');',
  274. 'title(''Original signal s and a3 d3 d2 d1.'');',
  275. 'h(2) = subplot(5,1,2); plot(1:ls,a3);',
  277. 'h(3) = subplot(5,1,3); plot(1:ls,d3);',
  278. 'h(4) = subplot(5,1,4); plot(1:ls,d2);',
  279. 'h(5) = subplot(5,1,5); plot(1:ls,d1);',
  280. 'set(h(1:4),''Xtick'',[],''Xticklabel'',[]);',
  281. 'set(h,''Xlim'',[1,ls]);',
  282. ''};
  284.   slide(idx).text = {
  285. ' % Reconstruct detail coefficients at levels 1,2 and 3,',
  286. ' % from the wavelet decomposition structure [c,l].',
  287. '',
  288. '        d3 = wrcoef(''d'',c,l,''db1'',3);',
  289. '        d2 = wrcoef(''d'',c,l,''db1'',2);',
  290. '        d1 = wrcoef(''d'',c,l,''db1'',1);',  
  291. ''};
  293.   slide(idx).info = 'wrcoef';
  295.   %========== Slide 12 ==========
  296.   idx = idx+1;
  297.   slide(idx).code = {
  298. 'a0 = waverec(c,l,''db1'');',
  300. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',   
  301. 'h(1) = subplot(2,1,1); plot(s,''r'');',
  302. 'title(''Original signal s .'');',
  303. 'h(2) = subplot(2,1,2); plot(1:ls,a0);',
  304. 'title(''Reconstructed signal a0.'');',
  305. 'set(h,''Xlim'',[1,ls]);',
  306. ''};
  308.   slide(idx).text = {
  309. ' % Reconstruct s from the wavelet decomposition structure [c,l].',
  310. '',
  311. '        a0 = waverec(c,l,''db1'');',
  312. ''};
  314.   slide(idx).info = 'waverec';
  316.   varargout{1} = slide;
  317. end