out.txt

博士师兄编写的电力系统状态估计程序

REST   9    6    3  0.000000
REST  10    7    3  1.000504
REST  11    2    3 -0.218079
REST  12    3    3 -0.727947
REST  13    4    4  0.000000
REST  14    5    4  0.946091
REST  15    0    1  0.687500
REST  16    1    3 -0.705566
迭代次数IT=8 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

****子系统量测数NOM=16

****形成雅克比矩阵的转置矩阵HT()  
残差的门槛值epsl=0.900000
 残差的门槛值epsl=0.900000  不良数据个数NOBD=5  
BD MBD -1.520760  2.000000  1.000000
BD MBD 1.298440  4.000000  1.000000
BD MBD -1.388542  6.000000  2.000000
BD MBD 1.847945  6.000000  3.000000
BD MBD 1.765148  4.000000  4.000000

****迭代结果输出
不包含不良数据的目标函数值jx=2.476327  
包含不良数据的目标函数值jb=14.933681 
自由矢量最大的节点AEMAX
外部结点号=4  自由矢量值a=0.000090
MaxResidual= 0.000000  量测类型=4  量测地点=4 

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1  0.000000
REST   2    3    1 -0.653603
REST   3    4    1  0.000000
REST   4    5    1  0.356197
REST   5    4    2  0.146217
REST   6    5    2  0.659818
REST   7    6    2  0.000000
REST   8    7    2  0.043331
REST   9    6    3  0.000000
REST  10    7    3  1.001434
REST  11    2    3 -0.034737
REST  12    3    3 -0.728951
REST  13    4    4  0.000000
REST  14    5    4  0.946701
REST  15    0    1  0.687477
REST  16    1    3 -0.705322
迭代次数IT=9 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

****子系统量测数NOM=16

****形成雅克比矩阵的转置矩阵HT()  
残差的门槛值epsl=0.900000
 残差的门槛值epsl=0.900000  不良数据个数NOBD=5  
BD MBD -1.524944  2.000000  1.000000
BD MBD 1.297173  4.000000  1.000000
BD MBD -1.385650  6.000000  2.000000
BD MBD 1.854289  6.000000  3.000000
BD MBD 1.770504  4.000000  4.000000

****迭代结果输出
不包含不良数据的目标函数值jx=2.475656  
包含不良数据的目标函数值jb=14.976867 
自由矢量最大的节点AEMAX
外部结点号=3  自由矢量值a=0.000032
MaxResidual= 0.000000  量测类型=4  量测地点=4 

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1  0.000000
REST   2    3    1 -0.653542
REST   3    4    1  0.000000
REST   4    5    1  0.355583
REST   5    4    2  0.142546
REST   6    5    2  0.660334
REST   7    6    2  0.000000
REST   8    7    2  0.042972
REST   9    6    3  0.000000
REST  10    7    3  1.001724
REST  11    2    3 -0.004848
REST  12    3    3 -0.728794
REST  13    4    4  0.000000
REST  14    5    4  0.946968
REST  15    0    1  0.687454
REST  16    1    3 -0.705276
迭代次数IT=10 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

****子系统量测数NOM=16

****形成雅克比矩阵的转置矩阵HT()  
残差的门槛值epsl=0.900000
 残差的门槛值epsl=0.900000  不良数据个数NOBD=5  
BD MBD -1.528231  2.000000  1.000000
BD MBD 1.295277  4.000000  1.000000
BD MBD -1.384435  6.000000  2.000000
BD MBD 1.855769  6.000000  3.000000
BD MBD 1.770571  4.000000  4.000000

****迭代结果输出
不包含不良数据的目标函数值jx=2.475687  
包含不良数据的目标函数值jb=14.984381 
自由矢量最大的节点AEMAX
外部结点号=4  自由矢量值a=0.000023
MaxResidual= 0.000000  量测类型=4  量测地点=4 

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1  0.000000
REST   2    3    1 -0.654030
REST   3    4    1  0.000000
REST   4    5    1  0.355152
REST   5    4    2  0.140978
REST   6    5    2  0.660499
REST   7    6    2  0.000000
REST   8    7    2  0.042805
REST   9    6    3  0.000000
REST  10    7    3  1.001564
REST  11    2    3  0.000000
REST  12    3    3 -0.729065
REST  13    4    4  0.000000
REST  14    5    4  0.946945
REST  15    0    1  0.687454
REST  16    1    3 -0.705246
迭代次数IT=11 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

****子系统量测数NOM=16

****形成雅克比矩阵的转置矩阵HT()  
残差的门槛值epsl=0.900000
 残差的门槛值epsl=0.900000  不良数据个数NOBD=5  
BD MBD -1.531286  2.000000  1.000000
BD MBD 1.293353  4.000000  1.000000
BD MBD -1.383538  6.000000  2.000000
BD MBD 1.856393  6.000000  3.000000
BD MBD 1.769872  4.000000  4.000000

****迭代结果输出
不包含不良数据的目标函数值jx=2.475838  
包含不良数据的目标函数值jb=14.986257 
自由矢量最大的节点AEMAX
外部结点号=4  自由矢量值a=0.000014
MaxResidual= 0.000000  量测类型=4  量测地点=4 

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1  0.000000
REST   2    3    1 -0.653847
REST   3    4    1  0.000000
REST   4    5    1  0.355370
REST   5    4    2  0.139822
REST   6    5    2  0.660617
REST   7    6    2  0.000000
REST   8    7    2  0.042687
REST   9    6    3  0.000000
REST  10    7    3  1.002121
REST  11    2    3  0.000568
REST  12    3    3 -0.728687
REST  13    4    4  0.000000
REST  14    5    4  0.947006
REST  15    0    1  0.687447
REST  16    1    3 -0.705231
迭代次数IT=12 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

