小波变换双正交小波讲稿.ppt

《小波变换双正交小波讲稿.ppt》由会员分享，可在线阅读，更多相关《小波变换双正交小波讲稿.ppt（45页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、小波变换课件双正交小波小波变换课件双正交小波第一页，讲稿共四十五页哦正交小波的性质对称性（），紧支撑（）对称性（），紧支撑（）对称性（），紧支撑（）光滑性（）Harr小波紧支撑且线性相位(对称性)？双正交小波！第二页，讲稿共四十五页哦5.1滤波器的相位特性在线性系统理论中,滤波器的传递函数可表达为 为幅频特性，为相频特性。如果 可以表示为 式中和为常数，那么称为 具有线性相位特性 第三页，讲稿共四十五页哦输出信号的相位特性，除了常数外，与延时为的输入信号 的相位特性完全一致hf(x)g(x)第四页，讲稿共四十五页哦当滤波器具有线性相位特性时,输出信号将不产生相位畸变。这一点对图像信号十分重要，

2、因为视觉对于相位畸变非常敏感 线性相位，振幅畸变非线性相位，振幅无畸变滤波器如何具有线性相位特性？第五页，讲稿共四十五页哦 定理5.1 滤波器 具有线性相位的必要与充分条件是它的脉冲响应函数具有如下关于 的共轭对称性：证明：必要性：充分性：第六页，讲稿共四十五页哦 推论5.1 如果限定脉冲响应 为实函数，那么由式(5.1.3)可知，这时 必为实数,即 所以实脉冲响函数具有线性相位的必要与充分条件是 任何实值脉冲响应的数字滤波器具有线性相位的必充条件是 第七页，讲稿共四十五页哦正交小波：紧支撑线性相位？定理5.2 紧支撑正交小波，除Haar小波之外，不可能是线性相位的。第八页，讲稿共四十五页哦5

3、.2 双正交小波的基本性质如果有两对函数 与 ，其中，尺度函数 和 分别生成MRA 和MRA ，而 和 则分别张成在下述意义上的补空间 和 ：并且它们之间还满足如下正交关系：那么这两对函数称为互为对偶的双正交小波。第九页，讲稿共四十五页哦5.2.2 双正交小波的二尺度关系二尺度关系第十页，讲稿共四十五页哦令则双正交基本条件若令则正交基本条件第十一页，讲稿共四十五页哦5.2.3紧支撑线性相位双正交小波的 和 之间的长度关系定理5.3 序列 ，除 之外，所有下标为偶数的元素取值为0。证明：利用双正交基本条件第十二页，讲稿共四十五页哦 和 的长度 和 之间的关系（1）两者长度均为奇数，并且长度相差的

4、奇数倍，因而两者不可能等长。小波 和 是一对偶对称双正交小波。第十三页，讲稿共四十五页哦（2）；两者长度均为偶数，并且长度相差的偶数倍，两者等长是可能的。小波 和 是一对反对称双正交小波。第十四页，讲稿共四十五页哦（3）L为奇数,为偶数或为L偶数,为奇数 序列 的终点下标和起点下标关于奇偶性出现矛盾，故此种情况不存在 第十五页，讲稿共四十五页哦5.3 构造双正交小波的CDF方法Step1 给定2M，根据下式计算Step2Step3 当 和 均为偶数，当 和 均为奇数，Step4 第十六页，讲稿共四十五页哦function y,t=biofilter1(n)t(1)=sym(1);%syms是定

5、义符号变量；sym则是将字符或者数字转换为字符。for i=1:n t(i+1)=t(i)*(n+1-i)/i;endfor i=1:n t=t/2;endy=sym(0);syms z;n2=floor(n/2);%朝负无穷方式舍入for i=-n2:n-n2 y=y+t(n2+i+1)*zi;end第十七页，讲稿共四十五页哦function y=biofilter2(n,m)k=(n+m)/2;t(1)=sym(1);for p=2:k t(p)=sym(1);for j=1:p-1 t(p)=t(p)*(k-1+p-1+1-j)/j;end endm0=sym(0);syms z;for

