MATLAB(离散傅里叶变换）

（1）fft(X):返回向量X的离散傅里叶变换。若X的长度为2的幂次则进行快速傅里叶变换，否则为普通较慢的算法

（2）fft(X,N) 计算N点离散傅里叶变换。

（3）fft(X,N,dim)

（4）ifft为傅里叶逆变换

例如：给定x(t)，取N=128,对t从0-1s采样，用fft函数作快速傅里叶变换，绘制振幅，频率图。

N=128;
T=1;
t=linspace(0,T,N);  %N个采样时间
x=12*sin(2*pi*10*t+pi/4)+5*cos(2*pi*40*t);
dt=t(2)-t(1);
f=1/dt;
X=fft(x);
F=X(1:N/2+1);
f=f*(0:N/2)/N;
plot(f,abs(F),'-*')
xlabel('Frequency');
ylabel('|F(k)|');

（可复制代码到matlab看绘图结果）

ix=real(ifft(X)); 求逆变换取实部

plot(t,x,t,ix,':') 逆变换与原函数的曲线