广义预测控制在机组控制系统中的仿真研究讲解.doc

广义预测控制在机组控制系统中的仿真研究 摘要：由于单元机组是一个复杂的系统，一般都存在着很大的延迟和惯性，被控对象具有不确定性和强耦合，常规的PID控制器在处理大延迟对象上很难获得较好的控制效果。我们将广义预测算法应用到机组控制中，利用MATLAB仿真，对响应速度和鲁棒性进行分析，通过比较仿真效果表明广义预测控制方法对系统具有良好的适应性和鲁棒性。 关键词：预测控制 控制系统 仿真 0 引言 广义预测控制（GPC）是在最小方差控制的基础上，汲取动态矩阵控制（DMC）和模型算法控制（MAC）中的多步预测优化策略，随着自适应控制的研究而发展起来的一种预测控制方法。此种控制算法的数学模型是通过推导控制律参数而获取的，在理论研究和实际应用中都得到了重视[1]。本文采用广义预测控制算法对500MW机组80%负荷点处定压运行进行仿真研究。 1 多变量广义预测控制算法 预测控制算法是在工业生产过程和相关的理论知识结合的基础上提出来的，与传统的PID控制方法相比较，预测控制算法具有自身的优越性。 广义预测控制具有预测控制算法的基本特征，因而具有优良的控制性能和鲁棒性，使系统的设计更加灵活，被认为是在复杂工业控制领域具有代表性的预测控制算法之一。广义预测控制算法采用的预测模型是具有一定结构和参数的最小化模型，因此对系统输出进行预测时需要已知模型结构，这给预测模型的建立增加了难度，但所需确定的参数较少，减少了计算量。 在GPC中，采用了最小方差控制中所用的受控自回归积分滑动平均模型来描述收到随机干扰的对象[2] （1） 其中 式中，是后移算子，表示后退一个采样周期的相应的量；为差分算子；是一个不相关的随机序列，表示一类随机噪声的影响。A、B、C都是的多项式，其中多项式的若干首项元素可以是零，以表示对象相应的时滞数。为了突出方法原理，这里假设。 式（1）实际上是用Z传递函数给出了对象的描述，由输入u到输出y间的Z传递函数为 （2） 为了利用模型（1）导出j步后输出y(t+j)的预测值，首先考虑下述丢翻图（Diophantine）方程 （3） 其中，，是由和预测长度j唯一确定的 这样，在式（1）两端乘以后可得 利用式（3），可以写出t+j时刻的输出量 注意到，的形式，可以知道与，有关，只与y(t)，y(t-1),...有关，而与有关[3]。由于在t时刻未来的噪声都是未知的，所以对y(t+j)最合适的预测值可以由下式得到 （4） 这样，根据已知的输入输出信息及未来的的输入值，就可以预测对象未来的输出，式（4）就是GPC的预测模型。 在GPC中，t时刻的优化性能指标具有以下形式： （5） 其中，E为数学期望；为对象输出的期望值；和分别为优化时域的识值和终值；NU为控制时域，即在NU步后控制量不再变化： 为控制加权系数，为简化设计一般常假设其为常数。 2 机组受控对象模型 在单元机组协调控制系统中，包括主控系统和子控系统，它们最终的控制对象为锅炉、汽轮机和发电机[4]，它具有大时滞、强耦合、非线性，同时难以平匹配的特点。在整个机组的发电过程中，锅炉和汽轮机的动态特性差别很大，它们对负荷的响应速度一慢一快，所以单元机组内外的两个能量平衡关系互相制约，而它内部的负荷响应性能又与内部参数稳定性存在固有固有的矛盾。因此，在单元机组实施负荷控制的时候，如何协调控制相当重要。为了便于对系统的控制进行分析，在保证满足控制要求的前提下，我们将受控对象的模型简化为两输入两输出的系统，如下图1所示： 图1 单元机组受控对象模型 上图由输入控制和传递函数以及过程输出组成，并且此模型的动态特性由线性常系数多变量传递函数式（6）和（7）来描述[5]： （6） 其中，， （7） 式中是被控制矢量，是控制矢量，是传递函数矩阵。 我们选取的动态模型是某厂300MW直流锅炉燃煤机组在100%负荷点上的数学模型[6]，模型如下[7]： （8） 由图1可以知道单元机组复杂的多输入多输出的控制系统简化成了有汽轮机调节门指令和锅炉的燃烧率指令为输入，实发电功率和机前压力为输出的模型。在公式（8）中，和对应于图1 中的和，和对应于图1中的y1(t)和y2(t)。 3 仿真 控制时域长度M越大，则控制的机动性变强，控制的灵敏度变高，但系统地稳定性和鲁棒性就会下降，为了提高系统的稳定性和鲁棒性，要求M小，M越小，远程跟踪控制能力虽有所削弱，但可得到一个稳定的控制[8]，因此，M的选择应兼顾快速性和稳定性两种，平衡考虑，本文仿真时选取M=6。预测时域P=10，仿真时间为500s，利用如下程序[9]： clear all g11=poly2t