00讲稿 自动控制原理实验指导书 00实验一 基于MATLAB实验平台的系统被控对象的建立与转换
- 格式:doc
- 大小:177.50 KB
- 文档页数:10
1 《自动控制原理》 实验指导书 2
实验一 基于MATLAB实验平台的系统被控对象的建立与转换 [实验目的] 1.了解MATLAB软件的基本特点和功能; 2.掌握线性系统被控对象传递函数数学模型在MATLAB环境下的表示方法及转换;
3.掌握多环节串联、并联、反馈连接时整体传递函数的求取方法; 4.了解在MATLAB环境下求取系统的输出时域表达式的方法。
[实验指导] 一、被控对象模型的建立 在线性系统理论中,一般常用的描述系统的数学模型形式有: (1)传递函数模型——有理多项式分式表达式 (2)传递函数模型——零极点增益表达式 这些模型之间都有着内在的联系,可以相互进行转换。
1、传递函数模型——有理多项式分式表达式 设系统的传递函数模型为
01110111......)()()(asasasabsbsbsbsRsCsGnnnnmmmm
对线性定常系统,式中s的系数均为常数,且an不等于零。 这时系统在MATLAB中可以方便地由分子和分母各项系数构成的两个向量唯一地确定,这两个向量常用num和den表示。 num=[bm,bm-1,…,b1,b0] den=[an,an-1,…,a1,a0]
注意:它们都是按s的降幂进行排列的。分子应为m项,分母应为n项,若有空缺项(系数为零的项),在相应的位置补零。 然后写上传递函数模型建立函数:sys=tf(num,den)。这个传递函数便在MATLAB平台中被建立,并可以在屏幕上显示出来。 举例1-1:已知系统的传递函数描述如下:
22642202412)(23423sssssssG 3
在MATLAB命令窗口(Command Window)键入以下程序: >> num=[12,24,0,20]; >> den=[2 4 6 2 2]; >> sys=tf(num,den) 回车后显示结果: Transfer function: 12 s^3 + 24 s^2 + 20 --------------------------------------- 2 s^4 + 4 s^3 + 6 s^2 + 2 s + 2
并同时在MATLAB中建立了这个相应的有理多项式分式形式的传递函数模型。 举例1-2:已知系统的传递函数描述如下:
)523()1()66)(2(4)(23322sssssssssG
其中,多项式相乘项,可借助多项式乘法函数conv来处理。 在MATLAB命令窗口(Command Window)键入以下程序: >> num=4*conv([1,2],conv([1,6,6],[1,6,6])); >> den=conv([1,0],conv([1,1],conv([1,1],conv([1,1],[1,3,2,5])))); >> sys=tf(num,den) 回车后显示结果: Transfer function: 4 s^5 + 56 s^4 + 288 s^3 + 672 s^2 + 720 s + 288 --------------------------------------------------------------------------- s^7 + 6 s^6 + 14 s^5 + 21 s^4 + 24 s^3 + 17 s^2 + 5 s
即同时在MATLAB中建立了这个有理多项式分式形式的传递函数模型。
2.传递函数模型——零极点增益模型 零极点增益模型为:
))...()(())...()(()(2121nmpspspszszszsKsG
其中:K为零极点增益,zi为零点,pj为极点。 该模型 在MATLAB中,可用[z,p,k]矢量组表示,即 z=[z1,z2,…,zm]; 4
p=[p1,p2,...,pn]; k=[K]; 然后在MATLAB中写上零极点增益形式的传递函数模型建立函数: sys=zpk(z,p,k)。这个零极点增益模型便在MATLAB平台中被建立,并可以在屏幕上显示出来。 举例1-3:已知系统的零极点增益模型:
)5)(2)(1()3(6)(sssssG
在MATLAB命令窗口(Command Window)键入以下程序: >> z=[-3];p=[-1,-2,-5];k=6; >> sys=zpk(z,p,k) 回车后显示结果: Zero/pole/gain: 6 (s+3) ----------------- (s+1) (s+2) (s+5)
则在MATLAB中建立了这个零极点增益的模型。
二、不同形式模型之间的相互转换 不同形式之间模型转换的函数: (1)tf2zp: 多项式传递函数模型转换为零极点增益模型。 格式为:[z,p,k]=tf2zp(num,den) (2)zp2tf: 零极点增益模型转换为多项式传递函数模型。 格式为:[num,den]=zp2tf(z,p,k)
三、环节串联、并联、反馈连接时等效的整体传递函数的求取 1、串联
这里: 在MATLAB中求取整体传递函数的功能,采用如下的语句或函数来实现。. ① G=G1*G2 5
② G=series(G1,G2) ③ [num,den]=series(num1,den1,num2,den2) 举例1-4 两环节G1、G2串联,求等效的整体传递函数G
32)(1ssG 127)(22sssG 解:①实现的程序: >> n1=2;d1=[1 3];n2=7;d2=[1 2 1];G1=tf(n1,d1);G2=tf(n2,d2);G=G1*G2 运行结果: Transfer function: 14 --------------------- s^3 + 5 s^2 + 7 s + 3
②实现的程序: >>n1=2;d1=[1 3];n2=7;d2=[1 2 1];G1=tf(n1,d1);G2=tf(n2,d2);G=series(G1,G2) 运行结果: Transfer function: 14 --------------------- s^3 + 5 s^2 + 7 s + 3
③实现的程序: >>n1=2;d1=[1 3];n2=7;d2=[1 2 1];G1=tf(n1,d1);G2=tf(n2,d2); >> [n,m]=series(n1,d1,n2,d2)
运行结果: n = 0 0 0 14 m = 1 5 7 3
举例1-5 四环节G1、G2、 G3、G4串联,求等效的整体传递函数G
32321sGGG 12722ssG 解:实现的程序: >> n1=2;d1=[1 3];n2=7;d2=[1 2 1];G1=tf(n1,d1);G2=tf(n2,d2); >>G=G1*G2*G1*G1 运行结果: Transfer function: 6
56 ------------------------------------------ s^5 + 11 s^4 + 46 s^3 + 90 s^2 + 81 s + 27
2、并联
两环节G1(s)与G2
(s)并联,则等效的整体传递函数为
G(s)=G1(s)+G2(s) 在MATLAB中求取整体传递函数的功能,采用如下的语句或函数来实现。 ① G=G1+G2 ② G= parallel (G1,G2) ③ [num,den]= parallel (num1,den1,num2,den2) 举例1-6 两环节G1、G2并联,求等效的整体传递函数G(s)
32)(1ssG 127)(22sssG 解:①实现的程序: >> n1=2;d1=[1 3];n2=7;d2=[1 2 1];G1=tf(n1,d1);G2=tf(n2,d2);G1+G2 运行结果: Transfer function: 2 s^2 + 11 s + 23 ---------------------------- s^3 + 5 s^2 + 7 s + 3
②实现的程序: >> n1=2;d1=[1 3];n2=7;d2=[1 2 1];G1=tf(n1,d1);G2=tf(n2,d2);G=parallel(G1,G2) 运行结果: Transfer function: 2 s^2 + 11 s + 23 --------------------------- s^3 + 5 s^2 + 7 s + 3
③实现的程序: >> n1=2;d1=[1 3];n2=7;d2=[1 2 1]; [n,d]=parallel(n1,d1,n2,d2) 运行结果: n = 0 2 11 23