[实验三用FFT对信号作频谱分析.docxVIP

实验三用FFT对信号作频谱分析一：实验目的学习用FFT对连续信号和时域离散信号进行频谱分析的方法，了解可能出现的分析误差及其原因，以便正确使用FFT.二：实验原理用FFT对信号作频谱分析是学习数字信号处理的重要内容。经常需要进行谱分析的信号是模拟信号和时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D和分析误差。频谱分辨率直接和FFT的变换区间N有关，因为FFT能够实现的频率分辨率是，因此要求。可以根据此式选择FFT的变换区间N。误差主要来自于用FFT作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N较大时离散谱的包络才能逼近于连续谱，因此N要适当选择大一些。周期信号的频谱是离散谱，只有用整数倍周期的长度作FFT，得到的离散谱才能代表周期信号的频谱。如果不知道信号周期，可以尽量选择信号的观察时间长一些。对模拟信号进行谱分析时，首先要按照采样定理将其变成时域离散信号。如果是模拟周期信号，也应该选取整数倍周期的长度，经过采样后形成周期序列，按照周期序列的谱分析进行。三：实验步骤及步骤（1）对以下序列进行谱分析。　　　　　　 选择FFT的变换区间N为8和16 两种情况进行频谱分析。分别打印其幅频特性曲线。 并进行对比、分析和讨论。（2）对以下周期序列进行谱分析。选择FFT的变换区间N为8和16 两种情况分别对以上序列进行频谱分析。分别打印其幅频特性曲线。并进行对比、分析和讨论。（3）对模拟周期信号进行谱分析选择 采样频率，变换区间N=16,32,64 三种情况进行谱分析。分别打印其频特性，并进行分析和讨论。四：试验程序及结果分析（1）对下面序列进行频谱分析X1(n)的程序代码：x1=[1,1,1,1];X16k=fft(x1,16);X8k=X16k(1:2:16);N1=1:16;N2=1:8;subplot(2,1,1);stem(2/16*N1, abs(X16k),.);title(x1(n)16点FFT);xlabel(w/pi)subplot(2,1,2);stem(2/8*N2, abs(X8k),.);title(x1(n)8点FFT);xlabel(w/pi)运行结果：X2(n)的程序代码：X2=[1,2,3,4,4,3,2,1];X16k=fft(X2,16);X8k=X16k(1:2:16);N1=1:16;N2=1:8;subplot(2,1,1);stem(2/16*N1, abs(X16k),.);title(x2(n)16点FFT);xlabel(w/pi)subplot(2,1,2);stem(2/8*N2, abs(X8k),.);title(x2(n)8FFT);xlabel(w/pi)运行结果：X3(N)的程序代码：x3=[4,3,2,1,1,2,3,4];X16k=fft(x3,16);X8k=X16k(1:2:16);N1=1:16;N2=1:8;subplot(2,1,1);stem(2/16*N1, abs(X16k),.);title(x3(n)的16点 FFT);xlabel(w/pi)subplot(2,1,2);stem(2/8*N2, abs(X8k),.);title(x3(n) 的8点FFT);xlabel(w/pi)运行结果：结果分析:、和 的8点DFT和16点DFT分别是、和的频谱函数的8点和16点采样；（2）周期序列进行谱分析X4(n)的程序代码：n=0:15;x4=cos(pi*n/4);X16k=fft(x4,16);n=0:7;x4=cos(pi*n/4);X8k=fft(x4,8);N2=0:7;N1=0:15;subplot(2,1,1);stem(2/16*N1, abs(X16k),.);title(x4(n)的16点 FFT)；xlabel(w/pi)subplot(2,1,2);stem(2/8*N2, abs(X8k),.);title(x4(n) 的8点FFT)；xlabel(w/pi)运行结果：X5(n)的程序代码：n=0:15;x5=cos(pi*n/4)+cos(pi*n/8);X16k=fft(x5,16);n=0:7;x5=cos(pi*n/4)+cos(pi*n/8);X8k=fft(x5,8);N2=0:7;N1=0:15;subplot(2,1,1);stem(2/16*N1, abs(X16k),.);title(x5(n)μ?16μ? FFT);xlabel(w/pi)subplot(2,1,2);stem(2/8*N2, abs(X8k),.);title(x5(n) μ?8μ?FFT);xlabel(w/pi)运行结果：实验结果分析：的周期为8，N=8和N=16均是其周期的整数倍，得到正确的单一频率