第三篇模型预测控制及其MATLAB实现.pptx

《第三篇模型预测控制及其MATLAB实现.pptx》由会员分享，可在线阅读，更多相关《第三篇模型预测控制及其MATLAB实现.pptx（228页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1v7.1 动态矩阵控制理论v7.2 广义预测控制理论v7.3 预测控制理论分析第7章预测控制理论第1页/共228页2模型预测控制（ModelPredictiveControl：MPC）是20世纪80年代初开始发展起来的一类新型计算机控制算法。该算法直接产生于工业过程控制的实际应用，并在与工业应用的紧密结合中不断完善和成熟。模型预测控制算法由于采用了多步预测、滚动优化和反馈校正等控制策略，因而具有控制效果好、鲁棒性强、对模型精确性要求不高的优点。第2页/共228页3实际中大量的工业生产过程都具有非线性、不确定性和时变的特点，要建立精确的解析模型十分困难，所以经典控制方法如PID控制以及现代控制

2、理论都难以获得良好的控制效果。而模型预测控制具有的优点决定了该方法能够有效地用于复杂工业过程的控制，并且已在石油、化工、冶金、机械等工业部门的过程控制系统中得到了成功的应用。第3页/共228页4目前提出的模型预测控制算法主要有基于非参数模型的模型算法控制（MAC）和动态矩阵控制（DMC），以及基于参数模型的广义预测控制（GPC）和广义预测极点配置控制（GPP）等。其中，模型算法控制采用对象的脉冲响应模型，动态矩阵控制采用对象的阶跃响应模型，这两种模型都具有易于获得的优点；广义预测控制和广义预测极点配置控制是预测控制思想与自适应控制的结合，采用CARIMA模型（受控自回归积分滑动平均模型），具有

3、参数数目少并能够在线估计的优点，并且广义预测极点配置控制进一步采用极点配置技术，提高了预测控制系统的闭环稳定性和鲁棒性。第4页/共228页57.1 动态矩阵控制理论动态矩阵控制是一种基于计算机控制的技术，它是一种增量算法，并基于系统的阶跃响应，它适用于稳定的线性系统，系统的动态特性中具有纯滞后或非最小相位特性都不影响该算法的直接应用。由于它直接以对象的阶跃响应离散系数为模型,从而避免了通常的传递函数或状态空间方程模型参数的辩识，采用多步预估技术从而能有效地解决时延过程问题，按使预估输出与给定值偏差最小的二次性能指标实施控制，因此是一种最优控制技术，动态矩阵控制算法的控制结构主要由预测模型、滚动

