离散系统分析.docx

本科学生实验报告 学号 姓名 学院 物电学院 专业、班级 1 实验课程名称 教师及职称 开课学期 2014 至 2015 学年 下 学期 填报时间 2015 年 5 月 14 日 云南师范大学教务处编印 实验序号 2 实验名称 离散系统分析 实验时间 2015,5,13 实验室 同析 3-312 实验预习 1．实验目的 深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意义，熟练掌握利用MATLAB 分析离散系统的时域响应、频响特性和零极点的方法。掌握利用 DTFT 和 DFT 确定系统特性的原理和方法。2．实验原理、实验流程或装置示意图 MATLAB 提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频域响应等分析函数。 离散系统的时域响应 离散时间 LTI 系统可用如下的线性常系数差分方程来描述： y[k ] ? a 1 y[k ? 1] ?? ? a  N ?1 y[k ? (N ? 1)] ? a N y[k ? N ] ? b x[k ] ? b x[k ? 1] ?? ? b x[k ? (M ? 1)] ? b x[k ? M ] 0 1 M ?1 M 已知输入信号 x[k]以及系统初始状态 y[-1]，y[-2] ...，就可以求出系统的响应。MATLAB 提供了计算差分方程的数值解的函数，可以计算上述差分方程描述的离散系统的单位脉冲响应、零状态响应和完全响应等。 在调用 MATLAB 函数时，需要利用描述该离散系统的系数函数。对差分方程进行 Z 变换即可得系统函数： Y (z) b ? b z?1 ?? ? b z?(M ?1) ? b z? M b(z) H (z) ? ? 0 1 M ?1 M ? X (z) 1? a 1 z?1 ?? ? a  N ?1 z?( N ?1) ? a N z? N a(z) 在 MATLAB 中可使用向量 a 和向量 b 分别保存分母多项式和分子多项式的系数： a ? [1, a , ? , a , a ] b ? [b , b , ? , b , b ] 1 N ?1 N 0 1 M ?1 M 这些系数均从 z0 按 z 的降幂排列。 离散系统的单位冲激响应 h[k]的计算 [h,k] = impz(b, a)：计算系统的单位脉冲响应 h[k]和相应的时间向量 k；也可简写为： h = impz(b, a)。其中： a ? [1, a , ? , a , a ] b ? [b , b , ? , b , b ] 1 N ?1 N 0 1 M ?1 M [h,k] = impz(b, a, n)：计算 n 点单位脉冲响应 h[k]； 也可简写为：h = impz(b, a, n)。 impz(b, a)：绘制单位脉冲响应 h[k]的图形。 离散系统响应 y[k]的计算 y = filter(b, a, x)： 计算系统在输入 x 作用下的零状态响应 y[k]； a ? [1, a , ? , a , a ] b ? [b , b , ? , b , b ] 1 N ?1 N 0 1 M ?1 M y = filter(b, a, x, zi)： 计算系统在输入 x 和初始状态作用下的完全响应 y[k]。 zi 是由系统的初始状态经过 filtic 函数转换而得到的初始条件：zi= filtic (b, a, Y0) ，Y0 为系统的初始状态， Y0= [ y[-1]，y[-2]，y[-3]，...]。 离散系统的系统函数零极点分析 离散 LTI 系统的系统函数 H(z)可以表示为零极点形式： Y (z) b ? b z?1 ?? ? b z ?(M ?1) ? b z? M (z ? z )(z ? z )...(z ? z ) H (z) ? ? 0 1 M ?1 M ? k 1 2 M X (z) 1? a 1 z ?1 ?? ? a N ?1 z?( N ?1) ? a N z ? N (z ? p 1 )(z ? p 2 )...(z ? p ) N 使用 MATLAB 提供的 roots 函数计算离散系统的零极点； 使用 zplane 函数绘制离散系统的零极点分布图。 注意：在利用这些函数时，要求 H(z)的分子多项式和分母多项式的系数的个数相等，若不等则需要补零。 离散系统的频率响应 当离散因果 LTI 系统的系统函数 H(z)的极点全部位于 z 平面单位圆内时，系统的频率响应可由 H(z)求出，即 H (ej? ) ? H (z)  z?ej? ? H (ej? ) ej?(?