多小波变换矩阵形式.ppt

{PAGE} 因此, 若m是变换次数，A是一维变换矩阵，p是预滤波方法，T是m次变换时的置换矩阵，则运算T*A*p*x*p’*A’*T’就是对二维信号x做m次2维多小波变换, 若设B=T*A*p, 则可知经过m次变换后的小波系数矩阵为y=B*x*B’。下面是对16X16的信号x=(x[i,j])=(i+j)做1次小波变换的例子。此程序先使用ghm小波对x做1次变换，然后对x使用ghm小波，应用我们提出的方法做1次变换，并对2种变换结果相减后取绝对值 。从x4的结果几乎为零可看出，2种计算方法是等价的。 %程序名：mdwt_test8.m n=16;for i=1:n for j=1:n x(i,j)=i+j; end; end; x1=prep2D_appe(x,ghmap);x2=dec2D_pe(x1,ghm,1); [p,q]=pmatrix(n,ghmap);[L1,H1]=mdwt_matrix(n,ghm); A=[L1;H1];x3=smatrix(n,1)*A*p*x*p*A*smatrix(n,1)’; x4=x3-x2; sum(sum(abs(x4))) mdwt_test8 ans = 7.6014e-013 上面的程序mdwt_test8.m中的x有对称性，并且x是16X16的尺寸太小，与实际应用的尺寸有差别，现在我们将上面的程序中的x改动一