电力系统课设报告

1、电力系统分析课程设计电力系统分析课程设计报告专业班级：电气116班 学生学号： 3110421175 学生姓名： 刘万成 指导教师：段建东教授 平时表现(10%)验收答辩(50%)设计报告(40%)总评成绩 2014 年 秋季 学期起止时间：2015年1月5日至2015年1月16日42 目 录 一、 设计目的. .2二、 设计任务. .2三、 基础理论3.1不对称短路故障概述.23.2不对称短路的计算条件及近似33.3系统标幺值计算.33.4不对称短路的计算方法 3.4.1对称分量法. .4 3.4.2简单不对称短路计算步骤.4 3.4.3正序等效定则.4 3.4.4节点导纳矩阵的计算方法.7

2、四、 程序设计4.1计算程序原理.84.2变量说明.94.3程序流程图.9五、 算例分析. .11六、 收获建议. .14 附件一：程序源代码.16 附件二：系统算例.31 附件三：输出结果.32附件四：参考文献.42一、 设计目的1、 加深理解并巩固电力系统基础知识；2、 掌握电力输电、配电系统的计算方法；3、 熟悉电力系统计算程序基本编写过程；4、 为今后从事与电力系统相关的设计、运行、管理及科学研究工作打下良好的基础。二、 设计任务1、 掌握电力系统不对称短路计算，理解专业理论、计算机算法；2、 用MATLAB或C语言编写一个计算程序并调试通过；3、 分别对教材例题与数个标准测试系统（I

3、EEE节点系统）解出算例结果并分析。三、 基础理论3.1不对称短路故障的概述在电力系统运行过程中，时常发生故障，其中大多数是短路故障。所谓短路：是指电力系统正常运行情况以外的相与相或相与地（或中性点）之间的连接。除中性点外，相与相或相与地之间都是绝缘的。电力系统可分为三相短路，单相短路，两相短路和两相接地短路等。三相短路的三相回路依然是对称的，故称为对称短路，其他的几种短路的三相回路均不对称，故称为不对称短路。电力系统运行经验表明，单相短路占大多数。产生短路的主要原因是电气设备载流量部分的相间绝缘或相对地绝缘被损坏，主要表现在如下几个方面：元器件损坏，例如绝缘材料的老化，设计，安装及维护不良所

4、带来的设备缺陷发展正短路；气象条件恶化，例如雷电造成的闪络放电，架空线路由于大风或冰雪灾害引起杆塔倒塌；违规操作，例如运行人员在线路检修后未拆除地线就加电压。不对称短路故障会引起很大的危害，主要表现在以下几个方面：短路回路中的短路电流值大大增加，电流的热效应会引起导体或其绝缘的损害；另一方面；导体也会受到很大的电动力的冲击；短路还会引起电网中电压降低，特别是靠近短路点处的电压下降得最多，结果可能使部分用户的供电受到破坏；破坏了系统的稳定性；不对称短路所引起的不平衡电流产生的不平衡磁通，造成对通信系统的干扰，甚至危及设备和人身安全。短路问题是电力技术方面的基本问题之一。短路计算用于合理选择电气接

5、线、选用有足够热稳定度和动稳定度的电气设备及载流导体、确定限制短路电流的措施、在电力系统中合理地配置各种继电保护并整定其参数等。实际电力系统中的短路故障大多数是不对称的，为了保证电力系统和它的各种电气设备的安全运行，必须进行各种不对称故障的分析和计算。3.2不对称短路的计算条件和近似（1）故障前，正常情况下系统视为空载；（2）故障点正常运行电压Uf=1，省去潮流计算；（3）系统中所有变压器均为Yd11接线；（4）不计各元件的电阻，即R=0；（5）各元件的正序和负序电抗相等，即X1=X2；（6）不计线路电容电流的影响；（7）系统三相参数相同；3.3系统元件标幺值计算先选取基准值则主要元件的标幺值

6、计算公式如下：（1） 发电机：（2） 变压器：（3） 输电线路：（4） 电抗器：3.4不对称短路的计算方法3.4.1对称分量法电力系统中不对称短路故障包括两相短路、单相接地短路以及两相接地短路三种类型。分析计算不对称短路的方法很对，目前实际中用的最多的和最基本的方法仍是对称分量法。对电力网络的进行不对称短路故障分析时，通常从短路点f处把网络分成两部分，对两部分分别进行讨论。用对称分量法吧该电路中任一组不对称相量分解成正序、负序、零序三组对称分量，这三组对称分量是各自独立的，即每一相序的电压只能产生本相序的电流。于是可得到三个各自独立的序网，以及三个序网方程；另一部分是短路点对地的不对称阻抗电路

