%**************************************************************************************************************************************************************
%   本程序?yàn)楹瘮?shù)功能，目的是對(duì)輸入信號(hào)進(jìn)行小波分解，并且對(duì)分解后的各層細(xì)節(jié)信號(hào)和最終的逼近信號(hào)進(jìn)行重構(gòu)，函數(shù)使用’db3’二進(jìn)制小波基函數(shù)
%      
%   中國(guó)北方發(fā)動(dòng)機(jī)研究所試驗(yàn)測(cè)試中心  蘇志霄   2006年10月
%**************************************************************************************************************************************************************
%
function R=xiaobofenxi(data,n);           %n為小波分解的尺度
len=max(size(data));
R=zeros(n+1,len);
[c,l]=wavedec(data,n,'db3');
%對(duì)分解結(jié)構(gòu)[c,l]中的第n逼近部分進(jìn)行重構(gòu)
R(1,:)=wrcoef('a',c,l,'db3',n);
%對(duì)分解結(jié)構(gòu)[c,l]中的各層細(xì)節(jié)部分進(jìn)行重構(gòu)
for i=1:n;
    R(i+1,:) = wrcoef('d',c,l,'db3',n+1-i);
end;