哈工大小波实验.docx

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