7、，反映该电路特点的是各种故障类型故障点的边界条件，即反映该不对称短路类型短路点的电压和电流的方程式。将两部分电路组合起来，即将三个序网方程与短路点边界条件联立起来，解出短路点的各序电流和各序电压，进而求得各节点的三相电压和各支路的三相电流。对称分量法如下：设为三相系统中的任意一组不对称相量，可以用一组对称的相量表示， 其中同样，由三个不对称的相量可以唯一地分解成三组对称的相量（即对称分量）：正序分量、负序分量和零序分量，公式为3.4.2简单不对称短路计算步骤电力系统发生不对称故障时，无论是单相接地短路，还是两相短路，两相短路接地，只是短路点的电压，电流出现不对称，利用对称分量法将不对称的电流电

8、压分解为三组对称的序分量，由于每一序系统中三相对称，则在选好一相为基准后，每一序只需要计算一相即可，用对称分量法计算电力系统的不对称故障。其大概步骤如下：（1） 针对正常部分，找正序、负序和零序网络，计算电力系统各个元件的序阻抗；（2） 针对故障点处，找边界条件（序分量），制定电力系统的各序网络；（3） 由各序网络和故障，联立3个序网方程与3个边界条件；（4） 从联立方程组解出故障点电流和电压的各序分量，将相对应各序分量进行Park反变换，以求得故障点的各相电流和电压；（5） 计算各序电流和各序电压在网络中的分布，进而求得各指定支路的各相电流和指定节点的各相电压。3.4.3正序等效定则依据上述

9、步骤，不难求出三种不对称短路的电流。分析结果，可以看出这三种情况下短路电流的正序分量的计算式和三相短路电流在形式上很相似，只是阻抗为为附加阻抗。在单相短路时附加阻抗为和（或）的串联；两相短路时附加阻抗为（或）；两相接地短路时为和（或）的并联。对于任一种不对称短路，其短路电流的正序分量可以利用图3-4示出的正序增广网络计算。 根据不同的故障，可分别利用表3-4所列公式计算故障处各序电流、电压，进而可合成得三相电流、电压。表3-4 各类故障处各序电流、电压计算公式故障种类故障端口各序电流公式故障端各序电流公式单相短路两相短路两相接地短路求解得出故障处的各序电流、电压后，可根据以下公式计算网络中任一

10、节点电压和任一支路电流。任一节点电压的各序分量为： （3-4-1）式中：为正常运行时该节点的电压；为各序网阻抗矩阵中与故障点f相关的一列元素。任一支路电流的各序分量为： （3-4-2） 若进过变压器（Yd11接线），则正序分量和负序分量的相位要发生变化，即对于正序分量三角形侧电压、电流较星形侧超前30°（即11点钟）或落后330°，负序分量三角形侧电压、电流较星形侧落后30°或超前330°，零序分量不流通。其关系式可表达为： （3-4-3） （3-4-4）3.4.4节点导纳矩阵的计算方法计算f节点短路电流的关键，就是求各序网形成阻抗矩阵的第f列元素。各序

11、阻抗第f列元素等于在f节点加单位电流，其他节点注入电流都为零时，各节点的电压值。故可用节点导纳矩阵为系数矩阵的导纳型网络节点电压方程来求阻抗矩阵的第f列元素。节点导纳矩阵的生成公式为：（1）对角线元素： （3-4-5）（2）非对角线元素： （3-4-6）求解下列方程 （3-4-7）可得即为。将式（3-4-7）简写成YU=I （3-4-8）Y阵是个非奇异的对称矩阵，按照矩阵的三角分解法，Y可表示为 （3-4-9）式中：D为对角阵；L为单位下三角阵；R为单位上三角阵；L和R互为转置阵。式（3-4-9）说明Y阵分解为单位下三角阵，对角阵和单位上三角阵（即单位下三角阵的转置）的乘积。这些因子矩阵元素的

