MATLAB模型预测控制工具箱函数

1、11MATLAB 模型预测控制工具箱函数8.2 系统模型建立与转换函数前面读者论坛了利用系统输入/输出数据进行系统模型辨识的有关函数及使用方 法，为时行模型预测控制器的设计，需要对系统模型进行进一步的处理和转换。MATLAB 的模型预测控制工具箱中提供了一系列函数完成多种模型转换和复杂系统 模型的建立功能。在模型预测控制工具箱中使用了两种专用的系统模型格式，即 MPC 状态空间模型和 MPC 传递函数模型。这两种模型格式分别是状态空间模型和传递函数模型在模 型预测控制工具箱中的特殊表达形式。这种模型格式化可以同时支持连续和离散系统模型的表达，在 MPC 传递函数模型中还增加了对纯时延的支持。表

2、 8-2 列出了模型预测控制工具箱的模型建立与转换函数。表 8-2 模型建立与转换函数函数名将通用状态空间模型转换为 MPC 状态空间模型 将 MPC 状态空间模型转换为通用状态空间模型 将通用传递函数模型转换为 MPC 传递函数模型 将 MPC 传递函数模型转换为 MPC 状态空间模型 将 MPC 状态空间模型转换为 MPC 阶跃响应模型 将 MPC 传递函数模型转换为 MPC 附跃响应模型 将通用状态空间模型转换为 MPC 附跃响应模型 改变 MPC 状态空间模型的采样周期 将 Theta 格式模型转换为 MPC 状态空间模型将两个开环 MPC 模型连接构成闭环模型，使其中一个模型输出叠加

3、 到另一个模型输入向 MPC 对象添加一个或多个测量扰动 向 MPC 对象添加一个或多未测量扰动 将两个 MPC 系统模型并型 将两个 MPC 系统模型串联用两个 MPC 系统模型构成增广系统模型8.2.1 模型转换在 MATLAB 模型预测工具箱中支持多种系统模型格式。 这些模型格式包括:1通用状态空间模型；2通用传递函数模型；3MPC 阶跃响应模型；4MPC 状态空间模型；5MPC 传递函数模型。在上述 5 种模型格式中，前两种模型格式是 MATLAB 通用的模型格式，在其他 控制类工具箱中，如控制系统工具箱、鲁棒控制工具等都予以支持；而后三种模型 格式化则是模型预测控制工具箱特有的。其中

4、，MPC 状态空间模型和 MPC 传递函ss2mod()mod2ss()addmod()addmd()addumd()paramod()sermod()appmod()功能11数模型是通用的状态空间模型和传递函数模型在模型预测控制工具箱中采用的增广 格式。模型预测控制工具箱提供了若干函数，用于完成上述模型格式间的转换功能。 下面对这些函数的用法加以介绍。1 通用状态空间模型与 MPC 状态空间模型之间的转换MPC 状态空间模型在通用状态空间模型的基础上增加了对系统输入/输出扰动和采样周期的描述信息，函数 ss2moc()和 mod2ss()用于实现这两种模型格式之间的 转换。1)通用状态空间模

5、型转换为 MPC 状态空间模型函数 ss2mod()该函数的调用格式为pmod= ss2moc(A,B,C,D)pmod=ss2moc(A,B,C,D,mi nfo)pmod=ss2moc(A,B,C,D,mi nfo,xO,uO,yO,fO)式中，代 B, C, D 为通用状态空间矩阵；minfo 为构成 MPC 状态空间模型的其他描述信息，为 7 个元素的向量，各元 素分别定义为：minfo(1)=dt，系统采样周期，默认值为 1;minfo(2)=n，系统阶次，默认值为系统矩阵 A 的阶次； minfo(3)=nu，受控输入的个数，默认值为系统输入的维数；minfo(4)=nd，测量扰的

