上次写了一篇【Simulink】基于FCS-MPC的三相并网逆变器控制（Matlab Function），控制部分用 Matlab Function 写的，因实验室一般用 Stateflow，所以这篇把 Matlab Function 改成 Stateflow.

原理

电路原理图：

直流电压源与三相桥式逆变器相连，经过RL滤波器并网。

三相桥式逆变器包含6个IGBT，由于同相上下臂的IGBT不能同时导通，因此总共有 $2^3=8$ 种状态，且有两个状态计算得到的代价函数相同。

FCS-MPC 步骤：
步骤1：根据系统结构建立合适的系统模型。
步骤2：采集系统三相电流、电网电压值。
步骤3：根据系统模型建立离散化预测函数，预测变换器输出电流或并网功率。
步骤4：建立合适的代价函数，将8种开关状态依次带入代价函数。
步骤5：通过枚举的方式选择使代价函数最小的矢量所对应的开关状态进行下一控制周期的控制。

状态0：
$cost0=abs(i\alpha (k+1)-i\alpha ref)+abs(i\beta (k+1)-i\beta ref)=abs((1-R\ast Ts/L)\ast i\alpha (k)+Ts/L\ast (0-e\alpha (k))-ir\alpha )+abs((1-R\ast Ts/L)\ast i\beta (k)+Ts/L\ast (0-e\beta (k))-ir\beta ):$

状态1：
$cost1=abs(i\alpha (k+1)-i\alpha ref)+abs(i\beta (k+1)-i\beta ref)=abs((1-R\ast Ts/L)\ast i\alpha (k)+Ts/L\ast (-Udc/3-e\alpha (k))-ir\alpha )+abs((1-R\ast Ts/L)\ast i\beta (k)+Ts/L\ast (-Udc/sqrt(3)-e\beta (k))-ir\beta ):$

状态2：
$cost2=abs(i\alpha (k+1)-i\alpha ref)+abs(i\beta (k+1)-i\beta ref)=abs((1-R\ast Ts/L)\ast i\alpha (k)+Ts/L\ast (-Udc/3-e\alpha (k))-ir\alpha )+abs((1-R\ast Ts/L)\ast i\beta (k)+Ts/L\ast (Udc/sqrt(3)-e\beta (k))-ir\beta ):$

状态3：
$cost3=abs(i\alpha (k+1)-i\alpha ref)+abs(i\beta (k+1)-i\beta ref)=abs((1-R\ast Ts/L)\ast i\alpha (k)+Ts/L\ast (-2\ast Udc/3-e\alpha (k))-ir\alpha )+abs((1-R\ast Ts/L)\ast i\beta (k)+Ts/L\ast (0-e\beta (k))-ir\beta ):$

状态4：
$cost4=abs(i\alpha (k+1)-i\alpha ref)+abs(i\beta (k+1)-i\beta ref)=abs((1-R\ast Ts/L)\ast i\alpha (k)+Ts/L\ast (2\ast Udc/3-e\alpha (k))-ir\alpha )+abs((1-R\ast Ts/L)\ast i\beta (k)+Ts/L\ast (0-e\beta (k))-ir\beta ):$

状态5：
$cost5=abs(i\alpha (k+1)-i\alpha ref)+abs(i\beta (k+1)-i\beta ref)=abs((1-R\ast Ts/L)\ast i\alpha (k)+Ts/L\ast (Udc/3-e\alpha (k))-ir\alpha )+abs((1-R\ast Ts/L)\ast i\beta (k)+Ts/L\ast (-Udc/sqrt(3)-e\beta (k))-ir\beta ):$

状态6：
$cost6=abs(i\alpha (k+1)-i\alpha ref)+abs(i\beta (k+1)-i\beta ref)=abs((1-R\ast Ts/L)\ast i\alpha (k)+Ts/L\ast (Udc/3-e\alpha (k))-ir\alpha )+abs((1-R\ast Ts/L)\ast i\beta (k)+Ts/L\ast (Udc/sqrt(3)-e\beta (k))-ir\beta ):$

主电路

250V直流电压源与三相桥式逆变器相连，经过RL滤波器并网。

控制电路

Clark变换

Stateflow

记得要在模型管理器定义变量

运行结果

当参考电流幅值为 6A 时：

当参考电流幅值为 6A→3A 时：

可以看到，和采用 Matlab Function 编写代码的效果是一样的。