matlab潮流计算.pdf

前言电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。

在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。

可靠性和经济性。

此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。

所以潮流计算是研究电力系统的一种很重要和很基础的计算。

随着科学技术的发展，电力系统变得越来越复杂，电气工程师掌握一种好的能对电力系统进行仿真的软件是学习和研究的需要。

文章简要介绍了MATLAB发展历史、组成和强大的功能，并用简单例子分别就编程和仿真两方面分析了MATIAB软件在电力系统研究中的具体应用。

采取等效电路法，能对特殊、复杂地电力系统进行高效仿真研究，因此，掌握编程和仿真是学好MATLAB的基础。

与众多专门的电力系统仿真软件相比，MATLAB软件具有易学、功能强大和开放性好，是电力系统仿真研究的有力工具。

1电力系统的基本概念电力系统：发电机把机械能转化为电能，电能经变压器和电力线路输送并分配到用户，在那里经电动机、电炉和电灯等设备又将电能转化为机械能、热能和光能等。

这些生产、变换、输送、分配、消费电能的发电机、变压器、变换器、电力线路及各种用电设备等联系在一起组成的统一整体称为电力系统。

电力网：电力系统中除发电机和用电设备外的部分。

动力系统：电力系统和“动力部分”的总和。

2潮流计算2.1潮流计算概述与发展电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。

利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。

在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。

对潮流计算的要求可以归纳为下面几点：1）计算方法的可靠性或收敛性；2）对计算机内存量的要求；3）计算速度；4）计算的方便性和灵活性。

matlab潮流计算仿真方法
MATLAB 是一种强大的编程语言和环境，可用于执行各种仿真和计算任务，包括电力系统潮流计算。

以下是一个简单的 MATLAB 潮流计算仿真方法的
示例：
1. 定义系统参数：首先，你需要定义电力系统的参数，如发电机、负荷、变压器等。

这些参数通常包括额定电压、额定功率、电抗、电阻等。

2. 建立系统模型：使用这些参数，你可以在 MATLAB 中建立电力系统的模型。

这通常涉及到定义节点和支路，以及为它们分配相应的参数。

3. 编写潮流计算函数：接下来，你需要编写一个用于执行潮流计算的函数。

这个函数应该能够接收系统的模型和参数，并返回计算出的潮流结果，如电压、电流、功率等。

4. 运行仿真：最后，你可以运行仿真并调用你编写的潮流计算函数。

这将返回计算出的潮流结果，你可以使用这些结果进行进一步的分析或可视化。

这只是一个简单的示例，实际上在编写 MATLAB 潮流计算仿真方法时可能
需要考虑更多因素，例如系统的约束条件、初始条件、迭代算法的收敛性等。

如果你需要具体的 MATLAB 代码示例或更详细的指导，我建议你查阅MATLAB 的官方文档或相关的教程和文献。

前言电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。

在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。

可靠性和经济性。

此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。

所以潮流计算是研究电力系统的一种很重要和很基础的计算。

随着科学技术的发展，电力系统变得越来越复杂，电气工程师掌握一种好的能对电力系统进行仿真的软件是学习和研究的需要。

文章简要介绍了MATLAB发展历史、组成和强大的功能，并用简单例子分别就编程和仿真两方面分析了MATIAB软件在电力系统研究中的具体应用。

采取等效电路法，能对特殊、复杂地电力系统进行高效仿真研究，因此，掌握编程和仿真是学好MATLAB的基础。

与众多专门的电力系统仿真软件相比，MATLAB软件具有易学、功能强大和开放性好，是电力系统仿真研究的有力工具。

1电力系统的基本概念电力系统：发电机把机械能转化为电能，电能经变压器和电力线路输送并分配到用户，在那里经电动机、电炉和电灯等设备又将电能转化为机械能、热能和光能等。

这些生产、变换、输送、分配、消费电能的发电机、变压器、变换器、电力线路及各种用电设备等联系在一起组成的统一整体称为电力系统。

电力网：电力系统中除发电机和用电设备外的部分。

动力系统：电力系统和“动力部分”的总和。

2潮流计算2.1潮流计算概述与发展电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。

利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。

在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。

对潮流计算的要求可以归纳为下面几点：1）计算方法的可靠性或收敛性；2）对计算机内存量的要求；3）计算速度；4）计算的方便性和灵活性。

