哈工大小波实验报告

1、小波理论实验报告院 （系）专 业 学 生 学 号 日 期 2015年12月实验报告一一、 实验目的1. 运用傅立叶变换知识对常用的基本函数做基本变换。2. 加深对因果滤波器的理解，并会判断因果滤波器的类型。3. 运用卷积公式对基本信号做滤波处理并分析，以加深理解。4. 熟悉Matlab中相关函数的用法。二、 实验原理1.运用傅立叶正、反变换的基本公式：及其性质，对所要处理信号做相应的傅里叶变换和逆变换。2.运用卷积的定义式：对所求信号做滤波处理。三、 实验步骤与内容1.实验题目：Butterworth滤波器，其冲击响应函数为 1. 求2. 判断是否因果；是低通、高通、带通还是带阻？3. 对于信

2、号,画出图形4. 画出滤波后图形，比较滤波前后图形，你会发现什么，这里取5. 取采用不同的变量值(初始设定A=10) 画出原信号图形与滤波后图形，比较滤波效果。2.实验步骤及分析过程：1.求由傅里叶变换的定义式可得： （1）故该滤波器的幅频特性为：，转折频率；假定，绘制该滤波器的幅频特性曲线如下：图1.1滤波器的幅频特性曲线2. 判断是否因果；是低通、高通、带通还是带阻？（1）观察滤波器响应函数可知，只有在输入信号到达后，该滤波器才会有输出响应，此外实际应用的滤波器均是因果滤波器，所以，题中滤波器是因果滤波器。（2）由图1可知，该滤波器为低通滤波器。3. 对于信号，画出图形编写matlab程序

3、：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.2 f(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=

4、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,

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(/f(t);grid on可以得到如下图形：图1.4 取不同A和a值后的f(t)比较以上图形中的曲线：可以看出随着A、a值逐渐增大，波形幅值增大，滤波后信号毛刺（高频波动信号）也随着增多，即对高频信号的抑制效果变差，同时也可以看出滤波器输出信号中的低频成分也呈增大趋势。由此可知，滤波器在A、a值较小

6、时对高频的抑制效果最好，但这种情况下低频信号也受到一定的削弱，滤波效果并不一定是最好，因此，需要根据实际使用需求设定参数。实验报告二一、 实验目的1 学习Haar小波的定义及性质，掌握Haar小波分解与重构的原理。2 通过例子学习小波分析在一维信号奇异性检测中的应用；3 学习并掌握信号处理的相关步骤。4 熟悉Matlab中相关函数的用法二、实验原理一般来说，噪声信号多包含在具有较高频率细节中，在对信号进行了小波分解之后，再利用门限阈值等形式对所分解的小波系数进行权重处理，然后对小信号再进行重构即可达到信号去噪的目的。具体步骤为：a.一维信号的小波分解，选择一个小波并确定分解的层次，然后进行分解

7、计算。b.小波分解高频系数的阈值量化，对各个分解尺度下的高频系数选择一个阈值进行软阈值量化处理。C一维小波重构，根据小波分解的最底层低频系数和各层高频系数进行一维小波的重构。利用小波分析检测信号突变点的一般方法是：对信号进行多尺度分析，在信号出现突变时，其小波变换后的系数具有模极大值，因而可以通过对模极大值点的检测来确定故障发生的时间点。通常情况下，信号的奇异性分两种情况，一种是信号在某一个时刻内其幅值发生突变，引起信号的非连续，幅值的突变处是第一种类型的间断点。另一种是信号外观上很光滑，幅值没有突变，但是，信号的一阶微分有突变发生，且一阶微分是不连续的，成为第二种类型的间断点。三、实验步骤与

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

9、，程序如下：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=0&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);grid onv1=wrcoef(a,C,L,Haar,7);subplot(2,4,2);plot(x,v1);title(V1);grid onv2=wrcoef(a,C,L,Haar,6);subplot(2,4,3);pl

10、ot(x,v2);title(V2);grid onv3=wrcoef(a,C,L,Haar,5);subplot(2,4,4);plot(x,v3);title(V3);grid onv4=wrcoef(a,C,L,Haar,4);subplot(2,4,5);plot(x,v4);title(V4);grid onv5=wrcoef(a,C,L,Haar,3);subplot(2,4,6);plot(x,v5);title(V5);grid onv6=wrcoef(a,C,L,Haar,2);subplot(2,4,7);plot(x,v6);title(V6);grid onv7=wrc

11、oef(a,C,L,Haar,1);subplot(2,4,8);plot(x,v7);title(V7);grid on画出的波形如下：图2.3 v中的各个分量将图2.3跟图2.1作比较可知，随着j的增大，得到的图形越接近原函数波形。（3）进行压缩比为80%的压缩，画出压缩后的图像与原图像比较；MATLAB程序如下：F= waverec(C,L,Haar);subplot(1,2,1);plot(x,F);title(原函数);grid onXC,CXC,LXC,PERF0,PERFL2 = wdencmp(gbl,f,Haar,8,0.8,h,1) ;subplot(1,2,2);plot

12、(x,XC);title(压缩比为80%的图形);grid on得到的波形如下：图2.4 压缩后的图像与原图像比较（4）选择合适的参数去噪，画出去噪后的图像与原图像比较；为了去除噪声，我们选择去除高频分量W7，MATLAB程序如下：F= waverec(C,L,Haar);subplot(1,3,2);plot(x,F);title(原信号);grid onw7=wrcoef(d,C,L,Haar,1);subplot(1,3,1);plot(x,w7);title(W7);grid onFl=f-w7;subplot(1,3,3);plot(x,Fl);title(滤波后的信号);grid

13、on得到的波形如下：图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);for

