导航:首页 > 五金知识 > 小波分析matlab工具箱

小波分析matlab工具箱

发布时间:2022-12-17 16:33:18

㈠ 用matlab软件怎么进行小波分析

它自己有小波工具箱,wav开头的,里面有很多相关函数,一般直接用就是了。初学建议看帮助里面小波工具箱的demo,可能更好理解。

㈡ 用matlab中工具箱进行小波去噪步骤

matlab读取excel文件比较方便,建议你把数据放到xls文件中保存,然后在matlab中用xlsread这个函数读取出来。版

读取出的数据应该是一权个一维数组了,用plot画出图的话,就是常见的曲线。

然后做小波分解:选用你觉得合适的小波基,例如haar,然后用这个小波基做小波分解,再把高频部分去掉,然后用低频部分还原,就得到了去噪后的信号。

其实你这个问题估计也可以用神经网络或者其它曲线拟合一类手段来解决。具体的情况要根据数据特征来判断。

以上。
专业路过的老狼

㈢ 小波分析在层序地层划分中的应用

1.小波分析简介

20世纪80年代后期至今,一种著名的、在各行各业有重要应用价值的数学理论和方法技术在科学技术界得到了广泛的重视和采用,它就是被誉为“数学显微镜”的小波分析(李世雄,1994)。小波分析的主要功能和特点是,它具有多分辨分析或多尺度分析功能,可以把信号分解成各种不同的尺度成分;它具有很强的局部分析功能,同时具有时间(或空间)域和频率域的局部分析性质,它可自动地通过伸缩、平移聚焦到信号的任一细节对其加以分析(侯遵泽,1998)

(1)小波分析基本原理。小波(wavelet),即小区域的波,是一种特殊的长度有限、平均值为0的波形。它有两个特点:一是“小”,即在时域都具有紧支集或近似紧支集;二是正负交替的波动性。如果用小波和构成傅里叶分析基础的正弦波做比较的话,傅里叶分析所用的正弦波在时间上没有限制,从负无穷到正无穷,但小波则倾向于不规则与不对称。

傅里叶分析是把信号分解到一组相互正交的正弦波上的,也就是基函数,我们可以把基函数看成是度量信号某些特征的一把“尺子”,傅里叶分析度量的就是信号的频谱特征,但是如果这把“尺子”过于规则,有时候就不能十分精确地表达信号蕴含的信息,而在小波分析中,“尺子”换成了规则程度更低的小波函数,从而可以更加有效地表达信号中信息的成分。

小波变换对不同频率在时域上的取样步长是调节性的,即在低频时小波变换的时间分辨率较差,而频率分辨率较高;在高频时小波变换的时间分辨率较高,而频率分辨率较低(图2-13),这正符合低频信号变化缓慢而高频信号变化迅速的特点(胡昌华,1999)。这就构成了利用小波变换进行信号分析的基础。

图2-13 数字信号的小波变换

(2)一维连续小波变换。小波变换实际上是求取信号在各小波函数上的投影值。每个小波函数均由一个母小波函数经过尺度伸缩与时间平移得来的。信号分析的一般思路就是分解与组合,寻找一组最能代表信号特征的函数形式,将信号用这些量来逼近,或者写成这些量的线性组合的形式。

小波分析的思想来源于伸缩和平移方法:对波形的尺度伸缩就是在时间轴上对信号进行压缩与伸展,而时间平移就是指小波函数在时间轴上的波形平行移动。

(3)离散小波变换。由于连续小波变换的伸缩和平移系数是相互独立的,所以通过伸缩和平移得到的各个小波函数之间有一定的相似性,但由于这两个系数之间的独立,就引入了信息的冗余。在分辨率一定的情况下,描述了多余的信息,使得反映信号特征的一些参数相互重叠,给我们的分析带来不便,但这些特点可以用在本身就有自相似性的信号上,可以让我们更清楚地看到信号自身的自相似性。

此外,由于冗余信息的存在,也使得小波逆变换的重构过程不唯一,也就是说,由同一母小波生成的不同的小波变换函数可能重构成同一个信号。为了减少冗余信息,就引入了离散小波变换的概念,其中的伸缩和平移系数是可数的,重构过程用求和的形式给出。如果伸缩和平移系数满足一定的对应关系,则称为二进小波变换(尺度以2的幂的形式给出)。离散小波变换主要是建立在二进制小波变换的基础上的。

