实验二用FFT进行谱分析.pptVIP

实验目的 （1）进一步加深DFT算法原理和基本性质的理解（因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质）。熟悉FFT程序结构及编程方法。 （2）熟悉应用FFT对确定信号进行谱分析方法，熟悉FFT算法原理和FFT子程序的应用。 （3）学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应有FFT。 实验原理及方法 xa(t)的频谱 取样 x（t） x（n）序列 t k T 1 2T 2 3T 3 4T 4 幅度量化信号 X(k)是X(ejω)在[0,2π]区间上的N个等间隔采样点 一个序列x(n)的离散时间傅里叶变换就是它的频谱函数。 实验内容及步骤(Matlab DFT) function y=mydft(x) % y=mydft(x) % x为给定时间序列 % y为x的离散傅立叶变换 N=length(x); % 输入序列的长度 n=0:N-1; k=n; % 确定时域位置序列n和频域位置序列k WN=exp(-j*2*pi/N); % 计算DFT所需的旋转因子 nk=n*k; WNnk=WN.^nk; % 构成旋转因子矩阵 Xk=x*WNnk;y=Xk; % 按DFT定义计算x的傅立叶变换 实验内容及步骤(Ma