【文档说明】计算机控制系统第七章课件.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)----控制量，零阶保持器的输出。v1c(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）其中：dTkTeeke

eFTcAdAT0)()(,（7）式中，T为采用周期，是等效的离散白噪声序列，其协方差为：)(kedTBBdeBBeQcTcTATAT220（8）显然有：TBDDDQcT,（9）从而式

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

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

一多项式。0c)(zC（15）其中有：)()(,)()(00keckeczCzCd（16）2202202)()(ckEeckEed（17）设首项系数为，则（14）式变为：)()()()()()()(kezAzCkuzAzBky从而有于是，C(z)也

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

onMovingAverage)随机过程模型。对于式（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越小越好。对于离散系统，性能指标可以表示为：Nk

NkyNEJ121)(1lim（19）对于连续系统，性能指标表示为：201lim()TcNJytdtT（20）具有上述性能指标的最优控制问题为最小方差控制。为对控制量进行限制，可在指标中对控制量进行加权，即)()(222kukyEJ（21）其中为加权系数。上述指

标为更一般形式的二次型性能指标，称为广义最小方差控制。设计问题：设计控制器D(z)，使二次型指标J1或J2最小。一、最优预报控制对象模型：)()()()()()()()()()()()()(1*1*1*1*kezAzCkuzzAzBkezAzCkuzAzBk

yd（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)(dkyu使J1最小的方法：根据k-d及以前的信息最好地估计出k时刻的干扰量，并使)/(ˆdkk)/(ˆ)(dkkdkyu从而)/(ˆ)()()()(dkkkdk

ykkyu（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*kezAzGkefdkefdkedkezAzGzzfzfdkezAzCdkdddd（8）（9）由图1得到：由式

（6），得到：（10）21*1*2212121*1*112)/(ˆ)()()()1()/(ˆ)()()()1()1()()/(ˆ)(kdkkezAzGE

ffkdkkezAzGkefdkefdkeEkdkdkEdd其中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），即221212)1()(dffkEy二、最小方差控

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

dkvkyu由得到)()()()()()()(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)(zzzCzzBzzzA已知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）)()()()()()()(ky

zFzBzGkyzDku由于故)()()()()()()()()(kezCkyzFzBzGzBkyzA（22）（23）于是)()()()]()()()[()()()()(kezHkezGzFzAzBzCzBzFky（24）其中)]()()()[()(

)()()(zGzFzAzBzCzBzFzH（25）由式（16），即)()()()(1zGzFzAzCzd得到)()()(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则得到)()(),()()(zGzS

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

zBkyzA（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、计算实例控制对象模型：)()()()()()(ke

zCkuzBkyzA)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方程（待定系数法）：)()()()()()(*1zGzBzFzAzBzCzd))(19.

0())(7.0)(1()9.0)(7.0(1010gzgzfzfzzzzz通过系数比较，得到7.0,1,1,11010ggff即7.0)(,1)(zzGzzF故最小方差控制为：17.0)()()()(zzzFzBz

GzD<2>计算输出方差：y(k)u(k)e(k)++)()(zAzC)()()(zFzBzG)()(zAzB求出e(k)到y(k)的闭环传递函数H(z):*1*1][)(BzFCBzBFCBGBAFBF

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

ekz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**********************************kuCFBkyCGkefdkefdke

kuFACABkyCGdkeFkuzGCCABkyCGdkeFkuzCAGBABkyCGdkeFkuzCBkyCAAGkuABdkeFdkydddd（4）根据e(k)是白噪声序列及的假设，可以得到：22)(kEe2*****221212)()()1()

(kuCFBkyCGEffdkEyd（5）若令上式取极小，便可以得到最小方差控制，即)()()()()()()()()(1*1*1*kyzFzBzGkyzFzBzGku（

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

)()()()1(22*****221212kukuCFBkyCGEffJd（8）为使J2最小，求J2对u(k)的导数。（9）)1

()()()()(1)1)(()()()()('101'1011)1(1111101*1*1*kubkubkuzbbkuzczczfzfzbzbbkuzCzFzBnnddmm（10）由于故0*****)()()(bkuCFBkyCGku

（11）（11）式代入（9）式，得到)(2)()(2)(0*****2kubkuCFBkyCGEkuJ（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()()()zGzz

DzzzBzFzCzb求输出方差)(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)(913.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(922.001.17.0)()()()()()(kezHzzzzkezzzzzzkezHzDkyzDkuu与前面类似，可以求得：2211264.5)()(21)(cuudz

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

7.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)()()(zDzB

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

)()()[()()()]()()()[()]()()()[()()()()()()()()]()()()[()()()()()()()()('''1''''''zAzBzzCzCzFzzBzCzCzAzCzzzBzCzFzzB

zCzCzAzGzFzAzzBzCzFzzBzCzGzzBzCzAzFzzBzAzCzFzzBzCzCzFzzBzzGzBzAzCzHdd（17）由此可见，系统的极点由两部分组成：<1>C

(z)的零点。这部分极点将被抵消，因此可看成是观测器的极点；<2>的零点。可以看作是闭环传递函数要求的极点。)()('zAzBzd与极点配置设计方法的比较：设)()()()(kyzRzSku（17）对比式（15），显然有)(

)()()()()('zzGzSzCzFzzBzR（18）参考式（17）的推导过程，可以求得)]()()[()()()()()()()()()()()(''zAzBzzCzzGzBzCzAzFzzBzAzSzBzRzAd（19）对比标准的Diop

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

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

极点来确定上式中的加权系数。)(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）式)()()()()()(0zzGzSzCbzFzzBzR广义最小方差控制（7）式

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

计算方法不变，即按（7）式计算。于是，R(z)，S(z)和T(z)计算方法归纳如下：)()()()()()()()()()(0zCbzFzzBzRzFzzBzRzFzzBzR（9）I型最小方差控制

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

zBzRzAzBzCzDzAzBzAzBzDzH)()()()()()()()()(1)()()(22zSzBzRzAzRzCzDzAzBzAzCzH（13）（14）由前述可知，R(z)和S(z)满足如下的Dio

phantine方程：)()()()()()()(0zAzAzBzSzBzRzAm（15）对于I型和II型最小方差控制，这里的R(z)和S(z)相当于前面的zR(z)和zS(z)，考虑到这个变化，结

合前面推导的结果得到：1)()()(zBzBzBI型最小方差控制II型最小方差控制广义最小方差控制（16）)()(0zCzA（17）)()()()(0*zAbzBzzBz

zzAdddm（18）I型最小方差控制II型最小方差控制广义最小方差控制将（16）（17）和（18）式代入（13）（14）式，得到：)()()()()(1)()()()(0*1zAbzBzzBzBzzBzzAzBzBzHdddm（19）I型最小方差控制II

型最小方差控制广义最小方差控制)()()()()()()()()()()()(00*112zAbzBzzCbzFzzBzBzzFzzFzAzBzRzHdddm（20）I型最小方差

控制II型最小方差控制广义最小方差控制第七章结束