模型参考自适应控制

1、10.自适应控制严格地说，实际过程中的控制对象自身与能所处的环境都是十分复杂的，其参数会由于种种外部与内部的原因而发生变化。女口，化学反响过程中的参数随环境温度和湿度的变化而变化外部原因，化学反响速度随催化剂活性的衰减而 变慢内部原因，等等。如果实际控制对象客观存在着较强的不确定，那么， 前面所述的一些基于确定性模型参数来设计控制系统的方法是不适用的。所谓自适应控制是对于系统无法预知的变化，能自动地不断使系统保持所希望的状态。因此，一个自适应控制系统，应能在其运行过程中，通过不断地测取 系统的输入、状态、输出或性能参数，逐渐地了解和掌握对象，然后根据所获得 的过程信息，按一定的设计方法，作出控

2、制决策去修正控制器的结构， 参数或控 制作用，以便在某种意义下，使控制效果达到最优或近似更优。目前比拟成熟的 自适应控制可分为两大类：模型参考自适应控制 Model Referenee Adaptive Control和自校正控制Self-Turning。10.1模型参考自适应控制模型参考自适应控制原理模型参考自适应控制系统的根本结构与图10.1所示：_.可调果婕二1飞5 误 z：1/1*为弋厂前愤训节器一 控对象1r 1 /Jd>储lae1 11厂a_1 _占自适应扒杓10.1模型参考自适应控制系统它由两个环路组成，由控制器和受控对象组成内环，这一局部称之为可调系 统，由参考模型和自适

3、应机构组成外环。 实际上，该系统是在常规的反响控制回 路上再附加一个参考模型和控制器参数的自动调节回路而形成。在该系统中，参考模型的输出或状态相当于给定一个动态性能指标，通常，参考模型是一个响应比拟好的模型，目标信号同时加在可调系统与参考模型上， 通过比拟受控对象与参考模型的输出或状态来得到两者之间的误差信息，按照一定的规律自适应律来修正控制器的参数参数自适应或产生一个辅助输入 信号信号综合自适应，从而使受控制对象的输出尽可能地跟随参考模型的输 出。在这个系统，当受控制对象由于外界或自身的原因系统的特性发生变化时，将导致受控对象输出与参考模型输出间误差的增大。于是，系统的自适应机构再次发生作用

4、调整控制器的参数，使得受控对象的输出再一次趋近于参考模型的输 出即与理想的希望输出相一致。这就是参考模型自适应控制的根本工作原理。模型参考自适应控制设计的核心问题是怎样决定和综合自适应律，有两类方法，一类为参数最优化方法，即利用优化方法寻找一组控制器的最优参数， 使与t o系统有关的某个评价目标，如：J=e2(t)dt，达到最小。另一类方法是基于稳o定性理论的方法，其根本思想是保证控制器参数自适应调节过程是稳定的。如基于Lyapunov稳定性理论的设计方法和基于 Popov超稳定理论的方法。系统设计举例以下通过一个设计举例说明参数最优化设计方法的具体应用。例10.1设一受控系统的开环传递函数为

5、W(s)=s(s 1)其中K可变，要求用一参考模型自适应控制使系统得到较好的输出。解：对于该系统，我们选其控制器为PID控制器，而PID控制器的参数由自 适应机构来调节，参考模型选性能综合指标良好的一个二阶系统：1W(d)=二s 1.414s 1自适应津决定的评价函数取tminj= e2(t)dt ，e(t)为参考模型输出与对象输出的误差。0由于评价函数不能写成PID参数的解析函数形式，因此选用单纯形法做为寻 优方法。参见有关优化设计参考文献。在上述分析与考虑下，可将系统表示具体结构表示如如下图所示图系统仿真结构图为了进展数字仿真连续系统离散相似法，图中的模型是采用典型环节积 分环节，惯性环节

6、型式表示。对于图5-17所示系统，可写出各类环节的差数方程，其中， PID控制器的 差分方程为：TTdX3(k) X3(k 1) KpX5(k) X5(k 1) T X5(k) #(X5(k) 2X5(k 1) x5(k 2)式中，K为比例系数，T为积分时间常数，Td为微分时间常数，它们是自适 应津的调整参数，也是单纯形法的寻优参数。T是采样周期。用MATLABS写的仿真程序如下，chp1O-1 .m和，程序中有关离散相似法 仿真与单纯形法寻优的原理请参见有关参考书，如参考文献。global R I9 L1 L2 L3 x p1 p2 p3global Kp E1 F1 G1 E5 F5 G5

