实验6离散时间系统的z域分析
- 格式:doc
- 大小:104.15 KB
- 文档页数:12
实验6 离散时间系统的z 域分析一、实验目的1.掌握z 变换及其反变换的定义,并掌握MAT LAB 实现方法。
2.学习和掌握离散时间系统系统函数的定义及z 域分析方法。
3.掌握系统零极点的定义,加深理解系统零极点分布与系统特性的关系。
二、实验原理1. Z变换序列x(n )的z变换定义为()()nn X z x n z+∞-=-∞=∑Z 反变换定义为11()()2n rx n X z z dzj π-=⎰在MA TLAB 中,可以采用符号数学工具箱的ztr ans 函数和iztr an s函数计算z 变换和z 反变换:Z=z trans(F) 求符号表达式F 的z 变换。
F=ilaplace(Z) 求符号表达式Z 的z 反变换。
2.离散时间系统的系统函数离散时间系统的系统函数H(z)定义为单位抽样响应h(n)的z 变换()()nn H z h n z+∞-=-∞=∑此外,连续时间系统的系统函数还可以由系统输入和输出信号的z变换之比得到()()/()H z Y z X z =由上式描述的离散时间系统的系统函数可以表示为101101()MM N N b b z b z H z a a z a z ----+++=+++……3.离散时间系统的零极点分析离散时间系统的零点和极点分别指使系统函数分子多项式和分母多项式为零的点。
在MATL AB 中可以通过函数roots 来求系统函数分子多项式和分母多项式的根,从而得到系统的零极点。
此外,还可以利用MATLAB 的zplane 函数来求解和绘制离散系统的零极点分布图,zpla ne 函数调用格式为:zp lane(b,a) b,a 为系统函数的分子、分母多项式的系数向量(行向量)。
zpl an e(z ,p) z,p 为零极点序列(列向量)。
系统函数是描述系统的重要物理量,研究系统函数的零极点分布不仅可以了解系统单位抽样响应的变化,还可以了解系统的频率特性响应以及判断系统的稳定性:①系统函数的极点位置决定了系统单位抽样响应h(n)的波形,系统函数零点位置只影响冲激响应的幅度和相位,不影响波形。
数字信号处理实验报告实验名称:离散系统的Z 域分析学号:姓名: 评语: 成绩: 一、实验目的1、掌握离散序列z 变换的计算方法。
2、掌握离散系统系统函数零极点的计算方法和零极点图的绘制方法,并能根据零极点图分析系统的因果性和稳定性。
3、掌握利用MATLAB 进行z 反变换的计算方法。
二、实验原理与计算方法1、z 变换离散序列x (n )的z 变换定义为:。
∑∞-∞=-=n n z n x Z X )()(在MATLAB 中可以利用符号表达式计算一个因果序列的z 变换。
其命令格式为:syms n; f=(1/2)^n+(1/3)^n;ztrans(f)2、离散系统的系统函数及因果稳定的系统应满足的条件一个线性移不变离散系统可以用它的单位抽样响应h (n )来表示其输入与输出关系,即y (n )= x (n )* h (n )对该式两边取z 变换,得: Y (z )= X (z )· H (z )则: )()()(z X z Y z H =将H (z )定义为系统函数,它是单位抽样响应h (n )的z 变换,即∑∞-∞=-==n n z n h n h Z z H )()]([)(对于线性移不变系统,若n <0时,h (n )=0,则系统为因果系统;若,则系统稳∞<∑∞-∞=n n h |)(|定。
由于h (n )为因果序列,所以H (z )的收敛域为收敛圆外部区域,因此H (z )的收敛域为收敛圆外部区域时,系统为因果系统。
因为,若z =1时H (z )收敛,即∑∞-∞=-=n n z n h z H )()(,则系统稳定,即H(z)的收敛域包括单位圆时,系统稳定。
∞<=∑∞-∞==n z n h z H |)(||)(1因此因果稳定系统应满足的条件为:,即系统函数H (z )的所有极点全部落在1,||<∞≤<ααz z 平面的单位圆之内。
3、MATLAB 中系统函数零极点的求法及零极点图的绘制方法MATLAB 中系统函数的零点和极点可以用多项式求根函数roots ()来实现,调用该函数的命令格式为:p=roots(A)。
《信号与系统》课程实验报告变换。
zz z z z z F 2112)(232+++-=一、实验原理的验证 1、离散系统零极点图实验原理如下:离散系统可以用差分方程描述:∑∑==-=-Mm m Ni i m k f b i k y a 0)()(Z 变换后可得系统函数:NN MM z a z a a z b z b b z F z Y z H ----++++++==......)()()(110110 可以用root 函数可分别求零点和极点。
例7-4 求系统函数零极点图131)(45+-+=z z z z H实验结果如下:2、离散系统的频率特性实验原理如下:离散系统的频率特性可由系统函数求出,既令ωj e z =,函数freqz 可计算频率特性,调用格式是:[H ,W]=freqz(b,a,n),b 和a 是系统函数分子分母系数,n 是π-0范围内n 个等份点,默认值为512,H 是频率响应函数值,W 是相应频率点; 例7-5 系统函数z z z H 5.0)(-=10个频率点的计算结果为幅频特性曲线相频特性曲线freqz语句直接画图例7-7已知系统函数114/11)1(4/5)(----=z z z H ,画频率响应和零极点图。
零极点图幅频特性曲线相频特性曲线二、已知离散系统的系统函数如下所示:1422)(232+-++=z z z z z H试用MATLAB 实现下列分析过程: (1)求出系统的零极点位置;(2)绘出系统的零极点图,根据零极点图判断系统的稳定性; (3)绘出系统单位响应的时域波形,并分析系统稳定性与系统单位响应时域特性的关系。
(1)由计算结果可知:系统的极点为p0=-3.3028、p1=1、p2=0.3028。
由计算结果可知:系统的零点为z0=1.4142i 、z1=-1.4142i 。
(2)系统的零极点图如下:程序清单如下: a=[1 2 -4 1]; b=[1 0 2]; ljdt(a,b)p=roots(a)q=roots(b)pa=abs(p)由图可知:第一个极点(p0)在单位圆外部,第二个极点(p1)在单位圆上,第三个极点(p2)在单位圆内部,因为有一个极点在单位圆外部,故该系统是不稳定的系统(稳定系统要求极点全部在单位圆内)。
第 六 章 离散信号与系统的 Z 域分析引言与线性连续系统的频域分析和复频域分析类似,线性离散系统的频域分析是输入信号分解为基本信号e jΩk 之和,则系统的响应为基本信号的响应之和。
这种方法的数学描述是离散时间傅里叶变换和逆变换。
如果把复指数信号e jΩk 扩展为复指数信号Z k ,Z=re jΩ ,并以Zk 为基本信号, 把输入信号分解为基本信号Z k 之和, 则响应为基本信号Z k 的响应之和。
这种方法的数学描述为Z 变换及其逆变换,这种方法称为离散信号与系统的Z 域分析法.如果把离散信号看成连续时间信号的 抽样值序列,则Z 变换可由拉普拉斯变换引入.因此离散信号与系统的Z 域分析 和连续时间信号与系统的复频域分析有许多相似之处.通过Z 变换,离散时间信 号的卷积运算变成代算,离散时间系统的差分方程变成Z 域的代数方程,因此可 以比较方便的分析系统的响应。
Z 变换从拉普拉斯变换到Z 变换对连续信号f(t)进行理想抽样,即f(t)乘以单位冲击序列δT (t),T 为 抽样间隔,得到抽样信号为f s (t)=f(t)δT (t)= =对fs(t)取双边拉普拉斯变换,得F s (s)=£[fs(t)]=令z=e sT , 则Fs(s)=F(z) ,得F(z)=因为T为常数,所以通常用f(k)表示f(kT),于是变为F(z)=称为f(k)的双边Z变换,z为复变量。
z和s的关系为:z=e sTs=(1/T)㏑z由复变函数理论,可以得到f(k)= ∮cF(z)z k-1 dz式(7.1-5)称为F(z)的双边Z逆变换(后面讨论).双边Z变换的定义和收敛域§双边 Z 变换的定义对于离散序列f(k)(k=0,±1,±2,┄),函数(z的幂级数)F(z)=称为f(k)的双边Z变换,记为F(z)=Z[f(k)].F(z)又称为f(k)的象函数,f(k)又 称为F(z)的原函数.为了表示方便,f(k)与F(z)之间的对应关系可表示为 f(k) F(z)§双边 Z 变换的收敛域f(k)的双边Z变换为一无穷级数,因此存在级数是否收敛的问题.只有当 (7.1-6)式的级数收敛,F(z)才存在.F(z)存在或级数收敛的充分条件是 ∞在f(k)给定的条件下,式(7.1-6)级数是否收敛取决于z的取值.在z复平面上, 使级数收敛的z取值区域称为F(Z)的收敛域。
第 6 章离散信号与系统的Z 域分析6.0 引言与拉氏变换是连续时间傅立叶变换的推广相对应,Z 变换是离散时间傅立叶变换的推广。
Z 变换的基本思想、许多性质及其分析方法都与拉氏变换有相似之处。
当然, Z 变换与拉氏变换也存在着一些重要的差异。
6.1 双边 Z 变换6.1.1双边Z变换的定义前面讨论过,单位脉冲响应为h[n] 的离散时间 LTI 系统对复指数输入z n的响应y[n]为y[ n]H ( z) z n(6.1)其中H ( z)h[ n] z n(6.2)n式 (6. 2) 就称为 h[n] 的双边 Z 变换。
当 z= e j时, Z 变换就转变为傅立叶变换。
因此一个离散时间信号的双边Z 变换定义为:X ( z)x[ n]z n(6.3)n式中 z 是一个复变量。
而x[n]与它的双边z 变换之间的关系可以记做zx[n]X (z)6.1.2双边Z变换的收敛域x[n] 的双边 Z 变换为一无穷级数,因此存在级数是否收敛的问题,即一方面并非所有信号的Z 变换都存在;另一方面即使某信号的Z 变换存在,但并非Z 平面上的所有点都能使X(z)收敛。
那些能够使X(z)存在的点的集合,就构成了X(z)的收敛域,记为ROC。
只有当式 (6.3) 的级数收敛,X (z) 才存在。
X ( z) 存在或级数收敛的充分条件是x[n]z n(6.4)n在 x[ n] 给定的条件下,式 (6.4)级数是否收敛取决于 z 的取值。
在 z 复平面上,使式 (6.4)级数收敛的 z取值区域就是 X(z)的收敛域。
6.1.3零极点图如果X(z) 是有理函数,将其分子多项式与分母多项式分别因式分解可以得到:N ( z)(z z i )X ( z)i(6.5)M(zD ( z)z p )p则由其全部的零极点即可表示出X ( z) ,最多相差一个常数因子。
在Z 平面上表示出全部的零极点,即构成X ( z) 的几何表示——零极点图。
实验6 离散时间系统的z 域分析一、实验目的1.掌握z 变换及其反变换的定义,并掌握MATLAB 实现方法。
2.学习和掌握离散时间系统系统函数的定义及z 域分析方法。
3.掌握系统零极点的定义,加深理解系统零极点分布与系统特性的关系。
二、实验原理1. Z 变换序列x(n)的z 变换定义为()()nn X z x n z+∞-=-∞=∑Z 反变换定义为11()()2n rx n X z z dzj π-=⎰在MATLAB 中,可以采用符号数学工具箱的ztrans 函数和iztrans 函数计算z 变换和z 反变换:Z=ztrans(F) 求符号表达式F 的z 变换。
F=ilaplace(Z) 求符号表达式Z 的z 反变换。
2.离散时间系统的系统函数离散时间系统的系统函数H(z)定义为单位抽样响应h(n)的z 变换()()nn H z h n z+∞-=-∞=∑此外,连续时间系统的系统函数还可以由系统输入和输出信号的z 变换之比得到()()/()H z Y z X z =由上式描述的离散时间系统的系统函数可以表示为101101()MM N N b b z b z H z a a z a z ----+++=+++……3.离散时间系统的零极点分析离散时间系统的零点和极点分别指使系统函数分子多项式和分母多项式为零的点。
在MATLAB 中可以通过函数roots 来求系统函数分子多项式和分母多项式的根,从而得到系统的零极点。
此外,还可以利用MATLAB 的zplane 函数来求解和绘制离散系统的零极点分布图,zplane 函数调用格式为:zplane(b,a) b,a 为系统函数的分子、分母多项式的系数向量(行向量)。
zplane(z,p) z,p 为零极点序列(列向量)。
系统函数是描述系统的重要物理量,研究系统函数的零极点分布不仅可以了解系统单位抽样响应的变化,还可以了解系统的频率特性响应以及判断系统的稳定性:①系统函数的极点位置决定了系统单位抽样响应h(n)的波形,系统函数零点位置只影响冲激响应的幅度和相位,不影响波形。
②系统的频率响应取决于系统的零极点,根据系统的零极点分布情况,可以通过向量分析系统的频率响应。
③因果的离散时间系统稳定的充要条件是H(z)的全部极点都位于单位圆内。
三、实验内容(1)已知因果离散时间系统的系统函数分别为:①23221()0.50.0050.3z z H z z z z ++=--+②324322()3331z z H z z z z z -+=+-+-试采用MATLAB 画出其零极点分布图,求解系统的冲激响应h(n)和频率响应H(),并判断系统是否稳定。
①MATLAB 代码如下: b=[1 2 1];a=[1 -0.5 -0.005 0.3]; zplane(b,a); b1=[1 2 1];a1=[1 -0.5 -0.005 0.3 0]; [r,p,k]=residue(b1,a1) r =-1.5272 - 2.2795i -1.5272 + 2.2795i -0.2790 + 0.0000i 3.3333 + 0.0000i p =0.5198 + 0.5346i 0.5198 - 0.5346i -0.5396 + 0.0000i 0.0000 + 0.0000i k = []Ωj e实验结果分析:由零极点分布可得 冲激响应:h(n)=((-1.5272 - 2.2795*i)*(0.5198 + 0.5346i)^n+(-1.5272 + 2.2795*i)*(0.5198 - 0.5346*i)^n+(-0.2790)*(-0.5396)^n)*heaviside(n) 频率响应:232()21()()0.5()0.0050.3jw jw jwjw jw jwe e H e e e e ++=--+ 由于该系统所有极点位于Z 平面单位圆内,故系统是稳定的。
②MATLAB 代码如下: b=[1 -1 0 2]; a=[3 3 -1 3 -1]; zplane(b,a); b1=[1 -1 0 2]; a1=[3 3 -1 3 -1 0]; [r,p,k]=residue(b1,a1) r =-0.1375 + 0.0000i 0.2628 + 0.3222i 0.2628 - 0.3222i 1.6119 + 0.0000i -2.0000 + 0.0000i p =-1.6462 + 0.0000i0.1614 + 0.7746i 0.1614 - 0.7746i 0.3234 + 0.0000i 0.0000 + 0.0000i k = [] 实验结果分析: 由零极点分布可得 冲激响应:h=((-0.1375)*(-1.6462)^n+(0.2628 + 0.3222*i)*(0.1614 + 0.7746*i)^n+(0.2628 - 0.3222*i)*(0.1614 - 0.7746*i)^n+(1.6119)*(0.3234)^n)*heaviside(n); 频率响应:32432()()2()3()3()()31jw jw jwjw jw jw jwe e H e e e e e -+=+-+- 由于该系统所有存在极点位于Z 平面单位圆外,故系统是不稳定的。
(2)已知离散时间系统系统函数的零点z 和极点p 分别为:①z=0,p=0.25 ②z=0,p=1 ③z=0,p=-1.25 ④z=0,p 1=60.8je π,p 2=60.8jeπ- ⑤z=0,p 1= 8jeπ,p 2= 8jeπ-⑥z=0,p 1= 341.2jeπ,p 2= 341.2jeπ-试用MATLAB 绘制上述6种不同情况下,系统函数的零极点分布图,并绘制相应单位抽样响应的时域波形,观察分析系统函数极点位置对单位抽样响应时域特性的影响和规律。
①MATLAB 代码如下:b=[1 0];a=[1 -0.25];subplot(121);zplane(b,a); %绘出零极点分布图subplot(122);impz(b,a,0:10); %绘出单位抽样响应得到图像如下:②MATLAB代码如下:b=[1 0];a=[1 -1];subplot(121);zplane(b,a);subplot(122);impz(b,a,0:10);得到图像如下:③MATLAB代码如下:b=[1 0];a=[1 1.25];subplot(121);zplane(b,a);subplot(122);impz(b,a,0:20);得到图像如下:④MATLAB代码如下:z=[0]';p=[0.8*exp(i*pi/6) 0.8*exp(-i*pi/6)]'; subplot(121);zplane(z,p);b=[1 0];a=[1 -1.6*cos(pi/6) 0.64];subplot(122);impz(b,a,0:30);得到图像如下:⑤MATLAB代码如下:z=[0]';p=[exp(i*pi/8) exp(-i*pi/8)]'; subplot(121);zplane(z,p);b=[1 0];a=[1 -2*cos(pi/8) 1];subplot(122);impz(b,a,0:30);得到图像如下:⑥MATLAB代码如下:z=[0]';p=[1.2*exp(3*i*pi/4) 1.2*exp(-3*i*pi/4)]'; subplot(121);zplane(z,p);b=[1 0];a=[1 -2.4*cos(-3*pi/4) 1.44];subplot(122);impz(b,a,0:30);得到图像如下:实验结果分析:由以上6种情况可以总结出:①当极点位于单位圆内时,h(n)为衰减序列; ②当极点位于单位圆上时,h(n)为等幅序列; ③当极点位于单位圆外时,h(n)为增幅序列; ④若h(n)有一阶实数极点,则h(n)为指数序列;⑤若h(n)有一阶共轭极点,则h(n)为指数振荡序列,并且当h(n)的极点位于虚轴左边时,h(n)按一正一负的规律交替变化。
(3)已知离散时间系统的系统函数分别为:①66(2)()(0.8)(0.8)jjz z H z z ez eππ-+=--②66(2)()(0.8)(0.8)jjz z H z z ez eππ--=--上述两个系统具有相同的极点,只是零点不同,试用MATLAB分别绘制上述两个系统的零极点分布图及相应单位抽样响应的时域波形,观察分析系统函数零点位置对单位抽样响应时域特性的影响。
①MATLAB代码如下:z=[0 -2]';p=[0.8*exp(i*pi/6) 0.8*exp(-i*pi/6)]';subplot(121);zplane(z,p);b=[1 2 0];a=[1 -1.6*cos(pi/6) 0.64];subplot(122);impz(b,a,0:30);得到图像如下:②MATLAB代码如下:z=[0 2]';p=[0.8*exp(i*pi/6) 0.8*exp(-i*pi/6)]';subplot(121);zplane(z,p);b=[1 -2 0];a=[1 -1.6*cos(pi/6) 0.64];subplot(122);impz(b,a,0:30);得到图像如下:实验结果分析:从图像看出,两个系统极点相同,零点互为相反数,得到的h(n)各值也对应相反,但收敛性一致,故在有相同极点的情况下,零点分布只影响系统时域响应的幅度,不影响响应模式。
四、实验心得体会。