[发明专利]一种基于Givens正交增量行变换的大电网状态估计方法

权利要求书

1.一种基于Givens正交增量行变换的大电网状态估计方法，其特征在于：所述方法包括以下步骤：

步骤1：进行Givens正交增量行变换；

步骤2：进行因子表局部修正；

步骤3：进行大电网状态估计。

2.根据权利要求1所述的基于Givens正交增量行变换的大电网状态估计方法，其特征在于：所述步骤1包括以下步骤：

步骤1-1：Givens正交行变换；

设实数c与s满足：

c²+s²＝1        （1）

其中，T_ij为Givens正交变换矩阵，由Givens正交变换矩阵所确定的线性变换即为Givens变换，且i＜j；

设x＝(a₁,a₂,…,a_i,…,a_j,…,a_n)^T，其中a_i，a_j不全为0，则选择Givens正交变换矩阵T_ij，使：

T_ijx＝y≡(a₁,a₂,…,a_i,,…,a_j,,…,a_n)^T（3）

则：

c=ai/ai2+aj2s=aj/ai2+aj2---(4)]]>

经过正交变换后的结果为：

ait=ai2+aj2ajt=0---(5)]]>

H_R为加权雅可比矩阵，该矩阵为实数矩阵且非奇异，则存在正交矩阵T₁，T₂，…，T_n-1一系列Givens正交变换矩阵使：

其中，T_k＝T_k,k,T_k,k-1,…,T_k,i,…T_k,1，实现对第k+1行下三角元素的消元，T_k,i实现对(k+1,i)元素消元；通过T₁,T₂,…,T_k,…,T_n-1一系列Givens正交变换矩阵，实现H_R逐行消元，并形成三角矩阵R，称三角矩阵R为因子表；

步骤1-2：Givens正交增量行变换；

当加权雅可比矩阵H_R发生局部变化，变为矩阵H_R′时，基于H_R的Givens正交行变换结果，在因子表R基础上构造增量矩阵d＝[λ₁d₁;…;λ_id_i;…]，通过对增量矩阵d继续分解形成新的分解结果R′，实现加权雅可比矩阵H_R发生变化时因子表R的快速修正；

加权雅可比矩阵H_R为：

HR=a11a12a13···a1na21a22a23···a2na31a32a33···a3na41a42a43···a4n···············am1am2am3···am4---(7)]]>

其中，m＞n；

加权雅可比矩阵H_R分解结果R为：

加权雅可比矩阵H_R发生局部变化变为矩阵H_R′，对应分解结果R变为R′，通过Givens正交增量行变换实现如下：

TdRd=R'---(9)]]>

其中，T_d为分解增量矩阵d所需的系列Givens正交变换矩阵。