function xk=myfft(xk,N,n)
//fft  by swell2005
//其中xk為源信號(hào)(列向量），2^n=N
/////////////////////進(jìn)行倒序
i=1;l=1;
while l~=(N-1)
if l<i
	t=a(i);
	a(i)=a(l);
	a(l)=t;
end
k=N/2;
while k<i
	i=i-k;
	k=k/2;
end
i=i+k;
l=l+1;
end
////////////////////// DFT  逆變換？
xk=xk/N;
m=1;
while m<=n
	u=1;
    i=1;
    k=2^(m-1);
    w=exp(-j*pi/k);
	while k>=i
		l=i;
		while l<=2^n;
			t=xk(l+k)*u;
			xk(l+k)=xk(l)-t;
			xk(l)=xk(l)+t; 
			l=l+2^m;
		end
		u=u*w;
		i=i+1;
	end
	m=m+1;
end