👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆
《《《《《《《《更多资源还请持续关注本专栏》》》》》》》
附带知网可参考论文,部分展示WORD注释:
部分代码展示:
clc;
clear;
close all;
clear all class;
profile on
addpath('./example');
addpath('./functions');
addpath('./HeatGasPowerCombination');
%%
% 读入
% casename = input('Please enter case name : ', 's');
% casename = 'case14mod_SCUC_parallel';
casename = 'HeatGasPowerSystem';
% casename = 'IEEE118_new';
% 安全系数,用于留一定的裕度,针对潮流安全约束
k_safe = 0.95;
% 初始化文件
initial;
%%
%导纳矩阵计算
% [Ybus, Yf, Yt] = makeYbus(baseMVA, bus, M_branch); % build admitance matrix
[Bbus, Bf, Pbusinj, Pfinj] = makeBdc(baseMVA, bus, branch); %直流潮流
%%
% 创建决策变量
%%
% 电网
% 火电发电机出力
gen_P = sdpvar(n_gen, n_T);
gen_P_upper = sdpvar(n_gen, n_T);
% 火电机组状态
u_state = binvar(n_gen, n_T);
% 电力系统各支路功率
PF_D = sdpvar(n_branch, n_T);
% 电力系统各节点相角
Va = sdpvar(n_bus,n_T);
%%
% 气网
GasFlow = sdpvar(n_GasBranch, n_T); %各管道气流量
GasPressure2 = sdpvar(n_GasBus, n_T); %各节点气压平方
GasSourceOutput = sdpvar(n_GasSource, n_T); %各天然气源节点出力
GasGenNeed = sdpvar(n_GasGen, n_T); %各天然气发电机耗气
%%
% 热网
TmprtrFromDir = sdpvar(n_HeatBranch, n_T); %正方向支路头结点温度
TmprtrToDir = sdpvar(n_HeatBranch, n_T); %正方向支路尾结点温度
TmprtrFromRev = sdpvar(n_HeatBranch, n_T); %逆方向支路头结点温度
TmprtrToRev = sdpvar(n_HeatBranch, n_T); %逆方向支路尾结点温度
TmprtrBusDir = sdpvar(n_HeatBus,n_T); %正方向系统各节点热水的温度
TmprtrBusRev = sdpvar(n_HeatBus,n_T); %逆方向系统各节点热水的温度
HeatSource = sdpvar(n_HeatBus, n_T); %热源供热,因为电炉和CHP连在同一个节点才写的这么诡异
HeatCHP = sdpvar(n_CHPgen,n_T); %chp机组热出力
HeatEBoiler = sdpvar(n_EBoiler,n_T); %电锅炉热出力
PowerEBoiler = sdpvar(n_EBoiler,n_T); %电锅炉耗电
C = []; %约束
% C = sdpvar(C)>=0;
SCUC_value = 0;
效果展示:
👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆
《《《《《《《《更多资源还请持续关注本专栏》》》》》》》