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

xueshoudaoke

最近在做电流预测控制算法，在本论坛找到了这个模型，想研究一下，有些问题需要大家给予解释一下，在此谢谢大家了！！！！！我把预测控制器的代码给贴在下面，谁能给我解释一下？
0 j7 G: e: M4 Z! K%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
/ }% z" J0 c1 ~3 b% This function contains the algorithm for the predictive current control
+ q6 U# ?0 D5 g* e$ i( U  @% of a two-level voltage source inverter.
+ \' E7 o) C, n- e# p% Inputs:
%   I_ref   := Two-element vector containing the reference current in
%              alpha-beta coordinates.
0 R$ w) }- ?+ z%   I_meas  := Two-element vector containing the measured current in
! C# B( l. z2 r, i2 e%              alpha-beta coordinates.
% Parameters (defined in the file parameters.m):
%   R       := Load resistance
%   L       := Load inductance
%   Ts      := Sampling time
# T& o) E/ y; x3 \4 a1 a' l0 k%   v       := Eight-element vector containing the voltage vectors that can
3 @+ v/ ?: \' O" [6 a2 T" ~%              be generated by the inverter, in alpha-beta coordinates.
%   states  := Eight-by-three array containing the switching states for
%              each voltage vector.
" Y  k/ N. I4 c8 ~% Outputs:
%   [Sa, Sb, Sc] := switching state corresponding to the optimum vector to
%                   be applied in the next sampling period.
' Z* G; q. p( H0 f# K%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
9 w& L, ~* Y; w
function [Sa,Sb,Sc] = control(I_ref,I_meas,R,L,Ts,v,states)
$ K; I% k" i' P) r0 N% Optimum vector and measured current at instant k-1
persistent x_old i_old
5 W: K7 F8 C1 D7 w% Initialize values
if isempty(x_old), x_old = 1; end
( F- u; \5 [, o# [if isempty(i_old), i_old = 0+1j*0; end
+ C% }* N* G# Og = zeros(1,8);
0 r. Y9 w& `" J: Y6 a+ B( z% Read current reference inputs at sampling instant k
ik_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
; ]1 N, |- B; J# V% p1 je = v(x_old) - L/Ts*ik - (R - L/Ts)*i_old;
% Store the measured current for the next iteration
i_old = ik;
for i = 1:8
# C- c" {9 e5 c% z/ _- h+ j" ^7 [. K    % i-th voltage vector for current prediction
, L; ?  ?) [- X1 g$ o    v_o1 = v(i);
  u+ Y$ z. @- m) R+ w6 l: ~    % Current prediction at instant k+1
    ik1 = (1 - R*Ts/L)*ik + Ts/L*(v_o1 - e);
1 O* x; H- V: |8 A6 L" |    % Cost function
8 J& R8 A+ W8 x    g(i) = abs(real(ik_ref - ik1)) + abs(imag(ik_ref - ik1));
. ]+ _$ H6 Y- g5 Iend
, P7 u. p) K0 z' v8 J  K% Optimization
[~,x_opt] = min(g);
% Store the present value of x_opt
x_old = x_opt;
1 ^) b% q% S; m7 t% Output switching states
Sa = states(x_opt,1);
$ x5 G* L6 O: c9 O9 hSb = states(x_opt,2);
Sc = states(x_opt,3);