复杂网络N-R法潮流分析与计算的设计

1、v1.0 可编辑可修改1电气工程及其自动化专业课程设计复杂网络 N-R 法潮流分析与计算的设计学生学号：_学生姓名：_班 级：_指导教师：_起止日期：_哈尔滨工程大学自动化学院v1.0 可编辑可修改-II -课程设计报告撰写内容一、设计要求（宋体，小四号字，加黑）用 matlab 编程，N_R 法计算潮流分布具体要求为：（1） 给出程序，并给出注释（2） 输出迭代次数，各节点电压，各支路电流（3） 在图中标明功率流向节点数据如下表所示（标幺值）123456P35Q1V10支路及变压器数据线路T1T2L2L3L4L5阻抗+v1.0 可编辑可修改-III-导纳/2变比:1:1精度要求:二、设计方案

2、（要求给出详细的设计思路及其必要的论证）（1.）潮流计算的方法（1） 高斯雅克比迭代法（2） 高斯-塞得尔法（对初值要求底，迭代次数多）（3） 牛顿-拉夫逊法（使用广泛）（4） PQ 快速分解法（提升运算速度）目前广泛应用的潮流计算方法都是基于节点电压法的，以节点导纳矩阵丫作为电力网络的数学模型。节点电压 Ui 和节点注入电流 Ii 由节点电压方程YV=I （ 1）根据 S=VI * （I *为 I 的共轭）可得非线性的节点方程YV=|=（S/V） * （2）在实际的电力系统中，已知的运行条件不是节点的注入电流，而是负荷和发电机的功率，而且这些功率一般不随节点电压的变化而变化。由于各节点注入功

3、率与注入电流的关系为 Si = Pi + jQi =ViIi *，因此可将式（2） 改写为Ii=Si/Vi=Pi+jQi/Vi （i= 1,2,3? n）（3）式中，Pi 和 Qi 分别为节点 i 向网络注入的有功功率和无功功率，当i 为发电机节点时 Pi 0;当 i 为负荷节点时 Pi 0;当 i 为无源节点 Pi =0，Qi = 0; Vi 和 Ii 分别为节点电压相量 Vi 和节点注入电流相量 Ii 的共 轭。式（3）亦即潮流计算的基本方程式，它可以在直角坐标也可以在极坐v1.0 可编辑可修改-IV -标上建立 2n 个实数形式功率方程式。发电机 Pi、Qi 为正，负荷 Pi、Qi 为负

4、。展开 YV=I 为Ii=艺 YijVj=YiiVi+ 艺 YijVi（ i=1 2 3? n）将式（4）代入式（3）,得 n 维的非线性复数的电压方程组潮流计算的基本方程为（Pi-jQi）/ Vi= YiiVi+ 艺 YijVi （i=1,2,? n）（5）（2.）变量的分类假设系统中有 n 个节点，构成 n 个复数方程，2n 个实数方程，变量总数为 6n 个。a）不可控变量（2n 个）：负荷消耗的有功功率 Li P 和无功功率 Li Q .由于该类变 量无法控制，取决于用户，而且出现事先没有预计的变动，使系统偏离原始运行状态，因此又称为不可控变量或扰动变量。b）控制变量（2n 个）：发电机

5、发出的有功功率 Gi P 和无功功率 Gi Q，因为该类 变量可控。也称独立变量。c）状态变量（2n 个）：母线电压或节点电压的幅值大小 i V 与相角大小 iS，又 称依从变量或因变量。并且 i V 受 Gi P 控制，iS受 Gi Q 控制。其中 2n 个扰动变量是给定的，2n 个控制变量和 2n 个状态变量中给定两个，求另外两个。（3.）变量的约束条件 a）扰动变量没有约束条件 b）控制变量约束条件：为满足发电机的技术经济特性指标。c）状态变量的 i V 的约束条件：保证良好的电能质量。d）状态变量的 iS的约束条件：保证系统的稳定运行v1.0 可编辑可修改-V -& 68,-

