哈工大小波实验报告.doc

小波理论实验报告院(系)专业学生学号日期2015年12月实验报告一实验目的运用傅立叶变换知识对常用的基本函数做基本变换。加深对因果滤波器的理解,并会判断因果滤波器的类型。运用卷积公式对基本信号做滤波处理并分析,以加深理解。熟悉Matlab中相关函数的用法。、反变换的基本公式:及其性质,对所要处理信号做相应的傅里叶变换和逆变换。:对所求信号做滤波处理。:Butterworth滤波器,;是低通、高通、带通还是带阻?,,比较滤波前后图形,你会发现什么,(初始设定A=α=10)画出原信号图形与滤波后图形,比较滤波效果。::(1)故该滤波器的幅频特性为:,转折频率;假定,绘制该滤波器的幅频特性曲线如下:;是低通、高通、带通还是带阻?(1)观察滤波器响应函数可知,只有在输入信号到达后,该滤波器才会有输出响应,此外实际应用的滤波器均是因果滤波器,所以,题中滤波器是因果滤波器。(2)由图1可知,该滤波器为低通滤波器。,画出图形编写matlab程序:t=linspace(0,pi,80000);f=exp(-t/3).*(sin(2*t)+2*sin(4*t)+*sin(2*t).*sin(40*t));plot(t,f);xlabel('时间/t');ylabel('信号值/f(t)');gridon绘制信号的图形如下:(t),比较滤波前后图形,你会发现什么,这里取。对f(t)进行卷积运算,编写MATLAB程序,如下:A=10;a=10;t=linspace(0,pi,80000);f=exp(-t/3).*(sin(2*t)+2*sin(4*t)+*sin(2*t).*sin(40*t));h=A*exp(-a*t);F=conv(f,h);plot(F);xlabel('时间/t');ylabel('滤波后信号值/f(t)');gridon运行程序,得到的图形如下:(t),可以看出:经滤波处理后,信号f(t)的幅值变大,高频成分得到了有效的抑制,信号的曲线特征变得平滑,而且持续分布相位并未失真,信号的基本信息得到无损传递。(初始设定A==10)画出原信号图形与滤波后图形,比较滤波效果。 对A和a分别取A=a=2,5,10,15,25,并将几个图形放在一起比较,MATLAB程序如下:A=25;a=25;%A==a,a=2,5,10,15,25%t=linspace(0,pi,80000);f=exp(-t/3).*(sin(2*t)+2*sin(4*t)+*sin(2*t).*sin(40*t));h=A*exp(-a*t);F=conv(f,h);holdon;plot(F);xlabel('ʱ¼ä/t');ylabel('Â˲¨ÐźÅÖµ/f(t)');gridon可以得到如下图形:(t)比较以上图形中的曲线:可以看出随着A、a值逐渐增大,波形幅值增大,滤波后信号毛刺(高频波动信号)也随着增多,即对高频信号的抑制效果变差,同时也可以看出滤波器输出信号中的低频成分也呈增大趋势。由此可知,滤波器在A、a值较小时对高频的抑制效果最好,但这种情况下低频信号也受到一定的削弱,滤波效果并不一定是最好,因此,需要根据实际使用需求设定参数。,掌握Haar小波分解与重构的原理。;。、实验原理一般来说,噪声信号多包含在具有较高频率细节中,在对信号进行了小波分解之后,再利用门限阈值等形式对所分解的小波系数进行权重处理,然后对小信号再进行重构即可达到信号去噪的目的。具体步骤为: ,选择一个小波并确定分解的层次,然后进行分解计算。 ,对各个分解尺度下的高频系数选择一个阈值进行软阈值量化处理。 ,根据小波分解的最底层低频系数和各层高频系数进行一维小波的重构。 利用小波分析检测信号突变点的一般方法是:对信号进行多尺度分析,在信号出现突变时,其小波变换后的系数具有模极大值,因而可以通过对模极大值点的检测来确定故障发生的时间点。通常情况下,信号的奇异性分两种情况,一种是信号在某一个时刻内其幅值发生突变,引起信号