用FFT对信号作频谱分析实验报告记录.docVIP

用FFT对信号作频谱分析实验报告记录 PAGE PAGE 1 ———————————————————————————————— 作者： ———————————————————————————————— 日期： 实验一报告、用FFT对信号作频谱分析 一、实验目的 学习用FFT对连续信号和时域离散信号进行频谱分析的方法，了解可能出现的分析误差及其原因，以便正确应用FFT。 二、实验内容 1．对以下序列进行频谱分析： 选择FFT的变换区间N为8和16两种情况进行频谱分析。分别打印其幅频特性曲线，并进行对比，分析和讨论。 2．对以下周期序列进行频谱分析： 选择FFT的变换区间N为8和16两种情况分别对以上序列进行频谱分析。分别打印其幅频特性曲线，并进行对比、分析和讨论。 3.对模拟信号进行频谱分析： 选择采样频率，对变换区间N=16，32，64 三种情况进行频谱分析。分别打印其幅频特性，并进行分析和讨论。 三、实验程序 1.对非周期序列进行频谱分析代码： close all;clear all; x1n=[ones(1,4)]; M=8;xa=1:(M/2);xb=(M/2):-1:1;x2n=[xa,xb]; x3n=[xb,xa]; X1k8=fft(x1n,8);X1k16=fft(x1n,16); X2k8=fft(x2n,8);X2k16=fft(x2n,16); X3k8=fft(x3n,8);X3k16=fft(x3n,16); subplot(3,2,1);mstem=(X1k8);title((1a)8点DFT[x_1(n)]); subplot(3,2,2);mstem=(X1k16);title((1b)16点DFT[x_1(n)]); subplot(3,2,3);mstem=(X2k8);title((2a)8点DFT[x_2(n)]); subplot(3,2,4);mstem=(X2k16);title((2b)16点DFT[x_2(n)]); subplot(3,2,5);mstem=(X3k8);title((3a)8点DFT[x_3(n)]); subplot(3,2,6);mstem=(X3k16);title((3b)16点DFT[x_3(n)]); 2.对周期序列进行频谱分析代码： N=8;n=0:N-1; x4n=cos(pi*n/4); x5n=cos(pi*n/4)+cos(pi*n/8); X4k8=fft(x4n); X5k8=fft(x5n); N=16;n=0:N-1; x4n=cos(pi*n/4); x5n=cos(pi*n/4)+cos(pi*n/8); X4k16=fft(x4n); X5k16=fft(x5n); figure(2) subplot(2,2,1);mstem(X4k8);title((4a)8点 DFT[x_4(n)]); subplot(2,2,2);mstem(X4k16);title((4b)16点DFT[x_4(n)]); subplot(2,2,3);mstem(X5k8);title((5a)8点DFT[x_5(n)]); subplot(2,2,4);mstem(X5k16);title((5a)16点DFT[x_5(n)]) 3.模拟周期信号谱分析 figure(3) Fs=64;T=1/Fs; N=16;n=0:N-1; x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); X6k16=fft(x6nT); X6k16=fftshift(X6k16); Tp=N*T;F=1/Tp; k=-N/2:N/2-1;fk=k*F; subplot(3,1,1);stem(fk,abs(X6k16),.);box on title((6a)16μ?DFT[x_6(nT)]);xlabel(f(Hz));ylabel(·ù?è); axis([-N*F/2-1,N*F/2-1,0,1.2*max(abs(X6k16))]); N=32;n=0:N-1; %FFTμ?±???????N=32 x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); X6k32=fft(x6nT); X6k32=fftsh