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

成 绩： 《数字信号处理》 作业与上机实验 （第三四章） 班 级： 13- 电信 学 号： 姓 名： 任课老师： 李宏民 完成时间： 15.11.08 信息与通信工程学院 2015— 2016 学年 第 1 学期 第 3-4 章 离散傅里叶变换及快速算法与应用 一、实验内容及要求： 1、对信号在 x(n)={1,2,3,4,5,6,7,8} ，n=0， 1，2....7, 求其 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，此时所得到的频谱称为高分辨率频谱。设有连续信号 xa (t ) cos(2 6.5 103 t) cos(2 7 103 t ) cos(2 9 103 t ) 以采样频率 f s 32kHz 对信号 xa (t) 采样，分析下列几种情况的幅频特性。 ① 采集数据长度 N=16点，做 N=16的 FFT；采集数据长度 N=16点， 补零到 256 点，做 256 点的 FFT； ② 采集数据长度 N=64点，做 N=64的 FFT；采集数据长度 N=64点， 补零到 256 点，做 256 点的 FFT； ③ 采集数据长度 N=256点，做 N=16点的 FFT。 ④ 观察以上几幅不同的幅频特性曲线，分析和比较它们的特点，并 说明形成的原因。 ⑤ 注：在 MATLAB中，可用 zeros 函数实现填零运算，例如 x1=[x(1:1:16),zeros(1,240)] 。所有频谱图的横坐标以 HZ为单 位，用 plot 函数画图。 ⑥ 一个连续信号含三个频率分量，经采样得以下序列 x(n)  sin 2  0.15n  cos2  (0.15  f )n  cos2  (0.15  2 f )n 0,1,....... N 1。已知 N =16， f 分别为 1/16 ，1/64 ，观察其频谱；当 N =64，128， f 不变，其结果有何不同，为什么？ 4、语音信号频谱分析 首先应用录音软件 ( 见下图 ) ，录音“我是 XXX（学生姓名）”，得 到语音数据“录音 .WAV”并将“录音 .WAV”文件复制到 MATLABWORK 目录下。其次在 MATLAB环境下， X(n)=wavread( 录音 .WAV,2048) ， 得到 2048 点的语音数据 X(n), 对 X(n) 进行 FFT 频谱分析，画出频谱 图（横坐标以 HZ为单位，用 plot 函数画图）。 5、双音多频（ DTMF）信号是将拨号盘上的 0-F 共 16 个数字，用音频范围的 8 个频率来表示的一种编码方式。 8 个频率分为高频群和低频群两组，分别作为列频和行频。 每个字符的信号由来自列频和行频的两个频率的正弦信号叠加而成。频率组合方式如表 8-2 所示。 双音多频（ DTMF）信号频率组合方式 频率 1209Hz 1336Hz 1477Hz 1633Hz 697Hz 1 2 3 A 770Hz 4 5 6 B 852Hz 7 8 9 C 941Hz  */E  0  #/F  D 找一部电话机，在免提状态下拨号  0～9  中的任意  7 位号码，并录制 双音多频拨号声音。采用快速傅里叶变换（ FFT）分析这 7 个号码拨 号时的频谱，并与理论值进行比较。 图 1 x(n ) 波形 三、波形图 1．DFT幅度谱与相位谱如图 1 可见 xn=[1,2,3,4,5,6,7,8]; M=nextpow2(length(xn)), N=2^M,