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