5step.txt

牛顿-拉夫逊法潮流计算的算法实现

**************************优化编号后的节点编号：*****************************
原节点-------->初始编号--------->新编号
1-------->1--------->1
2-------->2--------->2
3-------->3--------->3
4-------->4--------->4
5-------->5--------->5


**************************新编号并经排序后的非接地支路数据*************************************
1    2    0.04      0.25      0.25      
1    3    0.1       0.35      0         
2    3    0.08      0.3       0.25      
-2   4    0         0.015     1.05      
-3   5    0         0.03      1.05      

**************************新编号下的发电机数据*************************************
4    -1   5    0    1.05 
5    0    0    0    1.05 

**************************新编号下的负荷数据*************************************
1    1    -1.6 -0.8 1    
2    1    -2   -1   1    
3    1    -3.7 -1.3 1    

P-theta迭代过程中用于求解修正方程式的因子表B'为:
-0.152866 -0.596202 -0.403798                     
-0.0146673-0.0687442-0.931256           
          
-0.0276926-0.12087            
          
          
-0.26061  
Q-V迭代过程中用于求解修正方程式的因子表B''为:
-0.151351 -0.605405 -0.432432                     
-0.015413 -0.0780369                    
          
-0.0289484                    
          
          
0         
------------------------------------------------------------------
第0次迭代时计算得各节点注入有功功率为:
0         
0         
0         
0         
0         

第0次迭代时,用第一因子表计算的常数项delt_P/V为:
-1.6
-2
-3.7
4.7619

第0次迭代后，由第一因子表计算得delt_theta为：
0.0945502
-0.305801
0.07995
-0.380801

------------------------------------------------------------------
第0次迭代时计算得各节点注入无功功率为:
0.312918  
-6.60863  
-2.71235  
6.9411    
3.52314   

第0次迭代时,用第二因子表计算的常数项delt_Q/V为:
-1.11292
5.60863
1.41235
0

第0次迭代后，由第二因子表计算得delt_V为：
0.104117
-0.0790343
-0.0381016
0

------------------------------------------------------------------

第0次迭代后最大功率误差为: 5.60863
第0次迭代后各节点功率及电压值为:
电压值V   相角      有功P     无功Q     
0.895883  -0.09455020         0.312918  
1.07903   0.305801  0         -6.60863  
1.0381    -0.07995  0         -2.71235  
1.05      0.362668  0         0         
1.05      0         0         0         

------------------------------------------------------------------
第1次迭代时计算得各节点注入有功功率为:
-1.65606  
-1.03152  
-3.89478  
4.08853   
2.76359   

第1次迭代时,用第一因子表计算的常数项delt_P/V为:
0.0625753
-0.897544
0.187636
0.868067

第1次迭代后，由第一因子表计算得delt_theta为：
-0.0175149
-0.00862384
-0.00695312
-0.0222959

------------------------------------------------------------------
第1次迭代时计算得各节点注入无功功率为:
-0.665873 
-1.06859  
-1.33258  
1.74109   
2.23941   

第1次迭代时,用第二因子表计算的常数项delt_Q/V为:
-0.149715
0.0635704
0.0313859
0

第1次迭代后，由第二因子表计算得delt_V为：
0.0234046
0.00049732
0.00102674
0

------------------------------------------------------------------

第1次迭代后最大功率误差为: 0.968481
第1次迭代后各节点功率及电压值为:
电压值V   相角      有功P     无功Q     
0.872478  -0.0749997-1.65606  -0.665873 
1.07854   0.313794  -1.03152  -1.06859  
1.03707   -0.0732521-3.89478  -1.33258  
1.05      0.383902  0         0         
1.05      0         0         0         

------------------------------------------------------------------
第2次迭代时计算得各节点注入有功功率为:
-1.57213  
-2.04348  
-3.67789  
5.03685   
2.53      

第2次迭代时,用第一因子表计算的常数项delt_P/V为:
-0.0319385
0.0403166
-0.0213231
-0.0350979

第2次迭代后，由第一因子表计算得delt_theta为：
0.00815406
0.00446266
0.00151335
0.00501545

------------------------------------------------------------------
第2次迭代时计算得各节点注入无功功率为:
-0.757513 
-1.00336  
-1.31121  
1.77096   
2.27727   

第2次迭代时,用第二因子表计算的常数项delt_Q/V为:
-0.0486966
0.00311496
0.0108056
0

第2次迭代后，由第二因子表计算得delt_V为：
0.00778726
0.000434191
0.000356353
0

------------------------------------------------------------------

第2次迭代后最大功率误差为: 0.043483
第2次迭代后各节点功率及电压值为:
电压值V   相角      有功P     无功Q     
0.864691  -0.0843456-1.57213  -0.757513 
1.0781    0.309656  -2.04348  -1.00336  
1.03672   -0.0747113-3.67789  -1.31121  
1.05      0.379126  0         0         
1.05      0         0         0         

------------------------------------------------------------------
第3次迭代时计算得各节点注入有功功率为:
-1.60137  
-1.99616  
-3.69416  
4.98902   
2.57942   

