最近在进行电力系统最优潮流计算的过程中,想利用yalmip在MATLAB中编程实现,CSDN上有利用GAMS的优化方法,但是鉴于以后在求解电力系统优化问题时,yalmip使用的场景更加广泛,所以想自己编程利用yalmip调用gurobi求解最优潮流。(当然调用其他求解器进行求解也是可以的。)网上搜索很久都没有找到yalmip进行最优潮流计算的程序,只能自己想着写代码。
使用的算例是CSDN一位博主的,就是前面说的利用GAMS进行优化的那位。链接https://blog.csdn.net/qq_34194377/article/details/105303603?utm_medium=distribute.pc_relevant.none-task-blog-title-2&spm=1001.2101.3001.4242。
下面是我用MATLAB进行的编程代码:
clc;
clear all;
Gen_max = [40 170 520 200 600];%发电机最大出力
Gen_min = [0 0 0 0 0];%发电机最小出力
Branch_max = [400 400 400 400 400 240];
Branch_min = -Branch_max;
Q = [14 15 30 40 20];%成本系数
bus_num = 5;%节点个数
slackbus = 1;%平衡节点
branch_num = 6;%支路个数
Gen_num = 5;%发电机台数
GenBus = [1