百度文库 - 让每个人平等地提升自我
11
实验一 快速傅里叶变换之报告
一 、实验目的
1、在理论学习的基础上,通过本实验加深对快速傅立叶变换的理解;
2、熟悉并掌握按时间抽取FFT算法的程序;
3、了解应用FFT进行信号频谱分析过程中可能出现的问题,例如混淆、泄漏、
栅栏效应等,以便在实际中正确应用FFT。
二 实验内容
a) 信号频率F=50Hz,采样点数N=32,采样间隔T=
matlab程序代码为:
F=50; T=; N=32;
n=0:N-1;
t=n*T;
A=sin(2*pi*F*t);
figure;
Y = fft(A,N);
h = (abs(Y));
h=h/max(h(1:N));
for n=1:N;
string1=strcat('X(',num2str(n-1),
')=',num2str(h(n)));
disp(string1);
f=(n/T)/N;
end
stem([0:N-1]/N/T,h);
xlabel('ƵÂÊ/HZ');
ylabel('Õñ·ùX£¨ejw£©');
title('·ùƵÌØÐÔ');
上述代码命令中,将FFT变换后的数字变量K,在画图时转换成频域中的频率f。这主
要是根据数字频率与模拟域频率之间的关系:
T
其中
、
分别为数字和模拟域中的频率,且
N
k
2
f
2
于是有:
NT
k
f
运算结果:
X(1)=1 X(2)= X(3)= X(4)=