4、优化和误差校正及闭环控制形式构成。第5页/共228页6图7-1单位阶跃响应曲线预测模型从被控对象的阶跃响应出发，对象动态特性用一系列动态系数即单位阶跃响应在采样时刻的值来描述，p称为模型时域长度，ap是足够接近稳态值的系数。第6页/共228页7根据线性系统的比例和叠加性质(系数不变原理),若在某个时刻k-i(k=i)输入u(k-i),则对输出y(k)的贡献为:(7-1)若在所有时刻同时有输入,则跟据叠加原理有(7-2)利用上式容易得到y(k+j的n步预估(np)为:(7-3)第7页/共228页8由于只有过去的控制输入是已知的,因此在利用动态模型作预估时有必要把过去的输入对未来的输出贡献分离出来

5、,上式可写为(7-4)上式右端的后二项即为过去输入对输出n步预估,记为(7-5)将式(3-4)写成矩阵形式(7-6)第8页/共228页9为增加系统的动态稳定性和控制输入的可实现性，以及减少计算量可将向量减少为m维(m=n=m)个时刻的输出值尽可能接近期望值。为简单起见取控制加权系数 (常数)第11页/共228页12若令则式(7-9)可表示为(7-10)式中w(k+j)称为期望输出序列值,在预测控制这类算法中,要求闭环响应沿着一条指定的、平滑的曲线到达新的稳定值,以提高系统的鲁棒性.一般取其中为柔化系数 ；y(k)为系统实测输出值；yr为系统的给定值。第12页/共228页13用Y的最优预测值代替

6、Y,即将式(3-8)代入式(3-10)中并令得(7-11)式(7-11)与实际检测值无关，是算法的开环控制形式。由于模型误差,弱非线性特性等影响,开环控制式(7-11),不能紧密跟踪期望值,若等到经过m个时刻后,再重复式(7-11),必然造成较大的偏差,更不能抑制系统受到的扰动。故必须采用闭环控制算式,即仅将计算出来的m个控制增量的第一个值付诸实施,即现时的控制增量为(7-12)式中；第13页/共228页14如果A,都已确定,d 可事先离线解出,在线计算u(k)只需完成两个矢量的点积即可。可见，预测控制的控制策略是在实施了u(k)之后,采集k+1时刻的输出数据，进行新的预测、校正、优化，从而避

7、免在等待m拍控制输入完毕期间,由于干扰等影响造成失控。因此优化过程不是一次离线进行，而是反复在线进行的，其优化目标也是随时间推移的，即在每一时刻都提出一个立足于该时刻的局部优化目标，而不是采用不变的全局优化目标。第14页/共228页15误差校正由于每次实施控制，只采用了第一个控制增量u(k)，故对未来时刻的输出可用下式预测。(7-13)式中表示在t=kT时刻预测的有 作用时的未来p个时刻的系统输出；表示在t=kT时刻预测的无 作用时的未来p个时刻的系统输出；为单位阶跃响应在采样时刻的值。第15页/共228页16由于对象及环境的不确定性,在k时刻实施控制作用后,在k+1时刻的实际输出y(k+1)

8、与预测的输出未见得相等,这就需要构成预测误差并用此误差加权后修正对未来其它时刻的预测即(7-14)式中为t=(k+1)T时刻经误差校正后所预测的t=(k+1)T时刻的系统输出；为误差校正矢量，。第16页/共228页17经校正后的作为下一时刻的预测初值，由于在t=(k+1)T时刻的预测初值,应预测t=(k+2)T,(k+p+1)T时刻的输出值,故令(7-15)由式(3-14),式(3-15)得下一时刻的预测初值为(7-16)这一修正的引入,也使系统成为一个闭环负反馈系统,对提高系统的性能起了很大作用。由此可见,动态矩阵控制是由预测模型,控制器和校正器三部分组成的,模型的功能在于预测未来的输出值,

9、控制器则决定了系统输出的动态特性,而校正器则只有当预测误差存在时才起作用。第17页/共228页187.2 广义预测控制理论十多年来产生了许多自校正器,都成功地用于实际过程,但是对变时延,变阶次与变参数过程,控制效果不好。因此研制具有鲁棒性的自校正器成为人们关注的问题。Richalet等人提出了大范围预测概念,在此基础上,Clarke等人提出了广义预测自校正器,该算法以CARIMA模型为基础,采用了长时段的优化性能指标,结合辨识和自校正机制,具有较强的鲁棒性,模型要求低等特点,并有广泛的适用范围。第18页/共228页19这个算法可克服广义最小方差（需要试凑控制量的加权系数)、极点配置(对阶的不确

10、定性十分敏感)等自适应算法中存在的缺点,近年来,它在国内外控制理论界已引起了广泛的重视,GPC法可看成是迄今所知的自校正控制方法中最为接近具有鲁棒性的一种。广义预测控制是一种新的远程预测控制方法，概括起来具有以下特点基于CARIMA模型目标函数中对控制增量加权的考虑利用输出的远程预报控制时域长度概念的引入丢番图方程的递推求解第19页/共228页20预测模型在预测控制理论中，需要有一个描述系统动态行为的基础模型，称为预测模型。它应具有预测的功能，即能够根据系统的历史数据和未来的输入，预测系统未来输出值。采用CARIMA模型作为预测模型，模型CARIMA是ContrlledAuto-Regress

11、iveIntegratedMoving-Average的缩写，可以译为“受控自回归积分滑动平均模型”，这个模型可以写成(7-17)其中A(z-1),B(z-1),C(z-1)分别是n,m和n阶的 的多项式,=；y(k),u(k)和 分别表示输出、输入和均值为零的白噪声序列，如果系统时滞大于零；B(z-1)多项式领头的一项或几项系数等于零。第20页/共228页21Clarke等人在推导广义预测控制时,为了简单起见，令C(z-1)=1CARIMA模型具有下述优点模型中的积分作用可消除余差；可适用于一类非平稳随机噪声过程；可以描述能用ARMAX模型描述的过程；在大多数情况下,CARIMA模型比ARM

12、AX模型辩识效果更好；用CARIMA模型导出的控制器对未建模动态具有较好的鲁棒性。第21页/共228页22滚动优化1.目标函数为增强系统的鲁棒性(Robustness),在目标函数中考虑了现在时刻的控制u(k)对系统将来时刻的影响,采用下列目标函数(7-18)其中n称为最大预测长度，一般应大于B(z-1)的阶数；m表示控制长度(m=n)；(j)是大于零的控制加权系数。为简单起见取(j)（常数）。第22页/共228页23为了柔化控制作用,控制的目的不是使输出直接跟踪设定值,而是跟踪参考轨线,参考轨线由下式产生(7-19)式中yr,y(k),w(k)分别为设定值、输出和参考轨线，为柔化系数，01。

