基于MATLAB的电力系统潮流计算

本代码用于五节点系统（第五节点为平衡节点）的在直角坐标系下的潮流计算，也可以拓展为任意节点的情况。共含有四个子程序，分别为雅克比矩阵的计算、功率电压不平衡计算、中间好、过程以及最终结果计算。适用于课程设计，亲测有效。

%%%此函数为雅克比矩阵的算法
clcclear
%%电压赋值
U=input（‘输入电压各节点电压值‘）; 
U=U‘;                                      
e=real（U）;
f=imag（U）;
%%初始化
n=length（U）;                          %节点数量
Bus=[1 2 3 4 5]‘;                       %节点号
Vtype=[1 1 1 2 3]‘;                     %节点类型                              
swing=find（Vtype==3）;                   %平衡节点编号
PV=find（Vtype==2）;                      %PV节点编号
PQ=find（Vtype==1）;                      %PQ节点编号
Vset=[1 1 1 1.1 1.05];                  %电压设定值
PVPQ=find（Vtype<3）;                     %PV/PQ节点编号
%节点导纳矩阵及B，G提取
Y=[6.311-20.4022j -3.5587+11.3879j -2.7523+9.1742j 0 0;
    -3.5587+11.3879j 8.5587-31.5272j -5+15j 5.1948j 0;
    -2.7523+9.1742j -5+15j 8.5587-31.5272j 0 5.1948j
    0 5.1948j 0 -5.1948j 0;
    0 0 5.1948 0 -5.1948];
G=real（Y）

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

     文件       2140  2018-12-05 16:45  潮流计算\jacobi_matrix.m

     文件       4390  2018-12-06 12:02  潮流计算\Load_flow_calculation.m

     文件       1736  2018-12-11 09:57  潮流计算\Node_power_unbalance.m

     文件       3338  2018-12-12 18:26  潮流计算\zhongjianguocheng.m

     目录          0  2019-01-09 22:04  潮流计算

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

                11604                    5