【文档说明】计算机控制系统第七章课件.ppt，共(76)页，693.117 KB，由小橙橙上传

转载请保留链接：https://www.ichengzhen.cn/view-76804.html

以下为本文档部分文字说明：

基于传递函数模型的极点配置与最优化设计方法的比较：（1）极点配置设计方法：（a）考虑确定性的跟踪系统；（b）性能指标是闭环模型传递函数，给定较困难。（2）最优化设计方法：（a）主要讨论随机的调节系统；（b）性能指标为输出量或

输出量加控制量的二次函数，给定较容易。第7章基于传递函数模型的最优化设计方法第一节设计问题一、对象及干扰模型控制对象y(t)u(t)v1c(t)图1（a）图（a）中：u(t)----控制量，零阶保持器的输出。v

1c(t)---随机的过程干扰。y(t)-----输出量。Gp(s)y(t)u(t)v1c(t)图1（b）图（b）中：Gd1(s)++)()()(sGsusyp)()()(11sGsvsydcGp(s)y(t)u(t)v1c(t)图1（c）Gd1(s)

++Gd2(s)ec(t)图（c）中：v1c(t)----不是白噪声，为具有有理功率频谱密度函数的有色噪声。ec(t)----白噪声。Gd2(s)----假想的动态环节的传递函数。Gp(s)y(t)u(t)ec(t)图

1（d）图（d）中：Gd(s)++)()()(21sGsGsGddd控制对象离散化，对于Gp(s)，假设u(t)是零阶保持器的输出，采用零阶保持器法，将Gp(s)化为Gp(z)。ssGZzsGpp)()1()(1（1）Gd(

s)的离散化方法：CxBeAxxcc（2）即：将Gd(s)化为如下的状态方程：设ec(t)是均值为零、方差为的白噪声，进一步设：2cccBee（3）则（2）式变为：CxeAxxcc（4）2)]()([

cTTccBBteteE（5）可求得的协方差为：ce式（4）可离散化为：)()()()()1(kCxkkekFxkxcd（6）其中：dTkTeekeeFTcAdAT0)()(,（7）式中，T为采用周期

，是等效的离散白噪声序列，其协方差为：)(kedTBBdeBBeQcTcTATAT220（8）显然有：TBDDDQcT,（9）从而式（6）可以写为：)()()()()

1(kCxkkDekFxkxcd（10）其中ed(k)是均值为零、方差为1的白噪声序列。从而可以求得等效得离散传递函数为：DFzICzGd1)()(（11）其等效离散控制对象的结构图如下：Gp(z)y(k)u(k)ed(k)图2Gd(z)++设)()()

(,)()()(2111zAzCzGzAzBzGdp（12）由图2，有)()()()()()()(2111kezAzCkuzAzBkyd（13）经通分进一步化为：)()()()()()()(kezAzCkuzAzBkyd（14）其中A(z)为首一多

项式。0c)(zC（15）其中有：)()(,)()(00keckeczCzCd（16）2202202)()(ckEeckEed（17）设首项系数为，则（14）式变为：)()()()()()()(kezAzCku

zAzBky从而有于是，C(z)也变成首一多项式。则式（15）写为：)()()()()()(kezCkuzBkyzA（18）此式便是标准的控制对象及干扰模型。若令e(k)=0，则为确定性系统的传递函数模型。若令u(k)=0，则为ARMA(Autoregr

essionMovingAverage)随机过程模型。对于式（18）的标准模型，有如下条件：（1）A(z)和C(z)均为首一多项式；（2）degA>degB;（3）degC-degA=0；（4）C(z)的零点均在单位圆内。其中（3）（4）均由于C

(z)是噪声驱动模型的性质所决定。二、性能指标及容许控制0)(21kEyJ研究调节系统，即r(k)=0。B(z)/A(z)y(k)u(k)e(k)图3C(z)/A(z)++D(z)由于Ee(k)=0，则Ey(k)

=0但是希望J1越小越好。对于离散系统，性能指标可以表示为：NkNkyNEJ121)(1lim（19）对于连续系统，性能指标表示为：201lim()TcNJytdtT（20）具有上述性能指标的最优控制问题为最

小方差控制。为对控制量进行限制，可在指标中对控制量进行加权，即)()(222kukyEJ（21）其中为加权系数。上述指标为更一般形式的二次型性能指标，称为广义最小方差控制。设计问题：设计控制器D(z)，使二次型指标J1或J2最小。一、最优预报控制对象模型：)(

)()()()()()()()()()()()(1*1*1*1*kezAzCkuzzAzBkezAzCkuzAzBkyd（1）第二节最小方差控制其中)(deg)(deg)(deg)(degzBzAdzCzAn（2）y(k)u(k)e(

k)图1控制对象结构图++)()(1*1*zAzC)()(1*1*zAzBdz)(kyu)(dkyu)(k)()()(dkykkyu若使J1最小，最好使，则)()(kdkyu0)(2kEy（3）然而

，此控制不能实现。因为：)1()()()()()(kukkykukydkyuudz则)1()()()(dkudkdkydkyu因此，决不能完全抵消。)(k)(dky

u使J1最小的方法：根据k-d及以前的信息最好地估计出k时刻的干扰量，并使)/(ˆdkk)/(ˆ)(dkkdkyu从而)/(ˆ)()()()(dkkkdkykkyu（

4）于是221)/(ˆ)()(dkkkEkEyJ或221)/(ˆ)()(kdkdkEdkEyJ（5）（6）即：输出量的最小方差等于最优预报估计误差的方差。问题：求最优预报，使J1最小。)/(ˆkdk利用多项式

除法，得到：)()()()()(1*1*1*1*1*zAzGzzFzAzCd（7）其中)1(1111*1)(ddzfzfzF)()()()1()1()()()()()1()()()()(1*1*111*1*)1(1111*1*k

ezAzGkefdkefdkedkezAzGzzfzfdkezAzCdkdddd（8）（9）由图1得到：由式（6），得到：（10）21*1*2212121*1

*112)/(ˆ)()()()1()/(ˆ)()()()1()1()()/(ˆ)(kdkkezAzGEffkdkkezAzGkefdkefdkeEkdk

dkEdd其中22)(kEe要使（10）式取得最小值，必须有)()()()/(ˆ1*1*kezAzGkdk（11）22121222)1()]/(ˆ)([)()(dffddkdkEdkEykEy于是此即为最优

预报估计误差。（12）式（11）表示由e(k)，e(k-1)，…来获得最优预报，而我们希望由来获得，因此需要进行如下变换：)/(ˆkdk)/(ˆkdk),1(),(kk由图1，得到)()()()(1*1*kzCzAke（13）（的零点均在单位圆内）

)(1*zC代入（11）式，有)()()()/(ˆ1*1*kzCzGkdk（14）此即为最优预报公式。求最优预报计算步骤：（1）作多项式的带余除法运算，如式（7），即)()()()()(1*1*1*1*1*zAzGzzFzAzCd或写成：)()(

)()(1*1*1*1*zGzzFzAzCd（15）其中如式（8），即)(1*zF)1(1111*1)(ddzfzfzF（15）式两边同乘以，得到1dnz)()()()(1zGzFzAzCzd其中1211)(

dddfzfzzF（16）（17）从而可以求出和或F(z)和G(z)。)(1*zF)(1*zG（2）计算干扰量最优预报估计：)()()()()()()/(ˆ1*1*kzCzzGkzCzGkdk（18）最优预报估计误差如式（12），即22

1212)1()(dffkEy二、最小方差控制1、最小方差控制的实现：)/(ˆ)(dkkdkyu从而实现最小方差控制，控制器结构如图2所示。给定u(k)，使得图2最小方差控制系统结构图dz)(ky)()(1*1*zAzB)(1*zD)(ke)

(ku)(kyu)()(1*1*zAzC)(dkyu)(kv)/(ˆ)(kdkvkyu由得到)()()()()()()(1*1*11*1*kvzCzGkyzDzAzB于是有

*1*1*1*1*1()()()()()()()AzGzvkDzBzCzyk(19)由图2,得到)()()()()()(1*11kyzDzAzBzkvkyd于是得到)()()(1)()(1

*11zDzAzBzkykvd(20)将(20)代如(19),整理得到)()()()(1*1*1*1*zFzBzGzD或者)()()()(zFzBzGzD说明：（1）最小方差控制由两部分组成：（a）计算最优预报估计)/(ˆkdk（b）产生

最优控制，使得)/(ˆ)(kdkkyu()()()()()GzukykBzFz故与分离性原理相似。（2）D(z)将抵消对象的分子多项式B(z)，故最小方差控制只适用于B(z)零点均在单位圆内的情况。2、计算实例：控制对象模型：)()()(

)()()()(1*1*1*1*kezAzCkuzzAzBkyd211*11*211*6.05.11)(5.01)(7.07.11)(zzzCzzBzzz

A已知22)(kEe要求：计算最小方差控制器的传递函数及最小性能指标)(21kEyJ解：211121211*1*7.07.11)1.02.0(17.07.116.05.11)()(

zzzzzzzzzAzC于是有：11*1*1.02.0)(,1)(zzGzF从而5.0)5.0(2.0)(zzzD111*1*1*1*5.011.02.0)()()()(zzzFzBzGzD（1）延时

拍数d=1或最小方差为：221)(kEyJ（2）设延时拍数d=22112121211*1*7.07.11)14.024.0(2.017.07.116.05.11)()(zzzzzzzzzzAzC于是有：11*11*14.024.0)(,2.01)(

zzGzzF从而)2.0)(5.0(14.024.0)(2zzzzzD)2.01)(5.01(14.024.0)()()()(1111*1*1*1*zzzzFzBzGzD或最小方差为：222221104.1)2.01()1(

fJ可见，当延时增大时，最小方差也增大。3、与极点配置设计法的比较控制对象模型：)()()()()()(kezCkuzBkyzA（21）)()()()()()()(kyzFzBzGkyzDku由于故)()()()()()()()()(kezCkyzFzBzGzB

kyzA（22）（23）于是)()()()]()()()[()()()()(kezHkezGzFzAzBzCzBzFky（24）其中)]()()()[()()()()(zGzFzAzBzCzBzFzH（25）由式（16），即)()()()(1

zGzFzAzCzd得到)()()(1zCzBzzd)]()()()[()(zGzFzAzBz上式为e(k)到y(k)的闭环传递函数，其特征多项式为：（26）（27）（28）可见，系统的极点由三部分组成：（1）d-1个原极点；（2）B(z)的零点（n-d个

）；（3）C(z)的零点（n个）。极点配置设计法中，)()()()(kyzRzSku（29）则R(z)和S(z)满足如下的Diophantine方程：mAABBSAR0（30）选定,,,10dmzACABB则得到)()(),()()(zGzSzBzFzR（31）代入

（29）式，得到)()()()(zFzBzGzD（32）此即为最小方差控制。由此可见，最小方差控制是按极点配置设计方法的一个特例。三、对象具有单位圆外零点时的最小方差控制1、计算方法[定理1]给定控制对象的模型为：)()(

)()()()(kezCkuzBkyzA（33）将B(z)分解为)()()(zBzBzB（34）其中B+(z)包含所有单位圆内的零点（首一多项式），包含所有单位圆外和圆上的零点。假定C(z)的所有零点均在单位圆内，A(z)和互质，则最小方差控制为：

)(zB)(zB)()()()()(kyzFzBzGku（35）其中F(z)和G(z)满足如下的Diophantine方程：)()()()()()(*1zGzBzFzAzBzCzd（36）在上式中，求degG(z)<degA(z)的最小阶解

，F(z)和G(z)的阶次分别为：1)(deg)(degzBdzF1)(deg)(degzAzG（37）（38）)(*zB)(zB是的互反多项式。[定理证明略]定义：<1>互反多项式：设nnnnpzpzpzpzP1110)(称

01111*)()(pzpzpzpzPzzPnnnnn为P(z)的互反多项式。<2>I型最小方差控制：最小方差控制器抵消B(z)的全部零点。<3>II型最小方差控制：最小方差控制器只抵消B+

(z)而不抵消。)(zB说明：<1>I型最小方差控制使得输出方差达到极小值，但控制量可能趋于无穷大（抵消）。)(zB<2>II型最小方差控制使得输出方差达到有限制的极小值，而不是最小值，但它使得控制量是稳定的。2、计

算实例控制对象模型：)()()()()()(kezCkuzBkyzA)7.0()(19.0)()7.0)(1()(zzzCzzBzzzA已知22)(kEe要求：计算最小方差控制。BAddegdegB的零点在单位圆外，取9.0)(,19.0)

(,1*zzBzzBB11degdegBdF解：<1>控制器的设计11degdegAG于是，设1010)(,)(gzgzGfzfzF代入Diophantine方程（待定系数法）：)()()()()()(*1zGz

BzFzAzBzCzd))(19.0())(7.0)(1()9.0)(7.0(1010gzgzfzfzzzzz通过系数比较，得到7.0,1,1,11010ggff即7.0)(,1)(

zzGzzF故最小方差控制为：17.0)()()()(zzzFzBzGzD<2>计算输出方差：y(k)u(k)e(k)++)()(zAzC)()()(zFzBzG)()(zAzB求出e(k)到y(k)的闭环传递函数H(z):*1*1][)(

BzFCBzBFCBGBAFBFCBzHdd代入具体参数，得到从而111111()()()()()0.910.90.11()()()10.9zzykHzekekekzzzekekwkz其中110.1()()10.9zwkekz

9.01)(zzzH即()0.9(1)0.1(1)wkwkek由此可以看出，e(k)与不相关，从而有()wk222()()()EykEekEwk22222()0.9(1)0.1(1)EwkEwkEek假设系统处于平衡状态，有2222(1)()

(1)()EekEekEwkEwk于是22220.1()()10.9EwkEek所以222222053.11920)(9.011.01)(kEekEy3、与极点配置设计法比较闭环系统特征方程（y(k)/e(k)）：

)()()()(*1zCzBzBzzd（39）可见，系统的极点由四部分组成：（1）d-1个原极点；（2）对象中位于单位圆内的零点B+(z)；（3）对象中单位圆外零点关于单位圆周的镜象；（4）C(z)的零点（n个）。1dz)(*zB极点配置设计法中，)()()()(kyzRzSku

（40）则R(z)和S(z)满足如下的Diophantine方程：mAABBSAR0（41）选定,,*10BzACAdm则得到)()(),()()(zGzSzFzBzR（42）可得到)()()()()(kyzFzBzGku（43

）可见，II型最小方差控制器的设计也可以看成是按极点配置设计方法的一个特例。第三节广义最小方差控制最小方差控制的弱点：（1）性能指标中没有对控制量加以限制，因此控制量幅度大；为限制控制量幅度，需要取较大的采样周期，这常使系统的其他性能变差；（2）当对象包含有单位圆外的零点时，需要采取改进措施

，即采用II型最小方差控制。故采用如下更具一般性的二次型函数：)]()([222kukyEJ（1）以此作为性能指标的最优控制称为广义最小方差控制。由上节（1）式与（7）式，得到（2）)()()()()()()()()(

)()()()()()()()()()()()()(1*1*1*1*1*1*1*1*1*1*1*1*1*1*kezAzGdkezFkuzAzBdkezAzGzzFkuzAzBdkezAzCdkuzzAzBdkydd

一、广义最小方差控制的计算同理由上节（1）式，得到)()()()()()()(1*1*1*1*kuzzCzBkyzCzAked（3）（C(z)的零点均在单位圆内）（3）式代入（2）式，得到

：)()()1()1()()()()()()()()()()()()()()()(*****11******************************

****kuCFBkyCGkefdkefdkekuFACABkyCGdkeFkuzGCCABkyCGdkeFkuzCAGBABkyCGdkeFkuzCBkyCAAGkuABdkeFdkydddd（4）根据e(k)是白噪声

序列及的假设，可以得到：22)(kEe2*****221212)()()1()(kuCFBkyCGEffdkEyd（5）若令上式取极小，便可以得到最小方差控制，即)()()()()()()()()(1*1*1*kyzFz

BzGkyzFzBzGku（6）由于系统已处于平衡状态，从而（1）式可以表示为：)()()()(22222dkuEdkEykuEkEyJ（7）（5）式代入（7）式，得到：)(2)()()()()(2)(**********2k

ukuCFBkyCGkukuCFBkyCGEkuJ)()()()1(22*****221212kukuCFBkyCGEffJd（8）为使J2最小，

求J2对u(k)的导数。（9）)1()()()()(1)1)(()()()()('101'1011)1(1111101*1*1*kubkubkuzbbkuzc

zczfzfzbzbbkuzCzFzBnnddmm（10）由于故0*****)()()(bkuCFBkyCGku（11）（11）式代入（9）式，得到)(2)()(2)(0*****2kubkuCF

BkyCGEkuJ（12）使上式等于零即可求得使J2极小的控制为：)()()()()()()()()()()(01*01*1*1*kyzCbzFzzBzzGkyzCbzFzBzGku（13）

0此即为要求的广义最小方差控制。若令，则上式便变为最小方差控制（时D(z)不抵消控制对象的零点B(z)）。0计算广义最小方差控制的步骤如下：（1）计算和：)(1*zF)(1*zG)()()()(1*1*1*1*zGzzFzAzCd或计算F(z)和

G(z)：)()()()(1zGzFzAzCzd（2）代入广义最小方差控制器公式：)()()()()()(0kyzCbzFzzBzzGku二、计算实例控制对象模型：)()()()()()(kezCkuzBkyzA

)7.0()(19.0)()7.0)(1()(zzzCzzBzzzA已知1.0,)(22kEe要求：计算广义最小方差控制。设计：9.0,1degdeg0bBAd)7.0)(1(7.01)7.0)(1()7.0()()(1

zzzzzzzzAzCzd故有7.0)(,1)(zzGzF广义最小方差控制器为：0()0.7()1.010.922()()()zGzzDzzzBzFzCzb求输出方差)(2kEy（1）求出e(k)到y(k)得闭环传递函数)693.

0)(111.0()913.0(1)(zzzzBDACDABACzH（2）求输出方差cdzzzHzHjkEy2112)()(21)(其中c为单位圆周。)1693.0)(1111.0)(693.0)(111.0()1913.0)(91

3.0()693.0)(111.0()913.0()693.0)(111.0()913.0()()(1111111zzzzzzzzzzzzzzzzzHzH)()(1zHzH在单位圆内有两个极点，即z=-0.111和z=-0.693，利用计算留

数得方法求复变积分，有22693.0111.02069.1)1693.0)(1111.0)(693.0()1913.0)(913.0()1693.0)(1111.0)(693.0()1913.0)(913.0()

(zzzzzzzzzzzzkEy（3）求u(k)的方差：)()()693.0)(111.0()17.0(99.0)()693.0)(111.0()913.0(92

2.001.17.0)()()()()()(kezHzzzzkezzzzzzkezHzDkyzDkuu与前面类似，可以求得：2211264.5)()(21)(cuudzzzHzHjkEu（4）对于I型最小方差控制，有22)(kEyB(z)的

零点在单位圆外，因而从的传递函数是不稳定的，因此：)()(kuke)(2kEu（5）对于II型最小方差控制，有22053.1)(kEy（已求得）)(9.07.0)(9.0117.0)()()()(

)()(kezzkezzzzkezHzDkyzDku利用与前面类似的方法可以求得：2247.14)(kEu（6）结论<1>I型最小方差控制：输出方差最小（），控制量趋于无穷（）；<2>广义最小方差控制：输出方差最大（），控制量

方差较小（）；<3>II型最小方差控制：输出方差较大（），控制量方差较大（）。22069.1264.52053.1247.14三、与极点配置设计法的比较由前述可知，从e(k)到y(k)的闭环传递函数H(z)为：)()()()()()()(1)()()(zDzBzAzCzDzAzBzA

zCzH（14）由于)()()()()('zCzFzzBzzGzD其中0'b（15）（16）将（15）式代入（14）式，有)]()()[()]()()()[()()()()()]()()()[()()()]()()()[()]()(

)()[()()()()()()()()]()()()[()()()()()()()()('''1''''''zAzBzzCzCzFzzBzCzCzAzCzzzBzCzFzzBzCzCzAzGzFzAzzBzCzFzzBzCzGzzBzCzAzFzzBzAzCzFzz

BzCzCzFzzBzzGzBzAzCzHdd（17）由此可见，系统的极点由两部分组成：<1>C(z)的零点。这部分极点将被抵消，因此可看成是观测器的极点；<2>的零点。可以看作是闭环

传递函数要求的极点。)()('zAzBzd与极点配置设计方法的比较：设)()()()(kyzRzSku（17）对比式（15），显然有)()()()()()('zzGzSzCzFzzBzR（18）参考式（17）的

推导过程，可以求得)]()()[()()()()()()()()()()()(''zAzBzzCzzGzBzCzAzFzzBzAzSzBzRzAd（19）对比标准的Diophantine方程，即mAABBSAR0

（20）取)()()(),()(,1)('0zAzBzzAzCzAzBdm则（18）式，即)()()()()()('zzGzSzCzFzzBzR为标准Diophantine的解。因此，广义最小方差控制也可

看作是按极点配置设计方法的一个特例。控制器中加权系数的选择：（1）试凑法：<a>初选<b>设计广义最小方差控制器<c>仿真检验<d>修改<e>最终确定0（2）根据希望的闭环系统极点来确定。闭环系统的极点为：)()()(0zAbzBzzA

dm根据要求的闭环系统的极点来确定上式中的加权系数。)(zc第四节跟踪系统的设计1、跟踪系统控制器的设计y(k)u(k)e(k)++)()(zAzC)(2zD)()(zAzB)(1zD+_r(k)控制器图3跟踪系统的结构

图设计问题：设计出控制器传递函数D1(z)和D2(z)，以使系统具有满意的抗干扰和跟踪性能。设)()()(,)()()(21zRzSzDzRzTzD（1）则控制器方程为：)()()()()()(kyzSkrzTkuzR（2）从而设计问题变为设计R(z)，T(z)和S(z)

。与极点配置设计方法比较，C(z)相当于观测器多项式。按照方式1引入参考输入，即参考输入的引入不影响系统的状态重构或状态估计时，可选：)()(zCzT（3）由前述假设条件，有nzCzT)(deg)(deg（4）D2(z)的设计可直接利用前面关于调节系统

的设计结果，R(z)和S(z)选取如下：)()(),()()(zGzSzBzFzRI型最小方差控制（5）式)()(),()()(zGzSzFzBzRII型最小方差控制（6）式)()()()()()(0zzGzSz

CbzFzzBzR广义最小方差控制（7）式对于式（5）（6），有1)(deg)(degnzSzR（8）而由式（4）可知，degT(z)=n，从而degT(z)>degR(z)，这在物理上不能实现，因而对于I型和II型最小方差控制，按（5）（6）式计算完后，取R(z)=z

R(z)，S(z)=zS(z)；对于广义最小方差控制，计算方法不变，即按（7）式计算。于是，R(z)，S(z)和T(z)计算方法归纳如下：)()()()()()()()()()(0zCbzFzzBzRzFzzBzRzF

zzBzR（9）I型最小方差控制II型最小方差控制广义最小方差控制)()(zzGzS)()(zCzT（10）（11）2、跟踪系统的闭环传递函数由图3，得到)()()()()(21kezHkrzHky（12）其中)()()

()()()()()()(1)()()()(211zSzBzRzAzBzCzDzAzBzAzBzDzH)()()()()()()()()(1)()()(22zSzBzRzAzRzCzDzAzBzAzCzH（13）（14）由前述可知，R(z)和

S(z)满足如下的Diophantine方程：)()()()()()()(0zAzAzBzSzBzRzAm（15）对于I型和II型最小方差控制，这里的R(z)和S(z)相当于前面的zR(z)和zS(z)，考虑到这个变化，结合前面推导的结果得到：

1)()()(zBzBzBI型最小方差控制II型最小方差控制广义最小方差控制（16）)()(0zCzA（17）)()()()(0*zAbzBzzBzzzAdddm（18）I型最小方差控制II型最小方差控制广义最小方差控制将（16）（17）和（18）式代入（13）

（14）式，得到：)()()()()(1)()()()(0*1zAbzBzzBzBzzBzzAzBzBzHdddm（19）I型最小方差控制II型最小方差控制广义最小方差控制)()()(

)()()()()()()()()(00*112zAbzBzzCbzFzzBzBzzFzzFzAzBzRzHdddm（20）I型最小方差控制II型最小方差控制广义最小方差控制第七章结束