dwt
单一一维离散小波变换
描述
(返回单级DWT与指定的扩展模式cA,cD)= dwt (___“模式”,extmode)extmode。有关更多信息,请参见dwtmode。这个观点可以被添加到任何以前的输入语法。
请注意
为gpuArray输入,支持的模式金宝app“symh”(“符号”),“每”。所有“模式”选择除了“每”被转换成“symh”。看这个例子在GPU的单一离散小波变换。
例子
DWT使用小波的名字
获得单层DWT的嘈杂的多普勒信号利用小波的名字。
负载noisdoppgydF4y2Ba;(cA、cD) = dwt (noisdopp,“sym4”);
重建信号的平滑版本使用近似系数。情节和与原始信号。
xrec =得到(cA, 0(大小(cA)),“sym4”);情节(noisdopp)在网格在情节(xrec)传说(“原始”,“重建”)
DWT使用小波和扩展过滤器
获得单层DWT的嘈杂的多普勒信号利用小波高通滤波和扩展(低通滤波器)过滤器。
负载noisdoppgydF4y2Ba;(LoD、藏)= wfilters (“bior3.5”,' d ');[cA、cD] = dwt (noisdopp、LoD、藏);
创建一个可以应用于DWT滤波器组的多普勒信号使用相同的小波。获得的高通滤波和低通滤波器滤波器组。
len =长度(noisdopp);fb = dwtfilterbank (“SignalLength”兰,“小波”,“bior3.5”);(瞧,嗨)=过滤器(神奇动物);
为bior3.5小波,罗和嗨12-by-2矩阵。罗是低通滤波器,嗨高通滤波器。的第一列罗和嗨用于分析和第二列是用于合成。比较的第一列罗和嗨与LoD和藏分别。确认他们是平等的。
disp (“低通滤波器分析”)
低通滤波器分析过滤器
(罗(:1)LoD)
ans =12×2-0.0138 -0.0138 0.0414 0.0414 0.0525 0.0525 -0.2679 -0.2679 -0.0718 -0.0718 0.9667 0.9667 0.9667 0.9667 -0.0718 -0.0718 -0.2679 -0.2679 0.0525 0.0525⋮
disp (“高通滤波分析过滤器”)
高通滤波分析过滤器
(嗨(:1)藏)
ans =12×20 0 0 0 0 0 0 0 -0.1768 -0.1768 0.5303 0.5303 -0.5303 -0.5303 0.1768 0.1768 0 0 0 0⋮
情节片面级一级小波的频率响应和扩展过滤器。
[psidft f phidft] = freqz (fb);水平= 1;情节(f (len / 2 + 1:结束),abs (phidft(水平,len / 2 + 1:结束)))在情节(f (len / 2 + 1:结束),abs (psidft(水平,len / 2 + 1:结束)))网格在传奇(“扩展过滤器”,小波滤波器的)标题(“一级片面的频率响应”)包含(的归一化频率(周期/样本))ylabel (“级”)
在GPU的单一离散小波变换
指GPU的金宝app支持版本(并行计算工具箱)看看支持gpu。金宝app
负载的多普勒信号。把信号在GPU上使用gpuArray。保存当前扩展模式。
负载noisdoppgydF4y2Banoisdoppg = gpuArray (noisdopp);origMode = dwtmode (“状态”,“nodisp”);
使用dwtmode改变扩展模式,零填充。获得的单一信号的离散小波变换在GPU使用db2小波。
dwtmode (“zpd”,“nodisp”)(cA、cD) = dwt (noisdoppg,“db2”);
当前的扩展模式zpd不支持金宝appgpuArray输入。因此,DWT而是进行使用信谊扩展模式。通过DWT的证实了这一点noisdoppg设置为扩展模式信谊并与之前的结果。
[cAsym, cDsym] = dwt (noisdoppg,“db2”,“模式”,“符号”);(max (abs (cA-cAsym))马克斯(abs (cD-cDsym)))
ans = 0 0
一个不支持金宝app的扩展模式指定为转换为一个输入参数“符号”。确认的DWTnoisdoppg与“模式”设置为一个不受支持的模式金宝app也默认了信谊扩展模式。
(cA、cD) = dwt (noisdoppg,“db2”,“模式”,“社会民主党”);(max (abs (cA-cAsym))马克斯(abs (cD-cDsym)))
ans = 0 0
改变当前的扩展模式周期。获得的单一信号的离散小波变换在GPU使用db2小波。
dwtmode (“每”,“nodisp”)(cA、cD) = dwt (noisdoppg,“db2”);
确认当前的扩展模式每支持金宝appgpuArray输入。
[雀跃,cdp] = dwt (noisdopp,“db2”,“模式”,“每”);(max (abs (cA-cAper))马克斯(abs (cD-cDper)))
ans = 0 0
扩展模式恢复到原来的设置。
dwtmode (origMode“nodisp”)
输入参数
输出参数
算法
从一个信号年代的长度N两组系数计算:近似系数cA1,和细节系数cD1。卷积年代与扩展过滤器LoD,其次是二元大量毁灭,收益率近似系数。同样,卷积年代与小波滤波器藏收益率,终于被赶下台,紧随其后的是二元的细节系数。
在哪里
——卷积过滤器X——Downsample(保持偶数元素)
每个滤波器的长度等于2n。如果N=长度(年代),信号F和G的长度N+ 2n−1和系数cA1和cD1地板长度的 。
处理signal-end影响途径造成的算法,一个全局变量管理dwtmode定义所使用的信号扩展模式。可能的选项包括零填充和对称扩展,这是默认的模式。
请注意
对于相同的输入,dwt功能和DWT的块DSP系统工具箱™不能产生相同的结果。DWT块是专为实时实现,而小波工具箱™软件设计进行分析,所以产品处理边界条件和过滤状态不同。下载188bet金宝搏
为了使dwt函数输出匹配DWT块输出,输入函数边界条件设置为零填充dwtmode (“zpd”)在MATLAB®命令提示符。与DWT的延迟块,使用FIR滤波器实现,添加0的输入dwt函数。0你添加的数量必须等于滤波器长度的一半。
引用
[1]Daubechies,我。十个专题小波。在应用数学CBMS-NSF地区会议系列。费城,宾夕法尼亚州:工业与应用数学学会,1992年。
[2]Mallat, s . g .”小波多分辨率信号分解的一个理论:表示。“IEEE模式分析与机器智能。问题7卷。11日,1989年7月,页674 - 693。
[3]Meyer Y。小波和运营商。由d·h·塞林格翻译。英国剑桥:剑桥大学出版社,1995年。
