lifting_db53.m

来自「是学习时频分析比较基础的一本书」· M 代码 · 共 22 行

function [sn,dn]=db53_lifting(x)
%=====该程序对向量进行小波变换，所用的小波为利用提升方案实现的 Daubechies 9/7小波=====
% 这是一维小波分解
% cn和dn分别是分解得到的低频向量(概貌系数)和高频向量(细节系数)
% x为待分解的向量

L=floor(length(x)/2);
s0=zeros(1,L); d0=zeros(1,L);

for i=1:L                          %进行奇偶二抽取
    se(i)=x(2*i-1);
    so(i)=x(2*i);
end

    dn=so-(se+[se(2:end) se(1)])/2;
    sn=se+(dn+[dn(end) dn(1:end-1)])/4;
    sn=sn*sqrt(2);

if L < length(x)/2           %保持原向量的长度不变
    dn = [dn,dn(L)];
end