测井曲线数据也恰好是离散数据,符合离散变换的要求。在利用小波分析进行层序地层划分时,主要是对测井曲线进行多尺度分解,得到不同尺度下的小波变换图,利用其表现出来的特征来划分不同级次的层序。

2.利用小波分析进行层序地层划分

利用小波分析方法是层序地层划分方法上的一种新的尝试,其目的是尽量减少层序划分过程中的主观因素,依靠地层自身表现出来的客观特征来识别层序、准层序组以及准层序。在我们研究的沉积岩地层中,沉积物的特征可以反映沉积时水体的特征。随着沉积水深的变化,沉积物中多种特征都会相应的发生变化,如放射性物质含量、有机质含量或其他微量元素的含量等,这种变化就会在相应的测井曲线上反映出来。而沉积水深变化受到了多种因素的影响,有长期和短期的旋回,是多个不同周期旋回的叠加,因此测井曲线应该是沉积地层中某种随水深变化的特征的多种频率变化的响应的叠加。也就是说,测井曲线中包含着沉积水深不同周期变化的信息,是多个沉积水深变化周期相互叠加的响应。而小波分析能够将信号分解为不同频率不同周期的旋回,因此,小波分析的特点恰好可以和测井曲线的特点相对应,利用小波分析的方法可以比较准确地将测井曲线中相互叠加的反映水深变化的不同周期的信息分别识别出来,识别出的这些信息就可以用来进行沉积旋回的划分。

同时,小波分析方法还可以帮助解决传统研究方法所不能解决的一些难题,如大段单一岩性地层中的沉积旋回识别。大段单一岩性尤其是大段泥岩、页岩,并不是一个小的沉积旋回里沉积的产物,相反,应是一个相当长时期沉积下来的,但是通过传统的岩性划分方法却很难将其划分开,这就给准层序甚至准层序组的划分造成了困难。小波分析方法可以较好地解决这一问题,利用这种方法可以从测井曲线的细微变化中识别沉积间断和沉积旋回。

(1)测井曲线的选择。不同的曲线具有不同的地质含义,进行相同的变换可能会得到不同的结果。但在研究中通过对GR、AC、COND、电阻率等多条曲线进行小波变换后对比发现,不同测井曲线所得出的变换结果尽管形态上不完全一样,但在旋回的划分上却比较一致(图2-14)。图中曲线a是COND测井曲线经过db5小波变换后的结果,曲线b是同一井段AC曲线变换后的结果。出现这个结果是由于虽然不同的曲线代表着不同的地层响应,会呈现出不同的特征,但地层中各种参数的变化主要受沉积环境的影响,会随着沉积环境的旋回变化呈现出基本一致的旋回特征。这也从一个方面反映了小波变换在沉积旋回划分中的客观性。因此,可以选择目标井的测量精度较高、质量较好的曲线来进行小波变换,进而进行沉积旋回的划分。

图2-14 对COND和AC曲线进行小波变换结果对比

(2)小波的选择。同傅里叶分析不同,小波分析的基(小波函数)不是唯一存在的,所有满足小波条件的函数都可以作为小波函数,那么小波函数的选取就成了十分重要的问题,实际选取小波的标准主要有以下三种。

1)自相似性原则:对二进小波变换(因为在正交小波变换中,取样的方式就是按照小波函数取样的,所以不存在这个问题),如果选择的小波对信号有一定相似性,则变换后的能量就比较集中,可以有效减少计算量。

2)判别函数:针对某类问题,找出一些关键性的技术指标,得到一个判别函数,将各种小波函数代入其中,得到一个最优原则。

3)支集长度:大部分应用选支集长度在5~9之间的小波。因为支集太长会产生边界问题,支集太短不利于信号能量的集中。

但在实际应用中,因为大部分信号的信息量太大,很难找到相应的模式,因此主要是依靠经验来选取。根据前人研究经验及作者对不同函数所做结果的对比,本书采用的是Daubechies小波,阶数为5。

