牛顿拉夫逊算法潮流计算IEEE14节点

IEEE标准节点系统简化版，在理论分析中采用IEEE标准节点系统来进行模拟时耗时就较多，若仍采用此法来计算实际运行的复杂大电网系统，其计算量将达到难以想象的程度。而且，目前研究人员已经研发出了一系列较为成熟的开断计算法的，如补偿法或灵敏系数法，其原理均是在已知的正常潮流计算的基础上采用增加补偿功率或是节点注入功率增量来模拟支路开断。以上方法在本质上都是模拟近似方法，在开断潮流模拟计算的过程中均保持了雅克比矩阵不变，并能确保一定的计算速度与计算准确程度，从大类上也可以被归类为定常雅克比矩阵的开断潮流计算。

clear all
n=14 ;
m =20 ;
ph=1 ;
k=1;
B1=[120.01938+0.05917j2*0.0264j11;
    230.04699+0.1979j2*0.0219j12;
    240.05811+0.17632j2*0.0187j13;
    150.05403+0.22304j2*0.0246j14;
    250.05695+0.17388j2*0.017j15;
    340.06701+0.17103j2*0.0173j16;
    450.01335+0.04211j2*0.0064j17;
    560.25202j00.9328;
    470.20912j00.9789;
    780.17615j0110;
    490.55618j00.96911;
    790.11001j0112;
    9100.03181+0.0845j0113;
    6110.09498+0.1989j0114;
    6120.12291+0.25581j0115;
    6130.06615+0.13027j0116;
    9140.12711+0.27038j0117;
    10110.08205+0.19207j0118;
    12130.22092+0.19988j0119;
    13140.17093+0.34802j0120];
% 它以矩阵形式存贮支路的情况 每行存贮一条支路
% 第一列存贮支路的一个端点
% 第二列存贮支路的另一个端点
% 第三列存贮支路的阻抗
% 第四列存贮支路的对地导纳
% 第五列存贮变压器的变比 注意支路为1
% 第六列存贮支路的序号
B2= [001.063;
    0.40.2171.0452;
    00.9421.012;
    0（47.8-3.9j）/10011;
    0（7.6+1.6j）/10011;
    00.1121.072;
    0011;
    001.092;
    0（29.5+16.6j）/10011;
    0（9+5.8j）/10011;
    0（3.5+1.8j）/10011;
    0（6.1+1.6j）/10011;
    0（13.5+5.8j）/10011;
    0（14.9+5j）/10011];
% 第一列为电源侧的功率
% 第二列为负荷侧的功率
% 第三列为该点的电压值
% 第四列为该点的类型: 1为PQ节点2为PV节点3为平衡节点
B3=[2-0.40.5;
    300.4;
    6-0.060.24;
    8-0.060.24];
%B3为无功可调母线数据矩阵，第一列为节点号，第二列为下限值，第三列为上限值
B5=[10;20.127j;30.19j;60.075j;80];%%%%%%%%%%  PV节点跟