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

xueshoudaoke

最近在做电流预测控制算法，在本论坛找到了这个模型，想研究一下，有些问题需要大家给予解释一下，在此谢谢大家了！！！！！我把预测控制器的代码给贴在下面，谁能给我解释一下？
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
7 c# d. F! x6 K8 w3 C* l) p9 P6 k% This function contains the algorithm for the predictive current control
- U6 u6 \& v& `: I/ ]5 X: |) S% of a two-level voltage source inverter.
- Z/ W' T  p% L9 i% 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):
%   R       := Load resistance
% K2 ?# J2 U1 f% S8 J%   L       := Load inductance
6 y6 T7 a, r7 e; r; F%   Ts      := Sampling time
, b' L6 \) s* E  Z" A+ I%   v       := Eight-element vector containing the voltage vectors that can
- J% k' B' K# Z: l+ E- ?$ ^5 K%              be generated by the inverter, in alpha-beta coordinates.
& m4 |, z0 g5 k) O% `( ~( l" `%   states  := Eight-by-three array containing the switching states for
+ r& `1 t1 e8 [6 T& v%              each voltage vector.
- k( _1 r( n' e- K" W- h; Y5 z3 E% M% Outputs:
%   [Sa, Sb, Sc] := switching state corresponding to the optimum vector to
%                   be applied in the next sampling period.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
( g( ]0 I" c+ n  Q/ M
/ _& Z$ [; T" T0 E1 f# Rfunction [Sa,Sb,Sc] = control(I_ref,I_meas,R,L,Ts,v,states)
% Optimum vector and measured current at instant k-1
persistent x_old i_old
% Initialize values
if isempty(x_old), x_old = 1; end
2 V9 H; H: E+ q* \if isempty(i_old), i_old = 0+1j*0; end
g = zeros(1,8);
' R! m# s) w3 t! f8 {% Read current reference inputs at sampling instant k
2 }6 J1 K  T1 A) f/ Bik_ref = I_ref(1) + 1j*I_ref(2);
$ _; X  _# t8 f$ X0 ]# K% Read current measurements at sampling instant k
ik = I_meas(1) + 1j*I_meas(2);
% Back-EMF estimate
3 q" Z, v. e( U4 S  c& ]e = v(x_old) - L/Ts*ik - (R - L/Ts)*i_old;
0 K6 F7 F6 v5 @& M$ [% a6 S# ?9 s% Store the measured current for the next iteration
i_old = ik;
for i = 1:8
    % i-th voltage vector for current prediction
6 C( C* |# E$ `) i    v_o1 = v(i);
    % Current prediction at instant k+1
    ik1 = (1 - R*Ts/L)*ik + Ts/L*(v_o1 - e);
# M2 V5 B% c& u7 U! B7 l5 e; \    % Cost function
    g(i) = abs(real(ik_ref - ik1)) + abs(imag(ik_ref - ik1));
end
% Optimization
[~,x_opt] = min(g);
: W8 V& O6 H: h& @/ P. Y7 z' s% Store the present value of x_opt
& H- j: z) I8 ]/ _/ H- _x_old = x_opt;
% Output switching states
Sa = states(x_opt,1);
Sb = states(x_opt,2);
# S3 t; t8 y2 A7 l5 JSc = states(x_opt,3);