6、数目，默认值为0;minfo(5)=nw，未测量扰动的数目，默认值为0;minfo(6)=nym，测量输出的数目，默认值系统输出的维数；minfo(7)=nyu，未测量输出的数目，默认值为0;注：如果在输入参数中没有指定 minfo，则取默认值。x0, u0, y0, f0 为线性化条件，默认值均为 0；pmod 为系统的 MPC 状态空间模型格式。例 8-5 将如下以传递函数表示的系统模型转换为MPC 状态空间模型。解：MATLAB 命令如下：n um=1 3 1;:其输出结果为:AB)m(D=tf2ss( numde n);pmod=ss2mod(A，B,C,D)0 10NaN -2-2-

7、11000 1 0 0 0 0 0001 0 00001310002) MPC 状态空间模型转换为通用状态空间模型函数mod2ss()该函数的调用格式为11A, B, C, D= mod2ss(pmod)A, B, C, D, mi nfo= mod2ss(pmod)A,B,C,D,mi nfo,xO,uO,yO,fO= mod2ss(pmod)式中，pmod 为系统的 MPC 状态空间模型格式；A,B,C,D 为通用状态空间矩阵; minfo为构成 MPC 状态空间模型的其他描述信息，其说明参见函数ss2moc()。2 通用传递函数模型转换为 MPC 传递函数模型通用传递函数模型与 MPC

8、传递函数模型的转换函数 poly2tfd()的调用格式为g=poly2tfd(nu m,de n, delt,delay)式中，num 为通用传递函数模型的分子多项式系数向量；den 为通用传递函数模型的分母多项式系数向量；delt 为采样周期，对连续系统，该参数为0;delay 为系统纯时延，对于离散系统，纯时延为采样时间周期的整数倍；g 为被控对象的 MPC 传递函数模型。例 8-6 考虑如下的纯时延二阶对象，并将其转换为MPC 传递函数模型。解：MATLAB 命令如下： nu m=1 1;de n=1 4 4; g=poly2tfd( nu m,de n,0,0.5)结果显示：g =0

9、1.0000 1.00001.00004.00004.000000.500003. MPC 传递函数模型转换为 MPC 状态空间模型函数 tfd2mod() 该函数的调用格式为pmod=tfd2mod(delt ,n y,g1,g2,.,g25)式中，delt 为采样时间；ny 为输出个数；g1,g2,为 SISO 传递函数，对应多变量系统传递函数矩阵的各个元素按行 向量顺序排序构成的向量，其最大个数限制为25 ；pmod 为系统的 MPC 犬态空间模型。4MPC 阶跃响应模型与其他模型格式之间的转换函数 mod2step()、 tfd2step()和 ss2step()分别用于将 MPC 状

10、态空间模型、 MPC 传递函数模型和通用状态空间模型转换为 MPC 阶跃响应模型。下面对这个三函数的 用法进行说明。111)MPC 状态空间模型转换为 MPC 阶跃响应模型函数 mod2step()plant=mod2step(pmod,tfinal)plant, dplant=mod2step(pmod,tfinal,delt,nout)式中，pmod 为系统的 MPC 状态空间模型；tfinal 为阶跃响应模型的截断时间；delt 为采样周期，默认值由 MPC 状态空间模型的参数 minfo(1)决定； nout 为输出稳定性向量，用于指定输出的稳定性。对于稳定的系统，nout 等于输出的

11、个数；对于具有一个或多个积怨分输出的系统， nout 为一个长度等于输 出个数的向量，该向量对应积分输出的分量为 0，其余分量为 1；plant 为对象在受控变量作用下的阶跃响应系数矩阵； dplant 为对象在扰动作用下阶跃响应矩阵。2)MPC 传递函数模型转换为 MPC 阶跃响应模型函数 tfd2step()该函数的调用格式为plant=tfd2step(tfinal,delt,nout,g1)pla nt二tfd2step(tfi nal,delt ,n out,g1,.,g25)式中， tfinal 为阶跃响应的截断时间；delt 为采样周期；nout 为输出稳定性向量，参见函数mod

12、2step()的有关说明；g1,g25 为 SISO 传递函数，对应多变量系统传递函数矩阵的各个元素按 行向量顺序排列构成的向量，其最大个数限制为25；plant 为对象的阶跃响应系统矩阵。例 8-8 设系统传递函数为 将其转换为阶跃响应模型。解： MATLAB 命令如下： num=1 2;den=1 3 1; tf1=poly2tfd(num,den,0,0); plant=tfd2step(5,0.1,1,tf1); plotstep(plant) 由阶跃响应模型绘制的系统阶跃响应曲线如图 8-5 所示。图 8-5 系统阶跃响应曲线8.4 基于阶跃响应模型的控制器设计与仿真函数11基于系统

