CourseWavelet同济大学宣国荣老师小波分析的讲座PPS格式.pps
《CourseWavelet同济大学宣国荣老师小波分析的讲座PPS格式.pps》由会员分享,可在线阅读,更多相关《CourseWavelet同济大学宣国荣老师小波分析的讲座PPS格式.pps(60页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、小波分析及其应用小波分析及其应用 Wavelet Analysis andIts Applications同济大学 计算机系宣国荣2003年6月10日 星期二1研究生讲座:小波分析及其应用研究生讲座:小波分析及其应用1、小波的特点和发展2、小波分析在一维信号处理中的应用3、小波分析在图象分析中的应用 图象特征抽取 图象压缩 数据隐藏和图象水印21、小波的特点和发展小波的特点和发展 “小波分析”是分析原始信号各种变化的特性,进一步用于数据压缩、噪声去除、特征选择等。例如歌唱信号:是高音还是低音,发声时间长短、起伏、旋律等。从平稳的波形发现突变的尖峰。小波分析是利用多种“小波基函数”对“原始信号”
2、进行分解。3小波的时间和频率特性小波的时间和频率特性运用小波基,可以提取信号中的“指定时间”和“指定频率”的变化。时间:提取信号中“指定时间”(时间A或时间B)的变化。顾名思义,小波在某时间发生的小的波动。频率:提取信号中时间A的比较慢速变化,称较低频率成分;而提取信号中时间B的比较快速变化,称较高频率成分。时间A时间B4小波的小波的成就成就小波分析是纯数学、应用数学和工程技术的完美结合。从数学来说是大半个世纪“调和分析”的结晶(包括傅里叶分析、函数空间等)。小波变换是20世纪最辉煌科学成就之一。在计算机应用、信号处理、图象分析、非线性科学、地球科学和应用技术等已有重大突破,预示着小波分析进一
3、步热潮的到来。5多分辨度分析(多分辨度分析(MRA)1988年Mallat提出的多分辨度分析理论,统一了几个不相关的领域:包括语音识别中的镜向滤波,图象处理中的金字塔方法,地震分析中短时波形处理等。当在某一个分辨度检测不到的现象,在另一个分辨度却很容易观察处理。例如:6 7 参考:M.Vetterli,”WaveletsandSubbandCoding“,PrenticeHallPTR,1995p.118小波的小波的3 个特点个特点小波变换,既具有频率分析的性质,又能表示发生的时间。有利于分析确定时间发生的现象。(傅里叶变换只具有频率分析的性质)小波变换的多分辨度的变换,有利于各分辨度不同特征
4、的提取(图象压缩,边缘抽取,噪声过滤等)小波变换比快速Fourier变换还要快一个数量级。信号长度为M时,Fourier变换(左)和小波变换(右)计算复杂性分别如下公式:9小波基表示发生的时间和频率小波基表示发生的时间和频率“时频局域性”图解:Fourier变换的基(上)小波变换基(中)和时间采样基(下)的比较傅里叶变换(Fourier)基小波基时间采样基10 Haar小小波基母函数波基母函数(a)Haar“近似”基函数(b)Haar“细节”基函数低频滤波系数高频滤波系数H0=11qH1=1-1q=qq=q-q其中:11Haar小波的基函数小波的基函数第1行基函数是取平均(近似),第2-8行基
5、函数是取变化(细节)。细节包括变化速率和发生的时间。H0=11qH1=1-1q尺度函数近似基函数小波函数细节基函数12小波分析发展历史小波分析发展历史1807年Fourier提出傅里叶分析,1822年发表“热传导解析理论”论文1910年Haar提出最简单的小波1980年Morlet首先提出平移伸缩的小波公式,用于地质勘探。1985年Meyer和稍后的Daubeichies提出“正交小波基”,此后形成小波研究的高潮。1988年Mallat提出的多分辨度分析理论(MRA),统一了语音识别中的镜向滤波,子带编码,图象处理中的金字塔法等几个不相关的领域。13小波基可以通过给定滤波系数生成小波基可以通过
6、给定滤波系数生成小波基(尺度函数和小波函数)可以通过给定滤波系数生成。有的小波基是正交的,有的是非正交的。有的小波基是对称的,有的是非对称的。小波的近似系数和细节系数可以通过滤波系数直接导出,而不需要确切知道小波基函数,这是I.Daubechies等的重要发现,使计算简化,是快速小波分解和重建的基础。14小波基函数和滤波系数小波基函数和滤波系数(Haar-正交,对称正交,对称)“近似”基函数“反变换”低频和高频“滤波系数”“细节”基函数Haar小波“正变换”低频和高频“滤波系数”15小波基函数和滤波系数小波基函数和滤波系数(db 2-正交,不对称正交,不对称)“近似”基函数“细节”基函数db小
7、波“反变换”低频和高频“滤波系数”“正变换”低频和高频“滤波系数”16小波基函数和滤波系数小波基函数和滤波系数(db 4-正交,不对称正交,不对称)17小波基函数和滤波系数小波基函数和滤波系数(sym 4-正交,近似对称正交,近似对称)18小波基函数和滤波系数小波基函数和滤波系数(bior 2.4 双正交,对称双正交,对称)19小波基函数和滤波系数小波基函数和滤波系数(bior 6.8 双正交,对称双正交,对称)202 2、小波、小波分析分析在一维信号处理中的应用在一维信号处理中的应用小波变换小波变换就是将“原始信号s”变换成“小波系数w”,w=wa,wd包括近似(approximation)
8、系数wa与细节(detail)系数wd近似系数wa-平均成分(低频)细节系数wd-变化成分(高频)21小波原始信号分解过程:小波原始信号分解过程:原始信号s可分解成小波近似a与小波细节d之和。s=a+d小波系数w=wa,wd的分量,乘以基函数,形成小波分解:小波近似系数wa基函数A=近似分解a-平均小波细节系数wd基函数D=细节分解d-变化22小波分解和小波分解和小波基小波基小波基D小波基A原始信号小波系数wd小波系数wa正变换:原始信号在小波基上,获得“小波系数”分量反变换:所有“小波分解”合成原始信号例如:小波分解a=小波系数wa小波基A23离散小波变换公式离散小波变换公式正变换反变换其中
9、:是小波基函数参考“数字图象处理”英文版,电子工业出版社,2002年(R.C.Gonzalaz,”DigitalImageProcessing”,p.375)信号s有M个样本,J级小波变换:小波分解小波系数24一维信号小波变换例子一维信号小波变换例子Haar小波,例子:16点信号:6 5 9 8 3 7 8 5 6 5 9 8 3 7 8 5 6 5 9 8 1 3 3 9 6 5 9 8 1 3 3 9通过MATLAB实现(wavemenu)波形图小波正变换:小波系数:小波近似系数(加);小波细节系数(减)小波反变换:可以由分解信号恢复原始信号。有2种:近似分解;细节分解25一维信号的二级小
10、波变换系数一维信号的二级小波变换系数原始信号2级小波系数 w2=wa2,wd2,wd1*Haar是正交变换。除以常数,目的使变换后平方和不变。例如:16位2级近似系数2级细节系数1级细节系数16位26一维信号的二级小波变换分解一维信号的二级小波变换分解2级近似分解 (原始信号每4个平均值)2级细节分解 (原始信号每2个平均的差值)1级细节分解 (原始信号单数和双数的差值)恢复信号27一维信号的二级小波变换系数和分解一维信号的二级小波变换系数和分解原始信号2级小波系数w2=wa2,wd2,wd12级近似分解 (原始信号每4个平均值)2级细节分解 (原始信号每2个平均的差值)1级细节分解 (原始信
11、号单数和双数的差值)恢复信号28原始信号原始信号 16点点 16点原始信号点原始信号 6 5 9 8 3 7 8 5 6 5 9 8 1 3 3 9 29两级小波系数两级小波系数16点点原始信号小波系数原始信号(红)两级小波系数wd1wd2|wd2|wd1|3016点点 信号信号 的的Haar小波近似值和细节分解小波近似值和细节分解 两级分解31小波小波去噪声去噪声一般噪声特点:一般噪声特点:(1)高频成分(细节),(2)幅度小:用阈值;去噪声过程:去噪声过程:去除原始信号高频成分(细节)中幅度小于阈值部分。对2级小波,设定2个阈值,称“阈值2”和“阈值1”。去除1级噪声:去除1级小波细节分解
12、中小于“阈值1”部分。去除2级噪声:去除2级小波细节分解中小于“阈值2”部分。恢复:恢复:将小波近似分解,加上去噪声后小波细节分解,即获得去除噪声的信号32噪声去除噪声去除两级分解噪声去除,括号内保留部分数据原始信号(红),去噪后(黄)wd1两级小波系数wd233小波小波去噪声去噪声16点点 6 5 9 8 3 7 8 5 6 5 9 8 1 3 3 9|wd1|1级去噪前绝对值|wd1|1级去噪后绝对值|wd2|2级去噪后绝对值|wd2|2级去噪前绝对值原始信号(红),去噪后(黄)1级细节小波系数2级细节小波系数0.7071,1,-4,3,1,1,-2,-60.5-6,-3,-6,-8两级小
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- CourseWavelet 同济大学 宣国荣 老师 分析 讲座 PPS 格式
限制150内