电力系统分析大作业——潮流计算班级：XXX姓名：XXX学号：XXX程序清单：%直角坐标法求解潮流计算clear;close;clc;Branch=[1 2 0.1+0.4j 0.01528j 1 03 1 0.3j inf 1.1 11 4 0.12+0.50j 0.01920j 1 02 4 0.08+0.40j 0.01413j 1 0]%Branch矩阵：1、支路首端号；2、支路末端号；3、支路阻抗；4、支路对地导纳；5、支路的变化；6、支路首端处于K侧为1,1侧为0Y=zeros(4); %节点导纳矩阵for i=1:4if Branch(i,6)==0 %不含变压器的支路j=Branch(i,1);k=Branch(i,2);Y(j,k)=Y(j,k)-1/Branch(i,3);Y(k,j)=Y(j,k);Y(j,j)=Y(j,j)+1/Branch(i,3)+Branch(i,4);Y(k,k)=Y(k,k)+1/Branch(i,3)+Branch(i,4);elsej=Branch(i,1);k=Branch(i,2);Y(j,k)=Y(j,k)-Branch(i,5)/Branch(i,3);Y(k,j)=Y(j,k);Y(j,j)=Y(j,j)+1/Branch(i,3);Y(k,k)=Y(k,k)+Branch(i,5)^2/Branch(i,3);endenddisp('节点导纳矩阵');YG=real(Y);B=imag(Y);V=[1;1;1.1;1.05];%给定V的初始计算值e=real(V);f=imag(V);disp('节点电压的实部：')edisp('节点电压的虚部：')fdisp('节点注入有功功率：')Ps=[-0.3;-0.55;0.5;0]disp('节点注入无功功率：')Qs=[-0.18;-0.13;0;0]%由各节点电压向量（状态变量）可得各节点注入功率：P=e.*(G*e-B*f)+f.*(G*f+B*e);Q=f.*(G*e-B*f)-e.*(G*f+B*e);del_W=[-0.30-P(1);-0.18-Q(1);-0.55-P(2);-0.13-Q(2);0.5-P(3);1.1^2-e(3)^2-f(3)^2];n=0;%辅助循环计数变量while (any(del_W>1e-5))&(n<5)n=n+1;disp('迭代次数：');disp(n)%--------------------------------------------------------%修正方程式：% [-0.30-P(1) ] [del_e(1)]% [-0.18-Q(1) ] [del_f(1)]% [-0.55-P(2) ] [del_e(2)]% [-0.13-Q(2) ] + J * [del_f(2)] = 0 % [0.5-P(3) ] [del_e(3)]% [1.1^2-e(3)^2-f(3)^2] [del_f(3)]%雅克比矩阵为：J=[P1_e1 P1_f1 P1_e2 P1_f2 P1_e3 P1_e3;% Q1_e1 Q1_f1 Q1_e2 Q1_f2 Q1_e3 Q1_f3% P2_e1 P2_f1 P2_e2 P2_f2 P2_e3 P2_f3% Q2_e1 Q2_f1 Q2_e2 Q2_f2 Q2_e3 Q2_f3% P3_e1 P3_f1 P3_e2 P3_f2 P3_e3 P3_f3% V3_e1 V3_f1 V3_e2 V3_f2 V3_e3 V3_f3]%--------------------------------------------------------%-------------------求雅可比矩阵的参数-----------------------%GeBf=G*e-B*f;%辅助计算函数GfBe=G*f+B*e;%辅助计算函数for i=1:2for j=1:3if i==jP_e(i,j)=-GeBf(i)-G(i,j)*e(i)-B(i,j)*f(i);P_f(i,j)=-GfBe(i)-G(i,j)*f(i)+B(i,j)*e(i);Q_e(i,j)=GfBe(i)+B(i,j)*e(i)-G(i,j)*f(i);Q_f(i,j)=-GeBf(i)+G(i,j)*e(i)+B(i,j)*f(i);elseP_e(i,j)=-G(i,j)*e(i)-B(i,j)*f(i);Q_f(i,j)=-P_e(i,j);P_f(i,j)=B(i,j)*e(i)-G(i,j)*f(i);Q_e(i,j)=P_f(i,j);endendendfor j=1:2P_e(3,j)=-G(3,j)*e(3)-B(3,j)*f(3);P_f(3,j)=B(3,j)*e(3)-G(3,j)*f(3);V_e(3,j)=0;V_f(3,j)=0;endP_e(3,3)=-GeBf(3)-G(3,3)*e(3)-B(3,3)*f(3);P_f(3,3)=-GfBe(3)-G(3,3)*f(3)+B(3,3)*e(3);V_e(3,3)=-2*e(3);V_f(3,3)=-2*f(3);%----------------------------------------------------%%--------------------求出雅克比矩阵-------------------% J=zeros(6,6);for i=1:3for j=1:3J(2*i-1,2*j-1)=P_e(i,j);J(2*i-1,2*j)=P_f(i,j);endendfor i=1:2for j=1:3J(2*i,2*j-1)=Q_e(i,j);J(2*i,2*j)=Q_f(i,j);endendfor j=1:3J(6,2*j-1)=V_e(3,j);J(6,2*j)=V_f(3,j);enddisp('雅可比矩阵为：')J%-----------------------------------------------------%%------------解修正方程并得出修正后的电压向量-----------% del_V=-J\del_W;clear i;for k=1:3V(k)=V(k)+del_V(2*k-1)+i*del_V(2*k);enddisp('节点电压为：')disp(V)%-----------------------------------------------------%%------------------计算节点不平衡量--------------------%e=real(V);f=imag(V);P=e.*(G*e-B*f)+f.*(G*f+B*e);Q=f.*(G*e-B*f)-e.*(G*f+B*e);del_W=[-0.30-P(1);-0.18-Q(1);-0.55-P(2);-0.13-Q(2);0.5-P(3);1.1^2-e(3)^2-f(3)^2];disp('节点不平衡量为：')disp(del_W)end%------------------------最终结论-------------------------%disp('最终各节点电压幅值为：')disp(abs(V))disp('最终各节点电压相角(度)为：')disp(180*angle(V)/pi)disp('最终各节点注入功率：')S=P+i*Qdisp('最终各节点注入有功功率为：')Pdisp('最终各节点注入无功功率为：')Q运行结果：Branch =1.00002.0000 0.1000 + 0.4000i 0 + 0.0153i 1.0000 03.0000 1.0000 0 + 0.3000i Inf 1.1000 1.00001.0000 4.0000 0.1200 + 0.5000i 0 + 0.0192i1.0000 02.0000 4.0000 0.0800 + 0.4000i 0 + 0.0141i 1.0000 0节点导纳矩阵Y =1.0421 - 8.2429i -0.5882 +2.3529i 0 +3.6667i -0.4539 + 1.8911i-0.5882 + 2.3529i 1.0690 - 4.7274i 0 -0.4808 + 2.4038i0 + 3.6667i 0 0 - 3.3333i 0-0.4539 + 1.8911i -0.4808 + 2.4038i 0 0.9346 - 4.2616i节点电压的实部：e =1.00001.00001.10001.0500节点电压的虚部：f =节点注入有功功率：Ps =-0.3000-0.55000.5000节点注入无功功率：Qs =-0.1800-0.1300迭代次数：1雅可比矩阵为：J =-1.0194 -8.3719 0.5882 2.3529 0 3.6667 -8.1138 1.0648 2.3529 -0.5882 3.6667 00.5882 2.3529 -1.0450 -4.8770 0 02.3529 -0.5882 -4.5778 1.0930 0 00 4.0333 0 0 0 -3.66670 0 0 0 -2.2000 0节点电压为：0.9935 - 0.0088i0.9763 - 0.1078i1.1000 + 0.1267i1.0500节点不平衡量为：-0.0013-0.0028-0.0135-0.05470.0030-0.0160迭代次数：2雅可比矩阵为：J =-0.8091 -8.3613 0.6052 2.3325 0.0324 3.6429 -7.9992 1.4071 2.3325 -0.6052 3.6429 -0.03240.8280 2.2338 -1.0190 -4.6364 0 02.2338 -0.8280 -4.3641 2.0878 0 0-0.4644 4.0333 0 0 -0.0324 -3.64290 0 0 0 -2.2000 -0.2533节点电压为：0.9847 - 0.0086i0.9590 - 0.1084i1.0924 + 0.1289i1.0500节点不平衡量为：-0.0000-0.0000-0.0003-0.00110.0001-0.0001迭代次数：3雅可比矩阵为：J =-0.7940 -8.2936 0.5995 2.3120 0.0315 3.6107 -7.9228 1.4000 2.3120 -0.5995 3.6107 -0.03150.8191 2.1927 -0.9865 -4.6144 0 02.1927 -0.8191 -4.2210 2.0885 0 0-0.4728 4.0056 0 0 -0.0315 -3.61070 0 0 0 -2.1849 -0.2579。

前言电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。

在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。

可靠性和经济性。

此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。

所以潮流计算是研究电力系统的一种很重要和很基础的计算。

随着科学技术的发展，电力系统变得越来越复杂，电气工程师掌握一种好的能对电力系统进行仿真的软件是学习和研究的需要。

文章简要介绍了MATLAB发展历史、组成和强大的功能，并用简单例子分别就编程和仿真两方面分析了MATIAB软件在电力系统研究中的具体应用。

采取等效电路法，能对特殊、复杂地电力系统进行高效仿真研究，因此，掌握编程和仿真是学好MATLAB的基础。

与众多专门的电力系统仿真软件相比，MATLAB软件具有易学、功能强大和开放性好，是电力系统仿真研究的有力工具。

1电力系统的基本概念电力系统：发电机把机械能转化为电能，电能经变压器和电力线路输送并分配到用户，在那里经电动机、电炉和电灯等设备又将电能转化为机械能、热能和光能等。

这些生产、变换、输送、分配、消费电能的发电机、变压器、变换器、电力线路及各种用电设备等联系在一起组成的统一整体称为电力系统。

电力网：电力系统中除发电机和用电设备外的部分。

动力系统：电力系统和“动力部分”的总和。

2潮流计算2.1潮流计算概述与发展电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。

利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。

在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。

对潮流计算的要求可以归纳为下面几点：1）计算方法的可靠性或收敛性；2）对计算机内存量的要求；3）计算速度；4）计算的方便性和灵活性。

