数字信号处理第三四章作业.pdf

成 绩： 《数字信号处理》 作业与上机实验 （第三四章） 班 级： 13-电信 学 号： 姓 名： 任课老师： 李宏民 完成时间： 15.11.08 信息与通信工程学院 2015—2016 学年 第 1学期 第3-4章 离散傅里叶变换及快速算法与应用 一、实验内容及要求： 1、对信号在x(n) {1,2,3,4,5,6,7,8}，n 0，1，27,求其N 8 点的DFT。要求采用基于时间抽取算法编写FFT实现程序，画出DFT 幅度谱与相位谱,并将计算结果与用MATLAB 自带的FFT函数计算结果 进行比较。 2、一个由40Hz和100 Hz正弦信号构成的信号，受零均值随机噪声 的干扰（噪声服从标准正态分布，由randn函数产生，n(t) randn(m,n) 返回一个m*n的随机矩阵），即s(t) sin(2100t)sin(240t)n(t) 。数据 采样率为500Hz，试用FFT函数来分析其信号频率成分。 ① 求其幅度频谱，从频谱图 （横坐标以HZ为单位，用plot函数画 图）中能否观察出信号的2个频率分量？ ② 提高采样点数，再求该信号的幅度频谱图，此时幅度频谱发生了 什么变化？信号的2个模拟频率和数字频率各为多少？FFT频谱分析 结果与理论上是否一致？ 3、研究高密度频谱与高分辨率频谱。频率分辨率是指所用的算法能 将信号中两个靠得很近的谱峰分开的能力。信号末尾补零由于没有对 原信号增加任何新的信息，因此不能提高频率分辨率，但可以减小栅 栏效应，所得到的频谱称为高密度频谱。在维持采样频率不变的情况 下，为提高分辨率只能增加采样点数N，此时所得到的频谱称为高分 辨率频谱。设有连续信号 3 3 3 x (t) cos(26.510 t) cos(2710 t) cos(2910 t) a f 32kHz x (t) s a 以采样频率 对信号 采样，分析下列几种情况的幅频特 性。 1 采集数据长度N 16点，做N 16的FFT；采集数据长度N 16点， 补零到256点，做256点的FFT； 2 采集数据长度N 64点，做N 64 的FFT；采集数据长度N 64点， 补零到256点，做256点的FFT； 3 采集数据长度N 256点，做N 16点的FFT。 4 观察以上几幅不同的幅频特性曲线，分析和比较它们的特点，并 说明形成的原因。 5 注：在MATLAB 中，可用zeros函数实现填零运算，例如 x1 [x(1:1:16),zeros(1,240)]。所有频谱图的横坐标以HZ为单 位，用plot函数画图。 6 一个连续信号含三个频率分量，经采样得以下序列 x(n) sin2 0.15n cos2 (0.15 f )n cos2 (0.15 2f )n n 0,1,N  1。已知N 16， 分别为1/16，1/64，观察其频谱； f f 当N 64，128， 不变，其结果有何不同，为什么？ 4、语音信号频谱分析 首先应用录音软件(见下图)，录音 “我是XXX （学生姓名）”，得 到语音数据 “录音.WAV”并将 “录音.WAV”文件复制到 MATLABWORK 目录下。其次在M