科技公司网站建设方案书模板,娄底网站建设的公司,wordpress dedecms整合,做网站设计参考文献目录
0 概述
1 配电网重构的目标函数
2 算例
3 matlab代码实现 0 概述
配电系统中存在大量的分段开关和联络开关#xff0c;配电网重构正是通过调整分段开关和联络升大的组合状态来变换网络结构,用于优化配电网某些指标#xff0c;使其达到最优状态。正常运行时,则通…
目录
0 概述
1 配电网重构的目标函数
2 算例
3 matlab代码实现 0 概述
配电系统中存在大量的分段开关和联络开关配电网重构正是通过调整分段开关和联络升大的组合状态来变换网络结构,用于优化配电网某些指标使其达到最优状态。正常运行时,则通过两类开关的不同组合状态来提高系统的可靠性、稳定性或减小网损等。故障时﹐部分分段开关将打开以隔离故障支路,同时一些联络开关闭合将故障支路转移到其他馈线上或同一馈线的其余分支上以满足供电需求。
1 配电网重构的目标函数
配电网结构是闭环设计开环运行。配电网重构则是在保证网络拓扑结构成辐射状满足电压约
束、支路容量约束、变压器容量约束等要求的前提下确定使配电网某一指标或多个指标最佳的配电网运行方式。本文以配电网线路损耗最小为目标函数即: 式中f为所有支路总有功损耗;i表示支路编号;N表示配电网支路总数;Ki表示开关状态变量,即О表示打开状态,1表示闭合状态;Ri表示支路i上的电阻值;Pi、Qi、Ui分别表示支路的有功功率、无功功率以及支路i前端节点的电压幅值。 配电网重构约束条件包括潮流约束、电压约束、电流约束、线路容量约束和网络拓扑结构辐射状运行约束等如下式所示: 2 算例
仿真算例采用 IEEE33 节点配电网络拓扑结构如图 。该配电网包含 33 个节点和 37 条支路
5 个常开联络开关分别位于支路33、34、35、36、37上基准 电压为 12 66 kV。 3 matlab代码实现
部分代码
部分代码
r2rand(n,dim);% 创建随机矩阵大小( 20x3 )
iter0;
maxiter60;% 最大迭代次数
tap[8 9 10 11 21 33 35 0 02 3 4 5 6 7 18 19 2012 13 14 34 0 0 0 0 015 16 17 29 30 31 36 32 022 23 24 25 26 27 28 37 0];
tatap;
%% 建立关联矩阵
dataloadcase(case33);
docdata.branch;
nhanh37;
nut33;
matrixzeros(nhanh,nut);nutdaudoc(:,1);nutcuoidoc(:,2);
for i1:nhanh
matrix(i,nutdau(i))1;
matrix(i,nutcuoi(i))1;
end%% 计算个体最优的适应度函数fpbestzeros(1,n);for i1:nfpbest(i)50000;end
%% 主循环
while itermaxiteriteriter1;wwmax-(wmax-wmin)*iter/maxiter;% 指定权重系数c12*rand(1);c22*rand(1);%% 更新速度voldvnew;for i1:nfor j1:dimvnew(i,j)w*vnew(i,j)c1*r1(i,j)*(pbest(i,j)-x(i,j))c2*r2(i,j)*(gbest(j)-x(i,j));if abs(vnew(i,j))abs(vold(i,j))vnew(i,j)rand(1,1).*vnew(i,j);endendend-------
plot(dienap,-sr)
hold on
plot(gbestvolt,-^b)
ylabel(电压(p.u))
xlabel(Node)
title(电压分布)
legend(配电网重组前,配电网重组后)
hold off
disp( )
disp( )
disp( )
disp( )
disp()
disp(***************************** 33节点配电网仿真结果 ***************************)
disp()
disp( 重构前 重构后 )
disp(------------------------------------------------------------------------------------------)
disp([联络开关: , num2str(o), ,num2str(a)])
disp(------------------------------------------------------------------------------------------)
disp([功率损耗: ,num2str(tonthat), kW, ,num2str(fgbest), kW])
disp(------------------------------------------------------------------------------------------)
disp([减少的功率损耗: , _______, ,num2str(ploss), %])
disp(------------------------------------------------------------------------------------------)
disp([最小电压: ,num2str(dienapmin), pu, ,num2str(minvolt), pu])
disp(------------------------------------------------------------------------------------------)
toc 部分代码 iter0; maxiter60;% 最大迭代次数 tap[8 9 10 11 21 33 35 0 0 2 3 4 5 6 7 18 19 20 12 13 14 34 0 0 0 0 0 15 16 17 29 30 31 36 32 0 22 23 24 25 26 27 28 37 0]; tatap; %% 建立关联矩阵 dataloadcase(case33); docdata.branch; nhanh37; nut33; matrixzeros(nhanh,nut); nutdaudoc(:,1); nutcuoidoc(:,2); for i1:nhanh matrix(i,nutdau(i))1; matrix(i,nutcuoi(i))1; end
%% 计算个体最优的适应度函数 fpbestzeros(1,n); for i1:n fpbest(i)50000; end %% 主循环 while itermaxiter iteriter1; wwmax-(wmax-wmin)*iter/maxiter;% 指定权重系数 c12*rand(1); c22*rand(1); %% 更新速度 voldvnew; for i1:n for j1:dim vnew(i,j)w*vnew(i,j)c1*r1(i,j)*(pbest(i,j)-x(i,j))c2*r2(i,j)*(gbest(j)-x(i,j)); if abs(vnew(i,j))abs(vold(i,j)) vnew(i,j)rand(1,1).*vnew(i,j); end end end ------- plot(dienap,-sr) hold on plot(gbestvolt,-^b) ylabel(电压(p.u)) xlabel(Node) title(电压分布) legend(配电网重组前,配电网重组后) hold off disp( ) disp( ) disp( ) disp( ) disp() disp(***************************** 33节点配电网仿真结果 ***************************) disp() disp( 重构前 重构后 ) disp(------------------------------------------------------------------------------------------) disp([联络开关: , num2str(o), ,num2str(a)]) disp(------------------------------------------------------------------------------------------) disp([功率损耗: ,num2str(tonthat), kW, ,num2str(fgbest), kW]) disp(------------------------------------------------------------------------------------------) disp([减少的功率损耗: , _______, ,num2str(ploss), %]) disp(------------------------------------------------------------------------------------------) disp([最小电压: ,num2str(dienapmin), pu, ,num2str(minvolt), pu]) disp(------------------------------------------------------------------------------------------) toc
