baoluoxian.m

经改进的希尔伯特——黄分解

function [ envmax,envmin] = baoluoxian( m,t,indmax,indmin )
%BAOLUO Summary of this function goes here
%  Detailed explanation goes here
add=t(2)-t(1)

if  indmax(1)< indmin(1)           %left 
    k1=indmax(2)-indmax(1)
elseif indmax(1)> indmin(1)
    k1= indmin(2)- indmin(1)
elseif length(indmax)==length( indmin)==1
    k1=2*abs(indmax(1)- indmin(1))
end  
    tt1=t(indmax(1))-2*add*k1       %left extend points
    uu1=m(indmax(1))
    tt2=t(indmax(1))-add*k1
    uu2=m(indmax(1))
    ttt1=t( indmin(1))-2*k1*add
    vv1=m( indmin(1))
    ttt2=t( indmin(1))-k1*add
    vv2=m( indmin(1))
     
    
if  indmin(end)<indmax(end)          %right
    k2=indmax(end)-indmax(end-1)
elseif indmax(end)< indmin(end)
    k2= indmin(end)- indmin(end-1)
elseif length(indmax)==length( indmin)==1
    k2=2*abs(indmax(end)- indmin(end))
end
    
    tt3=t(indmax(end))+add*k2     %right extend points
    uu3=m(indmax(end))
    tt4=t(indmax(end))+2*add*k2
    uu4=m(indmax(end))
    ttt3=t( indmin(end))+add*k2
    vv3=m( indmin(end))
    ttt4=t( indmin(end))+2*add*k2
    vv4=m( indmin(end))
    
    %对数据进处理
    if m(1)>m(indmax(1))
        tt2=t(1)
        uu2=m(1)
    end
    if m(1)<m( indmin(1))
        ttt2=t(1)
        vv2=m(1)
    end
    if m(end)>m(indmax(end))
        tt3=t(end)
        uu3=m(end)
    end
    if m(end)<m( indmin(end))
        ttt3=t(end)
        vv3=m(end)
    end  
        ttmin=[ttt1 ttt2 t( indmin) ttt3 ttt4]
        ttmax=[tt1 tt2 t(indmax) tt3 tt4]
        ttxmin=[vv1 vv2 m( indmin) vv3 vv4]
        ttxmax=[uu1 uu2 m(indmax) uu3 uu4]
        envmax= interp1(ttmax,ttxmax,t,'spline');	

        envmin = interp1(ttmin,ttxmin,t,'spline')