13、目标函数中后一项的加入,主要用于压制过于剧烈的控制增量,以防止系统超出限制范围或发生剧烈振荡。第23页/共228页24广义预测控制问题,可以归结为求u(k),u(k+1),u(k+m-1)使得目标函数式(7-18)达到最小值的问题,这是一个优化问题。2.输出预测根据预测理论,为了预测超前j步输出,引入丢番图Dioaphantine方程：(7-20)其中第24页/共228页25将式(3-17)两边同时左乘以后与式(7-20)可得时刻k后j步的预测方程为(7-21)令注意到可知的前j项正是系统单位阶跃响应的前j项故有为书写方便将式(7-21)简记为(7-22)第25页/共228页26对未来输出值的

14、预报可忽略未来噪声的影响得到(7-23)上式中包括k时刻的已知量和未知量两部分,用表示已知量,即写成矩阵形式(7-24)其中第26页/共228页27根据式(7-23)可得最优输出预测值为(7-25)式中第27页/共228页283.最优控制律若令则式(7-18)可表示为(7-26)用Y的最优预测值代替Y,即将式(3-25)代入式(3-26)并令可得(7-27)实际控制时,每次仅将第一个分量加入系统即(7-28)式中 为 的第一行第28页/共228页29为了计算简单,通常选mm时,u(k+j-1)=0。这时,u变成了m1矩阵,(GTG+I)则变成mm方阵,降低了维数,减小了计算量。对于阶数较低的简

15、单系统,取m=1,则整个过程将不包括任何矩阵运算。与通常的最优控制不同，广义预测控制采用滚动优化，优化目标是随时间推移的。即在每一时刻都提出一个立足于该时刻的局部优化目标，不是采用不变的全局优化目标。因此，优化过程不是一次离线进行，而是反复在线进行的，这种滚动优化目标的局部性，使其在理想条件下，只能得到全局的次优解，然而当模型失配或有时变、非线性及干扰影响时，却能顾及这种不确定性，及时进行弥补，减小偏差，保持实际上的最优。第29页/共228页30反馈校正 在广义预测控制算法推导过程中，虽然没有明显给出反馈或闭环的表示,但它在进行滚动优化时，强调了优化的基点与实际系统一致。这意味着在控制的每一步

16、，都要检测实际输出并与预测值比较，并以此修正预测的不确定性。当实际系统存在非线性、时变、模型失配、干扰等因素时，这种反馈校正就能及时修正预测值，使优化建立在较准确的预测基础上。因此可降低对基础模型的要求，提高控制的鲁棒性，在实际工业应用中，这点具有十分现实的意义。第30页/共228页317.3 预测控制理论分析的性能分析1、引 言 近年来，广义预测控制无论是在工业应用界还是在控制理论界都得到了广泛的重视，然而由于优化的启发性质和算法的复杂性，对于这一算法的理论研究十分困难，本文通过系统的闭环方块图求出了闭环传递函数，并在内模控制结构的基础上，分析了系统的闭环动态特性、稳定性和鲁棒性。第31页/

