%Matlab file used for class demo conducted during Session 8
%on the effects of aliasing.
%A word utterance is effectively sampled at four different
%rates and played back.  A lowpass filtering is done
%prior to decimation to avert aliasing
clf
set(0,'defaultaxesfontsize',16);

%these commands read in the speech file
data=load('enf1s1t0');
Fs=8192;

%plot data to cut off silence
plot(data)
pause
%end1=input('end1?')
%end2=input('end2?')
%datar=data(end1:end2);
datar=data;
dsize=length(datar);

%compute 8192 point FFT to view spectrum of speech signal
deltaf=Fs/8192;
freq=-Fs/2:deltaf:Fs/2-deltaf;
plot(freq,abs(fftshift(fft(datar,8192))),'linewidth',2)
xlabel('Analog Frequency (Hz)')
ylabel('Spectral Magnitude')
domega=2*pi/8192;
omega=-pi:domega:pi-domega;
pause
input('Playback speech at original rate of 8192 Hz.');
soundsc(datar,Fs)
disp('Reduce sampling rate by a factor of 2 thru decimation');
input(' and playback speech at 4096 KHz.');
datar2=resample(datar,1,2);
soundsc(datar2,Fs/2)
disp('Reduce sampling rate by a factor of 3 thru decimation');
input(' and playback speech at 2731 KHz.')
datar3=resample(datar,1,3);
soundsc(datar3,Fs/3)
disp('Reduce sampling rate by a factor of 4 thru decimation');
input(' and playback speech at 2048 KHz.');
datar4=resample(datar,1,4);
soundsc(datar4,Fs/4)

%compute via the FFT the DTFT of each of the
%four sampled versions of the original speech
%Plot magnitude of each DTFT over -pi to pi
subplot(221)
dr1=abs(fftshift(fft(datar,8192)));
plot(omega,dr1,'Linewidth',3)
axis([-pi pi 0 max(dr1)])
title('Fs=12.5 KHz')
%xlabel('omega (radians)')
subplot(222)
dr2=abs(fftshift(fft(datar2,8192)));
plot(omega,dr2,'Linewidth',3)
axis([-pi pi 0 max(dr2)])
title('Fs=6.25 KHz')
%xlabel('omega (radians)')
subplot(223)
dr3=abs(fftshift(fft(datar3,8192)));
plot(omega,dr3,'Linewidth',3)
axis([-pi pi 0 max(dr3)])
title('Fs=4.167 KHz')
xlabel('\omega (radians)')
subplot(224)
dr4=abs(fftshift(fft(datar4,8192)));
plot(omega,dr4,'linewidth',2)
axis([-pi pi 0 max(dr4)])
title('Fs=3.125 KHz')
xlabel('\omega (radians)')