function [cn,dn]=lifting_db97(x)
%=====該程序對向量進行小波變換，所用的小波為利用提升方案實現的 Daubechies 9/7小波=====
% 這是一維小波分解
% cn和dn分別是分解得到的低頻向量(概貌系數)和高頻向量(細節系數)
% x為待分解的向量

L=floor(length(x)/2);
s0=zeros(1,L); d0=zeros(1,L);
alpha = -1.586134342; beta = -0.05298011854; gamma = 0.8829110762; delta = 0.4435068522;
K = 1.149604398;                   %設定參數值

for i=1:L                          %進行奇偶二抽取
    s0(i)=x(2*i-1);
    d0(i)=x(2*i);
end

s1=s0(2:L); 
s1=[s1,s0(L)]; 
d0=d0+alpha*(s0+s1);
d1=d0(1:L-1); 
d1=[d0(1),d1];
s0=s0+beta*(d0+d1);
s1=s0(2:L);
s1=[s1,s0(L)];
d0=d0+gamma*(s0+s1);
d1=d0(1:L-1);
d1=[d0(1),d1];
s0=s0+delta*(d0+d1);

cn=K*s0; dn=d0/K;            %求出 cn 和 dn

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