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

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

波变换

开发平台：

Matlab

dscedw1d.m：源码内容
  1. function varargout = dscedw1d(varargin)
  2. %DSCEDW1D Demonstrates typical wavelet 1-D scenarios in the Wavelet Toolbox (Auto play). 
  3. %
  4. % This is a slideshow file for use with wshowdrv.m
  5. % To see it run, type 'wshowdrv dscedw1d', 
  7. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
  8. %   Last Revision: 03-Feb-2003.
  9. %   Copyright 1995-2004 The MathWorks, Inc.
  10. % $Revision: 1.13.4.2 $
  12. % Initialization and Local functions if necessary.
  13. if nargin>0
  14. action = varargin{1};
  15. switch action
  16.   case 'auto'    , wshowdrv('#autoMode',mfilename,'close');
  17.   case 'gr_auto' , wshowdrv('#gr_autoMode',mfilename,'close');
  19.   case 'getFigParam'
  20. figName  = 'Short Scenario DW1D';
  21. showType = 'command';
  22. varargout = {figName,showType};
  24.   case 'slidePROC'
  25.   wait_mode = 0;
  26.   [figHandle,idxSlide]  = deal(varargin{2:end});
  27.   idxPREV = wshowdrv('#get_idxSlide',figHandle);
  28.   localPARAM = wtbxappdata('get',figHandle,'localPARAM');
  29.   switch idxSlide
  30.    case 1
  31. if ~isempty(localPARAM)
  32. win_dw1d = localPARAM{1};
  33. delete(win_dw1d);
  34. dmsgfun('close');
  35. wtbxappdata('del',figHandle,'localPARAM');
  36. end
  37. % wshowdrv('#bloc_BackBtn',figHandle);
  38. return
  40. case 2
  41. if ~isempty(localPARAM)
  42. win_dw1d = localPARAM{1};
  43. delete(win_dw1d);
  44. end
  45. %-------------------------------------------
  46. % Display main window (demo initialization).
  47. %-------------------------------------------
  48. win_dw1d = dw1dtool;
  49. %----------------------------------------------
  50. % Main window parameters.
  51. %----------------------------------------------
  52. tag_pop_viewm  = 'View_Mode';
  53. tag_valapp_scr = 'ValApp_Scr';
  54. pop_handles    = findobj(win_dw1d,'style','popupmenu');
  55. pop_viewm      = findobj(pop_handles,'tag',tag_pop_viewm);
  56. pop_app_scr    = findobj(pop_handles,'tag',tag_valapp_scr);
  57. %------------------------------------------------
  58. localPARAM = {win_dw1d,pop_viewm,pop_app_scr};
  59. wtbxappdata('set',figHandle,'localPARAM',localPARAM);
  60. wshowdrv('#modify_cbClose',figHandle,win_dw1d,'dw1dtool');
  61. msg = get_Message(idxSlide);
  62. wshowdrv('#disp_msg',figHandle,msg,'step',localPARAM{1});
  63. if idxPREV<idxSlide
  64. figMSG = dmsgfun('handle');
  65. wshowdrv('#modify_cbClose',figHandle,figMSG,'dmsgfun');
  66. end
  67. return
  68.   
  69.     otherwise
  70. [win_dw1d,pop_viewm,pop_app_scr] = deal(localPARAM{:});
  71.   end
  73.   switch idxSlide
  74.     case 3
  75. if idxPREV<idxSlide
  76. %--------------------------
  77. % Full decomposition mode.
  78. %--------------------------
  79. sig_nam = 'noisdopp'; wav_nam = 'sym4'; lev_dec = 5;
  80. dw1dmngr('demo',win_dw1d,sig_nam,wav_nam,lev_dec);
  81. wenamngr('Inactive',win_dw1d);
  82. else
  83. set(pop_viewm,'value',2);
  84. eval(get(pop_viewm,'Callback'));
  85. end
  87.     case 4
  88. if idxPREV<idxSlide
  89. %----------------
  90. % Separate mode.
  91. %----------------
  92. set(pop_viewm,'value',3);
  93. eval(get(pop_viewm,'Callback'));
  94. else
  95. win_mopt_Params = wtbxappdata('get',figHandle,'win_mopt_Params');
  96. [win_mopt,chk_det_on,chk_cfs_on,pop_ccfs] = deal(win_mopt_Params{:});
  97. figure(win_mopt);
  98. set(chk_det_on,'Value',1);
  99. pause(0.5);
  100. dw1ddisp('apply',win_dw1d,win_mopt);
  101. set(pop_app_scr,'Value',1);
  102. eval(get(pop_app_scr,'Callback'));
  103. delete(win_mopt)
  104. wtbxappdata('del',figHandle,'win_mopt_Params');
  105. end
  107.     case 5
  108. if idxPREV<idxSlide
  109. %--------------------------------------
  110. % Show and Scroll mode (app5+sig+cfs).
  111. %--------------------------------------
  112. set(pop_viewm,'value',1);
  113. eval(get(pop_viewm,'Callback'));
  114. pause(2);
  115. %----------------------------------------
  116. % More Display Option window parameters.
  117. %----------------------------------------
  118. win_mopt_Params = makeWinOpt(figHandle,win_dw1d);
  119. [win_mopt,chk_det_on,chk_cfs_on,pop_ccfs] = deal(win_mopt_Params{:});
  120. %------------------------------------------------
  121. pause(1);
  122. set(chk_det_on,'Value',0);
  123. pause(1);
  124. dw1ddisp('apply',win_dw1d,win_mopt);
  125. set(pop_app_scr,'Value',5);
  126. eval(get(pop_app_scr,'Callback'));
  127. else
  128. win_mopt_Params = wtbxappdata('get',figHandle,'win_mopt_Params');
  129. [win_mopt,chk_det_on,chk_cfs_on,pop_ccfs] = deal(win_mopt_Params{:});
  130. figure(win_mopt);
  131. pause(1)
  132. set(pop_ccfs,'Value',1);
  133. pause(1)
  134. dw1ddisp('apply',win_dw1d,win_mopt);
  135. end
  137.     case 6
  138. %--------------------------------------------
  139. % Show and Scroll mode (app5+sig+cfs ).
  140. %     coloration : init + all levels + abs
  141. %--------------------------------------------
  142. if idxPREV>idxSlide
  143. makeWinOpt(figHandle,win_dw1d);
  144. return;
  145. end
  146. win_mopt_Params = wtbxappdata('get',figHandle,'win_mopt_Params');
  147. [win_mopt,chk_det_on,chk_cfs_on,pop_ccfs] = deal(win_mopt_Params{:});
  148. figure(win_mopt);
  149. pause(1)
  150. set(pop_ccfs,'Value',3);
  151. pause(1)
  152. dw1ddisp('apply',win_dw1d,win_mopt);
  154.     case 7
  155. %-----------------------------------
  156. % Close More Display Option window.
  157. % And message for Histograms.
  158. %-----------------------------------
  159. if idxPREV<idxSlide
  160. win_mopt_Params = wtbxappdata('get',figHandle,'win_mopt_Params');
  161. win_mopt = win_mopt_Params{1};
  162. dw1ddisp('close',win_dw1d,win_mopt);
  163. wtbxappdata('del',figHandle,'win_mopt_Params');
  164. wenamngr('Inactive',win_dw1d);
  165. else
  166. win_hist_Params = wtbxappdata('get',figHandle,'win_hist_Params');
  167. win_hist = win_hist_Params{1};
  168. dw1dhist('close',win_hist);
  169. delete(win_hist)
  170. wtbxappdata('del',figHandle,'win_hist_Params');
  171. end
  173. case 8
  174. %--------------------------------------------------
  175. % Histograms window - all the detail coefficients.
  176. % And first message for De-noising.
  177. %--------------------------------------------------
  178. win_hist      = dw1dmngr('hist',win_dw1d);
  179. wenamngr('Inactive',win_hist);
  180. tag_det_sig   = 'Det_sig';
  181. tag_det_all   = 'Det_All';
  182. tag_sel_cfs   = 'Sel_Cfs';
  183. tag_show_hist = 'Show_Hist';
  184. uic_win_hist  = findobj(win_hist,'Type','uicontrol');
  185. chk_win_hist  = findobj(uic_win_hist,'Style','checkbox');
  186. pus_win_hist  = findobj(uic_win_hist,'Style','pushbutton');
  187. chk_det_sig   = findobj(chk_win_hist,'tag',tag_det_sig);
  188. pus_det_all   = findobj(pus_win_hist,'tag',tag_det_all);
  189. rad_sel_cfs   = findobj(uic_win_hist,'Style','radiobutton','tag',tag_sel_cfs);
  190. pus_show_hist = findobj(pus_win_hist,'tag',tag_show_hist);
  191. win_hist_Params = {win_hist,chk_det_sig,pus_det_all,rad_sel_cfs,pus_show_hist};
  192. wtbxappdata('set',figHandle,'win_hist_Params',win_hist_Params);
  193. %--------------------------------------------------
  194. set(chk_det_sig,'value',1)      
  195. eval(get(chk_det_sig,'Callback')); pause(1)
  196. eval(get(pus_det_all,'Callback')); pause(1)
  197. eval(get(rad_sel_cfs,'Callback')); pause(1)
  198. eval(get(pus_show_hist,'Callback'));
  199. wenamngr('Inactive',win_hist);
  201.     case 9
  202. if idxPREV<idxSlide
  203. %------------------------------------
  204. % Close Histograms window.
  205. % And first message for Compression.
  206. %------------------------------------
  207. deleteSubFig(figHandle,'win_hist_Params');
  208. wenamngr('Inactive',win_dw1d);
  209. else
  210. deleteSubFig(figHandle,'win_comp_Params');
  211. end
  213.     case 10
  214. if idxPREV>idxSlide
  215. deleteSubFig(figHandle,'win_comp_Params');
  216. deleteSubFig(figHandle,'win_deno_Params');
  217. end
  219. %----------------------------------------
  220. % Compression: Initialization of  window.
  221. %----------------------------------------
  222. win_comp = dw1dmngr('comp',win_dw1d);
  223. wenamngr('Inactive',win_comp);
  224. % UTTHRGBL handlesUIC (Compression).
  225. %---------------------------------------
  226. % ud.handlesUIC = ...
  227. %   [fra_utl;txt_top;pop_met;...
  228. %    txt_sel;sli_sel;edi_sel; ...
  229. %    txt_nor;edi_nor;txt_npc;...
  230. %    txt_zer;edi_zer;txt_zpc;...
  231. %    tog_res;pus_est];
  232. %---------------------------------------
  233. % Memory of stored values (Compression).
  234. %---------------------------------------
  235. n_misc_comp = ['MB1_dw1dcomp'];
  236. ind_status  = 2;
  237. ind_pop_mod = 8;
  238. pop_mod    = wmemtool('rmb',win_comp,n_misc_comp,ind_pop_mod);
  239. handlesUIC = utthrgbl('handlesUIC',win_comp);
  240. edi_zer    = handlesUIC(11);
  241. win_comp_Params = {win_comp,pop_mod,edi_zer,n_misc_comp,ind_status};
  242. wtbxappdata('set',figHandle,'win_comp_Params',win_comp_Params);
  243. %--------------------------------------------------
  245. case 11
  246. %----------------------------------------------------
  247. % Compression: Global mode.
  248. % After compression with the default threshold value.
  249. %----------------------------------------------------
  250. win_comp_Params = wtbxappdata('get',figHandle,'win_comp_Params');
  251. win_comp = win_comp_Params{1};
  252. dw1dcomp('compress',win_comp,win_dw1d);
  253. wenamngr('Inactive',win_comp);
  255. case 12
  256. %--------------------------------------
  257. % Compression: change the zeros value.
  258. %--------------------------------------
  259. win_comp_Params = wtbxappdata('get',figHandle,'win_comp_Params');
  260. [win_comp,pop_mod,edi_zer,n_misc_comp,ind_status] = deal(win_comp_Params{:});
  261. set(edi_zer,'String',num2str(96.4));
  262. eval(get(edi_zer,'Callback'));
  263. dw1dcomp('compress',win_comp,win_dw1d);
  264. wenamngr('Inactive',win_comp);
  266. case 13
  267. %----------------------------------------------------
  268. % Compression: By Level mode.
  269. %----------------------------------------------------
  270. win_comp_Params = wtbxappdata('get',figHandle,'win_comp_Params');
  271. [win_comp,pop_mod,edi_zer,n_misc_comp,ind_status] = deal(win_comp_Params{:});
  272. set(pop_mod,'Value',2)
  273. eval(get(pop_mod,'Callback'))
  274. wenamngr('Inactive',win_comp);
  276. case 14
  277. %-----------------------------------------------------
  278. % Compression: By Level mode.
  279. % After compression with the default threshold values.
  280. %-----------------------------------------------------
  281. win_comp_Params = wtbxappdata('get',figHandle,'win_comp_Params');
  282. win_comp = win_comp_Params{1};
  283. dw1dcomp('compress',win_comp,win_dw1d);
  284. wenamngr('Inactive',win_comp);
  286.     case 15
  287. %----------------------------
  288. % Close Compression window.
  289. % And message for De-noising.
  290. %----------------------------
  291. deleteSubFig(figHandle,'win_comp_Params');
  292. wenamngr('Inactive',win_dw1d);
  293.   
  294.     case 16
  295. if idxPREV>idxSlide
  296. deleteSubFig(figHandle,'win_deno_Params');
  297. lin_handles = findobj(win_dw1d,'type','line');
  298. try
  299. lin(1) = findobj(lin_handles,'tag','SSig_in_App');
  300. lin(2) = findobj(lin_handles,'tag','SSig_in_Det');
  301. delete(lin);
  302. end
  303. end
  305. %--------------------------------------------------
  306. % De-noising: 
  307. %   - Initialization of  window.
  308. %   - De-noising with the default threshold values,
  309. %     with Fixed form soft thresholding method.
  310. %--------------------------------------------------
  311. win_deno = dw1dmngr('deno',win_dw1d);
  312. wenamngr('Inactive',win_deno);
  313. % UTTHRW1D handlesUIC (De-noising).
  314. %-------------------------------------
  315. % ud.handlesUIC = ...
  316. %   [fra_utl;txt_top;pop_met; ...
  317. %    rad_sof;rad_har;txt_noi;pop_noi; ...
  318. %    txt_BMS;sli_BMS;txt_tit(1:4),...
  319. %    txt_nor;edi_nor;txt_npc; ...
  320. %    txt_zer;edi_zer;txt_zpc; ...
  321. %    tog_thr;tog_res;pus_est];
  322. %--------------------------------------------------
  323. handlesUIC = utthrw1d('handlesUIC',win_deno);
  324. pop_met    = handlesUIC(3);
  325. rad_har    = handlesUIC(5);
  326. win_deno_Params = {win_deno,pop_met,rad_har};
  327. wtbxappdata('set',figHandle,'win_deno_Params',win_deno_Params);
  328. %--------------------------------------------------
  329. pause(1)
  330. dw1ddeno('denoise',win_deno,win_dw1d);
  331. wenamngr('Inactive',win_deno);
  333. case 17
  334. %--------------------------------------------------
  335. % De-noising: 
  336. %   - De-noising with the default threshold values.
  337. %     with Fixed form hard thresholding.
  338. %--------------------------------------------------
  339. win_deno_Params = wtbxappdata('get',figHandle,'win_deno_Params');
  340. [win_deno,pop_met,rad_har] = deal(win_deno_Params{:});
  341. if idxPREV>idxSlide
  342. set(pop_met,'Value',1)
  343. eval(get(pop_met,'Callback'))
  344. end
  345. set(rad_har,'Value',1)
  346. eval(get(rad_har,'Callback'))
  347. dw1ddeno('denoise',win_deno,win_dw1d);
  348. wenamngr('Inactive',win_deno);
  350. case 18
  351. % De-noising: 
  352. %   - Change the method of thresholding to:
  353. %     the Penalize medium thresholding method.
  354. %--------------------------------------------------
  355. win_deno_Params = wtbxappdata('get',figHandle,'win_deno_Params');
  356. [win_deno,pop_met,rad_har] = deal(win_deno_Params{:});
  357. set(pop_met,'Value',6)
  358. eval(get(pop_met,'Callback'))
  360. case 19
  361. if idxPREV<idxSlide
  362. % De-noising: 
  363. %   - De-noising with the default threshold values,
  364. %     with Penalize medium thresholding method.
  365. %--------------------------------------------------
  366. win_deno_Params = wtbxappdata('get',figHandle,'win_deno_Params');
  367. win_deno = win_deno_Params{1};
  368. dw1ddeno('denoise',win_deno,win_dw1d);
  369. wenamngr('Inactive',win_deno);
  370. else
  371. dmsgfun('close');
  372. end
  374. case 20 
  375. if idxPREV>idxSlide
  376. deleteSubFig(figHandle,'win_deno_Params');
  377. win_deno = dw1dmngr('deno',win_dw1d);
  378. wenamngr('Inactive',win_deno);
  379. handlesUIC = utthrw1d('handlesUIC',win_deno);
  380. pop_met    = handlesUIC(3);
  381. rad_har    = handlesUIC(5);
  382. win_deno_Params = {win_deno,pop_met,rad_har};
  383. wtbxappdata('set',figHandle,'win_deno_Params',win_deno_Params);
  384. set(pop_met,'Value',6)
  385. eval(get(pop_met,'Callback'))
  386. dw1ddeno('denoise',win_deno,win_dw1d);
  387. wenamngr('Inactive',win_deno);
  388. else
  389. % DISPLAY MESSAGE ONLY!
  390. end
  392. case 21
  393. %-------------------------------
  394. % Close De-noising window
  395. % and update synthesized signal.
  396. %-------------------------------
  397. win_deno_Params = wtbxappdata('get',figHandle,'win_deno_Params');
  398. win_deno = win_deno_Params{1};
  399. % MemBloc1 of stored values.
  400. %---------------------------
  401. n_param_anal  = 'DWAn1d_Par_Anal';
  402. ind_ssig_type = 7;
  403. wwaiting('msg',win_deno,'Wait ... computing');
  404. lin_den = utthrw1d('get',win_deno,'handleTHR');
  405. wmemtool('wmb',win_dw1d,n_param_anal,ind_ssig_type,'ds');
  406. dw1dmngr('return_deno',win_dw1d,1,lin_den);
  407. wenamngr('Inactive',win_dw1d);
  408. wwaiting('off',win_deno);
  409. delete(win_deno)
  410. dmsgfun('close');
  412. case 22
  413. if idxPREV<idxSlide
  414. %----------------------------------------
  415. % More Display Option window parameters.
  416. %----------------------------------------
  417. win_mopt_Params = makeWinOpt(figHandle,win_dw1d);
  418. win_mopt = win_mopt_Params{1};
  419. wenamngr('Inactive',win_mopt);
  420. %-----------------------------------------
  421. else
  422. win_mopt_Params = wtbxappdata('get',figHandle,'win_mopt_Params');
  423. [win_mopt,chk_det_on,chk_cfs_on,pop_ccfs] = deal(win_mopt_Params{:});
  424. set(chk_cfs_on,'Value',1);
  425. dw1ddisp('apply',win_dw1d,win_mopt);
  426. end
  428.     case 23
  429. if idxPREV<idxSlide
  430. win_mopt_Params = wtbxappdata('get',figHandle,'win_mopt_Params');
  431. else
  432. win_mopt_Params = makeWinOpt(figHandle,win_dw1d);
  433. end
  434. [win_mopt,chk_det_on,chk_cfs_on,pop_ccfs] = deal(win_mopt_Params{:});
  435. set(chk_cfs_on,'Value',0);
  436. pause(1)
  437. dw1ddisp('apply',win_dw1d,win_mopt);
  439. case 24
  440. %-----------------------------------
  441. % Close More Display Option window.
  442. %-----------------------------------
  443. deleteSubFig(figHandle,'win_mopt_Params');
  444.   end
  445.   [msg,ok] = get_Message(idxSlide);
  446.   if ok , wshowdrv('#disp_msg',figHandle,msg,'step',win_dw1d); end
  447. end
  448. return
  449. end
  451. if nargout<1,
  452.   wshowdrv(mfilename)
  453. else
  454.   idx = 0; slide(1).code = {}; slide(1).text = {};
  455.   
  456.   %========== Slide 1 ==========
  457.   idx = idx+1;
  458.   slide(idx).code = {
  459.   'figHandle = gcf;',
  460.   [mfilename ,'(''slidePROC'',figHandle,',int2str(idx),');']};
  462.   %========== Slide 2 to Slide 24 ==========
  463.   for idx = 2:24
  464.     slide(idx).code = {[mfilename ,'(''slidePROC'',figHandle,',int2str(idx),');']};
  465.   end
  467.   %========== Set "previous" slide indices ==========
  468.   
  469.   % Compression slides.
  470.   %--------------------
  471.   slide(12).idxPrev = 10;
  472.   slide(13).idxPrev = 10;
  473.   slide(14).idxPrev = 10;
  474.   slide(15).idxPrev = 10;
  475.   
  476.   % De-noising slides.
  477.   %-------------------
  478.   slide(16).idxPrev = 10;
  479.   slide(21).idxPrev = 20;
  480.   slide(22).idxPrev = 20;
  481.   %===================================================
  483.   varargout{1} = slide;
  484.   
  485. end
  488. %------------------------------------------------------------------------------------------%
  489. function deleteSubFig(figHandle,win_Params_Name)
  491. win_Params = wtbxappdata('get',figHandle,win_Params_Name);
  492. if ~isempty(win_Params)
  493. win = win_Params{1};
  494. if ishandle(win) , delete(win); end
  495. wtbxappdata('del',figHandle,win_Params_Name);
  496. end
  497. %------------------------------------------------------------------------------------------%
  498. function win_mopt_Params = makeWinOpt(figHandle,win_dw1d)
  500. win_mopt = dw1ddisp('create',win_dw1d);
  501. wenamngr('Inactive',win_mopt);
  502. tag_cfs_on   = 'Chk_Cfs';
  503. tag_det_on   = 'Chk_Det';
  504. tag_pop_ccfs = 'Pop_CCfs';
  505. uic_win_mopt = findobj(win_mopt,'Type','uicontrol');
  506. chk_det_on   = findobj(uic_win_mopt,'Tag',tag_det_on);
  507. chk_cfs_on   = findobj(uic_win_mopt,'Tag',tag_cfs_on);
  508. pop_ccfs     = findobj(uic_win_mopt,'Tag',tag_pop_ccfs);
  509. win_mopt_Params = {win_mopt,chk_det_on,chk_cfs_on,pop_ccfs};
  510. wtbxappdata('set',figHandle,'win_mopt_Params',win_mopt_Params);
  511. %------------------------------------------------------------------------------------------%
  512. function [msg,ok_msg] = get_Message(idxSlide)
  514. ok_msg = 1;
  515. switch idxSlide
  516. case 2
  517. str1 = 'This demo shows the 1-D wavelet capabilities of the ';
  518. str2 = 'Wavelet GUI tools. The scenario supposes that we want to';
  519. str3 = 'de-noise or compress a 1-D signal. In this example, we ';
  520. str4 = 'decompose a simple synthetic noisy Doppler signal using ';
  521. str5 = 'the Symlet 4 wavelet (sym4) into five levels.';
  522. msg  = strvcat(str1,str2,str3,str4,str5);
  524. case 3
  525. str1 = 'The full decomposition appears, with the signal denoted s at the top.';
  526. str2 = 'This signal is equal to the sum of the coarsest approximation (a5)';
  527. str3 = 'and the details at levels 1 to 5 (d1 ... d5).';
  528. str4 = 'i.e.,  s = a5 + d5 + d4 + d3 + d2 + d1';
  529. str5 = 'Since sym4 is an orthogonal wavelet, the decomposition is orthogonal.';
  530. msg  = strvcat(str1,str2,str3,str4,str5,'$');
  532. str1 = 'How do we select the wavelet and decomposition level?';
  533. msg  = strvcat(msg,str1,'$');
  535. str1 = 'For de-noising or compression, both the decomposition level and ';
  536. str2 = 'wavelet are useful tuning parameters.';
  537. str3 = 'Recall that the de-noised or compressed signal is equal to';
  538. str4 = 'the coarsest approximation + some signal reconstructed using the ';
  539. str5 = 'thresholded detail coefficients.';
  540. msg  = strvcat(msg,str1,str2,str3,str4,str5,'$');
  542. str1 = 'Often the wavelet choice is not crucial, and any "short" wavelet';
  543. str2 = 'except Haar is a good starting point.';
  544. msg  = strvcat(msg,str1,str2,'$');
  546. str1 = 'Having selected a wavelet, the decomposition level needs to be selected.';
  547. str2 = 'To study this, let us change to the Separate display mode.';
  548. msg  = strvcat(msg,str1,str2);
  550. case 4
  551. str1 = 'As a rule of thumb, we can select the decomposition level in such';
  552. str2 = 'a way that the corresponding approximation is noise free. In this';
  553. str3 = 'example 4 or 5 are good candidates, since the approximations at levels';
  554. str4 = '1 to 3 are clearly noisy. In addition, the detail d5 contains';
  555. str5 = 'significant noise as it can be seen inspecting the second half of the';
  556. str6 = 'detail, so 5 should be a good value for the decomposition level.';
  557. msg  = strvcat(str1,str2,str3,str4,str5,str6,'$');
  559. str1 = 'The Wavelet GUI provides a flexible and configurable interface.';
  560. msg  = strvcat(msg,str1,'$');
  562. str1 = 'In the current view, one can see the wavelet coefficients, cfs. We';
  563. str2 = 'focus on the coefficients, selecting Show and Scroll display mode,';
  564. str3 = 'selecting the approximation at level 5 and using the More Display';
  565. str4 = 'Options and removing the Detail Axes.';
  566. msg  = strvcat(msg,str1,str2,str3,str4);
  568. case 6
  569. str1 = 'The time-scale pattern shows that the';
  570. str2 = 'deterministic part of the signal is concentrated in a few big coefficients,';
  571. str3 = 'and noise corrupts all the time-scale representation.';
  572. str4 = 'To reinforce this, use the coloration  mode init + all levels + abs';
  573. str5 = 'in order to rescale globally the absolute coefficients.';
  574. msg  = strvcat(str1,str2,str3,str4,str5);
  576. case 7
  577. str1 = 'Clicking on the Histograms button and selecting all the detail coefficients,';
  578. str2 = 'we can see that the detail coefficients are approximately Gaussian N(0,1)';
  579. str3 = 'except for extreme values outside the interval [-4 4] which are';
  580. str4 = 'linked to the decomposition of the deterministic part of the signal.';
  581. msg  = strvcat(str1,str2,str3,str4);
  583. case 8
  584. str1 = 'The de-noising or compression procedure acts on these coefficients killing';
  585. str2 = 'the lowest coefficients in absolute value terms and keeping large values.';
  586. msg  = strvcat(str1,str2);
  588. case 9
  589. str1 = 'Compression.';
  590. str2 = '------------';
  591. str3 = 'We press the Compress button, and the compression window appears.';
  592. str4 = 'This allows the user to select the threshold value with respect to';
  593. str5 = 'two design parameters:';
  594. msg  = strvcat(str1,str2,str3,str4,str5);
  596. case 10
  597. str1 = '- the retained energy in percent which is equal to 100 times';
  598. str2 = 'the square norm of the compressed signal over the square norm of the';
  599. str3 = 'original signal (magenta line on the left);';
  600. msg  = strvcat(str1,str2,str3,'$');
  601. str4 = '- the number of zeros in percent which is equal to the percentage of zeros of';
  602. str5 = 'the wavelet representation of the compressed signal (blue line on the left).';
  603. str6 = 'The threshold default value given by the toolbox in the 1-D case is the one';
  604. str7 = 'for which the two previously defined design parameters are approximately the';
  605. str8 = 'same.';
  606. msg  = strvcat(msg,str4,str5,str6,str7,str8,'$');
  607. str1 = 'Then compression using the threshold default value leads a smoothed signal,';
  608. str2 = 'and the explanation of how the compression works is clear from the';
  609. str3 = 'comparison of original and thresholded coefficients in the time-scale';
  610. str4 = 'representations.';
  611. msg  = strvcat(msg,str1,str2,str3,str4);
  613. case 12
  614. str1 = 'You can then iterate choosing the threshold value or the retained energy';
  615. str2 = 'or the number of zeros. A modification of each of the three parameters';
  616. str3 = 'automatically updates the other two.';
  617. msg  = strvcat(str1,str2,str3);
  619. case 13
  620. str1 = 'To learn more about the compression procedure we select';
  621. str2 = 'the By-Level thresholding mode.';
  622. str3 = 'The detail coefficients are thresholded according to the thresholds represented';
  623. str4 = 'by the yellow lines on the left. With respect to the previous mode, this one';
  624. str5 = 'allows level-dependent threshold values leading to additional flexibility.';
  625. str6 = 'The resulting compressed signal can be found in the following figure.';
  626. msg  = strvcat(str1,str2,str3,str4,str5,str6);
  628. case 15
  629. str1 = 'De-noising.';
  630. str2 = '-----------';
  631. str3 = 'The design objective is to recover the signal or to remove noise,';
  632. str4 = 'without reference to an original interesting signal which is supposed to be';
  633. str5 = 'unknown. It is assumed that additive noise corrupts the signal to be';
  634. str6 = 'recovered, and is supposed to be stationary (at least at each scale).';
  635. msg  = strvcat(str1,str2,str3,str4,str5,str6,'$');
  637. str1 = 'The main ideas for de-noising and for compression are the same. The differences';
  638. str2 = 'are located at two points:';
  639. str3 = '- soft thresholding (shrinkage) is often used instead of the more crude';
  640. str4 = 'hard thresholding.';
  641. str5 = '- the noise structure (scaled or unscaled white noise, coloured noise) is of';
  642. str6 = 'interest.';
  643. msg  = strvcat(msg,str1,str2,str3,str4,str5,str6,'$');
  645. str1 = 'The GUI tools for de-noising are similar to the Compression tools. ';
  646. str2 = ' ';
  647. str3 = 'When you push on De-noise, the de-noising window appears,';
  648. str4 = 'we apply the default selection.';
  649. msg  = strvcat(msg,str1,str2,str3,str4);
  651. case 16
  652. str1 = 'As can be seen, the procedure keeps roughly the highest coefficients.';
  653. str2 = 'Then using the same default, we can show the differences between';
  654. str3 = 'soft and hard thresholding.';
  655. msg  = strvcat(str1,str2,str3);
  656. ok_msg = 1;
  658. case 177
  659. str1 = ' ';
  660. str2 = ' ';
  661. msg  = strvcat(str1,str2);
  663. case 20
  664. str1 = 'Now we close the de-noising window, update the de-noised signal and compare';
  665. str2 = 'the two signals in full size.';
  666. msg  = strvcat(str1,str2);
  668. otherwise
  669. msg = '';
  670. ok_msg = 0;
  671. end
  672. %------------------------------------------------------------------------------------------%