最小方差调节器和自校正调节器.pptx

自适应控制篇自适应控制篇目录(1/2)自适应控制篇自适应控制篇q第10讲 自适应控制概述q第11讲 最优预报和自适应预报q第12讲 最小方差调节器和自校正调节器q第13讲 最小方差控制器与自校正控制器q第14讲 极点配置调节器与极点配置自校正调节器q第15讲 自校正PID调节器q第16讲 多变量自适应控制q第17讲 自适应信号处理与滤波q第18讲 模型参考自适应控制概述自适应控制篇自适应控制篇目录(2/2)自适应控制篇自适应控制篇(续续)q第19讲 模型参考自适应系统的数学模型表示q第20讲 基于李氏稳定性理论的状态空间模型参考自适应控制q第21讲 基于李氏稳定性理论的输入输出方程模型参考自适应控制q第22讲 基于Popov稳定性理论的状态空间模型参考自适应控制q第23讲 神经网络自适应控制第十二讲 最小方差调节器和STR(1/3)第十二讲第十二讲 最小方差调节器和自校正调节器最小方差调节器和自校正调节器q自 校 正 调 节 器(Self-tuning Regulator,STR)最早是由Astrom和Wittenmark于1973年首先提出来的,其结构如图1所示.第十二讲 最小方差调节器和STR(2/3)qSTR是以RLS参数估计方法在线估计最优预报模型,并在此基础上以输出方差最小为调节指标的一种可以适应参数未知或慢时变的自适应控制系统.当被估计参数收敛时,则根据估计参数而推得的输出方差最小调节律将收敛于被控系统参数已知时的输出方差最小调节律.这时,这种调节律就是渐近最优的了.q欲讨论参数未知时能调节系统输出方差至最小的STR,需先引入参数已知时调节系统输出方差最小的最小方差调节器.第十二讲 最小方差调节器和STR(3/3)q最小方差调节的基本思想是:由于系统称中信道存在着d步时滞,这就使得当前的控制作用u(k)要到d个采样周期后才能对输出产生影响.因此,要获得输出方差最小,就必须对输出量提前d步进行预报,然后根据预报值来计算适当的调节作用u(k).这样,通过不断的预报和调节,就能始终保持输出量的稳态方差为最小.q下面,我们将顺序讨论:最小方差调节律最小方差调节律、最小方差调节闭环系统的稳定性问题最小方差调节闭环系统的稳定性问题,STR,以及以及最小方差调节与自校正调节的计算机仿真最小方差调节与自校正调节的计算机仿真.1 最小方差调节器(1/6)1 最小方差调节器最小方差调节器q在最小方差调节器的研究中,所讨论的被控系统的模型为A(z-1)y(k)=B(z-1)u(k-d)+C(z-1)w(k)对该系统,有如下关于其最小方差调节律的定理.q定理定理1 对被控系统A(z-1)y(k)=B(z-1)u(k-d)+C(z-1)w(k),假设1.被控系统时滞时间d以及时滞算子z-1的多项式A、B和C的阶次的上界以及系数都已知;2.被控系统为逆稳定系统,即多项式B(z-1)的所有零点都在单位圆外;1 最小方差调节器(2/6)3.C(z-1)为稳定多项式,即它的所有零点都在单位圆外;4.w(k)为白色噪声序列,且Ew2(k)=2.那么,在最优指标函数J=Ey(k+d)2 (1)下,其最小方差调节律和最小方差调节误差分别为u(k)=-G/(BF)y(k)(2)y(k)=Fw(k)=w(k)+f1w(k-1)+.+fd-1w(k-d+1)(3)其中F和G满足当P(z-1)=1时的丢番图方程,即C=AF+z-dG1 最小方差调节器(3/6)q证明 设y(k+d/k)和y(k+d/k)分别为y(k+d)在k时刻的d步最优预报和最优预报误差.因此,被控系统输出量的方差为 J=Ey(k+d)2=Ey(k+d/k)+y(k+d/k)2 =Ey(k+d/k)2+Ey(k+d/k)2+2Ey(k+d/k)y(k+d/k)最优预报误差y(k+d/k)与最优预报y(k+d/k)统计无关且期望值为零 =Ey(k+d/k)2+Ey(k+d/k)2 Ey(k+d/k)2 (4)要使(4)式所示的输出量的方差为最小,即把上式的不等式取等式即可.因此,令y(k+d/k)=0可求得最优调节律.1 最小方差调节器(4/6)当P=1,j=d时,由第十三讲中的定理定理1可知,输出y(k)的d步最优预报和最优预报误差分别为y(k+d/k)=Gy(k)+BFu(k)/C (5)y(k+d/k)=Fw(k+d)(6)故,系统的最小方差调节律为u(k)=-G/(BF)y(k)(7)此时,最小方差调节误差为y(k)=y(k/k-d)=Fw(k)=w(k)+f1w(k-1)+.+fd-1w(k-d+1)(证毕).1 最小方差调节器(5/6)q对于Astrom的最小方差调节器,有两种实现方法:一为用数字器件实现的传递函数型控制器传递函数型控制器,如另一为可用数字计算机实现的在线递推计算型控制器在线递推计算型控制器,如1 最小方差调节器(6/6)q例1 求解被控系统(1-0.9z-1)y(k)=0.5u(k-2)+(1+0.7z-1)w(k)的最小方差调节律.q解 该算例与第十二讲中的算例为同一被控系统.因此,由定理1和第十一讲中求出的输出y(k)的2步最优预报式,可得如下最小方差调节律传递函数型控制器传递函数型控制器或在线递推计算型控制器在线递推计算型控制器u(k)=-2.88y(k)-1.6u(k-1)1 最小方差调节器(7/6)此时的输出误差的方差为Ey(k+2)2=EFw(k+2)2=E(1+1.6z-1)w(k+2)2=(1+1.62)2=3.5622 最小方差调节系统的闭环稳定性质(1/4)2 最小方差调节系统的闭环稳定性质最小方差调节系统的闭环稳定性质q由被控系统模型Ay(k)=Bu(k-d)+Cw(k)和最小方差调节律u(k)=-G/(BF)y(k)可得调节系统的闭环框图如图2所示.由图2可以导出最小方差调节系统的闭环方程2 最小方差调节系统的闭环稳定性质(2/4)因此,当B为稳定多项式时,上式中分子和分母中的多项式B可以对消,于是利用恒等式(8)可得y(k)=Fw(k)(9)2 最小方差调节系统的闭环稳定性质(3/4)q不难看出,最小方差调节系统的实质,就是利用调节器(2)的极点去对消被控系统的零点.实际上,由于存在建模误差使得系统模型的极点并不完全与实际系统的极点相等(即存在模型误差),且在控制过程中计算工具和计算方法存在计算误差(如舍入误差和截断误差)等原因,要能实现用调节器的极点去精确对消被控系统的零点而不至于产生不稳定的数值计算,这就要求B(z-1)的全部零点都在单位圆外,即要求系统具有逆稳定的性质.2 最小方差调节系统的闭环稳定性质(4/4)否则,闭环系统内部将隐含有不稳定极点的环节使得闭环系统在受到外界扰动时不稳定,或由计算机递推计算控制量u(k)时由于计算误差扩散使得系统的输出y(k)和控制量u(k)不稳定.要求被控系统为逆稳定系统,对许多离散化的工程系统达不到,这就在一定范围内限制了最小方差调节器的应用(对许多离散化过程,将增加临界不稳定的零点).3 STR(1/7)3 自校正调节器自校正调节器(STR)q前面我们讨论了被控系统在参数已知时的随机离散系统的最小方差调节规律,而STR主要解决被控系统参数未知或慢时变时的最小方差调节问题.对STR问题,与自适应预报类似,亦有直接法直接法和间接法间接法.所谓间接法间接法,即在每一控制(采样)周期先系统模型,然后基于实时辨识模型求解丢番图方程,计算最小方差调节律及相应的在线控制量.所谓直接法直接法,则直接辨识系统的输出预报模型,以避免在每一控制周期求解丢番图方程和计算最小方差调节律.3 STR(2/7)下面主要介绍STR的直接法.STR的基本思想如图1所示.下面,将讨论Astrom和Wittenmark最初提出的STR算法.q对被控系统Ay(k)=Bu(k-d)+Cw(k),Astrom和Wittenmark的STR算法考虑的是C(z-1)=1时，即系统所受到的扰动可用白噪声建模的输出调节问题。因此，系统输出y(k)的最优预报为：y(k+d)=y(k+d/k)+Fw(k+d)(10)3 STR(3/7)式中y(k+d/k)为如下最优预报y(k+d/k)=Gy(k)+BFu(k)(11)当被控系统的参数,即多项式A(z-1)和B(z-1)的系数未知时,则可根据(10)式的形式设定一个与k时刻系统输入u(k)和输出y(k)统计独立的随机扰动(k)=F(z-1)w(k).因此,由(10)式和(11)式,有如下预报模型y(k+d)=G(z-1)y(k)+(z-1)u(k)+(k+d)(13)其中(z-1)=B(z-1)F(z-1)=0+1z-1+.+,n=nb+d-1 (14)3 STR(4/7)因此,由定理1可得如下最小方差调节律u(k)=-G/(BF)y(k)其中G和F满足如下丢番图方程1=AF+z-dGq基于预报模型(13),类似于前一讲中讨论的自适应预报,我们可以递推估计预报模型(13)的未知参数.3 STR(5/7)q为了保证预报模型(13)在闭环下的参数可辨识性的要求,可以设定多项式(z-1)的首项系数0为一合理的估计值0,则由式(13)可列写出如下自回归方程y(k+d)-0u(k)=(k)+(k+d)(15)其中q自回归方程(15)的未知参数向量 可由如下带遗忘因子的渐消记忆ELS法和SA法来估计:3 STR(6/7)一、一、渐消记忆渐消记忆ELS法法二、二、SA法法3 STR(7/7)q基于参数向量 的在线估计值(k)和最小方差调节律(12),有如下自校正调节律q上述自校正调节器为保证闭环可辨识性，未辨识参数0，而是辨识其估计值。可以证明，若0的估计值满足:则上述自校正调节律一样可以收敛。3 STR(8/7)q综上所述,STR的基本算法步骤可总结如下:1.确定被控系统模型的结构,并选定预报模型;2.设定递推参数初值(0)和P(-1)(或r(-1)，以及0的估计值;3.采样获取新的观测数据y(k);4.组成观测数据向量(k)5.用渐消记忆的ELS法(或SA法)计算当前参数递推估计值;6.由自校正调节律(21)计算控制作用u(k),并输出控制作用u(k)到被控系统;7.将采样次数加1,转回到第3步继续循环.4 最小方差调节与最小方差调节与自校正调节的计算机仿真自校正调节的计算机仿真(1/1)4 4 最小方差调节与最小方差调节与自校正调节的计算机仿真自校正调节的计算机仿真q下面分别讨论最小方差调节最小方差调节与与自校正调节的计算机仿真自校正调节的计算机仿真4.1 最小方差调节最小方差调节的计算机仿真的计算机仿真(1/7)4.1 最小方差调节的最小方差调节的计算机仿真计算机仿真q下面给出针对随机线性离散系统(XARMA)A(z-1)y(k)=B(z-1)u(k-d)+C(z-1)w(k)的最小方差调节仿真程序伪代码./%第一步第一步 初始化初始化输入系统阶次na,nb和nc,以及时滞d,输入系统模型Az=1 a1 a2 和Bz=b0 b1 b2;输入噪声模型Cz=1 c1 c2 预报模型阶次nf=d-1;ng=max(nc,na+nf)-d;输入噪声w(k)的方差w设定系统变量初始值:yf1:na+1=0;uf1:d+nb+1=0;wf1:nc+1=0;设定预报变量初始值:yc1:d+na+1=0;uc1:d+nf+nb+1=04.1 最小方差调节最小方差调节的计算机仿真的计算机仿真(2/7)/%第二步第二步 求解求解Diophant方程和最优预报模型方程和最优预报模型PCz=Cz;Fz1=1;PCz(nc+2:n+1)=0;if nf0 for k=2:d PCz(i-1:i-1+na)=PCz(i-1:i-1+na)-Fz(i-1)*Az;Fz(i)=PCz(i);endPCz(d:d+na)=PCz(d:d+na)-Fz(d)*Az;Gz=PCz(d+1:ng+d+1);BFz=conv(Bz,Fz);4.1 最小方差调节最小方差调节的计算机仿真的计算机仿真(3/7)/%第三步第三步 最小方差调节仿真最小方差调节仿真for k=1:最大仿真步数/%被控对象模型仿真被控对象模型仿真(产生系统输入输出信号产生系统输入输出信号,即数据即数据)yf2:na+1=yf1:na;uf2:d+nb+1=uf1:d+nb;wf2:nc+1=wf1:nc;wf1=2*w*(rand()-0.5);yf1=-Az2:na+1*yf2:na+1+Bz1:nb+1*ufd+1:d+nb+1+Cz1:nc+1*wf1:nc+1;4.1 最小方差调节最小方差调节的计算机仿真的计算机仿真(4/7)%输出数据检测输出数据检测yc1=yf1;输出在线最小方差调节值y(k)=yc(1)；%计算在线最小方差调节控制量计算在线最小方差调节控制量uc1=(-Gz*yc(1:ng+1)-FBz(2:nb+nf+1)*uc(2:nb+nf+1)/FBz(1);%输出控制变量至被控对象输出控制变量至被控对象uf1=uc1;%数据移位数据移位yc2:d+ng+1=yc1:d+ng;uc2:d+nf+nb+1=uc1:d+nf+nb;4.1 最小方差调节最小方差调节的计算机仿真的计算机仿真(5/7)%计算输出y(k)的最小方差调节的统计分析结果输出y(k)的最小方差=norm(y)2/最大仿真步数4.1 最小方差调节最小方差调节的计算机仿真的计算机仿真(6/7)q仿真示例:模型:最优预报中所介绍的例例1 系统方程(1-0.9z-1)y(k)=0.5u(k-2)+(1+0.7z-1)w(k)白化后的噪声信号w(k)v均值 Ew(k)=0v方差 Ew2(k)=2.o参数可用于控制噪声信号的能量/噪信比.控制要求:试进行输出y(k)的最小方差调节,并对调节结果进行统计分析.4.1 最小方差调节最小方差调节的计算机仿真的计算机仿真(7/7)仿真结果表1 计算机仿真结果输出方差理论值仿真1仿真2仿真3仿真步数50003.563.50213.51113.58394.2 自校正调节自校正调节的计算机仿真的计算机仿真(1/13)4.2 自校正调节自校正调节的的计算机仿真计算机仿真q下面给出针对随机线性离散系统(XARMA)A(z-1)y(k)=B(z-1)u(k-d)+w(k)的自校正调节仿真程序伪代码.4.2 自校正调节自校正调节的计算机仿真的计算机仿真(2/13)/%第一步第一步 初始化初始化输入系统阶次na,nb和nc,以及时滞d,加权因子输入系统模型Az=1 a1 a2 和Bz=b0 b1 b2;预报模型阶次nf=d-1;ng=max(nc,na+nf)-d;辨识参数个数n=ng+nf+nb+1输入噪声w(k)的方差w输入b0的估计值b0g设定系统变量初始值:yf1:na+1=0;uf1:d+nb+1=0;wf1:nc+1=0;设定预报变量初始值:yc1:d+na+1=0;uc1:d+nf+nb+1=0;设定辨识的初始值1:n=0;Pn,n=106*I4.2 自校正调节自校正调节的计算机仿真的计算机仿真(3/13)%第二步第二步 自校正调节仿真自校正调节仿真for k=1:最大仿真步数%被控对象模型仿真被控对象模型仿真(产生系统输入输出信号产生系统输入输出信号,即数据即数据)yf2:na+1=yf1:na;uf2:d+nb+1=uf1:d+nb;wf2:nc+1=wf1:nc;wf1=2*w*(rand()-0.5);yf1=-Az2:na+1*yf2:na+1+Bz1:nb+1*ufd+1:d+nb+1+wf1;4.2 自校正调节自校正调节的计算机仿真的计算机仿真(4/13)%输出数据检测输出数据检测yc1=yf1;输出在线自校正调节值y(k)=yc(1)；%在线辨识自校正调节律在线辨识自校正调节律=yc(d+1:ng+d+1)uc(d+2:d+nf+nb+1);K=P*/(+*P*);=+K*yf(1)-b0g*u(d+1)-*;P=I-K*P/;修正矩阵P；4.2 自校正调节自校正调节的计算机仿真的计算机仿真(5/13)%在线计算自校正调节器控制量在线计算自校正调节器控制量=yc(1:ng+1)uc(2:nf+nb+1);uc(1)=-*/b0g;%输出控制变量至被控对象输出控制变量至被控对象uf1=uc1;%数据移位数据移位；yc2:d+ng+1=yc1:d+ng;uc2:d+nf+nb+1=uc1:d+nf+nb;%计算输出计算输出y(k)的自校正调节的统计分析结果的自校正调节的统计分析结果输出y(k)的最小方差=norm(y)2/最大仿真步数4.2 自校正调节自校正调节的计算机仿真的计算机仿真(6/13)q仿真示例:模型:最优预报中所介绍的例例1 系统方程(1-1.5z-1+0.7z-2)y(k)=(1+0.5z-1)u(k-d)+w(k)白化后的噪声信号w(k)v均值 Ew(k)=0v方差 Ew2(k)=2.o参数可用于控制噪声信号的能量/噪信比.自校正调节要求:试分别对d=1和d=2这2种情形,进行输出y(k)的自校正调节,并进行调节结果统计分析.4.2 自校正调节自校正调节的计算机仿真的计算机仿真(7/13)仿真结果表1 计算机仿真结果输出方差自适应调节律g0/b0 g1/b0FB1/b0 FB2/b0d=1数据组数5000真值11.5-0.70.5b0/b0=1 1.00601.5091-0.69290.5024b0/b0=1.61.01171.5159-0.69120.5001b0/b0=0.61.34201.4671-0.67750.5086d=2数据组数50000真值3.251.55-1.0852.00.75b0/b0=1 3.40611.5447-1.06781.99590.7483b0/b0=1.34.68091.6125-1.05022.06520.7823b0/b0=0.953.30871.5375-1.03931.98840.74264.2 自校正调节自校正调节的计算机仿真的计算机仿真(8/13)q自校正调节的仿真结果如下图所示d=1，b0/b0=1时的自校正调节结果参数估计误差的平方和系统输出4.2 自校正调节自校正调节的计算机仿真的计算机仿真(9/13)d=1，b0/b0=1.6时的自校正调节结果参数估计误差的平方和系统输出4.2 自校正调节自校正调节的计算机仿真的计算机仿真(10/13)d=1，b0/b0=0.6时的自校正调节结果参数估计误差的平方和系统输出4.2 自校正调节自校正调节的计算机仿真的计算机仿真(11/13)d=2，b0/b0=1时的自校正调节结果参数估计误差的平方和系统输出4.2 自校正调节自校正调节的计算机仿真的计算机仿真(12/13)d=2，b0/b0=1.3时的自校正调节结果参数估计误差的平方和系统输出4.2 自校正调节自校正调节的计算机仿真的计算机仿真(13/13)d=2，b0/b0=0.95时的自校正调节结果参数估计误差的平方和系统输出