from MELP source code by Oleg Basov
MELP encoder and decoder in Matlab

vbp=melp_bpva(melp_bands,melp_envelopes,p2)
function vbp=melp_bpva(melp_bands,melp_envelopes,p2)
%    
%  : 
%   melp_bands     -   
%   melp_envelopes -    
%   p2             -    
%  :
%   vbp            -  

p2=round(p2); %   
for j=1:4     %    
    k=j+1;    %   (2-5) 
    [p(1),r(1)]=fpr(melp_bands(k,:),p2);%    
                                        %     
    [p(2),r(2)]=fpr(melp_envelopes(j,:),p2);%   
                                            %    
    r(2)=r(2)-0.1;
    %    
    if r(2)>r(1) 
        temp=r(2); 
    else
        temp=r(1);
    end
    %     
    if temp>0.6
        vbp(j)=1; 
    else
        vbp(j)=0;
    end
end
%     ,   
% 
if vbp(1:3)==[0 0 0] 
    vbp(4)=0; 
end

Contact us at files@mathworks.com