6、8丿J?Jmax(4.)系统节点的分类，根据给定的控制变量和状态变量进行分类如下:(1)PQ 节点(即负荷节点)： 给定 Pi、Qi，求 Vi 和 iS( i i e , f )。通常变电所都是这一类型的节点，由于没有发电设备，因而发电功率为零电力系统中的绝大多数节点属于这一节点。其包含变电站节点(即联络节点或浮游节点)。G e - B f y J i/ J)y G f + B wU J U J)(6) )(2)PV 节点(即调节节点、电压控制节点)：给定 Pi 和 Vi,求 Qi 和 iS( i i e , f )。这类节点必须有足够的可调无功容量,用以维持给定的电压幅值。一般时选择有一顶武

7、功储备的发电厂和具有可调无功电源设备的变电所作为 PV 节点。在电力系统中，这类节点数很少。(3)平衡节点(即松弛节点、参考节点、基准节点)：给定 Vi 和 iS( iS=0 )，求 Pi 和 Qi。(只有一个)有功功率不能给定，这个节 点承担了系统的有功功率平衡。同时其电压幅值也是给定的，相位为零。v1.0 可编辑可修改-VI -(5. )P-Q 分解法是从改进和简化牛顿法潮流程序的基础上提出来的，它的基本思想是：把节点功率表示为电压向量的极坐标方程式，抓住主要矛盾，以有功功率误差作为修正电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，把有功功率和无功功率迭代分开来进行。牛顿法潮流

8、程序的核心是求解修正方程式，当节点功率方程式采取极坐标系统时，修正方程式展开为：P = HA+ NV/ VQ = JA+ LAV /V以上方程式是从数学上推倒出来的，并没有考虑电力系统这个具体对象的特点。电力系统中有功功率主要与各节点电压向量的角度有关，无功功率则主要受各节点电压幅值的影响。大量运算经验也告诉我们，矩阵N 及 J 中各元素的数值相对是很小的，因此对牛顿法的第一步简化就是把有功功率和无功功率分开来进行迭代，即将式(4)化简为：AP = HAAQ = LAV /V(5)这样，由于我们把 2n 阶的线性方程组变成了二个 n 阶的线性方程组，对牛顿法的第二个化简，也是比较关键的一个化简

9、，即把式(5)中的系数矩阵简化为在迭代过程中不变的对称矩阵。众所周知，一般线路两端电压的相角差是不大的(通常不超过 1020 度)，因此可以认为：COS0. - 1tr B.号号弓此外，与系统各节点无功功率相应的导纳 Li B 必定远远小于该节点自导纳的虚部，即：v1.0 可编辑可修改-VII -瓦=Siiv1.0 可编辑可修改-VIII -因此, 考虑到以上关系后，式（5）中系数矩阵中的元素表达式可以化简为:HR%疔彳巧场（8）将它代入（5）中，并利用乘法结合率，我们可以把修正方程式变为:勺0Z、KPl血九5华q：B21禺2尻1（心尤0瓷17就可得到Byv1.0 可编辑可修改-10 -以上两