Daubechies小波是由著名小波学者Ingrid Daubechies所创造,她发明的紧支集正交小波是小波领域的里程碑,使得小波的研究由理论转到可行。这一系列的小波简写成dbN,其中N表示阶数。

(3)工作流程。测井曲线能比较准确地反映井旁地层的电性、物性等特征,但往往会受到测井仪器、钻井泥浆等其他非地层因素的干扰,且不同频率的旋回相互叠加,对正确识别和划分沉积旋回造成一定的影响。小波分析能真正消除干扰信号,放大真实信息,按不同频率反映出测井曲线中包含的真正旋回特征,以此划分不同级别层序单元,同时还可以在划分高精度沉积旋回的基础上,与Fischer图解相结合划分出体系域。

MATLAB软件的小波分析工具箱是一种比较常用的工具。MATLAB是Math works公司于1982年推出的一套高性能的数值计算和可视化软件。MATLAB的推出得到了广大专家学者的广泛关注,其强大的扩展功能为各个领域的应用提供了基础。各个领域的专家学者相继推出了MATLAB工具箱,包括信号处理、神经网络、图像处理、小波分析等。其中小波分析工具箱可以满足对测井曲线进行小波变换的需要。

图2-15 小波分析流程图

在对测井曲线进行小波变换时,首先需要对所研究层段的顶底界面进行准确的标定,然后将需要变换的该深度段的测井曲线数值建立单独的文本文件作为原始文件。将原始文件导入后保存成.m格式的信号文件。选择MATLAB软件的wavelet(小波分析)工具箱进行离散一维小波变换,小波类型选择db,阶数为5,最大级数定为12,将上述参数选好后进行分析,即可得到一组12条不同级次的db5小波变换曲线(图2-15)。此外对其进行连续一维小波变换,可以得到小波的频谱分析图,选择合适的最大显示值,根据频谱图上图形的闭合方向可以区分出层序的界面和层序单元(图2-6,图2-7)。

(4)单井分析实例。牛100井位于牛庄洼陷西部,地层以砂泥岩互层为主,岩性变化较快(图2-16)。利用小波分析方法对AC、R25两条测井曲线进行了一维连续小波变换,分别得到其小波变换谱系图,对AC曲线进行了一维离散变换,得到不同阶数的小波,根据与地震、测井及录井岩性资料的对比,选用d11,d9,d7三个层的小波变换曲线分别进行层序、准层序组和准层序的划分。

将传统划分方法所得的结果与小波分析方法所得结果进行对比可以比较明显的看出,在层序和准层序组的划分上两种方法划分的层序单元基本一致,可以相互验证。在准层序级别上的划分,小波分析方法的优势就比较明显地体现了出来,划分的旋回数较多,精度也有提高。这也正是小波分析作为“数学显微镜”的特点所决定的。

从图2-16中小波分析得到的d11曲线可以看出,这一段地层可以划分为两个大的旋回,对应两个层序,谱系图上的特征也比较明显。其中每个大的旋回又可以分为三个次一级的旋回,在d9及谱系图上可以找到相关界面,相当于每个层序划分出三个准层序组,每个准层序组在测井曲线及录井资料上也有较明显的反旋回特征。在进行准层组的划分时,小波分析方法可以划分出肉眼不易识别的旋回,从而提高了划分精度。整段地层一共可以划分为21个准层序,代表不同的沉积旋回。其旋回特征在d7曲线上有较好体现,从谱系图上也可以找到各界面的标志。从测井曲线和岩性上看,基本上每一个准层序都是一个反旋回,代表着一期的水体变换,这也完全符合层序地层学的基本原理。

图2-16 牛100井小波分析资料的层序地层划分

王62井位于牛庄洼陷东部,与牛100井相比,划分出的各层序单元的厚度发生了明显的变化,但数目基本一致,这也证明了小波分析划分层序地层的结果是比较可靠的。通过对AC曲线的小波变换得到AC曲线的小波变换谱系图和小波变换曲线,如图2-17所示。从谱系图和d11曲线上可以将整段地层划分为两个大的旋回,分别对应层序Ⅲ和层序Ⅳ。其中每个层序又可以划分为3个准层序组,在d9曲线上可以看到相应的旋回出现,谱系图上可以找到界面的标志(图2-17)。王62井这一段地层一共可以划分成20个准层序,缺失第一个准层序。各准层序在岩石类型、颜色和测井曲线上基本上可以看出反旋回特征,符合层序地层划分方法。