17、共228页322、GPC系统闭环传递函数1)闭环方块图为推求系统的闭环方块图，须将GPC算法各部分稍加变换重写如下(1）预测模型式(7-17)（7-29）(2）预测向量f式(7-24)（7-30）(3）参考轨迹W式(7-19)（7-31）其中第32页/共228页33(4）控制增量式(7-28)（7-32）由以上4式得闭环方块图为图7-2GPC控制系统闭环方块图第33页/共228页342)控制结构图根据图7-2的闭环方块图，可画出GPC的控制结构图如图7-3所示。图7-3GPC控制结构图第34页/共228页35 图中粗线表示矢量信号流，细线表示标量信号流，由此可见，GPC是由柔化、调节和预测三部

18、分构成的。在每一时刻，给定值序列经过柔化作用后得到的期望输出向量，与预测输出相比较构成偏差向量，偏差向量与动态向量gT点积得到该时刻的控制增量 ，控制增量一方面通过数字积分运算求出控制量u(k)作用于对象，另一方面又与系统输出一起去预测新的系统输出值f。第35页/共228页363)闭环传递函数由图7-2可求出GPC系统的闭环传递函数为（7-33）闭环系统的输出响应为（7-34）由以上可得如下结论（1）GPC不是用控制器的极点去对消对象的零点，因此不存在因对消不精确所带来的不稳定性问题，所以GPC可用于非最小相位系统。第36页/共228页37（2）闭环传递函数中不包括，只要是稳定的多项式，至于的

19、精度，它只影响跟踪性能而不影响闭环稳定性和鲁棒性。（3）当时,而,由式（7-20）知即,从而代入式(7-34)可得，即系统在稳态时无差跟踪设定值，既使在有阶跃扰动情况下也是无差跟踪，这里由CARIMA模型内部积分作用所决定的。（4）GPC系统闭环传递函数形式较为复杂(涉及到向量F和H),不能明显看出GPC的设计参数对系统性能的关系，不便进行系统分析，为此将GPC算法结构变为内模结构形式，以便借助于内模原理研究系统的鲁棒性、稳定性和其它特性。第37页/共228页383.广义预测控制的内模控制描述1)GPC内模结构的推导过程由图7-2变换依次可得图7-4(a)、(b)和(c)。(a)第38页/共2

20、28页39(c)图7-4GPC系统内模结构的推导过程(b)第39页/共228页40最后由图7-4(c)可得GPC的内模结构图为图7-5GPC系统内模结构图第40页/共228页41图中被控过程的实际模型；被控过程的预测模型；为控制器；为滤波器；总的扰动量；为模型误差信号。其中（7-35）第41页/共228页422)GPC的闭环稳定性由系统的内模结构图7-5可得系统的控制量和输出量分别为:（7-36）（7-37）系统稳定的充要条件是闭环特征方程的零点位于单位园内第42页/共228页43由式（7-36），式（7-37）得特征方程（7-38）（7-39）当模型匹配时，由式（7-38），式（7-39）得

21、（7-40）（7-41）第43页/共228页44由此可得结论：（1）若对象稳定，控制器稳定，则闭环系统稳定；（2）既使控制器具有非线性特性，但只要保证输入输出的稳定性，对稳定的受控对象，一定能得到稳定的闭环响应，也就是说若系统输入u(k)受到约束时,不会影响整个系统的稳定性，这一优点深受广大实际工程人员的欢迎；（3）若对象开环不稳定，除非有准确的极点对消，否则不能得到稳定解。第44页/共228页45以上结论（3）说明,如果对消不完全,不稳定因素总会被激励出来。GPC解决这个问题的措施在于控制时域长度概念的引入，即经过一段时间m后,令控制增量为零，这相当在目标函数中对控制增量作无穷大加权,正是这

22、个无穷大加权，才使不稳定对消因子的作用受到抑制第45页/共228页463)GPC系统的鲁棒性这里我们所谈的鲁棒性是指当模型失配时，系统仍能维持一定的稳定性的程度。在讨论基于模型CARIMA的广义预测控制时，由于我们采用了下列措施，而增强了系统的鲁棒性。（1）参考轨迹的引入为了使控制过程平稳，我们不要求系统输出y(k)直接跟踪设定值yr，而使对象输出y(k)沿着参考轨迹到达给定值yr。由式(7-19)容易看出,小w(k)很快趋向yr，这时控制系统跟踪的快速性好，但鲁棒性差，增大 ，w(k)跟踪yr的过程较长，这时系统跟踪的快速性变差而鲁棒性提高，在实际应用中，值的选择是控制系统的快速性要求与鲁棒

23、性要求的折衷。第46页/共228页47（2）在目标函数中考虑了现在时刻的控制u(k)对系统将来输出的影响广义预测控制问题,可以归结为求,使得目标函数式(7-18)达到最小值的问题，这是一个最优化问题。下面我们借助内模原理对这些措施作一分析。a)由式(7-38)的特征方程可见当 时，通过对 中参数的调整,可使特征方程的根位于单位圆内。第47页/共228页48b)由图7-5得，在干扰d(k)作用下内模控制反馈为（7-42）当时，式（3-42）变为（7-43）由此可知，包含着模型失配的信息，改变的特性，就能获得不同的鲁棒性，而改变 的特性是由选择不同的 来实现的。第48页/共228页49由此可见，滤

24、波器在这里是一个关键因素，中的参数能直接对系统的鲁棒性进行调整，若有模型失配，则的作用将使失配的影响减少到尽量小的程度，而且还能镇定由于模型失配可能带来的不稳定性。若无模型失配,则又可补偿一类动态干扰,若存在模型失配，由式（7-43）直观地可知，要减小失配的影响，可使减小，而由式（7-35）知,增加柔化因子可使减小，故增大柔化作用（增大），可提高系统的鲁棒性。鲁棒性是预测控制突出的优点之一，现有的文献表明，预测控制在工业过程控制中获得成功的应用，首先应归功于预测控制的鲁棒性。第49页/共228页50与控制规律的等价性证明 1.的最优控制律由7.2知，当系统采用CARIMA模型式(7-17)对目