10、式就是 P-Q 分解法达到修正方程式，其中系数矩阵只不过是系统导纳 矩阵的虚部，因而是对称矩阵，而且在迭代过程中维持不变。它们与功率误差 方程式rAl耳久方口 %九.B吐B心-pn玖厂v1.0 可编辑可修改-XI -0= Q乜-P工VgjSU1薔-4 cos嗣(2= 123刃)构成了 P-Q 分解法迭代过程中基本计算公式，其迭代步骤大致是：(1) 给定各节点电压向量的电压初值 V i (0) ,9i (0);(2) 根据(12)计算各节点有功功率误差 Pi ,并求出； Pi/Vi(3) 解修正方程式(11)，并进而计算各节点电压向量角度的修正量 i修正各节点电压向量角度9i ;(5) 根据式(

11、16)计算各节点无功功率误差 i Q，并求出/ ; i i Q V(6) 解修正方程式(11)，求出各节点电压幅值的修正量 i V(7) 修正各节点电压幅值 i V (i ) (i 1) (i 1)i i i V =V - V - (18)(8) 返回(2)进行迭代，直到各节点功率误差及电压误差都满足收敛条件。P-Q 分解法与牛顿法潮流程序的主要差别表现在它们的修正方程式上。P-Q分解法通过对电力系统具体特点的分析，对牛顿法修正方程式的雅可比矩阵进 行了有效的简化和改进，有以下三个特点：(1) 在提高计算速度和减少内存方面的作用是明显的，不再叙述。(2) 使我们得到以下好处。首先，因为修正方程

12、式的系数矩阵就是导纳矩阵的虚部，因此在迭代过程中不必象牛顿法那样进行形成雅可比矩阵的 计算，这样不仅是仅减少了运算量，而且也大大简化了程序。其次，由于系数 矩阵在迭代过程中维持不变，因此在求解修正方程式时，可以迅速求得修正量， 从而显著提高了迭代速度。(3) 可以使我们减少形成因子表时的运算量，而且由于对称矩阵三角分解后，其上三角矩阵和下三角矩阵有非常简单的关系，所以在计算机中可以只储v1.0 可编辑可修改-XII -存上三角矩阵或下三角矩阵，从而也进一步节约了内存。三、设计内容%本程序的功能是用牛顿一一拉夫逊法进行潮流计算% B1 矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、线路对地

13、电纳（或变压器导 纳）；%5、支路的变比；6、支路首端处于 K 侧为 1，1 侧为 0；%7、线路/变压器标识（0/1）变压器参数当支路首端处于 K 侧标识为 1 时归算至末端侧,0 归算至首端侧% B2 矩阵：1、该节点发电机功率；2 、该节点负荷功率；%3、PQ 节点电压初始值，或平衡节点及 PV 节点电压的给定值%4、节点所接无功补偿并联电容（感）的电纳%5、节点分类标号：1 为平衡节点（应为 1 号节点）；2 为 PQ 节点；3 为 PV 节点；clear;isb=1; %i nput（请输入平衡母线节点号：isb=）;pr=; %input（请输入误差精度：pr=）;%-n=6;%i

14、nput（请输入节点数：n=）;nl=6;%input（请输入支路数：nl=）;B1=1 20+ 011:2 3+ 0+ 10 0;2 5+ 0+ 10 0;3 4+ 0+ 10 0;4 5+01 0 0;6 50+01 1v1.0 可编辑可修改-XIII -B2=0010 1;0 3+1i0 2;0 + 0 2;0 +0 2;0 +0 2;0 -5+0i0 3%input(请输入各节点参数形成的矩阵：B2=);%X=1 0;2 0;3 0;4 0;5 0;6 0%-%n=4;%input(请输入节点数：n=);nl=4;%input(请输入支路数：nl=);%B1=1 2 4+16i 0 1

15、 0 0;1 3 4+16i 0 1 0 0;2 3 2+8i 0 1 0 0;2 4 + 0 11/110 0 1 %i nput(请输入由支路参数形成的矩阵：B 仁)；%B2=0 0 115 0 1;0 0 110 0 2;0 20+4i 110 0 2;0 10+6i 10 0 2 %i nput(请输入各节点参数形成的矩阵：B2=);%-Y=zeros( n);e=zeros(1, n);f=zeros(1, n);V=zeros(1, n);sida=zeros(1, n);S1=zeros( nl);% % %- 求导纳矩阵-%for i=1: n% if X(i,2)=0;% p

