这个例子展示了如何使用双树复小波变换创建近似解析小波。为了获得近似解析小波变换并得到双树变换的好处,必须仔细构造两个滤波器组中的FIR滤波器。
创建长度为256个样本的零信号。获得0信号到5级的两个双树变换。默认情况下,dualtree
使用默认的近对称双正交滤波器对nearsym5_7
对于电平1,正交Q-shift Hilbert小波滤波器对长度为10,对于大于1的电平。
x = 0 (256 1);[a1, d1] = dualtree (x,“水平”5);(a2, d2) = dualtree (x,“水平”5);
在这两棵树中设置一个5级细节系数为1,并将变换反变换得到小波。
d1 {5} (5) = 1;d2{5}(5) = 1我;wav1 = idualtree (a1, d1);wav2 = idualtree (a2, d2);
用第一棵树作为实部,第二棵树作为虚部形成复小波。绘制小波的实部和虚部。
analwav = wav1 + 1我* wav2;情节(真实(analwav))在情节(图像放大(analwav),“r”)情节(abs (analwav),“k”,“线宽”2)轴紧传奇(实部的,“虚部”,“级”,“位置”,“西北”)
对解析小波进行傅里叶变换并绘制幅值。
zdft = fft (analwav);domega =(2 *π)/长度(analwav);ω= 0:domega:(2 *π)-domega;clf;情节(ω,abs (zdft))包含(“弧度/样本”)轴紧
小波的傅里叶变换实际上只支持频率轴的一半。金宝app