主要内容

wavemngr

小波经理

描述

使用wavemngr添加、删除、恢复或读取小波。

例子

wavemngr(“添加”,FN烟度WT全国矿工工会文件将小波族添加到工具箱中。这些参数定义了小波族:

  • FN——姓

  • 烟度-姓短

  • WT-小波族类型

  • 全国矿工工会——小波参数

  • 文件-小波定义文件

请注意

当你使用wavemngr为了添加小波族,在当前文件夹中创建了三个小波扩展文件:两个ASCII文件wavelets.ascwavelets.prv,以及mat -文件wavelets.inf.如果你添加了一个新的小波组,它只在这个文件夹中可用。

wavemngr(“添加”,FN烟度WT, {全国矿工工会TYPNUMS},文件加入带参数的小波族全国矿工工会输入格式类型TYPNUMS

wavemngr (___B增加小波族,其中B指定小波的有效支持。万博1manbetx的B输入参数仅对类型的小波有效WT= 3, 4, 5。您可以将此语法与前面的任何语法一起使用。

例子

wavemngr(“▽”,WN删除指定的小波族WN

例子

wavemngr(“恢复”)恢复以前的wavelets.asc美国信息交换标准代码文件

IN2 wavemngr(“恢复”)恢复最初的wavelets.ascASCII文件。在这里IN2是一个虚拟参数。

= wavemngr(阅读)返回字符数组中的所有小波族名称。

= wavemngr(“读”,IN2)返回字符数组中的所有小波名称。在这里IN2是一个虚拟参数。

= wavemngr(“read_asc”)读取wavelets.ascASCII文件,并返回所有小波信息。

例子

全部折叠

列出默认情况下可用的小波族。

wavemngr (“读”
ans =18 x35 char数组'===================================' ' 哈尔- > >哈雾“Daubechies - > > db“Symlets - > - >符号“Coiflets - > >头巾“BiorSplines - > > bior“ReverseBior - > > rbio“迈耶- > > meyr“DMeyer - > > dmey“高斯- > - >“Mexican_hat高斯- > > mexh“Morlet - > > morl“复杂的高斯- > > cgau“香农- > >山“频率b样条- > > fbsp”“复杂Morlet - > > cmor‘Fejer-Korovkin - > >颗  ' '==================================='

列出所有小波。

wavemngr (“读”, 1)
ans =71年x44 char数组'=================================== ' ' 哈尔- > >哈雾  ' '=================================== ' ' Daubechies - > - > db  ' '------------------------------ ' ' db1 - > db2 - > db4 - > db4 - >“db5跑车——> db6 -推出> db7 > db8 - >“db9 - > db10 - > db **-> ' '=================================== ' ' Symlets - > - >符号  ' '------------------------------ ' ' sym2 - > sym3 - > sym4 - > sym5——>“sym6 - > sym7 - > sym8 - >符号 **-> ' '=================================== ' ' Coiflets - > >的发型  ' '------------------------------ ' ' coif1 - > coif2 - > coif3 - > coif4 - >“coif5 -> ' '=================================== ' ' BiorSplines - > > bior  ' '------------------------------ ' ' bior1.1 - > bior1.3 - > bior1.5 - > bior2.2 - >“bior2.4 - > bior2.6 - > bior2.8 - > bior3.1——>“bior3.3 - > bior3.5 - > bior3.7 - > bior3.9——>“bior4.4 - > bior5.5 - > bior6.8 -> ' '=================================== ' ' ReverseBior - > > rbio  ' '------------------------------ ' ' rbio1.1 - > rbio1.3 - > rbio1.5 - > rbio2.2 - >“rbio2.4 - > rbio2.6 - > rbio2.8 - > rbio3.1 - >“rbio3.3 - > rbio3.5 - > rbio3.7 - > rbio3.9 - >“rbio4.4 - > rbio5.5 - > rbio6.8——>“'=================================== ' ' 迈耶- > > meyr  ' '=================================== ' ' DMeyer - > > dmey  ' '=================================== ' ' 高斯- > - >高斯  ' '------------------------------ ' ' gaus1 - > gaus2 - > gaus3 - > gaus4 - >“gaus5 - > gaus6 - > gaus7 - > gaus8 -> ' '=================================== ' ' Mexican_hat - > > mexh”'=================================== ' ' Morlet - > > morl  ' '=================================== ' ' 复杂的高斯- > > cgau  ' '------------------------------ ' ' cgau1 - > cgau2 - > cgau3 - > cgau4 - >“cgau5 - > cgau6 - > cgau7 - > cgau8 -> ' '=================================== ' ' 香农- > >山  ' '------------------------------ '“shan1 - 1.5 - > shan1-1 - > shan1 - 0.5 - > shan1 - 0.1 - >“shan2-3 - >山 **-> ' '=================================== ' ' 频率b样条- > > fbsp  ' '------------------------------ ' ' fbsp1 - 1 - 1.5 - > fbsp1-1-1 - > fbsp1 - 1 - 0.5 - > fbsp2-1-1 - >“fbsp2 - 1 - 0.5 - > fbsp2 - 1 - 0.1 - > fbsp **-> ' '=================================== ' ' 复杂Morlet - > > cmor”'------------------------------ ' ' cmor1 - 1.5 - > cmor1-1 - > cmor1 - 0.5 - > cmor1-1 - >“cmor1 - 0.5 - > cmor1 - 0.1 - > cmor **-> ' '=================================== ' ' Fejer-Korovkin - > >颗  ' '------------------------------ ' ' fk4 - > fk6 - > fk8 - > fk14 - >“fk18 - > fk22 -> ' '=================================== '

这个例子展示了如何将新的紧支持正交小波添加到工具箱中。万博1manbetx这些小波是对Daubechies小波的稍微推广,基于Bernstein多项式的使用,是由Kateb和Lemarié提供的。

添加一组新的正交小波。您必须定义:

  • 姓:Lemarie

  • 家庭短名称:登月舱

  • 小波的类型:1(奥尔特)

  • 小波数字:1 2 3 4 5 5

  • 文件驱动程序:lemwavf

的源代码lemwavf.m在示例的最后提供。的输入参数lemwavf是字符向量的形式吗登月舱N,在那里N= 1、2、3、4或5

wavemngr (“添加”“Lemarie”登月舱的,1,'1 2 3 4 5'“lemwavf”

ASCII文件wavelets.asc保存为wavelets.prv,然后将定义新家庭的信息添加到其中wavelets.asc,以及mat -文件wavelets.inf是生成的。

请注意,wavemngr在当前文件夹上工作。如果你添加了一个新的小波组,它只在这个文件夹中可用。

列出可用的小波族。确定添加了新的小波族。

wavemngr (“读”
ans =19 x35 char数组'===================================' ' 哈尔- > >哈雾“Daubechies - > > db“Symlets - > - >符号“Coiflets - > >头巾“BiorSplines - > > bior“ReverseBior - > > rbio“迈耶- > > meyr“DMeyer - > > dmey“高斯- > - >“Mexican_hat高斯- > > mexh“Morlet - > > morl“复杂的高斯- > > cgau“香农- > >山“频率b样条- > > fbsp”“复杂Morlet - > > cmor‘Fejer-Korovkin - > >颗”“Lemarie登月舱- > >  ' '==================================='

删除添加的家族。再生小波族列表。

wavemngr (“▽”“Lemarie”) wavemngr (“读”
ans =18 x35 char数组'===================================' ' 哈尔- > >哈雾“Daubechies - > > db“Symlets - > - >符号“Coiflets - > >头巾“BiorSplines - > > bior“ReverseBior - > > rbio“迈耶- > > meyr“DMeyer - > > dmey“高斯- > - >“Mexican_hat高斯- > > mexh“Morlet - > > morl“复杂的高斯- > > cgau“香农- > >山“频率b样条- > > fbsp”“复杂Morlet - > > cmor‘Fejer-Korovkin - > >颗  ' '==================================='

恢复之前的ASCII文件wavelets.prv,然后构建mat -文件wavelets.inf.列出恢复的小波。因为wavemngr读取当前工作目录中的ASCII文件,新的系列就会出现在列表中。

wavemngr (“恢复”) wavemngr (“读”, 1)
ans =76年x44 char数组'=================================== ' ' 哈尔- > >哈雾  ' '=================================== ' ' Daubechies - > - > db  ' '------------------------------ ' ' db1 - > db2 - > db4 - > db4 - >“db5跑车——> db6 -推出> db7 > db8 - >“db9 - > db10 - > db **-> ' '=================================== ' ' Symlets - > - >符号  ' '------------------------------ ' ' sym2 - > sym3 - > sym4 - > sym5——>“sym6 - > sym7 - > sym8 - >符号 **-> ' '=================================== ' ' Coiflets - > >的发型  ' '------------------------------ ' ' coif1 - > coif2 - > coif3 - > coif4 - >“coif5 -> ' '=================================== ' ' BiorSplines - > > bior  ' '------------------------------ ' ' bior1.1 - > bior1.3 - > bior1.5 - > bior2.2 - >“bior2.4 - > bior2.6 - > bior2.8 - > bior3.1——>“bior3.3 - > bior3.5 - > bior3.7 - > bior3.9——>“bior4.4 - > bior5.5 - > bior6.8 -> ' '=================================== ' ' ReverseBior - > > rbio  ' '------------------------------ ' ' rbio1.1 - > rbio1.3 - > rbio1.5 - > rbio2.2 - >“rbio2.4 - > rbio2.6 - > rbio2.8 - > rbio3.1 - >“rbio3.3 - > rbio3.5 - > rbio3.7 - > rbio3.9 - >“rbio4.4 - > rbio5.5 - > rbio6.8——>“'=================================== ' ' 迈耶- > > meyr  ' '=================================== ' ' DMeyer - > > dmey  ' '=================================== ' ' 高斯- > - >高斯  ' '------------------------------ ' ' gaus1 - > gaus2 - > gaus3 - > gaus4 - >“gaus5 - > gaus6 - > gaus7 - > gaus8 -> ' '=================================== ' ' Mexican_hat - > > mexh”'=================================== ' ' Morlet - > > morl  ' '=================================== ' ' 复杂的高斯- > > cgau  ' '------------------------------ ' ' cgau1 - > cgau2 - > cgau3 - > cgau4 - >“cgau5 - > cgau6 - > cgau7 - > cgau8 -> ' '=================================== ' ' 香农- > >山  ' '------------------------------ '“shan1 - 1.5 - > shan1-1 - > shan1 - 0.5 - > shan1 - 0.1 - >“shan2-3 - >山 **-> ' '=================================== ' ' 频率b样条- > > fbsp  ' '------------------------------ ' ' fbsp1 - 1 - 1.5 - > fbsp1-1-1 - > fbsp1 - 1 - 0.5 - > fbsp2-1-1 - >“fbsp2 - 1 - 0.5 - > fbsp2 - 1 - 0.1 - > fbsp **-> ' '=================================== ' ' 复杂Morlet - > > cmor”'------------------------------ ' ' cmor1 - 1.5 - > cmor1-1 - > cmor1 - 0.5 - > cmor1-1 - >“cmor1 - 0.5 - > cmor1 - 0.1 - > cmor **-> ' '=================================== ' ' Fejer-Korovkin - > >颗  ' '------------------------------ ' ' fk4 - > fk6 - > fk8 - > fk14 - >“fk18 - > fk22 -> ' '=================================== ' ' 登月舱Lemarie - > >  ' '------------------------------ 'lem1 - > lem2 - > lem3 - > lem4——>“lem5 -> ' '=================================== '

恢复初始小波。恢复初始ASCII文件wavelets.ini和初始的mat文件wavelets.bin.再生小波族列表。这个名单不包括这个新家庭。

wavemngr (“恢复”, 0) wavemngr (“读”
ans =18 x35 char数组'===================================' ' 哈尔- > >哈雾“Daubechies - > > db“Symlets - > - >符号“Coiflets - > >头巾“BiorSplines - > > bior“ReverseBior - > > rbio“迈耶- > > meyr“DMeyer - > > dmey“高斯- > - >“Mexican_hat高斯- > > mexh“Morlet - > > morl“复杂的高斯- > > cgau“香农- > >山“频率b样条- > > fbsp”“复杂Morlet - > > cmor‘Fejer-Korovkin - > >颗  ' '==================================='

所有的命令行功能都可以用于新的小波族。建立一个新的家庭。计算四个相关的滤波器和尺度和小波函数。

wavemngr (“添加”“Lemarie”登月舱的,1,'1 2 3 4 5'“lemwavf”);[Lo_D, Hi_D Lo_R Hi_R] = wfilters (“lem3”);[φ,psi, xval] = wavefun (“lem3”);情节(xval(φ;psi));传奇(“扩展功能”“小波”网格)

图中包含一个坐标轴。轴线包含2个线型对象。这些对象代表尺度函数、小波。

删除添加的家族。

wavemngr (“▽”“Lemarie”

lemwavf.m

函数F = lemwavf (wname)Lemarie小波滤波器。% F = LEMWAVF(W)返回与Lemarie相关联的缩放滤波器%小波由字符数组'lemN'指定。% N的取值为1、2、3、4或5。%这个函数只用于“添加小波族”的例子。它%可能会在未来的版本中更改或删除。版权所有2019 The MathWorks, Inc.TFlem = startsWith (wname,登月舱的);如果~ TFlem错误(小波的缩写是lem后跟滤波器号);结束fnum = regexp (wname,“(\ d +)”“匹配”“一次”);如果isempty (fnum)错误('指定过滤器编号为1,2,3,4,或5');结束如果~isempty(fnum) num = str2double(fnum);结束Tffilt = ismember(num,[1 2 3 4 5]);如果~ tffilt错误('过滤器编号必须为1、2、3、4或5');结束开关全国矿工工会情况下1 =中立...0.46069299844871 0.53391629051346 0.03930700681965 -0.03391629578182...];情况下2 .厨师...0.31555164655258 0.59149765057882 0.20045477817080 0.10034811856888...0.01528128420694 0.00846362066021 -0.00072514051618 0.00038684732960...];情况下3 .单词pilot联想记忆:...0.2310894231941 0.56838231367966 0.33173980738190 -0.09447000132310...-0.00209952890579 0.00001769381066...0.00128429679795 -0.00053703458679 0.00002283826072 -0.00000928544107...];情况下4 .单词conduct联想记忆:...0.17565337503255 0.52257484913870 0.42429244721660 -0.04601056550580...-0.11292720306517 0.03198741803409 0.00813124691980 -0.00743764392677...0.00548090619143 -0.00140066128481 -0.00054200083128 0.00025607264164...-0.00008795126642 0.00003025515674 -0.00000082014466 0.00000027569334...];情况下5 .选[f]。...0.13807658847623 0.47310642622099 0.48217097800239 0.02112933622031...-0.15081998732499 0.01935767268926...0.00671209165995 0.00120022744496 -0.00321203819186 0.00115266788547...-0.00018266213413 -0.00002953360842 0.00008433396295 -0.00002997969339...0.00000534552866 -0.00000159098026 0.00000003069431 -0.00000000895816...];结束

这个示例展示了如何使用与双正交小波相关的分析和合成滤波器,并使它们与小波工具箱™兼容。小波工具箱要求分析和合成低通和高通滤波器具有相等的偶数长度。这个例子使用了基于Burt和Adelson的拉普拉斯金字塔方案的近似正交双正交小波([1]第283页表8.4)。该例子还演示了如何检验双正交小波的性质。

定义了双正交小波的分析与综合滤波系数。

Hd =[-1 5 12 5 -1]/20*根号(2);Gd =[3 -15 -73 170 -73 -15 3]/280*根号(2);Hr =[-3 -15 73 170 73 -15 -3]/280*根号(2);Gr =[-1 -5 12 -5 -1]/20*根号(2);

高清Gd分别为低通和高通分析滤波器。人力资源Gr是低通和高通综合滤波器。它们都是有限脉冲响应(FIR)滤波器。确定低通滤波器系数和为√6 (2)高通滤波器系数和为0。

总和(Hd) /√(2)
ans = 1.0000
sum (Hr) /√(2)
ans = 1.0000
总和(Gd)
ans = -1.0061 e-16
总和(Gr)
ans = -9.7145 e-17

z-transform的FIR滤波器 h 是一个洛朗多项式 h z 给出的 h z k k b k e h k z - k .的程度 | h | 的定义为 | h | k e - k b .因此,过滤器的长度 h 1 + | h | .检验尺度和小波滤波器的洛朗展开式。

博士= laurpoly(高清“dmin”, 2)
博士(z) = - 0.07071 * z ^ (+ 2) + 0.3536 * z ^ (+ 1) + 0.8485 + 0.3536 * z ^ (1) - 0.07071 * z ^ (2)
PHr = laurpoly(人力资源,“dmin”3)
PHr (z) =…- 0.01515 * z ^ (+ 3) - 0.07576 * z ^ (+ 2) + 0.3687 * z ^ (+ 1) + 0.8586 + 0.3687 * z ^(1)……- z^(- 1) - z^(- 1)
PGd = laurpoly (Gd,“dmin”3)
PGd (z) =…+ 0.01515 * z ^ (+ 3) - 0.07576 * z ^ (+ 2) - 0.3687 * z ^ (+ 1) + 0.8586 - 0.3687 * z ^(1)……- z^(- 1) + z^(-2)
PGr = laurpoly (Gr,“dmin”, 2)
PGr (z) z ^ = - 0.07071 * 0.3536 (+ 2) * z ^ (+ 1) + 0.8485 - 0.3536 * z ^ (1) - 0.07071 * z ^ (2)

由于滤波器与双正交小波相关联,确定 P H d z P H r z + P G z P G r z 2

博士* PHr + PGd * PGr
ans (z) = 2

小波工具箱要求与小波相关联的滤波器具有相等的长度。要在工具箱中使用拉普拉斯小波滤波器,您必须将洛朗级数的缺失幂包含为零。

度的博士学位PHr分别是4和6。可容纳四个滤波器的最小偶长滤波器长度为8,对应于7次的洛朗多项式。策略是尽可能均匀地添加0,使所有过滤器的长度为8。给所有过滤器添加0,然后添加两个0高清Gr

高清=[0高清0 0];Gd = [0 Gd];Hr = [0 Hr];Gr = [0 Gr 0 0];

你可以通过创建DWT滤波器组来检查双正交小波的特性。使用这些过滤器创建两个自定义DWT滤波器组,一个用于分析,另一个用于合成。确认滤波器组是双正交的。

fb = dwtfilterbank (“小波”“自定义”...“CustomScalingFilter”(高清“人力资源”),...“CustomWaveletFilter”, (Gd的Gr));fb2 = dwtfilterbank (“小波”“自定义”...“CustomScalingFilter”(高清“人力资源”),...“CustomWaveletFilter”(Gd的Gr),...“FilterType”“合成”);流('fb: isOrthogonal = %d\tisBiorthogonal = %d\n'...isOrthogonal (fb) isBiorthogonal(神奇动物));
fb: is正交= 0 isbi正交= 1
流('fb2: isOrthogonal = %d\tisBiorthogonal = %d\n'...isOrthogonal (fb2) isBiorthogonal (fb2));
fb2: is正交= 0 isbi正交= 1

在最粗的尺度上绘制与滤波器组相关的尺度和小波函数。

[φ,t] = scalingfunctions (fb);[psi, ~] =小波(神奇动物);[phi2 ~] = scalingfunctions (fb2);[psi2 ~] =小波(fb2);次要情节(2,2,1)情节(t,φ(最终,:))网格标题(“标度函数-分析”子图(2,2,2)plot(t,psi(end,:))网格标题(“小波分析”)子plot(2,2,3) plot(t,phi2(end,:))网格标题(“缩放函数-综合”)子图(2,2,4)plot(t,psi2(end,:)标题(“小波合成”

图中包含4个轴。标题为“缩放函数分析”的轴1包含一个类型为line的对象。标题为小波分析的坐标轴2包含一个类型为线的对象。标题为“缩放函数-合成”的轴3包含一个类型为线的对象。标题为小波合成的轴4包含一个线型对象。

计算过滤器组框架边界。

[analysisLowerBound, analysisUpperBound] = framebounds(神奇动物)
analysisLowerBound = 0.9505
analysisUpperBound = 1.0211
[synthesisLowerBound, synthesisUpperBound] = framebounds (fb2)
synthesisLowerBound = 0.9800
synthesisUpperBound = 1.0528

输入参数

全部折叠

小波族名称,指定为字符向量或字符串标量。

小波族的短名称,指定为字符向量或字符串标量。中的字符数烟度必须小于或等于4。

小波族类型,指定为以下类型之一:

  • 1——正交小波

  • 2——双正交小波

  • 3.-带尺度函数的小波

  • 4-无尺度函数的小波

  • 5-无尺度函数的复小波

小波参数,具体为:

  • 如果这个族是由一个小波组成的,全国矿工工会为空字符串''.例如,mexhmorl每个族都包含一个小波。

  • 如果小波是有限小波族的成员,全国矿工工会包含用空格分隔的表示小波参数的项列表。例如,对于双正交小波族biorNums = '1.1 1.3 1.5 2.2 2.4 2.6 2.8 3.1 3.3 3.5 3.7 3.9 4.4 5.5 6.8'

  • 如果小波是无限小波族中的一员,全国矿工工会包含用空格分隔的表示小波参数的项列表,以特殊序列结束**.下表中列出了两个示例。

    小波家族 全国矿工工会
    db Nums = '1 2 3 4 5 6 7 8 9 10 **'
    Nums = '1-1.5 1-1 -0.5 1-0.1 2-3 **'

小波参数输入格式,具体为:

  • “整数”—当参数为整数时,使用此选项。例如,Daubechies小波族db使用整型参数。

  • “真实”的—当参数为“real”时,使用该选项。例如双正交小波族bior使用实参数。

  • “charactervector”—当参数为字符向量时,使用此选项。例如,香农小波族使用字符向量。

小波定义文件,指定为字符向量或字符串标量。文件是定义小波族的mat文件或代码文件的名称。

有效支持具有家族类型万博1manbetx的小波WT等于3,4或5,指定为二元素实值向量如果[标准答案][B],然后指定下界,并且乌兰巴托指定上限。

数据类型:

小波族,由字符向量或字符串标量指定WN.的价值WN是小波族名或小波族短名。

例子:wavemngr(“▽”、“Lemarie”)

限制

  • wavemngr允许您添加小波。您必须验证它确实是一个小波。不进行任何检查以确认添加的是小波或确认新小波的类型。您可以使用dwtfilterbank验证小波是正交还是双正交。

参考文献

[1] Daubechies,我。小波十讲.CBMS-NSF应用数学区域会议系列。费城:工业和应用数学学会,1992。

之前介绍过的R2006a