通过牛100井、王62井的划分可以看出,小波分析方法在砂泥岩互层的地层中有较好的应用效果,可以提高层序划分的精度和准确性。在层序划分中有比较好的可重复性,使得全区的划分结果比较客观和统一,减少了人为因素造成的干扰。

㈣ 求助matlab中的小波工具包

打开matlab软件,进入软件主界面在软件的左下方找到start按钮,点击选择toolbox,然后选择wavelet进入wavemenu界面,选择一维小波中的wavelet1-D并进入 7.选择要处理的信号,界面出现loaded信号,这就是没有去噪前的原9.分析后在左边栏目中出现s,a*,d*,其中s为原信号,a*为近似信号,d*为细节信号11.阈值方法常用的有4种fixed(固定阈值),rigorsure,heusure,minmax根据需要选择,一般情况下rigorsure方式去噪效果较好12.oft(软阈值),hard(硬阈值)一般选择软阈值去噪后的信号较为平滑13.在噪声结构中选择unscaled white noise,因为在工程应用中的噪声一般不仅仅含有白噪声14.在噪声结构下面的数值不要随意改,这是系统默认的去噪幅度16.在此窗口下点击file-save denoised singal,保存输出去噪后的信号18.去噪结束后,把去噪后信号(.mat格式)拖至matlab主界面的workspace中,与原信号一起打包,以便以后计算统计量不会用就查帮助文档啊!waverec函数是不需要你自己加零延拓的,上面的代码完全不知所谓,waverec函数的使用是要依赖wavedec函数得到的CL组构的,CL组构中存放小波系数的数组C本身就已经延拓了,而且你不知道它对数据延拓了多少,延拓的方式有多种根本不是你这样直接加零就行的。我发现你很有才,经常提问和编出一些匪夷所思的问题和代码,不耻下问的精神是好的,但我个人是很不提倡这种做法的,有时间在这打字提问,不如找几本基础参考书看看,不了解就查吗,不明白就往明白搞吗,但看你这些“新奇”的问题和代码真很抓狂,自己对于这些基本问题都懒得琢磨,打着勤奋好学,不耻下问的幌子,太没劲了!哦,看错了,waverec函数是可以用上面的代码的,我看成wrcoef函数了,wrcoef函数可以实现waverec、upwlev和upcoef三个函数的功能之和,所以比waverec函数应用简单,不需要你将其他分量置零,用它实现小波工具箱功能最方便。

㈤ matlab小波分析工具箱时信号加载出现了问题,请高手看下,谢谢~

这是多变量,也就是多个一维信号同时小波分析,如果是分析一个一维信号,直接使用工具箱的第一个按钮即可

㈥ matlab中小波分析工具箱中wrcoef和waverec的区别是什么

这两个函数是与wavedec这个函数紧密相连的。
wavedec就是小波分解,将一个信号分解成指定回层数n,并返回答各层的小波系数。

waverec——它的作用与wavedec相反,即将给定的小波系数一次性完全重建出信号。
wrcoef——这个也是输入小波系数,重建信号。但是它与上面有些区别,区别在于它重建的是原信号在指定层次的,高频或者低频分量。也就是说,这个信号不是原本的信号,而且某个层次上的逼近。

㈦ 如何使用matlab中的工具箱