13、的阶跃响应模型进行模型预测控制器设计的方法称为动态矩阵控制方 法。该方法是采用工程上易于获取的对象阶跃响应模型，算法较为简单，计算量较 少，鲁棒性较强，适用于纯时迟、开环渐近稳定的非最小相位系统，在工业过程控 制中得到成功应用。MATLAB 的模型预测控制工具箱提供了对动态矩阵控制方法的支持，有关的函 数能够基于阶跃响应模型的模型预测控制器设计与仿真，见表8-3。表 8-3 动态矩阵控制设计与仿真函数函数名功能cmpc() mpccon() mpccl ()mpcsim() nlcmpc()nlmpcsim()输入/输出有约束的模型预测控制器设计 输入/输出无约束的模型预测控制器设计 计算模型

14、预测控制系统的闭环模型模型预测控制系统的仿真(输入/输出无约束)Simulink 块 nIcmpc 对应的 S 函数Simulink 块 nImpcsim 对应的 S 函数例 8-12 考虑如下的双输入输出纯时延对象，其传递函数矩阵为解：MATLAB 程序如下：clc,clear;clos闭环系统的输出和控制量变化曲线如图1 所示。%各传递函数模型转换为阶跃响应模型g11=poly2tfd (12.8,16图11。,闭环系统的输出和控制量变化曲线8.41基于状态空间模型的预1测控制器设计函数g21在0MATLAB9模型预测控制工具箱中， 除了提供基于阶跃响应模型的预测控制器g22=poly2t

15、fd(-19.4,14.4 1,0,3);设就功能外采样还提供了 MPC 状态空间模型的预测控制器设计功能。有关的函数参见:Sny=4。tfin al=90;model=tfd2step%进行模型预函控制器设计pla nt=moselmpc()%页测时域长!度为0np=6;m=2smpccl()ywt=;uWSmpcsim()滋置输入约 1 束和参考轨迹等控制估参数器设计83.11输入/输出有约束的状态空间模型预测控制器设计tend=30;%仿真时间为 30ulim|f1-scmpb用于进行输入/输出有约束条件下的状态空间模型预测控制器设 Hylim该函数的调用格式为MuUcmflcSimpc

16、fplmOdmOShtyWtHWtiyM 忡,tend, r, ulim, ylim, kest,乙 v, w, wu) :式中皿湍蛰为 MPC 状态空间模型格式的对象状态空间模型，用于仿真；-imod 为 MPC 状态空间模型格式的对象内部模型，用于预测控制器设计;-ywt 为二次型性能指标的输出误差加权矩阵；表 8-4 基于 MPC 状态空间模型的预测控制器设计函数,(tfinal,delt,ny,g11,g12,g21,g22);功能输入/输出有约束的状态空间模型预测控制器设计 输入/输出无约束的状态空间模型预测控制器设计 计算输入/输出无约束的模型预测闭环控制系统模型 输入有约束的模型

17、预测闭环控制系统仿真11uwt 为二次型性能指标的控制量加权矩阵；M 为控制时域长度；P 为预测时域长度；tend 为仿真的结束时间；r 为输入设定值或参考轨迹；ulim=ulow uhigh delu ，式中， ulow 为控制变量的下界， uhigh 为控制变量的 上界，delu 为控制变量的变化率约束；ylim=ylow yhigh ，式中， ylow 为输出的下界， yhigh 为输出的上界；kest 为估计器的增益矩阵；z 为测量噪声； v 为测量扰动； w 为输出未测量扰动； wu 为施加到控制输入 的未测量扰动；y 为系统响应；u 为控制变量；ym 为模型预测输出。11例 8-15 考虑如下的双输入输出纯时延对象，其传递函数矩阵为解：MATLAB 程序如下：clc,clear;clos闭环系统的输出和控制量变化曲线如图2 所示。%各传递函数模型转换为状态空间模型g1 仁 poly2tfd(12.8,16.7图0,2)；闭环系统的输出和控制量变化曲线1g12李国勇tfd智能控制及其,7MATLAB实现M.北京：电子工业岀版社,2005.g21=poly2tfd(-18.9,21 1,0,3);g22=poly2tfd(-19.4,14.4 1,0,3);delt=3;ny=2;