一个小波变换实例Matlab实现.docxVIP

1 (t) ( ) ，使 (t k) k Z 为一组正交归一基； 、 选择 或 ? 2、 求 hn 。 hn (t), 1,n (t) 或H( ) ?(2 )/?( ) 3、 由 hn 求 gn 。 gn ( 1)n h1 n 或 G( ) e i H ( t ) 4、 由 gn ， (t) 构成正交小波基函数 (t ) (t) gn 1,n (t ) 或 ?( ) G( /2)?( /2) Haar 小波的构造 1）、选择尺度函数。 1 0 t 1 (t) 其他 易知 (t n) 关于 n 为一正交归一基。 2）、求 hn hn (t), 1,n (t) = 2 (t) （2t-n)dt 其中 n n 1 (2t n) 1 t 2 2 0 其他 当 n=0 时， 1 1 t (2t ) 2 其他 当 n=1 时， 1 1 t 1 (2 t 1) 2 0 其他 故，当 n=0， n=1 时 (t) ? (2t 1 n n n) 其他 0 当 n=0 时， 1 1 (t) ? (2t t n) 2 0 其他 当 n=1 时， (t) ? (2t 1 1 t 1 n) 2 0 其他 故 hn = 1/ 2 nn 2 (t )（2t-n)dt 0 其他 3）、求 gn 。 1/ 2 n 0 gn ( 1)n h1 n 1/ 2 n 1 0 其他 4）、求 (t) 。 (t) gn 1,n (t ) = g0 -1,0 (t) g1 1,1(t) = 1 ? 2 (2 t ) 1 2 (2t 1) 2 2 1 0 t 1 2 1 = 1 1 t 2 0 其他 其图形如下： t t 1、 Haar 尺度函数 (t) 1 0 t 1 0 其他 Haar 尺度函数空间： L , (2 j x 2), (2 j x 1), (2 j x), (2 j x 1), (2 j x 1),Lj 为非负的整数， 该空间又 称为 j 级阶梯函数空间 Vj 。则 V0 V1 V2 L Vj 1 V j V j 1 L 随 j 的增加，分辨更为精细。 2、 性质 函数集 2 j /2 (2 j x k) : k Z 是 Vj 的一个标准正交基。 f ( x) V0 当且仅当 f (2 j x) Vj 。 3、 Haar 小波函数 函数满足两点： （ 1） 是 V1 的成员；（ 2） 与 V0 正交。 ( x) (2 x) (2 x 1) 性质： (x)dx 0 ( x) 是对称的、局部支撑的函数； 小波函数空间 Wj : ak (2 j x k ), ak R k Z W j 是 V j 的正交互补，即 Vj 1 V j Wj 函数集 2 j /2 (2 j x k) : k Z 是 Wj 的一个标准正交基 4、 Haar 小波分解与重建 对 Haar 小波， 有 (2 j x) ( (2 j 1 x) (2 j 1 x)) /2(2 j x 1) ( (2 j 1 x)(2 j 1 x)) / 2 Haar 小波分解定理： 设： f j ( ) ak (2 j x k ), f j ( x ) V j x k Z 则它可以有如下分解： f j f j 1 j 1 f j 1 akj 1 (2 j 1 x k ), f j 1 (x) V j 1 k Z j 1 bkj 1 (2 j 1 x k ), j 1 (x) W j 1 k Z akj 1 a2k a2 k 1 2 bkj 1 a2 k a2 k 1 把函数 f 2 分解成一个小波空间与一个尺度空间的分量 f ( x) 2 (4 x) 2 (4 x 1) (4 x 2) (4 x 3) 解：按照分解定理，此 j=2, ；k=0,1,2,3 对应的系数是 2,2,1，-1；代入公式，得出分解后尺度 函数空间元素的系数是 2 2 2 1 1 2 2 2 ， 0 ；分解后小波函数空间元素的系数是 0 ， 1 ( 1) 2 2 1；从而 2 f ( x) 2 (2 x) (2 x 1) 2 2 f(x) 1 1/2 1/2 -1 1 1 1/2 Matlab 程序 image1=imread(512.jpg); image1=rgb2gray(image1); subplot(2,2,1); imshow(image1); title(original image); image1=double(image1); imagew=imread(shuiyin.bmp); imagew=