用Matlab计算潮流计算
潮流计算(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编写电力系统潮流计算程序。
首先，我们需要建立电力系统的节点模型。
节点模型是电力系统中各节点的电压、功率和电流等参数的数学表示。
在Matlab中，我们可以使用矩阵来表示节点模型。
假设电力系统有n个节点，我们可以定义一个n×n的复数矩阵Y来表示节点之间的导纳关系，其中Y(i,j)表示节点i和节点j之间的导纳。
同时，我们还需要定义一个n×1的复数向量V来表示各节点的电压，其中V(i)表示节点i的电压。
接下来，我们需要编写潮流计算的主程序。
主程序的主要功能是根据节点模型和潮流计算算法，计算出各节点的电压、功率和电流等参数。
在Matlab中，我们可以使用循环语句和矩阵运算来实现潮流计算。
具体的计算过程可以参考电力系统潮流计算的算法。
在编写主程序之前，我们还需要定义一些输入参数，如电力系统的节点数、发电机节点和负荷节点等。
这些参数可以通过用户输入或者读取文件的方式获取。
同时，我们还需要定义一些输出参数，如各节点的电压、功率和电流等。
这些参数可以通过矩阵运算和循环语句计算得到，并输出到文件或者显示在屏幕上。
最后，我们需要进行程序的测试和调试。
可以通过输入一些测试数据，运行程序并检查输出结果是否正确。
如果发现程序有错误或者结果不准确，可以通过调试工具和打印调试信息的方式进行调试。
总之，利用Matlab编写电力系统潮流计算程序可以提高计算效率和准确性，为电力系统的运行和规划提供有力的支持。
当然，编写一个完整的潮流计算程序需要考虑很多细节和特殊情况，这需要有一定的电力系统和编程知识。
matlab3节点牛顿拉夫逊潮流计算
1. Matlab3节点牛顿拉夫逊潮流计算简介Matlab是一种高度灵活的编程语言和数学工具，被广泛应用于科学计算和工程领域。
3节点牛顿拉夫逊潮流计算是一种电力系统分析方法，用于计算电力系统中各节点的电压和相角。
在本文中，我们将介绍如何使用Matlab进行3节点牛顿拉夫逊潮流计算，以及该方法的原理和应用。
2. 3节点牛顿拉夫逊潮流计算原理3节点牛顿拉夫逊潮流计算是一种基于潮流方程的迭代算法。
它通过不断迭代求解节点电压和相角，以达到系统在给定负荷下的稳态。
其核心原理是利用牛顿拉夫逊法迭代求解潮流方程，即功率平衡方程和节点电压方程，直至收敛得到结果。
3. Matlab在3节点牛顿拉夫逊潮流计算中的应用Matlab提供了丰富的数学工具和函数库，使其成为进行电力系统分析的理想工具。
在3节点牛顿拉夫逊潮流计算中，我们可以利用Matlab编写相应的算法和程序，对实际电力系统进行分析和计算。
通过Matlab的矩阵运算和迭代算法，可以高效地求解潮流方程，得到系统各节点的电压和相角。
4. 3节点牛顿拉夫逊潮流计算的应用3节点牛顿拉夫逊潮流计算在电力系统规划、运行和故障分析中具有重要的应用价值。
通过计算系统各节点的电压和相角，可以评估系统的电压稳定性和潮流分布，指导电力系统的规划和调度。
在系统发生故障时，可以利用3节点牛顿拉夫逊潮流计算分析系统的稳定性和可靠性，为故障处理提供依据。
5. 结语3节点牛顿拉夫逊潮流计算是一种重要的电力系统分析方法，Matlab作为一种强大的数学工具，为其提供了理想的支持和实现。
通过Matlab进行3节点牛顿拉夫逊潮流计算，可以高效地进行电力系统分析和计算，为电力系统的规划和运行提供科学依据。
希望本文可以帮助读者更加深入地了解3节点牛顿拉夫逊潮流计算及其在Matlab 中的应用。
6. Matlab3节点牛顿拉夫逊潮流计算的优势利用Matlab进行3节点牛顿拉夫逊潮流计算具有许多优势。
Matlab提供了丰富的数学函数和工具，能够快速高效地进行矩阵运算、迭代求解等操作，极大地简化了算法的实现。
运用MATLAB软件进行潮流计算论文
摘要本文运用MATLAB软件进行潮流计算，对给定题目进行分析计算，再应用DDRTS软件，构建系统图进行仿真，最终得到合理的系统潮流。
潮流计算是电力系统最基本最常用的计算。
根据系统给定的运行条件，网络接线及元件参数，通过潮流计算可以确定各母线的电压幅值和相角，各元件流过的功率，整个系统的功率损耗。
潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。
因此，潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有着广泛的应用。
首先，画出系统的等效电路图，在计算出各元件参数的基础上，应用牛顿—拉夫逊Newton-Raphson法以及MATLAB软件进行计算对给定系统图进行了四种不同负荷下的潮流计算，经过调节均得到符合电压限制及功率限制的潮流分布。
其次，牛顿—拉夫逊Newton-Raphson法具有较好的收敛性，上述计算过程经过四到五次迭代后均能收敛。
根据运算结果，分析各支路损耗和系统总损耗。
最后，应用DDRTS软件，构建系统图，对给定负荷重新进行分析，潮流计算后的结果也能满足相应的参数要求。
关键词：牛顿-拉夫逊法 MATLAB DDRTS 潮流计算目录1．摘要 (2)2．题目原始资料 (2)3．题目分析 (5)4.题目求解 (6)1）根据题意要求画出等值电路 (6)2）读程序画出拉夫逊法的流程图 (7)3）变电所负荷为题目所给数据进行求解 (8)4)编写程序并运行 (10)5）具体调压调损耗过程 (10)1.改变变压器变比调压 (10)2.改变发电机机端电压调压 (12)3.负荷按照一定比例变化的潮流计算分析 (15)4.轮流断开支路双回线中的一条的潮流计算 (19)5.仿真并比较 (26)6.设计心得 (28)7.参考文献 (29)一、题目原始资料：1．系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。
变电所1 变电所2母线2、发电厂资料：母线1和2为发电厂高压母线，发电厂一总装机容量为（ 300MW ），母线3为机压母线，机压母线上装机容量为（ 100MW），最大负荷和最小负荷分别为50MW和20MW；发电厂二总装机容量为（ 200MW ）。
用matlab潮流计算课程设计
用matlab潮流计算课程设计一、教学目标本课程的目标是使学生掌握Matlab在电力系统潮流计算方面的应用。
通过本课程的学习，学生将能够：1.理解电力系统潮流计算的基本原理；2.熟悉Matlab软件的操作和编程；3.利用Matlab进行电力系统潮流计算；4.分析计算结果，判断电力系统的运行状态。
二、教学内容本课程的教学内容主要包括以下几个部分：1.电力系统潮流计算的基本概念和原理；2.Matlab软件的基本操作和编程方法；3.Matlab在电力系统潮流计算中的应用；4.电力系统潮流计算案例分析。
三、教学方法为了提高学生的学习兴趣和主动性，本课程将采用多种教学方法，包括：1.讲授法：用于讲解电力系统潮流计算的基本原理和Matlab软件的基本操作；2.案例分析法：通过分析实际案例，使学生掌握Matlab在电力系统潮流计算中的应用；3.实验法：让学生亲自动手进行电力系统潮流计算，加深对知识的理解和记忆。
四、教学资源为了支持教学内容和教学方法的实施，丰富学生的学习体验，我们将准备以下教学资源：1.教材：《Matlab潮流计算入门》；2.参考书：《电力系统潮流计算》、《Matlab编程与应用》；3.多媒体资料：Matlab软件教程、电力系统潮流计算案例视频；4.实验设备：计算机、Matlab软件、电力系统仿真软件。
五、教学评估本课程的评估方式将包括以下几个方面：1.平时表现：通过课堂参与、提问、讨论等环节，评估学生的学习态度和积极性；2.作业：布置相关的Matlab编程和电力系统潮流计算作业，评估学生的理解和应用能力；3.考试：进行期中考试和期末考试，评估学生对课程知识的掌握程度。
六、教学安排本课程的教学安排如下：1.教学进度：按照教材的章节顺序进行教学，确保每个章节都有足够的教学时间；2.教学时间：安排每周两次课堂授课，每次90分钟；3.教学地点：教室和实验室交替使用，以便进行实验和实践操作。
七、差异化教学根据学生的不同学习风格、兴趣和能力水平，我们将采取以下差异化教学措施：1.学习风格：提供多种学习资源，如视频教程、实验操作等，以满足不同学习风格的学生；2.兴趣：引入与电力系统相关的实际案例和应用场景，激发学生的学习兴趣；3.能力水平：设置不同难度的作业和项目，以便适应不同能力水平的学生。
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的电力系统潮流计算
摘要潮流计算是电力系统最基本、最常用的计算。
根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值及相角），各元件中流过的功率、整个系统的功率损耗等。
潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。
因此潮流计算在电力系统的规划设计、生产运行、调度管理及科学研究中都有着广泛的应用。
本次设计的主要目的就是面向一般的电力网络，形成节点导纳矩阵，确定合适的算法，编写通用的计算程序，得到计算结果。
设计中主要介绍了牛顿拉夫逊和PQ分解两种算法，PQ分解法虽然在结构上比牛顿法更加简化，但是针对一般网络现代计算机在存储空间及计算速度上已经十分强大，鉴于对牛顿法的熟悉与其算法的直观性，本次设计在编程时采用了牛顿拉夫逊法的直角坐标形式。
解方程的过程利用Matlab的强大计算功能，编写M语言，合理设置变量，实现通用计算功能。
关键词: 电力系统，潮流计算，牛顿—拉夫逊法，Matlab。
AbstractPower system load flow calculation is the most basic and commonly used calculations. Given according to the system operating conditions, the network connection and device parameters can be determined by power flow calculation of the bus voltage (magnitude and phase angle), the power flowing through the components, overall system power consumption and so on. Flow calculation is to achieve economic development of power system supply the necessary means and important part of the work. Therefore flow calculation in power system planning and design, production and operation, scheduling management, and scientific research have a wide range of applications.The main purpose of this design is for the general electricity network, the formation of the node admittance matrix, determine the appropriate method, the preparation of general-purpose computer program to get results. Introduces the design and the PQ decomposition Newton Raphson two algorithms, PQ decomposition although the structure is more streamlined than the Newton method, but for the general network of modern computer storage space and computing speed has been very strong, in view of the Newton Familiar with its intuitive algorithm, this design in programming using Newton Raphson polar form. The process of solving equations using matlab powerful computing capabilities, the preparation of M language, a reasonable set variables, to achieve general-purpose computing functions.Keywords: power system, power flow calculation, Newton - Raphson method, Matlab.目录摘要 (I)Abstract (II)目录 .................................................................................................................................................... I II 1 引言 .. (1)1.1 潮流计算目的 (1)1.2 潮流计算意义 (1)1.3 潮流计算发展史 (1)1.4基于MATLAB 的电力系统潮流计算发展前景 (2)2简单电力系统潮流计算 (4)2.1简单辐射网络的潮流计算 (4)2.1.1简单支路的潮流分布和电压降落 (4)2.1.2 辐射型网络的手工潮流计算方法 (6)2.2 简单环网的潮流计算 (7)2.2.1两端电压相等 (7)2.2.2两端电压不相等 (8)3 复杂电力系统潮流计算的计算机算法 (10)3.1电力网络方程及等值电路 (10)3.2节点导纳矩阵形成及修改 (11)3.3节点的分类 (14)3.3.1 PQ节点 (14)3.3.2 PV节点 (14)3.3.3 平衡节点 (14)3.4潮流计算的约束条件 (15)3.5牛顿-拉夫逊法（直角坐标） (15)3.5.1牛顿-拉夫逊法的推导过程 (15)3.5.2潮流计算时的修正方程（直角坐标） (17)3.5.3雅可比矩阵的特点: (19)3.5.4牛顿-拉夫逊法计算步骤 (19)3.6 P-Q分解法潮流计算 (20)3.6.1 P-Q分解法潮流计算概述 (20)3.6.2 P-Q分解法的潮流计算步骤 (20)3.6.3 P-Q分解法的特点 (21)4 Matlab概述 (22)4.1Matlab简介 (22)4.2 Matlab中的变量 (22)4.3 Matlab编程 (23)4.3.1矩阵的输入 (23)4.3.2矩阵的运算 (24)4.3.3 MatLab的控制流 (24)5 牛顿法潮流计算程序设计及实例 (26)5.1手算 (26)5.2计算机算法的数据输入 (29)5.3潮流计算程序 (30)5.3 计算结果分析 (36)结论 (37)参考文献 (38)附录A 程序流程图 (39)附录B Matlab仿真 (40)致谢 (43)1 引言1.1 潮流计算目的电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。
matlab潮流计算程序
基于MATLAB的电力系统潮流计算%简单潮流计算的小程序，相关的原始数据数据数据输入格式如下：%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初始化为0if 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初始化为0if 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.000001B1(支路参数矩阵)【1 2 0.04+0.25i 0.5i 1 0;1 3 0.1+0.35i 0 1 0;2 3 0.08+0.30i 0.5i 1 0;4 2 0.015i 0 1.05 1;5 3 0.03i 0 1.05 1】B2(节点参数矩阵)【0 -1.6-0.8i 1 0 0 2;0 -2-1i 1 0 0 2;0 -3.7-1.3i 1 0 0 2;0 5+0i 1.05 1.05 0 3;0 0 1.05 1.05 0 1】X(节点号和对地参数)【1 0;2 0;3 0;4 0;5 0】亲爱的朋友，上文已完，为感谢你的阅读，特加送另一篇范文，如果下文你不需要，可以下载后编辑删除，谢谢！矿井水灾事故专项应急演练方案1 应急演练目的、意义和目标1.1应急演练目的①评估我矿水灾事故的应急准备状态，发现并修改我矿水灾事故专项应急预案和执行程序中存在的缺陷和不足；②评估我矿在发生水灾事故时的应急能力，识别处理水灾事故的资源需求，澄清相关单位和人员的应急职责，改善水灾事故应急救援中的组织协调问题；③检验应急响应人员对水灾事故应急预案及执行程序的了解程度和实际操作技能；同时，通过调整演练难度，进一步培训和提高应急响应人员的业务素质和能力；④提高全员安全意识。
matlab潮流计算
附录1使用牛顿拉夫逊法进行潮流计算的Matlab程序代码% 牛拉法计算潮流程序%-----------------------------------------------------------------------% B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳% 5、支路的变比；6、支路首端处于K侧为1，1侧为0% B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值% 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量% 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；3为PV节点；%------------------------------------------------------------------------clear all;format long;n=input('请输入节点数:nodes=');nl=input('请输入支路数:lines=');isb=input('请输入平衡母线节点号:balance=');pr=input('请输入误差精度:precision=');B1=input('请输入由各支路参数形成的矩阵:B1=');B2=input('请输入各节点参数形成的矩阵: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:nl %支路数if B1(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)+B1(i,4); %对角元K侧Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4); %对角元1侧end%求导纳矩阵disp('导纳矩阵Y=');disp(Y)%-------------------------------------------------------------------G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部for i=1:n %给定各节点初始电压的实部和虚部e(i)=real(B2(i,3));f(i)=imag(B2(i,3));V(i)=B2(i,4); %PV节点电压给定模值endfor i=1:n %给定各节点注入功率S(i)=B2(i,1)-B2(i,2); %i节点注入功率SG-SLB(i,i)=B(i,i)+B2(i,5); %i节点无功补偿量end%---------------------------------------------------------------------P=real(S);Q=imag(S); %分解出各节点注入的有功和无功功率ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %迭代次数ICT1、a；不满足收敛要求的节点数IT2 while IT2~=0 % N0=2*n 雅可比矩阵的阶数；N=N0+1扩展列IT2=0;a=a+1;for i=1:nif i~=isb %非平衡节点C(i)=0;D(i)=0;for j1=1:nC(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)endP1=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) %求i节点有功和无功功率P',Q'的计算值V2=e(i)^2+f(i)^2; %电压模平方%以下针对非PV节点来求取功率差及Jacobi矩阵元素-----------------------------if B2(i,6)~=3 %非PV节点DP=P(i)-P1; %节点有功功率差DQ=Q(i)-Q1; %节点无功功率差%以上为除平衡节点外其它节点的功率计算--------------------------------------%求取Jacobi矩阵----------------------------------------------------------for j1=1:nif j1~=isb&j1~=i %非平衡节点&非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de=-dQ/dfX2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df=dQ/deX3=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;J(p,N)=DQ;m=p+1; % X3=dQ/de J(p,N)=DQ节点无功功率差J(m,q)=X1;J(m,N)=DP;q=q+1; % X1=dP/de J(m,N)=DP节点有功功率差J(p,q)=X4;J(m,q)=X2; % X4=dQ/df X2=dp/dfelseif j1==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/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;J(p,N)=DQ;%扩展列△Qm=p+1;J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列△PJ(m,q)=X2;endendelse%下面是针对PV节点来求取Jacobi矩阵的元素-----------------------------------------DP=P(i)-P1; % PV节点有功误差DV=V(i)^2-V2; % PV节点电压误差for j1=1:nif j1~=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;J(p,N)=DV; % PV节点电压误差m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV节点有功误差J(m,q)=X2;elseif j1==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;J(p,N)=DV; % PV节点电压误差m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV节点有功误差J(m,q)=X2;endendendendend%以上为求雅可比矩阵的各个元素及扩展列的功率差或电压差---------------------------------------for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即N=2*n+1扩展列△P、△Q 或△Ufor k2=k1:N1 % 从k+1列的Jacobi元素到扩展列的△P、△Q 或△UJ(k,k2)=J(k,k2)./J(k,k);% 用K行K列对角元素去除K行K列后的非对角元素进行规格化endJ(k,k)=1; % 对角元规格化K行K列对角元素赋1%回代运算-------------------------------------------------------------------if k~=3 % 不是第三行k > 3k4=k-1;for k3=3:k4 % 用k3行从第三行开始到当前行的前一行k4行消去for k2=k1:N1 % k3行后各行上三角元素J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算（当前行k列元素消为0）end %用当前行K2列元素减去当前行k列元素乘以第k行K2列元素J(k3,k)=0; %当前行第k列元素已消为0endif k==N0 %若已到最后一行break;end%前代运算------------------------------------------------------------for k3=k1:N0 % 从k+1行到2*n最后一行for k2=k1:N1 % 从k+1列到扩展列消去k+1行后各行下三角元素J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算end %用当前行K2列元素减去当前行k列元素乘以第k行K2列元素J(k3,k)=0; %当前行第k列元素已消为0endelse %是第三行k=3%第三行k=3的前代运算----------------------------------------------------for k3=k1:N0 %从第四行到2n行（最后一行）for k2=k1:N1 %从第四列到2n+1列（即扩展列）J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算（当前行3列元素消为0）end %用当前行K2列元素减去当前行3列元素乘以第三行K2列元素J(k3,k)=0; %当前行第3列元素已消为0endendend%上面是用线性变换方式高斯消去法将Jacobi矩阵化成单位矩阵%-----------------------------------------------------------------------------------for k=3:2:N0-1L=(k+1)./2;e(L)=e(L)-J(k,N); %修改节点电压实部k1=k+1;f(L)=f(L)-J(k1,N); %修改节点电压虚部end%修改节点电压---------------------------for k=3:N0DET=abs(J(k,N));if DET>=pr %电压偏差量是否满足要求IT2=IT2+1; %不满足要求的节点数加1endendICT2(a)=IT2; %不满足要求的节点数ICT1=ICT1+1; %迭代次数end%用高斯消去法解"w=-J*V"disp('迭代次数：');disp(ICT1);disp('没有达到精度要求的个数：');disp(ICT2);for k=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)*1i; %将各节点电压用复数表示end%计算各输出量------------------------------------------------------disp('各节点的实际电压标幺值E为：');disp(E); %显示各节点的实际电压标幺值E用复数表示disp('-----------------------------------------------------');disp('各节点的电压大小V为：');disp(V); %显示各节点的电压大小V的模值disp('-----------------------------------------------------');disp('各节点的电压相角deg为：');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('各节点的功率S为：');disp(S); %显示各节点的注入功率disp('-----------------------------------------------------');disp('各条支路的首端功率Si为：');for i=1:nlp=B1(i,1);q=B1(i,2);if B1(i,6)==0Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))...-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));Siz(i)=Si(p,q);elseSi(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))...-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));Siz(i)=Si(p,q);enddisp(Si(p,q));SSi(p,q)=Si(p,q);ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];disp(ZF);disp('-----------------------------------------------------');enddisp('各条支路的末端功率Sj为：');for i=1:nlp=B1(i,1);q=B1(i,2);if B1(i,6)==0Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))...-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));Sjy(i)=Sj(q,p);elseSj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))...-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));Sjy(i)=Sj(q,p);enddisp(Sj(q,p));SSj(q,p)=Sj(q,p);ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];disp(ZF);disp('-----------------------------------------------------');enddisp('各条支路的功率损耗DS为：');for i=1:nlp=B1(i,1);q=B1(i,2);DS(i)=Si(p,q)+Sj(q,p);disp(DS(i));DDS(i)=DS(i);ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];disp(ZF);disp('-----------------------------------------------------');end附录2使用PQ分解法进行潮流计算的Matlab程序代码%PQ分解法潮流计算程序%本文中的实例数据如下：节点数为9；支路数为9；平衡母线节点号为1；误差精度为0.00001；PQ节点数为5；%主程序clear all;format long;n=input('请输入节点数：n=');nl=input('请输入支路数：nl=');isb=input('请输入平衡母线节点号：isb=');pr=input('请输入误差精度：pr=');B1=input('请输入由支路参数形成的矩阵：B1='); %输入B1B2=input('请输入由支路参数形成的矩阵：B2='); %输入B2na=input('请输入PQ节点数na=');Y=zeros(n);YI=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);for i=1:nlif B1(i,6)==0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));YI(p,q)=YI(p,q)-1./B1(i,3);Y(q,p)=Y(p,q);YI(q,p)=YI(p,q);Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4);YI(q,q)=YI(q,q)+1./B1(i,3);Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4);YI(p,p)=YI(p,p)+1./B1(i,3);end %求导纳矩阵disp('节点导纳矩阵为：');disp(Y);G=real(Y);B=imag(YI);BI=imag(Y);for i=1:nS(i)=B2(i,1)-B2(i,2);BI(i,i)=BI(i,i)+B2(i,5);endP=real(S);Q=imag(S);for i=1:ne(i)=real(B2(i,3));f(i)=imag(B2(i,3));V(i)=B2(i,4);endfor i=1:nif B2(i,6)==2V(i)=sqrt(e(i)^2+f(i)^2);O(i)=atan(f(i)./e(i));endendfor i=2:nif i==nB(i,i)=1./B(i,i);else IC1=i+1;for j1=IC1:nB(i,j1)=B(i,j1)./B(i,i);endB(i,i)=1./B(i,i);for k=i+1:nfor j1=i+1:nB(k,j1)=B(k,j1)-B(k,i)*B(i,j1);endendendendp=0;q=0;for i=1:nif B2(i,6)==2p=p+1;k=0;for j1=1:nif B2(j1,6)==2k=k+1;A(p,k)=BI(i,j1);endendendendfor i=1:naif i==naA(i,i)=1./A(i,i);else k=i+1;for j1=k:naA(i,j1)=A(i,j1)./A(i,i);endA(i,i)=1./A(i,i);for k=i+1:nafor j1=i+1:naA(k,j1)=A(k,j1)-A(k,i)*A(i,j1);endendendendICT2=1;ICT1=0;kp=1;kq=1;K=1;DET=0;ICT3=1;while ICT2~=0||ICT3~=0ICT2=0;ICT3=0;for i=1:nif i~=isbC(i)=0;for k=1:nC(i)=C(i)+V(k)*(G(i,k)*cos(O(i)-O(k))+BI(i,k)*sin(O(i)-O(k)));endDP1(i)=P(i)-V(i)*C(i);DP(i)=DP1(i)./V(i);DET=abs(DP1(i));if DET>=prICT2=ICT2+1;endendendNp(K)=ICT2;if ICT2~=0for i=2:nDP(i)=B(i,i)*DP(i);if i~=nIC1=i+1;for k=IC1:nDP(k)=DP(k)-B(k,i)*DP(i);endelsefor LZ=3:iL=i+3-LZ;IC4=L-1;for MZ=2:IC4I=IC4+2-MZ;DP(I)=DP(I)-B(I,L)*DP(L);endendendendfor i=2:nO(i)=O(i)-DP(i);endkq=1;L=0;for i=1:nif B2(i,6)==2C(i)=0;L=L+1;for k=1:nC(i)=C(i)+V(k)*(G(i,k)*sin(O(i)-O(k))-BI(i,k)*cos(O(i)-O(k)));endDQ1(i)=Q(i)-V(i)*C(i);DQ(L)=DQ1(i)./V(i);DET=abs(DQ1(i));if DET >=prICT3=ICT3+1;endendendelse kp=0;if kq~=0;L=0;for i=1:nif B2(i,6)==2C(i)=0;L=L+1;for k=1:nC(i)=C(i)+V(k)*(G(i,k)*sin(O(i)-O(k))-BI(i,k)*cos(O(i)-O(k)));endDQ1(i)=Q(i)-V(i)*C(i);DQ(L)=DQ1(i)./V(i);DET=abs(DQ1(i));endendendendNq(K)=ICT3;if ICT3~=0L=0;for i=1:naDQ(i)=A(i,i)*DQ(i);if i==nafor LZ=2:iL=i+2-LZ;IC4=L-1;for MZ=1:IC4I=IC4+1-MZ;DQ(I)=DQ(I)-A(I,L)*DQ(L);endendelseIC1=i+1;for k=IC1:naDQ(k)=DQ(k)-A(k,i)*DQ(i);endendendL=0;for i=1:nif B2(i,6)==2L=L+1;V(i)=V(i)-DQ(L);endendkp=1;K=K+1;elsekq=0;if kp~=0K=K+1;endendfor i=1:nDy(K-1,i)=V(i);endenddisp('迭代次数');disp(K);disp('每次没有达到精度要求的有功功率个数为');disp(Np);disp('每次没有达到精度要求的无功功率个数为');disp(Nq);for k=1:nE(k)=V(k)*cos(O(k))+V(k)*sin(O(k))*j;O(k)=O(k)*180./pi;enddisp('各节点的电压标幺值E为：');disp(E);disp('各节点的电压V大小为：');disp(V);disp('各节点的电压相角O为：');disp(O);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);enddisp('各节点的功率S为：');disp(S);disp('各条支路的首端功率Sj为：');for i=1:nlif B1(i,6)==0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endSi(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-conj(E(q)))*conj(1./( B1(i,3)*B1(i,5))));disp(Si(p,q));enddisp('各条支路的末端功率Sj为：');for i=1:nlif B1(i,6)==0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endSj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-conj(E(p)))*conj(1./( B1(i,3)*B1(i,5))));disp(Sj(q,p));enddisp('各条支路的功率损耗DS为：');for i=1:nlif B1(i,6)==0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endDS(i)=Si(p,q)+Sj(q,p);disp(DS(i));endfor i=1:KCs(i)=i;for j=1:nDy(K,j)=Dy(K-1,j);endend附录3进行三相短路容量计算的Matlab程序代码%短路容量计算程序%---------------------------------------------------------------------% B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳% 5、支路的变比；6、支路首端处于K侧为1，1侧为0% B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值% 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量% 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；%Yd为修改后的节点导纳矩阵%-----------------------------------------------------------------------clear all;format long;g1=input('300MW发电机数：g1=');g2=input('250MW发电机数：g2=');n=input('请输入节点数:n=');nl=input('请输入支路数:nl=');B1=input('请输入由各支路参数形成的矩阵:B1=');B2=input('请输入各节点参数形成的矩阵:B2=');Y=zeros(n);% Y为修改前节点导纳矩阵for i=1:nl %支路数if B1(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)+B1(i,4);%对角元K侧Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4); %对角元1侧end%----------------------------------------------------------%Y2-Y5为各PQ节点负荷的导纳Y1=0;Y2=conj(B2(2,2));Y3= conj(B2(3,2));Y4= conj(B2(4,2));Y5= conj(B2(5,2)); Xd300=0.51j/0.950413^2;XT300=0.033212j/0.950413^2;Xd250=0.714j/0.950413^2;XT250=0.038747j/0.950413^2;Y6=g1/(XT300+Xd300)+g2/(XT250+Xd250);%处理相应的负荷及机组部分的导纳%------------------------------------------------------------C=[Y1,Y2,Y3,Y4,Y5,Y6];Yd=Y;for i=1:nYd(i,i)=Yd(i,i)+C(i); %修改各节点自导纳enddisp(Yd);Z = inv(Yd); %求节点阻抗矩阵for j=1:nI(j) = 1/Z(j,j); %电压故障前电压标幺值为1S(j)=abs(I(j));Sn(j)=S(j)*100;%短路电流有名值end%计算完毕----------------------------------------------------disp('各节点短路时的短路电流幅值标幺值')disp(abs(I))disp('短路容量有名值Sn=');disp(Sn);。
潮流计算方法
Pi N ii Vi Vi V j (Gij cos ij Bij sin ij ) 2Vi 2 Gii Vi 2 Gii Pi ji Vi j i
J ii
Qi Vi V j (Gij cos ij Bij sin ij ) Vi 2 Gii初始值及线路数据输入excel文 件中，并通过xlsread这一函数输入到MATLAB 中。
形成导纳矩阵
1.形成一个零矩阵Y 2.检测线路参数前两列，将Y对应元素加上相应 参数。
形成导纳矩阵
1.非对角线元素：
for n=1:nbr Y(n1(n),n2(n))= Y(n1(n),n2(n))-Y1(n)/(a(k)^2); Y(n2(n),n1(n))= Y(n1(n),n2(n)); end
潮流计算方法 —计算机算法
————
需求分析：
使用MATLAB软件计算美国电气服务公司 IEEE30节点系统（图见课本P136页图6.16） 各节点的电压幅值和相角的潮流分布，并比较 各种算法的优劣。
总体思路：
1.直接使用“牛顿——拉夫逊”法，通过图像 显示每次迭代后的误差并记录程序运行时间。 2.直接使用“PQ分解法”，通过图像显示每次 迭代后的误差并记录程序运行时间。 3.首先使用直流潮流计算，计算各节点电压的 相角，并将计算结果作为 “牛顿——拉夫逊” 法的初值，通过图像显示每次迭代后的误差并 记录程序运行时间。
迭代步骤：
1.根据各节点电压幅值及相角初值，计算电压 相量。 2.计算各PQ节点的复功率。 3.计算各节点流出的电流初值。 4.根据S=VI*计算出各节点复功率，并与已知复 程序 功率相减得到误差。 5.当误差小于设定值时，迭代结束。
matlab实验电力系统潮流计算
matlab实验电力系统潮流计算电力系统潮流计算是电力系统运行分析的基础，它通过计算电力系统中各节点的电压幅值和相角，以及各支路的功率和电流，来研究电力系统的稳态工作状态。
本文将介绍潮流计算的原理及其在电力系统中的应用。
潮流计算的基本原理是基于电力系统节点间的功率平衡方程，即节点注入功率等于节点负荷消耗功率加上节点发电机注入功率和节点之间传输输电功率的代数和。
潮流计算通常分为直流潮流计算和交流潮流计算两种方法。
直流潮流计算是指忽略电流相位差的计算方法。
在直流潮流计算中，电力系统的节点电压幅值和相角可以用复数来表示，节点注入功率和节点负荷消耗功率也采用复功率的形式。
直流潮流计算的基本方程为：P+iQ = V*(Gcosθ+Bsinθ)其中，P和Q分别表示节点注入有功功率和无功功率，V表示节点电压幅值，θ表示节点电压相角，G和B分别表示节点导纳矩阵的实部和虚部。
交流潮流计算考虑了电流相位差的影响，是更为准确的潮流计算方法。
交流潮流计算通常采用牛顿-拉夫逊法（Newton-Raphson method）进行迭代求解。
该方法以功率不平衡最小为目标，通过迭代计算更新节点电压幅值和相角，直到收敛为止。
潮流计算在电力系统运行和规划中具有重要的应用价值。
首先，潮流计算可以用来评估电力系统的稳态工作状态，包括节点电压幅值和相角、支路功率和电流等信息。
通过分析潮流计算结果，可以发现电力系统中潜在的潮流瓶颈和潮流分布情况，为电网调度和运行提供参考依据。
其次，潮流计算可以用来优化电力系统的运行和规划。
通过分析潮流计算结果，可以确定潮流分布不均衡的节点和支路，进而优化电力系统的输电和变电容量配置，提高电力系统的可靠性、经济性和稳定性。
此外，潮流计算还可用于电力系统的故障分析和重构，对于故障点的电压幅值、相角以及故障后的支路功率和电流进行分析，有助于电力系统故障的定位和恢复。
总的来说，电力系统潮流计算是电力系统运行分析的重要工具，通过计算电力系统中各节点的电压幅值和相角，以及各支路的功率和电流，可以评估电力系统的稳态工作状态，优化电力系统的运行和规划，实现电力系统的安全、稳定和高效运行。
基于matlab的潮流计算
课程设计任务书第学期学院：信息商务学院专业：电气工程及其自动化学生姓名：学号课程设计题目：基于MATLAB的潮流计算起迄日期：月日～月日课程设计地点：电气工程系工程实验室指导教师：系主任：课程设计任务书课程设计任务书1 引言潮流计算是电力系统最基本最常用的计算。
根据系统给定的运行条件，网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值和相角），各支路流过的功率，整个系统的功率损耗。
潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。
因此，潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有着广泛的应用。
潮流计算在数学上是多元非线性方程组的求解问题，牛顿—拉夫逊Newton-Raphson法是数学上解非线性方程组的有效方法，有较好的收敛性。
运用电子计算机计算一般要完成以下几个步骤：建立数学模型，确定解算方法，制订计算流程，编制计算程序。
1.1 设计任务与要求利用MATLAB/SIMULINK仿真软件，搭建多电压等级电力网络潮流计算模型。
掌握各种潮流计算方法，掌握MATLAB/SPS工具箱使用方法及其在潮流计算中的应用。
MATLAB作为先进的仿真软件，在电力系统潮流计算方面，尤其是矩阵运算上具有突出的优秀性。
本设计要求利用MATLAB编写程序，输入数据并完成潮流分布计算。
1.2实用价值与理论意义在电力系统的正常运行中，随着用电负荷的变化和系统运行方式的改变，网络中的损耗也将发生变化。
要严格保证所有的用户在任何时刻都有额定的电压是不可能的，因此系统运行中个节点出现电压的偏移是不可避免的。
为了保证电力系统的稳定运行，要进行潮流调节。
随着电力系统及在线应用的发展，计算机网络已经形成，为电力系统的潮流计算提供了物质基础。
电力系统潮流计算是电力系统分析计算中最基本的内容，也是电力系统运行及设计中必不可少的工具。
根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线电压的幅值及相角、各元件中流过的功率、整个系统的功率损耗等。
matlab电力系统潮流计算
电力系统潮流计算是电力系统分析和运行中的重要环节。
潮流计算主要用于确定电力系统中各个节点的电压、功率和潮流方向，以便进行功率平衡、电压稳定和线路负荷等方面的分析和评估。
MATLAB作为一种强大的数学建模和仿真工具，被广泛应用于电力系统潮流计算的研究和实际工程中。
本文将介绍MATLAB在电力系统潮流计算中的应用，包括算法原理、建模步骤和实例分析等内容。
一、潮流计算的基本原理潮流计算是指在给定电力网拓扑结构、负荷信息和发电机功率的情况下，通过迭代计算求解节点电压的复数值，以确定各节点的电压幅值和相角，进而计算各支路和各节点上的有功和无功功率。
潮流计算的基本原理是基于电力系统的潮流方程和节点功率平衡等基本理论，通过建立节点电压的复数方程组，利用迭代计算方法求解该方程组，从而得到节点的电压和功率信息。
二、MATLAB在潮流计算中的应用MATLAB作为一种功能强大的数学建模和仿真工具，具有丰富的数学计算函数和图形显示功能，适合于电力系统潮流计算的建模和仿真。
在MATLAB环境下，可以利用其矩阵运算、方程求解和数据可视化等功能，实现电力系统潮流计算的数学模型和算法的实现。
下面将介绍MATLAB在电力系统潮流计算中的具体应用步骤。
1. 建立电力系统潮流计算的数学模型在MATLAB环境下，首先需要建立电力系统潮流计算的数学模型，包括节点电压方程、支路潮流方程、节点功率平衡方程等。
利用MATLAB的矩阵运算和符号计算工具，可以将电力系统的节点和支路参数、负荷信息、发电机功率等数据表示为矩阵形式，建立电力系统潮流计算的数学模型。
2. 编写潮流计算的求解算法在建立电力系统潮流计算的数学模型后，需要编写潮流计算的求解算法。
在MATLAB环境下，可以利用其丰富的数学计算函数和优化工具，实现潮流计算的迭代求解算法，包括高斯-赛德尔迭代法、牛顿-拉夫逊迭代法等。
通过编写求解算法，可以实现电力系统潮流计算的数值求解过程。
3. 进行潮流计算的仿真实验在完成潮流计算的求解算法后，可以利用MATLAB进行潮流计算的仿真实验。
MATLAB在潮流计算演示教学中的应用
MATLAB在潮流计算演示教学中的应用MATLAB在潮流计算演示教学中的应用潮流计算是电力系统分析的重要内容之一，它可以求解电网中各节点电压、电流、功率、潮流等参数。
在电力系统教学中，潮流计算是一个重要的课程内容，也是学生理论知识与实践运用的桥梁。
因此，如何提高潮流计算教学的效果是学术界和工业界共同关注的问题。
MATLAB是一款功能强大的计算软件，广泛应用于工程、科学和教育领域，在潮流计算教学中也有重要的应用。
一、MATLAB在潮流计算仿真中的应用MATLAB可以进行各种电力系统仿真，包括潮流计算仿真。
它可以计算无功补偿、阻抗、励磁等各种因素对电网的影响，以及计算故障流、负荷流、电压降等等电网的各种参数。
通过使用MATLAB，学生可以模拟各种操作和状态，了解电网在不同负载状态下的响应，并研究电路参数对电网的影响。
二、MATLAB在潮流计算可视化中的应用MATLAB提供了强大的图形功能，可以用来可视化潮流计算。
学生可以绘制节点电压、相角、功率等参数的图示，以及各个节点之间电流流向的关系图。
在课堂上，教师可以通过这些可视化图像演示潮流计算的过程，提高学生的理解能力。
三、MATLAB在潮流计算分析中的应用MATLAB还可以用于对潮流计算结果进行分析，发现其中的规律和特点。
学生可以借助MATLAB进行统计分析、回归分析、协方差分析等方法，对潮流计算的结果进行深入探究。
这些分析结果有助于学生深入理解电网潮流的本质，提高其对电力系统的认知。
四、MATLAB在潮流计算算法优化中的应用MATLAB还可以用于潮流计算算法的优化，参与设计各种算法的控制语句和流程，包括线性规划、非线性规划、最小二乘法等多种算法。
学生可以在课程中了解各种算法的基本原理和具体操作流程，进而掌握MATLAB以及其他编程工具的使用方法，提高自主研究和创新能力。
综上所述，MATLAB在潮流计算教学中发挥了重要的作用，为学生提供了一个高效和灵活的学习平台。
matlab潮流计算
前言电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。
在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。
可靠性和经济性。
此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。
所以潮流计算是研究电力系统的一种很重要和很基础的计算。
随着科学技术的发展，电力系统变得越来越复杂，电气工程师掌握一种好的能对电力系统进行仿真的软件是学习和研究的需要。
文章简要介绍了MATLAB发展历史、组成和强大的功能，并用简单例子分别就编程和仿真两方面分析了MATIAB软件在电力系统研究中的具体应用。
采取等效电路法，能对特殊、复杂地电力系统进行高效仿真研究，因此，掌握编程和仿真是学好MATLAB的基础。
与众多专门的电力系统仿真软件相比，MATLAB软件具有易学、功能强大和开放性好，是电力系统仿真研究的有力工具。
1电力系统的基本概念电力系统：发电机把机械能转化为电能，电能经变压器和电力线路输送并分配到用户，在那里经电动机、电炉和电灯等设备又将电能转化为机械能、热能和光能等。
这些生产、变换、输送、分配、消费电能的发电机、变压器、变换器、电力线路及各种用电设备等联系在一起组成的统一整体称为电力系统。
电力网：电力系统中除发电机和用电设备外的部分。
动力系统：电力系统和“动力部分”的总和。
2潮流计算2.1潮流计算概述与发展电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。
利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。
在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。
对潮流计算的要求可以归纳为下面几点：1）计算方法的可靠性或收敛性；2）对计算机内存量的要求；3）计算速度；4）计算的方便性和灵活性。
MATLAB潮流计算仿真
附录A MATLAB程序%本程序的功能是用牛顿——拉夫逊法进行潮流计算% B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳% 5、支路的变比；6、支路首端处于K侧为1，1侧为0% B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值% 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量% 6、节点分类标号clear;n=10;%input('请输入节点数：n=');nl=11;%input('请输入支路数：nl=');isb=1;%input('请输入平衡母线节点号：isb=');pr=0.00001;%input('请输入误差精度：pr=');B1=[1 2 1.755e-2+4.155e-2i 0.26i 1 0;1 4 3.159e-2+7.479e-2i 0.1215i 1 0;1 6 3.159e-2+7.479e-2i 0.1215i 1 0;2 3 3.68e-3+0.11135i 0 0.909 1;4 5 3.68e-3+0.11135i 0 0.909 1;4 6 2.808e-2+6.648e-2i 0.108i 1 0;6 7 3.0865e-3+0.0833i 0 0.909 1;6 8 3.159e-2+7.479e-2i 0.1215i 1 0;6 10 2.457e-2+5.817e-2i 0.0945i 1 0;8 9 3.0865e-3+0.0833i 0 0.909 1;8 10 2.808e-2+6.648e-2i 0.108i 1 0];%input('请输入由支路参数形成的矩阵： B1=');B2=[0 0 1.05 1.05 0 1;0 0 1 0 0 2;0 0.6+0.3718i 1 0 0 2;0 0 1 0 0 2;0 0.4+0.247i 1 0 0 2;0 0 1 0 0 2;0 0.35+0.2169i 1 0 0 2;0 0 1 0 0 2;0 0.5+0.3099i 1 0 0 2;0.8 0 1.05 1.05 0 3];%input('请输入各节点参数形成的矩阵： B2=');Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n) ;S1=zeros(nl);%－－－－－－－修改部分－－－－－－－－－－－－ym=0;SB=100;UB=220;%ym=input('您输入的参数是标么值？（若不是则输入一个不为零的数值）');if ym~=0%SB=input('请输入功率基准值:SB=');%UB=input('请输入电压基准值:UB=');YB=SB./UB./UB;BB1=B1;BB2=B2;for i=1:nlB1(i,3)=B1(i,3)*YB;B1(i,4)=B1(i,4)./YB;enddisp('B1矩阵B1=');disp(B1)for i=1:nB2(i,1)=B2(i,1)./SB;B2(i,2)=B2(i,2)./SB;B2(i,3)=B2(i,3)./UB;B2(i,4)=B2(i,4)./UB;B2(i,5)=B2(i,5)./SB;enddisp('B2矩阵B2=');disp(B2)end% % %---------------------------------------------------for i=1:nl %支路数if B1(i,6)==0 %左节点处于低压侧p=B1(i,1);q=B1(i,2);elsep=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)+B1(i,4)./2; %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元1侧end%求导纳矩阵disp('导纳矩阵 Y=');disp(Y)%----------------------------------------------------------G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部for i=1:n %给定各节点初始电压的实部和虚部e(i)=real(B2(i,3));f(i)=imag(B2(i,3));V(i)=B2(i,4); %PV节点电压给定模值endfor i=1:n %给定各节点注入功率S(i)=B2(i,1)-B2(i,2); %i节点注入功率SG-SLB(i,i)=B(i,i)+B2(i,5); %i节点无功补偿量end%================================================================ ===P=real(S);Q=imag(S);ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;while IT2~=0IT2=0;a=a+1;for i=1:nif i~=isb %非平衡节点C(i)=0;D(i)=0;for j1=1:nC(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)endP1=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)%求P',Q'V2=e(i)^2+f(i)^2; %电压模平方%========= 以下针对非PV节点来求取功率差及Jacobi矩阵元素 =========if B2(i,6)~=3 %非PV节点DP=P(i)-P1; %节点有功功率差DQ=Q(i)-Q1; %节点无功功率差%=============== 以上为除平衡节点外其它节点的功率计算 ================= %================= 求取Jacobi矩阵 ===================for j1=1:nif j1~=isb&j1~=i %非平衡节点&非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de=-dQ/dfX2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df=dQ/deX3=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;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;elseif j1==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/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;J(p,N)=DQ;%扩展列△Qm=p+1;J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列△PJ(m,q)=X2;endendelse%=============== 下面是针对PV节点来求取Jacobi矩阵的元素===========DP=P(i)-P1; % PV节点有功误差DV=V(i)^2-V2; % PV节点电压误差for j1=1:nif j1~=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;J(p,N)=DV;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2;elseif j1==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;J(p,N)=DV;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2;endendendendend%========= 以上为求雅可比矩阵的各个元素 =====================for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即 N=2*n+1扩展列△P、△Q for k2=k1:N1 % 扩展列△P、△QJ(k,k2)=J(k,k2)./J(k,k); % 非对角元规格化endJ(k,k)=1; % 对角元规格化if k~=3 % 不是第三行%============================================================k4=k-1;for k3=3:k4 % 用k3行从第三行开始到当前行前的k4行消去for k2=k1:N1 % k3行后各行下三角元素J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算endJ(k3,k)=0;endif k==N0break;end%==========================================for k3=k1:N0for k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算endJ(k3,k)=0;endelsefor k3=k1:N0for k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算endJ(k3,k)=0;endendend%====上面是用线性变换方式将Jacobi矩阵化成单位矩阵=====for k=3:2:N0-1L=(k+1)./2;e(L)=e(L)-J(k,N); %修改节点电压实部k1=k+1;f(L)=f(L)-J(k1,N); %修改节点电压虚部end%------修改节点电压-----------for k=3:N0DET=abs(J(k,N));if DET>=pr %电压偏差量是否满足要求IT2=IT2+1; %不满足要求的节点数加1endendICT2(a)=IT2;ICT1=ICT1+1;end%用高斯消去法解"w=-J*V"disp('迭代次数：');disp(ICT1);disp('没有达到精度要求的个数：');disp(ICT2);for k=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);EE=E*UB;disp(EE);disp('-----------------------------------------------------'); disp('各节点的电压大小V为(节点号从小到大排列)：');disp(V);VV=V*UB;disp(VV);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);enddisp('各节点的功率S为(节点号从小到大排列)：');disp(S);disp('~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~');SS=S*SB;disp(SS);disp('-----------------------------------------------------'); disp('各条支路的首端功率Si为(顺序同您输入B1时一致)：');for i=1:nlp=B1(i,1);q=B1(i,2);if B1(i,6)==0Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-con j(E(q)))*conj(1./(B1(i,3)*B1(i,5))));Siz(i)=Si(p,q);elseSi(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))-co nj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));Siz(i)=Si(p,q);enddisp(Si(p,q));SSi(p,q)=Si(p,q)*SB;ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];disp(ZF);%disp(SSi(p,q));disp('-----------------------------------------------------'); enddisp('各条支路的末端功率Sj为(顺序同您输入B1时一致)：');for i=1:nlp=B1(i,1);q=B1(i,2);if B1(i,6)==0Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-co nj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));Sjy(i)=Sj(q,p);elseSj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))-con j(E(p)))*conj(1./(B1(i,3)*B1(i,5))));Sjy(i)=Sj(q,p);enddisp(Sj(q,p));SSj(q,p)=Sj(q,p)*SB;ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];disp(ZF);%disp(SSj(q,p));disp('-----------------------------------------------------'); enddisp('各条支路的功率损耗DS为(顺序同您输入B1时一致)：');for i=1:nlp=B1(i,1);q=B1(i,2);DS(i)=Si(p,q)+Sj(q,p);disp(DS(i));DDS(i)=DS(i)*SB;ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))]; disp(ZF);%disp(DDS(i));disp('-----------------------------------------------------'); endfigure(1);subplot(2,2,1);plot(V);xlabel('节点号');ylabel('电压标幺值');grid on;subplot(2,2,2);plot(sida);xlabel('节点号');ylabel('电压角度');grid on;subplot(2,2,3);bar(real(S));ylabel('节点注入有功');grid on;subplot(2,2,4);bar(Siz);ylabel('支路首端无功');grid on;。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《电力系统潮流上机》课程设计报告院系：电气工程学院班级：电088班________学号：01学生姓名：刘东昇________指导教师：张新松设计周数：两周_________日期：2010年12月25日一、课程设计的目的与要求目的：培养学生的电力系统潮流计算机编程能力，掌握讣算机潮流讣算的相关知识要求：基本要求：1.编写潮流计算程序:2.在计算机上调试通过：3.运行程序并il•算岀正确结果：4.写出课程设计报告二、设计步骤：1.根据给泄的参数或工程具体要求（如图），收集和査阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。
勺节点1为平衡节点，节点2、3、4、5为PQ节点匚2.在给定的电力网络上画出等值电路图。
3.运用汁算机进行潮流计算。
4.编写设汁说明书。
三、设计原理1.牛顿•拉夫逊原理牛顿迭代法是取xO之后，在这个基础上，找到比xO更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。
牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f（x） = 0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。
电力系统潮流il•算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡肖点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。
为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知肖点赋电压初值，一般为额泄电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。
牛顿一拉夫逊迭代法的一般步骤：（1）形成各节点导纳矩阵Y。
（2）设个节点电压的初始值U和相角初始值e还有迭代次数初值为0。
（3） II-算各个节点的功率不平衡量。
（4〉根据收敛条件判断是否满足，若不满足则向下进行。
（5）计算雅可比矩阵中的各元素。
（6）修正方程式个节点电压（7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。
（8）汁算平衡节点输出功率和各线路功率2.网络节点的优化1）静态地按最少出线支路数编号这种方法由称为静态优化法。
在编号以前。
首先统汁电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。
当由n个节点的出线支路相同时，则可以按任意次序对这n个节点进行编号。
这种编号方法的根据是导纳矩阵中，出线支路数最少的肖点所对应的行中非零元素也2 ）动态地按增加岀线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统讣出来的，在编号过程中认为固泄不变的，事实上，在廿点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化（增加，减少或保持不变）。
因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个盯点进行编号，就可以预期得到更好的效果，动态按最少岀线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中％巧点出线支路数目的变动情况。
3.MATLAB编程应用Matlab是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。
由于使用Matlab编程运算与人进行科学汁算的思路和表达方式完全一致，所以不像学习髙级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一髙效的科研助手。
四、设计内容1.设计流程图2.程序clear;clc%重新编号，把原题中的节点1,2,3,4,5重新依次编号为54,2,3,4.其中1-4号为PQ节点，5号为平衡节点y 二0;%输入原始数据，求节点导纳矩阵y (i,2)=i/+； y(13)=1/+; y ⑴ 4)=i/+；y ⑴ 5)=1/+;y(2/3)=l/+；y(2/5)=l/+;y(3z4)=l/+;y(4,5)=0;for i=l:5for j=i:5y(j』=y(i,j);endendY=0;%求互导纳for i=l:5for j=l:5if r=jY(i,j)=-y(i,j)；endendend%求自导纳for i=l:5Y(i/i)=sum(y(i/:));endY%Y为导纳矩阵G=real(Y);B=imag(Y);%原始节点功率S(l)=+;S(2)=賦初值U 二ones(l,5);U(5)=;e=zeros(：U5);ox=on€s(8j);fx=ones(8j);count=0 %计算迭代次数while max(fx)>le-5for i=l:4for j=l:4H(i,j)=0;N(i』=0;M(i』=0;L(i,j)=0;oP(i)=0;oQ(i)=0；endendfor i=l:4for j=l:5oP(i)=oP(i)-U(i)#U(j)*(G(ij)*cos(e(i)-e(j))+B(iJ)#sin(e(i)-e(j))); oQ(i)=oQ(i)-U(i)*UO)*(G(i,j)*sin(e(i)-e(j))-B(i z j)*cos(e(i)-e(j)));endoP(i)=oP(i)+P(i); oQ(i)=oQ(i)+Q(i);endfx=[oP/oQ]';%求雅克比矩阵%当中时候求H,NML如下:for i=l:4for j=l:4 ifr=jH(i,j)=-U(i)*U0)*{G(i,j)*sin(e(i)-e(j))-B(i/j)*cos(e(i)-e(j)))；N(i/j)=-U(i)*U(j)*(G(ij)*cos{e(i)-e0))+B(iJ)*sin(e(i)-eO))); L(i,j)=H(i,j);M(iJ)=-N(iJ);endendendH,NML%当冃时HNM丄如下:for i=l:4for j=l:5if r=jH{i,i)=H(i z i)+U(i)*U0)*{G{i/j)*sin(e(i)-e(j))-B(i z j)*cos (e(i)-e(j)));N(i/i)=N(i,i)-U(i)*U0)*(G(i/j)*cos(e(i)-e(j))+B(i/j)*sin(e(i)-e(j))); M(i/i)=M(i,i)-U(i)*UO)*{G(i,j)*cos(e(i)-e(j))+B(i/j)*sin(e(i)-e(j))); L(i,i)=L(i/i)-U(i)*UO)*(G(i,j)*sin(e(i)-e(j))-B(i z j)*cos(e(i)-e(j))); endendN(i/i)=N(i z i)-2*(U(i))^2*G(i/i);L(i,i)=L(i/i)+2*(U(i))A2^B(i/i)；endJ=(H,N;M,L] %J为雅克比矩阵ox=-((inv(J))*fx);for i=l:4oe(i)=ox(i); oU(i)=ox(i+4)*U(i);endfor i=l:4e(i)=e(i)+oe(i); U(i)=U(i)+oU(i);endcount=co un t+1;endOX,U QCO unt%求节点注入的净功率i=5;for j=l:5P(i)=U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)))+P(i); Q(i)=U(i)*UO)*(G(i,j)*sin(e(i)-e(j))-B(iJ)*cos(e(i)-e(j)))+Q(i); endS(5)=P(5)+Q(5)*sqrt(-l);S%求节点注入电流l=Y*U'3.运行结果Y值：(p New to MATLAB? Vifatch this Video, see Demos, or read Gettina Started.10. 8333 -32.50001 -1.6667 + 5-OOOOi -1. 6667 + S.OOOOi -2.5000 + 7.5000i -5.0000 +15.0000i-1. 6667 + 5.0000i 12.9167 -38.7500i -10. 0000 +30. 00001 0 ・ L.2500 + 3.75001 -1. 6667 + 5.00001 -10.0000 +30.0000i 12. 9167 -38.75001 -1.2500 + 3. 7500i 0-2. 5000 + 7.5000i 0 -1.2500 + 3. 7500i 3.7500 -11.2500i 0-5.0000 +15.00001 -1.2500 + 3.7500i 0 0 6.2500 -18. 7500i 迭代过程:count =0 5. 0000 5. OOOO7. 50005. 0000 0 30. 0000 05. 0000 30.0000 0 3. 75007. 5000 0 3. 7500 00 1. 6667 1. 6667 2・ 50001.6667 0 10. 0000 01. 6667 10.0000 0 1. 25002. 5000 O 1.2500 00 -1. 6667 -1. 6667 -2. 50001.6667 0-10. 000001. 6667 -10.0000 0 一 1・25002. 5000 0 -1. 2500 00 5. 00005. 0000 7. 50005. 000030. 0000 05.0000 30.00000 3.75007.50003. 7500 0-33.4000 5. 0000 5. 0000 7. 5000 -10. 5333 1. 66671. 66672. 50005. 0000 -38. 9750 30. 0000 0 1. 6667 -12. 8417 10. 0000 0 5・ 0000 30. 0000 -38. 7500 3・ 7500 1. 6667 10. 0000-12. 9167 1. 2500 7. 5000 03. 7500 -11. 2500 2. 50000 1. 2500 -3. 7500 11. 1333-1.6667-1.6667 -2.5000 -31.60005.00005.0000 7. 5000 -1.6667 12.9917 -10.0000 05.0000 -38.5250 30. 0000 0-1.6667 -10.0000 12.9167 -1.2500 5.0000 30.0000 -38.7500 3. 7500 -2.5000-1.25003. 75007.50003.7500-11. 2500M =-1・ 5572-L ・ 5235-2. 1511 一1・ 97010 一 1O ・ 1 13*7 0一1・ 9985 -IO ・ 48030 -1. 2203一3・ 2290 -L ・ 3384 00 5. 35995. 36228. 03965. 2222 0 30. 95205・ 2038 30. 8298 0 3. 85767. 7223 0 3. 8183 0-35. 0648 5. 3599 5. 3622 8. 0396-12. 0310 1. 5572 1. 5235 2. 15115. 2222 -40. 0794 30. 9520 0 1. 9701 -12. S627 10. 1137 0 5. 2038 30. 8298 -39. 8913 3. 85761. 9985 10. 4803 -12. 86821. 2203 7. 7223 03. 8183 -11.5406 3. 1029 01. 3384. -3. 1996 1 1. 5373-1. 5572 -1. 5235 -2. 1511-35. 6400 5. 3599 5. 36228. 0396-1. 9701 L3. 7720 -10. 1 1375. 2222 -39. 8247 dO ・ 9520 0-1. 9986 -L0. 4S03 13. 6991 -1. 2203 5. 203S 30. 8298 -39. S105 3. S576-3. 1029 0 -1. 33844. 44127. 72230 3. 8183 -11. 38185 35995・ 2222 5・ 2038 7・ 722336 32985・ 3622 30・ 95203・ & 1838.03963. 8576O1・ 9701 1・ 9985 3・ 10291. 5572 4803L ・ 5235 10・ 1 137 0L ・ 33842. 1511O1. 22035. 29045・ 2J912 7・ 924(3峪.1585 O30・ 5426石.1397 30・ 42673・ 80L8T ・ 62 17 O 3. 7C44 O0 1・ 54361・ 5 1 122・ 13681・ 9393 O9・ 93761・ Q6DS 1 O. 3300O1 ・ 2OE5O3・ 0463O1・ 3X710 一1・ 5436一1・ 5 112一2・1368 一 1・ 0303O —Q. QS76O一 1・ 9658 -10・ 3366一1・ 2080一3・ 0453O 一1・ 3 1715. 29045. 2912 7. 92465. 158530. 54265. 1397了0. 42673. 80187. 62173. 7644 0-1・ 5436 —1.・ 511 1 -2-136S 一 1・ 93910 -9・ 9867 0 一1. 9655 -10.3343 0 一1. 2049一3・ 04.4S 0 -L ・ 3169 0行.28995・ 2907 92375・ 1S80 030. 5395O5・ 1392 30・ 42360 3. 80137・ 62100 3・ 7640 0-34.7163 5. 2904 5. 2912 7. 9246 5. 1585 -39.5849 30. 5426 05. 1397 30. 4267 -39. 36813. 8018 7・ 62173・ 76<14-LI ・ 3S611 1・ 4384 -1.5436-1. 5112 -2. 1368 -1. 9393 13.5950 -9. 9876 0-1. 9658 -10.335513. 5062-1. 2050-3. 0453-1. 31714. 3623H =5・ 2899 5・ 29077. 92375・ 1580 036 53955・ 1392 30・ 42363.80137・ 62103. 7640N =0 1・ 5436 L ・ 5111 2.13681・ 9391 0 9・ 9867 0 1. 9555 1O.3343 0 1-2049 3・ oaa8 0 L ・ 3169 0-11.8394 1. 5436 1. 51122. 13681.9393 -12. 6950 9・ 9876 0 1.9658 10.3355 -12. 70591. 20503・ 0<1531. 3171 7・ 1619 -36.1173 5. 29045. 2912 7.9246 氐 1585 -39. 2853 30.5426 0 5. 1397 30. 4267-39. 2683 3. 80187.62173. 7644 -1 1. 1866-34.7136 5. 2899 5.2907 7. 9237 -H.837S 1.5436 1. 5111 2.13685.1580 -39.5812 30.5395 0 1. 9391 -12.6937 9. 9867 05.1392 30. 4236 -39.3642 3.8013 1. 9655 10.3343 -12. 7047 1.20497.6210 0 3.7640 -11.3850 3. 0448 0 1.3169 -3.161711.4379 -1. 5436 -1.5111 -2.136S -35. 1136 5.2899 5・ 2907 7.9237-1.9391 13. 5937 -9.9867 0 5. 1580 ・ 39.2812 30. 5395 0-1.9655 -10. 3343 13.5047 -1.2049 5. 1392 30.4236-39. 2642 3.8013-3.04480 -1.3169 4.3617 7. 62100 3. 7640 -11.1850电压值:ox =1.0 e-OO8 *-O.0004O.08350.1153O-2L57-O.2645-O.3080-O-3215-0.39281.0365 1・ 0088 1・ 0073L・ 00 16 1・ 060。