小波上机实验报告-哈工大

《小波上机实验报告-哈工大》由会员分享，可在线阅读，更多相关《小波上机实验报告-哈工大（21页珍藏版）》请在金锄头文库上搜索。

1、Harbin Institute of Technology小波分析上机实验报告小波分析上机实验报告院院 系：系： 学学 科：科： 姓姓 名：名： 学学 号：号： 指导教师：指导教师：实验报告一实验报告一一、一、实验目的实验目的1、运用傅里叶变换知识对常用的基本函数做基本变换。2、加深对因果滤波器的理解，并会判断因果滤波器的类型。3、运用卷积公式对基本信号做滤波处理并作出分析，以加深理解4、熟悉 Matlab 中相关函数的用法二、二、实验原理实验原理1运用傅里叶正、反变换的基本公式： ( )( ) ( )( ),11( )( ),22i xi ti ti ti tff x edxf t edt

2、f t ef tfedf te及其性质，对所要处理信号做相应的傅里叶变换和逆变换。2运用卷积的定义式：1212( )( )( )()f tf tff td 对所求信号做滤波处理。三、三、实验步骤与内容实验步骤与内容实验题目：实验题目：Butterworth 滤波器，其冲击响应函数为,0( )0,0若若tAeth tt1. 求( )h 2. 判断是否因果；是低通、高通、带通还是带阻？3. 对于信号,画出图形3( )(sin22sin40.4sin2 sin40 )，tf tetttt0 t( )f t4. 画出滤波后图形，比较滤波前后图形，你会发现什么，这里取( )fh t10A5. 取采用不同

3、的变量值(初始设定 A=10) ( )(sin5sin3sinsin40 )，tf tettttA画出原信号图形与滤波后图形，比较滤波效果。实验步骤及分析过程：实验步骤及分析过程：1.求( )h 由傅里叶变换的定义式可得：（1） 0ti tti tAhAeedtAeedti故该滤波器的幅频特性为：，转折频率；假定 2221 (/) AAH，绘制该滤波器的幅频特性曲线如下：1,2A图 1.1 滤波器的幅频特性曲线2. 判断是否因果；是低通、高通、带通还是带阻？（1）观察滤波器响应函数可知，只有在输入信号到达后，该滤波器才会有输出响应，此外实际应用的滤波器均是因果滤波器，非因果不可用；所以，题中滤

4、波器是因果滤波器。（2）由图 1 可知，该滤波器为低通滤波器。3. 对于信号，画出图形3( )(sin22sin40.4sin2 sin40 )tf tetttt0t ( )f t编写 matlab 程序（见附录 1） ，绘制信号的图形如下：00.511.522.533.5-1.5-1-0.500.511.522.53信 信 f(t)信 信信 信 t信 信 信 f(t)图 1.2 信号 f(t) 曲线4. 画出滤波后图形，比较滤波前后图形，你会发现什么，这里取。( )fh t10A根据卷积运算的滤波特性，编写 matlab 程序，取时间离散点数 n=10000；可绘制该参数下的滤波信号如下图所

5、示：00.20.40.60.811.21.41.61.82 x 104-4000-200002000400060008000、 、 n、 、 、 、 、 f*h、 、 f*h、 、A=a=10图 1.3 f(t) 滤波后信号曲线比较图 1.2 和图 1.3 中，可以看出：经滤波处理后，信号 f（t）中的高频成分得到了有效的抑制，信号的曲线特征变得平滑，而且持续分布相位并未失真，信号的基本信息得到无损传递。5. 取采用不同的变量值(初始设定 A=10) ( )(sin5sin3sinsin40 )tf tettttA画出原信号图形与滤波后图形，比较滤波效果。（1）根据题意，绘制信号曲线；并取 A

6、=1、3、5、1、15 、20(初始设定 A=10)采用给定滤波器，进行滤波处理，结果如下：00.20.40.60.811.21.41.61.82 x 104-50005001000150020002500信 信 f*h信 信信 信 n信 信 信 信 信 f*hA=a=100.20.40.60.811.21.41.61.82 x 104-1000010002000300040005000信 信 f*h信 信信 信 n信 信 信 信 信 f*hA=a=500.20.40.60.811.21.41.61.82 x 104-2000-10000100020003000400050006000信 信

7、f*h信 信信 信 n信 信 信 信 信 f*hA=a=1000.20.40.60.811.21.41.61.82 x 104-2000-10000100020003000400050006000信 信 f*h信 信信 信 n信 信 信 信 信 f*hA=a=1500.20.40.60.811.21.41.61.82 x 104-2000-100001000200030004000500060007000信 信 f*h信 信信 信 n信 信 信 信 信 f*hA=a=20由 1 问可知：该滤波器的幅频特性为：转折频 2221 (/) AAH率；所以当 A=分别取 1、3、5、10、15、20

