资源简介
Locational marginal cost for pjm 5-bus
代码片段和文件信息
% %%case4
GSF=[];
B=[-2050/3010/3;50/3-155/62.520/3;02.5-75/1420/7;10/320/320/7-90/7];
p=[-1 0 0];
theta123=(B(1:31:3)^-1)*p‘;
GSF(11)=(theta123(1)-theta123(2))/0.06;
GSF(21)=0;
GSF(31)=(theta123(1)-0)/0.3;
GSF(41)=(theta123(2)-0)/0.15;
GSF(51)=(theta123(3)-0)/0.35;
GSF(61)=(theta123(2)-theta123(3))/0.4;
p=[0 -1 0];
theta123=(B(1:31:3)^-1)*p‘;
GSF(12)=(theta123(1)-theta123(2))/0.06;
GSF(22)=0;
GSF(32)=(theta123(1)-0)/0.3;
GSF(42)=(theta123(2)-0)/0.15;
GSF(52)=(theta123(3)-0)/0.35;
GSF(62)=(theta123(2)-theta123(3))/0.4;
p=[0 0 -1];
theta123=(B(1:31:3)^-1)*p‘;
GSF(13)=(theta123(1)-theta123(2))/0.06;
GSF(23)=0;
GSF(33)=(theta123(1)-0)/0.3;
GSF(43)=(theta123(2)-0)/0.15;
GSF(53)=(theta123(3)-0)/0.35;
GSF(63)=(theta123(2)-theta123(3))/0.4;
f=[10 15 20 30];
lb=[0 0 0 0];
ub=[300 50 60 200];
b=[135.94];
A=[GSF(11)GSF(12)GSF(12)GSF(13)];
Aeq=[1 1 1 1];
beq=[450];
[xfvalexitflagoutputlamda]=linprog(fAbAeqbeqlbub);
MP=lamda.eqlin;
LMP=-lamda.eqlin-lamda.ineqlin*GSF(11:3);
LMP=[LMP-lamda.eqlin];
gencost1=10*x(1)+15*x(2)+20*x(3)+30*x(4);
genrev=LMP(1)*x(1)+LMP(2)*(x(2)+x(3))+LMP(3)*x(4);
%%
%case1
- 上一篇:提取语音特征短时能量,平均幅度、平均过零率的matlab代码
- 下一篇:梯度投影法
评论
共有 条评论