[发明专利]一种基于对称CU三角分解求解电力系统节点阻抗矩阵的方法

权利要求书

1.一种基于对称CU三角分解求解电力系统节点阻抗矩阵的方法，其特征包括以下步骤：

步骤1：读取各支路数据文件；

步骤2：形成节点导纳矩阵Y；

步骤3：利用Y阵以及C、U阵元素的对称关系快速形成C、U因子阵的合成阵；

(1)建立CU合成阵，确定c、u元素与Y阵元素的关系；

(2)根据合成阵、对c_kj元素进行类消元计算，只需分步计算对角元c_kk及以右的所有u′_kj元素；

(3)在计算u_kj＝u′_kj/c_kk元素之前，先根据对称性将u′_kj元素直接赋值给相应的c_jk元素，再计算u_kj＝u′_kj/c_kk；

步骤4：根据方程CUZ_k＝E_k对方程CW_k＝E_k求取W_k阵；通过方程UZ_k＝W_k求取Z_k阵对角元Z_kk及以上的非对角元素；

(1)对方程UZ_k＝W_k求解Z_k阵元素的计算顺序为：Z_n,…,Z_k,…,Z₁，且在计算各个Z_k阵过程中，仅计算Z_k阵对角元Z_kk及其以上的非对角元素；

(2)对方程CW_k＝E_k中的W_k阵求解仅需求取W_k阵中的对角元素w_kk＝1/c_kk；

(3)对方程UZ_k＝W_k的求解仅计算Z_k阵对角元Z_kk及其以上的非对角元素；

所述E_k阵是第k行的元素为1，其余元素全部为零的矩阵；

步骤5：根据对称性求取对角元Z_kk以左的非对角元素；

所述步骤3-步骤5加快了求取电力系统Z阵元素的速度；

步骤6：将Z阵写入数据文件。