14、i = 1 : 6decmp = wrcoef(d,c,l,db5,7-i);subplot(4,2,i+2);plot(decmp);ylabel(d,num2str(7-i);end得到的图形如下：图2.6 用db5小波分解6层从图中可以看出，对于信号的小波分解的第一层高频系数d1和第二层高频系数d2，可观察到信号的不连续点，在d1图中估测，为x=114/128=0.891，在d2图中估测，为x=112/128=0.875，而实际值为x=8/9=0.889，用db1小波要比db2小波好。实验报告三一、 实验目的1. 掌握双线性插值方法的基本思想，通过实验了解其优缺点。2. 掌握二维多分辨分

15、析的知识和思想，通过对数字图片进行分解与重构操作，加深对理论知识的理解。3. 掌握图片处理的基本方法，并熟悉Matlab中相关函数的应用。二、 实验原理根据插值处理和小波变换的特点，运用一种基于小波分解和双线性插值相结合的图像超分辨率处理方法。首先将原图像进行小波分解，并把原图像作为低通部分，然后对小波分解后的相应高频子带进行双线性插值以近似高频的更多细节，通过小波逆变换获取比原图像分辨率更高的图像。设一幅图像f (x,y)经过一次小波分解后，被分成了四个部分，如图1所示。MH 1为水平方向上的高频细节信息，MV1为垂直方向上的高频细节信息，MD1为对角线方向上的高频细节信息。也就是说，小波分

16、解的过程就是将信号不断“剥落”的过程，随着逼近越来越粗，丢掉的信息越来越多，而被抛弃掉的信息可用小波的线性组合来表示。重建的过程就是将丢掉的细节加起来作为原始信号的近似表示，只要采用足够多的相同步骤，这种近似表示就可以达到足够精确。图3.1三、 实验步骤与内容 实验题目；要求：任意选择一张数字图片，利用二次线性插值将图片放大4倍，同时采用二维小波对原始图片进行分解与重构，实现对图片的超分辨率处理。 MATLAB程序如下：I1=imread(C:UsersAdministratorDesktopHIT.jpg); figure(1); imshow(I1),title(原始图像); Y1,map

17、1=imresize(I1,2,bilinear);c,s=wavedec2(Y1,1,db3); Xa11=appcoef2(c,s,db3,1);Xh11=detcoef2(h,c,s,1); Xv11=detcoef2(v,c,s,1); Xd11=detcoef2(d,c,s,1);Y=idwt2(Xa11,Xh11,Xv11,Xd11,db3); figure(2); imshow(uint8(Y);title(获得超分辨率图像);得到的图形如下： 图3.2 原始图像 图3.3 获得的超分辨率图像利用haar小波分解，MATLAB程序如下：clear all; close all; I1=imread(C:UsersAdministratorDesktopHIT.jpg); figure(1); imshow(I1),title(原始图像); %用双线性插值方法获得插值图像Y1 Y1,map1=imresize(I1,2,bilinear);c,s=wavedec2(Y1,1,haar); Xa11=appcoef2(c,s,haar,1);Xh11=detcoef2(h,c,s