实时电价机制下交直流混合微网优化运行Matlab代码实现

👨‍🎓个人主页：研学社的博客

💥💥💞💞欢迎来到本博客❤️❤️💥💥 🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。 ⛳️座右铭：行百里者，半于九十。 📋📋📋本文目录如下：🎁🎁🎁 目录 💥1 概述 📚2 运行结果 🎉3 参考文献 🌈4 Matlab代码实现

💥1 概述

文献来源：

摘要：为了提高交直流混合微网经济效益,提出实时电价机制下的交直流混合微网优化运行方法。首先针对交直流混合微网的结构和电气特性建立其优化运行模型;其次针对微网中负荷构成建立了基于负荷分类的负荷水平对实时电价的响应模型。所提出的日前调度模型将交直流混合微网制定的网内实时电价作为控制变量,给出实时电价制定策略和交直流混合微网日前调度方案。对所建立的模型采用 混沌粒子群 算法求解。最后通过一个算例验证了所提出方法可以提高风光消纳率,增加交直流混合微网收益并且减小用户平均电费支出。

关键词：

交直流混合微网;优化运行;实时电价;需求侧响应;混沌粒子群算法;

📚2 运行结果

部分代码“

function load=DR3(pload,detapr)

   %pload=ones(1,24);

   price=[0.4042,0.4042,0.4042,0.4042,0.4042,0.4042,0.4042,0.8084,0.8084,1.3339,1.3339,1.3339,1.3339,1.3339,1.3339,0.8084,0.8084,1.3339,1.3339,1.3339,1.3339,0.8084,0.4042,0.4042];

   %detapr=[-0.1.*ones(1,12) 0.1.*ones(1,12)];

   pl1=0.6.*pload;%刚性负荷

   pl2=0.2.*pload;%Ⅰ类

   pl3=0.2.*pload;%Ⅱ类

   [row,colz]=find(detapr>0);%电价上升时段

   [row,colf]=find(detapr<0);%电价下降时段

   num_colz= size (colz,2);%电价上升时段数量

   num_colf=size(colf,2);%%电价下降时段数量

   T=24;

   a=0.08;%死区阈值

   %Ⅰ类负荷

   linj=0;lini=0;

   for i=1:num_colf

       linj=linj+abs(detapr(colf(i)));%公式（2）

   end

   for i=1:num_colz

       lini=lini+abs(detapr(colz(i)));

   end

   for t=1:T
if ismember(t,colz)==1

           fdetap=fhzy(detapr(t),a);

           linp=0;

           for kj=1:num_colf

               linp=linp+fdetap*abs(detapr(colf(kj)))*pl2(t)/linj;

           end

            lam_pl2(t)=pl2(t)-linp;%公式（5）

          elseif ismember(t,colf)==1

           fdetap=fhzy(detapr(t),a);

           linp=0;

           for ki=1:num_colz

               linp=linp+fdetap*abs(detapr(colz(ki)))*pl2(t)/lini;

           end

            lam_pl2(t)=pl2(t)+linp;%公式（6）

       end

   end

   %Ⅱ类负荷

   e0=-1;e1=0.8;e2=0.6;

   Est=zeros(T,T);

   for i=1:T

       for j=1:T

           if  i==j

               Est(i,j)=-1;

           elseif i==j+1

               Est(i,j)=0.8;

           elseif i+1==j

               Est(i,j)=0.6;

           else

               Est(i,j)=0;

           end

       end

   end

   for t=1:T

       lam_pl3(t)=pl3(t)+pl3(t)*(Est(t,:)*(detapr./price)');%公式（10）

   end

   load=pl1+lam_pl2+lam_pl3;

🎉3 文献来源

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

[1]陈安伟,华浩瑞,李鹏等.实时电价机制下交直流混合微网优化运行方法[J].电力系统保护与控制,2017,45(07):13-20.

🌈4 Matlab代码实现

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删