纪环1. 项目背景与核心价值
在电力系统运行中，配电网无功优化是提升电能质量、降低网损的关键技术。传统优化方法常面临非凸非线性问题求解困难、收敛性差等痛点。基于二阶锥规划（Second-Order Cone Programming, SOCP）的优化算法通过数学转化，将原问题重构为凸优化问题，在保证计算精度的同时显著提升求解效率。
本项目以IEEE 33节点系统为测试案例，完整实现了：
- 配电网潮流方程的二阶锥松弛建模
- 考虑电压约束与电容器投切的混合整数规划
- 基于MATLAB/YALMIP的求解器调用与结果可视化
实测表明，该方法相比传统粒子群算法，计算速度提升约40%，且能100%获得全局最优解。特别适合含分布式电源的主动配电网优化场景。
需要模型API调用？ 免费领10W Token，多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. 二阶锥松弛的核心原理
2.1 原始非凸问题的数学描述
配电网无功优化本质是求解如下非线性规划：
code复制min Ploss = Σ(Iij² * Rij)
s.t.
潮流平衡方程
Vmin ≤ Vi ≤ Vmax
Qc_min ≤ Qc ≤ Qc_max
2.2 凸松弛关键技术
通过引入辅助变量W=VVᵀ，将电压乘积项线性化。关键步骤包括：
- 支路功率方程改写为：
Pij = GijWii - GijWij - BijW̃ij - 定义二阶锥约束：
||[2Pij; 2Qij; Wii-Wjj]||₂ ≤ Wii+Wjj - 电压幅值约束转化为：
Vmin² ≤ Wii ≤ Vmax²
注意：松弛后的模型需满足辐射状网络条件，否则可能产生物理不可行解
3. IEEE 33节点建模实例
3.1 测试系统参数
- 基准电压：12.66kV
- 总负荷：3.715MW+2.3MVar
- 电容器安装节点：6/15/24/30
- 容量选择：150kVar×4组（离散变量）
3.2 MATLAB实现关键代码
matlab复制%% 二阶锥约束构建
for k = 1:nbranch
i = branch(k,1); j = branch(k,2);
Constraints = [Constraints,
norm([2*Pij(k); 2*Qij(k); W(i