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

xueshoudaoke

最近在做电流预测控制算法，在本论坛找到了这个模型，想研究一下，有些问题需要大家给予解释一下，在此谢谢大家了！！！！！我把预测控制器的代码给贴在下面，谁能给我解释一下？
9 l1 M( ]. |( _$ p6 c9 @4 S  W%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This function contains the algorithm for the predictive current control
% of a two-level voltage source inverter.
8 M) t* J' b7 P% Inputs:
; \: F! _8 j9 `7 G+ B%   I_ref   := Two-element vector containing the reference current in
- E0 W  y+ S  P3 `  ]. ?%              alpha-beta coordinates.
%   I_meas  := Two-element vector containing the measured current in
/ J( @4 L6 t/ e* D5 ^* i5 _7 u%              alpha-beta coordinates.
! M5 ?- y+ C  i( \" b; b9 |% Parameters (defined in the file parameters.m):
%   R       := Load resistance
%   L       := Load inductance
%   Ts      := Sampling time
0 y! b# C! ?* a2 l1 v7 M0 A+ A8 M%   v       := Eight-element vector containing the voltage vectors that can
4 ]: G' P% G7 Y, p  y' g%              be generated by the inverter, in alpha-beta coordinates.
! J; r1 H/ p6 N% \%   states  := Eight-by-three array containing the switching states for
9 m$ T- T8 }  X7 O" h8 g$ D%              each voltage vector.
4 \/ r  p( l# `3 M/ w7 v. R% Outputs:
! C, c% {1 x5 R+ A2 o' y5 S+ M%   [Sa, Sb, Sc] := switching state corresponding to the optimum vector to
%                   be applied in the next sampling period.
9 w6 m- k+ _, m# n; U' Z1 S4 Y%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Q3 b+ {* E. W' ?
function [Sa,Sb,Sc] = control(I_ref,I_meas,R,L,Ts,v,states)
+ k5 W3 _5 z+ a& z% Optimum vector and measured current at instant k-1
- i% H! ~  k0 Kpersistent x_old i_old
0 N- D; g" D3 O0 r5 Z5 `; H- b% Initialize values
if isempty(x_old), x_old = 1; end
8 e9 I% M/ ?- R4 n  xif isempty(i_old), i_old = 0+1j*0; end
g = zeros(1,8);
% Read current reference inputs at sampling instant k
$ i, y- Q: V* Q5 l, S3 n. bik_ref = I_ref(1) + 1j*I_ref(2);
% Read current measurements at sampling instant k
+ b& L; B* l. b* {* Y" Lik = I_meas(1) + 1j*I_meas(2);
% Back-EMF estimate
3 M2 P! Z; ?8 Je = v(x_old) - L/Ts*ik - (R - L/Ts)*i_old;
" M: l  O4 r2 C1 K% Store the measured current for the next iteration
i_old = ik;
. q! {7 @/ O0 s( p' c# qfor i = 1:8
    % i-th voltage vector for current prediction
    v_o1 = v(i);
    % Current prediction at instant k+1
    ik1 = (1 - R*Ts/L)*ik + Ts/L*(v_o1 - e);
, t2 ~# |/ o7 k% x; R' B    % Cost function
    g(i) = abs(real(ik_ref - ik1)) + abs(imag(ik_ref - ik1));
end
% Optimization
3 o# r( l/ }4 E( ~" T' y[~,x_opt] = min(g);
% Store the present value of x_opt
, ]& m; n# l. R/ C* |x_old = x_opt;
% Output switching states
Sa = states(x_opt,1);
Sb = states(x_opt,2);
Sc = states(x_opt,3);