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

xueshoudaoke

最近在做电流预测控制算法，在本论坛找到了这个模型，想研究一下，有些问题需要大家给予解释一下，在此谢谢大家了！！！！！我把预测控制器的代码给贴在下面，谁能给我解释一下？
' |" L% w/ e( ~$ w5 q6 K+ f" e9 B+ e) m6 W%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0 G) F$ W! q- }. H# w- S: N% This function contains the algorithm for the predictive current control
# ]: a& [9 {  ?- f% of a two-level voltage source inverter.
3 U( s4 K4 g4 X; Y# V, t9 }% i% Inputs:
4 R3 ~5 l* H; @- n3 Q( `%   I_ref   := Two-element vector containing the reference current in
1 z# c3 ]4 n* n%              alpha-beta coordinates.
, f( r# v" J) x! B; H# s& R( i%   I_meas  := Two-element vector containing the measured current in
%              alpha-beta coordinates.
. O: v. T4 ?& _% Parameters (defined in the file parameters.m):
; q% C% Y2 L- ~5 _+ I' k%   R       := Load resistance
%   L       := Load inductance
) e1 y8 E7 X, e6 U# n; L; C8 R%   Ts      := Sampling time
%   v       := Eight-element vector containing the voltage vectors that can
( a/ K% \' ~. t. ?& A%              be generated by the inverter, in alpha-beta coordinates.
% Z+ t/ T$ w4 p" P' K9 j& X7 J' V%   states  := Eight-by-three array containing the switching states for
%              each voltage vector.
( r, L4 I9 ?7 \+ f2 d1 B  D% Outputs:
%   [Sa, Sb, Sc] := switching state corresponding to the optimum vector to
5 Z4 H1 p& }) |8 x# d1 T%                   be applied in the next sampling period.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

) o7 r) P  z6 t! J" x- gfunction [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
( t& u5 I' b5 w) G7 N% Initialize values
if isempty(x_old), x_old = 1; end
5 _! s1 J+ B& p+ qif isempty(i_old), i_old = 0+1j*0; end
$ @% v9 B5 u% r- s' Xg = zeros(1,8);
% Read current reference inputs at sampling instant k
ik_ref = I_ref(1) + 1j*I_ref(2);
2 e* z9 t: G8 G2 V; M% Read current measurements at sampling instant k
ik = I_meas(1) + 1j*I_meas(2);
0 ~  m/ S; J; Z. c# V% Back-EMF estimate
( N/ o$ H/ u  S" Ie = v(x_old) - L/Ts*ik - (R - L/Ts)*i_old;
- |. n/ L# [( i- B1 }4 y7 t: F( k# a% Store the measured current for the next iteration
3 Y: n( F  N$ mi_old = ik;
for i = 1:8
    % i-th voltage vector for current prediction
    v_o1 = v(i);
    % Current prediction at instant k+1
7 {; P# g7 O" G( O1 e3 a: ?8 ?. d9 N    ik1 = (1 - R*Ts/L)*ik + Ts/L*(v_o1 - e);
. Q/ i0 S+ S9 ]    % Cost function
0 x- z- z8 X# s3 H. O1 C    g(i) = abs(real(ik_ref - ik1)) + abs(imag(ik_ref - ik1));
* K" ^! @! o" O9 hend
/ x. k1 K: r/ p8 w1 n" l6 @% Optimization
4 k" S' b% f' H! \' A: d[~,x_opt] = min(g);
% Store the present value of x_opt
: G; T" N  c; p0 t9 hx_old = x_opt;
9 B; \6 M4 Z/ `/ ]: l% Output switching states
Sa = states(x_opt,1);
Sb = states(x_opt,2);
+ o8 |/ j8 z8 w' z3 ~; ]Sc = states(x_opt,3);