(文章复现)微电网两阶段鲁棒优化经济调度方法(含matlab代码)

文献中两阶段鲁棒优化的模型的构建是基于下面这个确定性优化：

但是我复现的时候发现，这个模型其实是有点小问题的，第三行约束条件写的是Ky=0，但是对应的约束条件中等式右边并不全为零，所以按照常理，这一个约束应该写成Ky=k。但不确定这个文章是做了线性变换还是笔误。我也在这里纠结了好久。

后面不知怎么的，我想着给文章通信作者发邮件询问一下，没想到真的收到回复了。刘一欣老师说这里确实是个笔误，等式右边应还有一个列向量。特别感谢刘老师可以回复我这个小菜鸡，不然我要在这上面纠结很久。

我询问的邮件：

刘一欣老师的回复：

我邮件中提到的另一个问题，是关于这个对偶问题的：

因为根据强对偶定理，对偶变量对应的约束为等式约束，所以化为对偶形式时，对偶变量不应该含有符号约束，也就是≥0这个约束要删去。

同时因为约束Ky=0这一项要改成Ky=k，所以子问题的目标函数也要加上这一项。

解决了上面这些问题，后面的建模就相对没有那么困难了。

决策变量x是一个48维的列向量：

决策变量y是一个240维的列向量： 

 c和y一样都是一个240维的列向量：

 D是一个144行240列的矩阵，d是一个144维的列向量，Dy≥0可以写成矩阵形式：

表示24阶的单位矩阵，表示下三角和对角线全为1的矩阵，上三角全为0的矩阵。

K是一个50行240列的矩阵，k是一个50维的列向量，Ky=k可以写成矩阵形式：

式中，表示元素全为1的24维行向量。

F是一个96行48列的矩阵，G是一个96行240列的矩阵，h为96维的列向量，Fx+Gy≥h可以写成矩阵形式：

这样就实现了完成了所有系数矩阵的输入。

其他没有在矩阵中涉及的约束条件，用y>0表示即可。