bdzjf 发表于 2009-11-14 11:38:18

夏道止老师《电力系统分析》N-R迭代过程

=========================导纳矩阵输出开始:=========================
2.9049 - j11.5030 0.0000 + j5.3182 -1.6600 + j3.1618 -1.2450 + j2.3714
0.0000 + j5.3182 0.0000 - j4.6638 0.0000 + j0.0000 0.0000 + j0.0000
-1.6600 + j3.1618 0.0000 + j0.0000 2.4900 - j4.7039 -0.8300 + j1.5809
-1.2450 + j2.3714 0.0000 + j0.0000 -0.8300 + j1.5809 2.0750 - j3.9092
#########################导纳矩阵输出完毕!#########################

=========================初始功率向量输出开始:=========================
0.0000 + j0.0000 -0.5000 + j-0.3000 0.2000 + j0.0000 0.0000 + j0.0000
#########################初始功率向量输出完毕!#########################

=========================初始电压向量输出开始:=========================
1.0000 + j0.0000 1.0000 + j0.0000 1.0500 + j0.0000 1.0500 + j0.0000
#########################初始电压向量输出完毕!#########################

===================Iterationk = 0 Begin=========================
P = -0.14525
Q = 0.37497
P = 0.00000
Q = -0.65434
P = 0.08715
Q = 1.10250
误差向量 p、q:
0.14524743, 0.14524743, -0.37496691, -0.50000000, 0.35433884, 0.11285154, 0.00000000,
电流 aii、bii:
aii = -0.1452474277   bii = -0.3749669147 aii = 0.0000000000   bii = 0.6543388430 aii = 0.0829985301   bii = 0.0000000000
========================雅克比矩阵========================
-2.7597000    -11.128100    0.00000000    5.31820000    1.66000000    3.16180000   
-11.878000    3.05020000    5.31820000    0.00000000    3.16180000    -1.6600000   
0.00000000    5.31820000    0.00000000    -5.3182000    0.00000000    0.00000000   
5.31820000    0.00000000    -4.0095000    0.00000000    0.00000000    0.00000000   
1.74300000    3.31990000    0.00000000    0.00000000    -2.6975000    -4.9391000   
0.00000000    0.00000000    0.00000000    0.00000000    -2.1000000    0.00000000   

所求矩阵的逆矩阵为:
-0.08369   -0.15942   -0.08369   -0.21145   0.00000   -0.30618   
-0.23406   0.07113   -0.23406   0.09434   -0.17374   0.14524   
-0.11101   -0.21145   -0.11101   -0.52987   0.00000   -0.40612   
-0.23406   0.07113   -0.42209   0.09434   -0.17374   0.14524   
0.00000   -0.00000   0.00000   0.00000   0.00000   -0.47619   
-0.18686   -0.00845   -0.18686   -0.01121   -0.31925   0.24965   
方程的解:
0.0145420.070185-0.0690860.1642020.0000000.029460
新的电压值(e+jf):
0.98546 - j0.07019, 1.06909 - j0.16420, 1.05000 - j0.02946,
===================Iterationk = 0 End=========================
===================Iterationk = 1 Begin=========================
P = -0.02835
Q = -0.00268
P = -0.46151
Q = -0.20794
P = 0.20354
Q = 1.10337
误差向量 p、q:
0.02834863, 0.02834863, 0.00267722, -0.03848582, -0.09206171, -0.00353918, -0.00086787,
电流 aii、bii:
aii = -0.0284292626   bii = 0.0047414768 aii = -0.3925561267   bii = 0.2547940506 aii = 0.1642347318   bii = 0.0000000000
========================雅克比矩阵========================
-3.6416000    -11.136600    0.37330000    5.24080000    1.85770000    2.99940000   
-11.127100    3.69850000    5.24080000    -0.3733000    2.99940000    -1.8577000   
0.87330000    5.68560000    -0.3733000    -5.2408000    0.00000000    0.00000000   
5.68560000    -0.8733000    -4.7313000    1.15840000    0.00000000    0.00000000   
1.83610000    3.27100000    0.00000000    0.00000000    -2.9173000    -4.8658000   
0.00000000    0.00000000    0.00000000    0.00000000    -2.1000000    0.05890000   

所求矩阵的逆矩阵为:
-0.11168   -0.15094   -0.13743   -0.16516   -0.01477   -0.29385   
-0.23803   0.07191   -0.22576   0.07869   -0.17259   0.13189   
-0.15534   -0.17860   -0.23173   -0.40317   -0.03179   -0.34835   
-0.26578   0.06559   -0.44212   0.08656   -0.18743   0.11893   
-0.00558   -0.00024   -0.00562   -0.00026   -0.00903   -0.46893   
-0.19882   -0.00847   -0.20026   -0.00927   -0.32170   0.25892   
方程的解:
0.017232-0.0046150.0415680.0022480.0005200.003815
新的电压值(e+jf):
0.96823 - j0.06557, 1.02752 - j0.16645, 1.04948 - j0.03328,
===================Iterationk = 1 End=========================
===================Iterationk = 2 Begin=========================
P = -0.00050
Q = 0.00005
P = -0.49877
Q = -0.29567
P = 0.20071
Q = 1.10251
误差向量 p、q:
0.00050286, 0.00050286, -0.00004686, -0.00122615, -0.00432794, -0.00071146, -0.00001483,
电流 aii、bii:
aii = -0.0005202507   bii = -0.0000131695 aii = -0.4275816134   bii = 0.3570182626 aii = 0.1577813702   bii = 0.0000000000
========================雅克比矩阵========================
-3.5664000    -10.947100    0.34870000    5.14920000    1.81460000    2.95250000   
-10.947100    3.56740000    5.14920000    -0.3487000    2.95250000    -1.8146000   
0.88520000    5.46450000    -0.3487000    -5.1492000    0.00000000    0.00000000   
5.46450000    -0.8852000    -4.4352000    1.20390000    0.00000000    0.00000000   
1.84730000    3.26310000    0.00000000    0.00000000    -2.9275000    -4.8538000   
0.00000000    0.00000000    0.00000000    0.00000000    -2.0990000    0.06660000   