如果是系统自带的,你可以直接用,如果是外部的或者是自编的你需要先把文件夹拷贝到tools文件夹下,再设置路径。x0dx0aMatlab常用工具箱介绍(英汉对照)x0dx0aMatlab Main Toolbox——matlab主工具箱x0dx0aControl System Toolbox——控制系统工具箱x0dx0aCommunication Toolbox——通讯工具箱x0dx0aFinancial Toolbox——财政金融工具箱x0dx0aSystem Identification Toolbox——系统辨识工具箱x0dx0aFuzzy Logic Toolbox——模糊逻辑工具箱x0dx0aHigher-Order Spectral Analysis Toolbox——高阶谱分析工具箱x0dx0aImage Processing Toolbox——图象处理工具箱x0dx0aLMI Control Toolbox——线性矩阵不等式工具箱x0dx0aModel predictive Control Toolbox——模型预测控制工具箱x0dx0aμ-Analysis and Synthesis Toolbox——μ分析工具箱x0dx0aNeural Network Toolbox——神经网络工具箱x0dx0aOptimization Toolbox——优化工具箱x0dx0aPartial Differential Toolbox——偏微分方程工具箱x0dx0aRobust Control Toolbox——鲁棒控制工具箱x0dx0aSignal Processing Toolbox——信号处理工具箱x0dx0aSpline Toolbox——样条工具箱x0dx0aStatistics Toolbox——统计工具箱x0dx0aSymbolic Math Toolbox——符号数学工具箱x0dx0aSimulink Toolbox——动态仿真工具箱x0dx0aSystem Identification Toolbox——系统辨识工具箱x0dx0aWavele Toolbox——小波工具箱x0dx0ax0dx0a例如:控制系统工具箱包含如下功能:x0dx0a连续系统设计和离散系统设计x0dx0a状态空间和传递函数以及模型转换x0dx0a时域响应(脉冲响应、阶跃响应、斜坡响应)x0dx0a频域响应(Bode图、Nyquist图)x0dx0a根轨迹、极点配置x0dx0ax0dx0a较为常见的matlab控制箱有:x0dx0ax0dx0a控制类:x0dx0ax0dx0a控制系统工具箱(control systems toolbox)x0dx0a系统识别工具箱(system identification toolbox)x0dx0a鲁棒控制工具箱(robust control toolbox)x0dx0a神经网络工具箱(neural network toolbox)x0dx0a频域系统识别工具箱(frequency domain system identification toolbox)x0dx0a模型预测控制工具箱(model predictive control toolbox)x0dx0a多变量频率设计工具箱(multivariable frequency design toolbox)x0dx0ax0dx0a信号处理类:x0dx0a信号处理工具箱(signal processing toolbox)x0dx0a滤波器设计工具箱(filter design toolbox)x0dx0a通信工具箱(communication toolbox)x0dx0a小波分析工具箱(wavelet toolbox)x0dx0a高阶谱分析工具箱(higher order spectral analysis toolbox)x0dx0ax0dx0a其它工具箱:x0dx0a统计工具箱(statistics toolbox)x0dx0a数学符号工具箱(symbolic math toolbox)x0dx0a定点工具箱(fixed-point toolbox)x0dx0a射频工具箱(RF toolbox)x0dx0ax0dx0a1990年,MathWorks软件公司为Matlab提供了新的控制系统模型化图形输入与仿真工具,并命名为Simulab,使得仿真软件进入了模型化图形组态阶段,1992年正式命名为Simulink,即simu(仿真)和link(连接)。matlab7.0里的simulink为6.0版本,matlab6.5里的simulink为5.0版本。x0dx0ax0dx0aMATLAB的SIMULINK子库是一个建模、分析各种物理和数学系统的软件,它用框图表示系统的各个环节,用带方向的连线表示各环节的输入输出关系。x0dx0a启动SIMULINK十分容易,只需在MATLAB的命令窗口键入“SIMULINK”命令,此时出现一个SIMULINK窗口,包含七个模型库,分别是信号源库、输出库、离散系统库、线性系统库、非线性系统库及扩展系统库。x0dx0a1.信号源库x0dx0a包括阶跃信号、正弦波、白噪声、时钟、常值、文件、信号发生器等各种信号源,其中信号发生器可产生正弦波、方波、锯齿波、随机信号等波形。x0dx0a2.输出库x0dx0a包括示波器仿真窗口、MATLAB工作区、文件等形式的输出。x0dx0a3.离散系统库x0dx0a包括五种标准模式:延迟,零-极点,滤波器,离散传递函数,离散状态空间。x0dx0a4.线性系统库x0dx0a提供七种标准模式:加法器、比例环节、积分环节、微分环节、传递函数、零-极点、状态空间。x0dx0a5.非线性系统库x0dx0a提供十三种常用标准模式:绝对值、乘法、函数、回环特性、死区特性、斜率、继电器特性、饱和特性、开关特性等。x0dx0a6.系统连接库包括输入、输出、多路转换等模块,用于连接其他模块。x0dx0a7.系统扩展库x0dx0a考虑到系统的复杂性,SIMULINK另提供十二种类型的扩展系统库,每一种又有多种模型供选择。x0dx0a使用时只要从各子库中取出模型,定义好模型参数,将各模型连接起来,然后设置系统参数,如仿真时间、仿真步长、计算方法等。SIMULINK提供了Euler、RungeKutta、Gear、Adams及专用于线性系统的LinSim算法,用户根据仿真要求选择适当的算法。x0dx0ax0dx0a当然,不同版本的Matlab/Simulink内容有所不同。x0dx0ax0dx0a另外,Simulink还提供了诸如航空航天、CDMA、DSP、机械、电力系统等专业模块库,给快速建模提供了很大的便利。

