function [cA,cD] = mydwt(x,lpd,hpd,dim);
% 函數(shù) [cA,cD]=MYDWT(X,LPD,HPD,DIM) 對輸入序列x進行一維離散小波分解，輸出分解序列[cA,cD]
% 輸入?yún)?shù)：x——輸入序列；
%          lpd——低通濾波器；
%          hpd——高通濾波器；
%          dim——小波分解層數(shù)。
% 輸出參數(shù)：cA——平均部分的小波分解系數(shù)；
%           cD——細節(jié)部分的小波分解系數(shù)。

cA=x;       % 初始化cA，cD
cD=[];
for i=1:dim
    cvl=conv(cA,lpd);   % 低通濾波，為了提高運行速度，調(diào)用MATLAB提供的卷積函數(shù)conv()
    dnl=downspl(cvl);   % 通過下抽樣求出平均部分的分解系數(shù)
    cvh=conv(cA,hpd);   % 高通濾波
    dnh=downspl(cvh);   % 通過下抽樣求出本層分解后的細節(jié)部分系數(shù)
    cA=dnl;             % 下抽樣后的平均部分系數(shù)進入下一層分解
    cD=[cD,dnh];        % 將本層分解所得的細節(jié)部分系數(shù)存入序列cD
end


% 以下部分是信號去噪的測試程序，與上述函數(shù)無關(guān)
%     if sym
%         if i<=dim-3
%             dnh=averlize(dnh,1);
%         else if i<=dim-1
%                 dnh=averlize(dnh,2);
%             end
%         end
%     end

% function y=averlize(x,fa);
% N=length(x);
% % aver=sum(x(1:N))/N;
% for i=2:N
%     if abs(x(i)-x(i-1))>=0.7
%         switch fa
%             case 1
%                 x(i)=x(i-1);
%             case 2
%                 x(i)=x(i-1);
%             otherwise
%                 x(i)=x(i);
%         end
%     end
% end
% y=x;