基于MATLAB的YALMIP+cplex的配电网两阶段可调鲁棒优化

提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档

---

前言

	本文以系统网损成本为目标，求解配电网两阶段鲁棒优化问题。在规划求解过程中，利用二阶锥松弛法、大M法将非线性非凸问题转化为易于求解的混合整数二阶锥规划问题。同时将模型分解为min-max-min的主子问题，利用强对偶理论将max-min 凸锥问题转化为单层max问题，再使用C&CG算法交替迭代求解。最后，通过在改进的IEEE33节点系统上得到的仿真结果，验证本文模型的有效性。

一、学习流程

1.写出基本的原问题pri，只有IEEE33的交流潮流约束，为二阶锥形式；
11.写出基本的对偶问题dual，只有IEEE33的交流潮流约束，为对偶锥形式，并且其值应与pri保持一致或误差极小；
2.在pri的基础上写出包含连续变量SVC的原问题pri1；
22.写出包含SVC的对偶问题dual1;
3.在pri1的基础上添加离散变量CB，为原问题pri2;
33.由于此时引入了离散变量，离散变量是由主问题传递而来，因此CB中的0-1变量由主问题给定，即pri2给定，子问题仅需在目标函数中作出变化即可；
4.在pri2的基础上，添加时间耦合变量ESS，由于涉及到前后时刻变化（即SOC约束），单纯的求偏导行不通，需要引入新变量代替下一时刻的SOC，去除SOC（t+1）=SOC(t)+充放电功率变化量，这条约束。而添加新的等式约束。由于这是配网鲁棒，无法用矩阵对偶来处理此时间耦合，所以显得较为复杂。
44.同样的，由于ESS也涉及到了离散变量，因此充放电状态由pri3给定，而连续变量部分仍需在dual3中体现。
5.在pri3的基础上，添加不确定性变量，成为pri4；
55.不确定性变量为盒式集，在保留原源荷约束的对偶约束的基础上，添加新变量，并引入可调系数，这部分在很多论文都有。
6.到此所有内容均已单独实现，接下来就是将其整合，引入CCG，MP由pri4变化而来，SP由dual4变化而来。

总结

提示：持续更新。。。