[e,state_pul,state_noi,T]=d_mix(fm,T,jt,vp,state_pul,state_noi)

function [e,state_pul,state_noi,T]=d_mix(fm,T,jt,vp,state_pul,state_noi) 
% Description:    
% Inputs:
%     fm -   
%     T -   
%     jt - 
%     vp -  
%     state_pul -     
%     state_noi -     
% Output: 
%     e -  
%     state_pul -     
%     state_noi -     
%     T -   

global melp_firs; %  5 -

%  
T=T*(1+jt*0.5*(rand-0.5)); %     1, 
                           %       
T=round(T);                %    
%       
if T>160
    T=160;
elseif T<20
    T=20; 
end 

%   
    %   -   
    m(1:10)=fm; 
    m(T-10:T-1)=fliplr(fm); 
    if T>21 
        m(11:T-11)=1; 
    end
m=[0,m]; %  

ep=real(ifft(m));           %  ()  
                            %   
ep=[ep(T-9:T),ep(1:T-10)];  %   10   
                            %  
% 
ep=ep*sqrt(T)*1000;      %   
en=(rand(1,T)-0.5)*3464; %   

efirs(1:31,1)=0; 
pfirs=efirs; 
for i=1:5 
    if vp(i)>0.5 
        pfirs=pfirs+melp_firs(:,i); %   
                                    %   
    else
        efirs=efirs+melp_firs(:,i); %   
                                    %   
    end 
end 

ep=[state_pul,ep]; 
en=[state_noi,en]; 
for i=1:T 
    e(i)=sum(ep(i:30+i)*pfirs+en(i:30+i)*efirs); %   
                                                 %  
end
if T<30
    state_pul=[zeros(1,30-T),ep(1:T)];
    state_noi=[zeros(1,30-T),en(1:T)];
else
    state_pul=ep(T-29:T);%    
    state_noi=en(T-29:T);%    
end