%
% Written by:
% --
% John L. Weatherwax                2009-04-21
%
% email: wax@alum.mit.edu
%
% Please send comments and especially bug reports to the
% above email address.
%
%-----

close all; clc; clear; 

% EPage 137

nt = 512;
T = 1;
dt = T/nt;
t = 0:dt:(T-dt);
df = 1/T;
fmax = nt/(2*T);
f = 0:df:df*(nt/2-1);

f1 = 30;
f2 = 400;

y = sin( 2*pi*f1*t ) + 2 * sin( 2*pi*f2*t );

Y = fft(y);
Ys = zeros(1,nt/2);
Ys(1:nt/2) = (2/nt)*Y(1:nt/2);

plot( f, abs(Ys));
xlabel('frequency Hz'); ylabel('abs(DFT)'); title('DFT');
grid('on');

%saveas( gcf, '../../WriteUp/Graphics/Chapter7/c_7_p_14_plot.eps', 'epsc' );