16、=X(i,1);% 丫(p,p)=1/X(i,2);%e nd%e ndfor i=1:nl%从 1 到 n1 (总支路数)if B1(i,7)=1%-如果是变压器支路-if B1(i,6)=0加节点(首端)处于 1 侧v1.0 可编辑可修改-XIV -p=B1(i,1);q=B1(i,2);v1.0 可编辑可修改-XV -end disp(导纳矩阵 Y=);disp(Y);给定各节点初始电压及给定各节点注入功率e(i)=real(B2(i,3);f(i)=imag(B2(i,3);end% 合定各节点注入功率else左节点（首端）处于 K侧endp=B1(i,2);q=B1(i,1);丫(p

17、,q)=Y(p,q)-1./(B1(i,3)*B1(i,5); %非对角元Y(q,p)=Y(p,q);非对角元Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)A2);%对角元K侧丫(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4);%寸角元1侧+励磁导纳else %否则为线路支路p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1./B1(i,3);%非对角元Y(q,p)=Y(p,q);非对角元丫(q,q)=Y(q,q)+1./B1(i,3)+B1(i,4)./; %对角元 j 侧+线路电纳的一半end丫(p,p)=Y(p,p)+1./B1(i,

18、3)+B1(i,4)./; %对角元 i 侧+线路电纳的一半S(i)=B2(i,1)-B2(i,2);%i 节点注入功率 SG-SLB(i,i)=B(i,i)+B2(i,4);%i 节点无功补偿量（电纳值）G=real(Y);B=imag(Y);%分解出导纳阵的实部和虚部for i=1:n%合定各节点初始电压的实部和虚部V(i)=abs(B2(i,3);%PV 平衡节点及 PQ 节点电压模值for i=1:nv1.0 可编辑可修改-16 -end=牛顿-拉夫逊法迭代求解非线性代数方程(功率方程)P=real(S);Q=imag(S);%分解出各节点注入的有功和无功功率ICT1=0;IT2=1;

19、N0=2*n;N 仁 N0+1;a=0; %迭代次数 ICT1、a；不满足收敛要求的节点数 IT2 whileIT2=0% N0=2*n雅可比矩阵的阶数；N 仁 N0+扩展列IT2=0;a=a+1;JZ=Jacobi 矩阵第(,num2str(a),)次消去运算;JZ 仁Jacobi矩阵第(,num2str(a),)次回代运算;JZ0=功率方程第(,num2str(a),)次差值：;%- 求取各个节点的功率及功率偏差及 PV 节点的电压偏差for i=1:n%n个节点 2n 行(每节点两个方程 P 和 Q 或 U)p=2*i-1;m=p+1;C(i)=0;D(i)=0;for j1=1:n%第

20、 i 行共 n 列(n 个节点间互导纳及节点电压相乘即电流)C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);% 工(Gij*ej -Bij*fj)D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%艺(Gij*fj+Bij*ej)end% 求 i 节点有功和无功功率 P,Q的计算值P1=C(i)*e(i)+f(i)*D(i);%节点功率 P 计算 ei 工(Gij*ej-Bij*fj)+fi工(Gij*fj+Bij*ej)Q 仁 C(i)*f(i)-e(i)*D(i);%节点功率 Q 计算 fi 工(Gij*ej-Bij*fj)-ei工(G

21、ij*fj+Bij*ej)V2=e(i)A2+f(i)A2;灿压模平方%=求取功率差及 PV 节点电压模平方差if i=isb%非平衡节点(PQ 或 PV 节点)if B2(i,5)=3%非 PV 节点(只能是 PQ 节点)J(m,N1)=P(i)-P1;%PQ点有功功率差 J(m,N1)扩展列 PJ(p,N1)=Q(i)-Q1;%PQ节点无功功率差 J(p,N1)扩展列 Qv1.0 可编辑可修改-XVII -else%PV节点-J(m,N1)=P(i)-P1;%PV 节点有功功率差 J(m,N1)扩展列 PJ(p,N1)=V(i)A2-V2;%P节点电压模平方差 J(p,N1)扩展列 Uen

22、dend %(if i=isb)非平衡节点(PQ 或 PV 节点)end %(for i=1: n) n个节点 2n 行(每节点两个方程 P 和 Q 或 U)for m=1:N0JJN1(m)=J(m,N1);enddisp(JZ0);disp(JJN1);%- 判断功率偏差量及 PV 节点的电压偏差量是否满足要求for k=3:N0%除去平衡节点 1、2 号以外的所有节点DET=abs (J( k,N1);if DET=pr; %PQ节点的功率偏差量及 PV 节点的电压偏差量是否满足要求IT2=IT2+1; %不满足要求的节点数加 1endendICT2(a)=IT2;%不满足要求的节点数；

23、a 为迭代次数ICT1=ICT1+1;%迭代次数if ICT2(a)=0;%当前不满足要求的节点数为零break %退出迭代运算end%-以上为求取各个节点的功率及功率偏差及PV 节点的电压偏差%=求取 Jacobi 矩阵形成修正方程=个节点 2n 行(每节点两个方程 P 和 Q 或 U)for i=2:n%nv1.0 可编辑可修改-18 -j1=2:n% 第 i 行共 n 列(2n 个 Jacobi 矩阵元素 dP/de 及 dP/df 或 dQ/deX1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/deX2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i

24、);% dP/dfX3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/deX4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/dfp=2*i-1;q=2*j1-1;J(p,q)=X3;%if i=isb%非平衡节点（PQ 或 PV 节点）if B2(i,5)=3%F 面是针对 PQ 节点来求取 Jacobi 矩阵的元素C(i)=0;D(i)=0;for j1=1: n %第 i 行共 n 列（n 个节点间互导纳及节点电压相乘即电流）C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);% 艺(Gij*ej -Bij*fj)D

25、(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%艺(Gij*fj+Bij*ej)endfor及 dQ/df)ifj1=is%非平衡节点&非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i);% X仁dP/de=-dQ/df=-X4X2=B(i,j1)*e(i)-G(i,j1)*f(i);% X2=dP/df=dQ/de=X3X3=X2;% X2=dp/df X3=dQ/deX4=-X1;% X 仁 dP/de X4=dQ/dfp=2*i-1;q=2*j1-1;J(p,q)=X3;m=p+1;% X3=dQ/de J(p,N)=DQ 节点无功功率

