Input Args: 

[filename_array, phonemename_array, ph_start_array, ph_end_array]=textread('files.txt','%s %s %u %u','delimiter',' '); 

for vueltas=1:45 

filename_0 = filename_array(vueltas);

phonemename = phonemename_array(vueltas);

ph_start = ph_start_array(vueltas);

ph_end = ph_end_array(vueltas);

filename = strcat(char(filename_0),'.wav');

filename = strcat('outputs\',filename);

x = wavread(char(filename));

fs = 16000;

n = 1:length(x);

 %Plotting the entire word:---------------

fig1 = figure(1);

set(fig1, 'Color','w');

subplot(4,1,1);

plot(n / fs,x);

axis([1/fs length(x)/fs  -1 1]);

title1 = regexprep('/phone/- file', 'phone', phonemename);

title1 = regexprep(title1, 'file', filename_0);

title(title1);

ylabel('x(t)');

xlabel('time(seconds)');

%Plotting the phoneme (in time):-------------------

 phoneme_length = ph_end - ph_start + 1;

subplot(4,2,7);

m = (ph_start:ph_end)/fs;

plot(m,x(ph_start:ph_end)');

title(regexprep('Time waveform for phoneme /phonemename/', 'phonemename', phonemename));

ylabel('s(t)');

xlabel('time(seconds)');

v = axis;

axis([ph_start/fs ph_end/fs  v(3) v(4)]);

 %PLotting Spectrograms:-------------------

subplot(4,1,2);

N = 512;

specgram(x,N,fs,hamming(N),round(0.97*N));

title({'Narrowband'; 'spectrogram'}, 'Position',[length(x)/fs * 1.04 , 4000]);

subplot(4,1,3);

specgram(x,N,fs,hamming(N/4),round(0.9*N/4));

title({'Wideband'; 'spectrogram'}, 'Position',[length(x)/fs * 1.06 , 4000]);

 %Plotting phoneme magnitude spectrum:--------------------------

 subplot(4,2,8);

w = Hamming(length(x));

x_wind = w .* x;

ph_DFT = fft(x_wind(ph_start:ph_end), 4096);

M = length(ph_DFT);

ph_DFT_Mag = abs(ph_DFT);

l = 1:M/2;

plot(l .* (fs /M),20*log10(ph_DFT_Mag(1:M/2)./M));

regexprep('Magnitude spectrum of phoneme /phonemename/', 'phonemename',phonemename)

title(regexprep('Magnitude spectrum of phoneme /phonemename/', 'phonemename',phonemename));

ylabel('|S(f)| (dB)');

xlabel('frequency(Hz)');

%Plotting 3D--------------------------------------------------

fig2 = figure(2);

set(fig2, 'Color','w');

[S,F,T,P] = spectrogram(x,hamming(N),round(0.97*N),N,fs);

surf(T,F,10*log10(abs(P)),'EdgeColor','none');

axis tight;

%set(gca,'Xscale','log');

%set(gca,'XDir','reverse');

set(gca,'YDir','reverse');

ylabel('Freq(Hz)');

xlabel('Time(sec)');

zlabel('Power Density (dB)');

title(regexprep('Power Spectral Density - filename','filename',filename_0)); 

%Saving figures automatically:---------------------------------------- 

saveas(fig1,strcat('bitmaps\',char(filename_0),'_'),'jpg');

saveas(fig2,regexprep('bitmaps\file3D', 'file', filename_0),'jpg'); 

end 

filename_0 = 'excelente';  %filename of the audio file to be analyzed.

filename = strcat(filename_0,'.wav');

filename = strcat('outputs\',filename);

x = wavread(filename);

fl = 300; %low cutoff freq.

fh = 3500; %high cutoff freq.

fs = 16000;

overlap = 75; % overlap percentage

w_length = 512; % Lw = window length (5)

N = 1024; %FFT length

fl_bin = round(fl * N /fs);

fh_bin = round(fh * N /fs);

N_DFT = N - w_length;                                                       

x_length = length(x); %Lx = length of the whole input signal

w = hamming(w_length);                                             

shift_length = ceil(((100 - overlap)/100) * w_length);

p_times = ceil((x_length - w_length)/ shift_length);

w_times = p_times + 1;

new_x_length = p_times * shift_length + w_length; 

x(x_length + 1:new_x_length) = zeros(1,new_x_length - x_length); 

for i=1:w_times

    w_start = ((i - 1) * shift_length) + 1;

    mat(:,i) = x(w_start:w_start + w_length - 1);

    mat(:,i) = mat(:,i) .* w(:,1);

end

mat_defi(1:N_DFT/2,1:w_times) = zeros(1:N_DFT/2,1:w_times); %appending zeros before.

mat_defi((N_DFT/2)+1:(N_DFT/2)+ w_length,1:w_times) = mat(1:w_length,1:w_times); %copying the nonzero values.

mat_defi(w_length+(N_DFT/2)+1:N_DFT + w_length,1:w_times) = zeros(1:N_DFT/2,1:w_times); %appending zeros after.

mat_defi_DFT = fft(mat_defi,N);

mat_defi_mag = abs(mat_defi_DFT);

mat_defi_phase = angle(mat_defi_DFT); 

%Filtering------------------------

mat_defi_mag(1:fl_bin,:) = zeros(fl_bin,w_times);  %zeros at the beggining

mat_defi_mag(N - fl_bin + 2:N,:) = zeros(fl_bin - 1,w_times);  %zeros at the end (mirror)

mat_defi_mag(fh_bin:(N - fh_bin + 2),:) = zeros(N + 3 - 2*fh_bin,w_times);

[OUT_Re, OUT_Imag] = pol2cart(mat_defi_phase,mat_defi_mag);

OUT_DFT = OUT_Re + j * OUT_Imag;

%Resynthesis:-----------------------------------------------------

 out_mat_1 = ifft(OUT_DFT,N)';

new_window_size = N_DFT + w_length;

out_mat_2 = zeros(w_times, (shift_length * p_times) + new_window_size);

 for m = 1:w_times

    out_mat_2(m,(1:new_window_size) + ((m - 1) * shift_length)) = out_mat_1(m,:);   

end

 out = sum(out_mat_2);

out = out ./ max(out);

 %Plots:-------------------------------------------------------------

 figure('Color','w');

%Plotting channel response:

OUT_mag = abs(fft(out));

length_DFT_OUT = (length(OUT_mag));

X_mag = abs(fft(x,length_DFT_OUT ));

r = 1:(length(OUT_mag)/2);

X_log = 20*log10(X_mag(1:(length(OUT_mag)/2))/length(OUT_mag));

plot(r.*(fs /length(OUT_mag)),X_log,'r');

hold on;

OUT_log = 20*log10(OUT_mag(1:(length(OUT_mag)/2))/length(OUT_mag));

plot(r.*(fs /length(OUT_mag)),OUT_log);

axis([1 fs/2 -120 -30]);

title('Spoken Sentence Magnitud Spectrum');

ylabel('Magnitude(dB)');

xlabel('Frequency(Hz)');

legend('original signal','filtered signal'); 

%Saving wavfiles:

wavwrite(out,16000,strcat('outputs\',filename_0,'_filt'));