第3次迭代时,用第一因子表计算的常数项delt_P/V为:
0.0015878
-0.00356358
-0.00562839
0.0104617

第3次迭代后，由第一因子表计算得delt_theta为：
-0.00133064
-0.00176378
-9.00124e-005
-0.00192855

------------------------------------------------------------------
第3次迭代时计算得各节点注入无功功率为:
-0.791869 
-0.999057 
-1.29733  
1.80084   
2.28889   

第3次迭代时,用第二因子表计算的常数项delt_Q/V为:
-0.00940349
-0.00087495
-0.00257257
0

第3次迭代后，由第二因子表计算得delt_V为：
0.00158382
0.000117386
0.000207024
0

------------------------------------------------------------------

第3次迭代后最大功率误差为: 0.0109848
第3次迭代后各节点功率及电压值为:
电压值V   相角      有功P     无功Q     
0.863107  -0.0828068-1.60137  -0.791869 
1.07799   0.311292  -1.99616  -0.999057 
1.03651   -0.0746245-3.69416  -1.29733  
1.05      0.380962  0         0         
1.05      0         0         0         

------------------------------------------------------------------
第4次迭代时计算得各节点注入有功功率为:
-1.59686  
-2.0063   
-3.69724  
5.00286   
2.57591   

第4次迭代时,用第一因子表计算的常数项delt_P/V为:
-0.00364219
0.00584085
-0.00266704
-0.00272518

第4次迭代后，由第一因子表计算得delt_theta为：
0.000537259
-0.000100854
0.000100597
-5.79322e-005

------------------------------------------------------------------
第4次迭代时计算得各节点注入无功功率为:
-0.79646  
-0.999644 
-1.29976  
1.80846   
2.29603   

第4次迭代时,用第二因子表计算的常数项delt_Q/V为:
-0.00410127
-0.000329791
-0.000227459
0

第4次迭代后，由第二因子表计算得delt_V为：
0.000677811
4.83686e-005
6.42793e-005
0

------------------------------------------------------------------

第4次迭代后最大功率误差为: 0.00629635
第4次迭代后各节点功率及电压值为:
电压值V   相角      有功P     无功Q     
0.862429  -0.0834292-1.59686  -0.79646  
1.07794   0.311385  -2.0063   -0.999644 
1.03645   -0.0747216-3.69724  -1.29976  
1.05      0.381017  0         0         
1.05      0         0         0         

------------------------------------------------------------------
第5次迭代时计算得各节点注入有功功率为:
-1.59975  
-2.00099  
-3.69915  
4.99989   
2.5791    

第5次迭代时,用第一因子表计算的常数项delt_P/V为:
-0.000285842
0.0009227
-0.000822831
0.000109263

第5次迭代后，由第一因子表计算得delt_theta为：
-6.34814e-005
-0.000181403
2.41638e-006
-0.000183124

------------------------------------------------------------------
第5次迭代时计算得各节点注入无功功率为:
-0.799096 
-0.999876 
-1.29965  
1.81171   
2.29817   

第5次迭代时,用第二因子表计算的常数项delt_Q/V为:
-0.00104777
-0.000115123
-0.000337895
0

第5次迭代后，由第二因子表计算得delt_V为：
0.00017737
1.34702e-005
2.45908e-005
0

------------------------------------------------------------------

第5次迭代后最大功率误差为: 0.000994612
第5次迭代后各节点功率及电压值为:
电压值V   相角      有功P     无功Q     
0.862252  -0.0833556-1.59975  -0.799096 
1.07792   0.311554  -2.00099  -0.999876 
1.03642   -0.0747239-3.69915  -1.29965  
1.05      0.381192  0         0         
1.05      0         0         0         

------------------------------------------------------------------
第6次迭代时计算得各节点注入有功功率为:
-1.59971  
-2.00067  
-3.69969  
5.00026   
2.57912   

第6次迭代时,用第一因子表计算的常数项delt_P/V为:
-0.000334482
0.000621155
-0.000300903
-0.000248829

第6次迭代后，由第一因子表计算得delt_theta为：
3.74225e-005
-2.86051e-005
8.28635e-006
-2.4686e-005

------------------------------------------------------------------
第6次迭代时计算得各节点注入无功功率为:
-0.79964  
-0.999957 
-1.29993  
1.81259   
2.29901   

第6次迭代时,用第二因子表计算的常数项delt_Q/V为:
-0.000417766
-4.03248e-005
-6.29273e-005
0

第6次迭代后，由第二因子表计算得delt_V为：
6.96658e-005
5.1217e-006
7.71379e-006
0

------------------------------------------------------------------

第6次迭代后最大功率误差为: 0.000669558
第6次迭代后各节点功率及电压值为:
电压值V   相角      有功P     无功Q     
0.862182  -0.083399 -1.59971  -0.79964  
1.07792   0.31158   -2.00067  -0.999957 
1.03641   -0.0747319-3.69969  -1.29993  
1.05      0.381215  0         0         
1.05      0         0         0         

------------------------------------------------------------------
第7次迭代时计算得各节点注入有功功率为:
-1.59995  
-2.00015  
-3.6999   
5.00002   
2.57937   