25、标函数式(7-18)的最优控制律为式(7-27)第50页/共228页51其中 矩阵G中的元素 为系统单位阶跃响应的前n项。第51页/共228页522、的最优控制律由7.1知，若已知对象的单位阶跃响应,则对用式(7-9)所表示的目标函数的最优控制律为式(7-11)其中矩阵A中的元素 为系统单位阶跃响应的前n项第52页/共228页533.与的等价性由上可得如下结论()中G矩阵和中A矩阵中的元素同为系统单位阶跃响应的前n项，故对同一系统应有A=G。()中的U与中的U完全一样，同为控制量序列。()中的W与中的W一样，同为柔化作用后的输出跟踪序列第53页/共228页54()式(7-27)的导出是基于由式

26、(7-17)表述的CARIMA参数模型而式(7-11)导出是基于由式(7-17)表述的同一系统的非参数模型。众所周知，对同一线性系统，在同一目标函数下最优控制解具有唯一性。由式(7-18)与式(7-9)相同，知式(7-27)和式(7-11)完全等价。()由的控制规律与的控制规律完全等价知,中的f向量相当于中的Y0向量，而Y0的物理意义是明确的，它是在k时刻预测的未来n个时刻未加控制增量u(k)的系统输出量。从而可知f就表示k时刻基于以往数据对未来输出的预测。第54页/共228页55与的比较 广义预测控制汲取了动态矩阵控制的多步预测、滚动优化策略,同时又保持了自校正控制器的优点。把它与动态矩阵控

27、制作一比较，对于理解其控制机理及优良性能是十分有益的。第55页/共228页561 预测模型采用了CARIMA模型，可用来描述包括不稳定系统在内的任意对象,而采用有限卷积模型作为预测模型，它只能适用于渐近稳定的对象，若不加以修改，适用范围是有局限性的。2 控制机理由于优化策略的一致，两者的控制机理是相同的，系统的动态响应不仅取决于对象参数及优化设计参数，同时也取决于参考轨迹。第56页/共228页573 反馈校正在模型失配或存在干扰时，和都可以通过滤波器F(z-1)抑制干扰或保持闭环稳定性，采用了误差校正修正预测值的策略，主要是通过选择误差校正系数，或者增加滤波器的零点抑制干扰，或者增加其极点改善

28、鲁棒性，但因两者的设计用了同一设计参数，往往难以兼顾。而在中，通过增加滤波器的零点抑制干扰，而对模型失配，主要是通过模型在线辨识和自校正来纠正的。因此，采用了不同的反馈机制分别对付干扰和模型失配，综合的控制效果会更好些。第57页/共228页58 由此可见，广义预测控制汲取了动态矩阵控制的优化策略，而在预测模型，反馈机制方面都保留了自校正控制的优点，它依靠多步预测及滚动优化获取良好的动态性能，利用在线辩识与校正增强控制系统的鲁棒性，而以反馈环节有力地抑止了干扰。第58页/共228页59小结本章首先系统地阐明了广义预测控制与动态矩阵控制的原理；然后在内模结构的基础上对广义预测控制进行了理论分析；最

29、后对广义预测控制与动态矩阵控制规律进行了等价性的证明。第59页/共228页60v8.1系统模型辨识函数v8.2系统模型建立与转换函数v8.3基于阶跃响应模型的控制器设计与仿真函数v8.4基于状态空间模型的预测控制器设计函数v8.5系统分析与绘图函数v8.6通用功能函数第8章MATLAB预测控制工具箱函数第60页/共228页61MATLAB的模型预测控制工具箱提供了一系列用于模型预测控制的分析、设计、和仿真的函数。这些函数的类型主要有：（1）系统模型辨识函数主要功能包括通过多变量线性回归方法计算MISO脉冲响应模型和阶跃响应模型、对测量数据的尺度化等；（2）模型建立和转换函数主要功能包括建立模型

30、预测工具箱使用的MPC状态空间模型以及状态空间模型与MPC状态空间模型、阶跃响应模型、脉冲响应模型之间的转换；第61页/共228页62（3）模型预测控制器设计和仿真工具一分为面向阶跃响应模型的预测控制器设计与仿真函数和面向MPC状态空间模型的设计和仿真函数两类；（4）系统分析工具一包括计算模型预测控制系统频率响应、极点和奇异值的有关函数；（5）其他功能函数一包括绘图和矩阵计算函数等。本节将对模型预测工具箱的主要函数的原理和使用方法按以上分类进行介绍，并给出若干设计应用例子。第62页/共228页638.1 系统模型辨识函数为进行模型预测控制器设计，需要根据系统的输入输出数据建立开环系统的脉冲响应