6、 j=1:k m0=m0+t(j)*(z(-1/2)-z(1/2)/(2*i)(2*(j-1);endif n/2=fix(n/2)%fix(n)的意义是取小于n的整数(是向零点舍入）m0=m0*(z(1/2)+z(-1/2)/2)m;else m0=m0*z(1/2)*(z(1/2)+z(-1/2)/2)m;end y=collect(m0);%返回系数整理后的多项式 第十八页，讲稿共四十五页哦 biofilter2(2,4)ans=1/4/z+1/2+1/4*z ans=3/128*z4-3/64*z3-1/8*z2+19/64*z+45/64+19/64/z-1/8/z2-3/64/z3

7、+3/128/z4 biofilter2(4,2)ans=1/16/z2+1/4/z+3/8+1/4*z+1/16*z2 ans=3/32*z3-3/8*z2+5/32*z+5/4+5/32/z-3/8/z2+3/32/z3第十九页，讲稿共四十五页哦 biofilter2(3,3)ans=1/8/z+3/8+3/8*z+1/8*z2 ans=3/64*z4-9/64*z3-7/64*z2+45/64*z+45/64-7/64/z-9/64/z2+3/64/z3第二十页，讲稿共四十五页哦 biofilter2(1,5)ans=1/2+1/2*z ans=3/256*z5-3/256*z4-11/

8、128*z3+11/128*z2+1/2*z+1/2+11/128/z-11/128/z2-3/256/z3+3/256/z4 biofilter2(5,1)ans=1/32/z2+5/32/z+5/16+5/16*z+5/32*z2+1/32*z3 ans=3/16*z3-15/16*z2+5/4*z+5/4-15/16/z+3/16/z2 第二十一页，讲稿共四十五页哦5.4 基于双正交小波的分解与重构任给重构基 分解基 在双正交小波情况下，信号的分解与重构采用不同的基 第二十二页，讲稿共四十五页哦 (a)分解 (b)重构第二十三页，讲稿共四十五页哦第二十四页，讲稿共四十五页哦第二十五页，讲

9、稿共四十五页哦5.5 小波函数的消失矩性质定义5.1 当小波函数 满足如下条件时 称 具有N阶消失矩 定理 5.4 具有N阶消失矩，则 以 为其N重零点。推论5.2 具有N阶消失矩，则 以 为其N重零点。第二十六页，讲稿共四十五页哦5.5.2 小波函数的光滑性与消失矩的关系数学上用函数 的频谱 在 足够大时的衰减快慢来刻画 的光滑程度 定义5.2 如果存在尽可能大的正常数 ，使 成立，则称 具有光滑指数 ，并记为 越大，则表明 衰减愈快，因而 的 频域定域性愈好 第二十七页，讲稿共四十五页哦5.6 提升方案1994年Wim Sweldens提出了一种新的小波构造方法提升方案(lifting s

10、cheme)，也叫第二代小波变换(second generation wavelet transform,SGWT)或整数到整数小波变换(integer-to-integer wavelet transform,ITIWT)。第二代小波变换构造方法的特点是：1、继承了第一代小波的多分辨率的特性；2、不依赖傅立叶变换,直接在时域完成小波变换；3、小波变换后的系数可以是整数；4、图象的恢复质量与变换时边界采用何种延拓方式无关。第二十八页，讲稿共四十五页哦定理5.6 对于一对给定的双正交尺度函数 和 ，其相应的二尺度系数序列为 和 ，那么所有满足如下条件 的序列 将与序列 也满足双正交关系。式中 是

11、一个三角函数多项式。提升方案的基本公式第二十九页，讲稿共四十五页哦简单的提升步骤：已知：任何一对双正交尺度函数 和 (对应于 和 )，步骤：选择合适的三角多项式 并根据下式获得新的 和新的 用户定制（用户定制（custom_design)第三十页，讲稿共四十五页哦从haar小波出发的提升已知：要求：提升后小波具有二阶消失矩提升：依据定理5.4，提升后的小波所对应的 在 应具有二重零点，由此获得三角多项式 的具体表达式，并代入提升方案的基本公式即可。将Lazy小波提升到具有二阶消失矩的过程与此相同第三十一页，讲稿共四十五页哦交替提升：第一次提升：（和 不变）第二次提升：（和 不变）第三十二页，讲

12、稿共四十五页哦 Sweldens分解和重构算法 1、分解算法推导第三十三页，讲稿共四十五页哦按照（5.6.3）提升 则对应的序列 根据标准分解算法有（5.6.3）（4.1.5）第三十四页，讲稿共四十五页哦一次提升的分解步骤 关于 （提升 和 ）：第一步：用 和 对输入信号按标准算法作一次分解得 第二步：更新 第三十五页，讲稿共四十五页哦一次提升的分解步骤 关于 （提升 和 ）第一步：用 和 对输入信号按标准算法作一次分解得 第二步：更新 第三十六页，讲稿共四十五页哦交替提升的分解步骤第三十七页，讲稿共四十五页哦2、重构算法针对图示的分解算法，推导其重构过程第三十八页，讲稿共四十五页哦按照（5.

13、6.4）提升 则对应的序列 根据标准重构算法有（5.6.4）（4.2.4）第三十九页，讲稿共四十五页哦一次提升的重构步骤 第一步：计算 （去提升）第二步：计算第四十页，讲稿共四十五页哦一次提升的Sweldens重构算法 交替提升的Sweldens重构算法 第四十一页，讲稿共四十五页哦作业利用Sweldens算法实现信号的分解与重构可参考第四章示例程序的信号和模式给出图形结果使用haar小波或lazy小波的一次提升第四十二页，讲稿共四十五页哦整型小波变换提升方案为扩展小波变换的应用领域提供了更多的灵活性。常规的小波变换都是采用浮点运算的，但利用提升方案所带来的便利，可十分方便地构造整数到整数的小波变换。将整数小波变换用于图像压缩就可以用小波变换进行无失真的图像压缩。第四十三页，讲稿共四十五页哦即位算法（In_Place）以Lazy小波为例，进行分解。在同一数组 中，偶数下标的单元存放 ，奇数下标的单元存放 ，即 进行交替提升时，先更新细节序列：如果同样将细节序列乘以 ，并存入数组y.的奇数下标元素中，则 细节序列更新细节序列更新近似序列更新近似序列更新第四十四页，讲稿共四十五页哦重构（同址运算）整型提升(3,5)小波变换 第四十五页，讲稿共四十五页哦