Script_E.m

f1 = 852;
f2=1477;
fs = 8000;
n = 1600;
t = (0:(n-1))/fs;
f = (-n/2:n/2-1) * (fs/n);

x1 =@(t)sin(2*pi*f1*t);
% subplot(5,3,1);
% plot(t,x1(t));
% title('x1 =sin(2*pi*f1*t)');
% xlabel('Time');
% ylabel('Amplitude');
% sound(x1(t));

x2 = @(t)sin(2*pi*f2*t);
% subplot(5,3,2);
% plot(t,x2(t));
% xlabel('Time');
% ylabel('Amplitude');
% title('x2 =sin(2*pi*f2*t)');
%  sound(x2(t));
 
x =@(t)((sin(2*pi*f1*t)+sin(2*pi*f2*t)));
subplot(5,3,3);
plot(t,x(t));
xlabel('Time');
ylabel('Amplitude');
title('x =(sin(2*pi*f1*t) + sin(2*pi*f2*t))');
sound(x(t)); 

%Fourier Transform

 s1 = abs(fft(x1(t))); 
 title('Spectrum of x1 ');
 xlabel('frequency');
 ylabel('Amplitude');
 s1shift = fftshift(s1);
 subplot(5,3,4);
 plot(f,s1shift);
 

 s2 = abs(fft(x2)); 
 title('Spectrum of x2 ');
 xlabel('frequency');
 ylabel('Amplitude');
 s2shift = fftshift(s2);
 subplot(5,3,5);
 plot(f,s2shift);
 
 
 s3 = abs(fft(x)); 
 title('Spectrum of x ');
 xlabel('frequency');
 ylabel('Amplitude');
 s3shift = fftshift(s3);
 subplot(5,3,6);
 plot(f,s3shift);