哈工大数字信号处理实验报告..docVIP

Harbin Institute of Technology 实验报告 课程名称： 数字信号处理 院 系： 电子与信息工程学院 班 级： 姓 名： 学 号： 指导教师： 冀振元 实验时间： 2014 年 11 月 哈尔滨工业大学 实验一： 用FFT作谱分析 实验目的： 1．进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一 种快速算法，所以FFT的运算结果必然满足DFT的基本性质)。 2．熟悉FFT算法原理和FFT子程序的应用。 3．学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出 现的分析误差及其原因，以便在实际中正确应用FFT。 实验步骤 (1) 复习DFT的定义、 性质和用DFT作谱分析的有关内容。 (2) 复习FFT算法原理与编程思想， 并对照DIT-FFT运算流图和程序框图， 读懂本实验提供的FFT子程序。 (3) 编制信号产生子程序， 产生以下典型信号供谱分析用： (4) 编写主程序。 (5) 按实验内容要求， 上机实验， 并写出实验报告。 实验程序与结果: 实验程序： %整个程序通过switch语句切换x1 ~ x6 的波形产生和FFT变换，波形。clc; clear; n0=-10:0.1:20; n1=1:4; x1=[1 1 1 1]; n2=0:7; x2=(n2=0 n2=3).*(n2+1)+(n2=4 n2=7).*(8-n2); n3=0:7; x3=(n3=0 n3=3).*(4-n3) + (n3=4 n3=7).*(n3-3); n4=0:19; x4=cos(pi/4*n4); n5=0:19; x5=sin(pi/8*n5); switch v case 1 subplot(311); stem(n1,x1);xlabel(n); ylabel(x1(n)); title(原时间序列); N=8; n1=0:N-1; X=fft(x1,N);X=abs(X); subplot(312); stem(n1,fftshift(X)); title(N=8); N=16; n1=0:N-1; X=fft(x1,N);X=abs(X); subplot(313); stem(n1,fftshift(X)); title(N=16); case 2 subplot(311); stem(n2,x2);xlabel(n); ylabel(x2(n)); title(原时间序列); N=8; n2=0:N-1; X=fft(x2,N);X=abs(X); subplot(312); stem(n2,fftshift(X)); title(N=8); N=16; n2=0:N-1; X=fft(x2,N);X=abs(X); subplot(313); stem(n2,fftshift(X)); title(N=16); case 3 subplot(311); stem(n3,x3);xlabel(n); ylabel(x3(n)); title(原时间序列); N=8; n3=0:N-1; X=fft(x3,N);X=abs(X); subplot(312); stem(n3,fftshift(X)); title(N=8); N=16; n3=0:N-1; X=fft(x3,N);X=abs(X); subplot(313); stem(n3,fftshift(X)); title(N=16); case 4 subplot(311); stem(n4,x4);xlabel(n); ylabel(x4(n)); title(原时间序列);