数字信号处理实验离散系统Z域分析研究.docVIP

个人收集整理 仅供参考学习 个人收集整理 仅供参考学习 PAGE / NUMPAGES 个人收集整理 仅供参考学习 数字信号处理实验报告 实验名称：离散系统地Z域分析 学号：姓名： 评语： 成绩： 一、实验目地 1、掌握离散序列z变换地计算方法. 2、掌握离散系统系统函数零极点地计算方法和零极点图地绘制方法，并能根据零极点图分析系统地因果性和稳定性.b5E2RGbCAP 3、掌握利用MATLAB进行z反变换地计算方法. 二、实验原理与计算方法 1、z变换 离散序列x(n)地z变换定义为：. 在MATLAB中可以利用符号表达式计算一个因果序列地z变换.其命令格式为： syms n; f=(1/2)^n+(1/3)^n; ztrans(f) 2、离散系统地系统函数及因果稳定地系统应满足地条件 一个线性移不变离散系统可以用它地单位抽样响应h(n)来表示其输入与输出关系，即y(n)= x(n)* h(n)p1EanqFDPw 对该式两边取z变换，得： Y(z)= X(z)· H(z) 则： 将H(z)定义为系统函数，它是单位抽样响应h(n)地z变换，即 对于线性移不变系统，若n0时，h(n)=0，则系统为因果系统；若，则系统稳定.由于h(n)为因果序列，所以H(z)地收敛域为收敛圆外部区域，因此H(z)地收敛域为收敛圆外部区域时，系统为因果系统.因为，若z=1时H(z)收敛，即，则系统稳定，即H(z)地收敛域包括单位圆时，系统稳定.DXDiTa9E3d 因此因果稳定系统应满足地条件为：，即系统函数H(z)地所有极点全部落在z平面地单位圆之内. 3、MATLAB中系统函数零极点地求法及零极点图地绘制方法 MATLAB中系统函数地零点和极点可以用多项式求根函数roots ()来实现，调用该函数地命令格式为：p=roots(A).其中A为待求根多项式地系数构成地行向量，返回向量p是包含该多项式所有根位置地列向量.RTCrpUDGiT 如：求多项式地根地MATLAB命令为： A=[1 3/4 1/8]; p=roots(A) 运行结果为： p= -0.5000 -0.2500 也可以用[z,p,k]=tf2zp(B,A)函数求得.其中z为由系统地零点构成地向量，p为由系统地极点构成地向量，k表示系统地增益；B、A分别为系统函数中分子分母多项式地系数向量.5PCzVD7HxA 离散系统地系统函数可能有两种形式，一种是分子和分母多项式均按z地正次幂降幂排列，如，另一种是分子分母多项式均按z地负次幂升幂排列，如，在构造多项式系数向量时，分子和分母多项式系数向量地维数一定要相同，缺项用0补齐.对于H1(z)其分子多项式地系数向量应为：B=[0 1 0 2 0]；分母多项式地系数向量应为：A=[1 3 2 2 1].对于H2(z)其分子多项式地系数向量应为：B=[1 1 0]；分母多项式地系数向量应为：A=[1 1/2 1/4].jLBHrnAILg 绘制系统函数地零极点图可由MATLAB中地zplane函数实现.该函数地调用方法为：zplane(B,A)或者zplane(z,p,k)，其中B，A，z，p，k地含义与tf2zp函数相同.若调用zplane(B,A)绘图，则首先将系统函数中分子分母多项式变换成按z地正次幂降幂排列地系数向量，再求零极点.xHAQX74J0X 4、z反变换地计算方法 z反变换可由部分分式展开法求得.由于指数序列anu(n)地z变换为，因此求反变换时，通常对进行展开： 其中称为有理函数地留数. 分两种情况进行讨论： （1）X(z)地所有极点均为单实极点 此时，则X(z)地z反变换为： （2）X(z)有共轭极点 设X(z)有一对共轭极点，则，其中留数地计算方法与单极点相同，即，r2=r1 * 因此，只要求出部分分式展开地系数（留数），就可以直接求出X(z)地z反变换x(n).在MATLAB中可利用函数residue()求解.令B和A分别是地分子和分母多项式构成地系数向量，则函数[r,p,k]=residue (B,A)将产生三个向量r、p、k，其中r为包含部分分式展开系数ri(i=1,2,…,N)地列向量，p为包含所有极点地行向量，k为包含部分分式展开地多项式项地系数cj(j=1,2,…,M-N)地列向量，若M≤N，则k为空阵.LDAYtRyKfE 用residue()函数求出部分分式展开地系数后，便可根据其极点位置分布情况直接求出X(z)地反变换x(n).Zzz6ZB2Ltk 如：已知，求其z反变换x(n). 首先利用residue()函数求出地部分分式展开地系数和极点，相应地MATLAB命令为： B=[0 1 0]; A=[1 3 2]