本程序代码参考EI期刊论文《基于多智能体系统一致性算法的电力系统分布式经济调度策略》,主要做的是电力系统分布式调度策略,一致性变量为发电机组的增量成本和柔性负荷的增量效益,程序中算例丰富,注释清晰,干货满满,下面对程序和文章作简要介绍!
集中式优化方法难以应对未来电网柔性负荷广泛渗透以及电力元件“即插即用”的技术要求。 区别于集中式经济调度,提出一种电力系统分布式经济调度策略。 应用多智能体系统中的一致性算法,以发电机组的增量成本和柔性负荷的增量效益作为一致性变量,设计一种用于电力系统经济调度的算法,通过分布式优化的方式求解经济调度问题。 算例仿真与分析验证了所提调度策略的有效性。
运行平台:MATLAB!
算例系统:10机19节点系统!
文章结果:
程序结果:
部分程序:
%发电机%
be=[13.065 5.295 11.37 3.36 12.795 11.775 3.375 9.435 6.45 12.39];
ga=[0.0046 0.0111 0.0099 0.0095 0.0104 0.0029 0.0021 0.0062 0.0077 0.0048];
pg=[135.88 214.92 108.04 127.69 232.56 240.00 44.628 234.48 74.600 172.09];
spg=zeros(1,301);
spg(1)=sum(pg);
pgmax=2.5*pg;
pgmin=-2.5*pg;
li=be+2.*ga.*pg;
%负荷%
b=[25.755 18.42 27.63 10.59 16.275 28.365 28.14 23.55 21.42 15.225 28.56 10.305 23.94 22.05 26.25 16.455 24.375 26.295 14.76];
c=[14.02 6.25 15.1 8.42 8.1 21.21 11.9 15.96 12.75 6.9 9.75 8.23 9.28 9.49 9.1 34.08 18.36 12.39 13.04]./(-100);
pd=[110.15 176.75 109.69 75.55 120.64 80.26 142.02 88.57 100.8 132.38 175.75 75.13 154.69 139.3 172.85 28.98 79.67 127.37 67.92];
spd=zeros(1,301);
spd(1)=sum(pd);
pdmax=2.5*pd;
pdmin=-2.5*pd;
lj=b+2.*c.*pd;
%迭代%
l=[li lj];
ll=zeros(301,29);
ll(1,:)=l;
dp=zeros(1,301);
dp(1)=spd(1)-spg(1);
pgg=zeros(301,10);
pgg(1,:)=pg;
pdd=zeros(301,19);
pdd(1,:)=pd;
for t=1:1:300
for n=1:1:29 %更新lambda的值
if n==1
ll(t+1,n)=sum(a(n,:).*ll(t,:))+0.005*dp(t); % 这里和矩阵相乘是一样的 a(n,:)*ll(t,:)'
elseif n==11
ll(t+1,n)=sum(a(n,:).*ll(t,:))+0.005*dp(t);
else
ll(t+1,n)=sum(a(n,:).*ll(t,:));
end
end
for i=1:1:10 %判断pg是否越限并赋值
if (ll(t,i)-be(i))/(2*ga(i))>=pgmax(i)
pgg(t+1,i)=pgmax(i);
elseif (ll(t,i)-be(i))/(2*ga(i))<=pgmin(i)
pgg(t+1,i)=pgmin(i);
else
pgg(t+1,i)=(ll(t,i)-be(i))/(2*ga(i));
end
end
for j=1:1:19 %判断pd是否越限并赋值
if (ll(t,j+10)-b(j))/(2*c(j))>=pdmax(j)
pdd(t+1,j)=pdmax(j);
elseif (ll(t,j+10)-b(j))/(2*c(j))<=pdmin(j)
pdd(t+1,j)=pdmin(j);
else
pdd(t+1,j)=(ll(t,j+10)-b(j))/(2*c(j));
end
end
spg(t+1)=sum(pdd(t+1,:));
spd(t+1)=sum(pgg(t+1,:));
dp(t+1)=sum(pdd(t+1,:))-sum(pgg(t+1,:));
end
%绘图%
以上就是本次介绍的主要内容,欢迎感兴趣的小伙伴关注下方公众号获取完整版代码,小编会继续推送更有质量的学习资料、文章和程序代码,为您的科研加油助力!