12、表达式为 （3-4-10） 式中：d、l和r为矩阵D、L和R的相应元素。由于L和R互为转置，只需算出其中一个即可。将式（3-4-9）代入式（3-4-8）得 （3-4-11）式（3-4-11）可分解为以下的三个方程，并以次求解，得 （3-4-12） （3-4-13） （3-4-14）即由已知的节点电流相量I求W，由W求X，最后由X求得节点电压相量U，即一列阻抗元素值。四、 程序设计4.1计算程序原理（1）如果要求准确计算故障前的运行情况，则需要进行潮流计算，在近似使用计算中，对短路故障可假设节点；（2）这里采用形成节点导纳矩阵的方法。发电机的正序电抗用，可计算故障后瞬时量。发电机的负序的电抗近似

13、等于，当计算中不计负何影响时，在正、负序网络中不接入负荷阻抗。如果计及负荷影响时，负荷的正序阻抗可通过其额定功率和电压计算。负何阻抗很难却确定，一般取（以负荷额定功率为基准）。负荷的中性点一般不接地，零序无通路；（3）形成三个序网的节点导纳矩阵后，利用公式（3-4-11）求得，即为。对于短路故障，只需令（其余节电电流均为零），分别应用三序因子表求解一次所得电压，即为三序网和f点有关节点阻抗；（4）根据不同短路故障，可分别应用表3-4中相应公式计算故障处各序电流、电压，进而合成得三相电流，电压；（5）计算网络中任一点电压时，负序和零电压只须计算由故障点电流引起的电压。对于正序则还需要加上正常运行

14、时的电压。4.2变量说明因为C语言不能只能直接进行复数计算，而不对称短路计算时要用到复数运算，故考虑把计算分成实部和虚部分别进行。下面以计算任意节点m的三相电压为例进行说明。令其中：、表示三序电压取实部，、表示三序电压取虚部 。4.3程序框图开始4.3.1主程序流程图 输入数据形成节点导纳矩阵分解成节点阻抗矩阵选择故障类型两相接地短路单相接地短路两相短路任一节点电压电流输出数据结束4.3.2部分子程序流程图 五、 算例分析例5-1：一简单电力系统如图5-1所示，系统等值电路也已给出如图5-2所示。又已知三台发电机中性点均不接地，三台变压器均为Ynd接线（发电机侧为三角形），经试验测得三条输电线

15、路的零序电抗均为0.20（以60MVA为基准值）。要求计算节点分别发生单相短路接地、两相短路和两相短路接地时故障处的短路电流和电压（基波分量初始值）。 图5-1 电力系统线路图图5-2 等值电路图分析：1、手算上述算例计算三个序网络对故障点的等值阻抗： 故障处电流：(1) a相短路接地：(2) b、c两相短路：(3) b、c两相短路接地： 故障处电压：（1） a相短路接地：序电压： 相电压： （2） b、c两相短路：序电压： 相电压： （3） b、c两相短路接地：序电压： 相电压： 2、计算机计算结果（1）a相接地短路： 短路点序电压值： Uf1=0.684758 Uf2=-0.315242

16、Uf0=-0.369516短路点相电压值： Ua3=-0.000000 Ub3=-0.554274-j0.866000 Uc3=-0.554274+j0.866000 短路点序电流值： If1=-j3.108021 If2=-j3.108021 If0=-j3.108021短路点故障相电流值： If=-j9.324062 （2）b、c相短路： 短路点序电压值： Uf1=0.500000 Uf2=0.500000 Uf0=0.000000短路点相电压值： Ua3=1.000000 Ub3=-0.500000+j0.000000 Uc3=-0.500000+j0.000000 短路点序电流值： I

17、f1=-j4.929577 If2=j4.929577 If0=j0.000000 短路点故障相电流值： If=-j8.538279 （3）b、c相短路接地： 短路点序电压值： Uf1=0.350493 Uf2=0.350493 Uf0=0.350493短路点相电压值： Ua3=1.051480 Ub3=0.000000+j0.000000 Uc3=0.000000-j0.000000 短路点序电流值： If1=-j6.403588 If2=j3.455567 If0=j2.948021 短路点故障相电流值： If=j1.628048比较手算结果和计算机计算结果，数据全部一致。六、 收获建议6

18、.1设计中遇到的问题及解决方法（1）问题：拿到题目不知道从何下手。 解决方法：先分析算例，找到已知条件和要计算的未知量，然后进行变量定义。通过看教材例题，初步确定程序设计方案，列出主程序流程图，再分模块编写程序。（2）问题：不能把数据通过文本文档的形式输入。 解决方法：查看C语言教材，上网查阅相关资料，了解C语言是如何进行文件的输入输出操作的。学习fprintf和fscanf语句的使用方法。（3）问题：生成节点导纳矩阵后，不能运用C语言直接进行矩阵的求逆来计算节点阻抗矩阵。 解决方法：查阅相关资料，了解到形成节点阻抗矩阵可以运用高斯消去法、三角分解法和因子表法。为了避免每次重复对节点导纳矩阵作

