—()=U 可得
=1.039528； =0.008050； =0.006386； =0.000072； 所以
=1.039528； =0.008050； =0.006386； =0.000072；
求得各节点电压的新值后，就开始第二次迭代。每次迭代所得适于表1~表3。由表1可见经6次迭代就可满足ε≤的要求。 表1 迭代过程中各节点功率的不平衡量
K 0 -0.500000 1 -0.411720 -0.136497 2 0.043223 0.014875 -0.005899 3 -0.004833 4 0.000537 -0.001658 0.000184 -0.000650 -0.000072 5 -0.000060 6 0.000007 -0.000021 0.000000 -0.000001 0.000000 -0.000001 -0.000002 0.000000 0.000008 0.000003 -0.000096 0.000011 -0.000033 0.000004 -0.000150 0.000016 -0.000051 0.000006 -0.001936 0.000215 -0.007856 0.000872 -0.002664 0.000295 -0.012588 0.001376 -0.004265 0.000466 1.211930 -0.375000 0.049802 -0.077279 0.018300 0.076248 0.400000 -0.191947 0.024849 -0.600000 0.120899 -0.319599 0.040218
表2 迭代过程中各节点电压的修正量
k 0 1 2 3 4 5 6 -0.036952 0.039528 -0.085524 0.008050 0.000757 -0.091810 0.006386 0.002461 0.000941 -0.108571 0.000072 0.004597 0.001618 -0.010063 -0.003406 0.001828 0.001054 0.000385 -0.000230 -0.000064 -0.000285 -0.000082 -0.000484 -0.000149 0.000006 -0.000028 0.000008 0.000049 0.000015 -0.000120 -0.000043 0.000023 0.000013 0.000005 -0.000002 -0.000001 -0.000003 -0.000001 -0.000005 -0.000002 0.000000 0.000000 0.000000 0.000001 0.000000 -0.000002 -0.000001 0.000000
表3 迭代过程中各节点电压
k 0 1 2 3 4 5 6 0.000000 1.000000 1.000072 1.001690 1.000000 0.000000 1.000000 0.000000 1.000000 0.000000 -0.036952 1.039528 -0.085524 1.008050 -0.091810 1.006386 -0.047016 1.036122 -0.083696 1.008808 1.007327 -0.045962 1.036507 -0.083926 1.008744 -0.089633 1.007245 -0.104458 1.001540 -0.046082 1.036463 -0.083903 1.008750 -0.089605 1.007253 -0.104410 1.001555 -0.046068 1.036468 -0.083906 1.008750 -0.089608 1.007252 -0.104415 1.001554 -0.046070 1.036468 1.008750 -0.089608 1.007252 -0.104414 1.001554
=1.036468； =1.008750 =1.007252； =1.001554 6.计算平衡节点功率和线路功率
迭代收敛后，就可以计算平衡节点功率和线路功率。结果如下
=1.298162+j0.244472 各线路功率如表4
表4 各线路功率 i 1 0.889505 0.408657 j 1 2 3 4 5 +j0.138662 +j0.105810 2 -0.875079 -j0.095385 3 -0.395969 -0.243108 -j0.067747 -j0.070133 4 -0.274598 -0.188726 0.246884 0.279319 0.548870 +j0.081457 +j0.080606 +j0.13321 0.189079 -j0.012120 0.063325 +j0.003267 -0.063008 -j0.002316 -j0.066445 +j0.013178 5 -0.536990 -j0.097684 计算完毕。
第三章 程序设计
3.1 设计流程
3．2 程序设计
针对两机五节点网络潮流计算基于MATLAB的程序如下：
Clc %清除当前屏幕上显示的内容，但并不清楚工作空间的数据 Clear %清除工作空间内的所有变量 disp('节点总数为：'); %显示内容 N=5
disp('平衡节点为：'); 1
Y=input('请输入导纳矩阵：Y=');
pr=input('请输入误差精度：pr=');
G=real(Y);B=imag(Y); %real实部，image虚部 B1=B(2:5, 2:5); %去掉平衡节点
B2=B1; b1=inv(B1); %取逆 b2=inv(B2);
ure=input('请输入平衡节点1电压初值的实部，Re='); dim=input('请输入平衡节点1电压初值的虚部，Im='); for m=2:5
dt(m)=0; %电压虚部初值 u(m)=1.0; %电压实部初值 end for m=2:5
p(m)=input('请输入PQ节点的有功功率：p='); q(m)=input('请输入PQ节点的无功功率：q='); end
k=0;wucha=1; %迭代次数， while wucha(1)>pr; u(1)=ure;
dt(1)=dim; %平衡节点电压实部 for m=2:5 %平衡节点电压虚部 for n=1:5
pt(n)=u(m)*u(n)*(G(m,n)*cos(dt(m)-dt(n))+B(m,n)*sin(dt(m)-dt(n))); % end