含DG的配网无功优化计算研究，包括风力发电和光伏电池微型燃气轮机和燃料电池

clear; clc; F=[0.7500 0.4000 0 0.7000 0.4500 0 0.6750 0.3500 0 0.6750 0.3500 0 0.6800 0.4500 0 0.7200 0.4000 0 0.8000 0.5000 0 0.7825 0.4500 0.1000 0.8425 0.4000 0.1500 0.8650 0.3500 0.1850 0.9000 0.4000 0.2300 0.8100 0.2500 0.2500 0.7800 0.2000 0.3000 0.8600 0.4000 0.2500 0.8750 0.3500 0.2000 0.8900 0.5000 0.1500 0.9050 0.5500 0.1000 0.8850 0.6000 0 0.9000 0.5500 0 0.9800 0.6500 0 0.9500 0.6500 0 0.9200 0.6000 0 0.8600 0.7000 0 0.8000 0.5500 0]; for T=1:24 f=F(T,:); x=[0 0]; [U,ws,S0,S1]=yue(x,f); d=0; for n=2:33 d=d+abs(U(n)-1); end csdy=d;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 优化前电压偏差 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% csws=ws;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 优化前网损 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% csgl=S0;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 优化前上级电网输入功率 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% popsize=30;%粒子数 MAXITER=100;%迭代次数 dimension=2;%维数 w_max = 0.9; w_min = 0.4; c1 = 2; c2 = 2; q1min=0;%微型燃气轮机和燃料电池有功出力下限 q1max=0.06;%微型燃气轮机和燃料电池有功出力上限 q2min=0;%微型燃气轮机和燃料电池无功出力下限 q2max=0.15;%微型燃气轮机和燃料电池无功出力上限 v=3/2*(rand(popsize,2));%初始化当前速度 Q1=(q1max-q1min)*rand(popsize,1);%初始化微型燃气轮机和燃料电池的无功出力 Q2=(q2max-q2min)*rand(popsize,1);%初始化微型燃气轮机和燃料电池的无功出力 x=[Q1 Q2]; pbest=x;%定义局部最优粒子 gbest=zeros(1,dimension);%定义全局最优粒子 Ulim=zeros(32,1); for i=1:popsize; [U,ws,S0,S1]=yue(x(i,:),f); if S0>0 b=0; elseif S0<0 b=abs(S0); end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%对连接点无功潮流进行约束处理%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 if S1>0 c=0; elseif S1<0 c=abs(S1); end d=0; for n=2:33 d=d+abs(U(n)-1); end a=0; for n=2:33 if U(n)>=1.05 Ulim(n)=1.05; elseif U(n)<=0.95 Ulim(n)=0.95; elseif 0.95<U(n)<1.05 Ulim(n)=U(n); end a=a+(U(n)-Ulim(n))^2; %对电压越限的节点采用罚函数处理 end f_pbest(i)=0.4*ws/csws+0.6*d/csdy+1000*(a+b+c);%初始状态时各粒子的适应度函数 end g=find(f_pbest==min(f_pbest(1:popsize)));%找出f_pbest中等于min(f_pbest(1:popsize))元素的下标 gbest=pbest(g,:);%得到初始状态时的全局最优粒子 f_gbest=f_pbest(g);%初始状态时种群的最小目标函数值 MINIUM=f_gbest; for t=1:MAXITER; %disp('MINIUM='); %disp(MINIUM); %disp('迭代次数t='); % disp(t); w = w_max-(w_max-w_min)*t/MAXITER;%惯性权重自适应调整 for i=1:popsize v(i,:) = w.*v(i,:)+c2.*rand.*(gbest(1:dimension)-x(i,1:dimension))+c1.*rand.*(pbest(i,1:dimension)-x(i,1:dimension));%速度更新规则 if abs(v(i,1))>0.006;%最大速度设置，粒子的范围宽度 v(i,1) = sign(v(i,1))*0.006; end if abs(v(i,2))>0.015;%最大速度设置，粒子的范围宽度 v(i,2) = sign(v(i,2))*0.015; end x(i,1:dimension) = x(i,1:dimension)+v(i,1:dimension);%更新后的粒子 if x(i,1)>0.06; x(i,1)=0.06; end if x(i,1)<0; x(i,1)=0; end if x(i,2)>0.15; x(i,2)=0.15; end if x(i,2)<0; x(i,2)=0; end [U,ws,S0,S1]=yue(x(i,:),f); if S0>0 b=0; elseif S0<0 b=abs(S0); end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%对连接点无功潮流进行约束处理%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 if S1>0 c=0; elseif S1<0 c=abs(S1); end d=0; for n=2:33 d=d+abs(U(n)-1); end a=0; for n=2:33 if U(n)>=1.05 Ulim(n)=1.05; elseif U(n)<=0.95 Ulim(n)=0.95; elseif 0.95<U(n)<1.05 Ulim(n)=U(n); end a=a+(U(n)-Ulim(n))^2; %对电压越限的节点采用罚函数处理 end f_x(i)=0.4*ws/csws+0.6*d/csdy+1000*(a+b+c);%更新后的粒子的目标函数 if f_x(i)<f_pbest(i) pbest(i,:)=x(i,:);%若粒子的目标函数小于更新后的局部最优值，则局部最优值就为这个粒子 f_pbest(i)=f_x(i); end if f_pbest(i)<f_gbest gbest=pbest(i,:);%若更新后的局部最优值小于全局最优值，则全局最优值就是局部最优值 f_gbest=f_pbest(i); end MINIUM=f_gbest; end end G(T,:)=gbest%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%得到的DG3和DG4的无功出力%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% S(T,1)=d%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%得到每个时段的电压偏差%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% M(T,1)=S0%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%得到每个时段的上级电网输入有功功率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% L(T,1)=ws%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%得到每个时段的网损%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% end SH=sum(S)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%得到24个时段电压偏差的平均值%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% LH=sum(L)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%得到全天网损之和%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%