26、差J(p,N)=DQ;J(m,q)=X1;q=q+1;% X 仁 dP/de J（m,N）=DP 节点有功功 率差J(m,N)=DP;J(p,q)=X4;J(m,q)=X2;% X4=dQ/df X2=dp/dfelseif j1=i&j1=isb%非平衡节点&对角元扩展列 Q J(p,N)=DQ;v1.0 可编辑可修改-XIX -m=p+1;J(m,q)=X1;q=q+1;J(p,q)=X4;%扩展列 P J(m,N)=DP;J(m,q)=X2;endendelse %if B2(i,5)=3 %否则(即为 PV 节点)%=下面是针对 PV 节点来求取 Jacobi 矩阵的元

27、素for j1=1: nif j1=isb&j1=i%非平衡节点&非对角元X仁-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/deX2=B(i,j1)*e(i)-G(i,j1)*f(i);% dP/dfX5=0;X6=0;p=2*i-1;q=2*j1-1;J(p,q)=X5; % PV节点电压误差 J(p,N)=DV;m=p+1;J(m,q)=X1;q=q+1;J(p,q)=X6; % PV节点有功误差 J(m,N)=DP;J(m,q)=X2;elseif j1=i&j1=isb%非平衡节点&对角元X仁-C(i)-G(i,i)*e(i)-B(i

28、,i)*f(i);% dP/deX2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/dfX5=-2*e(i);X6=-2*f(i);p=2*i-1;q=2*j1-1;J(p,q)=X5; % PV节点电压误差 J(p,N)=DV;m=p+1;J(m,q)=X1;q=q+1;J(p,q)=X6; % PV节点有功误差 J(m,N)=DP;J(m,q)=X2;v1.0 可编辑可修改-XX -endv1.0 可编辑可修改-XXI -endend %(if B2(i,5)=3 else)end %(if i=isb)end %(for i=1:n)n个节点 2n 行(每节点两个

29、方程 P 和 Q 或 U)JZO=形成的第(,num2str(a),)次 Jacobi 矩阵：;disp(JZO);disp(J);%=以 上 为形成 完整的 Jacobi 矩 阵%=下面用高斯消去法对由 Jacobi 矩阵形成的修正方程进行求解(按列消去、回代)按列消去运算=for k2=k+1:N0 %从 k+1 行到 2*n 最后一行for k3=k+1:N1%从 k2+1 列到扩展列消去 k+1 行后各行下三角元素J(k2,k3)=J(k2,k3)-J(k2,k)*J(k,k3);%消去运算end %用当前行 K3 列元素减去当前行 K 列元素乘以第 k 行 K3 列元素for k1=

