版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、实验6离散时间系统的z域分析、实验目的1掌握z变换及其反变换的定义,并掌握MATLAB实现方法。2学习和掌握离散时间系统系统函数的定义及z域分析方法。3掌握系统零极点的定义,加深理解系统零极点分布与系统特性的关系、实验原理1. Z变换序列x(n)的z变换定义为X(z)x( n)zZ反变换定乂为x(n)?X (z)zn 1dz在MATLAB中,可以采用符号数学工具箱的ztrans函数和iztrans函数计算z变换和z反变换:Z=ztrans(F)求符号表达式 F的z变换。F=ilaplace(Z)求符号表达式 Z的z反变换。2. 离散时间系统的系统函数离散时间系统的系统函数H(z)定义为单位抽样
2、响应h(n)的z变换H(z)nh(n)z此外,连续时间系统的系统函数还可以由系统输入和输出信号的z变换之比得到H(z) Y(z)/X(z)MN aNZ由上式描述的离散时间系统的系统函数可以表示为boaobiz 1i aiz3. 离散时间系统的零极点分析离散时间系统的零点和极点分别指使系统函数分子多项式和分母多项式为零的点。在MATLAB中可以通过函数 roots来求系统函数分子多项式和分母多项式的根,从而得到系统的零极点。此外,还可以利用 MATLAB的zplane函数来求解和绘制离散系统的零极点分布图,zplane函数调用格式为:zplane(b,a) b,a为系统函数的分子、分母多项式的系
3、数向量(行向量)。zplane(z,p)乙p为零极点序列(列向量)。系统函数是描述系统的重要物理量,研究系统函数的零极点分布不仅可以了解系统单位抽样响应的变化,还可以了解系统的频率特性响应以及判断系统的稳定性: 系统函数的极点位置决定了系统单位抽样响应h(n)的波形,系统函数零点位置只影响冲激响应的幅度和相位,不影响波形。 系统的频率响应取决于系统的零极点,根据系统的零极点分布情况,可以通过向量分析系统的频率响应。 因果的离散时间系统稳定的充要条件是H(z)的全部极点都位于单位圆内、实验内容(1)已知因果离散时间系统的系统函数分别为 H(z)z2 2z 1z3 0.5z20.005z 0.3
4、H(z)c 4 c 323z 3z z 3z 1试采用MATLAB画出其零极点分布图,求解系统的冲激响应 h(n)和频率响应 H(e'),并判断系统是否稳定。 MATLAB代码如下:b=1 2 1;a=1 -0.5 -0.005 0.3;zplan e(b,a);b1=1 2 1;a1=1 -0.5 -0.005 0.3 0;r,p,k=residue(b1,a1)s 9CIS0口讯I Part-1.5272 - 2.2795i-1.5272 + 2.2795i-0.2790 + 0.0000i3.3333 + 0.0000i0.5198 + 0.5346i0.5198 - 0.534
5、6i-0.5396 + 0.0000i0.0000 + 0.0000i实验结果分析:由零极点分布可得冲激响应:h(n )=(-1.5272- 2.2795*i)*(0.5198+ 0.5346i)A n+(-1.5272+2.2795*i)*(0.5198 0.5346© n+(-0.2790)*(-0.5396)A n)*heaviside (n)频率响应:H(ejw)(ejw)22ejw 1(ejw)3 0.5(ejw)2 0.005ejw0.3由于该系统所有极点位于z平面单位圆内,故系统是稳定的 MATLAB代码如下: b=1 -1 0 2;111 1: * 'r-.-
6、1 -0 -r-oBr=- -ik:X j 1 110 -1id-J.b-a a115-1心Fteal Puttt.51a=3 3 -1 3 -1;zplan e(b,a);b1=1 -1 0 2;a仁3 3 -1 3 -1 0;r,p,k=residue(b1,a1)-0.1375 + 0.0000i0.2628 + 0.3222i0.2628 - 0.3222i1.6119 + 0.0000i-2.0000 + 0.0000iP =-1.6462 + 0.0000i0.1614 + 0.7746i0.1614 - 0.7746i0.3234 + 0.0000i0.0000 + 0.0000
7、i实验结果分析:由零极点分布可得冲激响应:h=(-0.1375)*(-1.6462)A n+(0.2628+0.3222*i)*(0.1614+0.7746*i)A n+(0.26280.3222*i)*(0.1614 - 0.7746*i)An+(1.6119)*(0.3234)An)*heaviside(n);频率响应:H(eJw)(ejw)3 (ejw)223(ejw)43(ejw)3 (ejw)2 3ejw 1由于该系统所有存在极点位于Z平面单位圆外,故系统是不稳定的已知离散时间系统系统函数的零点z和极点p分别为:Jj z=0,p=0.25 z=0,p=1z=0,p=-1.25 z=0
8、,p 1= 0.8e 6 ,p2= 0.8e.3J_r1.2e4.3z=0,p 1= e 8 ,p2=Jje 8 z=0,p 1= 1.2e 4 ,p2=试用MATLAB绘制上述6种不同情况下,系统函数的零极点分布图,并绘制相应单位抽样 响应的时域波形,观察分析系统函数极点位置对单位抽样响应时域特性的影响和规律MATLAB代码如下:b=1 0;a=1 -0.25;subplot(121);zplan e(b,a);subplot(122);impz(b,a,0:10);%绘岀零极点分布图%绘岀单位抽样响应得到图像如下S10fi (sarplr's)6 £0 i- 3 afl.
9、Q.o.IonReal Pa<tMATLAB代码如下:b=1 0;a=1 -1;subplot(121);zplan e(b,a);subplot(122); impz(b,a,0:10);得到图像如下i心勺 D 詰 1RlPirt5 Q flfl -0 £rd E-nip.kp丽詐“抽 MATLAB代码如下:b=1 0;a=1 1.25;subplot(121);zplan e(b,a);subplot(122);impz(b,a,0:20);得到图像如下:Irnp-jIsE RF3pTiSeWO1I1£1-曹 ua a IJU.afse-GJ40mTnm-dkwv
10、20#71>D« D D.S 1Re il PinMATLAB代码如下:z=0:-6051015ntiampisfilp=0.8*exp(i*pi/6) 0.8*exp(-i*pi/6)'subplot(121);zplan e(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);zplan e(z,p);b=1 0;a=1 -2*cos(pi/8) 1;subpl
11、ot(122);impz(b,a,0:30);得到图像如下-1.S -4 -1.5 D 0.51C 20 n(8Mnpii«)nipkp R?pCir3flMATLAB代码如下:z=0:p=1.2*exp(3*i*pi/4) 1.2*exp(-3*i*pi/4)' subplot(121);zplan e(z,p);b=1 0;a=1 -2.4*cos(-3*pi/4) 1.44;subplot(122);impz(b,a,0:30);得到图像如下:Rf-iI Faitmp J&g: 口日诃n阳航H旨trip刊J实验结果分析:由以上6种情况可以总结出当极点位于单位圆内
12、时,h(n)为衰减序列;当极点位于单位圆上时,h(n)为等幅序列;当极点位于单位圆外时,h(n)为增幅序列;若h(n)有一阶实数极点,贝U h(n)为指数序列;若h(n)有一阶共轭极点,则h(n)为指数振荡序列,并且当h(n)的极点位于虚轴左边时,h(n)按一正一负的规律交替变化(3)已知离散时间系统的系统函数分别为 H(z) Z(Z 2)匕jT(z 0.8e6)(z 0.8e 6) H(z) Z(Z 2)叱j7(z 0.8e 6)(z 0.8e 6)上述两个系统具有相同的极点,只是零点不同,试用MATLAB分别绘制上述两个系统的零极点分布图及相应单位抽样响应的时域波形,观察分析系统函数零点位置对单位抽样响应时域特性的影响。MATLAB代码如下:z=0 -2'p=0.8*exp(i*pi/6) 0.8*exp(-i*pi/6)'subplot(121);zplan e(z,p);b=1 2 0;a=1 -1.6*cos(pi/6) 0.64;subplot(122);impz(b,a,0:30);得到图像如下:|Lipij 悶 ResponseurLLA-H-H&eEMATLAB代码如下:z=0 2'p=0.8*exp(i*pi/6) 0.8*exp(-i*pi/6)'subplot(12
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024年电风扇项目资金申请报告代可行性研究报告
- 广告牌课课程设计
- 农作物分类课程设计
- 吕岩写的牧童课程设计
- 2023年强震仪投资申请报告
- 《 儿童乐园 》教学设计-2024-2025学年北师大版二年级数学上册
- 【高效备课】人教版七(上) 1.4 有理数的乘除法 1.4.2 有理数的除法 第1课时 有理数的除法 教案
- 三年级语文上册看拼音写词语
- 足球 教学设计-2023-2024学年高一上学期体育与健康人教版必修第一册
- 《乘法运算定律》(教学设计)-2023-2024学年四年级下册数学人教版
- 中医养生小妙招课件
- 2023年法律职业资格《客观题卷二》真题及答案
- 华东师大版2024七年级数学上册 第一章 有理数 单元测试卷
- 2.1.1活泼的金属单质钠说课课件高一上学期化学人教版
- 湖北省高中名校联盟2025届高三第一次联合测评 历史试卷(含答案详解 )
- Unit 3 Amazing animals大单元备课教案
- 新苏教版三年级上册科学全册复习资料 (背诵用)
- 2024版煤粉采购合同范本
- 专利培训课件
- 清雪合同范本
- 《秋天的画报》课件
评论
0/150
提交评论