7、 mglobal K1 T4 T5x=zeros(3,7);Q=zeros(1,7);p 仁zeros(1000,1);k=zeros(3,1);N=3;E=0.001;R=1;L=0.012;p2=zeros(1000,1);p3=zeros(1000,1);L1=0.05;L2=5;L3=200;m=0;m n=0;Kp1=zeros(100,1);Ti=zeros(100,1);Td=zeros(100,1);Kp=2;E1=1;F 1=L1;G 1=L1*L1/2;E5=exp(-L1);F5=Kp*(1-E5);G5=Kp*(E5-1)+Kp*L1;I9=0;x(1,1)=1;x(2

8、,1)=0.5;x(3,1)=0.9;for j=2:N+1for i=1:Nif i=j-1 x(i,j)=x(i,1)+L;elsex(i,j)=x(i,1); endendendfor j=1:N+1Q(j)=suba(j);endEx=0;while Ex<1Q0=0;for j=1:N+1 if Q(j)>Q0 Q0=Q(j); j0=j;endendQ仁0;for j=1:N+1if j=j0if Q(j)>Q1 Q1=Q(j); j1=j;endendendQ2=Q0;for j=1:N+1 if Q(j)<Q2 Q2=Q(j); j2=j;endendm

9、n=mn+1;Kp1(m n,1)=x(1,j2); Ti(m n,1)=x(2,j2);Td(m n,1)=x(3,j2);if Q2<1if (Q0-Q2)<E Ex=1;endelseif (Q0-Q2)<(E*Q2)Ex=1;endendfor i=1:Nx(i,N+2)=0;for j=1:N+1x(i,N+2)=x(i,N+2)+x(i,j);endx(i,N+2)=2*(x(i,N+2)-x(i,jO)/N-x(i,jO); end j=N+2;Q(N+2)=suba(j);if Q(N+2)<Q(j1)if (1.5*Q(N+2)-0.5*Q0)<

10、Q2;for i=1:Nx(i,N+4)=1.5*x(i,N+2)-0.5*x(i,j0); end j=N+4;Q(N+4)=suba(j);if Q(N+4)<Q(N+2);for i=1:Nx(i,N+3)=x(i,N+4);endQ(N+3)=Q(N+4);elsefor i=1:Nx(i,N+3)=x(i,N+2);endQ(N+3)=Q(N+2)endelsefor i=1:Nx(i,N+3)=x(i,N+2);endQ(N+3)=Q(N+2);endfor i=1:Nx(i,jO)=x(i,N+3);endQ(j0)=Q(N+3);elseif Q(N+2)>Q(j0

11、)for i=1:Nx(i,N+3)=x(i,N+2);x(i,N+2)=x(i,jO);x(i,j0)=x(i,N+3);endendfor i=1:N x(i,N+3)=0.3*x(i,j0)+0.7*x(i,N+2);endj=N+3;Q(N+3)=suba(j);if Q(N+3)<Q(j1)for i=1:Nx(i,j0)=x(i,N+3);endQ(j0)=Q(N+3);elsefor j=1:N+1for i=1:Nx(i,j)=(x(i,j)+x(i,j2)/2;endendfor j=1:N+1Q(j)=suba(j);endendendend19=1;j=j2;Fin

12、=suba(j)x(1,j2),x(2,j2),x(3,j2)figure(1)m1=1:1:m;m1=m1' p1=p1(1:m,1);p2=p2(1:m,1);p3=p3(1:m,1);plot(m1,p1,'r',m1,p2,'B',m1,p3,'g')gridpausefigure(2)m2=1:1:m n;m2=m2'Kp1=Kp1(1:m n,1);Ti=Ti(1:m n,1);Td=Td(1:m n,1); plot(m2,Kp1,'r',m2,Ti,'B',m2,Td,'g&

13、#39;)gridfunction resulta=suba(j)global R 19 L1 L2 L3 x pl p2 p3global Kp E1 F1 G1 E5 F5 G5 mglobal K1 T4 T5K1=x(1,j);T4=x(2,j);T5=x(3,j);x1=0;x2=0;x3=0;x4=0;x6=0;w仁0;w2=0;w4=0;w5=0;u5=0;t=O;F=O;T6=L1*L2;for i1=1:L3x5=R-x6;x3=x3+K1*(x5-w5+T6*x5/T4+T5*(x5-2*w5+u5)/T6);u5=w5;w5=x5;for i2=1:L2xO=R-1.41

14、4*x1-x2;u2=x1;u4=x4;S=( xO-w1)/L1;x1= E1*x1+F1*x0+G1*v1;w仁 x0;v2=(u2-w2)/L1;x2=E1*x2+F1*u2+G1*v2;w2=u2;v3=(u4-w4)/L1;x4=E1*x4+F1*x3;x6=E5*x6+F5*u4+G5*v3;w4=u4;x7=x2-x6;t=t+L1;F=F+x7*x7*L1;if I9=1m=m+1;p1(m,1)=x2;p2(m,1)=x6;p3(m,1)=x7;endendendresulta=F;程序运行后，绘制出受控对象和参考模型的输出与两者间的误差曲线如图10.3，控制器的三个参数&a

15、mp;,Ti ,Td的变化曲线如下列图，图中可以看到，经 过一段时间的调整kd,Td自动调整后，受控对象输出能很好地接近参考模 型的输出。图10.4 K p,Ti ,T d的变化曲线-0.5 11020040060080010001.5图10.3受控对象和参考模型的输出与两者间的误差曲线读者可以改变系统的某些参数如让程序中的ka=1或ka=3等后再运行程 序观看仿真结果。10.2自校正控制自校正控制技术特别适宜用于结构而参数未知但恒定或缓慢变化的随机系统。由于大多数工业对象正好具有这种特征，因此，自校正控制技术在工业控制过程中得到了广泛的应用。自校正控制系统的根本结构图10.5自校正控制系统的

16、典型结构与模型参考自适应控 制系统一样，自校正控制 系统也由两个环路组成， 典型结构如图10.5所示。 内环与常规反响系统类 似，由对象和控制器组成。 外环由参数估计器和控制 器设计计算机构组成。参 数估计器的功用是根据受控对象的输入与输出信息， 连续不断地估计受控对象的参数，而控制器如此根据参数估计器不断送来的参数估计值， 通过一定的控制算 法，按某一性能指标，不断形成最优控制作用，由于存在着多种参数估计和控制 器设计算法，所以自校正控制的设计方法很多， 其中，以用最小二乘法进展参数 估计，按最小方差来形成控制作用的最小方差自校正控制器最为简单，并获得较多应用。1最小方差自校正控制器以一个单

17、输入、单输出、线性时不变的系统为讨论对象。假设系统经常处在随机扰动作用之下，如不加以控制、随机扰动将使系统输出相对其设定值产生 很大波动。或者说，系统输出的稳态方差将很大。最小方差自校正控制就是根据 给定的对象数学模型，综合一个最优控制津，使系统输出的稳态方差最小。设被控制对象的数学模型可表示为A(q-1)y(k)= (q -d)B(q-1)u(k)+C(q -1)e(k)(10.2-1)式中 A(q-1)=1+aiq-1+、+&q-n-1-1nB(q )=bo+biq +、+ boqb 0 0-1-1-nC(q )=1+Gq +、+ c 1qU(k),y(k)为系统的输入，e(k)为

18、均值为0的随机噪声，d为延时时间次数d 1 对于系统，设输出所跟踪的目标值yr(k)，要求求解u(k)使方差Ey(k)-y r(k) 2最小。如果系统中参数是的，这个冋题就是最小方差控制冋题。 而如果参数未知或不断缓慢变化，要在线递推估计和不断修正模型的参数值 且使系统输出方差保持最小，这就成为最小方差自校正控制器问题。为求解该问题，首先用 A(q-1-1 +an b°q-n并除C q-1=Hcq-1+、q-n，进展d 次除法后，C : q-1可写成。Cq-1=Aq-1R q-1+q-dS(q-1(10.2-2)式中 Rq-1=1+r1q-1+、+Sn-1q-d+1-1-1-n+1L

19、q=So+S1q +、Snq由(10.2-2).(1021)式，有，Cq-1y(k)=A q-1Rq-1y(k)+q -dS(q-1)y(K) =q-dR(q-1)B(q-1)u(k)+q -dS(q-1)y(k)+R(q -1)C(q-1)e(k)-1-1-1C(q )y(k)-y r(k)=R(q )B(q )u(k-d)+S(q-1)y(k-d)-C(q -1)y r(k)+C(q-1)R(q-1)e(k)(10.2-3)上式右边的前三项u()、y( )、yr()是能够测定的，此外因为11 1R(q )B(q )=b o,(b 汁rdq ,、b/dj=B0 B 1、Bn+d-1且令T =

20、s 0,S1,、Sn-1 ,B°,B 1,、,B n+d-1 , c 1,C2,、y(k d)y(k d n 1)u(k d)X k-d=u(k d n 1)yr(k 1)yr(k n)如此式可写成-1t-1_1C(q )y(k)-y r(k)=TX(k-d)-y r(k)+C(q )R(q )e(k)(10.2-4)y(k)-y r(k)=c-1(q-1)tTx(K-d)-y r(k)+R(q -1)e(k)2由(10.2-4)欲使 Ey(k)-y r(k) 最小，须使tTX(k-d)-y r(k)=o(10.2-6)这时 y(k)-y r(k)=R(q - )e(k)(10.2-

21、7)如果tT，从10.2-6丨中求u(k-d)可表示为1 u(k-d)= yr(k)-；xo(k d)Bo1 u(k)= y r(k+d)-J Xo (k) (10.2-8)bo式中，0、X)分别表示从、X中抽去B0和u(k-d)后的矢量。但如果 0、B0是未知的，便不能由10.2-8丨确定u(k),而须回到10.2-6，考虑由X和Yr来求的估计问题，这时可以用递推最小二乘方法来进展，其递推公式为PnXiPnX-M(k)1Xn1PnPn X n1 X n 1Pn 1Pn X n 1可求得的估计值，并由此求得o,如此1u(k)=丄y r(k+d)- T %(k)o由此而获得自校正控制输入例10.

22、2考虑一自校正控制问题，设受控对象为y(k)-0.5y(k-1)-0.14y(k-2)=u(k-2)+0.5u(k-3)+0.2e(k)解：用MATLAB编制设计程序chp10_2.m，程序中设置了一个变化的目标值 yr，程序运行结果如图10.6所示，该仿真结果表示，所设计的最小方差自校正 控制器可使系统在有随机干扰的情况下使系统输出很好地跟踪设定值的变化。% Chp10_2globalq1 pm1 id y1_vr1y1=zeros(26,1);u 1=zeros(26,1);v1=zeros(26,1);r1=zeros(26,1);q1=zeros(26,1);mu1= zeros(10

23、,1);my1=zeros(10,1);m n1=zeros(10,1);p m1= zeros(26,26);sf=zeros(2,1);f=zeros(10,1);py1= zeros(800,1);py2=zeros(800,1);c=1;% In itnf=1;di=1; n=2;m1=1;ld=0.99;um=80;mu=2;md=2;p0=1;for i=1:25pm1(i,i)=p0;endy1_v=0;c on tr=O;l=n+mu+md;f(1:5)=0.5 0.14 1.0 0.5 0.2'% Set poi nt valveR=zeros(800,1);for

24、i=1:200R(i)=20;endfor i=201:400R(i)=40;endfor i=401:650R(i)=30;endfor i=651:800R(i)=20;endfor ij=1:800% System_Modelfor j=9:-1:1my1(j+1)=my1(j);mu1(j+1)=mu1(j);mn 1(j+1)=mn 1(j);endmy1(1)=y1_v;mu1(1)=c ontr;mn 1(1)=nf*randn( 1);y1_v=0;for i=1: ny1_v=y1_v+f(i)*my1(i);endfor i=1:muy1_v=y1_v+f(i+n)*mu1(i+md);endfor i=1:m1y1_v=y1_v+f(i+n+mu)*m n1(i);end% X_vectorfor i=25:-1:1u1(i+1)=u1(i);endu1(1)=mu1(1);for i=1: nr1(