31、模型或阶跃响应模型，即进行系统模型的辨识。MATLAB模型预测工具箱提供的模型辨识函数如表3-1所示。第63页/共228页64表8-1系统模型辨识函数函数名功能autosc()矩阵或向量的自动归一化scal()根据指定的均值和标准差归一化矩阵rescal()由归一化的数据生成原数据wrtreg()生成用于线性回归计算的数据矩阵mlr()利用多变量线性回归计算MISO脉冲响应模型plsr()利用部分最小二乘回归方法计算MISO脉冲响应模型imp2step()由MISO脉冲响应模型生成MIMO阶跃响应模型validmod()利用新的数据检验MISO脉冲响应模型第64页/共228页65数据向量或矩阵

32、的归一化获得了系统输入输出的原始数据后，为进行参数估计和模型辨识，往往要求对数据进行归一化处理。模型预测工具箱提供的有关数据归一化处理的函数有：autosc()、scal()和，rescal()，这三个函数都是根据数据的均值和标准差来对数据进行归一化或反归一化处理。设有数据向量：其均值和标准差分别为进行归一化处理后得到的向量为第65页/共228页661.矩阵或向量的自动归一化函数autosc()该函数的调用格式为：ax,mx,stdx=autosc(x)例：x=1112;ax,mx,stdx=autosc(x)输出结果：ax=-0.5000-0.5000-0.50001.5000mx=1.25

33、00stdx=0.5000第66页/共228页672.根据指定的均值和标准差进行矩阵或向量的归一化函数scal()该函数的调用格式为：sx=scal(x,mx)sx=scal(x,mx,stdx)例：x=12;21;sx=scal(x,11,0.50.5)输出结果为：sx=0220第67页/共228页683.由归一化的矩阵或向量计算原矩阵或向量（反归一化）函数rescal()该函数的调用格式为：rx=resca(x,mx)rx=rescal(x,mx,stdx)例：x=02;20;rx=rescal(x,11,0.50.5)输出结果：rx=1221第68页/共228页69基于线性回归方法的脉冲

34、响应模型辨识在获得系统的输入输出数据之后，可以采用最小二乘方或部分最小二乘方等线性回归方法来计算系统的脉冲响应模型的系数。在进行线性回归计算之前，需要对原始数据进行预处理，函数wrtreg()用于完成这一预处理过程。函数mlr()和plsr()分别完成基于多变量最小二乘方和部分最小二乘方方法的脉冲响应模型辨识。第69页/共228页701.生成用于线性回归计算的输入输出数据矩阵函数wrtreg()该函数的调用格式为：krey,yrey=wrtrey(x,y,n)式中x为输入数据矩阵；y为输出数据向量；n为脉冲响应模型系数的个数；xreg为预处理后的输入数据矩阵；greg为预处理后的输出数据向量。

35、第70页/共228页712.基于多变量最小二乘方的脉冲响应模型辨识函数mlr()该函数的调用格式为：theta,yres=mlr(xreg,yreg,ninput)theta,yres=mlr(xreg,yreg,ninput,plotopt,wtheta,wdeltheta)式中xreg为预处理后的输入数据矩阵；yreg为预处理后的输出数据向量；ninput为输入变量的个数；plotopt为绘图选项：plotopt=0，缺省值，不绘制图形；plotopt=l，绘制实际输出和预测输出；plotopt=2，绘制实 际 输 出、预 测 输 出 以 及 输 出 误 差；wtheta，wdelthet

36、a为最小二乘方的加权向量，缺省值均为0；theta为脉冲响应模型的系数矩阵。该矩阵的每一列对应一个输入到输出的脉冲响应模型系数向量；yres为预测误差向量。第71页/共228页72例8-1 考虑一个单输入单输出的对象，其传递函数为解采用下面的MATLAB程序对该对象进行脉冲响应模型辨识，系统实际输出和模型的预测输出以及预测误差曲线如图8-1所示。Ex8_1第72页/共228页73图8-1脉冲响应模型预测输出与预测误差曲线第73页/共228页74例8-2 考虑一个双输入单输出的对象，其传递函数矩阵为采样时间为7秒。解 采用下面的MATLAB程序对该对象进行脉冲响应模型辨识，系统实际输出和模型的预

