求快速傅里叶变换后的函数表达式 ( 积分: 100 )

对采集的一组数据，求它的幅频特性的函数表达式和相频特性的函数表达式。

楼主有学过信号处理没？如果学过的话就看的懂下面的了。
下面是DIT-FIT算法：（在matlab里实现的）
function x=slowfft(x,N)
N=log2(N);
NV2=N/2;
NM1=N-1;
j=1;
for i=1:NM1
if i<j
T=x(j);x(j)=x(i);x(i)=T;
end
K=NV2；
while k<j
j=j-k;
k=k/2;
end
j=j+k;
end
for L=1:M
LE=2^L;
LE1=LE/2；
U=1；
w=exp(-j*pi/LE1);
for j=1:LE1
for i=j:LE:N
IP=I+LE1;
t=x(IP)*U;
x(IP)=x(i)-t;
x(i)=x(i)+t;
en
U=U*W;
end