牛顿-拉弗逊潮流计算matlab算法

基于牛顿拉弗逊法的潮流计算方法，通过多个子程序调用的新式实现，程序流程清晰

function Y = Admittance（busline）
[nbmb]=size（bus）;
[nlml]=size（line）;
Y=zeros（nbnb）;                        % 对导纳矩阵赋初值0
for k=1:nl   
    I=line（k1）;  
    J=line（k2）;   
    Zt=line（k3）+j*line（k4）;          % 读入线路参数
    if Zt~=0
        Yt=1/Zt;                       % 接地支路不计算Yt
    end
    Ym=line（k5）+j*line（k6）;   
    K=line（k7）;       
   
    if （K==0）&（J~=0）                   % 普通线路: K=0
        Y（II）=Y（II）+Yt+Ym;       
        Y（JJ）=Y（JJ）+Yt+Ym;       
        Y（IJ）=Y（IJ）-Yt;       
        Y（JI）=Y（IJ）;   
    end
   
    if （K==0）&（J==0）                   % 对地支路: K=0J=0R=X=0
        Y（II）=Y（II）+Ym;   
    end   
   
    if K>0                             % 变压器线路: Zt和Ym为折算到i侧的值K在j侧       
        Y（II）=Y（II）+Yt+Ym;       
      

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

     文件       1196  2018-12-26 09:33  牛顿拉弗逊法潮流计算\Admittance.m

     文件        292  2018-12-26 09:36  牛顿拉弗逊法潮流计算\bus_res_.m

     文件        679  2018-12-26 09:29  牛顿拉弗逊法潮流计算\CalculatedPQ.m

     文件       1995  2018-12-26 09:30  牛顿拉弗逊法潮流计算\Jacbi.m

     文件       2246  2018-12-26 10:23  牛顿拉弗逊法潮流计算\main.m

     文件       1052  2018-12-26 09:28  牛顿拉弗逊法潮流计算\NumSort.m

     文件        595  2018-12-26 09:30  牛顿拉弗逊法潮流计算\PQCalculate.m

     文件        815  2018-12-26 09:32  牛顿拉弗逊法潮流计算\ReNum.m

     文件         30  2018-12-25 19:39  牛顿拉弗逊法潮流计算\Result.m

     文件       1034  2018-12-26 10:13  牛顿拉弗逊法潮流计算\ShowBusLine.m

     文件        879  2018-12-26 10:08  牛顿拉弗逊法潮流计算\ShowResult.m

     文件        945  2018-12-26 09:35  牛顿拉弗逊法潮流计算\S_res_.m

     文件       1119  2018-12-26 09:32  牛顿拉弗逊法潮流计算\YtYm.m

     文件         63  2018-12-26 10:31  牛顿拉弗逊法潮流计算\使用说明.txt

     目录          0  2018-12-26 10:31  牛顿拉弗逊法潮流计算

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

                12940                    15