在matlab中通过yalmip平台调用cplex求解器，可用于求解MILP问题，适合于综合能源系统优化求解（完整程序分享）

综合能源系统优化求解

完整程序：

%% 请先确保YALMIP工具箱和CPLEX正确安装,MATLAB导入对应文件，否则无法运行程序！！
%CPLEX免费试用版对求解规模有限制，如出现规模过大无法求解，请购买正式版或申请教育版！！！
%本程序已在yalmip version：20171121 CPLEX version：：12.8教育版验证成功

%% 初始化
clc;
clear;
yalmip;
Cplex;
%% 各变量及常量定义
%------------------------变量定义-----------------------%
Pw=sdpvar(1,24,'full'); %风机出力
Ppv=sdpvar(1,24,'full');%光伏出力

Pbat=sdpvar(1,24,'full');%蓄电池出力
% Pde=Power(4,:);%柴油机组出力
% Lshift=zeros(1,24);
Pnet=sdpvar(1,24,'full');%交换功率
Pbuy=sdpvar(1,24,'full');%从电网购电电量
Psell=sdpvar(1,24,'full');%向电网售电电量
Temp_net=binvar(1,24,'full'); % 购|售电标志
Temp_cha=binvar(1,24,'full'); %充电标志
Temp_dis=binvar(1,24,'full'); %放电标志
Temp_static=binvar(1,24,'full'); %电池静置标志
Pcha=sdpvar(1,24);
Pdis=sdpvar(1,24);
Constraints = [];
%-------------------------常量定义-----------------------%
Load=[88.24    83.01    80.15    79.01    76.07    78.39    89.95    128.85    155.45    176.35    193.71    182.57    179.64    166.31    164.61    164.61    174.48    203.93    218.99    238.11    216.14    173.87    131.07    94.04];
%风机预测出力
Pw=[66.9   68.2   71.9   72   78.8   94.8   114.3   145.1   155.5   142.1   115.9   127.1   141.8   145.6...
145.3   150   206.9   225.5   236.1   210.8   198.6   177.9   147.2   58.7];
%光伏预测出力
Ppv=[0   0   0   0   0.06   6.54   20.19   39.61   49.64   88.62   101.59   66.78   110.46   67.41   31.53...
50.76   20.6   22.08   2.07   0   0   0   0   0];
%分时电价
C_buy=[0.25   0.25   0.25   0.25   0.25   0.25   0.25   0.53   0.53   0.53   0.82   0.82...
0.82   0.82   0.82   0.53   0.53   0.53   0.82   0.82   0.82   0.53   0.53   0.53];
C_sell=[0.22   0.22   0.22   0.22   0.22   0.22   0.22   0.42   0.42   0.42   0.65   0.65...
0.65   0.65   0.65   0.42   0.42   0.42   0.65   0.65   0.65   0.42   0.42   0.42];

