基于遗传算法二进制编码的混合发电系统配置优化问题求解matlab仿真,含仿真操作录像

%********混合发电系统的成本目标函数 function [CT]=CalCostTatal(Apv,Aw,Cb) %*********************************假设设备的置换成本与购入成本相同********************** %*********************************以下定义光伏组件的成本参数**************************** %**************光伏组件单位面积产生的功率为0.8 kW/m2,单位KW的买入成本为7000$/KW aPV=7000*0.8;%光伏组件单位面积的购入成本7000$/KW*0.8 kW/m2=5600$/m2 %**************光伏组件单位面积产生的功率为0.8 kW/m2,单位KW每年的运行和维护成本为50$/KW/year aOMPV=50*0.8;%为光伏组件单位面积每年的运行和维护成本(美元）50$/KW/year*0.8 kW/m2=40$/m2/year sPV=2000;%光伏组件单位面积的二手市场上可售出价格2000$/m2 %*********************************以下定义风力发电机的成本参数**************************** %**************风力发电机组件单位KW的买入成本为3000$/KW,风力机单位面积产生的功率为0.3KW/m2 aW=3000*0.3;%为项目初期购买风力机组件单位面积的购入和配件置换成本3000$/KW*0.3KW/m2=900$/m2 %**************风力发电机组件单位面积产生的功率为0.3 kW/m2,单位KW每年的运行和维护成本为500$/KW/year aOMW=50*0.8;%风力发电机组件单位面积每年的运行和维护成本(美元）500$/KW/year*0.3 kW/m2=150$/m2/year sW=350;%风力发电机组件单位面积的二手市场上可售出价格350$/m2 %*********************************以下定义蓄电池组件的成本参数**************************** %**************蓄电池组件单位KWH的买入成本为138.889$/KWH,置换成本为138.889$/KWH aB=138.889;%为项目初期购买风力机组件单位面积的购入和配件置换成本138.889$/KWH %**************蓄电池组件单位KWH运行和维护成本为9.26$/KWH/year aOMB=9.26;%蓄电池组件单位KWH运行和维护成本为9.26$/KWH/year sB=0;%蓄电池组件单位KWH的二手市场上可售出价格0$/m2,即忽略不计 es=0.08;%es为运行和维护的人工成本增长率 r=0.06;%r为金融存款利率 j=0.045;%j为通货膨胀率 N=20 ;%N为项目生命周期(年） Lw=15;%Lw为风力发电机的生命周期（年） Lb=10;%Lb为蓄电池组件的生命周期（年） years=N;%使用年限 Xw=round(N/Lw); Xb=round(N/Lb); %************************************以下计算成本系数c1~c6****************** %************************************* a1=(1+es)/(1+r);%等比数据首项 q=(1+es)/(1+r);%等比数据公比 if q~=1 Sn1=a1*(1-q^N)/(1-q);%等比数据求和 else Sn1=a1*N;%等比数据求和 end %************************************* a1=1;%等比数据首项 q=((1+es)/(1+r))^Lw;%等比数据公比 if q~=1 Sn2=a1*(1-q^Xw)/(1-q);%等比数据求和 else Sn2=a1*Xw;%等比数据求和 end %************************************* a1=((1+j)/(1+r))^Lw;%等比数据首项 q=((1+j)/(1+r))^Lw;%等比数据公比 if q~=1 Sn3=a1*(1-q^(Xw-1))/(1-q);%等比数据求和 else Sn3=a1*(Xw-1);%等比数据求和 end %************************************* a1=((1+j)/(1+r))^Lb;%等比数据首项 q=((1+j)/(1+r))^Lb;%等比数据公比 if q~=1 Sn4=a1*(1-q^(Xb-1))/(1-q);%等比数据求和 else Sn4=a1*(Xb-1);%等比数据求和 end %*********光伏组件成本系数 c1=aPV; c2=aOMPV*Sn1; c3=sPV*(1+j)^N/(1+r)^N; %*********风力发电机成本系数 c4=aW*Sn2; c5=aOMW*Sn2; c6=sW*Sn3+((sW-aW)/Lw*years+aW)*((1+j)/(1+r))^N; %*********蓄电池组件成本系数 c7=aB*Sn4; c8=aOMB*Sn1; CT=(c1+c2-c3)*Apv+(c4+c5-c6)*Aw+(c7+c8)*Cb;