37、测输出以及预测误差曲线如图3-7所示。Ex8_2第74页/共228页75图8-2脉冲响应模型预测输出与预测误差曲线第75页/共228页763.基于部分最小二乘方（PLS）方法的脉冲响应模型辨识函数plsr()该函数的调用格式为：theta,w,cw,ssqdif,yres=plsr(xreg,yreg,ninput,lv,plotopt)第76页/共228页77例8-3 考虑一个单输入单输出的对象，其传递函数为解 采用如下的MATLAB程序对该系统进行模型响应模型辨识，模型的预测输出和预测误差如图3-8所示。ex8_3第77页/共228页78图8-3基于PLS线性回归的模型预测输出与预测误差第

38、78页/共228页79脉冲响应模型转换为阶跃响应模型在模型预测控制工具箱中，模型预测控制器的设计主要基于系统阶跃响应模型。因此在完成系统的脉冲响应模型辨识后，往往需要转换为等价的阶跃响应模型。函数imp2step()用于完成这一转换过程，该函数的调用格式为plant=imp2step(delt,nout,theta1,theta2,)第79页/共228页80例8-4 考虑一个单输入单输出的对象，其传递函数为解 采用如下的MATLAB语句对该系统进行脉冲响应模型辨识，并转换为阶跃响应模型，系统的阶跃响应曲线如图3-9所示。ex8_4第80页/共228页81图8-4系统阶跃响应的预测曲线第81页/

39、共228页82模型的校验在根据系统的输入输出数据完成模型辨识以后，可以进一步利用新的数据对辨识模型进行校验，函数validmod()用于实现这一功能。该函数的调用格式为：yres=validmod(xreg,yreg,theta,plotopt)第82页/共228页838.2 系统模型建立与转换前面讨论了利用系统输入输出数据进行系统模型辨识的有关函数及使用方法，为进行模型预测控制器的设计，需要对系统模型进行进一步的处理和转换。MATLAB的模型预测工具箱中提供了一系列函数完成多种模型转换和复杂系统模型的建立功能。在模型预测工具箱中使用了两种专用的系统模型格式，即MPC状态空间模型和MPC传递函

40、数模型。这两种模型格式分别是状态空间模型和传递函数模型在模型预测工具箱中的特殊表达形式。这两种模型格式可以同时支持连续和离散系统模型的表达，在MPC传递函数模型中还增加了对纯时延的支持。表3-2列出了模型预测工具箱的模型建立与转换函数。第83页/共228页84表8-2模型建立与转换函数函数名功能ss2mod()将通用状态空间模型转换为MPC状态空间模型mod2ss()将MPC状态空间模型转换为通用状态空间模型poly2tfd()将通用传递函数模型转换为MPC传递函数模型tfd2mod()将MPC传递函数模型转换为MPC状态空间模型mod2step()将MPC状态空间模型转换为MPC阶跃响应模型

41、tfd2step()将MPC传递函数模型转换为MPC阶跃响应模型ss2step()将通用状态空间模型转换为MPC阶跃响应模型mod2mod()改变MPC状态空间模型的采样周期th2mod()将Theta格式模型转换为MPC状态空间模型addmod()将两个开环MPC模型连接构成闭环模型，使其中一个模型输出叠加到另一个模型输入addmd()向MPC对象添加一个或多个测量扰动addumd()向MPC对象添加一个或多个未测量扰动paramod()将两个MPC系统模型并联sermod()将两个MPC系统模型串联appmod()将两个MPC系统模型构成增广系统模型第84页/共228页85模型转换在MAT

42、LAB模型预测控制工具箱中支持多种系统模型格式，这些模型格式包括：通用状态空间模型；通用传递函数模型；MPC阶跃响应模型；MPC状态空间模型；MPC传递函数模型。第85页/共228页86 在上述模型格式中，前两种模型格式是MATLAB通用的模型格式，在其他控制类工具箱中如控制系统工具箱、鲁棒控制工具箱等都予以支持；而后三种模型格式则是模型预测工具箱特有的。其中MPC状态空间模型和MPC传递函数模型是通用的状态空间模型和传递函数模型在模型预测工具箱中采用的增广格式，为了方便计算。模型预测工具箱中提供了若干函数，完成上述模型格式之间转换功能。下面对这些函数的用法加以介绍。第86页/共228页871

43、.通用状态空间模型与MPC状态空间模型之间的转换MPC状态空间模型在通用状态空间模型的基础上增加了对系统输入输出扰动和采样周期的描述信息，函数ss2mod()和mod2ss()实现这两种模型格式之间的转换。（1）通用状态空间模型转换为MPC状态空间模型函数ss2mod()该函数的调用格式为：pmod=ss2mod(A,B,C,D)pmod=ss2mod(A,B,C,D,minfo)pmod=ss2mod(A,B,C,D,minfo,x0,u0,y0,f0)第87页/共228页88例8-5 下面的MATLAB语句计算具有如下传递函数的系统的MPC状态空间模型格式ex8_5输出结果为A=-2-2-

