以高斯塞德爾法和牛頓拉夫遜法為思想的基于Matlab的電力系統潮流計算

        電力系統的潮流計算常見的方法有兩種，一個是高斯賽德爾法，另一個是牛頓拉夫遜法，前者速度慢一些，但最終可以穩定的算出結果，而後者，雖然計算速度很快，但能否算出結果和初值有關。是以，這裡我采用了兩種方法結合的方式來計算潮流。首先通過高斯法算出一個大概的初值，然後再将其帶入牛頓拉夫遜法中解得最終的精确解。下面的圖為MATLAB的程式檔案，其中各個主要功能都做成了子函數，每個功能裡都有大量的注釋資訊，友善檢閱程，最下面是我的下載下傳位址，有問題歡迎留言讨論。

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

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