If I take the fft of the impulse response of
x(t)=(1/(m*wd))*exp(-zeta*wn*t)*sin(wd*t)
I should get
G(jw)=1/(m*(-w^2+wn^2+2*zeta*wn*w*i))]
in the frequency domain. I have not been able to accomplish this using the following code. Anyone have any ideas?
clear
k=25;
m=1;
c=1;
wn=sqrt(k/m)...
