小波实验报告双树复小波变换

一、题目：双树复小波变换 二、目的：双树复小波和实小波变换的比较 三、算法及其实现：提取阶梯型边界点 1． 算法。 幅值: 相位： 2． 代码实现。 我采用 Matlab 函数编程实现。具体程序见 shift_test_2D.m，drawcirc.m，setfig.m，dtwavexfm2.m，dtwaveifm2.m,waveifm2.m,wavexfm2.m SkelMap.m。 设 和 分别是双正交对偶尺度函数与对偶小波， ， ， 和 是相应的低通滤波器和高通滤波器，即它们满足 实部： 虚部： 双树复小波变换可以通过离散小波变换 DWT 实现：一个 DWT 产生实部，另一个产生虚部。 四、实现工具：Matlab 五、程序代码： （1 ）shift_test_2D.m： % shift_test_2D.m % % M-file to perform a 4-level wavelet transform on a circle using Q-shift % dual wavelet tree and DWT, and to compare shift invariance properties. % % Nick Kingsbury, Cambridge University, May 2002. clear all close all % Draw a circular disc. x = round((drawcirc(64,1,0,0,256) - 0.5) * 200); setfig(1); colormap(gray(256)) image(min(max(x+128,1),256)); set(gca,'position',[0.1 0.25 .25 .5]); ( )( )j( )chgttt,,j,chgfff22|( , ) |[( , )][( , )]chgdj ndj ndj n( , )( , )arctan( , )gchdj ndj ndj n,hh,hh  0h n 0hn 1h n 1h n0101( )2( ) (2)( )2( ) (2)( )2( ) (2)( )2( ) (2)hhnhhnhhnhhnth nt nth nt nth nt nth nt n0101( )2( )(2)( )2( )(2)( )2( )(2)( )2( )(2)ggnggnggnggntg nt ntg nt ntg nt ntg nt naxis('off'); axis('image'); % draw(xx); title('Input (256 x 256)','FontSize',14); drawnow % Do 4 levels of CWT. [Yl,Yh] = dtwavexfm2(x,4,'near_sym_b','qshift_b'); % Loop to reconstruct output from coefs at each level in turn. % Starts with the finest level. titl = ['1st';'2nd';'3rd';'4th';'Low']; yy = zeros(size(x) .* [2 3]); yt1 = 1:size(x,1); yt2 = 1:size(x,2); for mlev = 1:5, mask = zeros(6,5);...