⭐ 欢迎来到虫虫下载站! | 📦 资源下载 📁 资源专辑 ℹ️ 关于我们
⭐ 虫虫下载站

📄 out.txt

📁 博士师兄编写的电力系统状态估计程序
💻 TXT
📖 第 1 页 / 共 2 页
字号:
NOX=7

结点次序优化ORD()
****子系统结点数NOB=4
****子系统实际出现的最大外部结点号IMAY=4
****内部结点号-->外部结点号转换表IIO[]
IIO[0]=4
IIO[1]=2
IIO[2]=3
IIO[3]=1

****外部结点号-->程序内部结点号转换表IOI[]
IOI[0]=4
IOI[1]=2
IOI[2]=3
IOI[3]=1

****计算参考电(外部结点号)=1 


****形成导纳矩阵上三角非零元素的位置表JJY()
****合并并联支路后输出结果
****导纳矩阵上三角各行非零互导纳首址表IY[]
IY[0]=1
IY[1]=2
IY[2]=4
IY[3]=5
IY[4]=5

****导纳矩阵上三角各行非零互导纳的列号JYY[]
JYY[0]=3
JYY[1]=3
JYY[2]=4
JYY[3]=4
JYY[4]=0


****建立虚拟的导纳矩阵下三角非零元素位置表JYL()
****导纳矩阵各行下三角非零元素首址表IYL[]
IYL[0]=1
IYL[1]=1
IYL[2]=1
IYL[3]=3
IYL[4]=5

****导纳矩阵各行下三角非零元素的列号JYLL[][0]
****导纳矩阵各行下三角非零元素的导纳值的存放位置JYLL[][1]
JYLL[0][0]=1    JYLL[0][1]=1    
JYLL[1][0]=2    JYLL[1][1]=2    
JYLL[2][0]=2    JYLL[2][1]=3    
JYLL[3][0]=3    JYLL[3][1]=4    

****状态估计可计算性检查MSP()
****有功量测对各结点电压,相角覆盖次数的统计数组NPA[]
****无功或电压量测对各结点电压幅值覆盖次数的统计数组NQV[]
NPA[0]=2 NQV[0]=2 IIO[0]=4
 NPA[1]=4 NQV[1]=4 IIO[1]=2
 NPA[2]=6 NQV[2]=7 IIO[2]=3
 NPA[3]=5 NQV[3]=6 IIO[3]=1
 电压量测计数器NV=2  
有功量测计数器NP=7  
无功量测计数器NQ=7
能否进行状态估计标志NOP=0 

****形成导纳矩阵YY()
****导纳矩阵各结点的自电导GII[]  自电纳BII[]
gii[0]=0.000000,   bii[0]=-0.133333
gii[1]=0.156371,   bii[1]=-0.810317
gii[2]=0.124017,   bii[2]=-0.674202
gii[3]=0.109220,   bii[3]=-0.466871

****导纳矩阵各结点的互电导GIJ[]  互电纳BIJ[]
gij[0]=0.000000,   bij[0]=0.126984
gij[1]=-0.085584,   bij[1]=0.448796
gij[2]=-0.070787,   bij[2]=0.362102
gij[3]=-0.038433,   bij[3]=0.105470

****进行状态估计REST()
迭代收敛精度E=0.000010   

****形成量测量的稀疏存放的雅可比矩阵H()  

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

****形成雅克比矩阵的转置矩阵HT()  

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

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1 18.710026
REST   2    3    1 42.854752
REST   3    4    1 41.790001
REST   4    5    1 39.720867
REST   5    4    2 -19.100000
REST   6    5    2  4.803469
REST   7    6    2 17.549999
REST   8    7    2 -3.696531
REST   9    6    3 132.570007
REST  10    7    3 48.144154
REST  11    2    3 -49.490055
REST  12    3    3 46.702354
REST  13    4    4 -193.220001
REST  14    5    4  4.975098
REST  15    0    1  0.000000
REST  16    1    3 -1.160004
迭代次数IT=2 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

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

****形成雅克比矩阵的转置矩阵HT()  

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

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1 -1.364578
REST   2    3    1 -0.632225
REST   3    4    1  2.596371
REST   4    5    1  0.127388
REST   5    4    2  0.019020
REST   6    5    2  0.910180
REST   7    6    2 -1.662334
REST   8    7    2 -0.207407
REST   9    6    3  5.340813
REST  10    7    3 -0.182739
REST  11    2    3 -1.532787
REST  12    3    3 -14.199945
REST  13    4    4  1.178787
REST  14    5    4 -11.089199
REST  15    0    1  0.715225
REST  16    1    3 -0.692986
迭代次数IT=3 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

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

