小波分析考试题(附答案)

1、.word格式,小波分析试题适用范围：硕士研究生时 间：2013年6月一、名词解释(30分)1、线性空间与线性子空间解释：线性空间是一个在标量域(实或复)F上的非空矢量集合 V;设V1是数域K上的线性空间V的一个非空子集合，且对V已有的线性运算满足以下条件(1)如果x>y VI,则x+ y V1;(2)如果V1,K,则kx V1, 则称V1是V的一个线性子空间或子空间。2、基与坐标解释：在n维线性空间V中，n个线性无关的向量 .,；2,., ；n，称为V的一组基；设：是中任一向量，于是1, ；2,., ；n线性相关，因此可以被基；1, ；2，；n, 线性表出：1,；2，；n,二a1；1*

2、1；2. an；n其中系数a1,a1,., an是被向量和基1, ；2,., ；n唯一确定的，这组数就称为在基下的坐标，记为(a1, a1，an )。3、内积解释：内积也称为点积、点乘、数量积、标量积。x = (x!,x2,.,XnT, y = (%, y2,.,yn T, 令 x, = x1y1 x2y2 . xnyn，称 x,y 】为 x 与 y 的内积。4、希尔伯特空间解释：线性 完备的内积空间称为 Hilbert 空间。线性(linearity):对任意f , g H,a, b R, a*f+b*g 仍然 H。完备(completeness ):空间中的任何柯西序列都收敛在该空间之内。

3、内积(inner product ) ： <f, g>,它满足：f =：f，f?,.J , g =Ig1,g2,.,gn ? 时 x,y= X1y1 X22. xm。5、双尺度方程解释： (t)VoVi,(t)Wov所以(t)和*(t)都可以用Vi空间的一个基-(2t-k)z 线性表示：(t)hi (2t-k) () (t)八 gk':(2t-k)(2)kk()=h (与)忌2)(3= g (号)#冷)(41 1并且有 h ()二、hke 出一(5 ,g CO 二 '、gke*' (6)，其中(3 )、( 4)即为双2 k2 k尺度方程。二、简述小波的定义及

4、其主要性质(10分)答：小波(Wavelet)这一术语，顾名思义，“小波”就是小的波形。所谓“小”是指它具 有衰减性；而称之为“波”则是指它的波动性，其振幅正负相间的震荡形式。与Fourier 变换相比，小波变换是时间(空间)频率的局部化分析，它通过伸缩平移运算对信号(函数)逐步进行多尺度细化，最终达到高频处时间细分，低频处频率细分，能自动适应时频信号分析的要求，从而可聚焦到信号的任意细节，解决了Fourier变换的困难问题，成为继Fourier变换以来在科学方法上的重大突破。小波性能除了正交性以外还有光滑性、紧支性、衰减性、对称性以及消失矩和时频窗面积。三、简述小波理论的发展，并结合你所研究

5、的领域，对小波理论在该领域的应用 及发展进行综述。(10分)答：1807年，Fourier提出傅里叶分析 ，1822年发表“热传导解析理论” 论文；1910 年Haar提出最简单的小波；1980,年Morlet首先提出平移伸缩的小波公式，用于地质勘 探；1985年，Meyer和稍后的Daubeichies提出“正交小波基”，此后形成小波研究的高潮； 1988年，Mallat 提出的多分辨分析理论( MRA； Coifman, Meyer 等人在1989年引入了小 波包的概念。基于样条函数的单正交小波基由崔锦泰和王建忠在1990年构造出来。1992年A. Cohe n, I. Daubechhi

6、es等人构造出了紧支撑双正交小波基近年来，一种简明有效的构造小波基的方法-提升方案(LiftingScheme)得到很大的发展和重视，利用提升方案构造的小波被认为是第二代小波。Goodman, Lebrun等人提出的多小波 (Multi-wavelet) 理论，Can des和Donoho等提出的脊小波(Ridgelet )和曲小波(Curvelet) 理论，等等。四、简述连续小波变换的过程。(10分)答：可分成5个步骤， 步骤1:把小波和原始信号的开始部分进行比较；步骤 2:计 算系数c。该系数表示该部分信号与小波的近似程度。系数c的值越大表示信号与小波越相似，因此系数c可以反映这种波形的相

7、关程度；步骤3:把小波向右移，距离为 ，得到的小波函数为，然后重复步骤1和2。再把小波向右移，得到小波，重复步骤1和2。按上述步骤一直进行下去，直到信号结束；步骤4:扩展小波 ，例如扩展一倍，得到的小波函数为；步骤5:重复步骤14。五、阐述多分辨分析的思想并给出MALLA算法的表达式。(10分)答：Meyer于1986年创造性地构造出具有一定衰减性的光滑函数，其二进制伸缩与平移构成L2 (R )的规范正交基，才使小波得到真正的发展。1988年S.Mallat在构造正交小波基时提出了多分辨分析(Multi-Resolution An alysis)的概念，从空间的概念上形象地说明了小波的多分辨率

8、特性，将此之前的所有正交小波基的构造法统一起来，给出了正交小波的构造方法以及正交小波变化的快速算法，即Mallat算法。Mallat算法在小波分析中的地位相当于快速傅立叶变换算法在经典傅立叶分析中的地位。定义：空间L2 ( R)中的多分辨分析是指L2 ( R)满足如下性质的一个空间序列VjZ :( 1)单调性： V4 V0 V1；(2)逼近性：Vj=C, Vj=L2(R)；jriZj Z(3)伸缩性：f(t)Vj= f (2tbVj 1；(4)平移不变性：f(t)Vj=f(t-1)Vj，k Z ；( 5)存在函数g(t) V。，使得g(t -k)k Z构成V。的Riesz基。满足上述个条件的函

9、数空间集合成为一个多分辨分析，如果g(t)生成一个多 分辨分析，那么称 g(t)为一个尺度函数。 关于多分辨分析的理解，我们在这里以一个三层的分解进行 说明，其小波分解树如图所示。从图可以明显看出，多分辨分析只是对低频部分进行进一步分解，而高频部分则不予以考虑。分解的关系为 L2(R) =0 Vjl"：Vj ：-!：Vj 1。另外强调一点这 里只是以一个层分解 进行说明，如果要进行进一步的分解，则可以把低频部分分解成低频部分和高频部分，以下再分解以此类推。在理解多分解分析时，我们必须牢牢把握一点：其分解的最终目的是力求构造一个在频率上高度逼近l2(r)空间的正交小波基，这些频率分辨率

10、不同的正交小波基相当于带宽各异的带通滤波器。从上面的多分辨分析树型结构图可以看出，多分辨分析只对低频空间进行进一步的分解，使频率的分辨率变得越来越高。Mallat算法：通过下面公式 和(2)，可以很快计算出尺度系数和小波系数 cj,k,dj,k ，因此，只要确定Vj空间的初始序列Cj,kkZ，就可以算出任意空间系数和小波系数。公式(1 )和(2)称为离散小波变换的分解公式。Vj(j<J)的所有尺度7亠teZ这就是Mallat重构算法:专业.专注.word格式,DB3小波函数，进行六、(10分)基于MATLAB请自行选择一个一维信号，采用 3尺度的分解与重构。要求(1) 附上源程序；(2)

11、 绘出原始信号以及分解、重构的结果图。答：(1)源程序Load leleccum;S=leleccum(1:100);W= db3'Subplot(621);Plot(s);Title(原始程序)；Dwtmode;lcazpd,cdzpd=dwt(s,w);Lxtzpd=2*le ngth(cazpd) Xzpd=idwt(cazpd,cazpd,w,lx);Subplot(622);plot(xzpd); ITitle( zpd模式重构图);Dwtmode( sym' );casym,cdsym=dwt(s,w);Lxtzpd=2*le ngth(caspd) Xsym=id

12、wt(casym,cdsym,w,lx);Subplot(625);plot(xsym);Title( sym模式重构图);Dwtmode (' spd ');Lxtzpd=2*le ngth(caspd) Xsym=idwt(caspd,cdspd,w,lx);Subplot(626);plot(xspd); |(2)原始信号以及分解、重构的结果图,专业.专注原始图據粗造图像1七、给出一个小波分析用于图像压缩的应用实例。(10分)答：图像压缩可按如下程序进行处理clcclearX= imread ('5.jpg'); % 读入图像figure;image(X)

13、;title('原始图像');disp('压缩前图像X的大小：');whos(X)c,s=wavedec2(X,3,'db5');%寸图像用db5小波进行3层小波分解%取第二层低频高频系数ca仁 appcoef2(c,s,'db5',1);%提取低频系数他取小波分解结构中第一层低频系数和高频系数ch1=detcoef2('h',c,s,1);%水平方向cv仁 detcoef2('v',c,s,1);%垂直方向cd仁 detcoef2('d',c,s,1);%斜线方向.word格式,%

14、分别对各频率成分进行重构a仁 wrcoef2('a',c,s,'db5',1);h仁 wrcoef2('h',c,s,'db5',1);v仁 wrcoef2('v',c,s,'db5',1);d仁 wrcoef2('d',c,s,'db5',1);c1=a1,h1;v1,d1;%显示分解后第一层各频率成分的信息figure;c1=ui nt8(c1);image(c1);title('分解后低频和高频信息');%下面进行图像压缩处理%保留小波分解第一层

15、低频信息，进行图像的压缩%第一层的低频信息即为ca1，显示第一层的低频信息%首先对第一层信息进行量化编码ca1= appcoef2(c,s,'db5',1);ca仁 wcodemat(ca1,440,'mat',0);%改变图像的高度ca1= 0.25*ca1;figure;ca1= ui nt8(ca1*2.5);image(ca1);title('第一次压缩的图像');disp('第一次压缩图像的大小为：');whos('ca1')%保留小波分解第二层低频信息，进行图像的压缩，此时压缩比更大%第二层的低频信息即

16、为ca2，显示第二层的低频信息ca2=appcoef2(c,s,'db5',2);%首先对第二层信息进行量化编码ca2=wcodemat(ca2,440,'mat',0);液变图像的高度ca2=0.125*ca2;figure;ca2=ui nt8(ca2*4.5);image(ca2);title(' 第二次压缩后的图像');disp('第二次压缩图像的大小为：');whos('ca2')ca3=appcoef2(c,s,'db5',3);%首先对第二层信息进行量化编码ca3=wcodemat(c

17、a3,440,'mat',0);液变图像的高度ca3=0.125*ca3;figure;ca3=ui nt8(ca3*4.5);image(ca3);title(' 第三次压缩后的图像');disp('第三次压缩图像的大小为：');whos('ca3')MATLAB显示结果压缩前图像X的大小：Name SizeBytes Class AttributesX 768x1024x32359296 uint8第一次压缩图像的大小为：Name SizeBytes Class Attributesca1 388x516x3600624 uint8第二次压缩图像的大小为：Name SizeBytes Class Attributesca2 198x262x3155628 uint8第三次压缩图像的大小为：Name SizeBytes Class Attributesca3 103x135x341715 uint81002(X)30040050D600原始图傑1G0200300400500600T008009001000700专业.专注分解启低频和高频信息10)0120014002DD 4006009001000 1200 14W) 1600 1800