[发明专利]基于Matlab矩阵运算的快速分解法系数矩阵计算法

权利要求书

1.基于Matlab矩阵运算的快速分解法系数矩阵计算法，其特征在于：包括以下步骤：

C1、读支路首节点号数组I、末节点号数组J、电抗数组X；

所述的首节点号数组I、末节点号数组J、电抗数组X分别按顺序存放所有支路的首节点号i_m、末节点号j_m、电抗x_m；

支路包括输电线路和变压器支路，为了区分二者，变压器支路非标准变比k_m侧的节点号加个负号；

C2、用支路电抗计算支路电纳数组B_B；

计算支路电纳数组B_B的公式如下：

B_B＝-1·/X (1)

式中，“·/”表示两数组对应元素相除；电抗数组X按顺序存放所有支路的电抗x_m；

C3、由支路电纳数组B_B形成系数矩阵B'部分非对角元素组成的矩阵B1；

系数矩阵B'部分非对角元素组成的矩阵B1为：

B1＝sparse(abs(I),abs(J),-B_B,n,n) (2)

式中，sparse为Matlab形成稀疏矩阵函数，其参数分别为矩阵的行号数组、列号数组、元素值数组、行数、列数，abs为Matlab求数组元素绝对值函数，I、J分别为支路首、末节点号数组，其中变压器支路的非标准变比k_m侧的节点号为负数，n为节点数；

C4、形成系数矩阵B'的非对角矩阵；

系数矩阵B'的非对角矩阵B′_OD为：

B′_OD＝B1+B1^T (3)

式中，上标T表示矩阵的转置；

C5、形成系数矩阵B'的对角元素组成的n维列向量BD；

BD＝-B′_OD×1_n×1 (4)

式中，BD为系数矩阵B'的对角元素组成n维列向量，1_n×1为元素都为1的n维列向量；

C6、生成系数矩阵B'的对角矩阵B′_D；

B′_D＝sparse(1:n,1:n,BD,n,n) (5)

式中，1:n表示形成从1到n、步长为1的数组；

C7、计算系数矩阵B'；

修正方程系数矩阵B'为：

B′＝B′_OD+B′_D (6)

结束。