离散控制系统的数学描述与建模-计算机控制实验一(共7页)

1、精选优质文档-倾情为你奉上实验一 离散控制系统的数学描述与建模一、实验目的 学会使用Matlab建立离散控制系统的数学模型。二、实验指导1 传递函数模型传递函数模型由分子分母多项式来表达。传递函数模型分为SISO（单输入单输出）模型和MIMO（多输入多输出）模型。在这里我们只讨论SISO模型。离散时间SISO系统的脉冲传递函数模型为创建SISO传递函数模型，常用的方法是使用tf命令，用法是sys= tf (num，den)或sysd= tf (num，den ，Ts)式中的num和den分别是分子多项式和分母多项式的系数构成的向量，在向量中，系数按变量s的降幂排列；Ts为系统的采样周期。sys

2、为线性连续系统模型；sysd为线性离散系统模型。2零极点增益模型零极点模型实际上是传递函数模型的另一种表现形式，其原理是分别对原系统传递函数的分子分母进行分解因式处理，以获得系统的零极点表示形式。离散时间SISO系统的零极点模型的一般形式为创建SISO系统的零极点模型常用的方法是使用zpk命令，用法是sys = zpk ( z，p，k )或sysd = zpk ( z，p，k，Ts )式中的z，p，k分别是系统的零极点和增益向量。系统的零极点模型可以被直接用来判断系统的稳定性。3 状态空间模型状态空间模型是采用线性微分或差分方程来描述系统的动态行为。离散系统状态空间模型的一般形式为使用ss命令

3、可以创建系统的状态空间模型，其命令的调用格式为sys=ss ( A，B，C，D)或sysd=ss ( A，B，C，D，Ts)4 数学模型的转换各种数学模型适用于各类不同的场合，因而当研究的范围发生变化时，就需要对原有的数学模型进行转换，以适应工程实际的需要。Matlab提供了许多可以对同一控制系统的模型描述进行转换的函数，其中常用函数如下表10.5所示。表10.5 模型转换函数及说明函数说明tf2ss由传递函数模型转化为状态空间模型tf2zp由传递函数模型转化为零极点模型zp2ss由零极点模型转化为状态空间模型zp2tf由零极点模型转化为传递函数模型ss2tf由状态空间模型转化为传递函数模型s

4、s2zp由状态空间模型转化为零极点模型5 控制系统模型的典型连接控制系统是由多个环节组成的，每个环节又是由多个元件构成的。在控制系统设计中，有三种典型的连接方式。1） 串联环节如果有两个环节sys1和sys2串联，则其等效传递函数为sys=series(sys1,sys2)实际上，sys=series(sys1,sys2)命令现在很少用，它已由命令sys=sys1*sys2*sysn所取代，这样命令中不仅省掉了“series()”字符，而且可以实现多个传递函数模快的串联。2） 并联环节如果有两个环节sys1和sys2并联，则其等效传递函数为sys=parallel(sys1,sys2)或sys

5、=sys1+sys2+sysn3） 反馈连接反馈环节的连接，其等效传递函数可由命令G=feedback(G1,G2,sign)来计算，其中G1为闭环前向通道传递函数；G2位反馈通道传递函数；sign为反馈方式，sign=1为正反馈；sign=1或缺省，为负反馈系统6 脉冲传递函数的部分分式展开脉冲传递函数Matlab提供了函数residue()对两个多项式的比进行分部展开,其调用格式为r,p,k=residue(f,g) 。三、实验内容及结果1试用matlab语句编程，求下列脉冲传递函数的零点、极点、增益。 程序及结果： num=1,0,11,30,0den=1,9,45,87,10sysd=

6、 tf (num,den,1)z p k=tf2zp(num,den)num = 1 0 11 30 0den = 1 9 45 87 10 Transfer function: z4 + 11 z2 + 30 z-z4 + 9 z3 + 45 z2 + 87 z + 10 Sampling time: 1num = 1 0 11 30 0den = 1 9 45 87 10 Transfer function: z4 + 11 z2 + 30 z-z4 + 9 z3 + 45 z2 + 87 z + 10 Sampling time: 1 Zero/pole/gain: z (z+2) (z

7、2 - 2z + 15)-(z+3.157) (z+0.1225) (z2 + 5.72z + 25.85) Sampling time: 1z = 0 1.0000 + 3.7417i 1.0000 - 3.7417i -2.0000 p = -2.8602 + 4.2038i -2.8602 - 4.2038i -3.1572 -0.1225 k = 12已知二阶离散系统变换传递函数为：试求采样周期时系统的传递函数模型和零极点增益模型程序及结果：num=1.6,-5.8,3.9den=1,-0.7,1.4num = 1.6000 -5.8000 3.9000den = 1.0000 -0.

8、7000 1.4000 | >> sysd=tf(num,den,0.1) Transfer function:1.6 z2 - 5.8 z + 3.9- z2 - 0.7 z + 1.4 Sampling time: 0.1>> z p k=tf2zp(num,den)z = 2.7332 0.8918p = 0.3500 + 1.1303i 0.3500 - 1.1303ik = 1.6000 3 试用matlab语句编程将部分分式展开。程序及结果： >> num=10,0num = 10 0>> den=1,-3,2den = 1 -3 2

9、>> sysd=tf(num,den,1) Transfer function: 10 z-z2 - 3 z + 2 Sampling time: 1>> r,p,k=residue(num,den) r = 20 -10p = 2 1k = 4 试用matlab语句编程将部分分式展开。程序及结果：>> p=1,1,2p = 1 1 2>> z=0,1.5z = 0 1.5000>> k=2k = 2>> num,den=zp2tf(z,p,k) >> sysd = zpk ( z,p,k,1 ) Zero/pole/gain: 2 z (z-1.5)-(z-1)2 (z-2) Sampling time: 1>> num,den=zp2tf(z,p,k) >> num,den=zp2tf(z',p,k)num = 0