8、时，滤波器对于该频率成分的信号起到抑制作用。通过观察不同滤波曲线可知：可以看出随着 A、a 值逐渐增大，滤波后信号毛刺 （高频波动信号）增多，即对高频信号的抑制效果变差，同时也可以看出滤波器输 出信号中的低频成分也呈增大趋势。 总之，滤波器在 A、a 值为 1 时对高频的抑制效果最好，但必须指出在这种情况 下低频信号也受到一定的削弱，效果并不一定最好，需要根据实际使用需求设定参 数。程序附录程序附录%待求解信号程序段t=linspace(0,pi,30000);f=exp(-t/3).*(sin(2*t)+2*sin(4*t)+0.4*sin(2*t).*sin(4*t);plot(t,f);

9、grid onxlabel(时间/t);ylabel(信号值/f(t);附录 2%信号图绘制t=linspace(0,pi,30000);f=exp(-t/3).*(sin(2*t)+2*sin(4*t)+0.4*sin(2*t).*sin(4*t);plot(t,f);grid onxlabel(时间/t);ylabel(信号值/f(t);figure%滤波处理后信号A=10;a=10;h=A*exp(-a*t);fy=conv(f,h);plot(fy);grid onxlabel(时间/t);ylabel(滤波信号值/f(t);附录 3%信号图绘制t=linspace(0,pi,3000

10、0);f1=exp(-t).*(sin(5*t)+sin(3*t)+sin(t)+sin(40*t);plot(t,f1);grid onxlabel(时间/t);ylabel(信号值/f(t);figure%滤波处理后信号A=10;a=10;h=A*exp(-a*t);fy=conv(f1,h);plot(fy);grid onxlabel(时间/t);ylabel(滤波信号值/f(t);t=linspace(0,pi,10000); f=exp(-t/3).*(sin(2*t)+2*sin(4*t)+0.4*sin(2*t).*sin(40*t); A=10;a=10; h=A*exp(-

11、a*t); fy=conv(f,h); plot(fy); grid onxlabel(点数n);ylabel(滤波信号值f*h);legend(A=a=10);title(信号f*h波形);实验报告二实验报告二一、一、实验目的实验目的1 学习 Haar 小波的定义及性质，掌握 Haar 小波分解与重构的原理。2 通过例子学习小波分析在一维信号奇异性检测中的应用；3 学习并掌握信号处理的相关步骤。4 熟悉 Matlab 中相关函数的用法二、实验原理二、实验原理在理解了离散小波变换的基本原理和算法的基础上，通过设计 MATLAB 程序对简单的一维信号进行 Haar 小波变换，从而得到小波分解系数

12、；再通过改变分解得到的各层高频系数进行信号的小波重构达到消噪的目的。在这一程序实现的过程中能直观地理解信号小波分解重构的过程和在信号消噪中的重要作用，以及在对各层高频系数进行权重处理时系数的选取对信号消噪效果的影响。一般来说，噪声信号多包含在具有较高频率细节中，在对信号进行了小波分解之后，再利用门限阈值等形式对所分解的小波系数进行权重处理，然后对小信号再进行重构即可达到信号去噪的目的。具体步骤为： a.一维信号的小波分解，选择一个小波并确定分解的层次，然后进行分解计算。 b.小波分解高频系数的阈值量化，对各个分解尺度下的高频系数选择一个阈值进行软阈值量化处理。 C一维小波重构，根据小波分解的最

13、底层低频系数和各层高频系数进行一维小波的重构。 利用小波分析检测信号突变点的一般方法是：对信号进行多尺度分析，在信号出现突变时，其小波变换后的系数具有模极大值，因而可以通过对模极大值点的检测来确定故障发生的时间点。通常情况下，信号的奇异性分两种情况，一种是信号在某一个时刻内其幅值发生突变，引起信号的非连续，幅值的突变处是第一种类型的间断点。另一种是信号外观上很光滑，幅值没有突变，但是，信号的一阶微分有突变发生，且一阶微分是不连续的，成为第二种类型的间断点。三、实验步骤与内容三、实验步骤与内容实验题目：实验题目：1. 设信号将区间0,1实行 256 等分并得到信 24sin 32cos 50.2

14、sincos 55x fexxxx号在这些节点上的离散值（1）利用 Haar 小波对离散后的信号进行分解；（2）画出中的分量并与原信号进行比较；1,8,7,.,1jVjjf（3）进行压缩比为 80%的压缩，画出压缩后的图像与原图像比较；（4）选择合适的参数去噪，画出去噪后的图像与原图像比较；2.定义区间0,1上存在间断点的信号：21,8/91( )0,08/9 xxg xx去信号 g 在区间0,1上 128 等分节点的值，按照 Haar 小波分解与重构算法实现一次分解，描述第 6 层小波系数的值。确定最大的小波系数以及对应该系数的值。判断间断点的大致位置，输出你的理由。实验步骤及分析过程：实验

15、步骤及分析过程：1.第一题第一题（1）利用 Haar 小波对离散后的信号进行分解；信号原图：（2）画出中的分量并与原信号进行比较1,8,7,.,1jVjjf可以看出，与原信号相比，随着的增大，得到的分量越来越逼近原始信号。jfj（3）进行压缩比为 80%的压缩，画出压缩后的图像与原图像比较；从前后图像对比可以发现，当选择压缩比为 80%时，压缩后信号细节损失非常严重，仅 保留了完全重构信号中低频信号成分。（4）选择合适的参数去噪，画出去噪后的图像与原图像比较； 为了去除噪声，我们选择去除高频分量 W7滤波后的图像与原图对比从上下两幅图像的对比中，我们可以看出滤波后的信号曲线高频毛刺减少，使得信号更 加平滑。2.第二题第二题02040608010012014000.20.4信 db5信 信 信 信 信 信 信 s=a6+d6+d5+d4+d3+d2+d1s020406080100120140-0.0500.05a6020406080100120140-0.0500.05d6020406080100120140-0.0200.02d5020406080100120140-0.100.1d4020406080100120140-0.100.1d302