data = read('judahlevine.txt',-1,2);
t = data(:,1);
tilt = data(:,2);

tmax = 50;
dt = tmax/length(t);
time = 0:dt:tmax-dt;

subplot(2,1,1)
plot2d(time,tilt)

fmax = 1/(2 * dt);
df = fmax/floor(length(time));
freqs = 0:df:fmax-df;

subplot(2,1,2)
periodogram = 10*log10(abs(fft(tilt,-1)));
psd = periodogram(1:length(freqs));
plot2d(freqs,psd)