% Script:  aliasing_demo_audible

% Demonstration Script of Audible Aliasing
% Using Sinusoidal Signals
% Original: Professor Marple
% Modified: Samson Cheung

echo on
n = 0:15000; % time interval for samples >>> 3 SECONDS OF TONE <<<<<<
Fs = 5000;   % sampling rate in Hertz
Ts = 1/Fs;   % sampling interval in seconds (0.2 ms)
F1 = 500;    % sinusoidal frequency in Hertz <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
x1 = cos(2*pi*F1*n*Ts);
disp('Hit space bar to play sound')
pause
sound(x1,Fs)
pause
disp(' ')

n = 0:15000; % time interval for samples
Fs = 5000;   % sampling rate in Hertz
Ts = 1/Fs;   % sampling interval in seconds
F2 = 2000;   % sinusoidal frequency in Hertz <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
x2 = cos(2*pi*F2*n*Ts);
disp('Hit space bar to play sound')
pause
sound(x2,Fs)
pause
disp(' ')

n = 0:15000; % time interval for samples
Fs = 5000;   % sampling rate in Hertz
Ts = 1/Fs;   % sampling interval in seconds
F3 = 3000;   % sinusoidal frequency in Hertz <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
x3 = cos(2*pi*F3*n*Ts);
disp('Hit space bar to play sound')
pause
sound(x3,Fs)
pause
disp(' ')

n = 0:15000; % time interval for samples
Fs = 5000;   % sampling rate in Hertz
Ts = 1/Fs;   % sampling interval in seconds
F4 = 4500;   % sinusoidal frequency in Hertz <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
x4 = cos(2*pi*F4*n*Ts);
disp('Hit space bar to play sound')
pause
sound(x4,Fs)
pause
disp(' ')

n = 0:15000; % time interval for samples
Fs = 5000;   % sampling rate in Hertz
Ts = 1/Fs;   % sampling interval in seconds
F5 = 5500;   % sinusoidal frequency in Hertz <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
x5 = cos(2*pi*F5*n*Ts);
disp('Hit space bar to play sound')
pause
sound(x5,Fs)
pause
disp(' ')

n = 0:15000; % time interval for samples
Fs = 5000;   % sampling rate in Hertz
Ts = 1/Fs;   % sampling interval in seconds
F6 = 7000;   % sinusoidal frequency in Hertz <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
x6 = cos(2*pi*F6*n*Ts);
disp('Hit space bar to play sound')
pause
sound(x6,Fs)
echo off