Переключающиеся контроллеры на основе оптимальных затрат

В этом примере показано, как использовать "оптимальную стоимость" выходной порт блока MPC Controller, чтобы переключить между многоуровневой моделью прогнозирующие контроллеры, выходные параметры которых ограничиваются дискретными значениями.

Задайте модель объекта управления

Линейная модель объекта управления следующие::

plant = ss(tf({1,1},{[1 1.2 1],[1 1]}),'min');  % Plant with 2 inputs and 1 output
[A,B,C,D] = ssdata(plant);                      % Get state-space realization matrices, to be used in Simulink
x0 = [0;0;0];                                   % Initial plant state

Спроектируйте контроллер MPC

Задайте типы сигнала ввода и вывода.

plant = setmpcsignals(plant,'MV',1,'MD',2);     % First input is manipulated, second is measured disturbance

Спроектируйте два контроллера MPC с ограничениями мВ u=-1 и u=1, соответственно. Только u в текущее время квантуется. Последующие расчетные действия управления могут быть любым значением между-1 и 1. Диспетчер использует отступающий подход горизонта, таким образом, эти значения на самом деле не переходят к объектам.

Ts = 0.2;                                           % Sampling time
p = 20;                                             % Prediction horizon
m = 10;                                             % Control horizon
mpc1 = mpc(plant,Ts,p,m);                             % First MPC object
mpc2 = mpc(plant,Ts,p,m);                             % Second MPC object
% Specify weights
mpc1.Weights = struct('MV',0,'MVRate',.3,'Output',1);   % Weights
mpc2.Weights = struct('MV',0,'MVRate',.3,'Output',1);   % Weights
% Specify constraints
mpc1.MV = struct('Min',[-1;-1],'Max',[-1;1]);       % Constraints on the manipulated variable: u = -1
mpc2.MV = struct('Min',[1;-1],'Max',[1;1]);         % Constraints on the manipulated variable: u = 1
-->The "Weights.ManipulatedVariables" property of "mpc" object is empty. Assuming default 0.00000.
-->The "Weights.ManipulatedVariablesRate" property of "mpc" object is empty. Assuming default 0.10000.
-->The "Weights.OutputVariables" property of "mpc" object is empty. Assuming default 1.00000.
-->The "Weights.ManipulatedVariables" property of "mpc" object is empty. Assuming default 0.00000.
-->The "Weights.ManipulatedVariablesRate" property of "mpc" object is empty. Assuming default 0.10000.
-->The "Weights.OutputVariables" property of "mpc" object is empty. Assuming default 1.00000.

Симулируйте в Simulink®

Чтобы запустить этот пример, Simulink® требуется.

if ~mpcchecktoolboxinstalled('simulink')
    disp('Simulink(R) is required to run this example.')
    return
end
% Specify signals:
Tstop = 40;
ref.time = 0:Ts:(Tstop+p*Ts);
ref.signals.values = double(ref.time>10)';          % Step change in reference signal at time t=10
md.time = ref.time;
md.signals.values = double(md.time>30)';            % Step change in measured disturbance signal at time t=30

Откройте и симулируйте модель Simulink®:

mdl = 'mpc_optimalcost';
open_system(mdl);                     % Open Simulink(R) Model
sim(mdl,Tstop);                       % Start Simulation
-->Converting model to discrete time.
-->Assuming output disturbance added to measured output channel #1 is integrated white noise.
-->The "Model.Noise" property of the "mpc" object is empty. Assuming white noise on each measured output channel.
-->Converting model to discrete time.
-->Assuming output disturbance added to measured output channel #1 is integrated white noise.
-->The "Model.Noise" property of the "mpc" object is empty. Assuming white noise on each measured output channel.

Обратите внимание на то, что:

  • Со времени 0 ко времени 10, действие управления продолжает переключаться между MPC1 (-1) и MPC2 (+1). Это вызвано тем, что опорный сигнал 0, и он требует, чтобы контроллер выход в 0 достиг устойчивого состояния, которое не может быть достигнуто ни с одним контроллером MPC.

  • Со времени 10 - 30, управление MPC2 вывело (+1), выбран, потому что опорный сигнал становится +1, и это требует, чтобы контроллер выход в +1 достиг устойчивого состояния (усиление объекта равняется 1), который может быть достигнут MPC2.

  • Со времени 30 - 40, действие управления начинает переключаться снова. Это вызвано тем, что с присутствием измеренного воздействия (+1), MPC1 приводит к устойчивому состоянию 0, и MPC2 приводит к устойчивому состоянию +2, в то время как опорный сигнал все еще требует +1.

Симулируйте Используя команду MPCMOVE

Используйте mpcmove выполнить постепенную симуляцию и вычислить текущее действие MPC управления:

[Ad,Bd,Cd,Dd] = ssdata(c2d(plant,Ts));              % Discrete-time dynamics
Nsteps = round(Tstop/Ts);                           % Number of simulation steps

Инициализируйте матрицы для хранения результатов симуляции

YY = zeros(Nsteps+1,1);
RR = zeros(Nsteps+1,1);
UU = zeros(Nsteps+1,1);
COST = zeros(Nsteps+1,1);
x = x0;                                             % Initial plant state
xt1 = mpcstate(mpc1);                               % Initial state of the MPC controller #1
xt2 = mpcstate(mpc2);                               % Initial state of the MPC controller #2

Запустите симуляцию.

for td=0:Nsteps
    % Construct signals
    v = md.signals.values(td+1);
    r = ref.signals.values(td+1);
    % Plant equations: output update
    y = Cd*x + Dd(:,2)*v;
    % Compute both MPC moves
    options = mpcmoveopt;
    options.OnlyComputeCost = true; % the last input argument tells "mpcmove" to only compute the optimal cost
    [u1,Info1] = mpcmove(mpc1,xt1,y,r,v,options);
    [u2,Info2] = mpcmove(mpc2,xt2,y,r,v,options);
    % Compare the resulting optimal costs and choose the input value
    % corresponding to the smallest cost
    if Info1.Cost<=Info2.Cost
        u = u1;
        cost = Info1.Cost;
        % Update internal MPC state to the correct value
        xt2.Plant = xt1.Plant;
        xt2.Disturbance = xt1.Disturbance;
        xt2.LastMove = xt1.LastMove;
    else
        u = u2;
        cost = Info2.Cost;
        % Update internal MPC state to the correct value
        xt1.Plant = xt2.Plant;
        xt1.Disturbance = xt2.Disturbance;
        xt1.LastMove = xt2.LastMove;
    end
    % Store plant information
    YY(td+1) = y;
    RR(td+1) = r;
    UU(td+1) = u;
    COST(td+1) = cost;
    % Plant equations: state update
    x = Ad*x + Bd(:,1)*u + Bd(:,2)*v;
end

Постройте результаты mpcmove соответствовать результатам симуляции, полученным в Simulink®:

subplot(131)
plot((0:Nsteps)*Ts,[YY,RR]);                      % Plot output and reference signals
grid
title('OV and Reference')
subplot(132)
plot((0:Nsteps)*Ts,UU);                           % Plot manipulated variable
grid
title('MV')
subplot(133)
plot((0:Nsteps)*Ts,COST);                         % Plot optimal MPC value function
grid
title('Optimal cost')

Эти графики напоминают графики в осциллографах в модели Simulink®.

bdclose(mdl);

Смотрите также

Похожие темы

Для просмотра документации необходимо авторизоваться на сайте