//silence removal and end point detection

THRESHOLD=0.3; // adjust value yourself

TIME=;

%capture;

Fs = ;

y = wavrecord(TIME*Fs,Fs);

%plot(y)

%wavplay(y,Fs);

samplePerFrame=floor(Fs/);

bgSampleCount=floor(Fs/); %according to formula,  sample needed for  khz

%----------

%calculation of mean and std

bgSample=[];

for i=::bgSampleCount

    bgSample=[bgSample y(i)];

end

meanVal=mean(bgSample);

sDev=std(bgSample);

%----------

%identify voiced or not for each value

for i=::length(y)

   if(abs(y(i)-meanVal)/sDev > THRESHOLD)

       voiced(i)=;

   else

       voiced(i)=;

   end

end

% identify voiced or not for each frame

%discard insufficient samples of last frame

usefulSamples=length(y)-mod(length(y),samplePerFrame);

frameCount=usefulSamples/samplePerFrame;

voicedFrameCount=;

for i=::frameCount

   cVoiced=;

   cUnVoiced=;

   for j=i*samplePerFrame-samplePerFrame+::(i*samplePerFrame)

       if(voiced(j)==)

           cVoiced=(cVoiced+);

       else

           cUnVoiced=cUnVoiced+;

       end

   end

   %mark frame for voiced/unvoiced

   if(cVoiced>cUnVoiced)

       voicedFrameCount=voicedFrameCount+;

       voicedUnvoiced(i)=;

   else

       voicedUnvoiced(i)=;

   end

end

silenceRemovedSignal=[];

%-----

for i=::frameCount

    if(voicedUnvoiced(i)==)

    for j=i*samplePerFrame-samplePerFrame+::(i*samplePerFrame)

            silenceRemovedSignal= [silenceRemovedSignal y(j)];

        end

    end

end

%---display plot and play both sounds

    figure; plot(y);

    figure; plot(silenceRemovedSignal);

%%%play

wavplay(y,Fs);

wavplay(silenceRemovedSignal,Fs);