基于Matlab的小波信号的去噪仿真.docxVIP

基于Matlab的小波信号的去噪仿真1.引言小波分析是在傅里叶基础上发展起来的一个有效的调和分析工具，傅里叶分析作为最早的调和分析工具，在信号频谱分析中有着非常重要的作用，是调和分析发展的一个重要分支。但是在当需要对信号的某一时间点，或某一时间段的频率成分进行分析时，传统的傅里叶变换就显得无能为力了，所以在他得基础上发展了短时傅里叶变换，Gabor变换，以致到后来具有革命性意义的小波变换。小波变换是采用面积固定不变但形状不断变化的分析窗口对信号进行变换，其多分辨分析的特点，很适合于分析非平稳信号。语音信号是一个典型的非平稳信号，目前，小波变换已经成功运用于信号处理。2.小波定义（1）构造一个一元函数 ，生成函数系，由其产生的级数可以表示一类定义于R的非周期函数。（2）上述级数不同层次的部分和可以给出函数不同尺度的近似。满足以上要求的函数 即为一元小波函数，简称小波。3.构造小波的过程3.1多分辨分析理论定义令 是一个R-函数，那么ψ称为是一个正交小波，如果满足ψ称为是一个半正交小波，如果满足定义的闭子空间的嵌套序列称为形成一个多分辨分析，满足 1) 且 2)则，且 3) 存在使形成 一个Riesz基。由于是 的一个Riesz基，由条件(2)， 也是的一个Riesz基。所以，也称生成一个多分辨分析。设生成一个多分辨分析，由于而，所以，可用 的基底表示。由于是的一个Riesz基，所以存在唯一序列,使上式就称为的两尺度关系, 称为两尺度序列.对两尺度关系两边作Fourier变换，则得两尺度关系的Fourier变换形式前边我们引入了序列的符号,为了突出它与ω的联系,再写规范正交尺度函数的Fourier变换满足下述条件:几乎处处成立。几乎处处成立。对于生成的多分辨分析，由于现在考虑关于补空间，即至少要求满足满足上面两个性质的称为的直接和为了构造，像生成一样，生成的’小波’，由于，所以考虑序列及由它定义的中的函数。3.2构造小波的一般方法如果正交尺度函数生成的一个多分辨分析是中关于的正交补空间，那么，存在的小波基使对有并且,对于小波ψ构造的一种可能是：或等价地所以构造小波的方法是由方程求出，带入到式子中求出小波。3.3 Daubechies小波的构造给出了序列就完全确定了尺度函数与小波函数。假定正交尺度函数是紧支撑的,由容易得到得出(A)由正规性,则写为形如 （N1）(B)把(B)代入(A),由此得出(C)满足(D)(E)为方便,能够再写为的多项式(F)借助于R,记 ,约束式(D)变成 (G)对于所有成立,因此对所有成立。在式子(B)中如且仅如能够写为 而(S)且T 是一个奇多项式,使对 成立。 其中 (T)Daubechies构造了一簇正交小波,她是在(S)中令 这时。再注意式(C)与(E),得到(P)对每个N,相应的 有2N个非零系数,可以写(Q)求出则相应的Daubechies小波由下式给出：4.小波的分解与重构简介对于一个多分辨分析，是关于的补空间，则就可分解为空间 的直接和所以，对于都有唯一分解()令，则并且和的两尺度关系，和的分解关系其中记值得注意的是，在中,k代表分解水平.所以有唯一分解固定k,由求的算法称为分解算法由得分解算法固定k，由 求 的算法称为重构算法由 且得在实际计算中，假定取型值点所对应的f(t)的水平为N，即对于某个正整数，信号由水平N分解到水平N－M水平，即求及及而求得的 是中间过程。5.小波分解重构去噪仿真5.1任务的提出对含噪的信号去噪，实现信噪分离。含噪的正弦信号：5.2设置过滤噪声的阈值由分解的高频系数，可设阈值为根据阈值，可从高频系数中过滤噪声，从而得出目标信号。5.3仿真去噪代码% 生成含噪正弦信号N=1024;t=1:N;sig=sin(0.05*t);figure(1);subplot(211);plot(t,sig); title(正弦信号);% 叠加噪声x=sig+randn(1,N);subplot(212);plot(t,x); title(含噪正弦信号);% 一维小波分解，使用db3进行4层分解[c,l]=wavedec(x,4,db3);% 重构第1-4层逼近信号a4=wrcoef(a,c,l,db3,4);a3=wrcoef(a,c,l,db3,3);a2=wrcoef(a,c,l,db3,2);a1=wrcoef(a,c,l,db3,1);% 显示各层逼近信号figure(2);subplot(411);plot(a4);ylabel(a4);subplot(412);plot(a3);ylabel(a3);subplot(413);plot(a2);ylabel(a2);subplot(414);plot(a1);ylabel(a1);% 重构第1-4层细