预测控制课程报告.

预测控制

课程论文

题目DMC算法matlab编程与仿真学生姓名

学号

学院

专业

指导教师

二0一三年一月五日

DMC算法matlab编程与仿真

摘要: 预测控制在控制方面有重要作用，而动态矩阵控制是预测控制中一种重要的算法。本文分析了动态矩阵控制算法的原理以及算法包括的三个主要方面：预测模型、滚动优化、反馈控制。并通过仿真实例，来进一步表明动态矩阵控制算法的优越性。

关键词：预测控制，DMC，仿真

一、引言

随着科技的发展，人们对大型复杂和不确定性系统的控制品质要求逐渐提高，因此就需要提出一种新的计算机控制算法。利用状态空间法分析和设计系统，不仅提高了人们对控制对象的洞察能力，而且还提供了在更高层次上设计控制系统的手段。因为工业对象的结构、参数和环境都有很大的不确定性，按照理想模型得到的最优控制在实际中不能保持最优，有时会引起控制品质严重下降。预测控制的提出不要求对模型结构有先验知识，也不需要通过复杂的系统辨识，直接就可以设计控制系统。预测控制算法汲取了现代控制理论中的优化思想，滚动的在线优化，克服了不确定性，增强控制系统的鲁棒性。

预测控制算法一般分为三个部分，分别是预测模型、滚动优化和反馈校正。预测模型是展示系统未来动态行为的功能，任意给出未来的控制策略观察对象的在不同控制策略下的输出变化，为比较这些控制策略的优劣提供依据。滚动优化不是一次离线进行的，而是反复在线进行的。反馈校正可以在保持模型不变的基础上对未来的误差进行预测并加以校正，还可以在线辨识直接修改预测模型。因此预测控制能有效地应用于复杂系统，它在工业过程和其它领域有着诱人的应用前景。

二、动态矩阵控制算法( DMC)

预测控制是智能控制方法之一，目前提出的预测控制算法主要有基于非参数模型的模型算法控制( MAC) 、动态矩阵控制( DMC) 和基于参数模型的广义预测控制( GPC) 等。动态矩阵控制算法是一种基于被控对象非参数数学模型的控制算法，它是一种基于对象阶跃响应的预测控制算法，它以对象的阶跃响应离散系数为模型，避免了通常的传递函数或状态空间方程模型参数的辨识，又因为采用多步预估技术，能有效解决时延过程问题，并按预估输出与给定值偏差最小的二次性能指标实施控制，它适用于渐进稳定的线性对象，系统的动态特性中具有纯滞后或非最小相位特性都不影响改算法的直接应用，因此是一种最优控制技术［4］。20 世纪80 年代初，Culter 提出的动态矩阵预测控制( DMC) 算法是MPC 中的一种重要的算法，是一种利用被控对象的单位阶跃响应采样序列作为预测模型的预测控制算法。算法原理图1 所示。

图1 DMC 原理图

DMC 算法分为三部分:预测模型、滚动优化和反馈校正。

2.1预测模型

在DMC 算法中，首先需要测定对象单位阶跃响应的采样值)(iT a a i =,i=1,2,……。其中，T 为采样周期。对于渐进稳定的对象，阶跃响应在某一NT N =t ，后将趋于平稳，以至i a (i>N)与N a 的误差和量化误差及测量误差有相同的数量级，因而可认为，N a 已近似等于阶跃响应的稳态值)(∞=a a S 。这样，对象的动态信息就可以近似用有限集合{}N a a a ⋯⋯21,加以描述。这个集合的参数构成了DMC 的模型参数，向量a={}N a a a ⋯⋯21,称为模型向量，N 称为模型时域长度。

虽然阶跃响应是一种非参数模型，但由于线性系统具有比例和叠加性质，故利用这组模型参数{}i a ，已足以预测在任意输入作用下系统在未来时刻的输出值。在t=kT 时刻，假如控制量不再变化时系统在未来N 个时刻的输出值为

)(~),2(~),1(~000k N k y k k y k k y +⋯⋯++，,那么，在控制增量)(k u ∆作用后系统的输出

可由

)()(~)(~01k u a k y k y N N ∆+= （2-1）

预测，其中⎥⎥⎥

⎦⎤

⎢⎢⎢⎣⎡++=)(~)1(~)(~000k N k y k k y k y N M 表示在t=kT 时刻预测的尚无)(k u ∆作用时未来N 个时刻

的系统输出。⎥⎥⎥

⎦

⎤⎢⎢⎢⎣⎡++=)(~)1(~)(~111k N k y k k y k y N M 表示在t=kT 时刻预测的有控制增量)(k u ∆作用时未来

N 个时刻的系统输出。⎥⎥⎥

⎦

⎤⎢⎢⎢⎣⎡=N a a a M 1为阶跃响应模型向量，其元素为描述系统动态特性的N 个阶

跃响应系数。式中，上标~表示预测，k i k +表示在t=kT 时刻预测t=(k+i)T 时刻。 同样，如果考虑到现在和未来M 个时刻控制增量的变化，在t=kT 时刻预测在控制增量

)(k u ∆,⋯,)1(-+∆M k u 作用下系统在未来P 个时刻的输出为

)()(~)(~0k u A k y k y M P PM ∆+= （2-2）

式中⎥⎥⎥

⎦⎤

⎢⎢⎢⎣⎡++=)(~)1(~)(~000P k N k y k k y k y M 为t=kT 时刻预测的无控制增量时未来P 个时刻的系统输出。

⎥⎥⎥⎦

⎤⎢⎢⎢⎣⎡++=)P (~)1(~)(~M M PM k k y k k y k y M 为t=kT 时刻预测的有M 个控制增量)(k u ∆,……，)1(-+∆M k u 时未来P 个时刻的系统输出。⎥⎥

⎥⎦

⎤⎢⎢⎢⎣⎡-+∆∆=∆)1()()(M k u k u k u M M 为从现在起M 个时刻的控制增量。

⎥⎥⎥⎥⎥⎦

⎤

⎢⎢⎢⎢

⎢⎣⎡=+--11121000M p p p a a a a a a A ΛM ΛM M ΛΛ称为动态矩阵，其元素为描述系统动态特性的阶跃响应系数。

2.2滚动优化

DMC 是一种以优化确定控制策略的算法。在采样时刻t=kT 的优化性能指标可取为

∑=+-+=P i M i k i k y i k q k J 1

2)](~)([)(min ω+∑=-+∆M

1

2)1(j j j k u r （2-3）

即通过选择该时刻起M 个时刻的控制增量)(k u ∆,⋯，)1(-+∆M k u ，使系统在未来P （M P N ≥≥）个时刻的输出值)1(~k k y M +，…，)(~k P k y M + 尽可能接近其期望值

)1(+k ω，…)(P k +ω。性能指标中的第二项是对控制增量的约束，即不允许控制量的变

化过于剧烈。式中，i q ，j r 为权系数，P 和M 分别称为优化时域长度和控制时域长度。