****形成雅克比矩阵的转置矩阵HT()  

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

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1 -1.889677
REST   2    3    1 -0.679245
REST   3    4    1  1.824768
REST   4    5    1  0.357368
REST   5    4    2  0.265461
REST   6    5    2  0.633536
REST   7    6    2 -1.911633
REST   8    7    2  0.066994
REST   9    6    3  3.216370
REST  10    7    3  1.012749
REST  11    2    3 -1.346710
REST  12    3    3 -0.712997
REST  13    4    4  3.738586
REST  14    5    4  0.927597
REST  15    0    1  0.685814
REST  16    1    3 -0.708389
迭代次数IT=4 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

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

****形成雅克比矩阵的转置矩阵HT()  

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

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1 -1.889374
REST   2    3    1 -0.678085
REST   3    4    1  1.825233
REST   4    5    1  0.357914
REST   5    4    2  0.265371
REST   6    5    2  0.633672
REST   7    6    2 -1.911543
REST   8    7    2  0.066861
REST   9    6    3  3.217377
REST  10    7    3  1.012943
REST  11    2    3 -1.353752
REST  12    3    3 -0.713139
REST  13    4    4  3.730469
REST  14    5    4  0.927269
REST  15    0    1  0.685806
REST  16    1    3 -0.708405
迭代次数IT=5 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

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

****形成雅克比矩阵的转置矩阵HT()  
残差的门槛值epsl=0.900000
 残差的门槛值epsl=0.900000  不良数据个数NOBD=6  
BD MBD -1.453316  2.000000  1.000000
BD MBD 1.276400  4.000000  1.000000
BD MBD -1.481833  6.000000  2.000000
BD MBD 1.584877  6.000000  3.000000
BD MBD -0.902580  2.000000  3.000000
BD MBD 1.498186  4.000000  4.000000

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

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1  0.000000
REST   2    3    1 -0.678268
REST   3    4    1  0.000000
REST   4    5    1  0.357594
REST   5    4    2  0.265392
REST   6    5    2  0.633662
REST   7    6    2  0.000000
REST   8    7    2  0.066868
REST   9    6    3  0.000000
REST  10    7    3  1.012566
REST  11    2    3  0.000000
REST  12    3    3 -0.714043
REST  13    4    4  0.000000
REST  14    5    4  0.926949
REST  15    0    1  0.685814
REST  16    1    3 -0.708397
迭代次数IT=6 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

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

****形成雅克比矩阵的转置矩阵HT()  
残差的门槛值epsl=0.900000
 残差的门槛值epsl=0.900000  不良数据个数NOBD=5  
BD MBD -1.470633  2.000000  1.000000
BD MBD 1.267971  4.000000  1.000000
BD MBD -1.473689  6.000000  2.000000
BD MBD 1.601128  6.000000  3.000000
BD MBD 1.485373  4.000000  4.000000

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

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1  0.000000
REST   2    3    1 -0.656921
REST   3    4    1  0.000000
REST   4    5    1  0.368778
REST   5    4    2  0.254995
REST   6    5    2  0.643912
REST   7    6    2  0.000000
REST   8    7    2  0.056445
REST   9    6    3  0.000000
REST  10    7    3  0.992016
REST  11    2    3 -1.342308
REST  12    3    3 -0.727329
REST  13    4    4  0.000000
REST  14    5    4  0.944695
REST  15    0    1  0.687943
REST  16    1    3 -0.706764
迭代次数IT=7 
  
****形成量测量的稀疏存放的雅可比矩阵H()  

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

****形成雅克比矩阵的转置矩阵HT()  
残差的门槛值epsl=0.900000
 残差的门槛值epsl=0.900000  不良数据个数NOBD=5  
BD MBD -1.511317  2.000000  1.000000
BD MBD 1.295695  4.000000  1.000000
BD MBD -1.401211  6.000000  2.000000
BD MBD 1.812692  6.000000  3.000000
BD MBD 1.726799  4.000000  4.000000

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

****各量测量的残差

量测量序号 量测量类型 量测地点  残差
REST   1    2    1  0.000000
REST   2    3    1 -0.655125
REST   3    4    1  0.000000
REST   4    5    1  0.357613
REST   5    4    2  0.162441
REST   6    5    2  0.657282
REST   7    6    2  0.000000
REST   8    7    2  0.045568

⌨️ 快捷键说明

复制代码 Ctrl + C
搜索代码 Ctrl + F
全屏模式 F11
切换主题 Ctrl + Shift + D
显示快捷键 ?
增大字号 Ctrl + =
减小字号 Ctrl + -