文件

modwtmra

基于MODWT的多分辨率分析

句法

mra = modwtmra(w)
mra = modwtmra(w,wname)
MRA = MODWTMRA(W,LO,HI)
mra = modwtmra(___,'反射')

描述

例子

MRA.= modwtmra(W.的)返回最大重叠离散小波变换(MODWT)矩阵的多分辨率分析(MRA),W..MODWT矩阵,W.的输出modwt.功能。

例子

MRA.= modwtmra(W.瓦姆姆的)使用对应的小波构造MRA瓦姆姆.这瓦姆姆小波必须是用于获得MODWT的相同的小波。

例子

MRA.= modwtmra(W.你好的)使用缩放滤波器构建MRA和小波过滤器你好.这你好filters必须与获取MODWT所用的filters相同。

例子

MRA.= modwtmra(___,'反射')在构造MRA时使用前面语法中的任何参数使用反射边界条件。如果您指定“反射”modwtmra假设原始信号的长度是输入系数矩阵中列数的一半。

例子

全部收缩

获得了一个简单的时间序列信号的MODWTMRA,并证明了良好的重建。

创建时间序列信号

t = 1:10;x = sin(2 * pi * 200 * t);

获取MODWT和MODWTMRA,并将MODWTMRA行相加。

m = modwt (x);mra = modwtmra (m);xrec =总和(mra);

使用绝对值的最大值以显示原始信号与重建之间的差异非常小。最大的绝对值是按顺序排列$ 10 ^ { -  25} $,这表明了完美的重建。

max(abs(x-xrec))
ans = 5.5738e-25

使用“DB2”小波将ECG信号的MRA构建到四级。

加载wecg;lev = 4;wtecg = modwt(wecg,“db2”列弗);MRA = MODWTMRA(WTECG,“db2”);

绘制ECG波形和MRA。

次要情节(1,1)情节(wecg)为了kk = 2:lev + 2个子图(6,1,kk)plot(MRA(kk-1,:))结尾

构建南部振荡指数数据的多分辨率分析。抽样期是一天。绘制与规模相对应的八级细节$ 2 ^ 8 $天。这个尺度上的细节捕获了大约一年尺度上的振荡。

加载soiwtsoi = modwt(soi);MRASOI = MODWTMRA(WTSOI);plot(mrasoi(8,:))标题('级别8详细信息'的)

使用最小带宽缩放和带有四个系数的小波滤波器获取德国马克-美元汇率数据的MRA。

加载dm_usd.;Lo = [0.4801755, 0.8372545, 0.2269312, -0.1301477];你好= qmf (Lo);wdm = modwt (DM_USD,嗨);mra = modwtmra (wdm,嗨);

使用以下方法获取ECG信号的MRA“反射”边界处理。

加载wecg;wtecg = modwt(wecg,“反射”);MRA = MODWTMRA(WTECG,“反射”);

表明MRA中的列数等于原始信号中的元素数。

isequal(大小(mra), 2),元素个数(wecg))
ans =逻辑1

输入参数

全部收缩

最大重叠离散小波变换,指定为矩阵。W.是输出modwt.

输入W.是一个L.+ 1 -N.包含Modwt的矩阵N.- 点输入信号降至级别L..默认情况下,modwtmra假设你用含有四个消失矩的symlet小波得到了MODWT,'符号4',并使用定期边界处理。

例子:

数据类型:双倍的

合成小波,指定为字符向量。合成小波必须具有用于获得MODWT的相同小波modwt.功能。

缩放滤波器,指定为偶数长度实值向量。您可以指定除非你没有指定瓦姆姆必须是用于获得MODWT的相同缩放滤波器modwt.功能。

小波滤波器,指定为偶数长度实值向量。您可以指定你好除非你没有指定瓦姆姆你好必须是用于获得MODWT的相同小波滤波器modwt.功能。

输出参数

全部收缩

多分辨率分析,作为矩阵返回。默认情况下,MRA.与输入变换矩阵相同W..如果您指定反射边界处理,那么MRA.有一个半数作为输入矩阵的列W.

输出MRA.是一个L.+ 1 -N.矩阵。这K.MRA.的详细信息K.等级。这 (L.+ 1)MRA.包含L.水平流畅。

参考文献

珀西瓦尔,d。B。和瓦尔登,a。T。时间序列分析的小波方法.剑桥,U.K:剑桥大学出版社,2000年出版社。

[2] Whitcher,B.,P.Guttorp和D.B.Cercival。“对大气时间序列的应用对协方差的小波分析。”地球物理研究杂志,卷。105,2000,pp.14941-14962。

扩展能力

也可以看看

|

介绍了R2015b

这个话题有用吗?