数字信号处理MATLAB实验报告.docxVIP

实验报 课程名称： 数字信号处理实验 专业班级： 姓 名： 学 号： 实验名称 数字信号处理 实验地点 实验时间 实验成绩 、实验目的及任务 学会运用 MATLAB求离散时间信号的 z 变换和 z 反变换； 学会运用 MATLAB分析离散时间系统的系统函数的零极点； 学会运用 MATLAB分析系统函数的零极点分布与其时域特性的关 系； 学会运用 MATLAB进行离散时间系统的频率特性分析。 、 实验内容与步骤 z 正反变换 序列 x n 的 z 变换定义为 X z Z x n n xnz n (4-1) 其中，符号 Z 表示取 z 变换， z 是复变量。相应地，单边 z 变换定义为 X z Z x n n n xnz 0 (4-2) MATLAB符号数学工具箱提供了计算离散时间信号单边 ztrans 和 z 反变换函数 iztrans ，其语句格式分别为 z 变换的函数 Z=ztrans(x) x=iztrans(z) 上式中的 x 和 Z 分别为时域表达式和 z 域表达式的符号表示，可通过 sym 函数来定义。 如果信号的 z 域表示式 X(z)是有理函数，进行 z 反变换的另一个方法 是对 X ( z)进行部分分式展开， 然后求各简单分式的 z 反变换。设 X(z)的有 理分式表示为 X(z)b0 b1z X(z) b0 b1z 1 b2 z 2 1 a1z 1 a2 z 2 bm z m B(z) an z n A(z) 4-3) MATLAB信号处理工具箱提供了一个对 X(z) 进行部分分式展开的函数 residuez ，其语句格式为 [R,P,K]=residuez(B,A) 其中， B，A分别表示 X(z) 的分子与分母多项式的系数向量； R为部分分式 的系数向量； P为极点向量； K为多项式的系数。若 X(z) 为有理真分式， 则 K 为零。 系统函数的零极点分析 z 变换与激励的 z z 变换与激励的 z 4-4) 变换之比，即 H(z) YX((zz)) 如果系统函数 H(z) 的有理函数表示式为 H(z)b1zm H(z) b1zm b2zm 1 a1 zn a2zn 1 bmz anz bm 1a n 1 4-5) 那么，在 MATLAB中系统函数的零极点就可通过函数 roots 得到，也可借 助函数 tf2zp 得到， tf2zp 的语句格式为 [Z,P,K]=tf2zp(B,A) 其中，B与 A分别表示 H ( z)的分子与分母多项式的系数向量。它的作用是 将 H ( z) 的有理分式表示式转换为零极点增益形式，即 H(z) k (z z1)(z z2) (z zm) (4-6) (z p1)(z p2 ) (z pn ) 若要获得系统函数 H (z)的零极点分布图，可直接应用 zplane 函数， 其语句格式为 zplane(B,A) 其中，B与 A分别表示 H ( z)的分子和分母多项式的系数向量。它的作用是 在 Z 平面上画出单位圆、零点与极点。 系统函数的零极点分布与其时域特性的关系 与拉氏变换在连续系统中的作用类似，在离散系统中， z 变换建立了 时域函数 h(n) 与 z 域函数 H ( z)之间的对应关系。因此， z 变换的函数 H(z) 从形式可以反映 h(n) 的部分内在性质。 我们仍旧通过讨论 H(z) 的一阶极点 情况，来说明系统函数的零极点分布与系统时域特性的关系。 离散时间 LTI 系统的频率特性分析 对于因果稳定的离散时间系统，如果激励序列为正弦序列 x(n) A sin(n )u(n) ， 则系统 的 稳态响 应为 yss(n) A|H(ej ) |sin[n ( )]u(n) 。其中， H (e j ) 通常是复数。 离散时间 系统的频率响应定义为 H (ej ) |H (ej )|ej ( ) (14-7 ) 其中， |H(ej ) |称为离散时间系统的幅频特性； ( ) 称为离散时间系统的 相频特性； H(ej )是以 s( s 2 ，若零 T 1， s 2 )为周期的周期 函数。因此，只要分析 H(ej ) 在| | 范围内的情况，便可分析出系统的 整个频率特性。 MATLAB提供了求离散时间系统频响特性的函数 freqz ，调用 freqz 的 格式主要有两种。一种形式为 [H,w]=freqz(B,A,N) 其中，B与 A分别表示 H(z) 的分子和分母多项式的系数向量； N为正整数， 默认值为 512；返回值 w包含 [0, ]范围内的 N个频率等分点；返回值 H则 是离散时间系统频率响应 H(ej )在 0~ 范围内 N个频率处的值。另一种形 式为 [H,w]=freqz(B,A,N, whole ) 与第一种方式不同之处在于角频率的范围由 [0,