第六章系统模型建立与转换PPT讲稿.ppt

《第六章系统模型建立与转换PPT讲稿.ppt》由会员分享，可在线阅读，更多相关《第六章系统模型建立与转换PPT讲稿.ppt（46页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第六章系统模型建立与转换第1页，共46页，编辑于2022年，星期三6.1 6.1 控制系统的基本概念控制系统的基本概念输入输入偏差偏差反馈信号反馈信号控制系统要求有比较好的稳定性、准确性和动态特性。第2页，共46页，编辑于2022年，星期三6.2 6.2 反馈控制系统的数学模型反馈控制系统的数学模型为了对系统的性能进行分析首先要建立其数学模型为了对系统的性能进行分析首先要建立其数学模型数学模型：描述系统、输入和输出三者之间动态关系的数学数学模型：描述系统、输入和输出三者之间动态关系的数学表达式。表达式。微分方程微分方程、传递函数、频率特性、传递函数、频率特性第3页，共46页，编辑于2022年，

2、星期三 为为了了对对系系统统的的性性能能进进行行分分析析首首先先要要建建立立其其数数学学模模型型 ，在在MATLABMATLAB中提供了中提供了3 3种数学模型描述的形式：种数学模型描述的形式：（1 1）传递函数模型）传递函数模型tf(tf(）（2 2）零极点形式的数学模型）零极点形式的数学模型zpk()zpk()（3 3）状态空间模型）状态空间模型ss()ss()本本节节首首先先介介绍绍利利用用MATLABMATLAB提提供供的的3 3个个函函数数来来建建立立系系统统的的数数学学模型模型，然后在此基础上介绍各种数学模型之间的相互转换。，然后在此基础上介绍各种数学模型之间的相互转换。一、线性系

3、统的传递函数模型一、线性系统的传递函数模型第4页，共46页，编辑于2022年，星期三格式：格式：syssystftf（numnum，denden）功能：功能：建立系统的传递函数模型建立系统的传递函数模型说明：说明：numnum、den den分别为传递函数分子和分母多项式的系数向量分别为传递函数分子和分母多项式的系数向量1 1tf tf 传递函数模型传递函数模型第5页，共46页，编辑于2022年，星期三例：已知系统的传递函数为例：已知系统的传递函数为:试建立系统的传递函数模型。试建立系统的传递函数模型。num=29den=13246G=tf(num,den)第6页，共46页，编辑于2022年，

4、星期三例例 ：已知系统传递函数如下：已知系统传递函数如下应用应用MatlabMatlab语言建立系统的传递函数模型语言建立系统的传递函数模型若分子和分母不是完全展开形式若分子和分母不是完全展开形式,两个多项式的乘积两个多项式的乘积在在matlab中可以借用函数中可以借用函数conv()得出。得出。num=5*1,2.4den=conv(11,conv(11,conv(134,101)G=tf(num,den)第7页，共46页，编辑于2022年，星期三2 2线性系统的线性系统的零极点模型零极点模型格式：格式：syssyszpk(z,p,k)zpk(z,p,k)功能：功能：建立零极点形式的数学模型

5、建立零极点形式的数学模型说明：说明：系统的传递函数还可以表示成零极点形式，零极点模系统的传递函数还可以表示成零极点形式，零极点模 型一般表示为：型一般表示为：其中其中 ZiZi（i i1,21,2,m,m）和和 PiPi（i i1,21,2,n,n）分别为系统的零分别为系统的零点和极点，点和极点，K K为系统的增益。为系统的增益。z z、p p、k k分别为系统的零极点向量分别为系统的零极点向量和增益。和增益。第8页，共46页，编辑于2022年，星期三例例 ：已知系统传递函数如下：已知系统传递函数如下应用应用MatlabMatlab语言建立系统的语言建立系统的零极点形式零极点形式模型模型。p=

6、-1,-2,-3z=-4G=zpk(z,p,5)第9页，共46页，编辑于2022年，星期三3 3SS SS 状态空间模型状态空间模型格式：格式：syssysssss（A,B,C,DA,B,C,D）功能：功能：建立系统的状态空间模型建立系统的状态空间模型说明说明：状态方程是研究系统的最为有效的系统数学描述：状态方程是研究系统的最为有效的系统数学描述第10页，共46页，编辑于2022年，星期三6.3 系统的系统的组合和连接组合和连接所所谓谓系系统统组组合合，就就是是将将两两个个或或多多个个子子系系统统按按一一定定方方式式加加以以连连接接形形成成新新的的系系统统。这这种种连连接接组组合合方方式式主主

7、要要有有串串联联、并并联联、反反馈馈等等形形式式。MATLABMATLAB提供了进行这类组合连接的相关函数。提供了进行这类组合连接的相关函数。1.series 1.series 系统的串联系统的串联格式格式1 1：syssysseriesseries（sys1,sys2sys1,sys2），），功能：功能：用于将两个线性模型串联形成新的系统即用于将两个线性模型串联形成新的系统即syssyssys1*sys2 sys1*sys2 说明：说明：格式格式1 1：对应于：对应于SISOSISO系统的串联连接。系统的串联连接。)(1sG)(2sGU(s)Y(s)G=series(G1,G2)G=G1*G

8、2第11页，共46页，编辑于2022年，星期三2 2parallelparallel格式格式1 1：sys=parallel(sys1,sys2)sys=parallel(sys1,sys2)功能：功能：将两个系统以并联方式连接成新的系统，将两个系统以并联方式连接成新的系统，即即sys=sys1+sys2sys=sys1+sys2。说明：说明：并联连接时，输入信号相同，并联后其输出为并联连接时，输入信号相同，并联后其输出为sys1sys1和和 sys2 sys2这两个系统的输出之和。这两个系统的输出之和。所以总的传递函数为所以总的传递函数为G(s)=GG(s)=G1 1(s)+G(s)+G2

9、2(s)(s)。G=parallel(G1,G2)G=G1+G2第12页，共46页，编辑于2022年，星期三例例 已知两个线性系统已知两个线性系统 ，分别应用分别应用seriesseries和和parallelparallel函数进行系统的串并联连接。函数进行系统的串并联连接。第13页，共46页，编辑于2022年，星期三3 3feedback feedback 系统的反馈连接。系统的反馈连接。格式格式1 1：sys=feedback(sys1,sys2,sign)sys=feedback(sys1,sys2,sign)功能：功能：实现两个系统的反馈连接。实现两个系统的反馈连接。说明：说明：对于

10、对于SISOSISO系统，系统，sys1sys1表示前向通道传函，表示前向通道传函，sys2 sys2表示反馈通道表示反馈通道,sign=1,sign=1,正反馈正反馈.sign=-1,sign=-1,负反馈负反馈(默认值默认值,可省略可省略)（a）xo(s)xi(s)G(s)H(s)E(s)B(s)GB=feedback(G,H)GB=feedback(G,H,1)第14页，共46页，编辑于2022年，星期三例例 ：如图所示控制系统：如图所示控制系统求其总的传递函数求其总的传递函数第15页，共46页，编辑于2022年，星期三例例 ：如图所示控制系统，已知各个子传递函数如下：如图所示控制系统，

11、已知各个子传递函数如下求其总的传递函数求其总的传递函数G=tf(172424,110355024)Gc=tf(106,10)H=1GG=feedback(G*Gc,H)H(s)=1第16页，共46页，编辑于2022年，星期三 在进行系统分析时，往往根据不同的要求选择不同形式的数学模型，因在进行系统分析时，往往根据不同的要求选择不同形式的数学模型，因此经常要在不同形式数学模型之间相互转换，下面介绍三种模型之间的相互此经常要在不同形式数学模型之间相互转换，下面介绍三种模型之间的相互转换函数。转换函数。6.4 6.4 模型的转换模型的转换第17页，共46页，编辑于2022年，星期三1）tf模型模型z

12、pk模型，调用函数格式：模型，调用函数格式：zpk（sys）或）或tf2zp（sys）2）tf模型模型ss模型，调用函数格式：模型，调用函数格式：ss(sys)或或tf2ss（sys）3）zpk模型模型tf模型，调用函数格式：模型，调用函数格式：tf(sys)或或zp2tf(sys)4）zpk模型模型ss模型，调用函数格式：模型，调用函数格式：ss（sys）或）或zp2ss（sys）5）ss模型模型tf模型，调用函数格式：模型，调用函数格式：tf(sys)或或ss2tf(sys)6）ss模型模型zpk模型，调用函数格式：模型，调用函数格式：zpk(sys)或或ss2zp(sys)第18页，共4

13、6页，编辑于2022年，星期三应应用用MATLABMATLAB的的模模型型转转换换函函数数将将其其转转换换为为零零极极点点形形式式的的模模型型和和状状态态空空间间模型模型例：已知系统的传递函数为例：已知系统的传递函数为num=21840den=16116G=tf(num,den)GZPK=ZPK(G)%GZPK=tf2zp(G)GS=ss(G)%GS=tf2ss(G)第19页，共46页，编辑于2022年，星期三第七章第七章第七章第七章 控制系统的计算机辅控制系统的计算机辅控制系统的计算机辅控制系统的计算机辅助分析助分析助分析助分析7.1 7.1 系统的时域分析系统的时域分析7.2 7.2 系统

14、的频域分析系统的频域分析7.3 7.3 系统的稳定性分析系统的稳定性分析本章介绍进行系统分析和获得系统各项性能指标的方法本章介绍进行系统分析和获得系统各项性能指标的方法,通过本通过本章的学习，掌握利用章的学习，掌握利用MATLAB及其控制系统工具箱对各种控制及其控制系统工具箱对各种控制系统进行时域及频域分析的方法。系统进行时域及频域分析的方法。第20页，共46页，编辑于2022年，星期三7.1时域分析的一般方法1、拉氏逆变换法例1：已知某系统的微分方程如下，其中T=0.5s，k=10，设初始条件为零，求系统单位阶跃响应曲线。Symssy=ilaplace(10/(0.5*s+1)*s)Ezpl

15、ot(y,05)grid第21页，共46页，编辑于2022年，星期三2、直接计算法、直接计算法1 1）impulse impulse 求连续系统的单位脉冲响应。求连续系统的单位脉冲响应。格式：格式：impulse(sys)Y,X=impulse(sys)impulse(sys)Y,X=impulse(sys)说说 明明：syssys为为tf(),zpk(),ss()tf(),zpk(),ss()中中任任一一种种模模型型，对对于于不不带带返返回回参参数数的的该该函函数数在在当当前前窗窗口口中中绘绘制制出出响响应应曲曲线线。对对于于带带有有返返回回参参数数的将不绘制曲线。的将不绘制曲线。第22页，

16、共46页，编辑于2022年，星期三例例2 2：如图所示的典型反馈控制系统结构，：如图所示的典型反馈控制系统结构，H(S)H(S)G(S)G(S)Gc(S)Gc(S)其中，其中，求系统的开环和闭环单位脉冲响应。求系统的开环和闭环单位脉冲响应。_ _G=tf(4,1234)Gc=tf(1-3,13)H=tf(1,0.011)G1=G*GcGk=G1*HGB=feedback(G1,H)Impulse(Gk)pauseImpulse(GB)第23页，共46页，编辑于2022年，星期三2 2）step step 求连续系统的单位阶跃响应。求连续系统的单位阶跃响应。格式格式1 1：step(sys)Y,

17、X=step(sys)step(sys)Y,X=step(sys)说说 明明：step()step()中中的的参参数数意意义义和和implse()implse()函函数数相相同同。如如果果用用户户在在调调用用step()step()函函数数时时不不返返回回任任何何向向量量，则则将将自自动动地地绘绘出出阶阶跃跃响响应输出曲线。应输出曲线。例例3:3:考虑下面传递函数模型：考虑下面传递函数模型：试绘制其单位阶跃响应曲线。试绘制其单位阶跃响应曲线。第24页，共46页，编辑于2022年，星期三 1 1、假设将自然频率固定为、假设将自然频率固定为 1 1，0,0.1,0.5,0.7,1,50,0.1,0

18、.5,0.7,1,5。例例4：典型二阶系统传递函数为：典型二阶系统传递函数为：试分析不同参数下的系统单位阶跃响应试分析不同参数下的系统单位阶跃响应2 2、将阻尼比、将阻尼比的值固定在的值固定在0.550.55，自然频率自然频率 变化范变化范 围为围为0.1-1 0.1-1 wn=1forkic=0,0.1,0.5,0.7,1,5G=tf(wn*wn,12*kic*wnwn*wn)Step(G)holdonendkic=0.55forwn=0.1:0.1:1G=tf(wn*wn,12*kic*wnwn*wn)Step(G)holdonend第25页，共46页，编辑于2022年，星期三 3)lsi

19、m 3)lsim 求任意输入信号时系统的响应求任意输入信号时系统的响应 格式格式1 1：lsim(sys,u,t)Y,X=lsim(sys,u,t)lsim(sys,u,t)Y,X=lsim(sys,u,t)说明：说明：sys sys 为为tf()tf()或或zpk()zpk()模型模型。例例6 6：系统传递函数为：系统传递函数为：对任意输入信号对任意输入信号u(t)求系统的输出响应曲线。求系统的输出响应曲线。g=tf(55,14230)t=0:0.1:10u=sin(t+30/180*pi)Lsim(g,u,t)第26页，共46页，编辑于2022年，星期三系统的性能指标系统的性能指标性能指标

20、：评价系统性能好坏的标准，它以准确的定量的方式来描述系统的性能。第27页，共46页，编辑于2022年，星期三时间时间tr上上升升峰值时间峰值时间tpAB超调量超调量Mp%=AB100%调节时间调节时间ts动态性能指标定义动态性能指标定义第28页，共46页，编辑于2022年，星期三1、编程用公式计算2、利用LTI Viewer获得响应曲线和性能指标系统的性能指标的计算方法第29页，共46页，编辑于2022年，星期三7.27.2利用利用LTI ViewerLTI Viewer获得响应曲线和性能指标获得响应曲线和性能指标在MATLAB中提供了线性时不变系统仿真的图形化工具LTIViewer，可方便地

21、获得系统的各种响应曲线和性能指标。使用方法：ltiview第30页，共46页，编辑于2022年，星期三例7：利用LTIViewer对典型二阶系统进行分析。假设固有频率wn 10，0.1,0.3,0.7,1,3。（1）在工作空间中建立一组二阶系统的传递函数。（2）打开LTIViewer，对此系统进行分析。wn=10;s1=tf(wn*wn,12*0.1*wnwn*wn)s2=tf(wn*wn,12*0.3*wnwn*wn)s3=tf(wn*wn,12*0.7*wnwn*wn)s4=tf(wn*wn,12*1*wnwn*wn)s5=tf(wn*wn,12*3*wnwn*wn)ltiview第31页

22、，共46页，编辑于2022年，星期三例例8.已知系统开环传函为：已知系统开环传函为：试绘制该单位负反馈系统的单位响应曲线，并计算系统的性试绘制该单位负反馈系统的单位响应曲线，并计算系统的性能指标。能指标。解：解：g=tf(1.25,110)g1=feedback(g,1)ltiview第32页，共46页，编辑于2022年，星期三7.3 7.3 系统的频域响应分析系统的频域响应分析时域分析时域分析时域分析时域分析：以时间以时间以时间以时间t t为独立变量，通过阶跃或脉冲输入下系统的瞬态为独立变量，通过阶跃或脉冲输入下系统的瞬态为独立变量，通过阶跃或脉冲输入下系统的瞬态为独立变量，通过阶跃或脉冲输

23、入下系统的瞬态时间响应来研究系统的性能；依据的数学模型为时间响应来研究系统的性能；依据的数学模型为时间响应来研究系统的性能；依据的数学模型为时间响应来研究系统的性能；依据的数学模型为G(s)G(s)。频域分析频域分析频域分析频域分析：以频率以频率以频率以频率 为独立变量，通过系统在不同频率的谐波（正弦）为独立变量，通过系统在不同频率的谐波（正弦）为独立变量，通过系统在不同频率的谐波（正弦）为独立变量，通过系统在不同频率的谐波（正弦）输入作用下的稳态响应来研究系统的性能；依据的输入作用下的稳态响应来研究系统的性能；依据的输入作用下的稳态响应来研究系统的性能；依据的输入作用下的稳态响应来研究系统的

24、性能；依据的数学模型为数学模型为数学模型为数学模型为G(j)G(j)。频域分析法是利用系统开环的奈氏图、波特图分析系统的性能，如系频域分析法是利用系统开环的奈氏图、波特图分析系统的性能，如系统的动态性能、稳定性统的动态性能、稳定性第33页，共46页，编辑于2022年，星期三一、频率特性曲线频率特性的极坐标图频率特性的极坐标图(Nyquist图图)频率特性的对数坐标图频率特性的对数坐标图(Bode图图)第34页，共46页，编辑于2022年，星期三1 1、nyquist nyquist 求连续系统的求连续系统的NyquistNyquist曲线曲线格式格式1 1：nyquist(sys)re,im,

25、w=nyquist(sys)nyquist(sys)re,im,w=nyquist(sys)格式格式2 2：nyquist(sys,w)re,im,w=nyquist(sys,w)nyquist(sys,w)re,im,w=nyquist(sys,w)说明：说明：syssys为为tf(),zpk(),ss()tf(),zpk(),ss()中任一种模型中任一种模型。w w设定频率范围设定频率范围 省略时由机器自动产生。省略时由机器自动产生。对于不带返回参数的将对于不带返回参数的将绘制绘制NyquistNyquist曲线曲线。对于带有返回参数的将不对于带有返回参数的将不绘制绘制曲线曲线g=tf(1

26、,11)Nyquist(g)第35页，共46页，编辑于2022年，星期三2 2、bode bode 求连续系统的求连续系统的BodeBode（伯德）频率响应。（伯德）频率响应。格式格式1 1：bode(sys)mag,phase,w=bode(sys)bode(sys)mag,phase,w=bode(sys)说明：说明：bodebode函数的输入变量定义与函数的输入变量定义与nyquistnyquist相同相同BodeBode图图可可用用于于分分析析系系统统的的增增益益裕裕度度、相相位位裕裕度度、增增益益、带带宽宽以以及及稳稳定定性性等等特性。特性。magmag和和phasephase分别是

27、幅值和相位数组。分别是幅值和相位数组。Bode(g)w=logspace(-1,3,100)Bode(g)第36页，共46页，编辑于2022年，星期三两个性能指标:相角裕度-反映系统的相对稳定性截止频率c-反映系统的快速性.-110jrc二、频率响应分析二、频率响应分析第37页，共46页，编辑于2022年，星期三margin margin 求取给定线性定常系统的幅值裕量和相角的裕量。求取给定线性定常系统的幅值裕量和相角的裕量。格式格式1 1：marginmargin（syssys）格式格式2 2：Gm,Pm,Wcg,Wcp=marginGm,Pm,Wcg,Wcp=margin（syssys）说

28、说明明：marginmargin函函数数可可从从频频率率响响应应数数据据中中计计算算出出幅幅值值裕裕度度、相相角角裕裕度度和剪切频率。和剪切频率。格式格式1 1 画出画出bodebode图，并标注幅值裕度和对应频率，相角裕图，并标注幅值裕度和对应频率，相角裕 度和对应频率。度和对应频率。格格式式2 2：不不画画图图，返返回回幅幅值值裕裕度度GmGm和和对对应应频频率率Wcg Wcg,相相角角裕裕度度PmPm和和对对应应频频率率Wcp Wcp。g=tf(1,10.21)margin(g)gmpmwcgwcp=margin(g)第38页，共46页，编辑于2022年，星期三线性定常系统稳定的充要条件

29、线性定常系统稳定的充要条件线性定常系统稳定的充要条件线性定常系统稳定的充要条件：若系统的全部特征根（传递函数的全部极点）均具若系统的全部特征根（传递函数的全部极点）均具若系统的全部特征根（传递函数的全部极点）均具若系统的全部特征根（传递函数的全部极点）均具有有有有负实部负实部负实部负实部（位于（位于（位于（位于ssss平面的左半平面），则系统稳定。平面的左半平面），则系统稳定。平面的左半平面），则系统稳定。平面的左半平面），则系统稳定。7.4 系统的稳定性分析第39页，共46页，编辑于2022年，星期三1.直接判定方法直接判定方法首首先先求求出出系系统统的的所所有有极极点点，当当其其实实部部大

30、大于于零零，则则系系统统为为不不稳稳定定系系统统，否否则称为稳定系统。若极点的实部等于则称为稳定系统。若极点的实部等于0的，则系统称为临界稳定系统。的，则系统称为临界稳定系统。对于传递函数模型对于传递函数模型tf（num，den）,利用求根函数利用求根函数roots（den）来求极点。也可利用零极点图来判断系统的稳定性）来求极点。也可利用零极点图来判断系统的稳定性例：例：假设系统的传递函数模型为假设系统的传递函数模型为判断系统的稳定性。判断系统的稳定性。p=110355024roots(p)第40页，共46页，编辑于2022年，星期三 2 2 稳定性判据稳定性判据 直接求取高阶代数方程的根是一

31、件很困难的工作，所以直接求取高阶代数方程的根是一件很困难的工作，所以出现了一些判定给定系统稳定性的间接方法，出现了一些判定给定系统稳定性的间接方法，Nyquist稳定判据当当当当 由由由由 到到到到+时时时时，若若若若系系系系统统统统的的的的开开开开环环环环奈奈奈奈魁魁魁魁斯斯斯斯特特特特频频频频率率率率特特特特性性性性曲曲曲曲线线线线逆时针方向包围（逆时针方向包围（逆时针方向包围（逆时针方向包围（1 1 1 1，j0j0j0j0）点）点）点）点P P P P圈，则闭环系统稳定。圈，则闭环系统稳定。圈，则闭环系统稳定。圈，则闭环系统稳定。第41页，共46页，编辑于2022年，星期三 例例:系统

32、开环传递函数系统开环传递函数绘制系统的绘制系统的NyquistNyquist图，并讨论其稳定性图，并讨论其稳定性NyquistNyquist图图顺顺时时针针包包围围（1 1，j0j0）点点2 2次次，而而原原开开环环系系统统中中没没有有不不稳稳定定极极点点，从而可以得出结论，闭环系统有从而可以得出结论，闭环系统有2 2个不稳定极点。个不稳定极点。G=tf(1000,181510)Nyquist(G)第42页，共46页，编辑于2022年，星期三第43页，共46页，编辑于2022年，星期三例例:系统开环传递函数系统开环传递函数绘制系统的绘制系统的NyquistNyquist图确定图确定K=-0.5K=-0.5和和k=-2k=-2时闭环系统的稳定时闭环系统的稳定性性第44页，共46页，编辑于2022年，星期三g=tf(-0.5,11)G1=tf(-2,11)第45页，共46页，编辑于2022年，星期三End第46页，共46页，编辑于2022年，星期三