改進(jìn)多目標(biāo)粒子群算法計(jì)算微電網(wǎng)多目標(biāo)調(diào)度 Multisim仿真

ID:122912 · 發(fā)表于 2016-10-21 16:29

以微電網(wǎng)為模型，以成本和網(wǎng)絡(luò)損耗為優(yōu)化目標(biāo)，采用多目標(biāo)粒子群算法對(duì)電網(wǎng)進(jìn)行優(yōu)化。仿真結(jié)果表明，收斂性較好。代碼直接可運(yùn)行。

改進(jìn)多目標(biāo)粒子群算法計(jì)算微電網(wǎng)多目標(biāo)調(diào)度.zip (11.7 KB, 下載次數(shù): 22)

部分源碼：

function cost=shiyingzhi_zonghexiaoyi_all_2(PGT)
%適應(yīng)值函數(shù)
%--------------------------------------------------------------------------
% PGT=[ 2.7238 0.3086 0.3078 0.3076 0.6724 0.8129 0.8239 0.8293 0.8240 0.8131 0.6770 0.8127 0.8150 0.8182 0.8286 0.6529 0.8358 2.6527 4.6283 4.6611 1.6757 1.6930 0.8180 0.6920
% 4.0000 3.4152 3.3845 3.4100 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000
% 4.0000 0 0 0 0 0.8463 2.3894 2.8585 2.3982 0.8902 0 0.8038 1.3032 1.7781 2.8045 0 3.3269 4.0000 4.0000 4.0000 4.0000 4.0000 1.3820 0
% 0.2762 0.2762 0.3077 0.2824 0.3276 0.3408 0.2867 0.3122 0.2778 0.2967 0.3230 0.3836 0.3818 0.4037 0.3669 0.3471 0.3373 0.3473 0.3717 0.3389 0.3243 0.3070 0.3000 0.3080];
T=24;
N=4;
% 輸入原始數(shù)據(jù),包括各時(shí)段負(fù)荷大小，發(fā)電機(jī)有功輸出上下限，發(fā)電機(jī)耗量成本系數(shù)，各時(shí)段風(fēng)電場(chǎng)預(yù)測(cè)的平均輸出功率
pmax(1)=7;pmin(1)=0;
pmax(2)=4;pmin(2)=0;
pmax(3)=4;pmin(3)=0;
% a=[1000,970,700,680,450,370,480,660,665,670];
% b=[16.19,17.26,16.60,16.50,19.70,22.26,27.74,25.92,27.27,27.79];
% c=[0.00048,0.00031,0.002,0.00211,0.00398,0.00712,0.00079,0.00413,0.00222,0.00173];
% 柴油發(fā)電機(jī)的耗量參數(shù)
a=0.4333;b=0.2333;c=0.0071;
% 微型燃?xì)廨啓C(jī),P_mt表示燃?xì)廨啓C(jī)發(fā)出的功率，Xl_mt表示燃?xì)廨啓C(jī)的效率
Price_mt=0.1;
% Xl_mt=0.0753*(P_mt/65)^3-0.3095*(P_mt/65)^2+0.1068;
% 燃料電池
Xl_fc=0.4;
Price_fc=0.2;
% Cost_fc=Price_fc*P_fc/Xl_fc;
% 污染物排放成本
% aa表示NOX，bb表示FC,cc表示DG
aa=[21.8,0.03,0.44];
bb=[0.454,0.006,0.008];
cc=[1.432,1.078,1.596];
% PL=[11,4,4,4,5,6,7.5,8,7.5,6,5,6,6.5,7,8,5,8.5,11,13,13,10,10,6.5,5];
PL=[3,4,4,4,5,6,6.5,7,7.5,8.5,9,10,10.5,10,9,8.5,9,10,11,11.5,10,9,5.5,5];
W=[37.35,43.52,41.81,41.54,47.29,42.10,44.74,33.11,23.98,22.70,30.76,29.23,24.25,45.37,26.82,27.71,47.12,54.10,44.52,49.51,50.67,44.35,35.56,17.67];
%--------------------------------------------------------------------------
C=zeros(1,T);
Cost=0;
chengben=0;
wuranwu=0;
zonghexiaoyi=0;
s1=zeros(1,T);
s2=zeros(1,T);
s=zeros(1,T);
for jj=1:T
s1(jj)=a+b*PGT(1,jj)+c*PGT(1,jj)*PGT(1,jj)+Price_fc*PGT(2,jj)/Xl_fc+Price_mt*PGT(3,jj)/(0.0753*(PGT(3,jj)/65)^3-0.3095*(PGT(3,jj)/65)^2+0.4174*PGT(3,jj)/65+0.1068);
s2(jj)=4.2*(aa(1)*PGT(1,jj)+aa(2)*PGT(2,jj)+aa(3)*PGT(3,jj))/1000+0.99*(bb(1)*PGT(1,jj)+bb(2)*PGT(2,jj)+bb(3)*PGT(3,jj))/1000+0.014*(cc(1)*PGT(1,jj)+cc(2)*PGT(2,jj)+cc(3)*PGT(3,jj))/1000;
% s=a+b*PGT(1,jj)+c*PGT(1,jj)*PGT(1,jj)+Price_fc*PGT(2,jj)/Xl_fc+Price_mt*PGT(3,jj)/(0.0753*(PGT(3,jj)/65)^3-0.3095*(PGT(3,jj)/65)^2+0.4174*PGT(3,jj)/65+0.1068)+4.2*(aa(1)*PGT(1,jj)+aa(2)*PGT(2,jj)+aa(3)*PGT(3,jj))/1000+0.99*(bb(1)*PGT(1,jj)+bb(2)*PGT(2,jj)+bb(3)*PGT(3,jj))/1000+0.014*(cc(1)*PGT(1,jj)+cc(2)*PGT(2,jj)+cc(3)*PGT(3,jj))/1000;
% Cost=Cost+s;
s(jj)=s1(jj)+s2(jj);
Cost=Cost+s(jj);
chengben=chengben+s1(jj);
wuranwu=wuranwu+s2(jj);
end
xx=1:T;
figure(4);
plot(xx,s1,'-bdiamond',xx,s2,'-rsquare',xx,s,'-k*'); % 繪制24小時(shí)的各微電源的功率曲線圖
xlabel('t/h');ylabel('目標(biāo)成本/
);
cost=Cost;

復(fù)制代碼

ID:338675 · 發(fā)表于 2018-5-26 13:39

支持一下樓主

帳號(hào)		自動(dòng)登錄	找回密碼
密碼			立即注冊(cè)

久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

改進(jìn)多目標(biāo)粒子群算法計(jì)算微電網(wǎng)多目標(biāo)調(diào)度 Multisim仿真

評(píng)分