44、1100010B=100第88页/共228页89C=131D=0pmod=1310010NaN-2-2-1100010000000100000131000第89页/共228页90（2）MPC状态空间模型转换为通用状态空间模型函数mod2ss()该函数的调用格式为：A,B,C,D=mod2ss(pmod)A,B,C,D,minfo=mod2ss(pmod)A,B,C,D,minfo,x0,u0,y0,f0=mod2ss(pmod)第90页/共228页912.通用传递函数模型转换为MPC传递函数模型通用传递函数模型与MPC传递函数模型的转换函数ploy2tfd()的调用格式为：g=ploytf2(

45、num,den,delt,delay)第91页/共228页92例8-6 考虑如下的纯时延二阶对象：Ex8_6输出结果为g=01.00001.00001.00004.00004.000000.50000第92页/共228页933.MPC传递函数模型转换为MPC状态空间模型函数tfd2mod()该函数的调用格式为：pmod=tfd2mod(delt,ny,g1,g2,g3,g25)第93页/共228页94例8-7 考虑如下的二阶对象：Ex8_7输出结果为g=011136000mod1=0.10002.00001.0000001.00000NaN1.6892-0.74081.000000001.00

46、000000000.0900-0.08150000第94页/共228页954.MPC阶跃响应模型与其他模型格式之间的转换函数mod2step(),tfd2step()和ss2step()分别用于将MPC状态空间模型、MPC传递函数模型和通用状态空间模型转换为MPC阶跃响应模型。下面对这三个函数的用法进行说明。（1）MPC状态空间模型转换为MPC阶跃响应模型函数mod2step()该函数的调用格式为：plant=mod2step(pmod,tfinal)plant,dplant=mod2step(pmod,tfinal,delt,nout)第95页/共228页96（2）MPC传递函数模型转换为M

47、PC阶跃响应模型函数tfd2step()该函数的调用格式为：plant=tfd2step(tfinal,delt,nout,g1)plant=tfd2step(tfinal,delt,nout,g1,g25)式中tfinal为阶跃响应的截断时间；deft为采样周期；nout为输出稳定性向量，参见函数mod2step()的有关说明；g1,g2，为SISO传递函数，对应多变量系统传递函数矩阵的各个元数按行向量顺序排列构成的向量，其最大个数限制为25；plant为对象的阶跃响应系数矩阵。第96页/共228页97例8-8设系统传递函数为将其转换为阶跃响应模型。解MATLAB命令如下Ex8_8图8-5系

48、统阶跃响应曲线由阶跃响应模型绘制的系统阶跃响应曲线如图8-5所示。第97页/共228页98例8-9考虑如下的双输入双输出系统的传递函数矩阵将多变量系统传递函数模型转换为阶跃响应模型。解MATLAB命令如下Ex8_9第98页/共228页99系统在阶跃输入作用下的响应曲线如图8-6所示。图8-6多变量系统的阶跃响应曲线第99页/共228页100（3）通用状态空间模型转换为MPC阶跃响应模型函数ss2step()该函数的调用格式为：plant=ss2step(A,B,C,D,tfinal）plant=ss2step(A,B,C,D,tfinal,delt1,delt2,nout)第100页/共228

49、页101例8-10 考虑如下的系统状态空间模型：将该模型转换为阶跃响应模型。解利用以下MATLAB命令可绘制如图8-7所示的阶跃响应曲线模型转换，截断时间为5s，阶跃响应的采样周期为0.2。Ex8_10图8-7系统阶跃响应曲线第101页/共228页1025.其他的模型转换函数（1）改变 MPC 状态空间模型的采样周期函数mod2mod()该函数的调用格式为：newmod=mod2mod(oldmod,delt)式中oldmod为离散系统原有的MPC状态空间模型；delt为指定系统新的采样周期；newmod为改变采样周期后的系统MPC状态空间模型。第102页/共228页103例8-11 考虑二阶

50、对象解 下面的MATLAB语句用于计算并绘制采样周期为0.5s时系统的阶跃响应曲线，如图8-11所示。ex8_11第103页/共228页104图8-8二阶对象的阶跃响应第104页/共228页105（2）系统辨识工具箱的Theta格式模型转换为MPCmod格式函数th2mod()该函数的调用格式为：umod=th2mod(th)umod,emod=th2mod(thl,th2,th3,th4,th5,th6,th7,th8)第105页/共228页106模型建立函数1.计算闭环系统模型函数addmod()在模型建立工具中，函数addmod()用于计算闭环系统的模型。该函数的调用格式为：pmod=a