第7次迭代时,用第一因子表计算的常数项delt_P/V为:
-5.50233e-005
0.000140101
-9.55849e-005
-1.95583e-005

第7次迭代后，由第一因子表计算得delt_theta为：
-1.80242e-006
-1.77726e-005
9.47078e-007
-1.74645e-005

------------------------------------------------------------------
第7次迭代时计算得各节点注入无功功率为:
-0.799895 
-0.999986 
-1.29996  
1.81293   
2.29926   

第7次迭代时,用第二因子表计算的常数项delt_Q/V为:
-0.00012176
-1.32763e-005
-3.48064e-005
0

第7次迭代后，由第二因子表计算得delt_V为：
2.0549e-005
1.5537e-006
2.72833e-006
0

------------------------------------------------------------------

第7次迭代后最大功率误差为: 0.000151017
第7次迭代后各节点功率及电压值为:
电压值V   相角      有功P     无功Q     
0.862162  -0.0833969-1.59995  -0.799895 
1.07792   0.311597  -2.00015  -0.999986 
1.03641   -0.0747328-3.6999   -1.29996  
1.05      0.381232  0         0         
1.05      0         0         0         

------------------------------------------------------------------
第8次迭代时计算得各节点注入有功功率为:
-1.59997  
-2.00007  
-3.69997  
5.00002   
2.5794    

第8次迭代时,用第一因子表计算的常数项delt_P/V为:
-3.23344e-005
6.41402e-005
-3.32727e-005
-2.25956e-005

第8次迭代后，由第一因子表计算得delt_theta为：
2.86837e-006
-3.99285e-006
7.57972e-007
-3.63697e-006

------------------------------------------------------------------
第8次迭代时计算得各节点注入无功功率为:
-0.799961 
-0.999995 
-1.29999  
1.81303   
2.29936   

第8次迭代时,用第二因子表计算的常数项delt_Q/V为:
-4.48286e-005
-4.51884e-006
-8.47357e-006
0

第8次迭代后，由第二因子表计算得delt_V为：
7.50142e-006
5.56465e-007
8.7799e-007
0

------------------------------------------------------------------

第8次迭代后最大功率误差为: 6.91378e-005
第8次迭代后各节点功率及电压值为:
电压值V   相角      有功P     无功Q     
0.862154  -0.0834003-1.59997  -0.799961 
1.07792   0.3116    -2.00007  -0.999995 
1.03641   -0.0747335-3.69997  -1.29999  
1.05      0.381235  0         0         
1.05      0         0         0         

------------------------------------------------------------------
第9次迭代时计算得各节点注入有功功率为:
-1.59999  
-2.00002  
-3.69999  
5         
2.57942   

第9次迭代时,用第一因子表计算的常数项delt_P/V为:
-7.42581e-006
1.73254e-005
-1.08036e-005
-3.72115e-006

第9次迭代后，由第一因子表计算得delt_theta为：
1.13494e-007
-1.81229e-006
1.45696e-007
-1.75369e-006

------------------------------------------------------------------
第9次迭代时计算得各节点注入无功功率为:
-0.799988 
-0.999998 
-1.3      
1.81307   
2.29939   

第9次迭代时,用第二因子表计算的常数项delt_Q/V为:
-1.39398e-005
-1.49259e-006
-3.6201e-006
0

第9次迭代后，由第二因子表计算得delt_V为：
2.34721e-006
1.76625e-007
3.01735e-007
0

------------------------------------------------------------------

第9次迭代后最大功率误差为: 1.86753e-005
第9次迭代后各节点功率及电压值为:
电压值V   相角      有功P     无功Q     
0.862152  -0.0834004-1.59999  -0.799988 
1.07792   0.311602  -2.00002  -0.999998 
1.03641   -0.0747337-3.69999  -1.3      
1.05      0.381237  0         0         
1.05      0         0         0         

------------------------------------------------------------------
第10次迭代时计算得各节点注入有功功率为:
-1.6      
-2.00001  
-3.7      
5         
2.57942   

第10次迭代时,用第一因子表计算的常数项delt_P/V为:
-3.30244e-006
6.79603e-006
-3.68524e-006
-2.1819e-006

第10次迭代后，由第一因子表计算得delt_theta为：
2.43113e-007
-4.89615e-007
7.4767e-008
-4.5525e-007

------------------------------------------------------------------
第10次迭代时计算得各节点注入无功功率为:
-0.799996 
-0.999999 
-1.3      
1.81308   
2.2994    

第10次迭代时,用第二因子表计算的常数项delt_Q/V为:
-4.90376e-006
-5.03014e-007
-1.01977e-006
0

第10次迭代后，由第二因子表计算得delt_V为：
8.21955e-007
6.12167e-008
9.87503e-008
0

------------------------------------------------------------------

第10次迭代后最大功率误差为: 7.32555e-006
第10次迭代后各节点功率及电压值为:
电压值V   相角      有功P     无功Q     
0.862151  -0.0834007-1.6      -0.799996 
1.07792   0.311603  -2.00001  -0.999999 
1.03641   -0.0747337-3.7      -1.3      
1.05      0.381238  0         0         
1.05      0         0         0