基于蒙特卡洛随机潮流研究(Matlab代码实现)

简介: 基于蒙特卡洛随机潮流研究(Matlab代码实现)

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

💥💥💞💞欢迎来到本博客❤️❤️💥💥



🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。



⛳️座右铭:行百里者,半于九十。


📋📋📋本文目录如下:🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现



💥1 概述

当所求问题的解是某个事件的概率,或者是某个随机变量的数学期望,或者是概率、数学期望有关的量时,通过某种试验的方法,得到该事件发生的频率,或者该随机变量若干个具体观察值的算数平均值,通过它得到问题的解。这就是蒙特卡罗方法的基本思想。本文利用蒙特卡洛模拟法计算随机潮流,计算电力系统安全性

📚2 运行结果

image.gif

image.gif

image.gif

image.gif

image.gif

image.gif

部分代码:

for i=1:n    %对PI节点的处理

   h=h+1;

   if i~=isb&&B2(i,2)==4  

        Q(i)=sqrt(Ig^2*((B2(i,5))^2)-B(i,3)^2);        %e=B2(i,5),f=0,e^2+f^2=B2(i,5))^2,其中e和f为光伏发电系统接入节点电压的实部和虚部

        B2(i,4)=Q(i);

        B2(i,2)=1;

        PI=h;

   end  

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

h=0;

j=0;

for i=1:n            %对PQ节点的处理

   if i~=isb&&B2(i,2)==1  

       h=h+1;

       for j=1:n

           OrgS(2*h-1,1)=OrgS(2*h-1,1)+B2(i,5)*(G(i,j)*B2(j,5)-B(i,j)*B2(j,6))+B2(i,6)*(G(i,j)*B2(j,6)+B(i,j)*B2(j,5));   %Pi 书P57页11-45

           OrgS(2*h,1)=OrgS(2*h,1)+B2(i,6)*(G(i,j)*B2(j,5)-B(i,j)*B2(j,6))-B2(i,5)*(G(i,j)*B2(j,6)+B(i,j)*B2(j,5));       %Qi  同上

       end

   end

end

for i=1:n           %对PV节点的处理

   if i~=isb&&B2(i,2)==2

       h=h+1;

       for j=1:n

           OrgS(2*h-1,1)=OrgS(2*h-1,1)+ B2(i,5)*(G(i,j)*B2(j,5)-B(i,j)*B2(j,6))+B2(i,6)*(G(i,j)*B2(j,6)+B(i,j)*B2(j,5)); %同上

           OrgS(2*h,1)=OrgS(2*h,1)+ B2(i,6)*(G(i,j)*B2(j,5)-B(i,j)*B2(j,6))-B2(i,5)*(G(i,j)*B2(j,6)+B(i,j)*B2(j,5));    %同上

       end

   end

end

OrgS;

%创建PVU 用于存储PV节点的初始电压

PVU=zeros(n-H-1,2);

t=0;

for i=1:n

   if B2(i,2)==2

       t=t+1;

       PVU(t,1)=B2(i,5);

       PVU(t,2)=B2(i,6);

   end

end

%disp('PV节点初始值:电压、相位ei、fi:')   ;

PVU;

%创建DetaS,用于存储有功功率、无功功率和电压幅值的不平衡量 书p58 11-46和11-47

h=0;

for i=1:n           %对PQ节点的处理

   if i~=isb&&B2(i,2)==1

       h=h+1;

       DetaS(2*h-1,1)=B2(i,3)-OrgS(2*h-1,1);

       DetaS(2*h,1)=B2(i,4)-OrgS(2*h,1);

   end

end

t=0;

for i=1:n           %对PV节点的处理

   if i~=isb&&B2(i,2)==2

       h=h+1;

       t=t+1;

       DetaS(2*h-1,1)=B2(i,3)-OrgS(2*h-1,1);

       DetaS(2*h,1)=PVU(t,1)^2+PVU(t,2)^2-B2(i,5)^2-B2(i,6)^2;

   end

end

%disp('P、Q、V不平衡量:')  ;

DetaS;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%创建I,用于存储节点电流参数(计算雅克比矩阵用 具体见纸板公式推导)

I=zeros(n-1,1);

h=0;

for i=1:n

  if i~=isb

       h=h+1;

       I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,5)+B2(i,6)*sqrt(-1));                                          

   end                                                          

end

I;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%创建Jacbi(雅可比矩阵)     各元素是书中各元素求法乘个-1 所以△W=J△V

Jacbi=zeros(2*n-2);

h=0;

k=0;

for i=1:n       %对PQ节点的处理

🎉3 参考文献

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

[1]雒浪,蓝文军,王明霞.一种基于蒙特卡洛的模糊随机潮流计算方法[J].电子世界,2019(04):15-16.DOI:10.19353/j.cnki.dzsj.2019.04.008.

🌈4 Matlab代码实现

https://ttaozhi.com/t/p.html?id=0jWVcWcATX

相关文章
|
1月前
|
数据可视化 算法
MATLAB Simulink 交交变流电路性能研究
MATLAB Simulink 交交变流电路性能研究
21 2
|
1月前
|
数据可视化 算法
MATLAB Simulink 直流斩波电路性能研究
MATLAB Simulink 直流斩波电路性能研究
32 1
|
1月前
|
数据可视化 算法
MATLAB Simulink 逆变电路性能研究
MATLAB Simulink 逆变电路性能研究
20 1
|
1月前
|
数据可视化 算法
MATLAB Simulink 三相整流电路性能研究
MATLAB Simulink 三相整流电路性能研究
23 1
|
1月前
|
数据可视化 算法
MATLAB Simulink 单相桥式整流电路性能研究
MATLAB Simulink 单相桥式整流电路性能研究
23 2
|
1月前
|
数据可视化 算法
MATLAB Simulink 单相半波可控整流电路性能研究
MATLAB Simulink 单相半波可控整流电路性能研究
14 2
|
1月前
|
数据安全/隐私保护
地震波功率谱密度函数、功率谱密度曲线,反应谱转功率谱,matlab代码
地震波格式转换、时程转换、峰值调整、规范反应谱、计算反应谱、计算持时、生成人工波、时频域转换、数据滤波、基线校正、Arias截波、傅里叶变换、耐震时程曲线、脉冲波合成与提取、三联反应谱、地震动参数、延性反应谱、地震波缩尺、功率谱密度
|
1月前
|
数据安全/隐私保护
耐震时程曲线,matlab代码,自定义反应谱与地震波,优化源代码,地震波耐震时程曲线
地震波格式转换、时程转换、峰值调整、规范反应谱、计算反应谱、计算持时、生成人工波、时频域转换、数据滤波、基线校正、Arias截波、傅里叶变换、耐震时程曲线、脉冲波合成与提取、三联反应谱、地震动参数、延性反应谱、地震波缩尺、功率谱密度
基于混合整数规划的微网储能电池容量规划(matlab代码)
基于混合整数规划的微网储能电池容量规划(matlab代码)
|
1月前
|
算法 调度
面向配电网韧性提升的移动储能预布局与动态调度策略(matlab代码)
面向配电网韧性提升的移动储能预布局与动态调度策略(matlab代码)

热门文章

最新文章