灵璧县住房和城乡建设局网站,城乡建设部网站稽查执法专栏,上海工商查询系统官网,网站专题页面开发目录 1 主要内容
目标函数
程序模型
2 部分代码
3 程序结果 1 主要内容
该程序参考《热电联产在区域综合能源系统中的定容选址研究》#xff0c;主要针对电热综合能源系统进行优化#xff0c;确定热电联产机组的位置和容量#xff0c;程序以33节点电网和17节点热网为例…目录 1 主要内容
目标函数
程序模型
2 部分代码
3 程序结果 1 主要内容
该程序参考《热电联产在区域综合能源系统中的定容选址研究》主要针对电热综合能源系统进行优化确定热电联产机组的位置和容量程序以33节点电网和17节点热网为例具体见程序模型采用5个备选位置作为优化基础热网采用牛拉法潮流电网潮流采用matpower工具箱【发消息】-【免费软件】提供免费下载及相应教程优化算法采用遗传算法程序采用matlab编写注释清晰方便学习 目标函数 该程序以网损为目标函数通过遗传算法进行优化求解。 程序模型 2 部分代码
clc;clear all;
T60;%仿真代数 遗传算法
M40;% 群体规模
pm0.01;pc0.6;%交叉变异概率
pcmax0.9;pcmin0.6;pmmax0.1;pmmin0.01;
xmax[5 130];%选址位置、容量kw最大值
xmin[1 10];%选址位置、容量kw最小值
Long5;%单个参数字串长度总编码长度3L
bestv-inf;
bvalround(rand(M,2*Long));%初始种群 round 四舍五入取整函数 编码后的初始种群
%%%%%%%%%%%%%bestv-inf%最优适应度初值%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for ii1:Tf_ave0;
%
解码计算适应度 for a1:Myzeros(1,2); for gn1:2for j1:1:Long y(gn)y(gn)bval(a,gn*Long-j1)*2^(j-1) ; %%%解码 二进制转换为十进制 end xd(gn)round((xmax(gn)-xmin(gn))*y(gn)/(2^Long-1)xmin(gn));%%%%%%%%解码换算为约束条件范围内的数 endxx(a,:)xd;%
%%%
%%修正网络参数进行潮流计算%
%%%
%%%
%%%
obj(a)usi_pf(xd); %目标函数
end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%种群M%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% func1./obj;%目标函数转换为适应度函数
fmaxmax(func);
f_avefunc/M;
pfunc./sum(func);
qcumsum(p);%累加 分布采样区间 以便进行选择
[fmax,indmax]max(func);%求当代最佳个体 indmax 当前最佳个体代数 if fmaxbestv bestvfmax;%到目前为止最优适应度值 bvalxxbval(indmax,:);%到目前为止最佳位串 optxxxx(indmax,:);%到目前为止最优参数 end Bfit1(ii)bestv ;% 存储每代的最优适应度
%%%%遗传操作开始
%轮盘赌选择 for i1:(M-1) rrand; tmpfind(rq);%%找到适应度大于随机r的个体 newbval(i,:)bval(tmp(1),:); %
%保存第一个这样的个体因为r就落到第一个这样的数的区间内然后继续选择 end newbval(M,:)bvalxx;%最优保留 bvalnewbval;
%两点交叉
for i1:2:(M-1) %隔一个交叉一对
funmaxmax(func(i),func(i1));if funmaxelsepcpcmax-(pcmax-pcmin)*(func-f_ave)/(bestv-f_ave);
endccrand; if cc%%pc交叉率3 程序结果