91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

  • 大小: 10.41MB
    文件類型: .rar
    金幣: 2
    下載: 0 次
    發(fā)布日期: 2023-11-19
  • 語(yǔ)言: Matlab
  • 標(biāo)簽: matlab??IEEE30??

資源簡(jiǎn)介

壓縮包為matlab書寫的源代碼,文件數(shù)較多,具體完成任務(wù)就是利用IEEE30節(jié)點(diǎn)設(shè)置了負(fù)荷的變化在matpower環(huán)境下進(jìn)行潮流計(jì)算,主程序?yàn)閙ain文件中

資源截圖

代碼片段和文件信息

function?ABC
%%算法參數(shù)設(shè)置
NP=12;???%%NP=雇傭蜂數(shù)=跟隨蜂數(shù)
limit=24;
maxCycle=100;

%控制變量上下限
ub=[1.05?8?5];???%機(jī)端電壓、變壓器檔位、并聯(lián)電容器檔位
lb=[0.95?0?0];????%變壓器變比=檔位*0.025+0.9??并聯(lián)電容器容量=檔位*10(MW)

%系統(tǒng)參數(shù)
shuju=data_ieee30;
n1=shuju.gen(:1);??????%發(fā)電機(jī)節(jié)點(diǎn)
n2=find(shuju.branch(:9)~=0);??%變壓器支路
n3=find(shuju.bus(:6)~=0);?????%并聯(lián)電容器節(jié)點(diǎn)
d1=length(n1);
d2=length(n2);
d3=length(n3);??????????
dim=[d1?d2?d3];????????%控制變量維數(shù)
ND=sum(dim);???????????%總維數(shù)

%初始化
for?i=1:NP
????foods(i:)=foodGenerate(ublbdim);???%蜜源,即解
????shuju=shujuGenerate(shujufoods(i:)dimn1n2n3);
????objVal(i)=objFun(shuju);??????????????%目標(biāo)函數(shù)值
????fitness(i)=1/(objVal(i)+1);???????????%適應(yīng)度
end

[best?bestInd]=min(objVal);
globalMin=best;
globalSolution=foods(bestInd:);
globalMins=zeros(1maxCycle);

trial=zeros(1NP);????%解未更新代數(shù)?????

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);
????????%計(jì)算適應(yīng)度
????????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);
????
????%跟隨蜂開(kāi)始工作
????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);
????????????%計(jì)算適應(yīng)度
????????????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;
????%偵查蜂開(kāi)始工作
????[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)

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

????.......???????711??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opf_model\display.m

????.......???????434??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opf_model\get_mpc.m

????.......??????2226??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opf_model\opf_model.m

????.......??????8151??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\add_constraints.m

????.......??????9467??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\add_costs.m

????.......??????6060??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\add_vars.m

????.......??????6993??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\build_cost_params.m

????.......??????3145??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\compute_cost.m

????.......??????2293??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\describe_idx.m

????.......??????5177??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\display.m

????.......???????645??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\get.m

????.......??????1625??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\getN.m

????.......??????4465??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\getv.m

????.......??????2196??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\get_cost_params.m

????.......??????2168??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\get_idx.m

????.......??????4191??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\linear_constraints.m

????.......??????7901??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\opt_model.m

????.......??????1132??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\@opt_model\userdata.m

?????文件???????4559??2017-05-29?18:38??三點(diǎn)估計(jì)潮流計(jì)算matlab\ABC.m

?????文件???????5191??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\add_userfcn.m

?????文件??????17385??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\apply_changes.m

?????文件????????319??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\AUTHORS

?????文件???????2065??2016-12-17?01:34??三點(diǎn)估計(jì)潮流計(jì)算matlab\bustypes.m

?????文件??????33539??2016-12-17?01:33??三點(diǎn)估計(jì)潮流計(jì)算matlab\case118.m

?????文件?????214629??2016-12-17?01:33??三點(diǎn)估計(jì)潮流計(jì)算matlab\case1354pegase.m

?????文件????2258871??2016-12-17?01:33??三點(diǎn)估計(jì)潮流計(jì)算matlab\case13659pegase.m

?????文件???????4611??2016-12-17?01:33??三點(diǎn)估計(jì)潮流計(jì)算matlab\case14.m

?????文件??????60070??2016-12-17?01:33??三點(diǎn)估計(jì)潮流計(jì)算matlab\case145.m

?????文件?????290087??2016-12-17?01:33??三點(diǎn)估計(jì)潮流計(jì)算matlab\case1888rte.m

?????文件?????305188??2016-12-17?01:33??三點(diǎn)估計(jì)潮流計(jì)算matlab\case1951rte.m

............此處省略500個(gè)文件信息

評(píng)論

共有 條評(píng)論