预测控制模型代码

预测控制模型建立与实现开发技术matlab语言编程实现

function [s rho omega HPhiJ g] = GPCcoef（n p alpha beta）
% Copyright （c） 2014 Yiping Cheng ypcheng@bjtu.edu.cn
% This is to compute the GPC coefficients used by the algorithm
% For details see my paper “A novel DMC-like implementation of GPC“
% p is the prediction horizon
% alpha and beta must be of dimension 1-by-n
% the last element of alpha must be 0

s=zeros（n1）;
for i=1:n
  s（i） = sum（beta（1:i））;
end

rho = zeros（p 1）;
for j=1:p
  if j    rho（j）=-alpha（j）;
  end
  for i=1:min（n-1 j-1）
    rho（j） = rho（j） -alpha（i）*rho（j-i）;
  end
end

omega = zeros（p 1）;
for j=1:p
  omega（j） = 1+sum（rho（1:j））;
end

HPhiJ = zeros（p n）;
for k=1:n-1
  for j=2:min（p+1k-1）
    HPhiJ（j-1k） = 0;
  end
  if k>=2 && k<=p+1
    HPhiJ（k-1k） = 1;
  end
  fo

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件          82  2014-11-03 15:22  radical.m
     文件       20532  2014-11-16 17:14  TwoTank.mdl
     文件        1064  2014-11-16 17:19  readme_twotank.txt
     文件       86297  2014-11-16 16:51  TwoTank.jpg
     文件        3553  2014-11-16 13:02  T2TankControl.m
     文件        1633  2014-11-16 13:02  T2Tank.m
     文件        1232  2014-11-16 10:38  GPCcoef.m
     文件        1529  2014-11-16 04:27  license.txt