三点估计法概率潮流计算
标签:
•
文件类型: .zip
•
文件大小: 10.3MB
•
下载次数: 1
•
1.matlab代码,在matpower环境下进行潮流计算,用略微修改的ieee30节点,程序在main.m中
2.仅考虑了负荷的波动
如有原理性的错误望指出
代码片段和文件信息
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
目录 0 2017-09-30 21:36 matpower6.0
目录 0 2017-05-29 18:22 matpower6.0@opf_model
文件 711 2016-12-17 01:34 matpower6.0@opf_modeldisplay.m
文件 434 2016-12-17 01:34 matpower6.0@opf_modelget_mpc.m
文件 2226 2016-12-17 01:34 matpower6.0@opf_modelopf_model.m
目录 0 2017-05-29 18:22 matpower6.0@opt_model
文件 8151 2016-12-17 01:34 matpower6.0@opt_modeladd_constraints.m
文件 9467 2016-12-17 01:34 matpower6.0@opt_modeladd_costs.m
文件 6060 2016-12-17 01:34 matpower6.0@opt_modeladd_vars.m
文件 6993 2016-12-17 01:34 matpower6.0@opt_modeluild_cost_params.m
文件 3145 2016-12-17 01:34 matpower6.0@opt_modelcompute_cost.m
文件 2293 2016-12-17 01:34 matpower6.0@opt_modeldescribe_idx.m
文件 5177 2016-12-17 01:34 matpower6.0@opt_modeldisplay.m
文件 645 2016-12-17 01:34 matpower6.0@opt_modelget.m
文件 1625 2016-12-17 01:34 matpower6.0@opt_modelgetN.m
文件 2196 2016-12-17 01:34 matpower6.0@opt_modelget_cost_params.m
文件 2168 2016-12-17 01:34 matpower6.0@opt_modelget_idx.m
文件 4465 2016-12-17 01:34 matpower6.0@opt_modelgetv.m
文件 4191 2016-12-17 01:34 matpower6.0@opt_modellinear_constraints.m
文件 7901 2016-12-17 01:34 matpower6.0@opt_modelopt_model.m
文件 1132 2016-12-17 01:34 matpower6.0@opt_modeluserdata.m
文件 4559 2017-05-29 18:38 matpower6.0ABC.m
文件 319 2016-12-17 01:34 matpower6.0AUTHORS
文件 972 2016-12-17 01:34 matpower6.0Contents.m
文件 2840 2016-12-17 01:34 matpower6.0LICENSE
文件 19884 2016-12-17 01:34 matpower6.0README
文件 20337 2016-12-17 01:34 matpower6.0README.txt
文件 5191 2016-12-17 01:34 matpower6.0add_userfcn.m
文件 17385 2016-12-17 01:34 matpower6.0apply_changes.m
文件 2065 2016-12-17 01:34 matpower6.0ustypes.m
文件 33539 2016-12-17 01:33 matpower6.0case118.m
............此处省略496个文件信息
function ABC
%%算法参数设置
NP=12; %%NP=雇佣蜂数=跟随蜂数
limit=24;
maxCycle=100;
%控制变量上下限
ub=[1.05 8 5]; %机端电压、变压器档位、并联电容器档位
lb=[0.95 0 0]; %变压器变比=档位*0.025+0.9 并联电容器容量=档位*10(MW)
%系统参数
shuju=data_ieee30;
n1=shuju.gen(:1); %发电机节点
n2=find(shuju.branch(:9)~=0); %变压器支路
n3=find(shuju.bus(:6)~=0); %并联电容器节点
d1=length(n1);
d2=length(n2);
d3=length(n3);
dim=[d1 d2 d3]; %控制变量维数
ND=sum(dim); %总维数
%初始化
for i=1:NP
foods(i:)=foodGenerate(ublbdim); %蜜源,即解
shuju=shujuGenerate(shujufoods(i:)dimn1n2n3);
objVal(i)=objFun(shuju); %目标函数值
fitness(i)=1/(objVal(i)+1); %适应度
end
[best bestInd]=min(objVal);
globalMin=best;
globalSolution=foods(bestInd:);
globalMins=zeros(1maxCycle);
trial=zeros(1NP); %解未更新代数
for iter=1:maxCycle
%雇佣蜂工作 vid=xid+r(xid-xjd)
for i=1:NP
d=randi(ND);
j=randi(NP);
while j==i
j=randi(NP);
end
%得到新的蜜源
sol=foods(i:);
sol(d)=foods(id)+(foods(id)-foods(jd))*(rand-0.5)*2;
%限制蜜源位置
sol=bound(solublbdim);
%计算适应度
shuju=shujuGenerate(shujusoldimn1n2n3);
objValSol=objFun(shuju);
fitnessSol=1/(1+objValSol);
%贪婪选择
if fitnessSol>fitness(i)
foods(i:)=sol;
fitness(i)=fitnessSol;
objVal(i)=objValSol;
trial(i)=0;
else
trial(i)=trial(i)+1;
end
end
prob=fitness/sum(fitness);
%跟随蜂开始工作
i=1;
t=0;
while t if rand t=t+1;
d=randi(ND);
j=randi(NP);
while j==i
j=randi(NP);
end
%得到新的蜜源
sol=foods(i:);
sol(d)=foods(id)+(foods(id)-foods(jd))*(rand-0.5)*2;
%限制蜜源位置
sol=bound(solublbdim);
%计算适应度
shuju=shujuGenerate(shujusoldimn1n2n3);
objValSol=objFun(shuju);
fitnessSol=1/(1+objValSol);
%贪婪选择
if fitnessSol>fitness(i)
foods(i:)=sol;
fitness(i)=fitnessSol;
objVal(i)=objValSol;
trial(i)=0;
else
trial(i)=trial(i)+1;
end
end
i=i+1;
if i==NP+1
i=1;
end
end
[best bestInd]=min(objVal);
if best globalMin=best;
globalSolution=foods(bestInd:);
end
globalMins(iter)=globalMin;
%侦查蜂开始工作
[maxTrial ind]=max(trial);
if maxTrial>limit
foods(ind:)=foodGenerate(ublbdim);
shuju=shujuGenerate(shujufoods(ind:)dimn1n2n3);
objVal(ind)=objFun(shuju);
fitness(ind)=1/(1+objVal(ind));
trial(ind)=0;
end
end
globalMins
plot(globalMins)
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
目录 0 2017-09-30 21:36 matpower6.0
目录 0 2017-05-29 18:22 matpower6.0@opf_model
文件 711 2016-12-17 01:34 matpower6.0@opf_modeldisplay.m
文件 434 2016-12-17 01:34 matpower6.0@opf_modelget_mpc.m
文件 2226 2016-12-17 01:34 matpower6.0@opf_modelopf_model.m
目录 0 2017-05-29 18:22 matpower6.0@opt_model
文件 8151 2016-12-17 01:34 matpower6.0@opt_modeladd_constraints.m
文件 9467 2016-12-17 01:34 matpower6.0@opt_modeladd_costs.m
文件 6060 2016-12-17 01:34 matpower6.0@opt_modeladd_vars.m
文件 6993 2016-12-17 01:34 matpower6.0@opt_modeluild_cost_params.m
文件 3145 2016-12-17 01:34 matpower6.0@opt_modelcompute_cost.m
文件 2293 2016-12-17 01:34 matpower6.0@opt_modeldescribe_idx.m
文件 5177 2016-12-17 01:34 matpower6.0@opt_modeldisplay.m
文件 645 2016-12-17 01:34 matpower6.0@opt_modelget.m
文件 1625 2016-12-17 01:34 matpower6.0@opt_modelgetN.m
文件 2196 2016-12-17 01:34 matpower6.0@opt_modelget_cost_params.m
文件 2168 2016-12-17 01:34 matpower6.0@opt_modelget_idx.m
文件 4465 2016-12-17 01:34 matpower6.0@opt_modelgetv.m
文件 4191 2016-12-17 01:34 matpower6.0@opt_modellinear_constraints.m
文件 7901 2016-12-17 01:34 matpower6.0@opt_modelopt_model.m
文件 1132 2016-12-17 01:34 matpower6.0@opt_modeluserdata.m
文件 4559 2017-05-29 18:38 matpower6.0ABC.m
文件 319 2016-12-17 01:34 matpower6.0AUTHORS
文件 972 2016-12-17 01:34 matpower6.0Contents.m
文件 2840 2016-12-17 01:34 matpower6.0LICENSE
文件 19884 2016-12-17 01:34 matpower6.0README
文件 20337 2016-12-17 01:34 matpower6.0README.txt
文件 5191 2016-12-17 01:34 matpower6.0add_userfcn.m
文件 17385 2016-12-17 01:34 matpower6.0apply_changes.m
文件 2065 2016-12-17 01:34 matpower6.0ustypes.m
文件 33539 2016-12-17 01:33 matpower6.0case118.m
............此处省略496个文件信息
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件举报,一经查实,本站将立刻删除。
评论列表(条)