㈧ matlab中的小波工具箱怎么用,希望能详细介绍

将原始数据文件夹到装有matlab的电脑
打开matlab软件,进入软件主界面
在软件的左下方找到start按钮,点击选择toolbox,然后选择wavelet
进入wavemenu界面,选择一维小波中的wavelet1-D并进入
5.将数据文件(.Mat格式)托到matlab软件主界面的workspace
6.在wavemenu主界面中选择file-load signal或者import from workspace—import signal
7.选择要处理的信号,界面出现loaded信号,这就是没有去噪前的原
始信号
8.右上角选择用于小波分析的小波基以及分解层数并点击analyse开始分析
9.分析后在左边栏目中出现s,a*,d*,其中s为原信号,a*为近似信号,d*为细节信号
10.然后点击denoise去噪
11.阈值方法常用的有4种fixed(固定阈值),rigorsure,heusure,minmax根据需要选择,一般情况下rigorsure方式去噪效果较好
12.oft(软阈值),hard(硬阈值)一般选择软阈值去噪后的信号较为平滑
13.在噪声结构中选择unscaled white noise,因为在工程应用中的噪声一般不仅仅含有白噪声
14.在噪声结构下面的数值不要随意改,这是系统默认的去噪幅度
15.点击denoise开始正式去噪
16.在此窗口下点击file-save denoised singal,保存输出去噪后的信号
17.去噪结束
18.去噪结束后,把去噪后信号(.mat格式)拖至matlab主界面的workspace中,与原信号一起打包,以便以后计算统计量
19.Matlab编程计算相关统计量以及特征量
20.得出统计量和特征量后结束

㈨ MATLAB 小波分析

[c l]=wavedec(x,5,'db6'); %%%x是要分析的数据%%%
d1=wrcoef('d',c,l,'db6',1);
d2=wrcoef('d',c,l,'db6',2);
d3=wrcoef('d',c,l,'db6',3);
d4=wrcoef('d',c,l,'db6',4);
d5=wrcoef('d',c,l,'db6',5);
a5=wrcoef('a',c,l,'db6',5);
plot(d1); %%%输出d1的图像
合并的话,直接用x1=a5+d5+d4+d3就可以吧

阅读全文

与小波分析matlab工具箱相关的资料

热点内容
仪表盘黄色闪灯是什么意思 浏览:325
农村申请健身器材去哪里申请 浏览:372
五金机电黄页 浏览:731
铸钢和铸铁的阀门有什么区别 浏览:136
煤气灶阀门芯坏了怎么办 浏览:615
判断实验装置 浏览:442
小学仪器报告单怎么写 浏览:558
氮气阀门怎么使用 浏览:400
暖气总阀门可以改吗 浏览:266
三叉式万向传动装置的拆装 浏览:296
重庆熔炼设备哪里好 浏览:887
牵引传动装置包括什么 浏览:215
EGO电动工具电池协议 浏览:940
煤气阀门需要每顿关吗 浏览:27
5匹冷库r404加多少制冷剂 浏览:645
天津设备进口代理公司有哪些品牌 浏览:384
图纸管道阀门手柄角度W45N 浏览:569
初中化学实验装置自制 浏览:623
尼尔机械纪元如何用辅助机甩 浏览:471
铸造铝锭都有什么型号 浏览:870