19、消去运算，一般不采用高斯消去法求解节点阻抗矩阵。（4）问题：C语言不能直接进行复数计算。解决方法：把复数运算分为实部和虚部分别进行。但计算时要特别注意j*j=-1，而且计算比较复杂，必须考虑周全。（5）问题：程序只能进行A相单相接地短路，B、C相短路和B、C相短路接地，不能进行其他相的故障计算。 解决方法：计算时把A相看作特殊相，然后再通过故障相判定程序，输出显示时再考虑不同特殊相。6.2心得总结电力系统发生不对称短路后，由于短路点对地故障支路的不对称，使得整个网络电流电压三相不对称。本次课设解决不对称短路的问题核心是对称分量法。我们平时手算时，就是把该网络分解为正、负、零序三个对称序网，这三

20、组对称序分量可分别按对称的三相电路分解，然后将其结果叠加起来。求解不对称短路，首先应该计算各原件的序参数和画出等值电路，然后制定各序网络，根据不同的故障类型，确定出以相分量表示的边界条件，进而列出以序分量表示的边界条件，按边界条件将三个序网联合成复合网，由复合网求出故障处各序电流和电压，进而合成三相电流电压。运用计算机程序法，就是通过计算机形成三个序网的节点导纳矩阵，然后利用高斯消去法或三角分解法，通过相应公式对它们进行数据运算，即可求得故障端点的等值阻抗矩阵，最后根据故障类型选取相关公式计算故障处各序电流电压，进而合成三相电流电压。本次课程设计我选择C语言编程，主要是考虑到之前接触的C语言比

21、较多，感觉C语言有一定的功底。而在着手编写程序时才发现，由于本次课程设计用到的计算特别多，而且复杂，例如复数的计算和矩阵的计算，MATLAB对于这些计算都比较简单，但C语言却很难实现。此次编程时遇到的最大的困难就是节点阻抗矩阵的生成，指导书上给的方法是三角分解法，而我按照此思路编写程序，总是不能生成正确的节点阻抗矩阵。最后查阅资料，了解到还有一种高斯消去法，虽然高斯分解法会重复计算，但计算原理比较简单，最后采用高斯消去法。本次设计最大的特色就是用简单的程序解决了不同相短路故障的问题，不再拘泥于只能A相为特殊相。设计中，故障计算部分，始终把A相当做特殊相，最后在输出显示之前，先做一次故障相的判定

22、，找到特殊相，把故障计算得到的A相的数值重新赋给特殊相，在进行输出显示。另外，经过变压器的支路，程序还可以实现变压器两侧的相位变化。当然，本次设计还存在很多不足之处。例如自己编写的用三角分解法求节点阻抗矩阵的程序始终通不过，而且发现不了错误之处；还有，本次设计的程序只能让故障发生在各个节点处，不能实现发生在输电线路上，而现实中任何点发生故障都是可能的；再者，也没用做出漂亮的人机交互界面；最后，本次设计用了很多假设，省去了很多工作，如潮流计算。课程设计虽然接触，但学习工作远未结束，我还要继续努力，完善程序，学习更多的知识。最后，要感谢段老师以及各位学长两周来的辛勤辅导，还有各位同学的帮助。没有他

23、们的帮助，我很难在两周内完成本次设计。通过本次设计，不仅让我对原来学过的知识有了系统的复习，还加强了我的编程能力，使我对以后的研究生生涯有了更大信心。附件一/*不对称短路近似计算约定：1、故障前，正常情况为空载；2、Uf=1,省去潮流计算；3、系统所有变压器均为Yd11接线；*/#include<stdio.h>#include<math.h>#define M 100 /最大值路数#define N 100 /最大节点数/* 全局变量定义*/int h, n, /节点数 m, /正序支路数v, /零序支路数 L, /故障类型 F; /故障点char k1,k2,k3;

