电流预测控制simulink仿真，帮忙解释一下！！！！

xueshoudaoke

最近在做电流预测控制算法，在本论坛找到了这个模型，想研究一下，有些问题需要大家给予解释一下，在此谢谢大家了！！！！！我把预测控制器的代码给贴在下面，谁能给我解释一下？
$ G3 E1 P  f2 G6 Q! T  H%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This function contains the algorithm for the predictive current control
% of a two-level voltage source inverter.
8 a9 [3 a$ f( n" M# `/ W$ C; v% Inputs:
%   I_ref   := Two-element vector containing the reference current in
%              alpha-beta coordinates.
%   I_meas  := Two-element vector containing the measured current in
%              alpha-beta coordinates.
% Parameters (defined in the file parameters.m):
0 B: F7 F$ O7 Z%   R       := Load resistance
%   L       := Load inductance
%   Ts      := Sampling time
& Q& B* C7 l6 V: `: b: g6 L  }%   v       := Eight-element vector containing the voltage vectors that can
%              be generated by the inverter, in alpha-beta coordinates.
%   states  := Eight-by-three array containing the switching states for
%              each voltage vector.
0 M& J; z. W4 W+ _. V% Outputs:
%   [Sa, Sb, Sc] := switching state corresponding to the optimum vector to
! J* q! |1 z) K* o! R+ |%                   be applied in the next sampling period.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [Sa,Sb,Sc] = control(I_ref,I_meas,R,L,Ts,v,states)
% Optimum vector and measured current at instant k-1
( w* M% x9 v3 }persistent x_old i_old
# E; m- L; f  G4 h9 v3 K% Initialize values
if isempty(x_old), x_old = 1; end
if isempty(i_old), i_old = 0+1j*0; end
g = zeros(1,8);
3 f8 a8 @1 \, j- c( E- _% Read current reference inputs at sampling instant k
' B$ z( @9 Y4 p' mik_ref = I_ref(1) + 1j*I_ref(2);
% Read current measurements at sampling instant k
ik = I_meas(1) + 1j*I_meas(2);
% Back-EMF estimate
3 L; X- v- U0 L# me = v(x_old) - L/Ts*ik - (R - L/Ts)*i_old;
' M+ G1 {) E3 L; c6 ^: I% Store the measured current for the next iteration
i_old = ik;
# d: U* d; c$ Z: R" g0 H0 p8 ffor i = 1:8
    % i-th voltage vector for current prediction
    v_o1 = v(i);
% e# o+ ^2 H$ x% r    % Current prediction at instant k+1
$ w+ P5 f% l, b2 S) w' ^    ik1 = (1 - R*Ts/L)*ik + Ts/L*(v_o1 - e);
    % Cost function
    g(i) = abs(real(ik_ref - ik1)) + abs(imag(ik_ref - ik1));
end
+ l7 G# X2 [% b, \. L# X) v1 ^& t! N- N% Optimization
[~,x_opt] = min(g);
( J( F5 j3 L3 e& _. c8 O/ j% Store the present value of x_opt
x_old = x_opt;
% Output switching states
, e3 {$ o& @' k  D8 @6 ^' Y1 XSa = states(x_opt,1);
Sb = states(x_opt,2);
( k" m! w  \* z$ V4 @Sc = states(x_opt,3);