%   Jerome P. Lynch
%   CEE619 - Advanced Dynamics and Smart Structures
%   Question #4 - Homework #1

%   Generate Noise and Sine Signal
%   ==============================
%clear;
%for i=1:4000
%    x(i,1) = ((30000 + round(randn(1)*10))/(2^16-1))*5;
%    t(i,1) = i*0.01;
%    y(i,1) = 0.0004*sin(25.6637*t(i,1));
%    tot(i,1) = x(i,1) + y(i,1);
%end;

%a = tot;
%save RoofAccel t a

%   Load Signal
%   ===========
load RoofAccel
subplot(2,1,1)
plot(t,a);
subplot(2,1,2)
autocorr(a,150);