利用双线性变换求其离散传递函数

课程设计报告学生姓名:学号：学院:班级:题目:专业方向课程设计仿真组-利用双线性变换求其离散传递函数指导教师：孟杰姜文娟职称:副教授讲师2017年1月13日目录1.题目背景与意义................................................12.设计题目介绍..................................................12.1设计内容和要求..........................................12.2设计工作任务及工作量的要求..............................13双线性变换法..................................................13.1双线性变换的定义........................................13.2双线性变换法的优缺点....................................23.2.1双线性变换法的优点................................23.2.2双线性变换法的缺点................................23.2双线性变换的原理........................................23.3双线性变换的主要特性....................................24设计步骤......................................................35理论计算......................................................56结果分析......................................................6参考文献.......................................................10附录1程序清单................................................11附录2计算机实现程序框图......................................1711.题目背景与意义本课程设计以自动控制理论、现代控制理论、MATLAB及应用等知识为基础，利用双线性变换求连续系统对应的离散化的系统，目的是使学生在现有的控制理论的基础上，学会用MATLAB语言编写控制系统的离散化的程序，通过上机实习加深对课堂所学知识的理解，掌握一种能方便地对系统进行离散化的设计工具。2.设计题目介绍2.1设计内容和要求1在理论上对连续系统采用双线性变换求离散化推导出算法和计算公式2画出计算机实现算法的框图3编写程序并调试和运行4以下面的系统为例，进行计算已知系统闭环传递函数)2)(1(4)(ssssG，利用双线性变换求其离散传递函数。5分析运算结果（离散化步长对系统性能的影响）6程序应具有一定的通用性，对不同参数能有兼容性。2.2设计工作任务及工作量的要求1本次课程设计要求每周学生至少见指导教师2次，其中集中辅导答疑部不于3次。2设计说明书的格式按设计说明书格式要求，采用word软件排版，计算机打印。（具体包括：封皮、目录、正文、参考文献等）3程序清单用A4纸打印后，作为附录订装在说明书后面。4框图和其他图表放在正文中。3双线性变换法3.1双线性变换的定义双线性变换法又称突斯汀（Tustin）法，是一种基于梯形积分规则的数字积分变换方法。23.2双线性变换法的优缺点3.2.1双线性变换法的优点双线性变换的主要优点:靠频率的严重非线性关系得到S平面与Z平面的单值一一对应关系，整个jΩ轴单值对应于单位圆一周，其中ω和Ω为非线性关系。在零频率附近，Ω~ω接近于线性关系，Ω进一步增加时，ω增长变得缓慢，(ω终止于折叠频率处)，所以双线性变换不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象。3.2.2双线性变换法的缺点双线性变换法的缺点:Ω与ω的非线性关系，导致数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变，(使数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变)。例如，一个模拟微分器，它的幅度与频率是线性关系，但通过双线性变换后，就不可能得到数字微分器。3.2双线性变换的原理双线性变换公式可以从梯形积分公式中直接推导出来。按这种替换公式进行替换，可以保证)(zG的稳定性，而且，具有一定的仿真精度。已知梯形积分公式为：)(211kkkkxxTxx即：xzTxz)1(2)1（则有：1121zzTs即：112zzTs（3.1）（3.1）式称为双线性变换公式，也可写成为：2121sTsTz（3.2）3.3双线性变换的主要特性1、若)(sD稳定，则)(zD一定稳定；32、变换前后，稳态增益不变；3、双线性变换后的阶次不变，且分子、分母具有相同的阶次。4设计步骤双线性替换不仅可以精度极高的仿真模型，而且能利用计算机程序来实现这种替换。下面来介绍一种程序替换法。折线形系统的传递函数为nnnnnnnnbsbsbsbasasasasEsUsG111011100......)()()(（4.1）在双线性替换下得到的Z传递函数为：；nnnnnnnnezezezedzdzdzdzEzUzG111011100......)()()(（4.2）现需要由iiba，（i=0，1，......，n）确定iied，（i=0，1，......，n），若直接将双线性替换公式代入)(0sG，可得：nnnnnnnnnnnnbzzTbzzTbzzTbazzTazzTazzTazG)11)(2(...)11()2()11()2()11)(2(...)11()2()11()2()(11110111100（4.3）将其分子、分母同时乘以nz）（1，可得：)()()1(...)1()1(2)1(2)1(...)1()1(2)1(2)(111011100zBzAzbzzTbzTbzazzTazTazGnnnnnnnnnnnn）（）（）（）（（4.4）将)(zA，)(zB写成向量形式，nnnnnnnzzzzTTTTaaaazA)1(...)1()1()1()2(0)2(...)2(0)2(...)(11100114由于nz）（1，，......，nz）（1均为n阶多项式，可得到：1........................)1(...)1(1...)1(111101011110001001zzzxxxxxxxxxxxxzzzzzznnnnnniniinniiinnn）（）（）（矩阵ijx为(n+1)×（n+1）阶；其中第一行诸元素为nz）（1的展开式的各系数，第一列诸元素为1，阶次n确定后，这些元素均为已知，并可以证明，其余n×n个元素可由下式求得：1,11,,1,jijijijixxxxi,j=1，2，......，n从而可得：TnnjinnnzzzxTTdiagaaaazA1...)2(...21...)(1,011这样就得到了)(0zG分子分母的各系数的表达式：AHaaaadddnnn01110......（4.5）其中：nTTdiagH)2(...21jixA,55理论计算)2)(1(4)(ssssGssssG234)(23（5.1）将112zzTs代入（5.1）式得上下同时乘以31）（z将A(z)，B(z)写成向量形式：由于，，2)1)(1(zz，均为3阶多项式，可得到：这样就得到了G(z)的分子、分母的各系数的表达式：61331111111111331)2()2(210004323210TTTdddd32323232812424124241248124TTTTTTTTTTTT采样周期取T=1时，0832243210eeee6结果分析分析离散化步长对系统性能的影响：由于对于连续系统离散化的问题，采样周期和步长一样。所以在这里仅分析T对系统性能的影响。以为例，分别，通过MATLAB程序求其性能指标pt，rt，st和超调量pos。表7.1不同的T对系统性能的影响性能指标strtpospt连续系统37.55151.745870.01612.9975T=0.232.8000269.91093.2000T=0.5392.500066.21273.5000T=148362.55734T=52551029.714910T=107002014.713740T=1002001.555726007(a)连续系统8(b)T=0.2(c)T=0.59(d)T=1(e)T=510(f)T=10(g)T=100图7.1连续系统及T取不同值时单位阶跃响应曲线由以上结果可知，当T的取值越大时，系统的性能指标，st，rt都随之增大，超调量随之减小，超过某一值时，系统不在稳定；T的取值越小，与原系统的性能指标越接近。参考文献[1]孙增晰.计算机辅助设计.北京：清华大学出版社，1995[2]郑大钟.线性系统理论.北京：清华大学出版社，2004[3]姚全珠.软件技术基础.北京：电子工业出版社，2006[4]王小丹.基于MATLAB的系统分析与设计.西安：西安电子科技大学出版社，200711附录1程序清单1利用双线性变换求其离散传递函数的计算机实现num=input('请输入分子多项式num=');den=input('请输入分母多项式den=');T=input('请输入步长T=');sys=tf(num,den)n=length(den)-1%计算传递系统阶数fori=1:(n+1)/2d=num(i);num(i)=num(n-i+2);num(n-i+2)=d;e=den(i);den(i)=den(n-i+2);den(n-i+2)=e;end%颠倒原系统分子分母向量顺序disp(num)disp(den)a(1)=1;b=[1,1];c=[1,1];forj=2:n+1;a(j)=a(j-1)*(2/T);end%构造内部向量以1为首项2/T为倍数的等比序列H=diag(a)%显示对角矩阵%求出z+1的n次方的序列的系数矩阵fori=1:n-1;b=conv(b,c);endfori=1:n+1;A(i,1)=1;endforj=1:n+1;A(1,j)=b(j);endfori=2:n+1;forj=2:n+1;A(i,j)=A(i-1,j)-A(i,j-1)-A(i-1,j-1);endenddisp(A)%显示矩阵A12dnum=num*H*A%得到离散后的分子系数矩阵dden=den*H*A%得到离散后的分母系数矩阵sys=tf(dnum,dden,T)运行程序，显示如下：tustin请输入分子多项式num=[0004]请输入分母多项式den=[1320]请输入步长T=[1]sys=4-----------------s^3+3s^2+2sContinuous-timetransferfunction.n=340000231H=100002001300400008133111-1-11-1-111-33-1dnum