实验二___用FFT作谱分析.docVIP

华中师范大学信息技术系课程实验报告 姓名： 日期： 2012 年 11 月 18 日 课程名称： 数字信号处理 指导教师： 实验题目： 实验二 用FFT作谱分析 一、实验目的 （1）进一步加深DFT算法原理和基本性质的理解； （2）熟悉FFT算法原理和FFT子程序的应用； （3）学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。 二、实验仪器 计算机和MATLAB软件 三、实验内容、步骤、分析 1．对下列信号作谱分析： （1），； %N=8 n=0:7; x1=ones(1,4); Xk1=abs(fft(x1,8)); figure(1) subplot(2,1,1) stem(n,[x1 0 0 0 0]);grid on; subplot(2,1,2) stem(n,Xk1);grid on; %N=16 n=0:15; x1=ones(1,4); Xk1=abs(fft(x1,16)); figure(2) subplot(2,1,1) stem(n,[x1 zeros(1,12)]);grid on; subplot(2,1,2) stem(n,Xk1);grid on; （2）， x=[1,2,3,4,4,3,2,1]; subplot(2,1,1); y=abs(fft(x,8)); stem(y);grid on; subplot(2,1,2); y=abs(fft(x,16)); stem(y);grid on; x=[4,3,2,1,1,2,3,4]; subplot(2,1,1); y=abs(fft(x,8)); stem(y);grid on; subplot(2,1,2); y=abs(fft(x,16)); stem(y);grid on; （3）， x=[4,3,2,1,1,2,3,4]; subplot(2,1,1); y=abs(fft(x,8)); stem(y);grid on; subplot(2,1,2); y=abs(fft(x,16)); stem(y);grid on; x=[1,2,3,4,4,3,2,1]; subplot(2,1,1); y=abs(fft(x,8)); stem(y);grid on; subplot(2,1,2); y=abs(fft(x,16)); stem(y);grid on; （4）， %N=8 n=0:1:7; x=cos(0.25*pi*n); y=abs(fft(x,8)); subplot(4,1,1); stem(x);grid on; subplot(4,1,2); stem(y);grid on; %N=16 n=0:1:15; x=cos(0.25*pi*n); y=abs(fft(x,16)); subplot(4,1,3); stem(x);grid on; subplot(4,1,4); stem(y);grid on; ， %N=8 n=0:1:7; x=sin(0.125*pi*n); y=abs(fft(x,8)); subplot(4,1,1); stem(x);grid on; subplot(4,1,2); stem(y);grid on; %N=16 n=0:1:15; x=sin(0.125*pi*n); y=abs(fft(x,16)); subplot(4,1,3); stem(x);grid on; subplot(4,1,4); stem(y);grid on; ，， N=16,32,64 m=input(FFT点数=) n=0:1:(m-1) subplot(411) x6=cos(pi*n/8)+cos(pi*n/4)+cos(pi*n*5/16) stem(n,x6,.) subplot(412) xa=fft(x6,16) i=0:15 stem(i,abs(xa),.) subplot(413) xb=fft(x6,32) i=0:31 stem(i,abs(xb),.) subplot(414) xc=fft(x6,64) i=0:63 stem(i,abs(xc),.) FFT点数=16 FFT点数=32 FFT点数=64 （程序和结果分析） 令，用FFT计算8点和16点离散傅立叶变换，，并根据DFT的对称性，又求出和，并与（1）中所得结果比较。 %N=8 n=0:1:7 x1=cos(pi*n/4)+sin(pi*n/8) xk1=fft(x1,8)