DSP实验Matlab分析信号.docx

DSP 实验报告 1 实 验 要 求 1 .1 常用数字信号序列的产生： 熟 悉 Ma t l a b 产 生 数 字 信 号 的 基 本 命 令 ， 加 深 对 数 字 信 号 概 念 的 理 解 ， 并 能 够 用 Ma t l a b 产 生 和 绘 制 出 一 些 常 用 离 散 信 号 序 列 ，例 如 δ ( n) 、单 位 阶 跃 序 列 u( n) 、矩 形 序 列 R( n) 、 正 弦 序 列 S in ( n w ) 等 。 1 .2 数字信号的基本运算： 加 、 减 、 尺 度 （ 乘 除 ） 和 移 位 是 数 字 信 号 处 理 中 最 基 本 的 算 术 运 算 ， 将 上 述 基 本 序 列 进 行 这 些 基 本 运 算 ， 得 到 多 个 序 列 构 成 的 组 合 序 列 。 通 过 本 次 实 验 ，掌 握 Ma t l a b 中 这 些 基 本 运 算 命 令 ，对 数 字 信 号 处 理 有 一 个 基 本 概 念 ， 为 后 面 的 数 字 信 号 分 析 和 滤 波 打 下 基 础 。 数 字 信 号 处 理 的 一 个 重 要 分 支 就 是 信 号 分 析 ， 而 信 号 分 析 的 基 本 工 具 是 离 散 傅 立 叶 变 换 。利 用 傅 立 叶 变 换 和 级 数 所 形 成 的 频 谱 分 析 技 术 作 为 处 理 连 续 信 号 的 重 要 工 具 已 经 应 用 得 很 久 了 , 1956 年 库 力 （ Cool e y ） 和 图 基 （ Tuke y ） 所 发 展 的 近 似 频 谱 的 快 速 算 法 为 频 谱 分 析 的 数 字 信 号 的 谱 分 析 铺 平 了 道 路 。 因 此 ， DFT（ FFT） 得 到 广 泛 应 用 。 1 .3 假设信号 x( n ) 由下述信号组成： x ( n) 0 . 1 * c os( 0 . 4 5 n ) s i n( 0 . 3 n ) c os( 0 . 3 1n ) 4 这 个 信 号 有 两 根 主 谱 线 0 .3 pi 和 0 .31 pi 靠 的 非 常 近 ， 而 另 一 根 谱 线 0 .4 5 pi 的 幅 度 很 小 ， 请 选 择 合 适 的 长 度 N 和 窗 函 数 ， 用 DTFT 和 DFT 分 析 其 频 谱 ， 得 到 清 楚 的 三 根 谱 线 。 通 过 本 次 实 验 ， 应 该 掌 握 ： (a ) 用 傅 立 叶 变 换 进 行 信 号 分 析 时 基 本 参 数 的 选 择 。 (b) 经 过 离 散 时 间 傅 立 叶 变 换（ DTFT） 和 有 限 长 度 离 散 傅 立 叶 变 换（ DFT） 后 信 号 频 谱 上 的 区 别 ， 前 者 DTFT 时 间 域 是 离 散 信 号 ， 频 率 域 还 是 连 续 的 ， 而 DFT 在 两 个 域 中 都 是 离 散 的 。 ( c ) 获 得 一 个 高 密 度 频 谱 和 高 分 辨 率 频 谱 的 概 念 和 方 法 ，建 立 频 率 分 辨 率 和 时 间 分 辨 率 的 概 念 ， 为 将 来 进 一 步 进 行 时 频 分 析 （ 例 如 小 波 ） 的 学 习 和 研 究 打 下 基 础 。 2 试 验 程 序 及 截 图 ( 1 ) 常用数字信号序列的产生： 单 位 冲 激 函 数 δ ( n) c l e a r a l l t 1 =- 0 . 5 : 0 . 0 0 1 : 1 ; K= 3 ; % k n 1 = l e n g t h ( t 1 ) ; u = z e r o s ( 1 , n 1 ) ; % f r o m - 0 . 5 t o 1 t 2 =- 0 . 5 : 0 . 0 0 1 : 0 ; n 2 = l e n g t h ( t 2 ) ; u ( n 2 ) = K; p l o t ( t 1 , u ) 单 位 阶 跃 序 列 u( n) c l e a r a l l t =- 2 : 1 : 3 ; N= l e n g t h ( t ) ; x = z e r o s ( 1 , N) ; x ( 1 ) = 1 ; x ( 2 ) = 1 ; x ( 3 ) = 1 ; x ( 4 ) = 1 ; x ( 5 ) = 1 ; x ( 6 ) = 1 ; s t