潮流计算matlab代码

电力系统的潮流计算代码，内含39节点，22节点等

function dPQ=dPQ（Ybus）
%delP--有功偏差量
%delQ--无功偏差量
%Y--节点导纳矩阵
%bus--节点参数（PQU及相角）矩阵
global n;
global m;
delP=zeros（1n-1）;
delQ=zeros（1m）;
for i=1:（n-1）  %形成delP矩阵（PQ、PV节点）
    P1=0;T=0;
    for j=1:n
        T=bus（i2）*bus（j2）*（real（Y（ij））*cos（bus（i3）-bus（j3））+imag（Y（ij））*sin（bus（i3）-bus（j3）））;
        P1=P1+T;
    end
    delP（1i）=bus（i4）-P1;
end
for i=1:m %形成delQ矩阵（PQ节点）
    Q1=0;T=0;
    for j=1:n
        T=bus（i2）*bus（j2）*（real（Y（ij））*sin（bus（i3）-bus（j3））-imag（Y（ij））*cos（bus（i3）-bus（j3）））;
        Q1=Q1+T;
    end
    delQ（1i）=bus（i5）-Q1;
end
dPQ=[delPdelQ];


    

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

     文件        986  2008-08-26 13:13  loadflow（2）\y.m

     文件        669  2008-08-26 15:58  loadflow（2）\dPQ.m

     文件       5240  2008-08-28 20:10  loadflow（2）\main.m

     文件       1704  2008-08-27 14:21  loadflow（2）\loss.m

     文件        505  2008-08-22 16:26  loadflow（2）\test11.m

     文件        443  2008-08-22 16:26  loadflow（2）\test12.m

     文件       1469  2008-08-24 23:17  loadflow（2）\test14.m

     文件       1797  2008-08-22 16:26  loadflow（2）\test13.m

     文件       9564  2008-08-22 16:26  loadflow（2）\典型系统1的计算结果.txt

     文件       5348  2008-08-22 16:26  loadflow（2）\典型系统2的计算结果.txt

     文件      20946  2008-08-22 16:26  loadflow（2）\典型系统3的计算结果.txt

     文件      34525  2008-08-22 16:26  loadflow（2）\典型系统4的计算结果.txt

     文件       3155  2008-08-22 16:26  loadflow（2）\test15.m

     文件       3952  2008-08-22 16:26  loadflow（2）\test16.m

     文件       7057  2008-08-27 14:30  loadflow（2）\output1.dat

     文件       1314  2008-08-26 15:58  loadflow（2）\jac.m

     目录          0  2008-08-27 14:32  loadflow（2）

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

                98674                    17