用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应急演练目的①评估我矿水灾事故的应急准备状态，发现并修改我矿水灾事故专项应急预案和执行程序中存在的缺陷和不足；②评估我矿在发生水灾事故时的应急能力，识别处理水灾事故的资源需求，澄清相关单位和人员的应急职责，改善水灾事故应急救援中的组织协调问题；③检验应急响应人员对水灾事故应急预案及执行程序的了解程度和实际操作技能；同时，通过调整演练难度，进一步培训和提高应急响应人员的业务素质和能力；④提高全员安全意识。
- 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。