线性时不变离散时间系统的频域分析课件.docVIP

云南大学信息学院数字信号处理实验报告 实验四线性时不变离散时间系统的频域分析 实验室名称：格物楼 1204 实验时间： 2015 年11 月 6 日 姓名： 成绩 教师签名： 年 月 日 一、实验目的 （1）学习 MATLAB 软件及其在信号处理中的应用， 加深对常用连续时间信号的理解。 （2）连续时间信号在时域和频域中的抽样效果。 （3）熟悉 MATLAB 模拟滤波器的。 二、实验内容 Q4.1 修改程序 P3.1中，取三个不同的 M 值，当0 2 时计算并画出式 (2.13)所示 滑动平均滤波器的幅度和相位谱。证明由幅度和相位谱表现出的对称类型。它表示了 那种类型的滤波器？你现在能解释习题 Q2.1 的结果吗？修改后的程序 P3.1也可以用 于计算并画出线性时不变离散时间系统的频率响应， 该系统用形如式 (4.12)的传输函数 描述。 Q4.2 使 用 修 改 后 的 程 序 P3.1， 计 算 并 画 出 当 0 2 时 传 输 函 数 2 0.15(1 z ) H (z) 的因果线性时不变离散时间系统的频率响应。它表示那种类 1 2 1 0.5 z 0.7 z 型的滤波器？ Q4.3 对下面的传输函数重做习题 Q4.2 2 0.15(1 z ) G 4.2 和4.3 给出的两个 (z) 1 2 0.7 0.5z z 滤波器之间的区别是什么？你将选择哪一个滤波器来滤波，为什么？ Q4.6 使用 zplane分别生成 4.2和4.2 所确定的两个滤波器的零极点图。 讨论你的结果。 Q4.7 用程序 P4.1 计算并画出近似理想低通滤波器的冲激响应。低通有限冲激响应滤 波器的长度是多少？在程序 P4.1 中，那个语句确定滤波器的长度？那个参数控制截止 频率？ Q4.8 修改程序 P4.1 ，计算并画出式 (4.39)所示的长度为 20，截止角频率为 c 0.45的 有限冲激响应低通滤波器的冲激响应。 Q4.9 修改程序 P4.1 ，计算并画出式 (4.39)所示的长度为 15，截止角频率为 c 0.65的 有限冲激响应低通滤波器的冲激响应。 1 云南大学信息学院数字信号处理实验报告 Q4.10 编写一个 MATLAB 程序，计算并画出式 (4.39)所示有限冲激响应低通滤波器的 振幅响应。使用这个程序，选取几个不同的 N 值，画出振幅响应并讨论你的结果。 Q4.11 运行程序 P4.2，计算并画出一个长度为 2 的滑动平均滤波器的增益响应。从图 中验证 3dB 截止频率在 π/2 处。 Q4.23 用 MATLAB 产 生 如 下 两 个 因 果 系 统 传 输 函 数 的 零 极 点 图 1 H (z) 研究生成的零极点图，你可以推断它们的稳定性吗？ 1 1 1.848 0.85 1 2 z z 三、实验器材及软件 1. 微型计算机 1 台 2. MATLAB 7.0软件 四、实验原理 熟悉 MATLAB 中产生信号和绘制信号的基本命令， MATLAB 已被开发成能对数 据向量或矩阵运算的工具。序列以向量的形式储存，并且所有的信号被限定为因果的 和有限长的，采用何种步骤执行程序，要根据 MATLAB 所运行的平台来决定。 (1) 若{h[n]} 表示一个线性时不变离散时间系统的冲激响应，对 {h[n]} 做离散时间傅里 j 叶变换得到其频率响应 H (e ) ，即 j h n e j n H (e ) [ ] 。 n j (2) 通常， H (e ) 是一个周期为 2π的ω的复值函数，可以根据实部，虚部或者幅度相 位来表示该函数。 因此， j j H ej jH ej H ej ej ( ) H re (e ) H (e ) ( ) ( ) | ( ) | 其中 和 re im H j im (e ) j j j 分别是 H (e ) 的实部和虚部， 并且 ( ) arg{H (e )} 又|H (e )| 称为幅度响 ( ) 应，而 称为线性时不变离散时间系统的相位响应。 j (3) 线性时不变系统的增益函数 g(ω)定义为 g( ) 20 log10 | H (e ) | dB，增益函数的 相反数 a(ω)=-g(ω),称为衰减或损益函数。 j (4) 对于用实冲激响应 h[n]描述的离散时间系统，幅度函数是 ω的偶函数，即 |H (e )| j Hre (e ) j =|H (e ) | ；而相位函数是 ω的奇函数，即 ( ) ( ) 。同样， 是ω的偶 j Him (e ) 函数， 是ω 的奇函数。 j (5) 线性时不变离散时间系统的频率响应可以由输出序列 y[n] 的傅里叶变换 Y(e ) 与 2 云南大学信息学院数字信号处理实验报告 j j j H (e ) Y(e ) X(e ) j 输入序列 x[n] 的