24、 /故障相float zp; /过渡电阻float b1MM,b0MM; /非零导纳元素 static float Z1MM,Z0MM; /阻抗矩阵数组static float If1,If2,If0; /故障口电流static float Uf1,Uf2,Uf0; /故障口电压static float Zf1,Zf2,Zf0; /故障口阻抗static float If; /故障相电流static float U0_sM,U1_sM,U2_sM, /序电压 U0_xM,U1_xM,U2_xM, Ua_sM,Ua_xM,Ub_sM,Ub_xM,Uc_sM,Uc_xM,/相电压的实部、虚部 Ia

25、_sMM,Ia_xMM,Ib_sMM,Ib_xMM,Ic_sMM,Ic_xMM, I1_sMM,I2_sMM,I0_sMM,/序电流 I1_xMM,I2_xMM,I0_xMM;/* 结构体定义*/struct zheng int h; / h=1为正序 int d; / 支路类型(d=1表示变压器支路，d=0表示无变压器支路) int p1,p2; / 支路前后的两个节点 float x; / 支路的电抗zhengM; / 正序结构体 struct ling int h; / h=0为零序 / int d; / 支路类型 int p1,p2; / 支路前后的两个节点 float x; / 支路

26、的电抗 lingM; / 零序结构体 FILE *fp1,*fp2; /文件指针定义 /* 文件操作*/ void Read_data() int i,j; /* 读取原始数据 */ fp1=fopen("input.txt","r"); / 打开input.txt文件 if(fp1=NULL) printf(" 文件为空，重新检查输入！n"); fscanf(fp1,"%d,%d,%dn",&n,&m,&v); /输入节点数，正序支路数,零序支路数 j=1; do for(i=1;i<

27、;=m+v;i+) fscanf(fp1,"%d",&h); if(h=1) /正序阻抗 fscanf(fp1,",%d,%d,%d,%fn",&zhengi.d,&zhengi.p1,&zhengi.p2,&zhengi.x); j+; if(h=0) /零序阻抗 fscanf(fp1,",%d,%d,%fn",&lingi.p1,&lingi.p2,&lingi.x); j+; while(j<=m+v); fclose(fp1); if(fp2=fopen(&

28、quot;output.txt","w")=NULL) printf(" 文件为空，重新检查输入!n"); /*程序原始数据输出*/ fprintf(fp2," 原始数据n"); fprintf(fp2,"n"); fprintf(fp2," 节点数:%2d 正序支路数:%2d 零序支路数:%2dn", n,m,v); fprintf(fp2,"n正序数据：n"); fprintf(fp2,"支路类型 节点号 节点号 支路阻抗n"); for(i

29、=1;i<=m;i+) fprintf(fp2," %d %2d %2d X1=%f n",zhengi.d,zhengi.p1,zhengi.p2,zhengi.x); fprintf(fp2,"零序数据：n"); fprintf(fp2,"节点号 节点号 支路阻抗n"); for(i=m+1;i<=m+v;i+) fprintf(fp2," %2d %2d X0=%f n",lingi.p1,lingi.p2,lingi.x); /* 形成阻抗矩阵*/void Form_Y()int a,b,i,j

30、,k;for(i=0;i<=m;i+) for(j=0;j<=m;j+) b1ij=0;b0ij=0; Z1ij=0;Z0ij=0; for(i=1;i<=m;i+) a=zhengi.p1;b=zhengi.p2; b1ab=-1/zhengi.x; /正序非零元素b1ba=-1/zhengi.x;for(i=m+1;i<=m+v;i+)a=lingi.p1;b=lingi.p2;b0ab=-1/lingi.x; /零序非零元素b0ba=-1/lingi.x;/*节点电压法形成节点导纳矩阵*/ for(i=1;i<=n;i+)for(j=1;j<=n;j+

31、)Z1ij=0;for(i=1;i<=n;i+)for(j=0;j<=n;j+)Z1ii=Z1ii+b1ij; /正序自导纳for(i=1;i<=n;i+)for(j=1;j<=n;j+)if(i!=j)Z1ij=-b1ij; /正序互导纳for(i=1;i<=n;i+)for(j=1;j<=n;j+)Z0ij=0;for(i=1;i<=n;i+)for(j=0;j<=n;j+)Z0ii=Z0ii+b0ij; /零序自导纳for(i=1;i<=n;i+)for(j=1;j<=n;j+)if(i!=j)Z0ij=-b0ij; /零序互导

32、纳/*以文件的形式输出导纳矩阵*/ fprintf(fp2," nn正序导纳矩阵：n"); for(i=1;i<=n;i+)fprintf(fp2,"n ");for(j=1;j<=n;j+) if(Z1ij>=0)fprintf(fp2,"j%f ",Z1ij); elsefprintf(fp2,"-j%f ",-Z1ij); fprintf(fp2," n n零序导纳矩阵：n");for(i=1;i<=n;i+)fprintf(fp2,"n ");

33、for(j=1;j<=n;j+) if(Z0ij>=0)fprintf(fp2,"j%f ",Z0ij); else fprintf(fp2,"-j%f ",-Z0ij);fprintf(fp2,"nnn"); /* 导纳矩阵变换成阻抗矩阵(高斯消去法)*/void Form_Z()int i,j,k;float Z;float ZZ5050; /求正序阻抗矩阵 for(i=1;i<=n;i+) for(j=1;j<=n;j+) ZZij=0; for(k=1;k<=n;k+) ZZkk=-1/Z1kk;

34、 for(i=1;i<=n;i+) for(j=1;j<=n;j+) if(i!=k&&j!=k) ZZik=-Z1ik/Z1kk; ZZkj=-Z1kj/Z1kk; Z=Z1ik*Z1kj/Z1kk; ZZij=Z1ij-Z; for(i=1;i<=n;i+)for(j=1;j<=n;j+)Z1ij=ZZij; /求零序阻抗矩阵 for(i=1;i<=n;i+) for(j=1;j<=n;j+) ZZij=0; for(k=1;k<=n;k+) ZZkk=-1/Z0kk; for(i=1;i<=n;i+) for(j=1;j&l

35、t;=n;j+) if(i!=k&&j!=k) ZZik=-Z0ik/Z0kk; ZZkj=-Z0kj/Z0kk; Z=Z0ik*Z0kj/Z0kk; ZZij=Z0ij-Z; for(i=1;i<=n;i+) for(j=1;j<=n;j+) Z0ij=ZZij; fprintf(fp2,"正序阻抗矩阵：n");for(i=1;i<=n;i+)fprintf(fp2,"n ");for(j=1;j<=n;j+)if(Z1ij>=0)fprintf(fp2,"j%f ",Z1ij);els

36、efprintf(fp2,"-j%f ",-Z1ij);fprintf(fp2," n n零序阻抗矩阵：n");for(i=1;i<=n;i+)fprintf(fp2,"n ");for(j=1;j<=n;j+)if(Z0ij>=0)fprintf(fp2,"j%f ",Z0ij);elsefprintf(fp2,"-j%f ",-Z0ij);fprintf(fp2,"nnn"); /* 串联型故障*/void CHUAN() float zp;printf

37、(" 请输入发生单相接地短路故障的节点：");scanf("%d",&F);printf(" 请输入发生单相接地短路故障的相：");scanf("%c",&k1);k1=getchar(); printf(" 请输入过渡电阻zp："); scanf("%f",&zp); Zf1=Z1FF; Zf2=Z1FF; Zf0=Z0FF; If1=-1/(Zf1+Zf2+Zf0+3*zp); /阻抗有个j If2=If1; If0=If1;If=3*If1; U

38、f1=1+If1*Zf1; /本身有个j，因为j*j=-1,所以减号变加号 Uf2=If2*Zf2; Uf0=If0*Zf0;/* 并联型故障*/void BING() float zp; if(L=2|L=3)printf(" 请输入发生两相短路或两相接地短路故障的节点："); scanf("%d",&F);printf(" 请输入发生两相短路或两相接地短路故障相：");scanf("%c",&k2);k2=getchar();scanf("%c",&k3);k3=ge

39、tchar(); printf(" 请输入过渡电阻zp："); scanf("%f",&zp); Zf1=Z1FF; Zf2=Z1FF; Zf0=Z0FF;if(L=2)If1=-1/(2*Zf1+zp); If2=-If1; If0=0;If=sqrt(3)*If1; Uf1=1+If1*Zf1; Uf2=If2*Zf2; Uf0=0; if(L=3) If1=-1/(Zf1+Zf1*(Zf0+3*zp)/(Zf1+Zf0+3*zp); If2=-If1*(Zf0+3*zp)/(Zf1+Zf0+3*zp); /并联分流 If0=-If1*Zf1/(Zf1+Zf0+3*zp); If=sqrt(3)*sqrt(1-Zf1*(Zf0+3*zp)/(2*Zf1+Zf0+3*zp)*(2*Zf1+Zf0+3*zp); Uf1=1+Zf1*If1; Uf2=Zf1*If2;