N=4096;
A=1;
theta0=pi/6;
f0=10;
fs=64;
n=0:N-1;
xn=A*exp(j*(2*pi*f0*n/fs+theta0));
k=641.5;
ek=exp(-j*2*pi*k*n/N);
Xk=sum(x.*ek)
N=4096;
A=1;
theta0=pi/6;
f0=10;
fs=64;
n=0:N-1;
xn=A*exp(j*(2*pi*f0*n/fs+theta0));
k=641.5;
ek=exp(-j*2*pi*k*n/N);
Xk=sum(x.*ek)