牛顿-拉夫逊法的潮流计算程序

牛顿-拉夫逊法的潮流计算程序，并提供9节点、30节点的算例。

% Y = G + j*B;
% 导纳矩阵是对称矩阵；
function[Y G B Yfrom Yto busNumber lineNumber generatorNumber] = Admittance（baseMVA busData branchData generatorData）
    busNumber = size（busData 1）;                                          % 节点数量；
    lineNumber = size（branchData 1）;                                      % 线路数量；
    generatorNumber = size（generatorData 1）;                              % 发电机数量；
    Y = zeros（busNumber busNumber）;                                       % 定义导纳矩阵；

%% 自己编写    
    % 求取导纳矩阵的非对角元；
    % 这样即使相同两个节点之间追加了支路，也可以直接在原文件中添加，而不需要修正程序；
    for m = 1:1:lineNumber
        if branchData（m 8） == 1                                           % 该线路处于运行状态；
            if branchData（m 6） == 0                                       % Ratio = 0 ，求普通线路的导纳；
                temp = -1./（branchData（m 3） + 1i *  branchData（m 4））;
                Y（branchData（m 1） branchData（m 2）） = Y（branchData（m 1） branchData（m 2）） + temp;
                Y（branchData（m 2） branchData（m 1）） = Y（branchData（m 2） branchData（m 1）） + temp;
            else                                                           % Ratio = K ，求变压器支路对应的导纳；
                K = branchData（m 6）;
                temp = -1./（branchData（m 3） + 1i *  branchData（m 4））/K;
                Y（branchData（m 1） branchData（m 2）） = Y（branchData（m 1） branchData（m 2）） + temp;
                Y（branchData（m 2） branchData（m 1）） = Y（branchData（m 2） branchData（m 1）） + temp;
            end
        end
    end
    
    % 求取导纳矩阵的对角元；
    for m = 1:1:lineNumber
       if branchData（m 8） == 1                                            % 该线路处于运行状态；
           if branchData（m 6） == 0                                        % Ratio = 0 ，求普通线路的两端点的自导纳；
               % 对普通线路首端点；
               Y（branchData（m 1） branchData（m 1）） = Y（branchData（m 1） branchData（m 1）） + ...
                                                       1./（branchData（m 3） + 1i *  branchData（m 4）） + ...
                                                       1i *  branchData（m 5）;
                                                                           % 考虑线路自身充电电容；
                % 对普通线路末端点；
               Y（branchData（m 2） branchData（m 2）） = Y（branchData（m 2） branchData（m 2）） + ...
                                                       1./（branchData（m 3） + 1i *  branchData（m 4）） + ...
                                                       1i *  branchData（m 5）;
                                                                           % 考虑线路自身充电电容；
           else                                                            % Ratio = K ，求变压器支路的两端点的自导纳；
               K = branchData（m 6）; 
               % 考虑变压器支路的首端对地导纳；
               Y（branchData（m 1） branchData（m 1）） = Y（branchData（m 1） branchData（m 1）） + ...
                                                       1./（branchData（m 3） + 1i * branchData（m 4）） / K + ...
            

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

     文件       6487  2018-08-22 19:14  PowerFlow\Admittance.m

     文件       3383  2018-08-20 21:20  PowerFlow\DeltaP_Q_U2.m

     文件       1054  2018-08-02 20:53  PowerFlow\IEEE_30_Matpower\Branch.txt

     文件       1064  2018-08-22 20:45  PowerFlow\IEEE_30_Matpower\Bus.txt

     文件        197  2018-07-24 16:59  PowerFlow\IEEE_30_Matpower\Gen.txt

     文件         33  2018-07-23 16:39  PowerFlow\IEEE_30_Matpower\ReadMe.txt

     文件        267  2018-08-02 10:44  PowerFlow\IEEE_9\Branch.txt

     文件        290  2018-08-02 21:46  PowerFlow\IEEE_9\Bus.txt

     文件        114  2018-08-02 10:45  PowerFlow\IEEE_9\Gen.txt

     文件         61  2018-08-02 21:48  PowerFlow\IEEE_9\Readme.txt

     文件        216  2018-07-23 16:49  PowerFlow\IEEE_9\新建文本文档.txt

     文件      12215  2018-07-22 10:43  PowerFlow\Jacobi.m

     文件        673  2018-07-21 18:30  PowerFlow\JudgeConvergence.m

     文件       5640  2018-09-21 16:15  PowerFlow\main.m

     文件       2430  2018-08-20 21:22  PowerFlow\ReadData.m

     文件       2125  2018-09-21 16:02  PowerFlow\ReadMe.txt

     文件       1242  2018-07-21 18:22  PowerFlow\RenewV.m

     文件       4629  2018-08-22 19:40  PowerFlow\Results.m

     目录          0  2018-09-21 16:11  PowerFlow\IEEE_30_Matpower

     目录          0  2018-09-21 16:00  PowerFlow\IEEE_9

     目录          0  2018-09-21 16:11  PowerFlow

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

                42120                    21