Monte 随机潮流计算

随机潮流计算，采用蒙特卡洛法进行潮流计算，负荷和风电机采用正态分布 )

% B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳
%         5、支路的变比；6、支路首端处于K侧为1，1侧为0

%输入pq节点时采用列向量的形式，如[1;2;3;4;]
%节点初始功率S，发电机节点为正，负荷节点为负
clear;
v0=2;
vci=3;
vco=20;
vr=13.5;
pr=0.6;
a=input（‘请输入节点数:a=‘）;
nl=input（‘请输入支路数:nl=‘）;
isb=input（‘请输入平衡母线节点号:isb=‘）;
usb=input（‘请输入平衡母线电压初始值:usb=‘）;
pv=input（‘请输入PV节点号:pv=‘）;
upv=input（‘请输入PV节点电压初始值:upv=‘）;
epv=input（‘请输入PV节点相角初始值:epv=‘）;
pq=input（‘请输入PQ节点号:pq=‘）;
B1=input（‘请输入由各支路参数形成的矩阵:B1=‘）;
%S=input（‘请输入节点初始功率:S=‘）;
WV=input（‘请输入风电机节点:‘）;
W=input（‘输入有功无功的期望值和方差矩阵（维数a*4）:‘）;
Mo=input（‘输入发电机节点和出力:‘）
for ii=1:30000
 S=zeros（1a）;
 wind=wblrnd（112）+v0;
 if （wind>vci&&wind<=vr） pw=pr/（vr-vci）*wind-pr/（vr-vci）*vci;end
 if （wind>vr&&wind<=vco） pw=pr;end
 if （wind>vco||wind<=vci） pw=0;end
 num=binornd（300.5）;
 totpw=pw*num;
 S（WV）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4456  2014-05-10 23:36  Monte.m

----------- ---------  ---------- -----  ----

                 4924                    3