-------用matlab計算
《電子數字計算機的應用-電力系統計算》 西安交通大學、 清華大學、 浙江大學、 湖南大學、 成都工學院、 水利電力部電網排程研究所合編, 水利電力出版社， 1978。 計算第151頁例3-1。節點電壓計算結果在第156頁例表3-5中。
功率計算結果在第168頁例圖3-4中。
1. 使用雅克比迭代法進行潮流計算
為了程式設計簡便，第4節點由pv節點改為pq節點。
雅克比迭代法見第121頁。
節點方程i=yv或者v=zi, 由z=y-1得到y導納矩陣的逆矩陣,即是阻抗矩陣z。y導納矩陣見第152頁。i是節點注入電流，v是節點對地電壓。
計算步驟:
1) 求阻抗矩陣z=y-1
2) 設5個節點的電壓初值為（見第152頁例表3-1）
1, 1, 1, 1, 1.05
3) 計算各節點注入電流初值
已知節點1、2、3、4的注入功率pi、qi（見第168頁例圖3-4）
求節點1、2、3、4的注入電流初值 i=1,2,3,4
由求平衡節點的電流初值
4) 求pq節點電壓電流 i=1,2,3,4
5) 求平衡點節點電流
返回4)進行迭代
m檔案內容:
y=[1.37874-6.29166i, -0.62402+3.90015i, -0.75471+2.6415i, 0, 0;
-0.62402+3.90015i, 1.45390-66.98082i, -0.82987+3.11203i, 0+63.49206i, 0;
-0.75471+2.64150i, -0.82987+3.11203i, 1.58459-35.73786i, 0, 0+31.74603i;
0, 0+63.49206i, 0, 0-66.66667i, 0;
0, 0, 0+31.74603i, 0, 0-33.33333i]
z=inv(y)
…….2. 利用matlab中的fsolve()函式求節點電壓
matlab函式fsolve()是用於求解形如f(x)=0的非線性代數方程的。
節點導納方程i=yv；
注入功率
令節點導納元素，
令節點電壓
由上式展開得：
簡記為：
即： 下標i=1,2,3,4; j=1,2,3,4,5 因為平衡點節點電壓已知
上式就是f(x)=0形式，其中、是待求變數，，是已知的注入節點功率。用fsolve()解上面的8個方程，求出、，然後求出i=yv。
變數、用x=[e1; e2; e3; e4; f1; f2; f3; f4] 表示，其初值為x0=[1; 1; 1; 1; 0; 0; 0; 0]
計算步驟:
1. 求ai 和 bi
2. 求 f(i) = i=1, 2, 3, 4
求 f(i) = i=5, 6, 7, 8
先建立乙個m檔案，檔名為powerflow2011，包含下面的函式
function f=powerflow2008new(x)
y=[1.37874-6.29166i, -0.62402+3.90015i, -0.75471+2.6415i, 0, 0;
-0.62402+3.90015i, 1.45390-66.98082i, -0.82987+3.11203i, 0+63.49206i, 0;
-0.75471+2.64150i, -0.82987+3.11203i, 1.58459-35.73786i, 0, 0+31.74603i;
0, 0+63.49206i, 0, 0-66.66667i, 0;
0, 0, 0+31.74603i, 0, 0-33.33333i];
g=real(y);
b=imag(y);
......
潮流計算程式
b1是支路引數矩陣，第一列和第二列是起始節點編號和終點節點編號 第三列 第四列 第五列 第六列 第七列 第八列分別為 支路電阻 電抗 電導 電納 變壓器變比 是否有變壓器 1為有 0為無 b2為節點引數矩陣，其中第一列到第六列為節點編號 為節點型別 注入有功 注入無功 電壓幅值 電壓相位。1 為平衡...
潮流計算方法
本人乙個星期以來的汗水，弄清楚了計算機演算法計算潮流的基礎，如果有什麼不懂的可以發資訊到郵箱 接下來開始弄潮流的優化問題，吼吼！電力系統的潮流計算的計算機演算法 以matlab為環境 這裡理論不做過多介紹，推薦一本專門講解電力系統分析的計算機演算法的書籍電力系統分析的計算機演算法 邱曉燕 劉天琪編著...
什麼是潮流計算
潮流計算，電力學名詞，指在給定電力系統網路拓撲 元件引數和發電 負荷參量條件下，計算有功功率 無功功率及電壓在電力網中的分布。潮流計算是根據給定的電網結構 引數和發電機 負荷等元件的執行條件，確定電力系統各部分穩態執行狀態引數的計算。通常給定的執行條件有系統中各電源和負荷點的功率 樞紐點電壓 平衡點...