%?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