matlab小波分解与重构 -回复

1 小波变换的基本理论信号分析是为了获得时间和频率之间的相互关系。

小波变换（DWT ）是现代谱分析工具，他既能考察局部时域过程的频域特征，又能考察局部频域过程的时域特征，因此即使对于非平稳过程，处理起来也得心应手。

傅立叶变换提供了有关频率域的信息，但有关时间的局部化信息却基本丢失。

与傅立叶变换不同，小波变换能将图像变换为一系列小波系数，这些系数可以被高效压缩和存储，此外，小波的粗略边缘可以更好地表现图像，因为他消除了DCT 压缩普遍具有的方块效应。

通过缩放母小波（Mother wavelet ）的宽度来获得信号的频率特征， 通过平移母小波来获得信号的时间信息。

对母小波的缩放和平移操作是为了计算小波系数，这些小波系数反映了小波和局部信号之间的相关程度。

小波变换是当前应用数学中一个迅速发展的领域，是分析和处理非平稳信号的一种有力工具。

它是以局部化函数所形成的小波基作为基底展开的，具有许多特殊的性能和优点，小波分析是一种更合理的进频表示和子带多分辨分析。

2小波包变换的基本理论和原理概论：由于正交小波变换只对信号的低频部分做进一步分解，而对高频部分也即信号的细节部分不再继续分解，所以小波变换能够很好地表征一大类以低频信息为主要成分的信号，但它不能很好地分解和表示包含大量细节信息（细小边缘或纹理）的信号，如非平稳机械振动信号、遥感图象、地震信号和生物医学信号等。

与之不同的是，小波包变换可以对高频部分提供更精细的分解，而且这种分解既无冗余，也无疏漏，所以对包含大量中、高频信息的信号能够进行更好的时频局部化分析。

2.1小波包的定义：正交小波包的一般解释 仅考虑实系数滤波器.{}n n Z h ∈{}n n Zg ∈()11nn ng h -=-()()()()22k k Z kk Z t h t k t g t k φφψφ∈∈⎧=-⎪⎨=-⎪⎩为便于表示小波包函数，引入以下新的记号：通过,,h,g 在固定尺度下可定义一组成为小波包的函数。

第6章 小波变换与图像处理
2. 图像的小波分解实例 下面通过两个例子说明如何对图像进行单尺度分解和多 尺度分解，并提取多尺度分解的小波系数。 【例6-1】 对图像进行单尺度分解。 在本例中说明如何对图像进行单尺度分解。程序中调用 函数dwt2对图像进行分解，并画出图像分解的低频分量和水 平、垂直和斜线方向的三个高频分量，可以看出低频分量表 现了图像的轮廓，而高频分量表现了图像的细节。 程序代码如下：
第6章 小波变换与图像处理 subplot(231);image(wcodemat(chd2,nbc)); title('尺度2水平方向的高频系数'); subplot(232);image(wcodemat(cvd2,nbc)); title('尺度2垂直方向的高频系数'); subplot(233);image(wcodemat(cdd2,nbc)); title('尺度2斜线方向的高频系数');
第6章 小波变换与图像处理
2. 图像的平稳小波变换实例 下面举例说明函数swt2的用法。 程序代码如下：
%加载图像 load tire; nbc = size(map,1); colormap(pink(nbc)); cod_X = wcodemat(X,nbc); subplot(221)
第6章 小波变换与图像处理
第6章 小波变换与图像处理
C = [ A(N) | H(N) | V(N) | D(N) | ... H(N-1) | V(N-1) | D(N-1) | ... | H(1) | V(1) | D(1) ]
式中，A为低频系数；H为水平高频系数；V为垂直高频系 数；D为斜线高频系数；所有向量均以列向量存储在矩阵C中。

matlab小波能量
小波分析是一种信号处理方法，它可以将信号分解成不同频率的成分，并分析
这些成分在不同时间点的变化情况。

在MATLAB中，可以使用小波分析工具箱来进行小波变换和重构，并计算小波系数的能量。

以下是使用MATLAB进行小波能量计算的示例代码：
matlab
% 读取信号
x = wavread('signal.wav');
% 选择小波函数
wname = 'db4';
% 进行小波分解
[c,l] = wavedec(x,4,wname);
% 计算小波系数的能量
energy = abs(c).^2;
% 重构信号
x_re = waverec(c,l,wname);
% 绘制原始信号和小波重构信号的波形图
subplot(2,1,1);
plot(x);
title('原始信号');
subplot(2,1,2);
plot(x_re);
title('小波重构信号');
在上述代码中，首先使用wavread函数读取一个音频信号。

然后选择一个小波函数（这里选择的是db4），并使用wavedec函数对信号进行4层小波分解，得到小波
系数c和长度向量l。

