[發明專利]一種基于CU三角分解求取電力系統節點阻抗矩陣的方法

有效

• 申請號: 201410696600.6
• 申請日: 2014-11-27
• 公開（公告）號: CN104657337B
• 公開（公告）日: 2018-10-19
• 發明（設計）人: 陳懇;羅仁露;席小青;萬新儒
• 申請（專利權）人: 南昌大學
• 主分類號: G06F17/16
• 分類號: G06F17/16
• 代理公司: 南昌新天下專利商標代理有限公司 36115
• 代理人: 施秀瑾

1.一種基于的CU三角分解求取電力系統節點阻抗矩陣的方法，其特征包括以下步驟：

步驟1：讀取各支路數據文件；

步驟2：形成節點導納矩陣Y；

步驟3：利用Y陣以及C、U陣元素的對稱關系快速形成C、U因子陣的合成陣；

(1)建立CU合成陣，確定c、u元素與Y陣元素的關系；

(2)根據合成陣、按“Г”方式，確定對角元以右非零元素的位置，再根據對稱性確定對角元以下非零元素的位置，按類高斯消元的方式分步計算“Г”所包含的上三角中的所有元素，而下三角的所有元素均根據對稱性獲得；

(3)對c_kj元素進行類消元計算，只需分步計算對角元c_kk及以右的所有u_kj元素(c_kk～u_kn元素)，而對角元c_kk以左的元素均不計算，可由對稱性獲得；

步驟4：根據方程CUZ_k＝E_k對方程CW_k＝E_k求取W_k陣；通過方程UZ_k＝W_k求取Z_k陣對角元Z_kk及以上的非對角元素；

(1)Z_k陣的計算順序為：Z_n,…,Z_k,…,Z₁，且在計算各個Z_k陣中，僅計算對角元Z_kk及其以上的元素；

(2)對方程CW_k＝E_k求解W_k陣時僅求取其對角元及其以上元素；

(3)對方程UZ_k＝W_k求解Z_k陣時僅求取其對角元及其以上元素；

步驟5：根據對稱性求取對角元Z_kk以左的非對角元素；

這種根據Z陣元素的對稱性得到對角元Z_kk以左元素的計算方式可減少50％非對角元素的計算；也可不求Z陣下三角的非對角元素，以進一步加快計算速度；

步驟6：將Z陣寫入數據文件。