fft算法的matlab实现过程详解
(一) 实验目的: 理解离散傅立叶变换时信号分析与处理的一种重要变换,特别是FFT 在数字信号处理中的高效率应用。
(二) 实验原理:
1、有限长序列x(n)的DFT的概念和公式:
2、FFT 算法
调用格式是
X= fft(x) 或X=fft(x.N)
对前者,若x 的长度是2的整数次幂,则按该长度实现x 的快速变换,否则,实现的是慢速的非2的整数次幂的变换; 对后者,N应为2的整数次幂,若x的长度小于N,则补零,若超过N,则舍弃N 以后的数据。Ifft 的调用格式与之相同。
(三) 实验内容
l、题- ~:若x(n)=cos(n*pi/6)是- 一个N=l2 的有限序列,利用MATLAB 计算它的
DFT 并画出图形。
源程序:
clc;
N=12;
n=0:N- 1;
k=0:N- 1;
xn=cos( n *pi/6);
W=exp(-j* 2* pi/N );
kn=n‘*k
Xk=xn*(W.^kn)
stem( n,Xk );
xlabel(k<