close all;


N=8;
w=10;
m=1;
k=1;
c=1;
I=sqrt(-1);

T=2*pi/w;
delT=T/N;
t=0:delT:(T-delT);
sampledLoad = cos(w*(0:N-1)*delT);
[D,ws] = fft_easy(sampledLoad,t);
subplot(2,3,1);
stem([ws ; -ws(2:end)],abs([D ; conj(D(2:end))]),'o','fill');
subplot(2,3,2);
stem([ws ; -ws(2:end)],angle([D ; conj(D(2:end))]),'o','fill');


X=D ./ ( - ws.^2 * m  + I*c*ws + k);
subplot(2,3,3);
stem([ws ; -ws(2:end)],abs([X ; conj(X(2:end))]),'o','fill');
subplot(2,3,4);
stem([ws ; -ws(2:end)],angle([X ; conj(X(2:end))]),'o','fill');

xt=ifft_easy(X,ws);
subplot(2,3,5);
plot(t,xt);

s=tf('s');
sys=1/(m*s^2+s*c+k);
subplot(2,3,6);
bode(c2d(sys,delT));