****子系统量测数NOM=16

****形成雅克比矩阵的转置矩阵HT()  
残差的门槛值epsl=0.900000
 残差的门槛值epsl=0.900000  不良数据个数NOBD=5  
BD MBD -1.533963  2.000000  1.000000
BD MBD 1.291499  4.000000  1.000000
BD MBD -1.382758  6.000000  2.000000
BD MBD 1.856777  6.000000  3.000000
BD MBD 1.769180  4.000000  4.000000

****迭代结果输出
不包含不良数据的目标函数值jx=2.475635  
包含不良数据的目标函数值jb=14.986285 
自由矢量最大的节点AEMAX
外部结点号=4  自由矢量值a=0.000042
MaxResidual= 0.000000  量测类型=4  量测地点=4 

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1  0.000000
REST   2    3    1 -0.653179
REST   3    4    1  0.000000
REST   4    5    1  0.355564
REST   5    4    2  0.138874
REST   6    5    2  0.660951
REST   7    6    2  0.000000
REST   8    7    2  0.042511
REST   9    6    3  0.000000
REST  10    7    3  1.001976
REST  11    2    3  0.000717
REST  12    3    3 -0.728863
REST  13    4    4  0.000000
REST  14    5    4  0.947075
REST  15    0    1  0.687439
REST  16    1    3 -0.705231
迭代次数IT=13 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

****子系统量测数NOM=16

****形成雅克比矩阵的转置矩阵HT()  
残差的门槛值epsl=0.900000
 残差的门槛值epsl=0.900000  不良数据个数NOBD=5  
BD MBD -1.536780  2.000000  1.000000
BD MBD 1.289661  4.000000  1.000000
BD MBD -1.382044  6.000000  2.000000
BD MBD 1.857032  6.000000  3.000000
BD MBD 1.768555  4.000000  4.000000

****迭代结果输出
不包含不良数据的目标函数值jx=2.475914  
包含不良数据的目标函数值jb=14.987235 
自由矢量最大的节点AEMAX
外部结点号=4  自由矢量值a=0.000012
MaxResidual= 0.000000  量测类型=4  量测地点=4 

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1  0.000000
REST   2    3    1 -0.653667
REST   3    4    1  0.000000
REST   4    5    1  0.355118
REST   5    4    2  0.137951
REST   6    5    2  0.660982
REST   7    6    2  0.000000
REST   8    7    2  0.042479
REST   9    6    3  0.000000
REST  10    7    3  1.002228
REST  11    2    3  0.000618
REST  12    3    3 -0.728943
REST  13    4    4  0.000000
REST  14    5    4  0.947021
REST  15    0    1  0.687439
REST  16    1    3 -0.705208
迭代次数IT=14 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

****子系统量测数NOM=16

****形成雅克比矩阵的转置矩阵HT()  
残差的门槛值epsl=0.900000
 残差的门槛值epsl=0.900000  不良数据个数NOBD=5  
BD MBD -1.539365  2.000000  1.000000
BD MBD 1.287855  4.000000  1.000000
BD MBD -1.381380  6.000000  2.000000
BD MBD 1.857175  6.000000  3.000000
BD MBD 1.768028  4.000000  4.000000

****迭代结果输出
不包含不良数据的目标函数值jx=2.475974  
包含不良数据的目标函数值jb=14.987424 
自由矢量最大的节点AEMAX
外部结点号=3  自由矢量值a=0.000010
MaxResidual= 0.000000  量测类型=4  量测地点=4 

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1  0.000000
REST   2    3    1 -0.653973
REST   3    4    1  0.000000
REST   4    5    1  0.354980
REST   5    4    2  0.137096
REST   6    5    2  0.661092
REST   7    6    2  0.000000
REST   8    7    2  0.042370
REST   9    6    3  0.000000
REST  10    7    3  1.002178
REST  11    2    3  0.000404
REST  12    3    3 -0.728985
REST  13    4    4  0.000000
REST  14    5    4  0.946999
REST  15    0    1  0.687431
REST  16    1    3 -0.705200
最终迭代次数IT=14  

****状态估计结果输出

****结点电压结果输出IVE
外部结点号 内部结点号 电压幅值 电压相角（度）
1  4  111.462570   0.000000
2  2  110.339951  -0.399733
3  3  111.695198   0.797972
4  1  112.322144   7.923705

****支路电流结果输出
支路号 首端节点 尾端节点  电流实部  电流虚部
1 1 2  0.358400 -0.352985
2 1 3 -0.172588 -0.036392
3 2 3 -1.158896  0.405560
4 3 4 -1.857138 -0.057949
****子系统结点数NOB=4 

****结点注入结果输出
外部结点号  内部结点号  结点有功PH  结点无功QH
1  4  20.711267  34.703850
2  2  -167.998703  -89.589302
3  3  -49.490395  -41.200871
4  1  197.622375  101.350227

****支路功率结果输出
支路号 I  J   PIJ   QIJ   PJI   QJI
1 1 2 39.948383 37.605106 -39.816799 -40.375828
2 1 3 -19.237087 -2.901070 19.331976 -10.782376
3 2 3 -128.181885 -49.213066 128.799973 41.607925
4 3 4 -197.622375 -72.027023 197.622375 101.350159
形成回路时间 0.000000 秒