电力系统潮流（Matlab实现）

搭建一个基于PowerWorldd的复杂模型，如图所示，同时也表征着其正常运行状态下的潮流分布情况。

5.2 负荷增加

当节点Five上的负荷增加至530MW时，系统的潮流发生较大改变，如图所示，多条线路处于过负荷运行状态下，电压水平也降低了很多，说明在过负荷下会严重影响系统的电压水平。

5.3 发电机出力增加

当节点Six上的发电机有功出力增加至800MW时，系统的潮流发生较大改变，如图所示，多条线路处于过负荷运行状态下，电压水平也稍有降低。

6 MATLAB编程实例

---

6.1 潮流计算程序流程图

基于牛顿-拉夫逊法的潮流计算程序流程图如图所示。

6.2 潮流计算代码

function [node_result,s_result] = PowerSystem % 潮流计算主程序

%%

[node] = OpenNode;

[nn,mn] = size(node); % 打开数据文件.并返回node

%%

[line] = OpenLine;

[nl,ml] = size(line); % 打开数据文件.并返回line

%%

[node,line,nPQ,nPV,nodenum,PH,PV,PQ] = Num(node,line); % 对节点重新排序

%%

Y = sparse(Yij(node,line)) % 计算节点导纳矩阵

%%

[U] = abs(Gauss_Seidel(Y,node,nPQ,nPV)) % 返回GS算法的结果，作为初值

%%

[node_result,s_result] =Newton_Raphson(U,Y,node,nPQ,nPV,line,nodenum); % 用牛顿拉夫逊法计算潮流结果

%%

Result_Write(node_result,s_result,node,line); % 把结果写入文件

function [node] = OpenNode

[dfile,pathname]=uigetfile(‘*.m’,‘Select Node File’); % 数据文件类型为m文件，窗口标题为“Select Node File”

if pathname == 0

error(’ you must select a valid data file’) % 如果没有选择有效文件，则出现错误提示

else

lfile =length(dfile); % 去除文件后缀

eval(dfile(1:lfile-2)); % 打开文件

end