基于 MATLAB 的牛顿拉夫逊迭代法计算潮流(附加短路计算)
按后文中给出的代码建立.m 文件放在一个文件夹里面。
先运行 Powerflow_main.m 计算算例系统的潮流;
然后运行 ShortcircuitCalc.m 计算算例系统三相短路电流;
程序说明详见各.m 文件注释部分,写的已经很详细了,慢慢看吧。
Powerflow_main.m 文件代码如下:
clear
%牛顿拉夫逊迭代法计算潮流
format short %规定参数数据显示精度
%节点参数矩阵
%第八列表示节点类型(2 表示平衡节点,1 表示 PV 节点,0 表示 PQ 节点)
0, 1.04 ,
0, 1.025,
0, 1.025,
0, 1.0,
0, 0.3, 1.137, 2;
0, 0.3, 1.211, 1;
0, 0.3, 1.047, 1;
0, 0;
0, 0;
0, 0;
0, 0;
0, 0;
0, 0
5, -1.25, -0.5, 1.0,
6, -0.9, -0.3, 1.0,
0,
-1, -0.35, 1.0,
0, 0, 1.0,
for k=1:size(Node_p,1)
if Node_p(k,8)==1