30、k+1:N1%从 k+1 列的 Jacobi 元素到扩展列的厶 P、AQ 或AUJ(k, k1)=J(k,k1)./J (k, k);%行规格化用 K 行 K 列对角元素去除 K 行 K 列后的非对角元素进for k=3:N0% N0=2*n(从第三行开始，第一、二行是平衡节点)endJ(k,k)=1;%对角元规格化 K 行 K 列对角元素赋 1%=v1.0 可编辑可修改-22 -J(k2,k)=0; %当前行第 k 列元素已消为 0endendJZ=Jacobi 矩阵第(,num2str(a),)次消去运算;JZ1=Jacobi 矩阵第(,nu m2str(a),)次回代运算;v1.0 可编

31、辑可修改-XXIII -disp(JZ);disp(J);%=按列回代运算=for k=N0:-1:3for k仁k-1:-1:3J(k1,N1)=J(k1,N1)-J(k1,k)*J (k, N1);J(k1,k)=0;endendfor m=1:N0JJN1(m)=J(m,N1);enddisp(JZ1);disp(JJN1);%disp(J);%-for k=3:2:N0-1修改节点电压-L=(k+1)./2;e(L)=e(L)-J(k,N1);%k1=k+1;f(L)=f(L)-J(k1,N1);%U(L)=sqrt(e(L)A2+f(L)A2);end修改节点电压实部修改节点电压虚部