用matlab潮流计算(牛顿拉夫逊法)%主程序：[d]=uigetfile('ieee14.m','Select Data File');if pathname==0error('you must select a valid data file')elsel(dfile);% strip off .meval(d)); %打开数据文件endglobal n;global m;[nb,mb]=size(bus); %节点重新编号[nl,ml]=size(line);nSW=0; % 平衡节点数目nPV=0; % PV节点数目nPQ=0; % PQ节点数目for i=1:nb, % nb为总节点数type=bus(i,6);if type==3,nSW=nSW+1; % 统计平衡节点数目SW(nSW,:)=bus(i,:);elseif type==2,nPV=nPV+1; % 统计PV节点数目PV(nPV,:)=bus(i,:);elsenPQ=nPQ+1; % 统计PQ节点数目PQ(nPQ,:)=bus(i,:);endend;bus=[PQ;PV;SW];newbus=[1:nb]';f=bus(:,1);nodenum=[newbus bus(:,1)];bus(:,1)=newbus;for i=1:nlfor j=1:2for k=1:nbif line(i,j)==nodenum(k,2)line(i,j)=nodenum(k,1);breakendendendendY=y(bus,line); %形成节点导纳矩阵K=0; %迭代次数初值Kmax=10; %最大迭代次数eps1=1.0e-10;eps2=1.0e-10;m=nPQ;n=nb;Um=eye(m,m);myf=fopen('output1.dat','w');for K=1:Kmaxfor i=1:mfor j=1:mif i==jUm(i,j)=bus(i,2);endendendb=dPQ(Y,bus);C=jac(bus,Y);dX=C\b';dx=dX';[nx,mx]=size(dx);for i=1:n-1 %计算相角bus(i,3)=bus(i,3)-dX(i,1);endB=dx(nx,n:mx)*Um; %计算电压差bus(1:m,2)=bus(1:m,2)-B'; %计算电压值dx(nx,n:mx)=B;fprintf(myf,'--第%d次迭代时雅可比矩阵--',K)fprintf(myf, '\n');for i=1:(n+m-1)for j=1:(n+m-1)fprintf(myf,'%8.6f', C(i,j));fprintf(myf, ' ');endfprintf(myf, '\n');endfprintf(myf,'--第%d次迭代时dPQ的误差--',K)fprintf(myf, '\n');for i=1:(n+m-1)fprintf(myf,'%8.6e', b(1,i));fprintf(myf, '\n');endfprintf(myf, '\n');fprintf(myf,'--第%d次迭代时dx(误差)--',K)fprintf(myf, '\n');for i=1:(n+m-1)fprintf(myf,'%8.6e', dX(i,1));fprintf(myf, '\n');endfprintf(myf, '\n');fprintf(myf,'第%d次迭代后节点电压(仅PQ节点)',K)fprintf(myf, '\n');for i=1:mfprintf(myf,'%8.6f', bus(i,2));fprintf(myf, '\n');endfprintf(myf, '\n');fprintf(myf,'第%d次迭代后相角(角度)',K)fprintf(myf, '\n');for i=1:nfprintf(myf,'%8.6f', bus(i,3)*180/pi);fprintf(myf, '\n');endfprintf(myf, '\n');if (max(abs(dx))<eps1)&(max(abs(b))<eps2) %判断是否达到计算精度break;endend%计算功率P1=0;T=0; %计算平衡节点的有功for j=1:nT=bus(n,2)*bus(j,2)*(real(Y(n,j))*cos(bus(n,3)-bus(j,3))+imag(Y(n,j))*sin(bus(n,3)-bus(j,3)));P1=P1+T;endbus(n,4)=P1;for k=m+1:n %计算PV节点、平衡节点的无功Q1=0;T=0;for j=1:nT=bus(k,2)*bus(j,2)*(real(Y(k,j))*sin(bus(k,3)-bus(j,3))-imag(Y(k,j))*cos(bus(k,3)-bus(j,3)));Q1=Q1+T;endbus(k,5)=Q1;endbus(:,1)=f; %换回各节点、支路的初始编号r=zeros(1,mb);for t=1:nfor l=t+1:nif bus(t,1)>bus(l,1)r=bus(t,:);bus(t,:)=bus(l,:);bus(l,:)=r;endendendfor i=1:nlfor j=1:2for k=1:nbif line(i,j)==nodenum(k,1)line(i,j)=nodenum(k,2);breakendendendendfclose(myf);Pf=loss(bus,line); %计算支路潮流及损耗%将节点导纳矩阵、节点潮流计算结果写入文件output2myf=fopen('output2.dat','w');fprintf(myf, '--节点导纳矩阵--\n');for k=1:nfor j=1:nfprintf(myf,'%8.6f', real(Y(k,j)));fprintf(myf, '+i*');fprintf(myf,'%8.6f', imag(Y(k,j)));fprintf(myf, ' ');endfprintf(myf, '\n');endfprintf(myf, '------------------牛顿－拉夫逊法潮流计算结果-------------\n');fprintf(myf, '------------节点计算结果------------\n');fprintf(myf, '-- 节点节点电压节点相角注入有功功率（P）注入无功功率(Q) 类型--\n');for l=1:nbfor j=1:mbif j==1|j==6fprintf(myf, ' %8.1f ', bus(l,j));elseif j==3fprintf(myf, ' %8.6f ', bus(l,j)*180/pi);elsefprintf(myf, ' %8.6f ', bus(l,j));endendfprintf(myf, ' \n');endfprintf(myf, '--支路计算结果--\n');fprintf(myf, '-- 节点(I) 节点(J) 线路功率S(I,J) 线路功率S(J,I) 线路损耗dS(I,J)--\n');for k=1:nlfor j=1:5if j<=2fprintf(myf,'%8.1f ', Pf(k,j));fprintf(myf, ' ');elsefprintf(myf,'%8.6f', real(Pf(k,j)));fprintf(myf, '+i*');fprintf(myf,'%8.6f', imag(Pf(k,j)));fprintf(myf, ' ');endendfprintf(myf, '\n');endfclose(myf);%根据支路参数建立节点导纳矩阵程序：function Y=y(bus,line)%目的：根据支路参数建立节点导纳矩阵%输入：节点参数矩阵--bus；支路参数矩阵--line%输出：节点导纳矩阵--Y[nb,mb]=size(bus);[nl,ml]=size(line);Y=zeros(nb,nb);for k=1:nlI=line(k,1);J=line(k,2);Zt=line(k,3)+j*line(k,4);if Zt==0disp('对地支路');Yt=inf;elseYt=1/Zt;endYm=line(k,5)+j*line(k,6);K=line(k,7);if(K==0)&(J~=0) %普通线路Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt+Ym;Y(I,J)=Y(I,J)-Yt;Y(J,I)=Y(I,J);endif(K==0)&(J==0) %对地支路K=J=0，R=X=0 Y(I,I)=Y(I,I)+Ym;endif K>0 %K>0时变压器支路Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt/K^2;Y(I,J)=Y(I,J)-Yt/K;Y(J,I)=Y(I,J);endif K<0 %K<0时变压器支路Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt*K^2;Y(I,J)=Y(I,J)+Yt*K;Y(J,I)=Y(I,J);endend%形成雅克矩阵程序：function J=jac(bus,Y)%形成雅可比矩阵global n;global m;for i=1:(n-1) %计算PQ、PV节点的有功P1=0;T=0;for j=1:nT=bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)-bus(j,3))+imag(Y(i,j))*sin(bus(i,3)-bus(j,3)));P1=P1+T;endbus(i,4)=P1;endfor i=1:n-1 %计算PV、PQ节点的无功Q1=0;T=0;for j=1:nT=bus(i,2)*bus(j,2)*(real(Y(i,j))*sin(bus(i,3)-bus(j,3))-imag(Y(i,j))*cos(bus(i,3)-bus(j,3)));Q1=Q1+T;endbus(i,5)=Q1;endfor i=1:n-1 %计算Hfor j=1:n-1if i~=jH(i,j)=-bus(i,2)*bus(j,2)*(real(Y(i,j))*sin(bus(i,3)-bus(j,3))-imag(Y(i,j))*cos(bus(i,3)-bus(j,3)));N(i,j)=-bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)-bus(j,3))+imag(Y(i,j))*sin(bus(i,3)-bus(j,3)));K(i,j)=bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)-bus(j,3))+imag(Y(i,j))*sin(bus(i,3)-bus(j,3)));L(i,j)=H(i,j);endendendfor i=1:n-1 %计算Hfor j=1:n-1if j==iH(i,i)=bus(i,5)+imag(Y(i,i))*bus(i,2)^2;N(i,i)=-bus(i,4)-real(Y(i,i))*bus(i,2)^2;K(i,i)=N(i,i)+2*real(Y(i,i))*bus(i,2)^2;L(i,i)=-H(i,i)+2*imag(Y(i,i))*bus(i,2)^2;endendendN=N(1:n-1,1:m);K=K(1:m,1:n-1);L=L(1:m,1:m);J=[H,N;K,L];%计算dPQ的程序：function dPQ=dPQ(Y,bus)%delP--有功偏差量%delQ--无功偏差量%Y--节点导纳矩阵%bus--节点参数（P,Q,U及相角）矩阵global n;global m;delP=zeros(1,n-1);delQ=zeros(1,m);for i=1:(n-1) %形成delP矩阵(PQ、PV节点)P1=0;T=0;for j=1:nT=bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)-bus(j,3))+imag(Y(i,j))*sin(bus(i,3)-bus(j,3)));P1=P1+T;enddelP(1,i)=bus(i,4)-P1;endfor i=1:m %形成delQ矩阵(PQ节点)Q1=0;T=0;for j=1:nT=bus(i,2)*bus(j,2)*(real(Y(i,j))*sin(bus(i,3)-bus(j,3))-imag(Y(i,j))*cos(bus(i,3)-bus(j,3)));Q1=Q1+T;enddelQ(1,i)=bus(i,5)-Q1;enddPQ=[delP,delQ];%计算线路损耗、线路潮流程序：function Pf=loss(bus,line)%计算线路损耗、线路潮流[nl,ml]=size(line);Pf=zeros(nl,5);for k=1:nlI=line(k,1);J=line(k,2);Zt=line(k,3)+i*line(k,4);if Zt==0Yt=inf;elseYt=1/Zt;endYm=line(k,5)+i*line(k,6);K=line(k,7);if (K==0)&(J~=0) %普通线路潮流S(I,J)=bus(I,2)^2*(conj(Yt)+conj(Ym))-bus(I,2)*(cos(bus(I,3))+i*sin(bus(I,3)))*bus(J,2)*(cos(bus(J,3))-i*sin(bus(J,3))) *conj(Yt);S(J,I)=bus(J,2)^2*(conj(Yt)+conj(Ym))-bus(J,2)*(cos(bus(J,3))+i*sin(bus(J,3)))*bus(I,2)*(cos(bus(I,3))-i*sin(bus(I,3))) *conj(Yt);delS(I,J)=S(I,J)+S(J,I);endif(K==0)&(J==0) %对地支路潮流J=5;S(I,5)=bus(I,2)^2*conj(Ym);endif K>0 %变压器支路k>0时的潮流S(I,J)=bus(I,2)^2*(conj(Ym+Yt*(1-1/K))+conj(Yt/K))-bus(I,2)*(cos(bus(I,3))+i*sin(bus(I,3)))*bus(J,2)*(cos(bus(J,3))-i *sin(bus(J,3)))*conj(Yt/K);S(J,I)=bus(J,2)^2*(conj(Yt))/K^2-bus(J,2)*(cos(bus(J,3))+i*sin(bus(J,3)))*bus(I,2)*(cos(bus(I,3))-i*sin(bus(I,3)))*conj( Yt/K);delS(I,J)=S(I,J)+S(J,I);endif K<0 %变压器支路k<0时的潮流S(I,J)=bus(I,2)^2*(conj(Ym+Yt))+bus(I,2)*(cos(bus(I,3))+i*sin(bus(I,3)))*bus(J,2)*(cos(bus(J,3))-i*sin(bus(J,3)))*conj( Yt*K);S(J,I)=bus(J,2)^2*(conj(Yt))*K^2+bus(J,2)*(cos(bus(J,3))+i*sin(bus(J,3)))*bus(I,2)*(cos(bus(I,3))-i*sin(bus(I,3)))*conj (Yt*K);delS(I,J)=S(I,J)+S(J,I);endif J==5&Zt==0Sp=[line(k,1) line(k,2) S(I,5) 0 S(I,5)];elseSp=[line(k,1) line(k,2) S(I,J) S(J,I) delS(I,J)];endPf(k,:)=Sp;end%输入的参数数据：% data for test case%各节点参数：节点编号,注入有功,注入无功,(Sn=100MV A)电压幅值,电压相位,类型%类型：1=PQ节点,2=PV节点,3=平衡节点% (bus#)( volt )( ang )( p )( q )(bus type)bus=[1,1.0,0.0,-0.478,0.039,1;2,1.0,0.0,-0.076,-0.016,1;3,1.0,0.0,0.0,0.0,1;4,1.0,0.0,-0.295,-0.166,1;5,1.0,0.0,-0.09,-0.058,1;6,1.0,0.0,-0.035,-0.018,1;7,1.0,0.0,-0.061,-0.016,1;8,1.0,0.0,-0.135,-0.058,1;9,1.0,0.0,-0.149,-0.05,1;10,1.045,0.0,0.183,0.0,2;11,1.010,0.0,-0.942,0.0,2;12,1.70,0.0,-0.112,0.047,2;13,1.90,0.0,0.0,0.174,2;14,1.060,0.0,0.0,0.0,3];%各支路参数：起点编号,终点编号,电阻,电抗,电导,电纳line = [1,2,0.01335,0.04211,0.0,0.0,0;1,3,0.0,0.20912,0.0,0.0,0;1,4,0.0,0.55618,0.0,0.0,0;1,10,0.05811,0.17632,0.0,0.0340,0;1,11,0.06701,0.17103,0.0,0.0128,0;2,10,0.05695,0.17388,0.0,0.0346,0;2,12,0.0,0.25202,0.0,0.0,0;2,14,0.05403,0.22304,0.0,0.0492,0;3,4,0.0,0.11001,0.0,0.0,0;3,13,0.0,0.17615,0.0,0.0,0;4,5,0.03181,0.08450,0.0,0.0,0;4,9,0.12711,0.27038,0.0,0.0,0;5,6,0.08205,0.19207,0.0,0.0,0;6,12,0.09498,0.19890,0.0,0.0,0;7,8,0.22092,0.19988,0.0,0.0,0;7,12,0.12291,0.25581,0.0,0.0,0;8,9,0.17093,0.34802,0.0,0.0,0;8,12,0.06615,0.13027,0.0,0.0,0;10,11,0.04699,0.19797,0.0,0.0438,0;10,14,0.01938,0.05917,0.0,0.0528,0];输出结果数据1：--第1次迭代时雅可比矩阵---38.624033 21.578554 4.781943 1.797979 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 5.346051 5.119505 -0.000000 -0.000000 -10.417258 6.840981 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.00000021.578554 -38.240787 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 5.427654 -0.000000 6.745496 -0.000000 6.840981 -9.429913 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.0000004.781943 -0.000000 -24.658288 9.090083 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 10.786262 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.0000001.797979 -0.000000 9.090083 -24.282506 10.365394 -0.000000 -0.000000 -0.000000 3.029050 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -5.326055 3.902050 -0.000000 -0.000000 -0.000000 1.424005-0.000000 -0.000000 -0.000000 10.365394 -14.768338 4.402944 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 3.902050 -5.782934 1.880885 -0.000000 -0.000000 -0.000000-0.000000 -0.000000 -0.000000 -0.000000 4.402944 -11.362870 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 6.959926 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 1.880885 -2.467393 -0.000000 -0.000000 -0.000000-0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -7.651113 2.251975 -0.000000 -0.000000 -0.000000 5.399139 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -2.946815 2.489025 -0.000000-0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 2.251975 -14.941622 2.314963 -0.000000 -0.000000 10.374684 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 2.489025 -4.555697 1.136994-0.000000 -0.000000 -0.000000 3.029050 -0.000000 -0.000000 -0.000000 2.314963 -5.344014 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 1.424005 -0.000000 -0.000000 -0.000000 1.136994 -2.5610005.346051 5.427654 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -32.727644 5.047017 -0.000000 -0.000000 1.761905 1.777691 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.0000005.119505 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 5.047017 -10.166523 -0.000000 -0.000000 2.005835 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000-0.000000 6.745496 -0.000000 -0.000000 -0.000000 6.959926 5.399139 10.374684 -0.000000 -0.000000 -0.000000 -29.479246 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 3.323549 2.594145 5.268177 -0.000000-0.000000 -0.000000 10.786262 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -10.786262 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.00000010.608721 -6.840981 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -1.761905 -2.005835 0.000000 0.000000 -37.968631 21.578554 4.781943 1.797979 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000-6.840981 9.706123 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -1.777691 0.000000 0.000000 0.000000 21.578554 -31.542421 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 4.781943 -0.000000 -14.439724 9.090083 -0.000000 -0.000000 -0.000000 -0.000000 -0.0000000.000000 0.000000 0.000000 5.326055 -3.902050 0.000000 0.000000 0.000000 -1.424005 0.000000 0.000000 0.000000 0.000000 1.797979 -0.000000 9.090083 -24.282506 10.365394 -0.000000 -0.000000 -0.000000 3.0290500.000000 0.000000 0.000000 -3.902050 5.782934 -1.880885 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 -0.000000 10.365394 -14.768338 4.402944 -0.000000 -0.000000 -0.0000000.000000 0.000000 0.000000 0.000000 -1.880885 5.204433 0.000000 0.000000 0.000000 0.000000 0.000000 -3.323549 0.000000 -0.000000 -0.000000 -0.000000 -0.000000 4.402944 -5.631166 -0.000000 -0.000000 -0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 5.083169 -2.489025 0.000000 0.000000 0.000000 -2.594145 0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -3.204764 2.251975 -0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -2.489025 8.894195 -1.136994 0.000000 0.000000 -5.268177 0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 2.251975 -6.397765 2.3149630.000000 0.000000 0.000000 -1.424005 0.000000 0.000000 0.000000 -1.136994 2.561000 0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 -0.000000 3.029050 -0.000000 -0.000000 -0.000000 2.314963 -5.344014--第1次迭代时dPQ的误差---3.822688e-0016.210513e-0020.000000e+000-2.950000e-001-9.000000e-0021.333520e+0001.007177e+0002.034249e+000-1.490000e-0016.056626e-002-9.219354e-001-7.942109e+0000.000000e+0003.667009e-0013.333183e+0005.109282e+000-1.660000e-001-5.800000e-0022.847852e+0002.207175e+0004.213929e+000-5.000000e-002--第1次迭代时dx(误差)---7.699084e-001-8.544764e-001-1.189723e+000-1.410571e+000-1.585607e+000-1.994895e+000-2.196974e+000-2.162427e+000-1.721659e+000-4.249173e-001-6.178169e-001-2.246296e+000-1.189723e+000-5.568104e-001-5.586033e-001-1.299237e+000-1.208867e+000-1.285646e+000-1.499291e+000-2.011550e+000-1.901143e+000-1.488513e+000第1次迭代后节点电压(仅PQ节点)1.5568101.5586032.2992372.2088672.2856462.4992913.0115502.9011432.488513第1次迭代后相角(角度)44.11250048.95789268.16611180.81979290.848607114.299041125.877362123.89794698.64381324.34596735.398301128.70329568.1661110.000000--第2次迭代时雅可比矩阵---88.468596 50.770062 15.630499 4.956802 0.000000 0.000000 0.000000 0.000000 0.000000 8.760031 8.351203 0.000000 0.000000 -17.679024 20.962621 6.976678 3.695668 -0.000000 -0.000000 -0.000000 -0.000000 -0.00000053.574257 -70.163300 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 8.844924 -0.000000 1.871649 0.000000 12.117328 -29.191523 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.00000015.630499 -0.000000 -85.475262 45.044595 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 24.800168 -6.976678 -0.000000 3.136302 10.112981 -0.000000 -0.000000 -0.000000 -0.000000 -0.0000004.956802 -0.000000 45.044595 -111.557696 48.101358 0.000000 0.000000 0.000000 13.454941 -0.000000 -0.000000 0.000000 -0.000000 -3.695668 -0.000000 -10.112981 -24.720611 28.512426-0.000000 -0.000000 -0.000000 12.548251-0.000000 -0.000000 -0.000000 54.962689 -73.761205 18.798516 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 -0.000000 -0.000000 -0.000000 10.286004 -30.269750 19.866386 -0.000000 -0.000000 -0.000000-0.000000 -0.000000 -0.000000 -0.000000 27.350216 -42.131939 0.000000 0.000000 -0.000000 -0.000000 -0.000000 14.781722 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.152201 -35.701334 -0.000000 -0.000000 -0.000000-0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -36.269587 20.414751 -0.000000 -0.000000 -0.000000 15.854836 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -43.168982 21.053877 -0.000000-0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 18.912486 -66.242459 18.617657 -0.000000 -0.000000 28.712316 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 22.413069 -72.744607 0.293713-0.000000 -0.000000 -0.000000 18.246829 -0.000000 0.000000 0.000000 11.613550 -29.860379 -0.000000 -0.000000 0.000000 -0.000000 -0.000000 -0.000000 -0.000000 2.355263 -0.000000 -0.000000 -0.000000 14.554342 -14.8093936.904762 6.537084 0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 -35.851861 4.723753 -0.000000 0.000000 5.396002 6.042147 -0.000000 -0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.0000007.404987 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 5.183063 -12.588050 -0.000000 0.000000 4.294174 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 0.000000 -0.000000 -0.000000-0.000000 1.871649 -0.000000 -0.000000 -0.000000 18.914409 16.625167 31.272979 -0.000000 -0.000000 -0.000000 -68.684204 -0.000000 -0.000000 -10.345614 -0.000000 -0.000000 -0.000000 3.718215 7.001258 12.708639 -0.000000-0.000000 -0.000000 24.800168 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 -24.800168 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.00000033.280772 -20.962621 -6.976678 -3.695668 0.000000 0.000000 0.000000 0.000000 0.000000 0.233337 -1.879142 0.000000 0.000000 -97.165875 50.770062 15.630499 4.956802 0.000000 0.000000 0.000000 0.000000 0.000000-12.117328 17.294580 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.004159 0.000000 -10.345614 0.000000 53.574257 -99.357152 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.0000006.976678 0.000000 3.136302 -10.112981 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 15.630499 -0.000000 -121.215950 45.044595 0.000000 0.000000 0.000000 0.000000 0.0000003.695668 0.000000 10.112981 27.252028 -28.512426 0.000000 0.000000 0.000000 -12.548251 0.000000 0.000000 0.000000 0.0000004.956802 -0.000000 45.044595 -125.395530 48.101358 0.000000 0.000000 0.000000 13.4549410.000000 0.000000 0.000000 -10.286004 30.152390 -19.866386 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 -0.000000 54.962689 -80.543606 18.798516 0.000000 0.000000 0.0000000.000000 0.000000 0.000000 0.000000 0.152201 12.220377 0.000000 0.000000 0.000000 0.000000 0.000000 -12.372578 0.000000 -0.000000 -0.000000 -0.000000 -0.000000 27.350216 -64.020523 0.000000 0.000000 -0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 29.658409 -21.053877 0.000000 0.000000 0.000000 -8.604532 0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000-62.187049 20.414751 -0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -22.413069 40.458166 -0.293713 0.000000 0.000000 -17.751384 0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 18.912486 -113.363276 18.6176570.000000 0.000000 0.000000 -2.355263 0.000000 0.000000 0.000000 -14.554342 16.909605 0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 -0.000000 18.246829 -0.000000 0.000000 0.000000 11.613550 -36.327356--第2次迭代时dPQ的误差--7.322874e+000-6.024472e+0003.136302e+0009.707086e-001-1.486800e-001-1.177548e+001-6.816287e+000-1.627822e+0019.011063e-0011.442203e+0003.766432e-001-6.045480e+0000.000000e+000-4.309639e+000-1.461293e+001-1.787034e+001-7.084917e+000-3.449200e+000-1.096229e+001-1.297473e+001-2.361841e+001-3.283488e+000--第2次迭代时dx(误差)--7.463082e-0017.954075e-0011.282017e+0001.553903e+0001.801690e+0002.314056e+0002.532583e+0002.510944e+0002.031354e+0002.758084e-0014.877708e-0012.596165e+0001.282017e+000-1.022801e-001-1.281804e-0012.457598e-002-1.510861e-001-1.422363e-0015.472778e-0022.816086e-0012.267577e-001-7.420019e-002第2次迭代后节点电压(仅PQ节点)1.7160411.7583862.2427312.5425962.6107482.3625102.1634712.2432872.673161第2次迭代后相角(角度)1.3521903.384402-5.288058-8.212275-12.380604-18.286594-19.228934-19.968564-17.7441768.5433117.451092-20.045989-5.2880580.000000--第3次迭代时雅可比矩阵---107.449063 64.339505 18.280415 7.735891 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 8.723401 8.369851 -0.000000 -0.000000 -33.906343 22.938392 -2.128130 -1.303493 -0.000000 -0.000000 -0.000000 -0.000000 -0.00000065.803521 -93.903508 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 9.224175 0.000000 10.883156 -0.000000 18.320439 -40.148863 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.00000018.280415 0.000000 -94.238507 51.767412 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 24.190680 2.128130 -0.000000 -0.516231 -2.644362 -0.000000 -0.000000 -0.000000 -0.000000 -0.0000007.735891 0.000000 51.767412 -151.916589 70.507019 -0.000000 -0.000000 -0.000000 21.906266 0.000000 0.000000 -0.000000 0.000000 1.303493 -0.000000 2.644362 -37.947834 20.832273 -0.000000 -0.000000 -0.000000 6.1357450.000000 0.000000 0.000000 66.741532 -94.948096 28.206563 -0.000000 -0.000000 -0.0000000.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 -0.000000 30.834903 -39.252891 8.745203 -0.000000 -0.000000 -0.0000000.000000 0.000000 0.000000 0.000000 25.819129 -42.495346 -0.000000 -0.000000 0.000000 0.000000 0.000000 16.676217 0.000000 -0.000000 -0.000000 -0.000000 -0.000000 14.333924 -21.142642 -0.000000 -0.000000 -0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -22.844227 11.084502 0.000000 0.000000 0.000000 11.759725 0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -20.202132 11.937858 -0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 10.772631 -47.668940 13.606970 0.000000 0.000000 23.289338 0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 12.220026 -36.325917 7.3518450.000000 0.000000 0.000000 18.700765 0.000000 -0.000000 -0.000000 14.136240 -32.837005 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 -0.000000 12.954276 -0.000000 -0.000000 -0.000000 6.274231 -17.3722379.480361 9.786321 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -41.877577 5.068935 -0.000000 -0.000000 1.851316 2.255032 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.0000009.101262 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 5.023265 -14.124527 -0.000000 -0.000000 2.489221 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.0000000.000000 10.883156 0.000000 0.000000 0.000000 16.194071 11.599663 23.257398 0.000000 0.000000 0.000000 -61.934289 0.000000 -0.000000 4.716418 -0.000000 -0.000000 -0.000000 8.353052 5.778354 11.849469 -0.0000000.000000 0.000000 24.190680 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 -24.190680 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.00000028.010895 -22.938392 2.128130 1.303493 0.000000 0.000000 0.000000 0.000000 0.000000 -4.148120 -4.356006 0.000000 0.000000 -118.100775 64.339505 18.280415 7.735891 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000-18.320439 19.018229 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -3.971376 0.000000 4.716418 0.000000 65.803521 -121.860594 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000-2.128130 0.000000 -0.516231 2.644362 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 18.280415 0.000000 -102.418265 51.767412 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000-1.303493 0.000000 -2.644362 30.915873 -20.832273 0.000000 0.000000 0.000000 -6.135745 0.000000 0.000000 0.000000 0.000000 7.735891 0.000000 51.767412 -162.046258 70.507019 -0.000000 -0.000000 -0.000000 21.9062660.000000 0.000000 0.000000 -30.834903 39.580107 -8.745203 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 66.741532 -106.373981 28.206563 -0.000000 -0.000000 -0.0000000.000000 0.000000 0.000000 0.000000 -14.333924 21.677303 0.000000 0.000000 0.000000 0.000000 0.000000 -7.343379 0.000000 0.000000 0.000000 0.000000 0.000000 25.819129 -52.356079 -0.000000 -0.000000 0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 17.383077 -11.937858 0.000000 0.000000 0.000000 -5.445220 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -27.967881 11.084502 0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -12.220026 31.358440 -7.351845 0.0000000.000000 -11.786569 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 10.772631 -59.717994 13.6069700.000000 0.000000 0.000000 -12.954276 0.000000 0.000000 0.000000 -6.274231 19.228507 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 18.700765 0.000000 -0.000000 -0.000000 14.136240 -43.537423--第3次迭代时dPQ的误差---3.425724e+000-1.064132e+001-5.162313e-001-3.810980e+0007.360782e-0022.323305e-001-1.470528e+000-2.618738e+0007.791352e-001-2.042190e+000-3.425627e-0011.156931e+0010.000000e+000-5.286856e+000-1.399454e+001-4.089879e+000-5.230835e+000-5.770943e+000-4.948367e+000-2.577827e+000-6.082527e+000-5.400209e+000--第3次迭代时dx(误差)--9.674289e-0021.172732e-001-7.870880e-003-4.668014e-002-8.644519e-002-1.584698e-001-1.940180e-001-1.949970e-001-1.391568e-0011.392652e-0011.715593e-001-2.170152e-001-7.870880e-0032.010999e-0012.142528e-0012.135755e-0012.788091e-0012.825328e-0012.223360e-0011.748684e-0011.963064e-0012.880642e-001第3次迭代后节点电压(仅PQ节点)1.3709451.3816471.7637381.8336971.8731261.8372391.7851491.8029151.903119第3次迭代后相角(角度)-4.190769-3.334859-4.837089-5.537700-7.427659-9.206943-8.112522-8.796060-9.7710800.564000-2.378530-7.611937-4.8370890.000000--第4次迭代时雅可比矩阵---70.787576 40.675150 11.561949 4.518688 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 7.103700 6.928088 -0.000000 -0.000000 -20.202620 13.567039 -0.130429 -0.106247 -0.000000 -0.000000 -0.000000 -0.000000 -0.00000041.062278 -63.775374 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 7.314738 0.000000 9.293937 -0.000000 12.345919 -20.057751 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.00000011.561949 0.000000 -59.982694 29.396602 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 19.024142 0.130429 -0.000000 -0.229050 -0.359479 -0.000000 -0.000000 -0.000000 -0.000000 -0.0000004.518688 0.000000 29.396602 -80.849041 36.025140 -0.000000 -0.000000 -0.000000 10.908611 0.000000 0.000000 -0.000000 0.000000 0.106247 -0.000000 0.359479 -18.954787 12.221085 -0.000000 -0.000000 -0.000000 4.1755430.000000 0.000000 0.000000 35.141107 -50.486980 15.345872 -0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 -0.000000 14.569427 -20.011329 5.999251 -0.000000 -0.000000 -0.000000。

前言电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。

在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。

可靠性和经济性。

此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。

所以潮流计算是研究电力系统的一种很重要和很基础的计算。

随着科学技术的发展，电力系统变得越来越复杂，电气工程师掌握一种好的能对电力系统进行仿真的软件是学习和研究的需要。

文章简要介绍了MATLAB发展历史、组成和强大的功能，并用简单例子分别就编程和仿真两方面分析了MATIAB软件在电力系统研究中的具体应用。

采取等效电路法，能对特殊、复杂地电力系统进行高效仿真研究，因此，掌握编程和仿真是学好MATLAB的基础。

与众多专门的电力系统仿真软件相比，MATLAB软件具有易学、功能强大和开放性好，是电力系统仿真研究的有力工具。

1电力系统的基本概念电力系统：发电机把机械能转化为电能，电能经变压器和电力线路输送并分配到用户，在那里经电动机、电炉和电灯等设备又将电能转化为机械能、热能和光能等。

这些生产、变换、输送、分配、消费电能的发电机、变压器、变换器、电力线路及各种用电设备等联系在一起组成的统一整体称为电力系统。

电力网：电力系统中除发电机和用电设备外的部分。

动力系统：电力系统和“动力部分”的总和。

2潮流计算2.1潮流计算概述与发展电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。

利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。

在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。

对潮流计算的要求可以归纳为下面几点：1）计算方法的可靠性或收敛性；2）对计算机内存量的要求；3）计算速度；4）计算的方便性和灵活性。

％简单潮流计算的牛顿拉夫逊程序,相关的原始数据数据数据输入格式如下：％B1是支路参数矩阵，第一列和第二列是节点编号。

节点编号由小到大编写%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点%编号,将变压器的串联阻抗置于低压侧处理。

％第三列为支路的串列阻抗参数。

%第四列为支路的对地导纳参数。

%第五列为含变压器支路的变压器的变比％第六列为变压器是否含有变压器的参数,其中“1"为含有变压器，%“0”为不含有变压器。

％B2为节点参数矩阵,其中第一列为节点注入发电功率参数；第二列为节点负荷功率参数;第三列为节点电压参数；第六列为节点类型参数,其中“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。

％X为节点号和对地参数矩阵.其中第一列为节点编号，第二列为节点对地％参数。

n=input(’请输入节点数:n=');n1=input（’请输入支路数:n1=');isb=input（’请输入平衡节点号:isb='）；pr=input('请输入误差精度：pr=');B1=input（’请输入支路参数:B1=’);B2=input(’请输入节点参数：B2=’)；X=input('节点号和对地参数:X=');Y=zeros（n);Times=1; %置迭代次数为初始值％创建节点导纳矩阵for i=1:n1if B1(i,6)==0 %不含变压器的支路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(p，p)=Y（p，p)+1/B1(i，3)+0。

5＊B1（i,4）;Y(q，q）=Y(q，q)+1/B1(i，3)+0。

5＊B1(i，4)；else %含有变压器的支路p=B1(i,1)；q=B1(i，2);Y(p，q）=Y（p,q）—1/（B1(i,3）＊B1(i，5));Y（q，p）=Y（p,q)；Y（p，p)=Y（p,p）+1/B1（i，3)；Y（q,q)=Y(q，q）+1/(B1(i，5）^2*B1（i，3））;endendYOrgS=zeros(2＊n—2，1);DetaS=zeros(2＊n—2,1); %将OrgS、DetaS初始化％创建OrgS，用于存储初始功率参数h=0；j=0;for i=1：n ％对PQ节点的处理if i～=isb＆B2（i,6）==2h=h+1;for j=1：nOrgS(2*h-1，1）=OrgS(2*h—1，1）+real（B2（i，3））*（real（Y （i，j）)＊real(B2(j，3）)-imag（Y(i，j））*imag(B2(j,3）)）+imag(B2(i,3))＊(real（Y(i,j））*imag(B2(j，3）)+imag（Y(i，j)）*real(B2(j,3））)；OrgS（2＊h,1）=OrgS(2＊h,1）+imag（B2（i，3））＊（real(Y（i，j）)＊real(B2（j，3））—imag(Y(i，j））＊imag(B2（j,3)）)-real（B2（i,3））＊（real(Y(i，j)）＊imag(B2(j,3）)+imag（Y(i，j))＊real（B2（j，3））);endendendfor i=1:n ％对PV节点的处理，注意这时不可再将h初始化为0 if i～=isb&B2（i，6)==3h=h+1;for j=1:nOrgS（2*h—1,1）=OrgS(2*h—1,1)+real(B2（i，3）)＊(real(Y（i，j）)＊real（B2（j，3）)-imag(Y（i,j))＊imag(B2(j，3）））+imag(B2(i,3）)*(real(Y（i，j)）*imag（B2(j，3））+imag（Y(i，j))＊real(B2（j，3））);OrgS（2＊h,1)=OrgS(2*h，1）+imag（B2（i，3)）*(real（Y(i，j））*real(B2（j,3）)-imag（Y(i，j)）＊imag(B2（j,3)))—real（B2(i,3））*(real（Y（i,j)）*imag(B2(j,3)）+imag（Y（i，j))＊real(B2(j，3)))；endendendOrgS％创建PVU 用于存储PV节点的初始电压PVU=zeros(n-h—1，1）；t=0;for i=1:nif B2（i,6)==3t=t+1；PVU（t,1)=B2(i,3）;endendPVU%创建DetaS，用于存储有功功率、无功功率和电压幅值的不平衡量h=0；for i=1：n ％对PQ节点的处理if i～=isb&B2(i,6)==2h=h+1;DetaS(2*h—1，1)=real（B2（i,2））-OrgS(2*h—1，1）;DetaS(2＊h,1）=imag(B2(i，2)）—OrgS（2*h，1)；endendt=0；for i=1:n ％对PV节点的处理，注意这时不可再将h初始化为0 if i～=isb&B2（i,6）==3h=h+1;t=t+1;DetaS(2*h-1，1）=real(B2(i,2））-OrgS(2*h-1，1);DetaS（2＊h,1)=real（PVU（t，1））^2+imag（PVU（t,1）)^2-real(B2(i,3）)^2—imag(B2（i,3）)^2；endendDetaS%创建I,用于存储节点电流参数i=zeros(n—1，1）；h=0；for i=1：nif i～=isbh=h+1;I（h，1)=(OrgS(2＊h—1，1)-OrgS(2＊h，1)＊sqrt（-1))/conj（B2(i,3））;endendI％创建Jacbi（雅可比矩阵)Jacbi=zeros(2*n—2);h=0;k=0；for i=1:n ％对PQ节点的处理if B2（i,6）==2h=h+1;for j=1:nif j～=isbk=k+1;if i==j ％对角元素的处理Jacbi(2＊h-1,2＊k—1)=—imag(Y（i，j））*real(B2(i,3))+real（Y（i,j）)*imag（B2(i，3））+imag(I（h，1));Jacbi(2＊h-1,2＊k）=real（Y(i,j）)＊real(B2（i，3)）+imag(Y(i,j)）＊imag（B2(i，3））+real(I（h，1））；Jacbi（2*h，2＊k-1)=-Jacbi(2*h-1，2＊k）+2*real （I(h,1));Jacbi(2＊h,2*k)=Jacbi（2＊h-1,2*k-1)—2*imag （I（h,1）)；else %非对角元素的处理Jacbi（2＊h—1,2*k-1）=—imag(Y(i,j））＊real(B2(i,3）)+real(Y（i,j）)＊imag（B2(i，3)）;Jacbi(2*h—1,2＊k)=real(Y（i，j))*real(B2（i,3))+imag(Y(i，j)）＊imag(B2（i，3））；Jacbi（2＊h，2*k-1)=—Jacbi(2＊h—1，2*k);Jacbi(2＊h,2＊k)=Jacbi(2＊h—1,2＊k—1);endif k==（n-1） %将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0；endendendendendk=0;for i=1：n %对PV节点的处理if B2(i，6）==3h=h+1;for j=1：nif j～=isbk=k+1;if i==j ％对角元素的处理Jacbi（2＊h—1,2＊k-1）=—imag(Y(i，j））*real(B2(i，3）)+real（Y（i,j））*imag(B2（i，3)）+imag(I(h，1））；Jacbi(2＊h—1,2*k）=real（Y(i，j）)＊real(B2（i,3)）+imag(Y（i,j））*imag(B2(i，3)）+real(I(h,1））;Jacbi（2*h,2*k-1)=2＊imag(B2（i,3）)；Jacbi（2＊h，2*k)=2*real（B2（i,3）)；else %非对角元素的处理Jacbi（2*h-1，2*k—1)=-imag（Y(i,j））＊real(B2（i,3））+real(Y(i，j))＊imag(B2（i，3)）;Jacbi（2*h—1,2*k)=real(Y（i,j))＊real(B2（i，3))+imag(Y(i,j））*imag(B2（i,3））;Jacbi（2＊h，2*k—1)=0;Jacbi（2＊h,2＊k）=0;endif k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0；endendendendendJacbi％求解修正方程，获取节点电压的不平衡量DetaU=zeros(2＊n—2,1)；DetaU=inv（Jacbi)＊DetaS；DetaU％修正节点电压j=0；for i=1:n ％对PQ节点处理if B2（i,6)==2j=j+1;B2(i，3）=B2(i，3）+DetaU(2＊j,1)+DetaU(2*j—1,1）*sqrt（—1);endendfor i=1:n %对PV节点的处理if B2(i,6)==3j=j+1；B2(i，3)=B2(i,3)+DetaU（2＊j,1）+DetaU（2＊j-1，1)*sqrt(-1）;endendB2％开始循环＊**＊*******＊*******＊*＊*＊**＊＊**＊**＊＊＊＊****＊＊＊＊**＊＊***＊******＊＊＊*****＊*＊while abs（max(DetaU))>prOrgS=zeros（2＊n-2,1); ％初始功率参数在迭代过程中是不累加的，所以在这里必须将其初始化为零矩阵h=0;j=0;for i=1:nif i~=isb＆B2（i,6）==2h=h+1;for j=1：nOrgS(2*h—1,1)=OrgS(2＊h-1，1)+real（B2(i，3））＊（real(Y(i，j)）*real(B2(j，3）)—imag（Y（i,j))*imag（B2(j,3)）)+imag(B2(i,3))*（real（Y（i，j)）*imag（B2(j，3）)+imag(Y（i,j))＊real（B2(j,3）））；OrgS（2＊h,1）=OrgS（2＊h,1）+imag（B2(i,3))*（real（Y（i，j））＊real（B2(j，3））—imag(Y(i，j)）*imag(B2(j,3））)-real(B2（i，3）)*（real(Y（i,j））*imag(B2（j,3)）+imag(Y（i,j）)＊real（B2(j，3）)）；endendendfor i=1:nif i~=isb&B2(i,6）==3h=h+1；for j=1：nOrgS(2＊h-1,1)=OrgS(2*h-1,1)+real（B2（i,3））＊（real(Y(i，j））*real(B2(j，3))-imag（Y（i,j））＊imag（B2(j,3)))+imag(B2(i，3））*（real（Y（i，j））*imag(B2(j，3))+imag(Y（i，j)）*real(B2(j，3）））;OrgS(2＊h,1）=OrgS（2*h,1)+imag(B2（i，3）)＊（real（Y（i,j))*real （B2（j，3））-imag(Y（i，j）)＊imag(B2(j，3）））—real（B2（i,3）)＊(real(Y(i,j）)＊imag(B2(j，3)）+imag（Y(i，j）)*real（B2（j，3)））；endendendOrgS%创建DetaSh=0；for i=1:nif i~=isb&B2（i，6)==2h=h+1；DetaS(2＊h-1,1）=real（B2（i，2））—OrgS(2＊h—1，1）；DetaS（2＊h，1)=imag(B2（i,2））—OrgS(2＊h,1）;endendt=0;for i=1:nif i～=isb&B2（i,6)==3h=h+1;t=t+1；DetaS(2*h—1,1）=real(B2（i，2）)—OrgS（2＊h—1,1)；DetaS（2*h,1)=real（PVU（t，1）)^2+imag（PVU(t，1)）^2-real（B2（i,3）)^2-imag(B2(i，3）)^2；endendDetaS%创建Ii=zeros(n-1,1)；h=0;for i=1:nif i~=isbh=h+1;I(h,1）=（OrgS（2*h—1，1）-OrgS（2＊h,1)*sqrt（—1)）/conj（B2(i，3)）；endendI%创建JacbiJacbi=zeros(2*n—2）;h=0;k=0;for i=1:nif B2(i，6)==2h=h+1；for j=1：nif j～=isbk=k+1；if i==jJacbi（2＊h—1，2＊k—1）=-imag（Y(i，j）)＊real（B2(i，3）)+real(Y(i，j))＊imag(B2（i,3）)+imag(I(h，1)）;Jacbi(2*h-1,2*k)=real(Y（i,j））＊real（B2（i,3））+imag（Y(i,j））*imag(B2（i，3))+real(I(h,1))；Jacbi（2＊h,2＊k—1)=-Jacbi(2*h-1，2*k）+2*real(I(h，1））;Jacbi(2＊h，2*k)=Jacbi（2＊h-1，2＊k—1)—2*imag（I(h，1）)；elseJacbi（2*h—1，2*k—1)=-imag(Y（i，j))＊real(B2（i，3）)+real(Y(i,j）)＊imag（B2（i，3)）；Jacbi（2＊h—1，2＊k)=real（Y(i，j)）＊real （B2(i,3））+imag（Y（i，j）)＊imag(B2(i,3))；Jacbi（2＊h，2＊k—1）=-Jacbi（2＊h-1，2*k）；Jacbi（2＊h,2＊k)=Jacbi(2＊h—1，2*k-1)；endif k==（n-1)k=0;endendendendendk=0；for i=1：nif B2(i，6）==3h=h+1；for j=1:nif j～=isbk=k+1;if i==jJacbi(2*h—1,2*k-1）=—imag(Y（i，j))*real(B2(i，3）)+real（Y(i,j))＊imag（B2（i，3））+imag(I(h,1）);Jacbi（2*h-1，2*k）=real（Y（i，j)）*real(B2（i,3)）+imag(Y(i，j))＊imag(B2（i,3))+real（I（h,1))；Jacbi(2*h,2＊k-1)=2＊imag(B2（i，3）);Jacbi（2＊h，2*k)=2＊real(B2(i,3））；elseJacbi(2*h-1,2＊k—1)=-imag（Y（i，j））＊real （B2(i，3））+real（Y（i，j）)＊imag（B2(i，3）)；Jacbi（2＊h-1,2*k)=real(Y(i,j））＊real（B2(i,3)）+imag（Y(i,j）)*imag（B2（i，3））;Jacbi(2＊h，2*k-1）=0；Jacbi（2＊h，2＊k）=0;endif k==（n—1）k=0；endendendendendJacbiDetaU=zeros（2*n-2，1）；DetaU=inv(Jacbi）＊DetaS；DetaU％修正节点电压j=0;for i=1:nif B2(i,6）==2j=j+1;B2（i,3)=B2(i，3)+DetaU(2*j，1)+DetaU(2＊j—1,1)＊sqrt(-1）;endendfor i=1：nif B2（i，6)==3j=j+1；B2(i,3)=B2(i,3)+DetaU（2*j,1）+DetaU（2＊j-1,1）*sqrt(—1)；endendB2Times=Times+1; %迭代次数加1endTimes一个原始数据的例子节点数 5支路数 5平衡节点编号 5精度pr 0。

MATLAB极坐标求解潮流计算极坐标是一种坐标系，用于描述平面上的点，其坐标由极径和极角组成。

在电力系统中，潮流计算是一种用于估计电力系统各节点电压幅值和相角的方法，以确定电力系统的负荷和发电机出力分布情况。

1.建立潮流计算的数学模型：潮流计算问题可以通过组成节点功率方程和支路潮流方程来描述。

节点功率方程描述了节点电压与注入和抽取的功率之间的关系，而支路潮流方程描述了支路接线和节点之间的相对功率转移。

2.收集输入数据：潮流计算需要收集电力系统的拓扑结构信息（节点和支路的连接关系），以及节点功率信息（负荷和发电机出力）。

这些数据可以从现场测量、电力系统数据库或网络模拟软件中获取。

3.构建节点导纳矩阵：节点导纳矩阵描述了电力系统各节点之间的电气连接关系。

它可以通过支路导纳矩阵和节点导纳矩阵的顺序乘法来计算。

4.进行潮流计算：利用节点功率方程和支路潮流方程，可以建立一个非线性方程组，其未知数是各节点的电压幅值和相角。

可以使用迭代的方法（如牛顿-拉夫逊法）或者线性化的方法（如高斯-赛德尔法）来求解这个方程组。

MATLAB提供了丰富的数值计算函数和算法，可以用来求解非线性方程组。

5.计算结果分析：潮流计算会得到电力系统各节点的电压幅值和相角，以及支路上的功率。

可以使用MATLAB绘制极坐标图来显示电力系统节点的电压幅值和相角分布情况。

此外，还可以对潮流计算结果进行故障分析、潮流分布优化等。

总结来说，使用MATLAB进行极坐标求解潮流计算可以帮助电力系统运营人员了解电力系统的电压幅值和相角分布情况，进而优化电力系统的运行和规划。

基于MATLAB进行潮流计算学生：王仕龙2011148213指导老师：李咸善摘要：电力系统潮流计算方法有两类，即手算潮流和计算机潮流计算。

手算潮流主要借助于形成简化的等值电路来实现，这种方法尤其适用于规模不大的辐射型电力潮流计算。

计算机潮流计算的实现有两种途径：其一是编程实现网络方程的迭代求解；其二是借助与电力系统分析仿真软件，搭建系统模型来完成潮流计算。

MATLAB具有强大的矩阵运算功能，同时其具有电力系统仿真平台也为直观地实现潮流计算提供了更便捷的手段[1]。

本文是基于MATLAB软件，采用极坐标形式牛顿─拉夫逊法进行潮流计算，为其他形式的潮流计算有借鉴的作用。

关键词：电力系统；计算机潮流计算；MATLAB ；牛顿─拉夫逊法Abstract:The power flow calculation method has two kinds,which are the hand calculation of tidal current and computer power flow calculation.Hand calculation tidal current is mainly realized by means of the formation of simplified equivalent circuit.This method is especially suitable for small scale radiation power flow calculation.There are two ways to realize the computer power flow calculation.The first one is through the programming iteration for solving network equation,the second one is with the help of analysis of power system simulation software to build the system model complete the power flow calculation.The software of MATLAB has strong matrix function,.At the same time,It’s power system simulation platform provides a more convenient means to realize power flow calculation intuitively[1].This paper is based on the software of MATLAB to calculate the power flow calculation by adopting the form of Newton-Raphson method of power flow calculation of polar coordinates.And it can be the role of reference of other forms of power flow calculation.Key words: power system computer; power flow calculation;MATLAB;Newton-Raphson1.计算原理电力系统潮流是指系统中所有运行参数的总体，包括各个母线电压的大小和相位，各个发电机和负荷的功率及电流，以及各个变压器和线路等元件所通过的功率,电流和其中的损耗。