课 程 设 计 报 告学生姓名:学 号：学 院:班 级:题 目:专业方向课程设计仿真组-利用双线性变换求其离散传递函数指导教师： 孟杰 姜文娟 职称: 副教授 讲师 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 计算机实现程序框图171.题目背景与意义本课程设计以自动控制理论、现代控制理论、MATLAB及应用等知识为基础，利用双线性变换求连续系统对应的离散化的系统，目的是使学生在现有的控制理论的基础上，学会用MATLAB语言编写控制系统的离散化的程序，通过上机实习加深对课堂所学知识的理解，掌握一种能方便地对系统进行离散化的设计工具。2.设计题目介绍2.1 设计内容和要求1 在理论上对连续系统采用双线性变换求离散化推导出算法和计算公式2 画出计算机实现算法的框图3 编写程序并调试和运行4 以下面的系统为例，进行计算已知系统闭环传递函数，利用双线性变换求其离散传递函数。5 分析运算结果（离散化步长对系统性能的影响）6 程序应具有一定的通用性，对不同参数能有兼容性。2.2设计工作任务及工作量的要求1 本次课程设计要求每周学生至少见指导教师2次，其中集中辅导答疑部不于3次。2 设计说明书的格式按设计说明书格式要求，采用word软件排版，计算机打印。（具体包括：封皮、目录、正文、参考文献等）3 程序清单用A4纸打印后，作为附录订装在说明书后面。4 框图和其他图表放在正文中。3 双线性变换法3.1 双线性变换的定义双线性变换法又称突斯汀（Tustin）法，是一种基于梯形积分规则的数字积分变换方法。3.2 双线性变换法的优缺点3.2.1 双线性变换法的优点双线性变换的主要优点:靠频率的严重非线性关系得到S平面与Z平面的单值一一对应关系，整个j轴单值对应于单位圆一周，其中和为非线性关系。在零频率附近，接近于线性关系，进一步增加时，增长变得缓慢，(终止于折叠频率处)，所以双线性变换不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象。3.2.2 双线性变换法的缺点双线性变换法的缺点:与的非线性关系，导致数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变，(使数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变)。例如，一个模拟微分器，它的幅度与频率是线性关系，但通过双线性变换后，就不可能得到数字微分器。3.2 双线性变换的原理双线性变换公式可以从梯形积分公式中直接推导出来。按这种替换公式进行替换，可以保证的稳定性，而且，具有一定的仿真精度。已知梯形积分公式为： 即： 则有： 即： （3.1）（3.1）式称为双线性变换公式，也可写成为： （3.2） 3.3双线性变换的主要特性1、 若稳定，则一定稳定；2、 变换前后，稳态增益不变；3、 双线性变换后的阶次不变，且分子、分母具有相同的阶次。4 设计步骤双线性替换不仅可以精度极高的仿真模型，而且能利用计算机程序来实现这种替换。下面来介绍一种程序替换法。折线形系统的传递函数为 （4.1）在双线性替换下得到的Z传递函数为：； （4.2）现需要由（i=0，1，.，n）确定（i=0，1，.，n），若直接将双线性替换公式代入，可得： （4.3）将其分子、分母同时乘以，可得： （4.4）将，写成向量形式， 由于，.，均为n阶多项式，可得到：矩阵为(n+1)（n+1）阶；其中第一行诸元素为的展开式的各系数，第一列诸元素为1，阶次n确定后，这些元素均为已知，并可以证明，其余nn个元素可由下式求得： i,j=1，2，.，n从而可得： 这样就得到了分子分母的各系数的表达式： （4.5）其中： 5 理论计算 （5.1）将代入（5.1）式得上下同时乘以将A(z)，B(z)写成向量形式：由于，均为3阶多项式，可得到：这样就得到了G(z)的分子、分母的各系数的表达式：采样周期取T=1时，6 结果分析分析离散化步长对系统性能的影响：由于对于连续系统离散化的问题，采样周期和步长一样。所以在这里仅分析T对系统性能的影响。以为例，分别，通过MATLAB程序求其性能指标，和超调量pos。表7.1 不同的T对系统性能的影响性能指标连续系统37.55151.745870.01612.9975T=0.232.8000269.91093.2000T=0.5392.500066.21273.5000T=148362.55734T=52551029.714910T=107002014.713740T=1002001.55572600(a) 连续系统(b) T=0.2(c) T=0.5(d) T=1(e) T=5(f) T=10(g) T=100图7.1 连续系统及T取不同值时单位阶跃响应曲线由以上结果可知，当T的取值越大时，系统的性能指标，都随之增大，超调量随之减小，超过某一值时，系统不在稳定；T的取值越小，与原系统的性能指标越接近。参考文献1 孙增晰.计算机辅助设计.北京：清华大学出版社，19952 郑大钟.线性系统理论.北京：清华大学出版社，20043 姚全珠.软件技术基础.北京：电子工业出版社，20064 王小丹.基于MATLAB的系统分析与设计.西安：西安电子科技大学出版社，2007附录1 程序清单1 利用双线性变换求其离散传递函数的计算机实现num=input(请输入分子多项式num=);den=input(请输入分母多项式den=);T=input(请输入步长T=);sys=tf(num,den)n=length(den)-1 %计算传递系统阶数for i=1:(n+1)/2 d=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;for j=2:n+1; a(j)=a(j-1)*(2/T);end %构造内部向量以1为首项2/T为倍数的等比序列H=diag(a) %显示对角矩阵%求出z+1的n次方的序列的系数矩阵for i=1:n-1; b=conv(b,c);endfor i=1:n+1; A(i,1)=1;endfor j=1:n+1; A(1,j)=b(j);endfor i=2:n+1; for j=2:n+1; A(i,j)=A(i-1,j)-A(i,j-1)-A(i-1,j-1); endenddisp(A) %显示矩阵Adnum=num*H*A %得到离散后的分子系数矩阵dden=den*H*A %得到离散后的分母系数矩阵sys=tf(dnum,dden,T) 运行程序，显示如下： tustin请输入分子多项式num=0 0 0 4请输入分母多项式den=1 3 2 0请输入步长T=1sys = 4 - s3 + 3 s2 + 2 s Continuous-time transfer function.n = 3 4 0 0 0 0 2 3 1H = 1 0 0 0 0 2 0 0 0 0 4 0 0 0 0 8 1 3 3 1 1 1 -1 -1 1 -1 -1 1 1 -3 3 -1dnum = 4 12 12 4dden = 24 -32 8 0sys = 4 z3 + 12 z2 + 12 z + 4 - 24 z3 - 32 z2 + 8 z Sample time: 1 secondsDiscrete-time transfer function.2 分析离散化步长对系统性能的影响的计算机实现num=input(请输入分子多项式num=);den=input(请输入分母多项式den=);G0=tf(num,den);bnum,bden=cloop(num,den,-1); %计算闭环系统的分子分母多项式sys=tf(bn