实验一FFT与谱答案.docx

《数字信号处理实验报告》姓名：朱道萍班级：通信1401学号：201403090129实验一 FFT与谱分析1.实验目的加深对DFT算法原理和物理意义的理解。(2)熟悉FFT算法，增强对FFT结果的分析能力。(3)掌握用FFT对连续时间信号和离散时间信号进行谱分析的方法，了解误差及其产生的原因。2.实验原理(1)复习 DFT 的定义、性质和用 DFT 作谱分析的有关内容。(2)复习 FFT 算法原理与编程思想。(3)查询Matlab中fft函数的用法，学会通过改变fft的参数计算不同DFT。(4)编制信号产生子程序，产生一下信号：编写实验程序： 参考图1-1所示的程序流程图，编写实验程序。3.实验内容 对 2 中所给出的11个信号逐一进行谱分析。各信号的FFT变换区间N、连续时间信号x6(t)的采样频率fs等参数的参考取值如下： a) x1[n],x2[n], x3[n], x4[n], x5[n] , x7[n], x8[n]: N=8, 16 b) x6(t): fs=64Hz，N=16, 32, 64 c) x9[n], x10[n], x11[n]: N=16 x1[n]:N=8n=0:N-1;x1=[1 1 1 1 0 0 0 0];m=0:2*N-1;x2=[1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0];f1=fft(x1,N);f2=fft(x2,2*N);subplot(2,2,1)stem(n,x1)subplot(2,2,2)stem(n,abs(f1))subplot(2,2,3)stem(m,x2)subplot(2,2,4)stem(m,abs(f2))x2[n]:N=8 n=0:N-1; x1=[1 2 3 4 4 3 2 1];m=0:2*N-1; x2=[1 2 3 4 4 3 2 1 0 0 0 0 0 0 0 0];f1=fft(x1,N);f2=fft(x2,2*N);subplot(2,2,1)stem(n,x1) subplot(2,2,2)stem(n,abs(f1))subplot(2,2,3)stem(m,x2)subplot(2,2,4)stem(m,abs(f2))x3[n]:N=8n=0:N-1; x1=[4 3 2 1 1 2 3 4];m=0:2*N-1; x2=[4 3 2 1 1 2 3 4 0 0 0 0 0 0 0 0];f1=fft(x1,N);f2=fft(x2,2*N)subplot(2,2,1)stem(n,x1)subplot(2,2,2)stem(n,abs(f1))subplot(2,2,3)stem(m,x2) subplot(2,2,4)stem(m,abs(f2))x4[n]:N=8;n=0:N-1; x1=cos(0.25*pi*n);m=0:2*N-1; x2=cos(0.25*pi*m); f1=fft(x1,N);f2=fft(x2,2*N); subplot(2,2,1)stem(n,x1) subplot(2,2,2)stem(n,abs(f1)) subplot(2,2,3)stem(m,x2) subplot(2,2,4)stem(m,abs(f2))x5[n]:N=8;n=0:N-1; x1=sin(0.125*pi*n);m=0:2*N-1; x2=sin(0.125*pi*m);f1=fft(x1,N);f2=fft(x2,2*N);subplot(2,2,1)stem(n,x1) subplot(2,2,2)stem(n,abs(f1))subplot(2,2,3)stem(m,x2)subplot(2,2,4)stem(m,abs(f2))x6[n]:N=16;fs=64;n=0:N-1; x1=cos(8*pi*n/fs)+cos(16*pi*n/fs)+cos(20*pi*n/fs);m=0:2*N-1; x2=cos(8*pi*m/fs)+cos(16*pi*m/fs)+cos(20*pi*m/fs);l=0:4*N-1; x3=cos(8*pi*l/fs)+cos(16*pi*l/fs)+cos(20*pi*l/fs);f1=fft(x1,N);f2=fft(x2,2*N);f3=fft(x3,4*N);subplot(3,2,1)stem(n,x1)subplot(3,2,2)stem(n,abs(f1))subplot(3,2,3)stem(m,x2)subplot(3,2,4)stem(m,abs(f2))subplot(3,2,5)stem(l,x3) subplot(3,2,6)stem(l,abs(f3))