clear?all
%讀入數(shù)據(jù)文件
n=load（‘節(jié)點數(shù)n.txt‘）;
nl=load（‘支路數(shù)nl.txt‘）;
pr=load（‘迭代精度.txt‘）;
kmax=load（‘最大迭代次數(shù).txt‘）;
X=load（‘接地支路參數(shù).txt‘）;%?矩陣X[na]各個元素意義：n--節(jié)點號a--電納;
B1=load（‘支路參數(shù)陣.txt‘）;%?矩陣B1[ijabcde]各個元素意義：i--支路首端號Pj--支路末端號Qa--支路電阻b--支路電抗c--1/2充電電容電納;
?????????????????????????%?d--支路變比Ke--折算到哪一側(cè)的標志，如果支路的首端位于高壓側(cè)則為1，否則輸入0;
PG=load（‘發(fā)電機輸出有功功率.txt‘）;
QG=load（‘發(fā)電機輸出無功功率.txt‘）;
PL=load（‘負荷有功功率.txt‘）;
QL=load（‘負荷無功功率.txt‘）;
V=load（‘節(jié)點電壓.txt‘）;%V[i，j]各個元素意義：i--節(jié)點電壓幅值，j--節(jié)點類型:1-平衡節(jié)點2-PQ節(jié)點3-PV節(jié)點;
%--------------------------------------------------------------------------
Y=zeros（n）;C1=zeros（n1）;%初始電壓的相角；
%求節(jié)點導納矩陣
for?j=1:n???????????????????????%將接地支路參數(shù)加到節(jié)點導納陣中
????if?X（j2）~=0;
????????p=X（j1）;
????????Y（pp）=X（j2）*i;
????end
end
for?j=1:nl								
????if?B1（j7）=

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????????55??2011-10-28?10:07??發(fā)電機輸出無功功率.txt

?????文件?????????43??2011-10-28?10:06??發(fā)電機輸出有功功率.txt

?????文件?????????67??2011-10-28?10:10??負荷無功功率.txt

?????文件?????????67??2011-10-28?10:08??負荷有功功率.txt

?????文件?????????77??2011-10-28?08:27??接地支路參數(shù).txt

?????文件????????156??2011-10-30?12:12??節(jié)點電壓.txt

?????文件??????????2??2011-10-28?08:03??節(jié)點數(shù)n.txt

?????文件????????592??2011-10-28?08:51??支路參數(shù)陣.txt

?????文件??????????2??2011-10-28?08:11??支路數(shù)nl.txt

?????文件??????????2??2011-10-28?09:56??最大迭代次數(shù).txt

?????文件???????4888??2011-11-02?00:08??jzbcl.asv

?????文件???????4687??2011-11-04?10:56??jzbcl.m

?????文件?????134656??2016-07-11?08:55??潮流數(shù)據(jù)格式.doc

?????文件??????????8??2011-10-28?08:24??迭代精度.txt

-----------?---------??----------?-----??----

???????????????145302????????????????????14