资源描述:
《用FFT做谱分析实验报告》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、........实验二用FFT做谱分析一、实验目的1.进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法,所以FFT的运算结果必然满足DFT的基本性质)。2.熟悉FFT算法原理和FFT子程序的应用。3.学习用FFT对连续信号和时域离散信号进行谱分析的方法二、实验原理如果给出的是连续信号xa(t),则首先要根据其最高频率确定抽样频率fs以及由频率分辨率选择抽样点数N,然后对其进行软件抽样(即计算x(n)=xa(nT),0≤n≤N-1),产生对应序列x(n)。再利用MATLAB所提供的库函数fft(n,x)进行FFT计算三、
2、实验内容①实验信号:x1(n)=R4(n)x2(n)=x3(n)=x4(n)=cos(πn/4)x5(n)=sin(πn/8)x6(t)=cos8πt+cos16πt+cos20π.专业学习资料.........FFT变换区间及x6(t)抽样频率fsx1(n),x2(n),x3(n),x4(n),x5(n):N=8,16x6(t):fs=64(Hz),N=16,32,64②MATLAB程序代码N1=8;N2=16;x1=ones(1,4);x2=[1:4,4:-1:1];x3=[4:-1:1,1:4];n=0:1:16;x4=cos(pi*n/4
3、);x5=sin(pi*n/8);X11=fft(x1,N1);X11=abs(X11);X21=fft(x2,N1);X21=abs(X21);X31=fft(x3,N1);X31=abs(X31);X41=fft(x4,N1);X41=abs(X41);X51=fft(x5,N1);.专业学习资料.........X51=abs(X51);X12=fft(x1,N2);X12=abs(X12);X22=fft(x2,N2);X22=abs(X22);X32=fft(x3,N2);X32=abs(X32);X42=fft(x4,N2);X42=
4、abs(X42);X52=fft(x5,N2);X52=abs(X52);figure(1);subplot(3,1,1);stem(x1);grid;%x1时域波形xlabel('n');ylabel('x1(n)')title('N=8的时域图')subplot(3,1,2);stem(X11);grid;%x1在N=8的FFT变换频谱图xlabel('Hz');ylabel('
5、X11(k)
6、')title('N=8的频谱图')subplot(3,1,3);stem(X12);grid;%x1在N=16的FFT变换频谱图xlabel('Hz
7、');ylabel('
8、X12(k)
9、')title('N=16的频谱图')figure(2);.专业学习资料.........subplot(3,1,1);stem(x2);grid;%x2时域波形xlabel('n');ylabel('x2(n)')title('N=8的时域图')subplot(3,1,2);stem(X21);grid;%x2在N=8的FFT变换频谱图xlabel('Hz');ylabel('
10、X21(k)
11、')title('N=8的频谱图')subplot(3,1,3);stem(X22);grid;%x2在N=16的F
12、FT变换频谱图xlabel('Hz');ylabel('
13、X22(k)
14、')title('N=16的频谱图')figure(3);subplot(3,1,1);stem(x3);grid;%x3时域波形xlabel('n');ylabel('x3(n)')title('N=8的时域图')subplot(3,1,2);stem(X31);grid;%x3在N=8的FFT变换频谱图xlabel('Hz');ylabel('
15、X31(k)
16、')title('N=8的频谱图')subplot(3,1,3);stem(X32);grid;%x3在N=16的
17、FFT变换频谱图xlabel('Hz');ylabel('
18、X32(k)
19、')title('N=16的频谱图')figure(4);subplot(3,1,1);stem(x4);grid;%x4时域波形xlabel('n');ylabel('x4(n)').专业学习资料.........title('N=8的时域图')subplot(3,1,2);stem(X41);grid;%x4在N=8的FFT变换频谱图xlabel('Hz');ylabel('
20、X41(k)
21、')title('N=8的频谱图')subplot(3,1,3);stem(X4
22、2);grid;%x4在N=16的FFT变换频谱图xlabel('Hz');ylabel('
23、X42(k)
24、')title('N=16的