哈工大小波实验报告

小波理论实验报告院（系）专业学生学号日期2015年12月实验报告一一、实验目的1.运用傅立叶变换知识对常用的基本函数做基本变换。2.加深对因果滤波器的理解，并会判断因果滤波器的类型。3.运用卷积公式对基本信号做滤波处理并分析，以加深理解。4.熟悉Matlab中相关函数的用法。二、实验原理1.运用傅立叶正、反变换的基本公式：ˆ()()()(),11ˆ()(),22ixititititffxedxftedtfteftfedfte及其性质，对所要处理信号做相应的傅里叶变换和逆变换。2.运用卷积的定义式：1212()()()()ftftfftd对所求信号做滤波处理。三、实验步骤与内容1.实验题目：Butterworth滤波器，其冲击响应函数为,0()0,0若若tAethtt1.求()h2.判断是否因果；是低通、高通、带通还是带阻？3.对于信号3()(sin22sin40.4sin2sin40)，tftetttt0t,画出图形()ft4.画出滤波后图形()fht，比较滤波前后图形，你会发现什么，这里取10A5.取()(sin5sin3sinsin40)，tftetttt采用不同的变量值A(初始设定A=α=10)画出原信号图形与滤波后图形，比较滤波效果。2.实验步骤及分析过程：1.求()h由傅里叶变换的定义式可得：0ˆtittitAhAeedtAeedti（1）故该滤波器的幅频特性为：2221(/)AAH，转折频率；假定1,2A，绘制该滤波器的幅频特性曲线如下：图1.1滤波器的幅频特性曲线2.判断是否因果；是低通、高通、带通还是带阻？（1）观察滤波器响应函数可知，只有在输入信号到达后，该滤波器才会有输出响应，此外实际应用的滤波器均是因果滤波器，所以，题中滤波器是因果滤波器。（2）由图1可知，该滤波器为低通滤波器。3.对于信号3()(sin22sin40.4sin2sin40)tftetttt0t，画出图形()ft编写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)');gridon绘制信号的图形如下：图1.2f(t)波形图4.画出滤波后图形()fht，比较滤波前后图形，你会发现什么，这里取10A。对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)');gridon运行程序，得到的图形如下：图1.3滤波后的f(t)比较图1.2和图1.3中，可以看出：经滤波处理后，信号f（t）的幅值变大，高频成分得到了有效的抑制，信号的曲线特征变得平滑，而且持续分布相位并未失真，信号的基本信息得到无损传递。5.取()(sin5sin3sinsin40)tftetttt采用不同的变量值A(初始设定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);holdon;plot(F);xlabel('ʱ¼ä/t');ylabel('Â˲¨ÐźÅÖµ/f(t)');gridon可以得到如下图形：图1.4取不同A和a值后的f(t)比较以上图形中的曲线：可以看出随着A、a值逐渐增大，波形幅值增大，滤波后信号毛刺（高频波动信号）也随着增多，即对高频信号的抑制效果变差，同时也可以看出滤波器输出信号中的低频成分也呈增大趋势。由此可知，滤波器在A、a值较小时对高频的抑制效果最好，但这种情况下低频信号也受到一定的削弱，滤波效果并不一定是最好，因此，需要根据实际使用需求设定参数。实验报告二一、实验目的1．学习Haar小波的定义及性质，掌握Haar小波分解与重构的原理。2．通过例子学习小波分析在一维信号奇异性检测中的应用；3．学习并掌握信号处理的相关步骤。4．熟悉Matlab中相关函数的用法二、实验原理一般来说，噪声信号多包含在具有较高频率细节中，在对信号进行了小波分解之后，再利用门限阈值等形式对所分解的小波系数进行权重处理，然后对小信号再进行重构即可达到信号去噪的目的。具体步骤为：a.一维信号的小波分解，选择一个小波并确定分解的层次，然后进行分解计算。b.小波分解高频系数的阈值量化，对各个分解尺度下的高频系数选择一个阈值进行软阈值量化处理。C．一维小波重构，根据小波分解的最底层低频系数和各层高频系数进行一维小波的重构。利用小波分析检测信号突变点的一般方法是：对信号进行多尺度分析，在信号出现突变时，其小波变换后的系数具有模极大值，因而可以通过对模极大值点的检测来确定故障发生的时间点。通常情况下，信号的奇异性分两种情况，一种是信号在某一个时刻内其幅值发生突变，引起信号的非连续，幅值的突变处是第一种类型的间断点。另一种是信号外观上很光滑，幅值没有突变，但是，信号的一阶微分有突变发生，且一阶微分是不连续的，成为第二种类型的间断点。三、实验步骤与内容实验题目：1.设信号24sin32cos50.2sincos55xfexxxx将区间[0,1]实行256等分并得到信号在这些节点上的离散值（1）利用Haar小波对离散后的信号进行分解；（2）画出1,8,7,...,1jVj中的分量并与原信号jf进行比较；（3）进行压缩比为80%的压缩，画出压缩后的图像与原图像比较；（4）选择合适的参数去噪，画出去噪后的图像与原图像比较；2.定义区间[0,1]上存在间断点的信号：21,8/91()0,08/9xxgxx去信号g在区间[0,1]上128等分节点的值，按照Haar小波分解与重构算法实现一次分解，描述第6层小波系数的值。确定最大的小波系数以及对应该系数的值。判断间断点的大致位置，输出你的理由。实验步骤及分析过程：1.第一题（1）利用Haar小波对离散后的信号进行分解；先画出原信号的波形，程序如下：x=linspace(0,1,256);f=exp(-x.^2/4).*(sin(3*x)+2*cos(5*x)+0.2*sin(x).*cos(55*x)).*(x=0&x=1);plot(x,f)；得到图形如下：图2.1原信号波形下面运用Haar小波对已知函数进行分解，程序如下所示：x=linspace(0,1,256);j7=linspace(0,1,128);f=exp(-x.^2/4).*(sin(3*x)+2*cos(5*x)+0.2*sin(x).*cos(55*x)).*(x=0&x=1);[a,b]=dwt(f,'haar');subplot(2,1,1);plot(j7,a);xlabel('尺度空间');gridonsubplot(2,1,2);plot(j7,b);xlabel('小波空间');gridon得到的波形如下：图2.2haar小波对函数进行分解（2）画出1,8,7,...,1jVj中的分量并与原信号jf进行比较；MATLAB程序如下：x=linspace(0,1,256);f=exp(-x.^2/4).*(sin(3*x)+2*cos(5*x)+0.2*sin(x).*cos(55*x)).*(x=0&x=1);[C,L]=wavedec(f,8,'Haar');v0=wrcoef('a',C,L,'Haar',8);subplot(2,4,1);plot(x,v0);title('V0');gridonv1=wrcoef('a',C,L,'Haar',7);subplot(2,4,2);plot(x,v1);title('V1');gridonv2=wrcoef('a',C,L,'Haar',6);subplot(2,4,3);plot(x,v2);title('V2');gridonv3=wrcoef('a',C,L,'Haar',5);subplot(2,4,4);plot(x,v3);title('V3');gridonv4=wrcoef('a',C,L,'Haar',4);subplot(2,4,5);plot(x,v4);title('V4');gridonv5=wrcoef('a',C,L,'Haar',3);subplot(2,4,6);plot(x,v5);title('V5');gridonv6=wrcoef('a',C,L,'Haar',2);subplot(2,4,7);plot(x,v6);title('V6');gridonv7=wrcoef('a',C,L,'Haar',1);subplot(2,4,8);plot(x,v7);title('V7');gridon画出的波形如下：图2.3v中的各个分量将图2.3跟图2.1作比较可知，随着j的增大，得到的图形越接近原函数波形。（3）进行压缩比为80%的压缩，画出压缩后的图像与原图像比较；MATLAB程序如下：F=waverec(C,L,'Haar');subplot(1,2,1);plot(x,F);title('原函数');gridon[XC,CXC,LXC,PERF0,PERFL2]=wdencmp('gbl',f,'Haar',8,0.8,'h',1);subplot(1,2,2);plot(x,XC);title('压缩比为80%的图形');gridon得到的波形如下：图2.4压缩后的图像与原图像比较（4）选择合适的参数去噪，画出去噪后的图像与原图像比较；为了去除噪声，我们选择去除高频分量W7，MATLAB程序如下：F=waverec(C,L,'Haar');subplot(1,3,2);plot(x,F);title('原信号');gridonw7=wrcoef('d',C,L,'Haar',1);subplot(1,3,1);plot(x,w7);title('W7');gridonFl=f-w7;subplot(1,3,3);plot(x,Fl);title('滤波后的信号');gridon得到的波形如下：图2.5去噪后的图像与原图像比较可以看出，去噪后的图形与原图基本相同。毛刺变少，信号更加平滑。第二题MATLAB程序如下：x=linspace(0,1,128);s=0.*(x8/9)+(1-x.*x).*(8/9=x&x=1);ls=length(s);[c,l]=wavedec(s,6,'db5');subplot(4,2,1);plot(s);title('用db5分解6层：s=a6+d6+d5+d4+d3+d2+d1');ylabel('s');a6=wrcoef('a',c,l,'db5',6);subplot(4,2,2);plot(a6);ylabel('a6');fori=1:6decmp=wrcoef('d',c,l,'db5',7-i);subplot(4,2,i+2);plot(decmp);ylabel(['d