Tp=64/1000;
Fs=1000;T=1/Fs;
N=Tp*Fs;
n=0:N-1;
A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;
xnt=A*exp(-alph*n*T).*sin(omega*n*T);
xk=T*fft(xnt,N);
wk=2*pi*n/N;
fk=n/Tp;
subplot(2,1,1);
stem(fk,xnt);
k=0:N-1;fk=k/Tp;
subplot(2,1,2);
plot(fk,abs(xk))
Fs=1000;T=1/Fs;
N=Tp*Fs;
n=0:N-1;
A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;
xnt=A*exp(-alph*n*T).*sin(omega*n*T);
xk=T*fft(xnt,N);
wk=2*pi*n/N;
fk=n/Tp;
subplot(2,1,1);
stem(fk,xnt);
k=0:N-1;fk=k/Tp;
subplot(2,1,2);
plot(fk,abs(xk))
