Matlab电力系统潮流计算程序

Matlab电力系统潮流计算程序，潮流数据由系统参数、线路参数、接地支路参数、变压器参数、运行参数、PV节点参数等6大块参数组成，各块参数间以0作为区分标志。

clear all
%读入数据文件
n=load（‘节点数n.txt‘）;
nl=load（‘支路数nl.txt‘）;
pr=load（‘迭代精度.txt‘）;
kmax=load（‘最大迭代次数.txt‘）;
X=load（‘接地支路参数.txt‘）;% 矩阵X[na]各个元素意义：n--节点号a--电纳;
B1=load（‘支路参数阵.txt‘）;% 矩阵B1[ijabcde]各个元素意义：i--支路首端号Pj--支路末端号Qa--支路电阻b--支路电抗c--1/2充电电容电纳;
                         % d--支路变比Ke--折算到哪一侧的标志，如果支路的首端位于高压侧则为1，否则输入0;
PG=load（‘发电机输出有功功率.txt‘）;
QG=load（‘发电机输出无功功率.txt‘）;
PL=load（‘负荷有功功率.txt‘）;
QL=load（‘负荷无功功率.txt‘）;
V=load（‘节点电压.txt‘）;%V[i，j]各个元素意义：i--节点电压幅值，j--节点类型:1-平衡节点2-PQ节点3-PV节点;
%--------------------------------------------------------------------------
Y=zeros（n）;C1=zeros（n1）;%初始电压的相角；
%求节点导纳矩阵
for j=1:n                       %将接地支路参数加到节点导纳阵中
    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  发电机输出无功功率.txt

     文件         43  2011-10-28 10:06  发电机输出有功功率.txt

     文件         67  2011-10-28 10:10  负荷无功功率.txt

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

     文件         77  2011-10-28 08:27  接地支路参数.txt

     文件        156  2011-10-30 12:12  节点电压.txt

     文件          2  2011-10-28 08:03  节点数n.txt

     文件        592  2011-10-28 08:51  支路参数阵.txt

     文件          2  2011-10-28 08:11  支路数nl.txt

     文件          2  2011-10-28 09:56  最大迭代次数.txt

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

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

     文件     134656  2016-07-11 08:55  潮流数据格式.doc

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

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

               145302                    14