以高斯塞德尔法和牛顿拉夫逊法为思想的基于Matlab的电力系统潮流计算

        电力系统的潮流计算常见的方法有两种,一个是高斯赛德尔法,另一个是牛顿拉夫逊法,前者速度慢一些,但最终可以稳定的算出结果,而后者,虽然计算速度很快,但能否算出结果和初值有关。所以,这里我采用了两种方法结合的方式来计算潮流。首先通过高斯法算出一个大概的初值,然后再将其带入牛顿拉夫逊法中解得最终的精确解。下面的图为MATLAB的程序文件,其中各个主要功能都做成了子函数,每个功能里都有大量的注释信息,方便检阅程,最下面是我的下载地址,有问题欢迎留言讨论。

以高斯塞德尔法和牛顿拉夫逊法为思想的基于Matlab的电力系统潮流计算_第1张图片

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

function [node_result,s_result] = PowerSystem                                                         

 

 

[node] = OpenNode; 
[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.csdn.net/download/qq_24253277/9849699                                

你可能感兴趣的:(Matlab程序)