主要内容

qmf

尺度与小波滤波器

    描述

    Y= qmf (X更改反向向量过滤器系数的偶数索引元素的符号X

    例子

    Y= qmf (X,P更改反向向量过滤器系数的偶数索引元素的符号X如果P0.如果P1,奇数索引元素的符号相反。改变P以π弧度改变所得小波滤波器的傅里叶变换相位。

    例子

    全部折叠

    属性关联的正交镜像过滤器db10小波。

    的关联的缩放过滤器db10小波。

    科幻小说= dbwavf (“db10”);

    dbwavf对滤波系数进行归一化,使其范数等于 1 / 2 .对系数进行归一化,使滤波器的范数等于1。

    G = sqrt(2) *科幻小说;

    利用,得到小波滤波系数qmf.情节的过滤器。

    H = qmf (G);次要情节(2,1,1)干细胞(G)标题(“缩放(低通)滤波器G”网格)次要情节(2,1,2)干细胞(H)标题(小波(高通)滤波器H网格)

    图中包含2个轴对象。过滤器G包含一个类型为stem的对象。标题为小波(高通)滤波器H的轴对象2包含一个类型为stem的对象。

    保存当前的扩展模式。设置扩展模式为“周期化”。生成长度为64的随机信号。执行信号的单级小波分解使用GH.为了重现性,将随机种子设置为默认值。

    origmode = dwtmode (“状态”,“nodisplay”);dwtmode (“每”,“nodisplay”n = 64;rng默认的sig = randn (1, n);[d] = dwt (sig, G, H);

    近似系数和细节系数的长度均为32。确认过滤器保存能量。

    [总和(sig。^ 2)和(a . ^ 2) +和(d ^ 2)]。
    ans =1×292.6872 - 92.6872

    计算的频率响应GH.对滤波器进行傅里叶变换时进行零运算。

    n = 128;F = 0:1 / n: 1 - 1 / n;Gdft = fft (G, n);Hdft = fft (H, n);

    绘制每个频率响应的幅度。

    图绘制(F (1: n / 2 + 1), abs (Gdft (1: n / 2 + 1)),“r”)举行情节(F (1: n / 2 + 1), abs (Hdft (1: n / 2 + 1)),“b”网格)标题(“频率响应”)包含(“归一化频率”) ylabel (“级”)传说(“低通滤波器”,“高通滤波器”,“位置”,“东”)举行

    图中包含一个轴对象。标题为“频率响应”的轴对象包含两个类型为line的对象。这些对象代表低通滤波器,高通滤波器。

    确定的频率响应的大小的平方和GH每个频率都等于2。

    sumMagnitudes = abs (Gdft)。^ 2 + abs (Hdft)。^ 2;(最小(sumMagnitudes)最大(sumMagnitudes))
    ans =1×22.0000 - 2.0000

    确认滤波器是标准正交的。

    df = (G; H);id = df * df”
    id =2×21.0000 -0.0000 -0.0000 1.0000

    恢复原来的扩展模式。

    dwtmode (origmode“nodisplay”

    这个例子展示了设置相位参数的效果qmf函数。

    获得与Daubechies小波相关的分解低通滤波器。

    lowfilt = wfilters (“db4”);

    使用qmf函数,得到小波的分解低通滤波器。然后,比较值的符号qmf相位参数设置为0或1。相反的符号表示相移 π 弧度,也就是DFT乘以 e π

    p0 = qmf (lowfilt, 0)
    p0 =1×80.2340 -0.7148 0.6309 0.0280 -0.1870 -0.0308 0.0329 0.0106
    p1 = qmf (lowfilt, 1)
    p1 =1×80.2340 0.2340 -0.6309 -0.0280 0.1870 0.0308 -0.0329 -0.0106

    计算大小并显示它们之间的差异。与相位不同,幅度不受符号反转的影响。

    abs (p0) abs (p1)
    ans =1×80 0 0 0 0 0 0

    输入参数

    全部折叠

    过滤器系数,指定为向量。

    数据类型:|

    相位参数,指定如下。

    • 0-改变反转向量的偶数索引元素的符号X

    • 1-改变反转向量的奇数索引元素的符号X

    数据类型:|

    更多关于

    全部折叠

    正交镜像滤波器

    x做一个有限能量信号。两个过滤器F0F1是否有正交镜滤光片(QMF)x,

    y 0 2 + y 1 2 x 2

    在哪里y0是信号的抽取版本吗x过滤与F0,所以y0被定义为x0F0x),y0n) =x0(2n).同样的,y1被定义为x1F1x),y1n) =x1(2n).这个特性确保了相关的双通道滤波器组方案的完美重构(见[1]p . 103)。

    例如,如果F0是一个范数等于1的Daubechies缩放过滤器吗F1qmfF0),然后是传递函数F0z),F1z)的过滤器F0F1满足条件:

    | F 0 z | 2 + | F 1 z | 2 2.

    参考文献

    [1] Strang, Gilbert和Truong Nguyen。小波和滤波器组.韦尔斯利,马萨诸塞州:韦尔斯利-剑桥出版社,1997年。

    扩展功能

    C / c++代码生成
    使用MATLAB®Coder™生成C和c++代码。

    之前介绍过的R2006a