用Matlab计算潮流计算
MATpower潮流计算软件的使用
MATpower软件的使用一、MATpower软件的使用方法在MATLAB软件中的命令窗口输入runpf（‘程序名’）就可以通过MATpower已经编好的程序计算潮流，而函数runpf的参数是相应需计算潮流的数据文件。
数据文件主要用来定义和返回一下四个变量。
1、baseMVA baseMVA是一个标量，用来设置基准容量。
对于计算中采用有名值，可以根据实际情况设置，如设置100MV A；对于计算中采用标幺值，一般设置为12、bus bus变量是一个矩阵，用来设置电网中各节点参数，该矩阵内的参数如下：%% bus data%bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin其中，第1列参数即bus_i用来设置母线编号，范围为1~29997；第2列参数type用来设置母线类型，1为PQ节点，2为PV节点，3为平衡节点；第3列参数Pd用来设置母线注入负荷的有功功率；第4列参数Qd用来设置母线注入负荷的无功功率；第5列参数Gs用来设置与母线并联的电导；第6列参数Bs用来设置与母线并联的电纳；第7列参数area用来设置电网断面号，可设置范围为1~100，一般设置为1；第8列参数Vm用来设置母线电压的幅值初值；第9列参数Va用来设置母线电压的相角初值；第10列参数baseKV用来设置该母线的基准电压；第11列参数zone用来设置省耗分区号，可设置范围为1~999，一般设置为1；第12列参数Vmax用来设置工作时母线电压最高幅值；第13列参数Vmin用来设置工作时母线电压最低幅值。
3、gen gen变量是一个矩阵，用来设置接入电网的发电机参数，该矩阵的参数如下：%% generator data%bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin 其中，第1列参数bus用来设置接入发电机的母线编号；第2列参数Pg用来设置接入发电机的有功功率，注意功率输入的是有名值；第3列参数Qg用来设置接入发电机的无功功率；第4列参数Qmax用来设置接入发电机的无功功率的最大允许值；第5列参数Qmin用来设置接入发电机的无功功率的最小允许值；第6列Vg用来设置接入发电机的工作电压，注意输入的是标幺值；第7列mBase用来设置接入发电机的功率基准；第8列status用来设置发电机的工作状态，1表示投入运行，2表示投出运行；第9列Pmax用来设置接入发电机的无功功率的最大允许值；第10列参数Pmin用来设置接入发电机的无功功率的最小允许值。
潮流计算(matlab)实例计算
潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。
2.在给定的电力网络上画出等值电路图。
3.运用计算机进行潮流计算。
4.编写设计说明书。
一、设计原理1．牛顿-拉夫逊原理牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。
牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。
电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。
为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。
牛顿—拉夫逊迭代法的一般步骤：（1）形成各节点导纳矩阵Y。
（2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。
（3）计算各个节点的功率不平衡量。
（4）根据收敛条件判断是否满足，若不满足则向下进行。
（5）计算雅可比矩阵中的各元素。
（6）修正方程式个节点电压（7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。
（8）计算平衡节点输出功率和各线路功率2．网络节点的优化１）静态地按最少出线支路数编号这种方法由称为静态优化法。
matlab潮流计算仿真方法
matlab潮流计算仿真方法
MATLAB 是一种强大的编程语言和环境，可用于执行各种仿真和计算任务，包括电力系统潮流计算。
以下是一个简单的 MATLAB 潮流计算仿真方法的
示例：
1. 定义系统参数：首先，你需要定义电力系统的参数，如发电机、负荷、变压器等。
这些参数通常包括额定电压、额定功率、电抗、电阻等。
2. 建立系统模型：使用这些参数，你可以在 MATLAB 中建立电力系统的模型。
这通常涉及到定义节点和支路，以及为它们分配相应的参数。
3. 编写潮流计算函数：接下来，你需要编写一个用于执行潮流计算的函数。
这个函数应该能够接收系统的模型和参数，并返回计算出的潮流结果，如电压、电流、功率等。
4. 运行仿真：最后，你可以运行仿真并调用你编写的潮流计算函数。
这将返回计算出的潮流结果，你可以使用这些结果进行进一步的分析或可视化。
这只是一个简单的示例，实际上在编写 MATLAB 潮流计算仿真方法时可能
需要考虑更多因素，例如系统的约束条件、初始条件、迭代算法的收敛性等。
如果你需要具体的 MATLAB 代码示例或更详细的指导，我建议你查阅MATLAB 的官方文档或相关的教程和文献。
电力系统潮流计算matlab程序
电力系统潮流计算matlab程序电力系统潮流计算是电力系统运行和规划中的重要环节，它用于计算电力系统中各节点的电压、功率和电流等参数。
随着电力系统规模的不断扩大和复杂性的增加，传统的手工计算方法已经无法满足需求，因此，利用计算机编程进行潮流计算成为了一种必要的选择。
Matlab是一种功能强大的科学计算软件，它提供了丰富的数学函数和工具箱，可以方便地进行电力系统潮流计算。
下面我将介绍一下如何使用Matlab编写电力系统潮流计算程序。
首先，我们需要建立电力系统的节点模型。
节点模型是电力系统中各节点的电压、功率和电流等参数的数学表示。
在Matlab中，我们可以使用矩阵来表示节点模型。
假设电力系统有n个节点，我们可以定义一个n×n的复数矩阵Y来表示节点之间的导纳关系，其中Y(i,j)表示节点i和节点j之间的导纳。
同时，我们还需要定义一个n×1的复数向量V来表示各节点的电压，其中V(i)表示节点i的电压。
接下来，我们需要编写潮流计算的主程序。
主程序的主要功能是根据节点模型和潮流计算算法，计算出各节点的电压、功率和电流等参数。
在Matlab中，我们可以使用循环语句和矩阵运算来实现潮流计算。
具体的计算过程可以参考电力系统潮流计算的算法。
在编写主程序之前，我们还需要定义一些输入参数，如电力系统的节点数、发电机节点和负荷节点等。
这些参数可以通过用户输入或者读取文件的方式获取。
同时，我们还需要定义一些输出参数，如各节点的电压、功率和电流等。
这些参数可以通过矩阵运算和循环语句计算得到，并输出到文件或者显示在屏幕上。
最后，我们需要进行程序的测试和调试。
可以通过输入一些测试数据，运行程序并检查输出结果是否正确。
如果发现程序有错误或者结果不准确，可以通过调试工具和打印调试信息的方式进行调试。
总之，利用Matlab编写电力系统潮流计算程序可以提高计算效率和准确性，为电力系统的运行和规划提供有力的支持。
当然，编写一个完整的潮流计算程序需要考虑很多细节和特殊情况，这需要有一定的电力系统和编程知识。
基于MATLAB进行潮流计算
基于MATLAB进行潮流计算本文介绍了基于MATLAB软件的潮流计算方法。
电力系统潮流计算方法分为手算潮流和计算机潮流计算两类。
手算潮流主要适用于规模较小的辐射型电力潮流计算，而计算机潮流计算有两种途径：编程实现网络方程的迭代求解和借助电力系统分析仿真软件搭建系统模型完成潮流计算。
MATLAB具有强大的矩阵运算功能和电力系统仿真平台，可以为实现潮流计算提供更便捷的手段。
本文采用极坐标形式牛顿─拉夫逊法进行潮流计算，为其他形式的潮流计算提供借鉴。
Abstract: The power flow n method can be divided into two categories: hand n of tidal current and computer power flow XXX simplified equivalent circuits。
making it XXX: programming XXX ns。
or using power system XXX system model for power flow n。
MATLAB are has strong matrix ns and its power system XXX-Raphson method of power flow n in polar coordinates with MATLAB are。
and can serve as a reference for other forms of power flow n.1.电力系统中的牛顿法潮流计算是一种常用的电力系统分析方法。
该方法基于节点电压的相等条件和潮流方程的等式条件，通过迭代求解电压和相位的不平衡量，最终得到各节点的电压、相位和功率等参数。
2.牛顿法潮流计算的步骤包括输入系统原始数据、形成节点导纳矩阵、给定各节点电压初值、计算功率偏差向量、判断收敛条件、计算雅克比矩阵、解修正方程、计算节点电压和相位的修正值、迭代计算直至满足收敛条件、计算各节点功率等参数并输出计算结果。
基于MATLAB的电力系统潮流计算
%输出计算结果
disp('节点电压为：');
通过这个程序，我们可以得到电力系统的节点电压向量。同样地，我们也可 以用节点电流法或迭代算法来求解潮流计算问题。
对于不同的算法，它们的优缺点也不尽相同。节点电压法具有计算量小、收 敛速度快等优点，但需要已知各节点的电压初始值。节点电流法相对于节点电压 法而言，收敛速度较慢，但不需要知道电压初始值。迭代算法具有普适性，可以 处理各种复杂的
基于MATLAB的电力系统潮流计算
目录
01 引言
03 Matlab工具
02 背景 04 潮流计算方法
05 结果分析
07 参考内容
目录
06 结论
引言
电力系统潮流计算是电力工程领域中非常重要的分析工具之一，用于研究电 力系统中电压、电流、功率等参数的分布和分配情况。准确地进行电力系统潮流 计算能够为电力系统的设计和运行提供重要的参考依据。本次演示将介绍使用 Matlab进行电力系统
2、利用Matlab的仿真功能，设置计算迭代的步长和算法类型等参数。
3、调用电力系统潮流计算函数， 开始计算并得到潮流结果。
4、对潮流结果进行分析和优化，为电力系统的设计和运行提供参考。
潮流计算方法
电力系统潮流计算的方法主要包括以下几个步骤：
1、网络拓扑分析：根据电力系统的结构，分析其网络拓扑关系，确定电力 系统的运行状态。
电力系统，但需要设定合适的迭代步长和初始值。
在未来研究中，我们可以进一步探索混合潮流计算方法，将不同的算法进行 组合，以获得更好的计算性能。此外，随着智能电网技术的发展，我们可以考虑 将潮流计算与优化、控制相结合，实现电力系统的智能化运行。
综上所述，基于MATLAB的电力系统潮流计算在电力工程领域具有广泛的应用 前景。通过深入研究和不断创新，我们可以进一步提高潮流计算的精度和效率， 为电力系统的稳定和经济运行提供更好的支持。
Matlab实现潮流计算程序
程序代码如下:111111.%读入数据clcclearfilename='123.txt';a=textread(filename)n=a(1,1);pinghengjd=a(1,2);phjddianya=a(1,3);jingdu=a(1,4);b=zeros(1,9);j1=0;[m1,n1]=size(a);for i1=1:m1if a(i1,1)==0j1=j1+1;b(j1)=i1;endendb;%矩阵分块a1=a(b(1)+1:b(2)-b(1)+1,1:n1);a2=a(b(2)+1:b(3)-1,1:n1);a3=a(b(3)+1:b(4)-1,1:n1);a4=a(b(4)+1:b(5)-1,1:n1);a5=a(b(5)+1:b(6)-1,1:n1);%设置初值vcz=1;dcz=0;kmax=20;k1=0;%求节点导纳矩阵a11=zeros(4,6);for i0=1:3for j0=1:6a11(i0,j0)=a1(i0,j0);a11(4,j0)=a2(1,j0);endenda11;linei=a11(1:4,2);linej=a11(1:4,3);liner=a11(1:4,4);linex=a11(1:4,5);lineb=a11(1:4,6);branchi=0;branchj=0;branchb=0;G=zeros(4,4);B=zeros(4,4);for k=1:4i2=linei(k,1);j2=linej(k,1);r=liner(k,1);x=linex(k,1);b=0;GIJ=r/(r*r+x*x);BIJ=-x/(r*r+x*x);if k>=4 & lineb(k)~=0k0=lineb(k);G(i2,j2)=-GIJ/k0;G(j2,i2)=G(i2,j2);B(i2,j2)=-BIJ/k0;B(j2,i2)=B(i2,j2);G(i2,i2)=G(i2,i2)+GIJ/k0/k0; B(i2,i2)=B(i2,i2)+BIJ/k0/k0;elseG(j2,i2)=-GIJ;G(i2,j2)=G(j2,i2);B(j2,i2)=-BIJ;B(i2,j2)=B(j2,i2);G(i2,i2)=G(i2,i2)+GIJ;b=lineb(k);B(i2,i2)=B(i2,i2)+BIJ+b;endG(j2,j2)=G(j2,j2)+GIJ;B(j2,j2)=B(j2,j2)+BIJ+b;endG;B;B=B.*i;Yf=G+BY=abs(Yf);alf=angle(Yf);%赋Jacobian矩阵参数P=zeros(n,1);Q=zeros(n,1);Pd=zeros(1,n);Qd=zeros(1,n);dP=zeros(1,n);dQ=zeros(1,n);PG=a4(:,3);PD=a4(:,5);QG=a4(:,4);QD=a4(:,6);i8=a4(:,2);for j8=1:length(i8)P(i8(j8))=PG(i8(j8))-PD(i8(j8));Q(i8(j8))=QG(i8(j8))-QD(i8(j8));enddelt=zeros(n,1);V=ones(n,1);V(3)=1.10;V(4)=1.05;ddelt=zeros(n,1);dV=zeros(n,1);A=zeros(2*n,2*n);B=zeros(2*n,1);Jacobian=Jaco(V,delt,n,Y,alf)%求取矩阵功率for j5=1:kmaxdisp(['第' int2str(j5) '次计算结果'])if k>=kmaxbreakendfor i10=1:4Pd(i10)=0;Qd(i10)=0;for j10=1:nPd(i10)=Pd(i10)+V(i10)*Y(i10,j10)*V(j10)*cos(d elt(i10)-delt(j10)-alf(i10,j10));Qd(i10)=Qd(i10)+V(i10)*Y(i10,j10)*V(j10)*sin(d elt(i10)-delt(j10)-alf(i10,j10));endendfor i4=1:3dP(i4)=P(i4)-Pd(i4);endfor j4=1:2dQ(j4)=Q(j4)-Qd(j4);endA=Jaco(V,delt,n,Y,alf)for i14=1:nB(i14*2-1)=-dP(i14);B(i14*2)=-dQ(i14);endif max(abs(B))>jingduX=A\B;for i16=1:nddelt(i16)=X(2*i16-1);dV(i16)=X(2*i16)*V(i16);endV=V+dVdelt=delt+ddeltelsebreakenddisp('----------------')end%流氓算法% for ii=1:2% V(ii)=V(ii)+dV(ii);% end% V222222.function A=Jaco(V,delt,n,Y,alf)%计算Jacobian矩阵for i7=1:nHd1(i7)=0;Jd1(i7)=0;for j7=1:nHd1(i7)=Hd1(i7)+V(i7)*Y(i7,j7)*V(j7)*sin(delt(i7)-delt(j7)-alf(i7,j7));Jd1(i7)=Jd1(i7)+V(i7)*Y(i7,j7)*V(j7)*cos(delt(i7)-delt(j7)-alf(i7,j7));endendfor i6=1:nfor j6=1:nif i6~=j6H(i6,j6)=-V(i6)*Y(i6,j6)*V(j6)*sin(delt(i6)-delt(j6)-alf(i6,j6));N(i6,j6)=-V(i6)*Y(i6,j6)*V(j6)*cos(delt(i6)-delt(j6)-alf(i6,j6));J(i6,j6)=-N(i6,j6);L(i6,j6)=H(i6,j6);elseH(i6,i6)=Hd1(i6)-V(i6)*Y(i6,i6)*V(i6)*sin(delt(i6)-delt(j6)-alf(i6,j6));J(i6,j6)=-Jd1(i6)+V(i6)*Y(i6,j6)*V(j6)*cos(delt(i6)-delt(j6)-alf(i6,j6));N(i6,j6)=-Jd1(i6)-V(i6)*Y(i6,i6)*V(i6)*cos(alf(i6,i6));L(i6,i6)=-Hd1(i6)+V(i6)*Y(i6,i6)*V(i6)*sin(alf(i6,i6));endendend%修正Jacobian矩阵for j9=3for i9=1:nN(i9,j9)=0;L(i9,j9)=0;J(j9,i9)=0;L(j9,i9)=0;endendL(j9,j9)=1;for j9=4for i9=1:nH(i9,j9)=0;N(i9,j9)=0;J(i9,j9)=0;L(i9,j9)=0;H(j9,i9)=0;N(j9,i9)=0;J(j9,i9)=0;L(j9,i9)=0;endendH(j9,j9)=1;L(j9,j9)=1;%Jaco=[H N;J L];%Jaco=zeros(2*n,2*n);for i11=1:nfor j11=1:nJaco(2*i11-1,2*j11-1)=H(i11,j11); Jaco(2*i11-1,2*j11)=N(i11,j11); Jaco(2*i11,2*j11-1)=J(i11,j11);Jaco(2*i11,2*j11)=L(i11,j11);endendA=Jaco;33333.数据:4 4 1.05 0.000011 12 0.1 0.40 0.015282 1 4 0.12 0.50 0.019203 24 0.08 0.40 0.014131 1 3 0 0.3 0.909090911 1 0 0 0.30 0.182 2 0 0 0.55 0.133 3 0.5 0 0 01 3 1.10 0 0。
matlab电力系统潮流计算程序
matlab电力系统潮流计算程序电力系统潮流计算是电力系统分析的关键步骤之一，用于确定电力系统各节点的电压和相角分布。
以下是一个简单的MATLAB电力系统潮流计算的基本步骤和代码示例：1.定义电力系统参数：-定义系统节点数量、支路数据、发电机数据、负荷数据等电力系统参数。
```matlab%电力系统参数busdata=[1,1.05,0,0,0,0,0,0;2,1.02,0,0,0,0,0,0;%...其他节点数据];linedata=[1,2,0.02,0.06,0.03;%...其他支路数据];gendata=[1,2,100,0,999,1.05,0.95;%...其他发电机数据];loaddata=[1,50,20;%...其他负荷数据];```2.构建潮流计算矩阵：-利用节点支路导纳、节点负荷和发电机功率等信息构建潮流计算的阻抗矩阵。
```matlabYbus=buildYbus(busdata,linedata);```3.迭代求解潮流方程：-利用迭代算法（如牛顿-拉夫森法）求解潮流方程，更新节点电压和相角。
```matlab[V,delta]=powerflow(Ybus,gendata,loaddata,busdata);```4.结果分析和可视化：-分析计算结果，可视化电压和相角分布。
```matlabplotVoltageProfile(busdata,V,delta);```这只是一个简单的潮流计算示例。
具体的程序实现可能涉及更复杂的算法和工程细节，取决于电力系统的复杂性和精确性要求。
您可能需要根据实际情况和数据格式进行调整和改进。
在实际工程中，也可以考虑使用专业的电力系统仿真软件。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
用Matlab计算潮流计算-电力系统分析《电力系统潮流上机》课程设计报告院系：电气工程学院班级：电088班_______学号：0812002221学生姓名：刘东昇________指导教师：张新松________设计周数：两周_________日期：2010年12月25日一、课程设计的目的与要求目的：培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识要求：基本要求：1.编写潮流计算程序；2.在计算机上调试通过；3.运行程序并计算出正确结果；4.写出课程设计报告二、设计步骤：1.根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。
Alt ；' T = r、二戶土旳「亠2.在给定的电力网络上画出等值电路图3.运用计算机进行潮流计算。
4.编写设计说明书。
三、设计原理1.牛顿-拉夫逊原理牛顿迭代法是取X0之后，在这个基础上，找到比x0更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。
牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f（x） = 0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。
电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。
为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。
牛顿一拉夫逊迭代法的一般步骤：（1）形成各节点导纳矩阵丫。
（2）设个节点电压的初始值U和相角初始值e还有迭代次数初值为0。
（3）计算各个节点的功率不平衡量。
（4）根据收敛条件判断是否满足，若不满足则向下进行。
（5）计算雅可比矩阵中的各元素。
（6）修正方程式个节点电压（7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。
（8）计算平衡节点输出功率和各线路功率2.网络节点的优化1）静态地按最少出线支路数编号这种方法由称为静态优化法。
在编号以前。
首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。
当由n个节点的出线支路相同时，则可以按任意次序对这n个节点进行编号。
这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也2）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化（增加，减少或保持不变）。
因此, 如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。
3. MATLAB^程应用Matlab 是“Matrix Laboratory "的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。
由于使用Matlab编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。
四、设计内容1.设计流程图设非平僑节点电JI 令迭代次对PQ 节点计算R（k ）.Q （k ）（对J计算节点的新电压e（k 1） e（k ） e （k ） ii iI 增加迭代次数coun t=co un t+12-程序clear;clc%重新编号，把原题中的节点 1,2,3,4,5重新依次编号为5,1,2,3,4 ，其中1-4号为PQ 节点,5号为平衡节点 y=0;%输入原始数据，求节点导纳矩阵 y (1,2)=1/(0.06+0.18i); y (1,3)=1/(0.06+0.18i);y (1,4)=1/(0.04+0.12i);y(1,5)=1/(0.02+0.06i);y(2,3)=1/(0.01+0.03i);y(2,5)=1/(0.08+0.24i); y(3,4)=1/(0.08+0.24i); y(4,5)=0; for i=1:5 for j=i:5 y(j,i)=y(i,j); end end*计算平 l衡节点丿输J卄Y=0;%求互导纳for i=1:5for j=1:5if i~=jY(i,j)=-y(i,j);endendend%求自导纳for i=1:5丫(i,i)=sum(y(i,:));end丫%Y为导纳矩阵G=real(Y);B=imag(Y);%原始节点功率S(1)=0.2+0.2i;S(2)=-0.45-0.15i;S(3)=-0.4-0.05i;S(4)=-0.6-0.1i;S(5)=0;P=real(S);Q=imag(S);%赋初值U=o nes(1,5);U(5)=1.06;e=zeros(1,5);ox=on es(8,1);fx=on es(8,1);cou nt=0 %计算迭代次数while max(fx)>1e-5for i=1:4for j=1:4 H(i,j)=0;N(i,j)=0;M(i,j)=0;L(i,j)=0;oP(i)=0;oQ(i)=0;endendfor i=1:4for j=1:5oP(i)=oP(i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*s in (e(i)-e(j))); oQ(i)=oQ(i)-U(i)*U(j)*(G(i,j)*si n(e(i)-e(j))-B(i,j)*cos(e(i)-e(j))); endoP(i)=oP(i)+P(i); oQ(i)=oQ(i)+Q(i);endfx=[oP,oQ]';%求雅克比矩阵%当i~=j 时候求H,N,M丄如下：for i=1:4for j=1:4if i~=jH(i,j)=-U(i)*U(j)*(G(i,j)*si n(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)))； N(i,j)=-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*s in (e(i)-e(j)));L(i,j)=H(i,j);M(i,j)=-N(i,j);endendendH,N,M 丄%当i=j 时H,N,M丄如下：for i=1:4for j=1:5if i~=jH(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*si n(e(i)-e(j))-B(i, j)*cos (e(i)-e(j)));N(i,i)=N(i,i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*s in (e(i)-e(j)));M(i,i)=M(i,i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*s in (e(i)-e(j)) );L(i,i)=L(i,i)-U(i)*U(j)*(G(i,j)*si n( e(i)-e(j))-B(i,j)*cos(e(i)-e(j)))；endendN(i,i)=N(i,i)-2*(U(i))A2*G(i,i);L(i,i)=L(i,i)+2*(U(i)f2*B(i,i);endJ=[H,N;M,L] %J为雅克比矩阵ox=-((i nv(J))*fx);for i=1:4oe(i)=ox(i); oU(i)=ox(i+4)*U(i);endfor i=1:4e(i)=e(i)+oe(i); U(i)=U(i)+oU(i);endcou nt=cou nt+1;endox,U,e,cou nt%求节点注入的净功率i=5;for j=1:5P(i)=U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*s in (e(i)-e(j)))+P(i); Q(i)=U(i)*U(j)*(G(i,j)*si n(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)))+Q(i); end S(5)=P(5)+Q(5)*sqrt(-1);S%求节点注入电流匸Y*U'3.运行结果Y值:¥ =10. 8333 -32.50001-1. 6667 + 5.00001-1.6067 + 5.OOOOi -2.5000 + 7. 5000i-5.0000 H5. OOQOi-L 6667 + 5+OOOCi12* 9167 -38. 7500i-10. 0000 +30. 000010-1.3500 + 3. 75001 -1. 6667 + 5.000(Ji-10. 0000 +30.0000112. 9167 -38. 75001-1.2500十3. 750010 -2.5000 + 7.500010-1. 2500 + 3. 7500i 3. 7500 -11. 260010-5. 0000 +15.00001-1. 2500 + 3.7500i00 6.2500 -18.7500i 迭代过程:courit —H =0 5. 0000 5. 0000 7. 5000H 000030_ 0000氐000030. 00003. 7SOOC 50003. 7500N -0U 5S67 1, 5567 2.50001.0007 010. 0000 01.6667 10. 00000 1. 25002. 5000 0 1.2500 0M =0—1- 6667 —1. 666? -2.5000-1. fififi? 0-10.00000-1.6667 -10.00000~1* 2500-2,5000 0-1.2500 00 5. 0000 5., 0000 7.500030000□o.oooo05.000030+0000Q3・75007. 50000 3.7500 a'33. 100Q 5. 0000 5.0000 7.5000 -10* 5333 1.i.. «667 2. 5000 氐 0000 -381. 975030, 0000 0 1. 6667 13. 3417 LO.0000 05. 0000 30. 0000 -38. 7500 3. 7500 1. 6€fl710. 0000 -1Z- 9107 1. 25007. BQQQ0工75QQ-11. 2500 2. 5000 a U 2500-工75oa1 k 1333-1, S66 : -1* 6667 -2. 5000-SLfiOOO E. OOOD E. 00007. 5000-1. 6067 12. 9917 -10. 0000 a 5.0000 -38.5250 30. 0000 0-LS667-10. 0000 12-91S7 -1.2500 5. O0U03d 0000-38. 7500 3. 7500-2.5000 0-L 250Q 3.7500 二5COO Q 3.7S00 -11. 250005, 35995* 302 28,0^96吕.2222 030.9520 OS, 2033 30.S2S3 O3* 8576了. 7Z2303,8183O01・S572 1. 5235 2.1E11L.970J010*11370U1i■■曲】町・-4SO30】■3.102&01-3384O0-1- 5572一L.5235-2* 11511-1. 97010-10, 1 1370一丄.目當曰吕-lDa 2S03O—1. 2^037・ 10 2&0-L.33840L -I0 5. 3E9D 5. 3622S. 03965. 2222 0廿Q・9520£ 2038 30, S298 03, 85767. 722303-8183J =-35. 0648 5. 3590 氐 3622 8> 0396 -12, 0310 1.55721+ 5235 2. 1EU5. 2222 -40. 0 794 30. 9520 01, 9701 -12.8027 10” 113706. 2038 30. S 299 -39・ 8&13 3. 8576 1.驱E 10. 4帥3 -12. Q6S2 1, 22037. 7223 0 3. 8183-11. 5406 3. 1029 0 1. 33SJ：—3. 199611. 5373 -h 5572 =1.52357 1511-皈6400 辰3594? 乩 3632-L 970113a7 720 ・ia 1157 05+2222-39+ 8247 30+ &520 a -E 99S5 -10- 4 903 13. SQQl-1. 2203 E” 203530.82DS -3&, 9105 3, 8576*3. 1029 0=1・3384 4・ 44127223 0 3. 3183 -11. 3818M ■O 6™ 15S5 13977.02175_ 230ao30* 4267QE. 2S1230™ 5 42 63. 77.92aao3七 SOL aQW ■Q1・93 93 1. Q0 QS3.01531« 5430o1 O- 3300oI. &1 139. i9ST6oI* 3 1T121.o1 _2 0 ESOOwr ■=-1 x S3Q3 -1. P廿耳甘-3. 0*53=1. 506o-10T 3^55O^■1,5 112—4 QS7S—1 . 3 l 7 1Y・ 13G9O-1 * 20&0O0 5.2904 5+ 2912 7+924^ 乩1顶030+ &4300 S. 139730_ 42B70 3. 8018 7.62170■ 7644034.71 63 E.2904 5, 2912 7* D246 -11* 3394 1.E436 1. 51125.1505 -39-5S49 30. 5426 01,9393 -12. &950 9, »S760, 139730-42A7-39, 3A8I 3. 80131.9038 10. 33B5-12, 70 as 札6Z 1T0 3. Toi 丄-1 lx 3K61 3. 0^53O 1. 317111. 4384 -1_ 6436-1. 61 12 -2.1368 -35.1173 6. 2S04 B. 2012 _ h 9393 13.5950 一乩 9076 0 5.1535-39.2653 3Q* 5426 -L 阴阳-1Q- 335513* 2T. 20!i0 5.1397节4 4707 -g ZQ83-九 0453 0-1.3171亠36237. 82170人 7^442.13 6Sr 20 50-3. 15 197- 92463.aois -]1L0 3 2SQ& 5. 29077.9237野.1压E0030. 53D50号.139230,42300 3.80)3二621003- 764000 1・E433 L. Bl 1 1 2.L36S 1・ 93910玄980 7 0 1- 955&JICL 3343O1-Z C49 3. 0148 0L-3169 O0 -1, '6436-1.5111一2・一1・ 939109867O -1a 9SSB-1 0_ 33430-J,20^^-1. 04480一—313900FL 2899 5.Z9077-9237 刁.1巨曲00□ 0 ™ 5*19 5Q3时.42J003-80137. 021003* 76100■34.7136 5. 2S995,29077.9237-tL S379 1.5436 1. 51115J580-39.581230.5395Q 1.9391-12.69379.9367 5.139230. 4236-39.3642 3.S013 1.965510.3343-12. 7047 7-52100 3.7640-11.3850 3.04480L3169 11.437&-L 5436-1.5111-Z136S-35.1136 5.2S99 5. 2907 -1.939113. 5937-9.98670 5. 1580-39.281230, 5395-1.9655-10.334313.5047 -1.2049 5. 139230.4236-39, 2642 -3.04430-1.3169 4.36177.621003*76402.13680 1.2049 -3.1617 7. ft2370 3.SD13 -11.1850氛 CIO 7 3 -O. Oifit -□., QS3 9 -Q ・ DS96 -Q a> 1044 c ouiit =平衡节点注入功率及电流:0. 2000 + 0. 2000i -0. 4500 - 0. 150Q1 -0. 4OQ0 - 0. 05001-0.5000 - 0. 100011*2932 + Q.24J0-0645 - 6 19351-9.0953 十 0.38581 -0.0565 + 0.16961*0.0944 + 0.2S321 0. 1817 - D. 54521五、课程设计心得与体会在电力系统分析课程中，我们学习了电力网络的 数学模型和确定解算方法。