最小方差控制.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,最小方差自校正调节器,自校正调节器,(Self-Tuning Regulator,STR),最早是由,Astrom,和,Wittenmark,于,1973,年首先提出来的,其结构如图,1,所示,.,STR,是以,RLS,参数估计方法在线估计最优预报模型,并在此基础上以输出方差最小为调节指标的一种可以适应参数未知或慢时变的自适应控制系统,.,欲讨论参数未知时能调节系统输出方差至最小的,STR,需先引入参数已知时调节系统输出方差最小的最小方差调节器,.,最小方差调节的基本思想是,:,由于系统中信道存在着,d,步时滞,这就使得当前的控制作用,u,(,k,),要到,d,个采样周期后才能对输出产生影响,.,因此,要获得输出方差最小,就必须对输出量提前,d,步进行预报,然后根据预报值来计算适当的调节作用,u,(,k,).,这样,通过不断的预报和调节,就能始终保持输出量的稳态方差为最小,.,在最小方差调节器的研究中,所讨论的被控系统的模型为,A,(,q,-1,),y,(,k,)=,B,(,q,-1,),u,(,k,-,d,)+,C,(,q,-1,),e,(,k,)(1),其中：,对该系统,有如下假设,:,1.,被控系统时滞时间,d,以及时滞算子,q,-1,的多项式,A,、,B,和,C,的阶次及系数都已知,;,2.,被控系统为最小相位系统,即多项式,B(,q,-1,),的所有零点都在单位圆内,;,3.A(q,-1,),、,C(q,-1,),所有零点都为稳定的,即所有零点都在单位圆内,;,4.,e,(,k,),为零均值白色噪声序列,且,E,e,2,(,k,)=,2,.,1,最小方差预测,设在,k,时刻已观测到输出值,y,(,k,),y,(,k,-1),等，希望由此得到预,测值 。,右边为,u,(,k,),u,(,k,-1),；,e,(,k,+,d,),e,(,k,+,d,-1),e,(,k,+1),e,(,k,),e,(,k,-1),等变量的线性组合，,u,(,k,),u,(,k,-1),为系统输入，可直接测量。,随机变量,e,(,k,),e,(,k,-1),可根据在,k,时刻为止的系统输入输出,值计算得到。,e,(,k+,1),e,(,k+d,),为系统在,k,时刻以后的干扰输入，与直到,k,时刻为止的系统输入输出观测值无关。,为了有效利用直到,k,时刻为止的系统输入输出观测值进行预测，,须将这两类变量区分开。为此可将,C,(,q,-1,)/,A,(,q,-1,),分成两部分：,(2),由式,(1),有：,F,和,G,可通过长除法得到，,F,为商，而,q,-d,G,(,q,-1,),为余因子。,也可通过将式,(3),写成,(3),然后比较两边系数得到。由,(3),可将,(2),右边的噪声项写成：,(4),(5),代入,(2),可得：,而由式,(1),有：,(6),代入式,(6),有：,(7),利用式,(4),可将式,(8),化简为：,(8),(9),记基于,k,时刻的观测值对,y,(,k,+,d,),的预报为：则它是,k,时,刻及以前的输入输出的函数。若对预测的要求是使预测的误差平,方即系统误差的方差为最小，则损失函数可表示为：,上式中,F,(,q,-1,),e,(,k,+,d,),与其它项均不相关，且由于,e,(,k,),为零均值,白噪声序列，式,(10),可写为,(10),与 的选择无关,因此当上式中第,2,项为,0,时，可使,J,最小。因此最小方差预测为：,(11),(12),最小方差预测估计的误差 的方差为,(13),2,最小方差控制,最小方差控制的目的是要确定,u,(,k,),，使得输出的方差为最,小，由于,u,(,k,),最早只能影响到,y,(,k+d,),，因此选择性能指标为,上式可改写为：,(14),(15),预测误差，,e,(,k,+1),e,(,k+d,),的线性组合。,显然，使式,(15),中性能指标取最小值的充要条件是：,(16),因此最小方差控制律为：,此时系统输出的方差为：,(17),(18),由式,(16),可见，最小方差控制律可以通过先求出输出提前,d,步的,预测值 ，然后令 等于理想输出值,y,r,(,这里,y,r,=0),而得到，因此最小方差控制问题可分离成两个问题，一个是预测,问题，另一个是控制问题。,例,1,求解被控系统,(1-1.7,q,-1,+0.7,q,-2,),y,(,k,)=(1+0.5,q,-1,),u,(,k,-,d,)+(1+1.5,q,-1,+0.9,q,-2,),e,(,k,),的最小方差控制,.,解,首先考虑时滞,d,=1,的情况，这时显然有,F,(,q,-1,)=1,。,设,G,(,q,-1,)=,g,0,+,g,1,q,-1,则由式,(4),可得,(1+1.5,q,-1,+0.9,q,-2,)=(1-1.7,q,-1,+0.7,q,-2,)+,q,-1,(,g,0,+,g,1,q,-1,),则由式,(17),可得最小方差控制：,而：,其次考虑时滞,d,=2,的情况，这时设,G,(,q,-1,),与前面一致。而,设,F,(,q,-1,)=1+,f,1,q,-1,则通过比较系数可得,f,1,=3.2,g,0,=5.64,g,1,=-2.24.,最小方差控制：,而：,最小方差调节系统的闭环稳定性质,由被控系统模型,A,y,(,k,)=,B,u,(,k,-,d,)+,C,e,(,k,),和最小方差调节律,u,(,k,)=-,G,/(,BF,),y,(,k,),可得调节系统的闭环框图如图,2,所示,.,由图,2,可以导出最小方差调节系统的闭环方程,因此,当,B,为稳定多项式（即系统,(1),为最小相位系统）时,上式中分子和分母中的多项式,B,可以对消,于是,y,(,k,)=,F,e,(,k,)(20),不难看出,最小方差调节系统的实质,就是利用调节器,(17),的极点去对消被控系统的零点,.,3,自校正调节器,(STR),前面我们讨论了被控系统在参数已知时的随机离散系统的最小方差调节规律,而,STR,主要解决被控系统参数未知或慢时变时的最小方差调节问题,.,对,STR,问题,有,直接法,和,间接法,.,所谓,间接法,即在每一控制,(,采样,),周期先辨识系统模型,然后基于实时辨识模型求解丢番图方程,计算最小方差调节律及相应的在线控制量,.,所谓,直接法,则直接辨识系统的输出预报模型,以避免在每一控制周期求解丢番图方程和计算最小方差调节律,.,由最小方差控制的原理可知，最小方差是通过置输出的,d,步预测值 为,0,而实现的。因此最小方差控制的核心是预测。,令,(,q,-1,)=,G,(,q,-1,),=,0,+,1,q,-1,+.+,n,-,1,q,-(,n,-,1),(,q,-1,)=,B,(,q,-1,),F,(,q,-1,)=,0,+,1,q,-1,+.+,n,+,d,-,1,q,-(,n,+,d,-,1),(21),则最小方差控制律式,(17),可写成：,(,q,-1,),u,(,k,)=-,(,q,-1,),y,(,k,),(22),或：,显然，若能直接估计出参数,i,和,i,，则可由上式立即得到最小方差控制，为了能估计参数,i,和,i,，我们设定一个具有白噪声干扰的预测模型：,y,(,k,+,d,)=,(,q,-1,),y,(,k,)+,(,q,-1,),u,(,k,)+,(,k,+,d,)(24),(23),由式,(12),，系统的,d,步最优预测为,预测误差,因此,(25),(26),(27),式,(27),与设定的预测模型式,(24),一致。根据系统辨识原理可知，,这时采用最小二乘法等即可得到,i,和,i,的无偏估计。,将消失,为了保证预报模型在闭环下的参数可辨识性的要求,可以设定多项式,(,q,-1,),的首项系数,0,为一合理的估计值,0,则可列写出如下自回归方程,y,(,k,+,d,)-,0,u,(,k,)=,T,(,k,),+,(,k,+,d,)(28),其中,自回归方程,(28),的未知参数向量,可由带遗忘因子的渐消记忆,ELS,法或,SA,法来估计。,基于参数向量,的在线估计值,(,k,),和最小方差调节律有如下自校正调节律,上述自校正调节器为保证闭环可辨识性，未辨识参数,0,，而是辨识其估计值。,可以证明，若,0,的估计值满足,:,则上述自校正调节律一样可以收敛。,广义最小方差自校正控制,在前一讲讨论最小方差调节器和,STR,的稳定性时已指出,STR,仅适用于最小相位系统。,再者,最小方差调节器和,STR,中仅考虑使输出的方差为最小,没有考虑使系统的输出跟踪给定的伺服输入项,也没有考虑对控制作用加以约束以避免使控制作用变化剧烈或频繁以降低控制仪表和测量仪表的损耗率,.,而许多实际控制系统的目的是使系统的输出跟踪给定的伺服输入项,并对系统的输入及输入的变化量加以约束以提高系统的稳定性,系统运行的平稳性,以及降低仪表的损耗,.,为了克服最小方差调节器和,STR,的上述缺陷,Clake,和,Gawthrop,等人于,1975,年提出了广义最小方差控制算法,.,这种算法仍然采用二次型的指标函数,但在指标函数中引入了对伺服输入项的跟踪和对控制作用的约束,.,由于引入了对控制作用的约束,不仅限制了控制作用的不适当的变化,同时使得可以通过适当选择对控制项约束的权因子的大小来使得该算法能适用于非最小相位系统,.,考虑到使系统的输出能跟踪给定的伺服输入项,以及对控制作用加以约束以避免使控制作用变化剧烈并使该控制方案能适用于非最小相位系统,最小方差控制器的指标函数为,J,=E,P,(,q,-1,),y,(,k,+,d,)-,R,(,q,-1,),w,(,k,),2,+,Q,(,q,-1,),u,(,k,),2,(31),式中,w,(,k,),为已知的参考输入量,;,P,R,和,Q,分别为对实际输出、参考输入和控制输入的加权多项式。,对指标函数,J,的意义有如下说明,:,该指标函数中的第一项的作用为跟踪伺服输入,w,(,k,).,该指标函数中的第二项的作用为约束控制量,u,(,k,),幅值和变化量,.,加权多项式,Q,的作用为约束控制量的幅值和变化,因此,Q,的选取要兼顾对控制量的幅值和变化量的约束,.,该指标函数中若,Q,=0,即不对控制项加以约束,则由该指标函数所定义的控制问题即为输出跟踪问题,.,若,R,=,Q,=0,即不考虑跟踪伺服输入,则该指标函数所定义的控制问题即为上一讲中讨论的最小方差调节器问题,.,也就是说,最小方差调节器可视为广义最小方差控制器的特例,.,由上一节的结果可求出,y,(,k+d,),的,d,步最优预报为：,(32),预报误差：,(33),将以上两式代入,J,中可得：,不相关于,y,(,k-i,),u,(,k-i,);,i,0,。记,则：,(34),(35),(36),由：,其中,b,0,q,0,分别是多项式,B,(,q,-1,),和,Q,(,q,-1,),中的常数项，可得使,J,为最小的广义最小方差控制,式中：,也可写成：,上式中,d,步最优预测值 由式,(32),给出。,(37),(38),(36),另一种表示方法是直接求变量,的,d,步最优预测 。采用类似于求 的方法。,首先，,其中,n,p,为多项式,P,(,q,-1,),的次数。,(39),(40),(41),类似上一节的式,(6),到式,(8),推导过程可得,与其他项不相关,(42),再类似上一节由式,(10),到式,(12),的推导过程可得,(43),(44),(45),(46),利用 可以将广义最小方差控制律式,(37),写为：,辅助系统方法,(Clark,等,),首先定义如下辅助系统,其中 称为广义输出。,利用 ，立即可以得到,其最优预测误差 。与求解最小方差控制的思想一样，令关于广义输出的,d,步最优预测值,即可得广义最小方差控制律：,与直接通过优化过程而得到的式,(45),完全一致。,(47),(48),(49),将式,(43),代入式,(48),：,(50),式,(50),可简写成：,(51),由 ，立即可得到与式,(49),等价的广义最小方差控制律,(52),也可写成：,(53),将上式代入受控对象模型式,(1),，可得到采用广义最小方差控制后的闭环系统模型：,(54),闭环系统的特征方程：,(55),选择,Q,和,P,，可以配置闭环系统的极点，即,P,Q,的选择将影响到系统的稳定性。当,Q,=0,时，闭环特征方程为,B,(,q,-1,),的零点成为闭环特征方程的根。,因此当受控对象为非最小相位系统时，闭环系统就不稳定。从性能指标式,(31),来看，,Q,=0,相当于在性能指标中不包括对控制作用的约束，广义最小方差控制退化为最小方差控制，因此它不适用于非最小相位系统。而适当地 选择,Q,不仅可使广义最小方差控制适用于非最小相位的受控对象，保证闭环系统稳定，还可保证控制,u,(,k,),不至于过大。,例,2,已知受控对象为,(1-1.5,q,-1,+0.7,q,-2,),y,(,k,)=,q,-1,(1+0.5,q,-1,),u,(,k,)+(1-0.5,q,-1,),e,(,k,),性能指标为,J,=,E,y,(,k,)-,w,(,k,),2,+0.5,u,2,(,k,),。,试确定广义最小方差控制律,解,A,=(1-1.5,q,-1,+0.7,q,-2,),，,C,=(1-0.5,q,-1,),，,d,=1.,将,A,C,d,代入式,(40),(1-0.5,q,-1,)=(1-1.5,q,-1,+0.7,q,-2,),E,(,q,-1,)+,q,-1,D,(,q,-1,),可解得,E,(,q,-1,)=1,，,D,(,q,-1,)=1-0.7,q,-1,代入式,(53),即可得广义最小方差控制律。,(40),广义最小方差自校正控制（两种算法）：,第一种算法根据式,(45),，利用 的最优,d,步预测值来实现。式,(45),中，,R,和,Q,均是加权多项式，一旦性能确定，他们就是确定的，而,w,(,k,),为外加输入，因此只要确定了最优预测值,即可确定控制量,u,(,k,),。将 的,d,步预测方程式,(43),改写成线性回归的形式：,基于参数向量,的在线估计值，得到 ，代入式,(46),即可得到广义最小方差控制律。,