✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。
🍎个人主页:Matlab科研工作室
🍊个人信条:格物致知。
更多Matlab仿真内容点击👇
⛄ 内容介绍
提出了基于人工鱼群优化算法(AFSA)的最优潮流(OPF)计 算方法;算法结合动态调整罚函数的方式,将最优潮流问题转化为一个无约束求极值问题,有效提高了算法的全局收敛能力和计算精度.应用此算法对标准 IEEE6节点的电力系统进行最优潮流计算,仿真结果表明,该算法能够更好地获得全局最优解,具有实用意义.
⛄ 部分代码
function [operate,overflowvalue]=dcflow(N,L,BL,Pi,nbl)
isb=2; %默认的平衡节点号(该程序平衡节点不能为1)
Y=zeros(N);
for i=1:L
p=BL(i,1);q=BL(i,2);
Y(p,q)=-1./BL(i,4);
Y(q,p)=Y(p,q);
Y(q,q)=Y(q,q)+1./BL(i,4);
Y(p,p)=Y(p,p)+1./BL(i,4);
end
for i=1:N
j=N+1;
Y(i,j)=Pi(i,1)-Pi(i,2);
end
Y(:,isb)=[];Y(isb,:)=[];%去掉参考节点对应的行和列
for j=2:N %高斯法解线性方程组
Y(1,j)=Y(1,j)/Y(1,1);
end
Y(1,1)=1;
for i=2:N-1
for k=1:i-1
for j=k+1:N
Y(i,j)=Y(i,j)-Y(i,k)*Y(k,j);
end
Y(i,k)=0;
end
for j=i+1:N
Y(i,j)=Y(i,j)/Y(i,i);
end
Y(i,i)=1;
end
thta(N-1)=Y(N-1,N);
for i=N-2:-1:1
S=0;
for j=i+1:N-1
S=S+Y(i,j)*thta(j);
thta(i)=Y(i,N)-S;
end
end
for i=N:-1:isb
thta(i)=thta(i-1);
end
thta(isb)=0;
for i=1:L %计算支路功率
PL(i)=(thta(BL(i,1))-thta(BL(i,2)))/BL(i,4);
end
for i=1:L
if nbl(i)~=0
PL(i)=PL(i)/(nbl(i)+1);
end
end
operate=0;
for i=1:L %计算运行费用
if nbl(i)==0
operate=operate+BL(i,3)*PL(i)^2; %operat记录方案的运行费用
else
operate=operate+(nbl(i)+1)*BL(i,3)*PL(i)^2;
end
end
operate=operate*0.35;
overflownumber=0;
overflowvalue=0;
for i=1:L
if abs(PL(i))>BL(i,5)
overflownumber=overflownumber+1;
overflow(overflownumber,1)=BL(i,1);%overflow矩阵储存过负荷支路首末节点号
overflow(overflownumber,2)=BL(i,2);
overflowvalue=overflowvalue+abs(PL(i))-BL(i,5);%计算过负荷费用
end
end
overflowvalue=overflowvalue*100;
⛄ 运行结果
⛄ 参考文献
[1] 洪彬倬,聂一雄,盛丹,et al.基于改进人工鱼群算法的含风电场电力系统最优潮流计算[J].电力科学与技术学报, 2012, 27(3):6.DOI:10.3969/j.issn.1673-9140.2012.03.009.
[2] 刘耀年,李迎红,张冰冰,等.基于人工鱼群算法的最优潮流计算[C]//电工电能新技术.CNKI;WanFang, 2006:5.DOI:10.3969/j.issn.1003-3076.2006.04.007.
⛳️ 代码获取关注我
❤️部分理论引用网络文献,若有侵权联系博主删除
❤️ 关注我领取海量matlab电子书和数学建模资料