接着计算小波系数的能量，这里使用了abs函数计算系数的绝对值，并使用乘方运算得到系数的平方。

最后使用waverec函数对小波系数进行重构，得到重构信号x_re，并绘制原始信号和小波重构信号的波形图。

基于小波变换的信号降噪研究2 小波分析基本理论设Ψt ∈L 2 R L 2 R 表示平方可积的实数空间,即能量有限的信号空间 , 其傅立叶变换为Ψt;当Ψt 满足条件4,7：2()Rt dw wCψψ=<∞⎰1时,我们称Ψt 为一个基本小波或母小波,将母小波函数Ψt 经伸缩和平移后,就可以得到一个小波序列：,()()a bt bt aψ-=,,0a b R a ∈≠ 2 其中a 为伸缩因子,b 为平移因子;对于任意的函数ft ∈L 2 R 的连续小波变换为：,(,),()()f a b Rt bW a b f f t dt aψψ-=<>=3 其逆变换为：211()(,)()fR R t b f t W a b dadb C a aψψ+-=⎰⎰ 4 小波变换的时频窗是可以由伸缩因子a 和平移因子b 来调节的,平移因子b,可以改变窗口在相平面时间轴上的位置,而伸缩因子b 的大小不仅能影响窗口在频率轴上的位置,还能改变窗口的形状;小波变换对不同的频率在时域上的取样步长是可调节的,在低频时,小波变换的时间分辨率较低,频率分辨率较高：在高频时,小波变换的时间分辨率较高,而频率分辨率较低;使用小波变换处理信号时,首先选取适当的小波函数对信号进行分解,其次对分解出的参数进行阈值处理,选取合适的阈值进行分析,最后利用处理后的参数进行逆小波变换,对信号进行重构;3 小波降噪的原理和方法小波降噪原理从信号学的角度看 ,小波去噪是一个信号滤波的问题;尽管在很大程度上小波去噪可以看成是低通滤波 ,但由于在去噪后 ,还能成功地保留信号特征 ,所以在这一点上又优于传统的低通滤波器;由此可见 ,小波去噪实际上是特征提取和低通滤波的综合 ,其流程框图如图所示6：小波分析的重要应用之一就是用于信号消噪 ,一个含噪的一维信号模型可表示为如下形式：(k)()()S f k e k ε=+* k=…….n-1其中 ,f k 为有用信号,sk 为含噪声信号,ek 为噪声,ε为噪声系数的标准偏差;假设ek 为高斯白噪声,通常情况下有用信号表现为低频部分或是一些比较平稳的信号,而噪声信号则表现为高频的信号,下面对 sk 信号进行如图结构的小波分解,则噪声部分通常包含在Cd1、Cd2、Cd3中,只要对 Cd1,Cd2,Cd3作相应的小波系数处理,然后对信号进行重构即可以达到消噪的目的;降噪方法一般来说, 一维信号的降噪过程可以分为 3个步骤进行5,6：1一维信号的小波分解,选择一个小波并确定一个小波分解的层次N,然后对信号进行N 层小波分解计算;2） 小波分解高频系数的阈值量化,对第1层到第N 层的每一层高频系数, 选择一个阈值进行软阈值量化处理．3） 一维小波的重构;根据小波分解的第 N 层的低频系数和经过量化处理后的第1层到第N 层的高频系数,进行一维信号的小波重构;在这 3个步骤中,最核心的就是如何选取阈值并对阈值进行量化,在某种程度上它关系到信号降噪的质量．在小波变换中,对各层系数所需的阈值一般根据原始信号的信号噪声比来选取,也即通过小波各层分解系数的标准差来求取,在得到信号噪声强度后,可以确定各层的阈值;这里着重讨论了信号在两种不同小波恢复后信号质量的不同和对信号中的信号与噪声进行分离;4．仿真实验本文采用Mtalab 本身程序提供的noissin 信号函数及初设原始信号fx 为例进行Matlab 分析1,3,其中：()sin(0.03)f x t =e = noissin + randnsizee1;首先对noissin 函数上叠加上随机噪声信号得到e,分别对比采用db10小波和sym8小波对信号e 进行5层分解,并且细节系数选用minimaxi 阈值模式和尺度噪声db10以及选用sure阈值模式和尺度噪声sym8;在进行噪声消除后,还对原信号进行进一步分析,将原始信号和噪声信号分离开来,仿真结果如图所示：图1图2图3图1-1为原始信号图形,1-2为叠加随机噪声后的图形,而1-3和1-4为利用db10和sym8小波默认阈值降噪后的信号图形;从图1-3和1-4可以看出利用db10和sym8小波降噪后的信号基本上恢复了原始信号,去噪效果明显;但是滤波后的信号与原始信号也有不同,从图中可以很直观地看到采用阈值消噪后信号特征值较少无法准确还原原始信号这是由于为降噪过程中所用的分析小波和细节系数的阈值不恰当所致,如需要更好的恢复信号,还可以采用其它种类小波对其进行分析,通过选取不同的阈值,分析结果,得到一个合适的阈值;从图2和图3中看出,在经过用db10对信号进行5层分解,然后分别对分解的第5层到第1层的低频系数和高频系数进行重构;可以得出其主要基波函数和高频噪声函数的图形,其中小分波分解的细节信号是有白噪声分解得到的,而正弦信号可以在图2中的近似信号a5得到;因为在这一层的影响已经可以忽略了,所以获得的信号就是初始信号的波形,从而把淹没在噪声中的有用信号有效地分离出来;5 总结小波变换对平稳信号的去噪声,要比传统的滤波去噪声得到的效果好．用小波变换进行信号降噪处理, 既降低了噪声同时又提高了信噪比,这说明小波降噪方法是切实可行的方案, 但是由于小波函数很多,采用不同的小波进行分解, 得到的结果可能相差很大, 而变换前并不能预知哪一种小波降噪效果更好,需反复试验比较才能得到良好的效果,这也是小波变换的困难之处之一;另外信号降噪过程中阀值的选取是十分重要的;本文利用两个小波sym8 ,db 10 以及将信号中的信噪分离开来,更加直观可行,通过分别进行信号降噪处理对所得结果与原始信号进行比较可以得出Sym8小波以及默认阈值处理后的重构信号与原始信号最为接近,与分离的结果相同;小波分析是一种信号的视频分析方法,它具有多分辨率分析的特点 ,很适合探测正常信号中夹带的瞬态反常现象并展示其成分,有效区分信号中的突变部分和噪声;通过MATLAB编制程序进行给定信号的噪声抑制和非平稳信号的噪声消除实验表明:基于小波分析的消噪方法是一种提取有用信号、展示噪声和突变信号的优越方法 ,具有广阔的实用价值;在这个越来月信息化的社会中,基于小波分析的应用前景必将越来越广泛;N=10;t=1:10;f=sint.expt+20sint.expt+5sint.expt;plott,f;f=sint.expt+20sint.expt+5sint.expt;输出数据fid=fopen'E:','wt';>> fprintffid,'%f\n',L;C,L=wavedecf,5,'db10';>> fid=fopen'E:','wt';>> fprintffid,'%f\n',L;>> fprintffid,'%f\n',C;>> C,L=wavedecf,1,'db10';>> fid=fopen'E:','wt';>> fprintffid,'%f\n',C;>> C,L=dwtf,'db10';>> fid=fopen'E:','wt';>> fprintffid,'%f\n',C;>> fprintffid,'%f\n',L;参考文献1徐明远,邵玉斌.MATALAB仿真在通信与电子工程中的应用M.西安：西安电子科技大学出版社,2010.2张志涌,杨祖樱等编著.MATLAB教程R2006a-R2007aM.北京：北京航空航天出版社,2006. 3张德丰.详解MATLAB数字信号处理M北京：电子工业出版社,2010.4杨建国.小波分析及其工程应用M北京：机械工业出版社,2005.5冯毅,王香华.小波变换降噪处理及其MATLAB实现J.数字采集与处理,2006,,2112：37-39. 6禹海兰,李天云.基于小波理论的噪声信号分析J.东北电力学院学报.3:36-40.7潘泉,张磊,孟晋丽,张洪才著,小波滤波方法及应用M.北京：清华大学出版社,2005.附仿真源码如下：N=1000;t=1:1000;f=sint;load noissin;e1=noissin;init=66;randn'seed',init;e = e1 + randnsizee1;subplot2,2,1;plott,f;xlabel'1 样本序列'; //x轴标记ylabel'原始信号幅值'; //y轴标记grid ;subplot2,2,2;plote ;xlabel'2 测试样本序列' ;ylabel'含有已加噪声的信号幅值' ;grid ;s1=wdene,'minimaxi','s','one',5,'db12'; subplot2,2,3;plots1;xlabel'3 db10降噪后信号' ;ylabel 'db10小波降噪后的信号幅值';grid;s2=wdene,'heursure','s','one',5,'sym8'; subplot2,2,4;plots2;xlabel'4 sym降噪后信号';ylabel'sym8小波降噪后的信号幅值';grid;figure;subplot6,1,1;plote;ylabel'e';C,L=wavedece,5,'db10';for i=1:5a=wrcoef'a',C,L,'db10',6-i;subplot6,1,i+1; plota;ylabel'a',num2str6-i;endfigure;subplot6,1,1;plote;ylabel'e';for i=1:5d=wrcoef'd',C,L,'db10',6-i;subplot6,1,i+1;plotd;ylabel'd',num2str6-i;end。