以高斯塞德尔法和牛顿拉夫逊法为思想的基于Matlab的电力系统潮流计算
电力系统的潮流计算常见的方法有两种,一个是高斯赛德尔法,另一个是牛顿拉夫逊法,前者速度慢一些,但最终可以稳定的算出结果,而后者,虽然计算速度很快,但能否算出结果和初值有关。所以,这里我采用了两种方法结合的方式来计算潮流。首先通过高斯法算出一个大概的初值,然后再将其带入牛顿拉夫逊法中解得最终的精确解。下面的图为MATLAB的程序文件,其中各个主要功能都做成了子函数,每个功能里都有大量的注释信息,方便检阅程,最下面是我的下载地址,有问题欢迎留言讨论。
---------------------------------------------------------------------------------------------------
function [node_result,s_result] = PowerSystem
[nn,mn] = size(node);
[line] = OpenLine;
[nl,ml] = size(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))
[node_result,s_result] =Newton_Raphson(U,Y,node,nPQ,nPV,line,nodenum);
Result_Write(node_result,s_result,node,line);
---------------------------------------------------------------------------------------------------
下载地址:https://download.****.net/download/qq_24253277/9849699