32、disp(各个节点电压模);disp(U);=结束一次迭代=endQ%*disp(迭代次数：)；disp(ICTI-l);F 面为迭代计算结束后的有关输出过程*v1.0 可编辑可修改-24 -v1.0 可编辑可修改-25 -disp(没有达到精度要求的个数：)；disp(ICT2);for k=1: nV(k)=sqrt(e(k)A2+f(k)A2); %计算各节点电压的模值sida(k)=ata n(f(k)./e(k)*180./pi; %计算各节点电压的角度E(k)=e(k)+f(k)*j;%将各节点电压用复数表示end%=计算各输出量=disp(各节点的电压复数值 E 为(节点号从小到

33、大排列)：)；disp(E); % 显示各节点的实际电压值 E 用复数表示disp(- );disp(各节点的电压模值大小 V 为(节点号从小到大排列)：);disp(V); % 显示各节点的电压大小 V 的模值disp(- );disp(各节点的电压相角 sida 为(节点号从小到大排列)：)；disp(sida); %显示各节点的电压相角for p=1: nC(p)=0;for q=1: nC(p)=C(p)+conj(Y(p,q)*conj(E(q); %计算各节点注入电流的共轭值endS(p)=E(p)*C(p);%计算各节点的功率 S =电压 X 注入电流的共轭值enddisp(各节

34、点的功率 S 为(节点号从小到大排列)：)；disp(S); %显示各节点的注入功率v1.0 可编辑可修改-26 -disp(- );disp(各条支路的首端功率 Si 为(顺序同您输入 B1 时一致)：)；for i=1: nlp=B1(i,1);q=B1(i,2);if B1(i,7)=0Si(p,q)=E(p)*conj(E(p)*B1(i,4)./2+(E(p)-E(q)./B1(i,3);Siz(i)=Si(p,q);elseif B1(i,6)=0Si(p,q)=E(p)*(conj(E(p)*B1(i,4) +(E(p)*B1(i,5)-E(q)*(1./(B1(i,3)*B1(

35、i,5);Siz(i)=Si(p,q);elseSi(p,q)=E(p)*conj(E(p)-E(q)*B1(i,5)*(1./(B1(i,3)*B1(i,5)A2);Siz(i)=Si(p,q);endendZF=S(, nu m2str(p), ,n um2str(q),)= ,n um2str(Si(p,q);disp(ZF);disp(- );enddisp(各条支路的末端功率 Sj 为(顺序同您输入 B1 时一致)：)；for i=1: nlp=B1(i,1);q=B1(i,2);if B1(i,7)=0Sj(q,p)=E(q)*conj(E(q)*B1(i,4)./2+(E(q)-

36、E(p)./B1(i,3);Sjy(i)=Sj(q,p);v1.0 可编辑可修改-xxi -elseif B1(i,6)=0Sj(q,p)=E(q)*conj(E(q)-E(p)*B1(i,5)*(1./(B1(i,3)*B1(i,5F2);Sjy(i)=Sj(q,p);elseSj(q,p)=E(q)*(conj(E(q)*B1(i,4) +(E(q)*B1(i,5)-E(p)*(1./(B1(i,3)*B1(i,5);Sjy(i)=Sj(q,p);endendZF=S(, nu m2str(q), nu m2str(p),)= ,n um2str(Sj(q,p);disp(ZF);disp

37、(- );enddisp(各条支路的功率损耗 DS 为(顺序同您输入 B1 时一致)：)；for i=1: nlp=B1(i,1);q=B1(i,2);DS(i)=Si(p,q)+Sj(q,p);ZF=DS(, nu m2str(p), nu m2str(q),)=, num2str(DS(i);disp(ZF);disp(- );endzws=0;JDDY=;JDP=;JDQ=;JDDYJD=;for i=1:n%总网损为所有节点注入功率的代数和zws=zws+S(i);JDDYJD=strcat(JDDYJD, num2str(i),(, nu m2str(sida(i),),);v1.0

38、 可编辑可修改-XXVIII -JDDY=strcat(JDDY, nu m2str(i),(, nu m2str(V(i),),);JDP=strcat(JDP, nu m2str(i),( ,n um2str(real(S(i),),);JDQ=strcat(JDQ ,n um2str(i),(, nu m2str(imag(S(i),),); endJDDYJD=strcat(节点电压角度：JDDYJD);JDDY=strcat(节点电压模值：,JDDY);JDP=strcat(节点有功:,JDP);JDQ=strcat(节点无功：JDQ);ZFsubplot(4,1,2);plot(s

39、ida);title(JDDYJD);ylabel(节点电压角度);grid on;subplot(4,1,3);P=real(S);Q=imag(S);bar(P);title(JDP);ylabel(节点注入有功);grid on;subplot(4,1,4);bar(Q);title(JDQ);xlabel(ZF);ylabel(节点注入无功);grid on*figure(2);v1.0 可编辑可修改-XXIX -subplot(3,2,1);JDH=;JDH 仁;for i=1: nlJDH=strcat(JDH, num2str(i),(, num2str(B1(i,1), nu

40、m2str(B1(i,2),),);JDH 仁 strcat(JDH1, num2str(i),(, nu m2str(B1(i,2), ,n um2str(B1(i,1),),); endP 仁 real(Siz);Q 1= imag(Siz);bar(P1);title(JDH);ylabel(支路首端注入有功);%xlabel(支路号)；grid on;subplot(3,2,2);bar(Q1);title(JDH);ylabel(支路首端注入无功);grid on;subplot(3,2,3);P2=real(Sjy);Q2=imag(Sjy);bar(P2);title(JDH1)

41、;ylabel(支路末端注入有功);grid on;subplot(3,2,4);bar(Q2);title(JDH1);ylabel(支路末端注入无功);v1.0 可编辑可修改-XXX -grid on;subplot(3,2,5);P3=real(DS);Q3=imag(DS);bar(P3);xlabel(JDH);ylabel(支路有功损耗);grid on;subplot(3,2,6);bar(Q3);xlabel(JDH);ylabel(支路无功损耗);grid on;四、结果分析及结论B1 =Columns 1 through 50 +0+0 +0 +0 +00 +0Column

42、s 6 through 700000000v1.0 可编辑可修改Colum ns 11 through 12-XXXI -B2 =0 0 00+00+00+00+00 0导纳矩阵丫=Columns 1 through 500 +0000 +0+0+- +000+ - +0+0+00000 +Colum n 600000 +0功率方程第(1)次差值：Colum ns 1 through 100 0v1.0 可编辑可修改-XXXII -形成的第次 Jacobi 矩阵:0 0 0 0 00 0 0 0 00 00 00 00 0000 00 00 0 0 0 0 0Colum ns 11 throu

43、gh 130000000000Colum ns 1 through 100000000000000000000000000000v1.0 可编辑可修改Colum ns 11 through 13-XXXIII -0 00 0000 0Jacobi 矩阵第(1)次消去运算Colum ns 1 through 100 0 0 0 00 0 0 0 00 00 000 00 0 00 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 000000000000000000000000000000000000000000000v1.0 可编辑可修改Colum ns

44、 1 through 10-XXXIV -0 0 00 00 00 00Jacobi 矩阵第次回代运算Colum ns 1 through 100 0Colum ns 11 through 120各个节点电压模0功率方程第(2)次差值：Colum ns 1 through 100 0Colum ns 11 through 12形成的第次 Jacobi 矩阵:v1.0 可编辑可修改-XXXV -0 00 00 00 00 0 0 0 0 0 0 0 0 0 0 0 0 0Colum ns 11 through 130000000 00 00 00 0v1.0 可编辑可修改-36 -Jacobi

45、矩阵第(2)次消去运算Colum ns 1 through 100 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 00 0Colum ns 11 through 1300000000000 0v1.0 可编辑可修改-XXXVII -0 00 00 0Jacobi 矩阵第次回代运算Colum ns 1 through 100 0Colum ns 11 through 12各个节点电压模0功率方程第(3)次差值：Colum ns 1 through 10 0 0Colum ns 11 through 12形成的第次 Jacobi 矩阵:Colum ns 1 through 100

46、0 0 00 00 0 0 000000000v1.0 可编辑可修改-XXXVIII -0 0 0 00 0 0 0000000000 0 0 00 0 0 0000 0 0 0 0 00 0000 0 0 0 0 0Colum ns 11 through 13000000000000000000Jacobi 矩阵第(3)次消去运算Colum ns 1 through 100 0 0 0 0 0 0 00 0v1.0 可编辑可修改-XXXIX -0 0 0 0 0 0 0 00 0Colum ns 11 through 130000000000v1.0 可编辑可修改-40 -Jacobi 矩阵第次回代运算Colum ns 1 through 10 0 0Colum ns 11 through 12各个节点电压模0功率方程第(4)次差值：Colum ns 1 through 100 0Colum ns 11 through 12形成的第次 Jacobi 矩阵:Colum ns 1 through 10v1.0 可编辑可修改-XLI-0000000000000000000000000000000