%% 约束
for k = 1:24
Constraints = [Constraints, -160<=Pnet(1,k)<=160,0<=Pbuy(1,k)<=160, -160<=Psell(1,k)<=0]; %主网功率交换约束
Constraints =[Constraints,Pnet(1,k)+Pw(1,k)+Ppv(1,k)==Load(1,k)+Pbat(1,k)]; %功率平衡约束
Constraints = [Constraints, implies(Temp_net(1,k),[Pnet(1,k)>=0,Pbuy(1,k)==Pnet(1,k),Psell(1,k)==0])]; %购电情况约束
Constraints = [Constraints, implies(1-Temp_net(1,k),[Pnet(1,k)<=0,Psell(1,k)==Pnet(1,k),Pbuy(1,k)==0])]; %售电情况约束
%----------------------蓄电池约束--------------------%
sum_bat=zeros(1,24);
Constraints = [Constraints, -40<=Pbat(1,k)<=40,0<=Pcha(1,k)<=40,-40<=Pdis(1,k)<=0];%电池充放电约束
Constraints = [Constraints, implies(Temp_cha(1,k),[Pbat(1,k)>=0,Pcha(1,k)==Pbat(1,k),Pdis(1,k)==0])];%充电情况约束
Constraints = [Constraints, implies(Temp_dis(1,k),[Pbat(1,k)<=0,Pdis(1,k)==Pbat(1,k),Pcha(1,k)==0])];%放电情况约束
Constraints = [Constraints, implies(Temp_static(1,k),[Pbat(1,k)==0,Pdis(1,k)==0,Pcha(1,k)==0])];%静置情况约束
Constraints = [Constraints,Temp_cha(1,k)+Temp_dis(1,k)+Temp_static(1,k)==1];
% sum_bat(1,k+1)=sum_bat(1,k)+Pcha(1,k)+Pdis(1,k);%计算SOC
% Constraints=[Constraints,-30<=sum(Pdis(1,1:k)+Pcha(1,1:k))<=165] ;%SOC约束，电池容量300kwh，初始S0C为0.4，0.3<=SOC<=0.95
end
Constraints=[Constraints,sum(Pdis+Pcha)==0] ;%ST=S0，始末SOC相等约束
%% 目标函数
F=0;
%------------------总费用--------------------%
for k = 1:24
F = F + 0.52*Pw(:,k)+0.72*Ppv(:,k)+C_buy(:,k)*Pbuy(:,k)+C_sell(:,k)*Psell(1,k)+0.2*Pdis(1,k);
end
ops=sdpsettings('solver', 'cplex');%参数指定程序用cplex求解器
optimize(Constraints,F,ops)
value(F)%费用
value(Pcha)
value(Pdis)
%% 画图
x=1:24;
PP=[Pbuy;-Pdis;Pw;Ppv];
PP_neg=[Psell;-Pcha];
figure
bar(PP','stack');
h=legend('交换功率','蓄电池','风机出力','光伏出力','Location','NorthWest');
set(h,'Orientation','horizon')
hold on
bar(PP_neg','stack');
plot(x,value(Load),'r','linewidth',2);
xlabel('时段');ylabel('功率/kW');
hold off
for k=1:24
s(k)=value(sum(Pdis(1,1:k)+Pcha(1,1:k)))/300+0.4;
soc(k+1)=s(k);
end
soc(1)=0.4;
xx=0:24;
figure
plot(xx,soc,'r');
xlabel('时段');ylabel('SOC值');
title('蓄电池SOC状态');

在matlab中通过yalmip平台调用cplex求解器，可用于求解MILP问题，适合于综合能源系统优化求解（完整程序分享）相关推荐

MATLAB中安装YALMIP及CPLEX详细步骤
一.MATLAB中安装YALMIP (YALMIP压缩包下载链接:链接:百度网盘请输入提取码密码:ujfp) 1.下载YALMIP-master工具压缩包,YALMIP-master工具压缩包解压 ...
路径调度问题（CVRP）后续之如何在MATLAB中安装YALMIP及CPLEX包
订阅博主专栏的用户可私信博主获取安装包,不用额外再付费下载安装包. 安装包下载地址:https://download.csdn.net/download/wenyusuran/14954633 一.M ...
MATLAB中安装Yalmip及Cplex
注意点!! Yalmip及Cplex安装对版本匹配要求很高推荐使用与自己安装的MATLAB相匹配的安装版本 Yalmip安装各版本Yalmip下载地址 1. 下载所需版本后,将解压后的文件夹保存至 ...
matlab 中的 yalmip安装方法
matlab 中的 yalmip安装方法在 https://yalmip.github.io/download/ 网站上下载 yalmip ,根据自己的需求下载想要的版本 . 将下载好的文件解压后, ...
【微电网优化】基于matlab粒子群算法求解综合能源系统优化问题【含Matlab源码 1969期】
一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[微电网优化]基于matlab粒子群算法求解综合能源系统优化问题[含Matlab源码 1969期] 点击上面蓝色字体,直接付费下载,即可. ...
【智能优化求解】基于粒子群算法实现综合能源系统优化附matlab代码
1 简介为了解决现有冷热电联供型综合能源系统大多只单一考虑系统机组投资成本或系统环境污染,影响系统整体优化运行的问题,以系统经济性和环保性为目标,对冷热电联供系统进行研究分析.构建含燃气轮机.燃气锅 ...
MATLAB代码：电-气-热综合能源系统耦合优化调度代码主要做的是一个考虑电网、热网以及气网耦合调度的综合能源系统优化调度模型
MATLAB代码:电-气-热综合能源系统耦合优化调度关键词:综合能源系统优化调度电气热耦合仿真平台:MATLAB YALMIP+cplex/gurobi 主要内容:代码主要做的是一个考虑电网. ...
MATLAB程序：综合能源系统优化调度，考虑了阶梯型碳交易机制和氢能，具有一定的创新。
MATLAB程序:综合能源系统优化调度,考虑了阶梯型碳交易机制和氢能,具有一定的创新. 采用CPLEX+Yalmip求解,基本复现. ID:6990667995938854
MATLAB代码：考虑用户舒适度的冷热电多能互补综合能源系统优化调度
MATLAB代码:考虑用户舒适度的冷热电多能互补综合能源系统优化调度关键词:用户舒适度综合能源 PMV 优化调度仿真平台:MATLAB+yalmip+cplex 主要内容:代码主要做的是考虑用户 ...

在matlab中通过yalmip平台调用cplex求解器，可用于求解MILP问题，适合于综合能源系统优化求解（完整程序分享）

在matlab中通过yalmip平台调用cplex求解器，可用于求解MILP问题，适合于综合能源系统优化求解（完整程序分享）相关推荐

最新文章

热门文章