小波分析实验：二维离散小波变换(Mallat快速算法).docxVIP

小波分析实验： 实验 小波分析实验： 实验2二维离散小波变换（Mallat快速算法） 小波分析实验： 实验 小波分析实验： 实验2二维离散小波变换（Mallat快速算法） 实验目的： 在理解离散小波变换原理和 Mallat快速算法的基础上，通 过编程对图像进行二维离散小波变换，从而加深对二维小波分解 和重构的理性和感性认识，并能提高编程能力，为今后的学习和 工作奠定基础。 实验工具： 计算机，matlab6.5 分解算法: 重构算法: “二工必（刃- 2上*［十三g（刃- 2k）d［ * |g[0] g[l] g[2] g[3] 0 0 0 I 0 0 g[0] g[l]g[2] S[3] - 0 ? ????? ? ? ?????■ ? ? g[2] g[3] 0 0 0 0 ...g[0] g[l]J |_勺4-1[ ?(2?分解算法写成矩阵的形式！ （lb g的长度为4） ?(2? ?勺[0] _ ^[0] A[l] h[2] h[3] 0 0 0 1 r订]] 勺?卩] ■ ■ ■ 二 0 0 h[0] h[}]h[2] h[3] … 0 ? ? ? ?? ? ? ? 申⑵ 拭3] 0 0 0 0 …凤0] /?[!]] 叨］ I 二 4[0]如]力⑵ h[3] 0 0 0 ［勺【0】? 记 0 0 h[0] h[\]h[2]山⑶ … 0 ? ????? ? ? Cj = 勺【1] ? ? ? 申[2] h[3] 0 0 0 0 -.^[0] ^[1]_ .勺［乃-1】_ [0] g[l] g⑵ g[3] 0 ? ? ? 0 0 e= 0 ? 0 ? g[0] g[l]g ⑵ ? ? g[3] ■ ? ?? ■ 0 ? DJ = [i] ■ ? ■ 目2] ■ g[3] 0 0 0 0 …茎0] 畀］ 于是Mallat分解公式为矩阵变换？丄 TOC \o 1-5 \h \z Cj- = PC^ ⑶卩 \o Current Document Dj = Qdj-l ⑷ 重构算法写成矩阵变换：- Cj_i = C$ + Dj (5) 4 M N \o Current Document PPq. 一片『 峰值信噪比计算公式：P沁沁逻竺皿e卢H耿V 屈E M ｛皿，00分别表示原始图像和重建图像，且 本实验采取的一些小技乐P (I) 分SW法… 加0]方[1]农］ 加0] 方[1] 农］ 舛3] 0 0 ? ? 0 o T = 0 纽0] 列1] 風2］疏3］ 0 ■ 町] ■ 列2] ■ 扯3] ? ? D 0 ? 0 …0 sm g[l] 列2] 或习0 ? ? ? 0 0 11 = 0 烈0] 畀］ 02]烈 3] 0 :.° ■ -g[l] ■ g[2] ■ g[3] ? ? 0 0 ■ 0 ???0 g[0]J[ %[1] 勺4-1] 怕［0】 编程时用如下思想：（h, g的长度为4）“ 今［1］ 勺［刀-1］ ■ V ■ ■丐⑼ ￡［1］ 4刀-1】 将数据。』和低通〔高通）滤波器进行添零到数据长度诳，再2抽祥，相加得到 （2） 科瞎 小波 系数丄 后一半 系数门 高通. 重构 滤波 器卩 重构结果 编程吋用如下思想：2 ii[0] 山1] 爪)] 川］ /〔0] 0 00 0 00 ? ? ? ? ? ? ? ? 0 00 ?托3] 42] 0??托3] ? ? ? ? ? ? 0 ??/(!] 0 0 00…灶3] 夙0] 0 00 0 ??鼻 ? ? ? ? ? 0 0 00 0 ■心[0]「 心[1] ■ ? ■ 十 咕[0「 吆[1] ? ? ? 弓4-1]_ /j-Jn-l] 5?1[0] 5 di] 川] ■ ??炷1] 址0] 0 勺[1] 0 曲］ 「77+1 实验结果：“ 多尺度分堺E3像 SC1）三级小波分解图Q 50 100 150 200 250 原始图條 逅构图像 原始图條 逅构图像 50 15C 200 250 50 100 130 200 250 EC2 ）原始凰慷和重构图像 峰值信噪比（psm—252.8923db^ 附录： 二维小波分解函数 %二维小波分解函数 fun cti on Y=mallatdec2(X,w name,level) %输入:X 载入的二维图像像数值； % level 小波分解次(级)数设定值(如果设定值超过最高可分解次数 ,按最高分解次 数分解) % wname 小波名字 wavelet name %输岀：丫 多极小波分解后的小波系数矩阵 [h,g]=wfilters(wname,d); %h,g分别为低通和高通滤波器 X=double(X); t=1; hh=size(X,2); while t=level %先进行行小波变换 fo