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

xueshoudaoke

最近在做电流预测控制算法，在本论坛找到了这个模型，想研究一下，有些问题需要大家给予解释一下，在此谢谢大家了！！！！！我把预测控制器的代码给贴在下面，谁能给我解释一下？
) R# {0 k4 v' e2 D" ]%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1 C& y' I3 {$ F% This function contains the algorithm for the predictive current control
% of a two-level voltage source inverter.
% Inputs:
%   I_ref   := Two-element vector containing the reference current in
# m% ?5 @, j+ B%              alpha-beta coordinates.
%   I_meas  := Two-element vector containing the measured current in
8 F2 Y1 G+ W! v# Q+ h%              alpha-beta coordinates.
$ ?+ I' u! C0 [6 L7 ^# m! W% Parameters (defined in the file parameters.m):
%   R       := Load resistance
4 P6 A- Q' o/ a# E# J7 k- I) ^%   L       := Load inductance
; S/ Y& a6 Y5 I9 X8 R% p: |# r%   Ts      := Sampling time
. [+ S4 @+ }  X' F2 `4 o/ l%   v       := Eight-element vector containing the voltage vectors that can
5 V9 o* Z! ^# {- w- ^%              be generated by the inverter, in alpha-beta coordinates.
%   states  := Eight-by-three array containing the switching states for
/ @: Q& F. L1 D# \: W1 b$ q% x2 |%              each voltage vector.
; C  S/ \+ O; H' Y8 N% Outputs:
%   [Sa, Sb, Sc] := switching state corresponding to the optimum vector to
: z8 i% @) w4 F- J% X0 m%                   be applied in the next sampling period.
: B/ d- t7 A( ~* C  v( Q%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [Sa,Sb,Sc] = control(I_ref,I_meas,R,L,Ts,v,states)
% Optimum vector and measured current at instant k-1
2 {  k- x3 h% W3 L* b7 ?# Cpersistent x_old i_old
% Initialize values
if isempty(x_old), x_old = 1; end
* s, c5 r# F; uif isempty(i_old), i_old = 0+1j*0; end
g = zeros(1,8);
% Read current reference inputs at sampling instant k
ik_ref = I_ref(1) + 1j*I_ref(2);
' G( S: B& r4 P. G, L- D$ k% Read current measurements at sampling instant k
ik = I_meas(1) + 1j*I_meas(2);
. n/ t1 G4 L0 |1 q' D# N% Back-EMF estimate
e = v(x_old) - L/Ts*ik - (R - L/Ts)*i_old;
8 k6 z* h& @1 ]' f2 R% Store the measured current for the next iteration
i_old = ik;
for i = 1:8
2 {" K( B3 h+ D: ^% H- L, K    % i-th voltage vector for current prediction
7 f0 x* k% }: W& s& V    v_o1 = v(i);
    % Current prediction at instant k+1
    ik1 = (1 - R*Ts/L)*ik + Ts/L*(v_o1 - e);
& ~% r  b/ _+ K    % Cost function
: O; y& n; q+ @$ V; o' A    g(i) = abs(real(ik_ref - ik1)) + abs(imag(ik_ref - ik1));
end
% Optimization
[~,x_opt] = min(g);
' ]& N* b& Y, c0 X( O5 P9 K% Store the present value of x_opt
x_old = x_opt;
% Output switching states
Sa = states(x_opt,1);
  L( b9 v) Q# u  o0 ~# T8 oSb = states(x_opt,2);
Sc = states(x_opt,3);