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

实用文案电气工程及其自动化专业课程设计复杂网络N-R法潮流分析与计算的设计学生学号：学生姓名：标准实用文案班 级：指导教师：起止日期：哈尔滨工程大学自动化学院标准实用文案课程设计报告撰写内容一、设计要求（宋体，小四号字，加黑）用matlab编程，N_R法计算潮流分布具体要求为：（1）给出程序，并给出注释（2）输出迭代次数，各节点电压，各支路电流（3）在图中标明功率流向34L4L2L52 L3 5T1T216节点数据如下表所示（标幺值）123456P31.80.63.55Q10.50.81.3V11.050支路及变压器数据标准实用文案线路T1T2L2L3L4L5阻抗j0.04j0.020.06+j0.0250.01+j0.20.06+j0.50.05+j0.3导纳/2j0.25j0.25j0.25j0.25变比1.05:11.05:1精度要求：0.0001二、设计方案（要求给出详细的设计思路及其必要的论证）1.）潮流计算的方法1）高斯雅克比迭代法2）高斯-塞得尔法（对初值要求底，迭代次数多）3）牛顿-拉夫逊法（使用广泛）4）PQ快速分解法（提升运算速度）目前广泛应用的潮流计算方法都是基于节点电压法的，以节点导纳矩阵Y作为电力网络的数学模型。节点电压Ui和节点注入电流Ii由节点电压方程YV=I（1）根据S=VI﹡(I﹡为I的共轭)可得非线性的节点方程YV=I=(S/V)﹡（2）在实际的电力系统中，已知的运行条件不是节点的注入电流，而是负荷和发电机的功率，而且这些功率一般不随节点电压的变化而变化。由于各节点注入功率与注入电流的关系为 Si＝Pi＋jQi=ViIi ﹡，因此可将式（2）改写为Ii=Si/Vi=Pi+jQi/Vi(i=1,2,3 ? n) （3）式中，Pi和Qi分别为节点i向网络注入的有功功率和无功功率，当为发电机节点时Pi﹥0；当i为负荷节点时Pi﹤0；当i为无源节点Pi＝0，Qi＝0；Vi和Ii分别为节点电压相量Vi和节点注入电流相量Ii的共轭。式（3）亦即潮流计算的基本方程式 ,它可以在直角坐标也可以在极坐标准实用文案标上建立2n 个实数形式功率方程式。发电机 Pi、Qi 为正，负荷Pi、Qi为负。展开YV=I 为Ii=ΣYijVj=YiiVi+ ΣYijVi(i=123 ?n)(4)将式（4）代入式（3），得n维的非线性复数的电压方程组潮流计算的基本方程为(Pi-jQi)/Vi=YiiVi+ ΣYijVi （i=1,2, ? n） (5)（2.）变量的分类假设系统中有n个节点，构成n个复数方程，2n 个实数方程，变量总数为 6n 个。不可控变量（2n个）：负荷消耗的有功功率LiP和无功功率LiQ.由于该类变量无法控制，取决于用户，而且出现事先没有预计的变动，使系统偏离原始运行状态，因此又称为不可控变量或扰动变量。控制变量（2n个）:发电机发出的有功功率GiP和无功功率GiQ，因为该类变量可控。也称独立变量。状态变量（2n个）：母线电压或节点电压的幅值大小iV与相角大小iδ，又称依从变量或因变量。并且 iV 受GiP 控制，iδ受GiQ 控制。其中2n 个扰动变量是给定的，2n 个控制变量和2n 个状态变量中给定两个，求另外两个。（3.）变量的约束条件扰动变量没有约束条件。控制变量约束条件：为满足发电机的技术经济特性指标。状态变量的iV的约束条件：保证良好的电能质量。状态变量的iδ的约束条件：保证系统的稳定运行。(4.)系统节点的分类，根据给定的控制变量和状态变量进行分类如下：标准实用文案（1）PQ节点（即负荷节点）：给定Pi、Qi，求Vi和iδ（iie,f）。通常变电所都是这一类型的节点，由于没有发电设备，因而发电功率为零电力系统中的绝大多数节点属于这一节点。其包含变电站节点（即联络节点或浮游节点）。（2）PV节点（即调节节点、电压控制节点）：给定Pi和Vi，求Qi和iδ（iie,f）。这类节点必须有足够的可调无功容量，用以维持给定的电压幅值。一般时选择有一顶武功储备的发电厂和具有可调无功电源设备的变电所作为PV节点。在电力系统中，这类节点数很少。（3）平衡节点（即松弛节点、参考节点、基准节点）：给定Vi和iδ（iδ=0），求Pi和Qi。（只有一个）有功功率不能给定，这个节点承担了系统的有功功率平衡。同时其电压幅值也是给定的，相位为零。(5.)P-Q 分解法是从改进和简化牛顿法潮流程序的基础上提出来的，它的基本思想是：把节点功率表示为电压向量的极坐标方程式，抓住主要矛盾，以有功功率误差作为修正电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，把有功功率和无功功率迭代分开来进行。牛顿法潮流程序的核心是求解修正方程式，当节点功率方程式采取极坐标系标准实用文案统时，修正方程式展开为：P=HΔΘ+N V/VQ=JΔΘ+L V/V以上方程式是从数学上推倒出来的，并没有考虑电力系统这个具体对象的特点。电力系统中有功功率主要与各节点电压向量的角度有关，无功功率则主要受各节点电压幅值的影响。大量运算经验也告诉我们，矩阵 N及J中各元素的数值相对是很小的，因此对牛顿法的第一步简化就是把有功功率和无功功率分开来进行迭代，即将式(4)化简为：P=HΔΘQ=L V/V (5)这样，由于我们把2n 阶的线性方程组变成了二个 n阶的线性方程组，对牛顿法的第二个化简，也是比较关键的一个化简，即把式 (5)中的系数矩阵简化为在迭代过程中不变的对称矩阵。众所周知，一般线路两端电压的相角差是不大的(通常不超过10~20 度)，因此可以认为：(6)此外，与系统各节点无功功率相应的导纳LiB必定远远小于该节点自导纳的虚部，即：因此， (7)考虑到以上关系后，式(5)中系数矩阵中的元素表达式可以化简为：（8）这样，式(5)中系数矩阵可以表示为：标准实用文案（9）进一步可以把它们表示为以下矩阵的乘积：（10）将它代入(5)中，并利用乘法结合率，我们可以把修正方程式变为：将以上两式的左右两侧用以下矩阵左乘标准实用文案就可得到以上两式就是P-Q分解法达到修正方程式，其中系数矩阵只不过是系统导纳矩阵的虚部，因而是对称矩阵，而且在迭代过程中维持不变。它们与功率误差方程式构成了P-Q 分解法迭代过程中基本计算公式，其迭代步骤大致是：(1)给定各节点电压向量的电压初值Vi(0),θi(0);(2)根据(12)计算各节点有功功率误差Pi，并求出;Pi/Vi(3)解修正方程式(11)，并进而计算各节点电压向量角度的修正量iΔθ标准实用文案修正各节点电压向量角度θi；(5)根据式(16)计算各节点无功功率误差iQ，并求出/;iiQV(6)解修正方程式(11)，求出各节点电压幅值的修正量iV(7)修正各节点电压幅值iV(i)(i1)(i1)iiiV=V--V-(18)返回(2)进行迭代，直到各节点功率误差及电压误差都满足收敛条件。P-Q 分解法与牛顿法潮流程序的主要差别表现在它们的修正方程式上。 P-Q分解法通过对电力系统具体特点的分析，对牛顿法修正方程式的雅可比矩阵进行了有效的简化和改进，有以下三个特点：（1）在提高计算速度和减少内存方面的作用是明显的，不再叙述。（2）使我们得到以下好处。首先，因为修正方程式的系数矩阵就是导纳矩阵的虚部，因此在迭代过程中不必象牛顿法那样进行形成雅可比矩阵的计算，这样不仅是仅减少了运算量，而且也大大简化了程序。其次，由于系数矩阵在迭代过程中维持不变，因此在求解修正方程式时，可以迅速求得修正量，从而显著提高了迭代速度。（3）可以使我们减少形成因子表时的运算量，而且由于对称矩阵三角分解后，其上三角矩阵和下三角矩阵有非常简单的关系，所以在计算机中可以只储存上三角矩阵或下三角矩阵，从而也进一步节约了内存。三、设计内容本程序的功能是用牛顿——拉夫逊法进行潮流计算%B1矩阵：1、支路首端号； 2、末端号； 3、支路阻抗； 4、线路对地电纳(或变压器导纳)；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; %input(' 请输入平衡母线节点号：isb=');pr=0.0001; %input(' 请输入误差精度：pr=');%---------------------------------------------------n=6;%input(' 请输入节点数：n=');nl=6;%input('请输入支路数：nl=');B1=[120+0.04i01.0511;230.06+0.025i0+0.5i100;250.01+0.2i0+0.5i100;340.06+0.50i0+0.5i100;450.05+0.3i0100;650+0.02i01.0511]B2=[00101;03+1i1.0002;01.8+0.50i1.0002;00.6+0.8i1.0002;03.5+1.3i1.0002;0-5+0i1.0503]%input('请输入各节点参数形成的矩阵：B2=');%X=[10;20;30;40;50;60]%-------------------------------------------------------------%n=4;%input('请输入节点数：n=');nl=4;%input('请输入支路数：nl=');标准实用文案%B1=[124+16i0100;134+16i0100;232+8i0100;241.49+48.02i011/11001]%input(' 请输入由支路参数形成的矩阵： B1=');%B2=[0011501;0011002;020+4i11002;010+6i1002]%input(' 请输入各节点参数形成的矩阵：B2=');%-------------------------------------------------------------Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);%%%----------- 求导纳矩阵------------------------%fori=1:nifX(i,2)~=0;p=X(i,1);Y(p,p)=1/X(i,2);%end%endfori=1:nl %从1到n1（总支路数）ifB1(i,7)==1 %----------- 如果是变压器支路--------ifB1(i,6)==0 %左节点(首端)处于1侧p=B1(i,1);q=B1(i,2);else %左节点(首端)处于K侧p=B1(i,2);q=B1(i,1);endY(p,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)^2); %对角元K侧Y(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); %非对角元Y(q,q)=Y(q,q)+1./B1(i,3)+B1(i,4)./2.0000; %对角元j侧+线路电纳的一半Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2.0000; %对角元i侧+线路电纳的一半endenddisp('导纳矩阵Y=');disp(Y);%-----------给定各节点初始电压及给定各节点注入功率--------------------------G=real(Y);B=imag(Y);%分解出导纳阵的实部和虚部fori=1:n%给定各节点初始电压的实部和虚部e(i)=real(B2(i,3));f(i)=imag(B2(i,3));V(i)=abs(B2(i,3));%PV、平衡节点及PQ节点电压模值endfori=1:n%给定各节点注入功率S(i)=B2(i,1)-B2(i,2);%i节点注入功率SG-SLB(i,i)=B(i,i)+B2(i,4);%i节点无功补偿量(电纳值)end%==================用牛顿-拉夫逊法迭代求解非线性代数方程（功率方程）=======================P=real(S);Q=imag(S);%分解出各节点注入的有功和无功功率ICT1=0;IT2=1;N0=2*n;N1=N0+1;a=0;% 迭代次数ICT1、a；不满足收敛要求的节点数IT2标准实用文案whileIT2~=0 %N0=2*n 雅可比矩阵的阶数；N1=N0+1 扩展列IT2=0;a=a+1;JZ=['Jacobi 矩阵第(',num2str(a),') 次消去运算'];JZ1=['Jacobi 矩阵第(',num2str(a),')次回代运算'];JZ0=[' 功率方程第(',num2str(a),') 次差值：'];%---------------- 求取各个节点的功率及功率偏差及PV节点的电压偏差--------------------fori=1:n %n个节点2n行(每节点两个方程P和Q或U)p=2*i-1;m=p+1;C(i)=0;D(i)=0;forj1=1:n %第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)Q1=C(i)*f(i)-e(i)*D(i);% 节点功率Q计算fiΣ(Gij*ej-Bij*fj)-ei Σ(Gij*fj+Bij*ej)V2=e(i)^2+f(i)^2; %电压模平方%===求取功率差及PV节点电压模平方差=========ifi~=isb %非平衡节点(PQ或PV节点)ifB2(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)扩展列△Qelse %PV节点==================J(m,N1)=P(i)-P1; %PV节点有功功率差J(m,N1)扩展列△PJ(p,N1)=V(i)^2-V2; %PV节点电压模平方差J(p,N1)扩展列△Uend标准实用文案end %(ifi~=isb) 非平衡节点(PQ或PV节点)end %(fori=1:n) n个节点2n行(每节点两个方程P和Q或U)form=1:N0JJN1(m)=J(m,N1);enddisp(JZ0);disp(JJN1);%------------- 判断功率偏差量及PV节点的电压偏差量是否满足要求-----------------fork=3:N0 %除去平衡节点1、2号以外的所有节点DET=abs(J(k,N1));if DET>=pr; %PQ节点的功率偏差量及PV节点的电压偏差量是否满足要求IT2=IT2+1;% 不满足要求的节点数加1endendICT2(a)=IT2; %不满足要求的节点数；a为迭代次数ICT1=ICT1+1; %迭代次数ifICT2(a)==0; %当前不满足要求的节点数为零break %退出迭代运算end%-------------------- 以上为求取各个节点的功率及功率偏差及 PV节点的电压偏差-------------%================= 求 取 Jacobi 矩 阵 形 成 修 正 方 程===================fori=2:n %n个节点2n行(每节点两个方程P和Q或U)标准实用文案ifi~=isb %非平衡节点(PQ或PV节点)if B2(i,5)~=3 %下面是针对 PQ 节点来求取 Jacobi 矩阵的元素===========C(i)=0;D(i)=0;forj1=1:n %第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)endfor j1=2:n %第i行共n列(2n个Jacobi矩阵元素dP/de 及dP/df 或dQ/de 及dQ/df)ifj1~=isb&j1~=i %非平衡节点&非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); %X1=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; %X1=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 节点无功功率差J(p,N)=DQ;J(m,q)=X1;q=q+1; %X1=dP/de J(m,N)=DP 节点有功功率差J(m,N)=DP;J(p,q)=X4;J(m,q)=X2; %X4=dQ/df X2=dp/dfelseifj1==i&j1~=isb %非平衡节点&对角元X1=-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);%dP/df标准实用文案X3=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;% 扩展列△QJ(p,N)=DQ;m=p+1;J(m,q)=X1;q=q+1;J(p,q)=X4;% 扩展列△PJ(m,N)=DP;J(m,q)=X2;endendelse %ifB2(i,5)~=3 %否则（即为PV节点）%=============== 下面是针对 PV节点来求取 Jacobi 矩阵的元素===========forj1=1:nifj1~=isb&j1~=i %非平衡节点&非对角元X1=-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;elseifj1==i&j1~=isb %非平衡节点&对角元X1=-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);%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;endendend %(ifB2(i,5)~=3 else)end %(ifi~=isb)end %(fori=1:n)n 个节点2n行(每节点两个方程P和Q或U)JZ0=['形成的第(',num2str(a),') 次Jacobi矩阵：'];disp(JZ0);disp(J);%=============================== 以上为形成完整的 Jacobi 矩阵================================%==== 下面用高斯消去法对由 Jacobi 矩阵形成的修正方程进行求解(按列消去、回代)==========fork=3:N0 %N0=2*n （从第三行开始，第一、二行是平衡节点）fork1=k+1:N1 %从k+1列的Jacobi元素到扩展列的△P、△Q或△UJ(k,k1)=J(k,k1)./J(k,k);% 用K行K列对角元素去除K行K列后的非对角元素进行规格化endJ(k,k)=1; %对角元规格化K行K列对角元素赋1%================== 按 列 消 去 运 算==================================标准实用文案fork2=k+1:N0 %从k+1行到2*n最后一行fork3=k+1:N1 %从k2+1列到扩展列消去k+1行后各行下三角元素J(k2,k3)=J(k2,k3)-J(k2,k)*J(k,k3);% 消去运算end%用当前行K3列元素减去当前行K列元素乘以第k行K3列元素J(k2,k)=0;%当前行第k列元素已消为0endendJZ=['Jacobi 矩阵第(',num2str(a),') 次消去运算'];JZ1=['Jacobi 矩阵第(',num2str(a),')次回代运算'];disp(JZ);disp(J);%==================== 按 列 回 代 运 算=======================================fork=N0:-1:3fork1=k-1:-1:3J(k1,N1)=J(k1,N1)-J(k1,k)*J(k,N1);J(k1,k)=0;endendform=1:N0JJN1(m)=J(m,N1);enddisp(JZ1);disp(JJN1);%disp(J);%---------------------------------- 修改节点电压-------------------------------fork=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)^2+f(L)^2);enddisp('各个节点电压模');disp(U);%============================== 结 束 一 次 迭 代==============================end%********************** 下面为迭代计算结束后的有关输出过程 *****************disp('迭代次数：');disp(ICT1-1);disp('没有达到精度要求的个数：');disp(ICT2);fork=1:nV(k)=sqrt(e(k)^2+f(k)^2); %计算各节点电压的模值sida(k)=atan(f(k)./e(k))*180./pi; %计算各节点电压的角度E(k)=e(k)+f(k)*j; %将各节点电压用复数表示end%=============== 计算各输出量===========================disp('各节点的电压复数值E为(节点号从小到大排列)：');disp(E); %显示各节点的实际电压值E用复数表示disp('-----------------------------------------------------');标准实用文案disp('各节点的电压模值大小V为(节点号从小到大排列)：');disp(V); %显示各节点的电压大小V的模值disp('-----------------------------------------------------');disp('各节点的电压相角sida为(节点号从小到大排列)：');disp(sida); %显示各节点的电压相角forp=1:nC(p)=0;forq=1:nC(p)=C(p)+conj(Y(p,q))*conj(E(q)); %计算各节点注入电流的共轭值endS(p)=E(p)*C(p); %计算各节点的功率S=电压X注入电流的共轭值enddisp('各节点的功率S为(节点号从小到大排列)：');disp(S); %显示各节点的注入功率disp('-----------------------------------------------------');disp('各条支路的首端功率Si为(顺序同您输入B1时一致)：');fori=1:nlp=B1(i,1);q=B1(i,2);ifB1(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);elseifB1(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(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)^2)));Siz(i)=Si(p,q);endendZF=['S(',num2str(p),',',num2str(q),')=',num2str(Si(p,q))];disp(ZF);disp('-----------------------------------------------------');enddisp('各条支路的末端功率Sj为(顺序同您输入B1时一致)：');fori=1:nlp=B1(i,1);q=B1(i,2);ifB1(i,7)==0Sj(q,p)=E(q)*conj(E(q)*B1(i,4)./2+(E(q)-E(p))./B1(i,3));Sjy(i)=Sj(q,p);elseifB1(i,6)==0Sj(q,p)=E(q)*conj((E(q)-E(p)*B1(i,5))*(1./(B1(i,3)*B1(i,5)^2)));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(',num2str(q),',',num2str(p),')=',num2str(Sj(q,p))];disp(ZF);disp('-----------------------------------------------------');enddisp('各条支路的功率损耗DS为(顺序同您输入B1时一致)：');fori=1:nlp=B1(i,1);q=B1(i,2);DS(i)=Si(p,q)+Sj(q,p);ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DS(i))];disp(ZF);disp('-----------------------------------------------------');endzws=0;JDDY=[];JDP=[];JDQ=[];JDDYJD=[];fori=1:n %总网损为所有节点注入功率的代数和zws=zws+S(i);JDDYJD=strcat(JDDYJD,num2str(i),'(',num2str(sida(i)),'),');JDDY=strcat(JDDY,num2str(i),'(',num2str(V(i)),'),');JDP=strcat(JDP,num2str(i),'(',num2str(real(S(i))),'),');JDQ=strcat(JDQ,num2str(i),'(',num2str(imag(S(i))),'),');endJDDYJD=strcat(' 节点电压角度：',JDDYJD);标准实用文案JDDY=strcat(' 节点电压模值：',JDDY);JDP=strcat(' 节点有功：',JDP);JDQ=strcat(' 节点无功：',JDQ);ZFsubplot(4,1,2);plot(sida);title(JDDYJD);ylabel('节点电压角度');gridon;subplot(4,1,3);P=real(S);Q=imag(S);bar(P);title(JDP);ylabel('节点注入有功');gridon;subplot(4,1,4);bar(Q);title(JDQ);xlabel(ZF);ylabel(' 节点注入无功');gridon;%***********************************figure(2);subplot(3,2,1);JDH=[];JDH1=[];fori=1:nl标准实用文案JDH=strcat(JDH,num2str(i),'(',num2str(B1(i,1)),',',num2str(B1(i,2)),'),');JDH1=strcat(JDH1,num2str(i),'(',num2str(B1(i,2)),',',num2str(B1(i,1)),'),');endP1=real(Siz);Q1=imag(Siz);bar(P1);title(JDH);ylabel('支路首端注入有功');%xlabel(' 支路号');gridon;subplot(3,2,2);bar(Q1);title(JDH);ylabel('支路首端注入无功');gridon;subplot(3,2,3);P2=real(Sjy);Q2=imag(Sjy);bar(P2);title(JDH1);ylabel('支路末端注入有功');gridon;subplot(3,2,4);bar(Q2);title(JDH1);ylabel('支路末端注入无功');gridon;标准实用文案subplot(3,2,5);P3=real(DS);Q3=imag(DS);bar(P3);xlabel(JDH);ylabel(' 支路有功损耗');gridon;subplot(3,2,6);bar(Q3);xlabel(JDH);ylabel(' 支路无功损耗');gridon;四、结果分析及结论B1=Columns1through51.00002.00000+0.0400i01.05002.00003.00000.0600+0.0250i0+0.5000i1.00002.00005.00000.0100+0.2000i0+0.5000i1.00003.00004.00000.0600+0.5000i0+0.5000i1.00004.00005.00000.0500+0.3000i01.00006.00005.00000+0.0200i01.0500Columns6through71.00001.000000000000标准实用文案1.0000 1.0000B2=001.000001.000003.0000+1.0000i1.000002.000001.8000+0.5000i1.000002.000000.6000+0.8000i1.000002.000003.5000+1.3000i1.000002.00000-5.00001.050003.0000导纳矩阵Y=Columns1through50-22.6757i0+23.8095i0000+23.8095i14.4506-35.4047i-14.2012+5.9172i0-0.2494+4.9875i0-14.2012+5.9172i14.4378-7.3888i-0.2366+1.9716i000-0.2366+1.9716i0.7771-4.9649i-0.5405+3.2432i0-0.2494+4.9875i0-0.5405+3.2432i0.7899-57.9808i00000+47.6190iColumn6000标准实用文案0+47.6190i-45.3515i功率方程第(1)次差值：Columns1through1000-1.6905-3.0000-0.0000-1.8000-0.5500-0.6000-1.0500 -3.5000Columns11through120 5.0000形成的第(1)次Jacobi矩阵：Columns1through100000000000000000000000-36.095214.45065.9172-14.2012004.9875-0.249400-14.4506-34.714214.20125.9172000.24944.9875005.9172-14.2012-6.888814.43781.9716-0.2366000014.20125.9172-14.4378-7.88880.23661.971600标准实用文案00001.9716-0.2366-4.71490.77713.2432-0.540500000.23661.9716-0.7771-5.21490.54053.2432004.9875-0.2494003.2432-0.5405-57.73080.7899000.24944.9875000.54053.2432-0.7899-58.2308000000000000000000050.0000Columns11through1300000000-1.690500-3.000000-0.000000-1.800000-0.550000-0.600047.61900-1.0500047.6190-3.5000-2.100000标准实用文案0-47.61905.0000Jacobi矩阵第(1)次消去运算Columns1through1000000000000000000000001.0000-0.4003-0.16390.393400-0.13820.00690001.0000-0.2922-0.2865000.0431-0.125600001.0000-0.9301-0.21030.0252-0.14170.1629000001.00000.0875-0.1200-0.0124-00000-0.2070-0.79780.139300000001.00000.0256-0.7122000000001.0000-0.00600000000001.0000标准实用文案00000000000000000000Columns11throug0468000.057400-0.0428000.1919000.1739000.1652-0.880100.03550.0052-0.86740.09061.00000001.0000-0.1107Jacobi矩阵第(1)次回代运算Columns1through100 0 0.0730 0.1638 0.1855 0.1901 0.2361 0.16040.0355 -0.0055Columns11through120 -0.1107各个节点电压模标准实用文案00.94140.83640.78050.96451.0558功率方程第(2)次差值：Columns1through1000-0.62960.0415-0.3413-0.2986-0.2337-0.0345-0.0185 -0.1300Columns11through12-0.0123 0.1870形成的第(2)次Jacobi矩阵：Columns1through100000000000000000000000-29.504722.30853.1592-14.1342004.5828-1.048100-16.0823-31.403914.13423.1592001.04814.5828002.1200-12.6919-2.681014.86931.5609-0.5675000012.69192.1200-11.4594-3.86640.56751.56090000001.4681-0.4970-2.80871.95002.3906-0.9332标准实用文案00000.49701.4681-0.8303-4.52680.93322.3906004.8120-0.2133003.1311-0.5037-54.61943.9470000.21334.8120000.50373.13113.0557-57.2371000000000000000000-5.272750.0000Columns11through1300000000-0.6296000.041500-0.341300-0.298600-0.233700-0.034545.92980.2597-0.0185-0.259745.9298-0.1300-2.1000-0.2215-0.01230.2597-45.92980.1870Jacobi矩阵第(2)次消去运算标准实用文案Columns1through1000000000000000000000001.0000-0.7561-0.10710.479000-0.15530.03550001.0000-0.2849-0.2494000.0333-0.118300001.0000-1.9754-0.27810.1011-0.12440.2471000001.00000.0644-0.1101-0.0363-0.12790000001.0000-0.8086-1.04120.387000000001.0000-0.0399-0.5826000000001.0000-0.08010000000001.00000000000000标准实用文案0000000000Columns11throug021300-0.0088000.086300-0.0058000.135600-0.0104-0.9174-0.00520.0131-0.0643-0.84900.00541.00000.10550.005801.00000.0009Jacobi矩阵第(2)次回代运算Columns1through10000.05400.01630.1003-0.01440.1480-0.00590.0189 0.0065Columns11through120.0057 0.0009各个节点电压模0 0.8914 0.7355 0.6350 0.9457 1.0500功率方程第(3)次差值：标准实用文案Columns1through1000-0.04210.0138-0.0504-0.0538-0.0708-0.0116-0.0025 -0.0032Columns11through12-0.0000 0.0010形成的第(3)次Jacobi矩阵：Columns1through100000000000000000000000-26.571422.08592.6084-13.4636004.3093-1.115900-15.8976-30.042313.46362.6084001.11594.3093001.7315-11.1823-1.580413.76931.3666-0.5153000011.18231.7315-9.4503-3.90160.51531.36660000001.1777-0.4505-1.59801.86531.9138-0.834200000.45051.1777-0.6269-4.27700.83421.9138标准实用文案004.7162-0.2412003.0664-0.5146-53.45254.5052000.24124.7162000.51463.06642.8872-56.2049000000000000000000-5.230149.7265Columns11through1300000000-0.0421000.013800-0.050400-0.053800-0.070800-0.011645.0309-0.0509-0.00250.050945.0309-0.0032-2.0885-0.2197-0.0000-0.0509-45.03090.0010Jacobi矩阵第(3)次消去运算Columns1through10标准实用文案00000000000000000000001.0000-0.8312-0.09820.506700-0.16220.04200001.0000-0.2752-0.2465000.0338-0.115100001.0000-2.5639-0.33400.1260-0.14910.2918000001.00000.0616-0.0994-0.0316-0.11480000001.0000-1.4473-1.59280.648000000001.00000.0025-0.4904000000001.0000-0.09650000000001.000000000000000000000000标准实用文案Columns11throug001600-0.0009000.015100-0.0009000.063300-0.0046-0.95630.00110.0047-0.0761-0.84270.00111.00000.10520.000001.00000.0085Jacobi矩阵第(3)次回代运算Columns1through10000.01130.00580.0252-0.00390.0645-0.00060.0047 0.0082Columns11through12-0.0009 0.0085各个节点电压模0 0.8816 0.7102 0.5724 0.9410 1.0500功率方程第(4)次差值：Columns1through10标准实用文案0 0 -0.0010 0.0006 -0.0028 -0.0037 -0.0167 -0.0014-0.0003 0.0003Columns11through12-0.0001 -0.0022形成的第(4)次Jacobi矩阵：Columns1through100000000000000000000000-25.997022.12382.4585-13.3379004.2516-1.142300-15.9475-29.648213.33792.4585001.14234.2516001.6376-10.8022-1.319813.50301.3179-0.5017000010.80221.6376-8.9331-3.90230.50171.31790000001.0505-0.4340-1.01831.83201.7048-0.797400000.43401.0505-0.5540-4.21670.79741.7048004.6909-0.2808003.0468-0.5387-53.13374.9857标准实用文案000.28084.6909000.53873.04682.4260-55.9691000000000000000000-4.826549.7682Columns11through1300000000-0.0010000.000600-0.002800-0.003700-0.016700-0.001444.8086-0.4407-0.00030.440744.80860.0003-2.0903-0.2027-0.0001-0.4407-44.8086-0.0022Jacobi矩阵第(4)次消去运算Columns1through100 0 0 0 0 0 0 00 0标准实用文案0000000000001.0000-0.8510-0.09460.513100-0.16350.04390001.0000-0.2737-0.2462000.0339-0.114600001.0000-2.7663-0.35240.1342-0.15690.3075000001.00000.0607-0.0968-0.0304-0.11180000001.0000-2.4188-2.43691.057600000001.00000.0702-0.4464000000001.0000-0.11020000000001.000000000000000000000000Columns11through13000标准实用文案0 0 00 0 0.00000 0 -0.00000 0 0.00080 0 -0.00000 0 0.02200 0 -0.0019-1.0086 0.0099 0.0014-0.0855 -0.8393 0.00021.0000 0.0970 0.000001.00000.0022Jacobi矩阵第(4)次回代运算Columns1through10000.00220.00120.0048-0.00110.0207-0.00110.00140.0021Columns11through12-0.0002 0.0022各个节点电压模0 0.8796 0.7052 0.5522 0.9396 1.0500功率方程第(5)次差值：Columns1through100 0 -0.0000 0.0000 -0.0000 -0.0001 -0.0019 -0.0002-0.0000 0.0000标准实用文案Columns11through12-0.0000-0.0002形成的第(5)次Jacobi矩阵：Columns1through100000000000000000000000-25.889922.13672.4282-13.3134004.2402-1.147700-15.9556-29.562613.31342.4282001.14774.2402001.6248-10.7272-1.285813.44511.3087-0.4984000010.72721.6248-8.8356-3.89700.49841.30870000001.0099-0.4270-0.82561.81501.6381-0.782800000.42701.0099-0.5281-4.20510.78281.6381004.6833-0.2907003.0411-0.5445-53.03925.1057000.29074.6833000.54453.04112.3103-55.8956标准实用文案0 0 0 0 0 0 0 00 00 0 0 0 0 0 0 0-4.7229 49.7766Columns11through1300000000-0.0000000.000000-0.000000-0.000100-0.001900-0.000244.7409-0.5384-0.00000.538444.74090.0000-2.0906-0.1984-0.0000-0.5384-44.7409-0.0002Jacobi矩阵第(5)次消去运算Columns1through10标准实用文案00000000000000000000001.0000-0.8550-0.09380.514200-0.16380.04430001.0000-0.2735-0.2461000.0339-0.114500001.0000-2.7964-0.35490.1352-0.15810.3095000001.00000.0604-0.0965-0.0304-0.11130000001.0000-3.1204-3.05631.353600000001.00000.1136-0.4328000000001.0000-0.11520000000001.000000000000000000000000标准实用文案Columns11throug000000-0.0000000.0000000.0000000.003100-0.0002-1.04430.01260.0002-0.0918-0.83840.00001.00000.09490.000001.00000.0003Jacobi矩阵第(5)次回代运算Columns1through10000.00030.00010.0006-0.00020.0028-0.00020.0002 0.0