IEEE3机9节点潮流计算matlab程序

压缩包里有IEEE3机9节点的潮流计算matlab程序，并附上节点描述表和支路描述表，只要修改这两个就行

branch_dscb=load（‘branch_desc.txt‘）;%导入支路描述表
node_dscb=load（‘node_desc.txt‘）;%导入节点描述表
branch_num=size（branch_dscb1）;%支路数
node_num=size（node_dscb1）;%节点数
pq_num=0;
for k=1:node_num    %统计pq节点数
    if node_dscb（k6）==1
         pq_num=pq_num+1;
    end
end

Y=zeros（node_num）;%定义节点导纳矩阵

for k=1:branch_num
    head=branch_dscb（k1）;%首节点写入
    tail=branch_dscb（k2）;%末节点写入
    r=branch_dscb（k3）;%电阻写入
    x=branch_dscb（k4）;%电抗写入
    ban=branch_dscb（k5）;%对地电纳写入
    K=branch_dscb（k6）;%变比写入
    
    %==节点导纳写入==
    Y（headhead）=Y（headhead）+1/（r+1j*x）+1j*ban/2;
    Y（headtail）=Y（headtail）-1/（K*（r+1j*x））;
    Y（tailhead）=Y（tailhead）-1/（K*（r+1j*x））;
    Y（tailtail）=Y（tailtail）+1/（K*K*（r+1j*x））+1j*ban/2;

end
G=real（Y）;%实部为电导
B=imag（Y）;%虚部为电纳
precision=1;%精度初始值设

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

     文件        642  2019-04-15 16:12  chaoliujisuan\branch_desc.txt

     文件        396  2019-04-15 16:12  chaoliujisuan\jiedian.txt

     文件       4042  2019-04-15 16:11  chaoliujisuan\loadflowcal.m

     目录          0  2019-04-15 16:13  chaoliujisuan

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

                 5080                    4