实验二 控制系统的数学模型及其转换.doc

实验二 控制系统的数学模型及其转换 【实验2.1】求传递函数的分子和分母多项式，并求传递函数的特征根。 解：运行实验shiyan2_1.m代码 r = 0 0 -3.0000 -1.6506 -0.1747 + 1.5469i -0.1747 - 1.5469i 【实验2.2】某一以微分方程描述系统的传递函数，其微分方程描述如下： ，试用MATLAB建立其模型。 解：运行实验shiyan2_2.m代码 num/den = s^2 + 4 s + 8 ------------------------ s^3 + 11 s^2 + 11 s + 10 【实验2.3】已知某系统的传递函数，求其部分分式表示形式，并求其状态空间模型。 解：运行shiyan2_3.m代码，根据r、p、k结果可得部分分式形式为 r = 0.0000 - 0.2500i 0.0000 + 0.2500i -2.0000 p = -0.0000 + 2.0000i -0.0000 - 2.0000i -1.0000 k = 2 A = -1 -4 -4 1 0 0 0 1 0 B = 1 0 0 C = -2 1 -7 D = 2 【实验2.4】给定RLC网络如图2-1所示，其中为输入量，为输出量。求解这个系统的传递函数模型、零极点增益模型以及状态空间模型（假设），并求取系统的阶跃响应。 图2-1 RLC网络 解：从数学上求出系统的传递函数。 根据电路基本定理，列出该电路的微分方程： 同时还有 整理以上方程，并在零初始条件下，取拉普拉斯变换，可得： 代入具体数值可得 运行shiyan2_4.m代码。 Zero/pole/gain: 1 --------------- (s^2 + 2s + 2) a = x1 x2 x1 -2.00000 -1.41421 x2 1.41421 0 b = u1 x1 1.00000 x2 0 c = x1 x2 y1 0 0.70711 d = u1 y1 0 【练习2.1】控制系统结构如图4-2所示 图2-2 （1） 利用MATLAB函数代码对以上单位负反馈系统建立传递函数模型； （2） 将第一问中求得的传递函数模型转化为零极点增益形式和状态空间形式。 num1=[2];den1=[1]; num2=[1 1];den2=[1 2 1]; [num3,den3]=series(num1,den1,num2,den2); [num,den]=cloop(num3,den3,-1); printsys(num,den) [A,B,C,D]=tf2ss(num,den) [Z,P,K]=ss2zp(A,B,C,D) num/den = 2 s + 2 ------------- s^2 + 4 s + 3 A = -4 -3 1 0 B = 1 0 C = 2 2 D = 0  Z = -1.0000 P = -3 -1 K = 2.0000【练习2.2】某控制系统结构如图2-3所示 图2-3 利用MATLAB函数代码对以上系统建立传递函数模型； num1=[2];den1=[1];num2=[1];den2=[1 1]; num3=[1 1];den3=[1 2 1]; num4=[1 1];den4=[1 0];num5=[2];den5=[1]; [num6,den6]=parallel(num3,den3,num4,den4); [na,da]=series(num1,den