所求矩阵的逆矩阵为:
-0.11450   -0.15942   -0.14639   -0.18258   -0.01421   -0.30342   
-0.23498   0.07099   -0.22078   0.08130   -0.16768   0.13058   
-0.16419   -0.19402   -0.25476   -0.44360   -0.03241   -0.36965   
-0.25794   0.06106   -0.43642   0.08493   -0.17820   0.11145   
-0.00627   -0.00040   -0.00635   -0.00046   -0.01009   -0.46835   
-0.19777   -0.01271   -0.20031   -0.01456   -0.31808   0.25478   
方程的解:
0.000934-0.0000850.0021870.0001600.0000210.000432
新的电压值(e+jf):
0.96729 - j0.06549, 1.02533 - j0.16661, 1.04946 - j0.03371,
===================Iterationk = 2 End=========================
===================Iterationk = 3 Begin=========================
P = -0.00000
Q = 0.00000
P = -0.50000
Q = -0.29999
P = 0.20008
Q = 1.10250
误差向量 p、q:
0.00000048, 0.00000048, -0.00000006, -0.00000179, -0.00001164, -0.00008305, -0.00000019,
电流 aii、bii:
aii = -0.0000005025   bii = -0.0000000245 aii = -0.4287817842   bii = 0.3622513393 aii = 0.1567496678   bii = 0.0000000000
========================雅克比矩阵========================
-3.5632000    -10.936600    0.34830000    5.14420000    1.81270000    2.94970000   
-10.936600    3.56320000    5.14420000    -0.3483000    2.94970000    -1.8127000   
0.88610000    5.45290000    -0.3483000    -5.1442000    0.00000000    0.00000000   
5.45290000    -0.8861000    -4.4197000    1.20580000    0.00000000    0.00000000   
1.84860000    3.26230000    0.00000000    0.00000000    -2.9284000    -4.8527000   
0.00000000    0.00000000    0.00000000    0.00000000    -2.0989000    0.06740000   

所求矩阵的逆矩阵为:
-0.11472   -0.15989   -0.14693   -0.18356   -0.01423   -0.30393   
-0.23481   0.07095   -0.22052   0.08145   -0.16742   0.13050   
-0.16472   -0.19489   -0.25605   -0.44590   -0.03248   -0.37084   
-0.25751   0.06086   -0.43612   0.08491   -0.17772   0.11109   
-0.00635   -0.00042   -0.00643   -0.00048   -0.01021   -0.46826   
-0.19773   -0.01296   -0.20034   -0.01488   -0.31788   0.25453   
方程的解:
0.0000040.0000130.0000080.0000140.0000010.000027
新的电压值(e+jf):
0.96729 - j0.06550, 1.02532 - j0.16662, 1.04946 - j0.03373,
===================Iterationk = 3 End=========================
===================Iterationk = 4 Begin=========================
P = -0.00000
Q = -0.00000
P = -0.50000
Q = -0.30000
P = 0.20001
Q = 1.10250
误差向量 p、q:
0.00000000, 0.00000000, 0.00000000, 0.00000000, -0.00000000, -0.00000516, -0.00000000,
电流 aii、bii:
aii = -0.0000000000   bii = 0.0000000000 aii = -0.4287787078   bii = 0.3622711617 aii = 0.1566485901   bii = 0.0000000000
========================雅克比矩阵========================
-3.5634000    -10.936500    0.34830000    5.14420000    1.81280000    2.94970000   
-10.936500    3.56340000    5.14420000    -0.3483000    2.94970000    -1.8128000   
0.88610000    5.45290000    -0.3483000    -5.1442000    0.00000000    0.00000000   
5.45290000    -0.8861000    -4.4197000    1.20590000    0.00000000    0.00000000   
1.84870000    3.26220000    0.00000000    0.00000000    -2.9284000    -4.8526000   
0.00000000    0.00000000    0.00000000    0.00000000    -2.0989000    0.06750000   

所求矩阵的逆矩阵为:
-0.11473   -0.15989   -0.14693   -0.18357   -0.01423   -0.30393   
-0.23481   0.07095   -0.22052   0.08145   -0.16742   0.13050   
-0.16473   -0.19489   -0.25606   -0.44590   -0.03248   -0.37084   
-0.25751   0.06086   -0.43611   0.08491   -0.17772   0.11108   
-0.00636   -0.00042   -0.00644   -0.00048   -0.01022   -0.46825   
-0.19773   -0.01297   -0.20034   -0.01489   -0.31788   0.25452   
方程的解:
0.0000000.0000010.0000000.0000010.0000000.000002
新的电压值(e+jf):
0.96729 - j0.06550, 1.02532 - j0.16663, 1.04946 - j0.03374,
===================Iterationk = 4 End=========================
===================Iterationk = 5 Begin=========================
P = 0.00000
Q = -0.00000
P = -0.50000
Q = -0.30000
P = 0.20000
Q = 1.10250
误差向量 p、q:
-0.00000000, -0.00000000, 0.00000000, 0.00000000, 0.00000000, -0.00000032, -0.00000000,
=================k=5时,结果收敛=======================
================牛顿法计算共耗时:0秒:609毫秒=========

seanchio 发表于 2009-11-14 12:16:05

很详细呵呵
页: [1]
查看完整版本: 夏道止老师《电力系统分析》N-R迭代过程

招聘斑竹