电子信息工程数字号实验报告.docx

数字信号处理实验报告实验一：用 FFT 做谱分析实验目的 1、进一步加深 DFT 算法原理和基本性质的理解。 2、熟悉 FFT 算法原理和 FFT 子程序的应用。 3、学习用FFT 对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用 FFT。二、实验原理用FFT对信号作频谱分析是学习数字信号处理的重要内容。经常需要进行谱分析的信号是模拟信号和时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D和分析误差。频谱分辨率直接和FFT的变换区间N有关，因为FFT能够实现的频率分辨率是2π/N≤D。可以根据此时选择FFT的变换区间N。误差主要来自于用FFT作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N较大时离散谱的包络才能逼近于连续谱，因此N要适当选择大一些。周期信号的频谱是离散谱，只有用整数倍周期的长度作FFT，得到的离散谱才能代表周期信号的频谱。如果不知道信号周期，可以尽量选择信号的观察时间长一些。对模拟信号的频谱时，首先要按照采样定理将其变成时域离散信号。如果是模拟周期信号，也应该选取整数倍周期的长度，经过采样后形成周期序列，按照周期序列的谱分析进行。三、实验内容和步骤对以下典型信号进行谱分析： 对于以上信号，x1(n)~x5(n) 选择FFT的变换区间N为8和16 两种情况进行频谱分析。分别打印其幅频特性曲线。 并进行对比、分析和讨论;；x6(t)为模拟周期信号，选择 采样频率，变换区间N=16,32,64 三种情况进行谱分析。分别打印其幅频特性，并进行分析和讨论。四、实验代码如下所示：clear; %清除变量 close all; %关闭全部绘图窗口 b=menu(请选择信号x1(n)--x6(n),x1(n)=R4(n),x2(n)=[1 2 3 4 4 3 2 1],... x3(n)=[4 3 2 1 1 2 3 4],x4(n)=cos(npi/4),x5(n)=sin(npi/8),... x6(n)=cos(8pit)+cos(16pit)+cos(20pit),Exit); i=0; A=[8,16,32,64]; while(b~=7) %当选择EXIT时，返回值7，则退出循环 if b==6 m=menu(请选择FFT变换区间长度N,N=16,N=32,N=64); N=A(m+1); fs=64; n=0:(N-1); x=cos(8*pi*n/fs)+cos(16*pi*n/fs)+cos(20*pi*n/fs); else m=menu(请选择FFT变换区间长度N,N=8,N=16,N=32); N=A(m); n=0:(N-1); if b==1 x=[1,1,1,1,0,0,0,0,zeros(1,N-8)]; elseif b==2 x=[1,2,3,4,4,3,2,1,zeros(1,N-8)]; elseif b==3 x=[4,3,2,1,1,2,3,4,zeros(1,N-8)]; elseif b==4 x=cos(n*pi/4); elseif b==5 x=sin(n*pi/8); end end %先画出信号源图 i=i+1; figure(i); %创建绘图窗口 subplot(2,2,1); %指定1号子图 xlabel(n); %标记X坐标 stem(n,x,.r); ylabe