多微网优化调度(风机、光伏、蓄电池、燃料电池、柴油机、电网交互)(Matlab代码实现)

简介: 多微网优化调度(风机、光伏、蓄电池、燃料电池、柴油机、电网交互)(Matlab代码实现)

💥1 概述

以含有光伏(PV)、风机 (W)、柴油发电机(de)、燃料电池(fc)、蓄电池(bat)的子微网1、子微网2 构成微网群进行算例分析, 微网群算例系统拓扑结构如下图。


微网群采用并联式结构,其中子网1中分布式电源包括:PV1, W1,fc1,bat1;子网2中分布式电源包括PV2,W2,de2,bat2。微网群与配电网之间的公共连接点群PCC保持闭合,网PCC1与网PCC2也均保持闭合。微网1,2和配电网之间可存在功率交换Pbuy1,Psell1,Pbuy2,Psell2,微网1与微网2之间可存在功率交换Pex1,Pex2。


本文建立的微网群以多目标协调优化为目的,其中各个子微网优化目标与总微网群的优化目标一致。微网群优化调度周期取一天,以1h作为一个调度时段,全天共分为24个时段,且算例中采用分时电价a1,a2,b1,b2的方式。


数学模型及详细文章讲解见第4部分。

📚2 运行结果

 

部分代码:


[best index]=min(fitness);
gtsite=pop;             %个体最佳位置
gtfit=fitness;           %个体最佳适应度
popsite=pop(index,:);   %全局最佳位置
popfit=fitness(index);  %全局最佳适应度
trace(1)=popfit;
for k=1:maxgen
    w=wmax-(wmax-wmin)/maxgen*j;
    for i=1:sizepop
        %速度更新
        v(i,:)=w*v(i,:)+c*rand*(gtsite(i,:)-pop(i,:))+c*rand*(popsite-pop(i,:));
        %限制微粒速度
        v(find(v>vmax))=vmax;
        v(find(v<vmin))=vmin;
        %位置更新
        pop(i,:)=pop(i,:)+v(i,:);
        %限制微粒位置
        pop1=pop(i,1:dim);
        pop2=pop(i,dim+1:2*dim);
        pop3=pop(i,2*dim+1:3*dim);
        pop4=pop(i,3*dim+1:4*dim);
        pop5=pop(i,4*dim+1:5*dim);
        pop1(find(pop1>popmax(1)))=popmax(1);
        pop1(find(pop1<popmin(1)))=popmin(1);
        pop2(find(pop2>popmax(2)))=popmax(2);
        pop2(find(pop2<popmin(2)))=popmin(2);
        pop3(find(pop3>popmax(3)))=popmax(3);
        pop3(find(pop3<popmin(3)))=popmin(3);
       pop4(find(pop4>popmax(4)))=popmax(4);
        pop4(find(pop4<popmin(4)))=popmin(4);
         pop5(find(pop5>popmax(5)))=popmax(5);
        pop5(find(pop5<popmin(5)))=popmin(5);
        pop(i,:)=[pop1 pop2 pop3 pop4 pop5];
        for j=1:14
            pop(i,dim+j)=0;
            pop(i,4*dim+j)=0;
        end
            pop(i,dim+23)=0;
            pop(i,4*dim+23)=0;    
            pop(i,dim+24)=0;
            pop(i,4*dim+24)=0;  
        for j=1:24
            if abs(pop(i,j))>0&&abs(pop(i,j))<10
        pop(i,j)=0;
    end
        %if abs(pop(i,3*dim+j))>0&&abs(pop(i,3*dim+j))<10
      % pop(i,3*dim+j)=0;
       % end
       % if abs(pop(i,dim+j))>0&&abs(pop(i,dim+j))<10
        %xudianpin=xudianpin+10000;
       % pop(i,dim+j)=0;
        %end
        if abs(pop(i,4*dim+j))>0&&abs(pop(i,4*dim+j))<10
        %xudianpin=xudianpin+10000;
        pop(i,4*dim+j)=0;
        end
        if j<=7|j>=16
            pop(i,2*dim+j)=0;
        end
        end
    end
    for i=1:sizepop
        %计算适应度
%其他参数求解pbuy psell
pex2=-pop(i,2*dim+1:3*dim);
    bs1=Load1-pop(i,1:dim)-Ppv1-pop(i,dim+1:2*dim)-Pw1-pop(i,2*dim+1:3*dim);
   bs2=Load2-pop(i,3*dim+1:4*dim)-Ppv2-pex2-Pw2-pop(i,4*dim+1:5*dim);

🎉3 参考文献

部分理论来源于网络,如有侵权请联系删除。

[1]汤雯博. 面向配电网的多微电网协同优化调度研究[D].南京邮电大学,2020.DOI:10.27251/d.cnki.gnjdc.2020.001421.


[2]何力,吕红芳.考虑经济性的多微电网优化调度研究[J].发电技术,2018,39(05):397-404.

🌈4 Matlab代码及文章讲解

相关文章
|
8月前
|
机器学习/深度学习 传感器 算法
【红外图像】利用红外图像处理技术对不同制冷剂充装的制冷系统进行性能评估(Matlab代码实现)
【红外图像】利用红外图像处理技术对不同制冷剂充装的制冷系统进行性能评估(Matlab代码实现)
|
8月前
|
机器学习/深度学习 传感器 算法
【视频去噪】基于全变异正则化最小二乘反卷积是最标准的图像处理、视频去噪研究(Matlab代码实现)
【视频去噪】基于全变异正则化最小二乘反卷积是最标准的图像处理、视频去噪研究(Matlab代码实现)
|
8月前
|
机器学习/深度学习 传感器 算法
【大规模 MIMO 检测】基于ADMM的大型MU-MIMO无穷大范数检测研究(Matlab代码实现)
【大规模 MIMO 检测】基于ADMM的大型MU-MIMO无穷大范数检测研究(Matlab代码实现)
|
4月前
|
Perl
【MFAC】基于全格式动态线性化的无模型自适应控制(Matlab代码)
【MFAC】基于全格式动态线性化的无模型自适应控制(Matlab代码)
|
4月前
【数值分析】迭代法求方程的根(附matlab代码)
【数值分析】迭代法求方程的根(附matlab代码)
|
4月前
【数值分析】Jacobi、Seidel和Sor迭代法求解线性方程组(附matlab代码)
【数值分析】Jacobi、Seidel和Sor迭代法求解线性方程组(附matlab代码)
|
4月前
【数值分析】二分法求方程的根(附matlab代码)
【数值分析】二分法求方程的根(附matlab代码)
|
23天前
|
存储 人工智能 机器人
【Matlab】Matlab电话拨号音合成与识别(代码+论文)【独一无二】
【Matlab】Matlab电话拨号音合成与识别(代码+论文)【独一无二】
|
28天前
|
存储
【Matlab】Matlab电话拨号音合成与识别(代码+论文)【独一无二】
【Matlab】Matlab电话拨号音合成与识别(代码+论文)【独一无二】
|
3月前
|
机器学习/深度学习 算法 计算机视觉
霍夫变换车道线识别-车牌字符识别代码(matlab仿真与图像处理系列第5期)
霍夫变换车道线识别-车牌字符识别代码(matlab仿真与图像处理系列第5期)
30 2

热门文章

最新文章