




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
第一章 控制系统的状态空间表达式1.状态空间表达式x=Ax+Bun阶 w:rxly.mxlA:nxnB:nxrC:mxnD:mxry=Cx+DuA称为系统矩阵,描述系统内部状态之间的联系;B为输入(或控制)矩阵,表示输入对每个状态变量的作用情况;C输出矩阵,表示输出与每个状态变量间的组成关系,D直接传递矩阵,表示输入对输出的直接传递关系。2.状态空间描述的特点①考虑了“输入一状态一输出”这一过程,它揭示了问题的本质,即输入引起了状态的变化,而状态决定了输出。②状态方程和输出方程都是运动方程。③状态变量个数等于系统包含的独立贮能元件的个数,n阶系统有n个状态变量可以选择。④状态变量的选择不唯一。⑤从便于控制系统的构成来说,把状态变量选为可测量或可观察的量更为合适。⑥建立状态空间描述的步骤:a选择状态变量;b列写微分方程并化为状态变量的一阶微分方程组:c将一阶微分方程组化为向量矩阵形式,即为状态空间描述。⑦状态空间分析法是时域内的一种矩阵运算方法,特别适合于用计算机计算。.模拟结构图(积分器加法器比例器)已知状态空间描述,绘制模拟结构图的步骤:积分器的数目应等于状态变量数,将他们画在适当的位置,每个积分器的输出表示相应的某个状态变量,然后根据状态空间表达式画出相应的加法器和比例器,最后用箭头将这些元件连接起来。.状态空间表达式的建立①由系统框图建立状态空间表达式:a将各个环节(放大、积分、惯性等)变成相应的模拟结构图:b每个积分器的输出选作天,输入则为X;c由模拟图写出状态方程和输出方程。②由系统的机理出发建立状态空间表达式:如电路系统。通常选电容上的电压和电感上的电流作为状态变量。利用KVL和KCL列微分方程,整理。③由描述系统的输入输出动态方程式(微分方程)或传递函数,建立系统的状态空间表达式,即实现问题。实现是非唯一的。方法:微分方程一系统函数f模拟结构图f状态空间表达式。熟练使用梅森公式。注意:a如果系统函数分子累次等于分母幕次,首先化成真分式形式,然后再继续其他工作。b模拟结构图的等效。如前馈点等效移到综合反馈点之前。p28c对多输入多输出微分方程的实现,也可以先画出模拟结构图。.状态矢量的线性变换。也说明了状态空间表达的非唯一性。不改变系统的特征值。特征多项式的系数也是系统的不变量。特征矢量Pi的求解:也就是求(4/-A)尤=0的非零解。状态空间表达式变换为约旦标准型(A为任意矩阵):主要是要先求出变换矩阵。a互异根时,各特征矢量按列排。b有重根时,设3阶系统,4=4,4为单根,对特征矢量乃,P3求法与前面相同,”2称作4的广义特征矢量,应满足(41-A)p2=-P]。系统的并联实现:特征根互异;有重根。方法:系统函数f部分分式展开一模拟结构图f状态空间表达式。.由状态空间表达式求传递函数阵W(s)W(s)=C(sI-A)1+B+D6xr的矩阵函数[1%]叱)表示第j个输入对第i个输出的传递关系。子系统的并联、串联、反馈连接时,对应的状态空间表达及传递函数阵W(s)。方法:画出系统结构图,理清关系,用分块矩阵表示。.离散系统的状态空间表达式及实现(模拟结构图)x(k+1)=Gx(k)+Huy(k)=Cx(k)+Du.时变系统:四个矩阵是时间t有关的。非线性系统:各微分方程组的右端含有状态变量的非线性项。利用泰勒级数可以线性化。第二章控制系统状态空间表达式的解--线性定常系统齐次状态方程(*=Ax)的解:x(t)=eA'x0二.矩阵指数函数——状态转移矩阵.。⑴=e"表示x(0)到x(f)的转移。5个基本性质。. 的计算:a定义;b变换为约旦标准型A(^J)=T'ATyeA,=TeA'T-'^TeJ,T-'c用拉氏反变换e"=L 记忆常用的拉氏变换对\ 1 -at 1/1 -at 1 . (0 So(t)<->l;l(r)c-;r<->—<-> ;r———-\te- -;coscotc- -ss s+。s” (s+a) s+co" s4-cod应用凯莱・哈密顿定理三.线性定常系统非齐次方程(*=Ax+3“)的解:xQ)=0(r)x(O)+[/(r-7)3〃(7)d7。可由拉氏变换法证明(当然给出拉氏变换法的求解思路)。求解步骤:先求然后将B和u⑴代入公式即可。特殊激励下的解。四.线性时变系统的解.状态转移矩阵用。来表示。.。(此)的计算:当A")/。人(/必7=J。A(/W9(')时,。(,,,0)=exp[j:4加7];通常不等。不满足乘法可交换条件时,一般采用级数近似法:。(//o)=/+J。A(r)t/r+J。A(r0)j:4%)d/d7+••.解为:xQ)=0Q/o)x«o)+f。(/,7)8(丁)〃(7)]不五.离散时间系统状态方程的解(递推法和Z变换法)1.递推法0(。)=G*为状态转移矩阵;满足队k+1)=G0(Q;0(O)=I1 4一1解为,X(k)=。(幻x(0)+>(k-j-DHu(j)d曲(k)=。伏)x(0)+£。(忆-j-1)”〃(力dr;=o j=o
直接计算。(Z)=G"有一定困难,可采用这样的步骤:先将原状态方程化为约旦标准型,求变换矩阵T,x(k)=7x(k),再求出工(幻,再得到当然$(k)=A*,楸)=6卜=T$(kyr'.2.Z变换法 公式不用记忆,现推最好。x(Z)=Z-'[(z/-G)-1zr(O)]+Z-'[(z/-G)~'Hu(z)];可见。(Q=Gk=Z-'[(zZ-G)-1z];| 7计算x(Z)的用到的内容:部分分式展开(先除z后乘z);ZT对ak -= ;k>01-uzz-a六.连续时间状态空间表达式的离散化.定常系统的离散化G(T)=eAT;H(T)=「eA,dt-BJo即G(T)xTA+I;H(T)=TBx=Ax+Bux(k4-1)=G(T)=eAT;H(T)=「eA,dt-BJo即G(T)xTA+I;H(T)=TBb.近似离散化x-k+1)T)=(7X+I)x(kT)+TBu(kT)
y(k)=Cx(k)+b.近似离散化.时变系统的离散化略第三章线性控制系统的能控性和能观性一.能控性及能观性定义(线性连续定常、时变系统,离散时间系统)二.线性定常系统的能控性判别(具有一般系统矩阵的多输入系统)判别方法(一):通过线性变换x=Ax+Bu-^z=T-'ATz+T-'Bu.若A的特征值互异,线性变换(x=Tz)为对角线标准型,八=7-"丁,能控性充要条件:丁一力没有全为0的行。变换矩阵T的求法。.若A的特征值有相同的,线性变换(x=Tz)为约当标准型,J=T'AT,能控性充要条件:①对应于相同特征值的部分,每个约当块对应的7^8中最后一行元素没有全为0的。②7^3中对应于互异特征根部分,各行元素没有全为0的。变换矩阵T的求法。这种方法能确定具体哪个状态不能控。但线性变换比较复杂,关键是求T、T'.T'B.判别方法(二):直接从A,B判别x^Ax+Bu能控的充要条件是能控性判别矩阵M=(8,A8,A2b,…从一砂的秩为储在单输入系统中,M是一个〃x〃的方阵;而多输入系统,M是一个〃x〃「的矩阵,可通过ra〃0W=ra〃左(MM')三.线性定常系统的能观性判别x=Axz=T~'ATz判别方法(一):通过线性变换 cTy=Cxy=TCz.若A的特征值互异,线性变换(》=改)为对角线标准型,A=TTat,能观性充要条件:TC中没有全为0的列。变换矩阵T的求法。.若A的特征值有相同的,线性变换(x=Tz)为约当标准型,J=T'AT,能控性充要条件:①对应于相同特征值的部分,每个约当块对应的TC中第一列元素没有全为0的。②对应于互异特征根部分,对应的TC中各列元素没有全为0的。变换矩阵T的求法。这种方法能确定具体哪个状态不能观。但线性变换比较复杂,关键是求T、T-'.TC。判别方法(二):直接从A,C判别’C、能观性的充要条件是能观性判别矩阵N=.的秩为n。3,在单输入系统中,N是一个〃x〃的方阵:而多输入系统,N是一个〃机x〃的矩阵,可通过ra〃攵Af=ra〃攵(MM')四.离散时间系统的能控性与能观性x(k+l)=Gx(k)+Hu(k) , „.," , 能控性充要条件加=(”,6”,62”「..6向”)的秩为11。y(k)=Cx(k)+Du(k)'c、_ CG能控性充要条件N= : 的秩为n。、CG"L五.时变系统的能控性与能观性(与定常系统不同).比=A(r)x+8(r)”在「0,。]上状态能控的充要条件是格拉姆矩阵叱.(品。)非奇异。叱.4"f)=f'Q),(to,t)dt 。仇,。与。。也)一样么?这种方法要求先计算出状态转移矩阵,如果无法写成闭解,则失去工程意义。.使用A⑴B⑺信息Q,⑺=(B,⑺,一⑴,…,纥⑺),其中与⑺=B(t),与⑺=-A(f)%⑴+Ry⑺如果存在某个时刻>0,使得m〃2Q.(0)=〃,则系统在[0,。]上是状态完全能控的。.能观性判别与能控性类似,也可以使用格拉姆矩阵叱,(、//),但工作量太大。可使用c(r)信息:匕⑺、RQ)=2:"),其中G(f)=C«),2⑺=A(f)C,.T⑺+Ct⑺如果存在某个时刻。>0,使得加成/?(0)=〃,则系统在[0,。]上是状态完全能观测的。六.能控性与能观性的对偶原理1.若A2=A1,B]=C:,C2=B[,则Z](A,8],G)与£2(42,82,。2)对偶。
2.与1对偶,则%能控性等价于Z能观性,Z能观性等价于[能控性。时变系统的对偶原理????七.能控标准型和能观标准型对于状态反馈,化为能控标准型比较方便;对于观测器的设计及系统辨识,能观标准型比较方便。.能控标准I型(如果已知系统的状态空间表达式)①判别系统的能控性。②计算特征多项式即可写出彳。③求变换矩阵c=cTc],也可以验证是否有无=7;JA.能控标准n型①判别系统的能控性。②计算特征多项式|〃一A|=犬I+a“T/rl+…%/l+ao,即可写出彳。II③求变换矩阵,2=g,Ab,…,。④求nJ,计算5=7;2一%=.,c=cTc2,也可以验证是否有=小配。.能观标准【型①判别系统的能观性。②计算特征多项式…。伉+%,即可写出彳。③求变换矩阵CcA。CcA。④求与,计算5 c=cT<a=[10••o],也可以验证是否有Z=w.能观标准n型①判别系统的能观性。②计算特征多项式|2/一4|=/1"+。“_|/11+~。伉+。0,即可写出彳。③求变换矩阵7;2=h7;2=h,A7;「、A"Z],④求7;2,计算5=7;2-%,C=CTO2=[00…1],也可以验证是否有A=T„^ATo2..如果已知传递函数阵,可直接写出能控标准I型和能观标准n型的状态空间表达。w(s)= 3+必产+…+。"风丁y+a”〔si+a„ysn~2+---+a.s+ann-i i u
-0 1 0 …0一°〕0 0 1 •••0o能控标准I型:A=・・・bc=[/?()A•••A-il0 0 0 •••10一4 -q—6Z-) •••~an-\i00…0—4A.'10•••0一4A能观标准n型:a=01…0—h=♦c=[00…1]••・・・・•图200…1 —CL._ n—i_A-._八.线性系统的结构分解.按能控性分解(状态不完全能控,即依〃内0=勺<〃),通过非奇异变换x=R/完成。&=(&R2…%…R”),前多个列矢量是M中/个线性无关的列,其他列矢量保证R,非奇异的条件下是任意的。.按能观性分解(状态不完全能观,即m成N="<〃),通过非奇异变换工=凡戊完成。,\R;R,“,前〃।个行矢量是N中〃।个线性无关的行,其他行矢量保证此「非奇异的条件下是任意的。R,“.按能控性和能观性分解(系统是不完全能控和不完全能观的),采用逐步分解法,虽然烦琐,但直观。步骤:①首先按能控性分解(尤,能控状态,与不能控状态)。②对不能控子系统按能观性分解(x&不能控能观状态,X)不能控不能观状态)。③将能控子系统按能观性分解(X,”能控能观状态,能控不能观状态)。④综合各步变换结果,写出最后的表达式。另一种方法:化为约当标准型,判断各状态的能控性能观测性,最后按4种类型分类排列。九.传递函数阵的实现问题.实现的定义:由W(s)写出状态空间表达式,甚至画出模拟结构图,称为传递函数阵的实现问题。条件:①传递函数阵中每个元的分子分母多项式都是实常数:②元是S的真有理分式。注意:如果不是有理分式,首先求出直接传递矩阵。=iimW(s)。S—>00.能控标准型和能观标准型实现单入单出系统,W(s)是有理分式,可直接根据分子分母多项式系数写出能控标准1型和能观标准2型实现。多输入多输出系统,W(s)是矩阵,将W(s)整理成和单入单出系统传递函数相类似的形式,即W(s)= :+%丁2+…+南+自 ;此时的饱0、...\是mXr维常数阵。其能控标准S+an-\S+an-2S+…+。1S+。0型和能观标准型实现与单入单出系统类似,只是各矩阵中的0变为全零矩阵,1变为单位矩阵1,常数变为常数乘单位矩阵,即一即一>-《)/。注意:能控标准型实现的维数是〃xr;能观标准型实现的维数是〃X%。.最小实现(维数最小的实现)x=Ax+Bu八 为W(s)最小实现的充要条件是2(AB,C)是完全能控能观的。y=Cx步骤:对给定的W(s),初选一种实现(能控标准型或能观标准型),假设选能控标准型,判断是否完全能观测,若完全能观测则就是最小实现;否则进行能观性分解,进一步找出能控能观部分,即为最小实现。注意:传递函数阵W(s)的实现不是唯一的,最小实现也不是唯一的。十.传递函数W(s)中零极点对消与能控性和能观性之间的关系对单输入系统、单输出系统或者单输入单输出系统,系统能控能观的充要条件是传递函数没有零极点对消。而对多输入多输出系统,传递函数阵没有零极点对消只是最小实现的充分条件,也就是说,即使存在零极点对消,系统仍有可能是能控能观的(pl47例3-19)。对单输入单输出系统,若传递函数出现了零极点对消,还不能判断到底是不能控还是不能观,还是既不能控又不能观。第四章稳定性与李雅普诺夫方法--稳定性的定义李雅普诺夫给出了对任何系统都普遍适用的稳定性定义。.平衡状态£=/(xj)为齐次状态方程。满足对所有t,都有/(演/)三0成立的状态矢量互称为系统的平衡状态。稳定性问题都是相对于某个平衡状态而言的。通常只讨论坐标原点处的稳定性。.稳定性的几个定义①李雅普诺夫意义下稳定,(相当于自控里的临界稳定);②渐近稳定,(相当于自控里的稳定);③大范围渐近稳定,大范围渐近稳定的必要条件是整个状态空间只有一个平衡状态;④不稳定。二.李雅普诺夫第一法(间接法).线性定常系统的稳定判据状态稳定性:平衡状态%=0渐近稳定的充要条件是A的所有特征值具有负实部。输出稳定性:充要条件是传递函数的极点位于s的左半平面。.非线性系统的稳定性线性化处理。= 4=或,若A的所有特征值具有负实部,则原非线性系统在平衡状态覆渐dxU-4x=Xe近稳定。若A的所有特征值至少有一个具有正实部,则原非线性系统在平衡状态兀不稳定。若若A的所有特征值至少有实部为零,则稳定性不能有特征值的符号来确定。三.李雅普诺夫第二法(直接法)借助于一个李雅普诺夫函数来直接对平衡状态的稳定性做出判断。.预备知识V(x)是由n维矢量x定义的标量函数,且在x=0处,恒有V(x)=0,对任何非零矢量x,如果V(x)>0,则称之为正定;如果V(x)<0,则称之为负定:如果V(x)20则称之为半正定或非负定;如果V(x)W0则称之为半负定或非正定;如果V(x)>0或V(x)<0,则称之为不定。V(尤)=/分为二次型标量函数,P为实对称阵。要判别V(x)的符号只要判别P的符号即可。P的定号判据(希尔维特斯判据):首先求出产的各阶顺序主子式△,,若所有的Aj〉。,则P(V(x))正定:若,=偶数的>0,i=奇数的<0则尸(V(x))负定;.李雅普诺夫函数对于一个给定系统,如果能找到一个正定的标量函数丫(x),而V(x)是负定的,则这个系统是渐近稳定的,这个标量函数V(x)叫做李雅普诺夫函数。李雅普诺夫第二法的关键问题就是寻找李雅普诺夫函数V(x)的问题。.稳定性判据①设£=/(x),平衡状态为乙=0,如果存在标量函数V(x)是正定的,即x=0时,有丫(x)=0,xhO时,有V(x)>0,且满足U(x)<0,则称原点平衡状态是渐近稳定的;如果当Mlf8时,V(X)f8,则系统是大范围渐近稳定的。②设£=/(x),平衡状态为x«=0,如果存在标量函数V(x)是正定的,即x=0时,有丫(x)=0,x#0时,有V(x)>0,且满足廿(x)WO,但除x三。外,即xhO,R(x)不恒等于0,则称原点平衡状态是渐近稳定的;如果当H->8时,V(X)T8,则系统是大范围渐近稳定的。③设£=f(x),平衡状态为x,=0,如果存在标量函数V(x)是正定的,即x=0时,有V(x)=0,x#0时,有V(x)>0,且满足廿(x)W0,但任意的xhO,廿(x)恒等于0,则称原点平衡状态是李雅普诺夫意义下稳定的。④设£=/(x),平衡状态为乙=0,如果存在标量函数V(x)是正定的,即x=0时,有丫(x)=0,xhO时,有V(x)>0,且满足廿(x)>0,,则称原点平衡状态是不稳定的。需要注意:①这些判据定理知识充分条件,也就是说,没有找到合适的李雅普诺夫函数来证明原点的稳定性,不能说明原点一定是不稳定的。②如果V(x)是可找到的,那么通常是非唯一的,但不影响结论。③V(x)最简单的形式是二次型标量函数,但不一定都是简单的二次型。④构造V(x)需要较多技巧。四.李雅普诺夫方法在线性系统中的应用.线性定常连续系统渐近稳定判据定理:£=Ax,若A是非奇异的,原点%=0是唯一的平衡点。原点大范围渐近稳定的充要条件是对任意对称实正定矩阵Q,李雅普诺夫方程A7p+P4=-Q,存在唯一的对称正定解P。该定理等价于A的特征值具有负实部。但高阶系统求解特征值复杂。步骤:选定正定矩阵。,通常为Q=/,代入李雅普诺夫方程,确定出P,判断是否正定,进而做出系统渐近稳定的结论。.线性时变连续系统渐近稳定判据定理:x=A(t)x,在平衡点%=0大范围渐近稳定的充要条件是对任意对称实正定矩阵QQ),李雅普诺夫方程户⑺=-A(t)TP⑺-尸⑴A⑺-Q"),存在唯一的对称正定解尸⑴。.线性定常离散系统渐近稳定判据定理:%优+1)=6双外在平衡点七=0渐近稳定的充要条件是,对任意对称实正定矩阵。,离散李雅普诺夫方程GTpG-P=—Q,存在唯一的对称正定解P。该定理等价于G的特征值均在单位圆内。步骤:选定正定矩阵Q,通常为Q=/,代入离散李雅普诺夫方程,确定出P,判断是否正定,进而做出系统渐近稳定的结论。五.非线性系统的李雅普诺夫稳定性分析.雅可比矩阵法步骤:X=/(X),写出/(X),计算雅可比矩阵J(x)=j,对给定正定矩阵P(通常尸=/),dxQ(x)=~[JMTP+PJ(x)]为正定的。并且V(x)=fT(x)Pf(x)为系统的一个李雅普诺夫函数。.变量梯度法第五章线性定常系统的综合综合:常规综合,使系统性能满足某种笼统指标要求;最优综合,使系统性能指标在某种意义下达到最优。一.线性反馈控制系统的基本结构及其特性.状态反馈将系统的每一个状态变量乘以相应的反馈系数,然后反馈到输入端与参考输入相加,作为受控系统的控制输入。K称为状态反馈增益阵,rx〃。设原受控系统Zo=(A8,C),D=0»x=(A+BK)x+Bv状态反馈闭环系统的状态空间表达式 简称Zk=(A+8K,8,C)y=Cx与原受控系统=(A,B,C)比较,状态反馈增益阵K的引入,并不增加系统的维数,但可以通过K的选择改变闭环系统的特征值,从而使获得所要求的性能。.输出反馈 由输出端y引入输出反馈增益阵H(rxa),然后反馈到输入端与参考输入相加,作为受控系x=(A+BHC)x+Bv 「统的控制输入。状态空间表达式为 简称=(A+B”C,5,C)y=Cx通过H的选择也可以改变闭环系统的特征值,从而改变性能,但可供选择的自由度远比K小(通常加<〃)。.从输出到状态变量导数上的反馈 从输出y引入反馈增益阵G(〃x/n)到状态变量的导数比,所得状态空间表达式为 简称Zh=(a+gc,b,c)y=Cx通过G的选择也可以改变闭环系统的特征值,从而改变性能。以上三种反馈的共同点是,不增加新的状态变量,系统开环与闭环同维,其次,反馈增益阵都是常数矩阵,反馈为线性反馈。.闭环系统的能控性与能观性a状态反馈不改变受控系统=(A,8,C)的能控性,但不保证系统的能观性不变。b输出反馈不改变受控系统=(A,B,C)的能控性和能观性。二.极点配置问题 就是通过选择反馈增益矩阵,将闭环系统的极点恰好配置在根平面所期望的位置,以获得所希望的动态性能。 只讨论单输入单输出系统.采用状态反馈 对系统X。=(A4c)任意配置极点的充要条件是X。完全能控。给定Zo=(A*,c),给定期望的极点,设计状态反馈控制器的方法:⑴能控规范型法,适合于〃23。①首先判断是否完全能控,是,则存在状态观测器。②通过线性变换尤=7;3化为能控标准1型,得到5=(彳石Z)。③加入状态反馈增益矩阵匠=[%,无,…,匕_/,得到闭环系统&=(X+5k,瓦O状态空间表达式,求出对应的闭环特征多项式/(团=|〃一(才+5/。|。④由给定的期望极点,求出期望的闭环特征多项式/*(/i)=nu-4*)。⑤将/U)与/*(㈤比较,即可得到匠= …⑥把对应与F的冗,通过k=N&t=%,占,…次⑦进一步画出模拟结构图。⑵当阶次较低时,n<3,可直接由反映物理系统的A,b矩阵求状态反馈增益矩阵K=[A0,匕,…,£_J,不通过非奇异变换,使设计工作简单。①首先判断是否完全能控,是,则存在状态观测器。②加入状态反馈增益矩阵K=[%,匕,…,匕,_/,得到闭环系统=(A+bK,b,c)状态空间表达式,求出对应的闭环特征多项式/(2)^Ai-(A+bK)\,③由给定的期望极点,求出期望的闭环特征多项式/*(丸)=n(A-V).④将/(2)与/*(之)比较,即可得到K=[%«,…次⑤进一步画出模拟结构图。注意,如果给定的是传递函数,则先画出其要求的模拟结构图,写出状态空间描述,然后做其他工作。.采用输出反馈不能任意极点配置,正是输出线性反馈的基本弱点。.采用从输出到土的反馈对系统=(A,"C)任意配置极点的充要条件是完全能观。设计Eo从输出到土的反馈阵G的问题就是其对偶系统£o设计状态反馈阵K的问题。方法:(1)能观标准型法,适合于"23。①首先判断是否完全能观,是,则存在输出反馈G。②通过线性变换x=To2x化为能观标准2型,得到£=(A,h,c)0③加入输出反馈增益矩阵G=[藐商”…商“_|了,得到闭环系统Ng=(•+]}石1)状态空间表达式,求出对应的闭环特征多项式/(㈤=|〃一(印+&0|。④由给定的期望极点,求出期望的闭环特征多项式/*(/1)=口(/1一4*)。⑤将/(㈤与/*(㈤比较,即可得到
G=[go,否,…,匹T]、⑥把对应与5的G,通过G=Ts1=醇0*°-,8“_/。©进一步画出模拟结构图。⑵当阶次较低时,n<3,可直接由反映物理系统的A,c矩阵求状态反馈增益矩阵G=[g0,g],…,g“_J,不通过非奇异变换,使设计工作简单。①首先判断是否完全能观,是,则存在输出反馈G。②加入从输出到文的反馈增益矩阵6=30,81,…,g“_/,得到闭环系统Eg=(A+Gc,瓦c)状态空间表达式,求出对应的闭环特征多项式f(/l)=|R-(A+Gc)|。③由给定的期望极点,求出期望的闭环特征多项式/*(/1)=0(/1—4*)。④将/(㈤与/*(/1)比较,即可得到6=岛.,…,g“_J。⑤进一步画出模拟结构图。三.系统镇定问题所谓系统镇定,是对受控系统=(A,8,C)通过反馈使其极点均具有负实部,保证系统为渐近稳定。镇定问题是极点配置问题的一种特殊情况,它只要求把闭环极点配置在根平面的左侧,而并不要求将闭环极点严格地配置在期望极点上。状态反馈能镇定的充要条件是其不能控子系统为渐近稳定。输出反馈能镇定的充要条件是结构分解中能控能观子系统是输出反馈能镇定的,其余子系统是渐近稳定的。输出到上的反馈实现镇定的充要条件是不能观子系统为渐近稳定。四.系统解藕问题.目的是寻求适当的控制规律,使输入输出相互关联的多变量系统实现每一个输出仅受相应的一个输入控制,每一个输入也仅能控制相应的一个输出,这样的问题称为解藕问题。.定义:若系统Z=(A"c)m维输入m维输出,其传递函数矩阵是一个对角线有理多项式矩阵,则称该系统是解藕的。匕⑸W(s)匕⑸W(s)=C(sl-A)-'B=0W")」.方法:①前馈补偿器解耦:待解耦系统Z°=(A,6,c)的传递函数阵W0(s),在其前面串接一个前馈补偿器传递函数为叼(S),使整个系统的传递函数阵为W(s)=%(s)Wo(s),满足对角线有理多项式特点。其中也(s)=%t(s)W(s)。②状态反馈解藕。如何设计K和F,使系统从v到y是解藕的。设计步骤。五.状态观测器作用:闭环极点的任意配置、系统解藕以及最优控制系统都离不开状态反馈。但状态变量并不是都能直接检测,有些根本无法检测,这就提出状态观测或状态重构问题。龙伯格提出的状态观测器理论,解决的状态重构问题,使状态反馈成为一种可实现的控制律。.定义:动态系统名以的输入u和输出y作为输入量,产生一组输出量£逼近于x,即lim|x—£|=O,则称£为z。的一个状态观测器。构造原则:必须是完全能观或不能观子系统是渐近稳定的:£的输出£应以足够快的速度渐近于X;七在结构上尽可能简单(具有尽可能低的维数),以便于物理实现。.等价性指标«x=Ax+Bu _x=Ax+Bu动态系统z 原系统工y=cx y=cxx-x=A(x-x)得到x-£=e*'(x()-£())只要系统是稳定的,即A的特征值具有负实部,就可做到£与力是稳态等价的。3.重构状态方程原因:①系统的状态是不能直接量测的,因此很难判断是否有£逼近于x;②不一定能保证A的特征值均具有负实部。克服这个困难,用对输出量的差值y-亍的测量代替对状态误差的测量,当lim|x-£|=O,有lim|y-9|=lim-戌|=lim|c(x-£)|=O。同时,引入反馈阵G,使系统的特征值具有负实部。/—>00 /—>00 /—>00状态重构方框图为p2135.16(a)要求熟练记忆,这种状态观测器称为渐近观测器。状态观测器方程为f=至+8“+G(y-9)=(A-GC).i+Gy+Bu记为£=(4_GC氏G)y=Cx这里的G称为输出误差反馈矩阵。可以证明,如果A-GC的特征值具有负实部,那么状态误差将逐渐衰减到0,即估计状态£逼近于实际的状态X。逼近的速度取决于G的选择,即A—GC的特征值的配置。4.观测器的存在性对于完全能观测的线性定常系统,其观测器总是存在的。观测器存在的充要条件是Eo不能观子系统是渐近稳定的。5.观测器的极点配置定理:线性定常系统%=(A8,C),其观测器£=(A—GC,3,G)可以任意配置极点,即具有任意逼近速度的充要条件是Eo= 完全能观测。极点配置方法:(1)能观标准型法,适合于〃23。①首先判断是否完全能观,是,存在观测器可以任意极点配置。②通过线性变换》=笈化为能观标准2型,得到5=(不,5足)。③加入输出误差反馈阵G=[蔡,耳,…,匹得到闭环系统状态空间表达式2=(A-G3)2+8"+Gy),求出对应的闭环特征多项式/(团=|〃一(•一13)④由给定的期望极点,求出期望的闭环特征多项式=口(丸一4*)。⑤将/U)与/*(九)比较,即可得到3=[疏,⑥,…商⑥把对应与5的3,通过G=7后=[go,gc、g“T]。©得观测器方程,x=(A-Gc)x+Bu+Gy^Ji=Ai+Bu+GCy-y),进一步画出模拟结构图。⑵当阶次较低时,n<3,可由特征值不变原理求状态反馈增益矩阵G=[g0,g「…,g”_J,不通过非奇异变换,使设计工作简单。①首先判断是否完全能观,是,则存在观测器可以任意极点配置。②引入输出误差反馈矩阵6=旧0,8。…,8"--,得到观测器系统£=(4—6。,5,6)状态空间表达式(=(4—6(')£+8〃+6、。③求出对应的闭环特征多项式f(/l)=|〃-(A-Gc)|。④由给定的期望极点,求出期望的闭环特征多项式/*(4)=n(”4*)。⑤将/(㈤与/*⑷比较,即可得到6=/0闾,…,⑤得观测器方程,进一步画出模拟结构图。
5.降维观测器观测器维数与受控系统相同,称为全维观测器。如果有些状态变量能由输出y直接获得,那么仅对其余的状态变量用降维观测器进行重构即可。步骤:①通过线性变换把状态按能检测性分解。(n-m)维状态变量看需要重构,m维状态变量另由直接获得。②对可构造(n-m)维观测器。 详细步骤通过实例熟悉。六.利用状态观测器实现状态反馈的系统(带观测器的状态反馈闭环系统).系统的结构与状态空间表达结构框图要非常熟悉p221图5.21前提:受控系统完全能控能观,状态反馈闭环系统和观测器都可以任意极点配置。受控系统=(A,5,C)受控系统=(A,5,C)x=Ax+Buy=exm+TTi'mtwbzx„„nx=Ax+Bu+G(y-y)=(A-GC)x+Gy+Bu,n.状态观测器Zg=(A_GC,8,G) z . *2式y=Cx反馈控制率〃=U+依 *3式x=Ax+BKjc+Bv整理得整个闭环系统的状态空间表达式i=GCx+(A-GC)x+Gy+Bv也可写成矩阵形式y=Cx显然,这是一个2n维的闭环控制系统。.闭环系统的基本性质(1)分离性复合系统(由观测器构成的状态反馈闭环系统)其特征多项式等于矩阵A+BK和A-GC特征多项式的乘积。即闭环系统的极点等于直接状态反馈(A+BK)的极点和状态观测器(A-GC)的极点总和,且相互独立。所以输出误差反馈阵G和状态反馈阵K可以分别进行设计。(2)传递函数矩阵的不变性可以推出复合系统的传递函数为W(s)=C[sI-(A+BK)ylB,等于直接状态反馈闭环系统的传递函数。或者说它与采用观测器反馈无关。(3)观测器反馈与直接状态反馈的等效性稳态时,两者等价。选择K,可以改变闭环系统的极点到期望极点,从而改善系统性能。选择G,可以改变观测器的极点,从而加速使状态误差衰减到0。一般取观测器的极点比闭环系统的期望极点((A+6K)的极点)略负,既保证状态误差有较快的衰减速度,又不致引人更多的噪声干扰。.设计步骤(只给出低阶系统的设计步骤):①判断原受控系统的能控性能观性,是完全能控能观,则状态反馈阵K和观测器输出误差反馈阵G存在,且闭环系统和观测器极点可以任意配置。②设计状态反馈阵K:求A+BK的特征多项式力(丸),由期望的闭环极点得期望的特征多项式//(4),比较系数,从而得到K。③设计观测器输出误差反馈阵G:求A-GC的特征多项式上。),由观测器期望的配置极点得期望的特征多项式九*(九),比较系数,从而得到G。④给出观测器方程即*2式。⑤结合*1式和*3式,画出相应的模拟结构图。第六章最优控制三种设计最优控制系统的方法:古典变分法、极小值原理、动态规划一.概述在最优控制系统中,由于受控对象是一个动态系统,所有变量都是时间的函数,所以这是动态最优化问题。这时目标函数不再是普通的函数,而是时间函数的函数,称为泛函。在[品。]上目标泛函为,=『〃⑺/W,基本约束条件是受控对象的状态方程%=/[x(r),«(OJ],J是标量泛函数,L标量函数(是矢量u⑴,x⑴的函数),x⑴是n维状态矢量,u(t)是r维控制矢量。二.研究最优控制的前提条件(1)给出受控系统的动态描述,即状态方程连续£=/Ix(f),〃(r),r]离散x(Z+l)=/[x(&),〃(Z),Z](2)明确控制作用域(3)明确初始条件通常书给定,若x«o)给定,称为固定始端。若Ml。)任意,则称谓自由始端。(4)明确终端条件固定终端自由终端可变终端(5)给出目标泛函即性能指标, N-I连续J=0[x(")]+1 离散 J=0[x(N)]+Z〃x(r),“⑺用,0 k=k0等式右边第一项反映对终端性能的要求,称为终端指标函数。第二项L为状态控制过程中对动态品质及能量或燃料消耗的要求等,称为动态指标函数。若不考虑终端指标函数,仅有第二项则称为拉格郎日型(或积分型)。若仅有第一项,则称为终端型(梅耶型)。最优控制问题就是在约束条件下寻求最优控制u(t),受控系统在上,从初始状态Nr。)转移到终端状态x(。)时,性能指标J取极值。满足条件的U⑴称为最优控制“*«),这时状态方程的解称为最优轨线X"(r),此时的性能指标J称为最优指标J*。三.静态最优化问题的解(1)多元函数的极值•/=/(“)这里〃=(%,%,…,勺)'df d2f取极值的必要条件是乜■=(),取极小值还需满足汩•>() 海赛矩阵正定。du 8u-(2)具有等式约束的极值a嵌入法 先从目标函数解出一个变量,代入目标函数,即成为没有目标约束的函数。b拉格朗日乘子法 将约束条件乘以;I,与目标函数相加,构成一个新的可调整的没有约束的多元函数。目标函数)=/(x,〃)约束条件g(x,〃)=O,新函数”=/(x,w)+Xg(x,a)九是与g同维的列矢量。目标函数存在极值的必要条件是迫=0 —=0 —=0dxdu52四.泛函及其极值——变分法动态最优控制中的目标函数是一个泛函数,因此动态最优化问题可以归结为求泛函极值问题。1.变分法概念在控制系统中,自变量是3宗量函数是状态矢量x(f),因此J=『力,而£= ,所以,j可以写成,=r 是积分型泛函。j的值取决于函数以力,所以j是〃⑺的泛函。求最优控制〃*«),就是寻求使性能泛函J取极值的N”)oa泛函的变分: 泛函4y(x)]的变分定义为,d/=—J[y(x)+ady(x)]多元函数的变分:8J=—J[y,+a8yx,y2+a私,…,y“+a6yn]ba a=o多元函数取极值的必要条件是8=02.泛函极值的必要条件一一欧拉方程求泛函j=r 的极小值,就是确定X”)使j达到极小值。定理:设曲线X(r)的始点为x(fo)=x0,终点为X(O)=X/,则使性能泛函1/=[*〃乂£"]^取极值的必要条件是:X(r)是二阶微分方程更一4(更)=0的解。 池一4(更)=0称为欧拉方程。dxdtdx dxdtdx_.d/BL、 ddL[x,x,t] ddLdxddLdxddLdtdtdx dtdx dxdxdtdxdxdtdtdxdt实例熟悉步骤。欧拉方程是二阶微分方程,求解时有两个常数待定。对固定端点问题,给定x«o)=%,x(0)=与边界条件,可以确定常数。对于自由端点问题,应有横截条件来补足。—=0—=0/ dx,dx\,zo ,r/3.多元泛函的极值条件一一欧拉方程组1/(X],%,…,X”)=['L[x],x2,---,xn,xi,x2,---,xn,t}dt取极值的必要条件是:为⑴是二阶微分方程组以一旦(以)=0的解。 实例熟悉步骤。p252dXjdtdXj4.可变终端问题和综合型性能泛函的情况 略五.用变分法求解连续系统最优控制问题一一有约束条件的泛函极值前面讨论的是没有约束的泛函极值问题。有约束条件的泛函极值,解决思路:应用拉格朗日乘子法,构造增广泛函,转化为没有约束条件的极值问题。1.问题描述受控系统的状态方程为£=/[宜。,〃(力用,给定初始状态》依)=%,求最优控制〃*(f),在上使目标性能泛函为J=0(x(。)]+『〃x(r),“(f)j]力取得极值。最优控制是一类条件泛函极值问题,求解时首先要定义哈密顿函数"=L(x,u,t)+^f(x,u,t)2.不同边界状态下的最优控制必要条件(1)末态无约束,但末态时间。固定。此时最优控制的必要条件(由欧拉方程推出的)为dHn -3H SH. 3、Q,、-^—=0 4=--— ——=x=f(x
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025年志愿者工作方案
- 2025年卖场活动促销方案
- 汽车使用与维护 课件 项目二 行驶系统的使用与维护2-3 四轮定位综合检测维修
- 2025年电子式电动套筒调节阀项目可行性研究报告
- 2025年电吉他袋项目可行性研究报告
- 2025年玻纤纱窗项目可行性研究报告
- 2025年爪型螺帽项目可行性研究报告
- 内蒙古百校联盟2025届高三下学期生物试题(月考)独立作业1含解析
- 江苏理工学院《输油管道设计》2023-2024学年第二学期期末试卷
- 永城职业学院《食品安全卫生学》2023-2024学年第二学期期末试卷
- 小径分岔的花园
- 超星尔雅学习通《孙子兵法》与执政艺术(浙江大学)网课章节测试答案
- 勾头作业施工方案
- 《叩问师魂》观后感3篇
- 出版专业基础知识中级
- GB/T 9575-2013橡胶和塑料软管软管规格和最大最小内径及切割长度公差
- GB/T 9163-2001关节轴承向心关节轴承
- GB/T 4857.19-1992包装运输包装件流通试验信息记录
- GB/T 37573-2019露天煤矿边坡稳定性年度评价技术规范
- ICU院内感染控制
- GB/T 24528-2009炭素材料体积密度测定方法
评论
0/150
提交评论