电力系统及其自动化专业毕业设计

1、电力系统及其自动化专业设计人： 指导老师目录一、任务书34二、网络参数的计算59三、电力系统潮流的计算1°17附：pq分解法潮流计算源程序及结果183°四、短路电流的计算3146附：1、对称短路程序36372、不对称短路程序4346五、线路保护的整定4764六、附录：a：对称短路计算的系统接线图65b：零序网络图66c：潮流分布图67d：线路保护配置图68某220kv电网潮流计算及输电线路继电保护配置150mva 1-i(一) ：已知(1)系统最大运行方式为四台发电机满发和系统投入运行;系统最小运行方式为停两台发电机(f1,f3), 各负荷减半.系统各负荷及线路参数如图所示

2、，各变压器及发电机型号分别为：f 1f4: sf-100-40/8540pc= 1 oom wue= 13.8kv cos (p = 0.9r|=98.16%blb3:sfp7-150000/220220±2x2.5%/13.8ud=14%yn,dllb4:sfp7-150000/1001 10 ± 2 x 2.5%/13.8 ud=13%yn,dllb5b&sfpsz7-120000/220220±8x1.25%/115/10.5 yn,yno,dllui2=23.5%u23=13.3ui3=7.7%b9: rt = 3g, xt =1100(3)线路参

3、数如图中所示.(二) ：设计任务(1) 计算各元件阻抗标么值(sj=100mva,vj=vn)f：画出正序，负序，零序等效网络图；(2) 对系统进行潮流计算;(用c语言)(3) 对5、6、7、8点进行各种类型的短路电流计算;(qb)(4) 为线路67选择保护方式；(5) 对所选保护进行整定计算；(6) 对保护进行选型；(7) 画出保护原理图；(8) 书写设计说明书及准备答辩.(三) ：设计成果(1) 系统潮流分布图一张及短路电流表一份；(2) 线路保护配直图一-张；(3) 保护原理图一张；(4) 设计说明书一份。网络参数计算一、双绕组变压器的参数计算变压器的参数一般是指其等值电路中的电阻rt，

4、电抗xt，电导gt和电纳bt，变压器的变比k。根据铭牌 上所给的短路损耗aps，短路电压vs%,空载损耗ap。，空载电流1（）%。前两个数据由短路试验得到，用以确定 rt和xt；后两个数据由空载试验得到，用以确定gt和1. 电阻rt：变压器作短路试验时，将一侧绕组短接，在另一侧绕组丿施加电压，使短路绕组的电流达到额定 值由于此时外加电压较小湘应的铁耗也小，可以认为短路损耗即等于变压器通过额定电流时原、副方绕 组电阻的总损耗在电力系统计算中，常用变压器三相额定容量和额定线电压进行参数计算，则公式 为:rtq=apskwvn2kvx103/ sn2kva2. 电抗x：当变压器通过额定电流时，在电抗

5、上产生的电压降的人小，可以用额定电压的百分数表示, 对于大容量变压器，其绕纽电阻比电抗小得多贝j公式：xtq=vs%x vn2kvx 103/100/sn kva3. 电导gt：变压器的电导是用來衣示铁芯损耗的.由于空载电流相对额定电流來说是很小的，绕组小的 铜耗也很小，所以近似认为变压器的铁耗就等于空载损耗，则公式为：gts=4pokwx 10'3/vn2kv4. 电纳bt:变压器的电纳代表变压器的励磁功率.变压器空载电流包含有功分量和无功分量，与励磁功 率对应的是无功分量.由于有功分量很小，无功分量和空载电流在数值上几乎相等.bts=io%x snkvax 10-3/100/vn2

6、kv5. 变压比心：在三相电力系统计算中，变压器的变压比通常是指两侧绕组空载线电压的比值对于星形和三角形接法的变压器,变压比与原副方绕组匝数比相等;对于星三角形接法的变压器，变压比为原副方绕组匝数比的松倍.根据电力系统运行调节的耍求，变压器不一定丄作在主抽头上，因此，变压器运行中 的实际变比，应是工作时两侧绕组实际抽头的空载线电压z比.二、三绕组变压器的参数计算三绕组变圧器等值电路中的参数计算原则与双绕组变压器的相同，下面分別确定各参数的计算公式.1. 电阻r1r2r3：为了确定三个绕组的等值阻抗,耍有三个方程，为此，需耍有三种短路试验的数据三绕组变 压器的魚路试验是依次让一个绕组开路,按双绕

7、组变压器来作.通过查手册可得短路损耗分別为，则冇p s1 = 1/2（aps（1 -2）+ap s（3-1 ）-ap s（2-3）p s2= 1/2(ap s(1 2)+ap s(23) ap s(3-1)ap s3= 1/2( ap s(2-3)+ap s(3-1 厂 ap s(1 -2)求出各绕组的短路损耗后，便可导出双绕组变压器计算电阻相同形式的算式，bp：riq=psikwvn2kvx103/sn2kva2. 电抗x】,x2,x3：和双绕组变压器一样，近似地认为电抗上的电压降就等于短路电压.在给出短路电压力后, 与电阻的计算公式相似,各绕组的短路电压为v s 1 %= l/2(vs(l

8、 2)%+v s(3-1 )%-v s(23)%)v s2%= 1 /2( vs( 1 -2)%+v s(2-3)% v s(3-1)%)v s3%= 1/2(v s(23)%+v s(3-1 )%-v s(1 2)%)各绕组的等值电抗为:xia=vsi%x vn2kvx 103/100/ snkva3.导纳gt-jbt和变比ki2、ki3、k23：计算与双绕组变压器相同.三. 线路参数的计算.输电线路的参数有四个:反映线路通过电流吋产生有功功率损失效应的电阻;反映载流导线周围产牛磁场效应 的电感;反映线路帯电吋绝缘介质中产生泄漏电流及导线附近空气游离而产生有功功率损失的电导仮映带电导线 周围

9、电场效应的电容.输电线路的这些参数通常可以认为是沿全长均匀分布的,每单位长度的参数为!、x、g及b. 当线路长为l(km)时，r=rl;x=xl;g=gl;b=bl由于沿绝缘子的泄漏很小，可设g=0.四. 标么值的折算.建立电力网络和电力系统的数学模型,需解决标么值的折算问题.进行电力系统计算时，除采用有单位的阻抗、导纳、电压、电流、功率等进行运算外,还可采用没有单位的阻抗、 导纳、电压、电流、功率等的相对值进行运算.前者称有名制,后者称标么制.标么制之所以能在相当宽广的范围内 取代有名制,是由于标么制具有计算结果清晰、便于迅速判断计算结果的正确性、可大量简化计算筹优点。标么值=实际有名值(任

10、意单位)/基准值(与有名值同单位)对于直接电气联系的网络,在制订标么值的等值电路时，各元件的参数必须按统一的基准值进行归算由于各元件 的额定值可能不同，因此，必须把不同基准值的标么阻抗换处成统一基准值的标么值.现统一选定的棊准电压和基 准功率分别为v和s,则电抗的实际有名值换算成标么值，即在工程计算屮规电各个电压等级都以其平均额定电压v作为基准电压.根据我国现行的电压等级，各级平均额定 电压规定为 3.15, 6.3,10.5,15.75,37,115, 230,345, 525kv五具体计算及结果根据原始数据及资料显示，通过公式计算nj得：（选择基准t：sb=100mva,vb=vav） 1

11、、双绕组变压器:bbb2,b3 (型号 sfp7-150000/220 )(参数 ap()=14okw4ps=45okw,vs%=14%,io%=o.8,yndl 1) rt-aps x un2 x 103/sn2=450 x 2202 x 103/l 500002=0.968q rt*=rtxsb/vb2=0.968 x 100/2302=0.0018 xt-vx%xvn2x103/100/sn=14x2202x103/100/150000=45.173q xt*=xtxsbab2=45.1 73 x 100/2302=0.0854 bt=io%snx1o'3/1oo/vn2=o.8

12、x15oooox1o-3/1oo/23o2=o.o227x1o-3s bt*=btxvb2/sb=0.0227x10-3x2302/100=0.012b4(型号 sfp7150000/110)(参数 apo=lo7kw,aps=547kw?vs%=13%,io%=o.6,yndl 1) rt=547 x1102x103/1 500002=0.2942firt*=0.2942x 100/1152=0.0022xt=13x1102x103/100/1 50000= 10.487qxt*=10.487x 100/1152=0.0793bt=0.6x 150000x 1 0_3/100/1102=0.

13、744x 10_4sbp=0744x 10故 1152/l 00=0.00982、三绕组变压器:b5b8 (型号 sfpsz7-120000/220220±8xl.25%/115/10.5yn,yno,dll u2 =23.5%u23 =13.3%=7.7% apv = 425 vv = d0 )us% = -x (23.5 4-13.3-7.7) = 14.55%/ % =丄x (7.7 + 13.3 -23.5) = -1.25% « 0/$1 % =丄 x (23.5 + 7.7-13.3) = 8.95%rtq=apskwvn2kvx103/sn2k=425 x 2

14、302x103/1200002=1.5613=1.5613xl00/2302 =0.002951= ±55x 2202 x1q3 二5&6850100 120000xm = -x 22° xlo3 = 36.098q100 120000x厂= 58.685x100 = 0 11()9 =230201= 36.098x100 =0()68 =23020,11x|$ + x严=0.1109“242/1211k = 1230/1153、联络变压器 b9： r“ =艾孚=0.00567x厂=甌小。=0.2079 k = 231/110 = 0.95723022302230/

15、1154、线路l1l7:（l1,l5,l6为双回路,l2,l3,l4,l7为单回路）l1:型号为 2 x lgj -400/75km.参数为:r=0.08 q /km,x=0.397 q /km,b=2.92 x 10-6s/kmri=l/2lxr=75x0.08/2=3qxi=l/2lx=0.397x75/2=14.888qb=2xlb=2x2.92xl0-6x75=4.38xl0-4srp=rxsb/vb2=3x100/2302=0.00567xi*=x1xsb/vb2=14.888x100/2302=0.0281bp=b1xvb2/sb=4.38x10-4x2302/100=0.2317

16、 bp/2=0 .1158l2：型号为 lgj -400/50km.参数为:r=0.08 q /km,x=0.397 q /km,b=2.92 x 10-6s/kmr2=lxr=50x0.08=4qx2=lx=0.397x50=19.85qb2=lb=2.92xl0-6x50=1.46xl0-4sr2*=rxsb/vb2=4x 100/1152=0.0302x2*=x2>< sb/vb2=1 9.85 x 100/1152=0.1501b2*=b2xvb2/sb= 1.46x 10经 1152/l 00=0.0193b2*/2=0.00965l3:型号为 lgj-185/70km.

17、参数为:r=0.17 q/km,x=0.395 q/km,b=2.79x 10-6s/kmr3=lxr=70x0.17=11.9qx3=lx=0.395x70=27.65qb3=lb=2.79x 10_6x70=l .953x 10'4sr3*=rxsb/vb2=1 1.9x100/1152=0.09x3*=x3xsb/vb2=27.65x 100/1152=0.2091b3*=b3xvb2/sb=1.953x10-4x1152/100=0.0258b3，/2=0.0129l4：型号为 lgj -150/74km. 参数为:r=0.21 q/km,x=0.403 q/km，b=2.74

18、x 10-6s/kmr4=lxr=74x0.21=15.54。x4=lx=0.403x74=29.822qb4=lb=2.74x 10'6x74=2.03xl0_4sr4*=rxsb/vb2=l 5.54x 100/1152=0.1175x4*=x4xsb/vb2=29.822x100/1152=0.2255b4*=b4xvb2/sb=2.03 xloxll 52/l 00=0.0268b4*/2=0.0134l5:型号为 2 x lgj -120/70km.参数为:r=0.27 q /km,x=0.409 q /km,b=2.69 x 10-6s/kmr5=l/2lxr=70x0.2

19、7/2=9.45qx5=l/2lx=0.409x70/2=14.315qb5=2xlb=2x2.69xl0_6x70=3.766xl0_4sr5*=r5 x sb/vb2=9.45 x 100/1152=0.0714x5*=x5xsbab2=14.315x 100/1152=0.1082b5*=b5xvb2/sb=3.766x10_4x1152/100=0.0498b5*/2=0.0249l6：型号为 2xlgj -400/73km.参数为:r=0.08 q/km,x=0.397 q/km,b=2.92x 106s/kmr6=l/2lxr=73x0.08/2=2.92qx6= 1 /2lx=0

20、.397x 73/2= 14.491qb6=2xlb=2x2.92xl0-6x73=4.26xl0_4sr6*=r6xsb/vb2=2.92x1 00/2302=0.0055x6*=x6xsb/vb2= 14.491 x 100/2302=0.0274b6*=b6xvb2/sb=4.26x 10'4x2302/l 00=0.2254b6*/2=0.1127l7:型号为 lgj -400/135km.参数为:r=0.08 q/km,x=0.397 q/km,b=2.92x 106s/kmr7=lxr=135x0.08=10.8qx?=lx=0.397x 135=53.595qb7=lb=

21、2.92x 10'6x 135=7.844x10_4sr7*=r7x sb/vb2= 10.8x1 00/2302=0.0204x7*=x7xsb/vb2=53.595x100/2302=0.1013b7*=b7xvb2/sb=7.844x10-4x2302/100=0.4149b7*/2=0.207455、发电机 f1-f4：(型号：sf-100/40-854 pe=100mw cos = 0.9 ue =13.8kvx 八=0.20346. pq节点pq标么值:4：s4=180+jl00s4*=s4/sb=1.8+jl.o5：s5=72+j40s5*=s5/sb=0.72+j0.4

22、8：s8=120+j63s8*=sg/sb=l 2+j0.639：s9=36+j20s9 尸s9/sb=0.36+j0.210： si0=40+j25sio*=sio/sb=o.4+jo.257零序电抗标么值:x0(l1) =4.5xll =4.5x0.0281 = 02645 x0(l2) = 3xs =3x0.1501 = 0.4503 x(u3)= 3x3 =3x 0.2091 = 0.6273 x0(/4) =3x l4 = 3 x 0.2255 = 0.6765 x0(l5) =3xl, =3x0.1082 = 0.3246 xg6)=4.5x“ =4.5x0.0274 = 0233

23、 x0(l7) =4.5xl7 =4.5x0.1013 = 0.45585电力系统潮流计算p-q分解法潮流计算主要步骤1导纳矩阵的形成2因子表的形成（三角分解法）3给定电压初值4.计算不平衡功率 r（k）除以"仮）5计算不平衡功率a q/k）除以vi（k）6判断收敛性7回代，修正屮）,a vj（k）8计算全线路功率电力系统的潮流计算机算法潮流计算的任务：对给定的运行条件确定运行状态，如果各母线上的电压（幅值及相角）网络 小的功率分布以及功率损耗等，几个节点电力系统的潮流方程的一般形式根据电力系统的实际运行条件，按给定变量的不同，一般节点可分三种：1. pq节点：有功功率p和无功功率q

24、给定的，节点电压（v, 8 ）是待求量。2. pv节点：有功功率p和电压幅值v给定的，节点的无功功率q电压的相位角§是待求量。3. 平衡节点：网络中至少有一个节点的有功功率p不能给定，此节点承担系统的自功功率平 衡。网络方程是线性方程，但由于在定解条件小不能给定节点电流，只能给出节点功率，从而使潮 流方程变为非线性方程，由于平衡节点的电压已经给定，假设系统屮有n个节点，其屮有m个pq 节点，n-（m+l）个pv节点和一个平衡节点，平衡节点不参加求解。形成导纳矩阵的方法及框图1. 形成导纳矩阵的程序框图为了形成导纳矩阵，必须个计算机输入必要的原始数据，一条支路的原始数据应包括两端的节点

25、号和支路阻 抗，对变压器支路的原始数据应包括他的变比，对于电力线路还应包括它的容纳，即：i , j , r , x , k（-b/2）o对以上数据说明如下：（1）当支路为接地支路时，规定节点i处置零，j处填写接地支路所在的节点号，支路参数用阻抗填写, 对于k处置零，但必须填写。（2）对于纯阻抗支路，k处置零。（3）对于具有容纳的线路，电容电纳不作单独支路处理，而且把总容纳的一半负值填写在k处，其符 号用以区别支路的性质。（4）对于变压器支路，采用下图所示等值电路，即非标准变比在j侧，变压器阻抗在1侧，k处填写 实际变化yb/kzbl:k2.框图注:bib2不记充电电容和非标准变比时，导纳矩阵的

26、虚部 导纳矩阵的虚部二、形成因子表的方法及框图1.用行消去过程形成因子表的程序框图2. 因子表的形成（三角分解法）消去法求非线性方程组的一种常川算法是对方程式的系数矩阵a进行三角分解，在本次电力系统潮流计算中 采用的三角分解是将非奇方阵a分解为单位下三角矩阵l和上三角矩阵r的乘积。a=lr,非奇方阵a被表示为矩阵l和r的乘积：这两个三角矩阵称为a的因子矩阵，两个因子矩阵的元 素计算公式：z = l=s e ggk = 勺=1 ,2 厂nj = i , i + 1 , 丿将a=lr代入线性方程组, lf=b便得lrx=b,这个方程又可以分解为以下两个方程: rx=f展开为厶21厶32厶42ni1

27、21323nn先由方程组自上而下地依次算出f、f2>行，其计算通式为x - y lufjj = i这一步演算相当于消元过程中对原方程式右端常数向量所作的变换，只顺用到下三角因子矩阵。方程组的 求解属于回代过程，只顺用到上三角因子矩阵以及经过消元变换的右端常数向量，方程组可以自下而上地逐步 算出待求量，其计算通式为：nv1 (，一1) n ? n -1，1 )八-1)rijl ru xj j = i+三、节点电压的表示：1. 极坐标表示法：匕=%（sin$）节点功率方程表示为：rni=<j qi =匕 s 匕 5 sin 叭- bij cos 叭）；=ig = c-j （两节点电压的

28、相位角）方程式把节点功率表示为节点电压的幅值和相角的函数在有n个节点的系统屮，假定第1m号节点为pq节点，第m+1nl号节点为平衡节点。v.和片是给定 的。pv节点的电压幅值vm+1vn也是给定的。因此,只剩下n-1个节点的电压相角§1. 6和m个节点的电 压幅值w.vm是未知的量。对于每一个pq节点或每一个pv节点都可以列写一个有功功率不平衡方程：遊厂klxgcos爲+dsi叫)=0 (i = l,2,-n-1)对于每冃一个pq节点可以列写一个无功功率不平衡方程式：aq 二qis-qi 二qis -vy.(gjj si叫-b. cosq) = 0 (i = 12m)7=1所以可以写

29、出修正方程式_ap_hna<2klvd'2 av厂ap =w np.3 */=av =vp2 =匕_ a匕_-冬四、pq分解法1.在交流高压电网中，输电线路的电抗耍比电阻大得多，系统中母线有功功率的变化主耍受电压相位的影 响，无功功率的变化规则主耍受母线电压幅值变化的影响，在修止方程式的系数矩阵中，偏导数乞竺和竺 dv 58 的数值相当于偏导数乞竺和是相当小的，作为简化第一步，可将方程中n、k略去不计，既认为它 ds 5v们的元素都等于零，这样，便可分解为ml阶和m阶的两个方程：- av所谓p-q分解法，节点有功功率不平衡量又用于修正电压的相位。节点的无功功率不平衡最值用于修正电

30、 压幅值，方程、分别轮流迭代。2. h、l元素都是节点电压幅值和相角差函数，其数值在迭代过程中是不断变化的，从而，最关键的一 步是把系数矩阵h、l简化为常数矩阵。方法：在一般情况下，在线路两端电压的相角差是不人的(不超过0。20。)，因此，以为 cos站=1 ,此外，与系统各节点无功功率相适应的导纳bg必须小于该节点白导纳虚部，即bld严寻i或0"利矩阵h、l可简化为h = vp1 - bwdi ®l = vd1b"%2将分别带入得:= -vdbvdi3a0 = -vd2bav= -bf a5vdq =即：v2b、b2b3b*21b22b屮-1坊1,1eb r%i

31、/l1v262y r41攻 2 43 by”a%片a2221222mv2直泌avm利用计算节点功率的不平衡量，用修疋方程解出修疋量a §及av,并换下述条件: max ”*)< s p maxaq$)< cq校验收敛这就完成分解法的计算了.3p-q分解法计算潮流的程序框图pq分解法源程序include hrnath.hn#include "stdio.h"#include mstdlib.hm main ()int a=l,k,i j,d=l,m,n,ls,nl,kk;float z 1005,pq1003,p,q,r,x,kb,b,g;float y

32、g 100 100,yb 100 100,b 1 100 100,b2 100 100,v 100,w 100;float pp 100,qq 100,e 100,fl 00,gg 1002,tt 1002,pw 100,qw 100,ww 100;float t,tl ,tp,tq,p 1 ,p2,q 1 ,q2,t2;printfcxn");printf(u电力系统潮流计算n”)；printfc'nn);prinlff指导老师：设计：2002年12月5“)；printfn”)；printtr'vn");printf("注：各参数用标么值表示n&

33、quot;);prin(fcn”)；begin:printfce);prinlf(“请选择:1健立电网;2 査看数据;3 修改数据;4 运行计算;5 退出程序:”);scanff%d",&a);if (a=l)goto creat;else if (a=2)goto view;else if (a3)goto modify; else if (a=4)goto run;elseexit(0);creat:printf("请输入系统的节点数 n:");scanf("%d”,&n);prints请输入系统的支路数ls:”);scanfc%d”

34、,&ls);prin垠”请输入系统的pq节点数m:“);scanfc%d”,&m); fbr(k=l;k<=ls;k+)printf(“请输入支路的(1的数拥1(接地支路为0) :“,k);scantr%d”,&i); print”请输入支路的d 的数据 j:",k);scanf(”d”,&j);printft请输入支路的d 的数据 r:”,k);scanf(“f',&r); printf(“请输入支路的d 的数据 x:”,k);scanf(“f',&x);print"请输入支路的d 的数据 kb:”,k

35、);scanfc%f；&kb); zki=i;zk2=j;zk3=r;zk4=x;zk5=kb;for(k-l ;k<m;k+)pqkl=k;print”请输入pq节点d的注入有功p:”,k);scanfc%f；&p);print”请输入pq节点d的注入无功q:“,k);scanf(“玖&q); pqk=p；pqk3=q;fbr(k=m+1 ;k<=n-l ;k-h-)pqkl=k;printft请输入 pv 节点d 的注入有功 p:”,k);scanff%f；&p); printf(”请输入 pv 节点d 的电压模:”,k);scanfc%f；&a

36、mp;q); pqk2=p;pqkh3=q；pqkl=n；printnnin输入平衡节点d的电压模:”,n);scanfc%f；&p);printf(”请输入平衡节点d的电压角:”,n);scanfr%u&q); pqn2=p;pqn3=q;goto begin; vicw:printf(“n”);printfcn=%dtls=%dtm=%d”,n,ls,m);printfcn”)； printfc 支路数据 “);printffn”)；printff't节点it节点j电阻t电抗t变比或电纳”);printffn”);for(i=l;i<=ls;i+)prints

37、支路叫);for(h;j<=5;j+)if(jv=2)printf(” t %1.0r,zi|j);elseprint"%14.6 化 zi|jd；pringw);注入有功tprintfvtpq 节点t for(i=1 ;i<=m;i+)f0r(j=lij<=3;j+)if(j<=l)printf(ut%1.0f',pqi|j)；elseprintf(ht%14.6r,pqiu)；printf(恤”);printf(ntpv 节点t注入有功t电圧模“);printf("n”);fbr(i=m+1 ;i<=n-l ;i+)for(j=l;

38、j<=3;j 卄)if(j<=l)printfct%1.0fpqiu)；elseprintf(%14.6f"qij)；print”*”)；printnnnm);printf(”t平衡节点t电压模t电压角");printf(mn");for(j=l;j<=3;j 卄)if(j<=dprints %1.0f,pqnj);elseprintf("t%14.6r,pqn|j);printf(nnm); goto begin;modify: prints 1-修改支路数据;2修改节点注入功率;3返回主菜单t);scanfc%出,&d

39、); if (d=l)goto ml;else if (d=2)goto m2;else if (d=3)goto begin;m 1: printft请输入要修改的支路号:”);scan" %d”,&k); prints请输入支路的d的新数据i:”,k);scanf(”d”,&i); printf(uis输入支路的d 的新数据 j:”,k);scanf(“d“,&j); printf(”请输入支路的d 的新数据 r:”,k);scanft%f',&r); printf("iff输入支路的d 的新数据 x:",k);scan

40、f("%f',&x); printft请输入支路的d 的新数据 kb:”,k);scanf(”f',&kb); zkl=i;zk2=j;zk3=r;zk4=x;zk5=kb;goto modify;m2:k=0;print” 请输入需修改的节点号:”);scanf(”d“,&k);pqk1=k;prints请重新输入节点的注入有功或电压参数:”);scanf(” f；&pqk2);printft请重新输入节点的注入无功或电压参数:“);scanf(”f,&pqk 3); goto modify;run : fbr(i=l;i&l

41、t;=n;i+)for(j=l;j<=n;j-h-)ybiu=o；形成bl,不计线路电容和变压器非标准变比*/for(k=l;k<=ls;k-i)i=zk 1 y =zk 2 ;r=zk 3 ;x=zk 4;b=-x/(r*r+x*x);if(i=0)goto po;ybii=ybii+b;ybiu=ybiu-b;ybui=ybiu；po：ybuju=ybuu+b;fbr(i=l;i<n;i-h-)fbr(j=l;j<=n;j+)bliu=ybiuj;/*形成导纳矩阵/for (i=l;i<=n;i-h-)for (j=l;j<=n;j+)ygiu=o； y

42、biu=o；for (k=l ;k<=ls;k-h-)i=zkll;j=zkj2;r=zkj ;x=zk4;kb=zk；b=-x/(r*r+x*x);g=r/(r*r+x*x);if(i=0)goto p 1 ；if(kb>0)goto p2;ybii=ybii-kb;ybjj-ybjjkb；p=l;q=l;goto p3；pl: p=l;goto j-)4;p2:p=kb*kb;q=kb;p3：ygii=ygii+g;ybii=ybii+b;ygiu=ygiu-g/q;ybiu=ybiu-b/q；yglii=ygiu;ybji=ybiu;p4：ygjlj=ygjj+g/p;ybj

43、j=ybjj+b/p；/*形成b2*/fbr(i=l;i<=m;i+)fbr(j=1;j<=m;j+)b2ij=ybij;printf (”矩阵 b:nn);for(i=l;i<=n;i+)for(j=l；j<=n；j+)printfc%14.7fybij)； printfc'n");printf (”矩阵 g:n”);fbr(i=l;i<=n;i+)for(j-l;j<=nu-h-)printfc%14.7f,ygij);printfc'");/*形成bl因子表*/fbr(i=2;i<=n-l ;i+)fbr(k=

44、l;k<=il;k+)blik=-blik/blkk;for (j=k+l;j<=n-l;j+)biili=bliu+blik*blku；printfv'b 1 因子表:n”)；fbr(i=l;i<=n-l;i+)for(j=l;j<=n-ly-h-)printf("%14.7f,blij);printf("nh);/*形成b2因子表:*/for(i=2;i<=m;i-h-)for(k-1 ;kv=i-1 ;k+)b2ik=.b2ik/b2kk;for (j=k+l;j<=m;j+)b2iuj=b2iu+b2ik*b2ku;pri

45、nttr'b2 因子表:n”)；fbr(i=l ;i<=m;i-h-)for(j=l;j<=m;j+)print” 14.7 化b2ij);printf("nh);print”®”)；/*给定电压初值*/for(i=l;i<=m;i+)vi=l；wi=o;for(i=m+1 ;i<=n-lvi=pqi3；wi=o;vn-pqn2;wn=pqn3;kk=o;严计算n-1个节点不平衡有功并除英电压引al:for(i=1;i<=n-l;i+)t=0;fdr(j=l;j<=n;j+)if(ybiu!=o)t=t+vj*(ygij*cos(

46、wi 卜 wj)+ybij*sin(wi 卜 wj);严计算m个节点不平衡无功并除其电压卄fbr(i=l;i<=in;i+)t=0;for(j=l；j<=n;j+)if(ybi|j!=o) t=t+vj*(ygi|j*sin(wi-wj)-ybij*cos(wi-w|j);qqi=pqi3/vi-t;for(i=l ;i<=n-l ;i+)if(fabs(ppi)>0.00001)goto a4 ;fbr(i=l;i<=m;i-h-)if(fabs(qqij)>0.00001)goto a2;goto a3;/*利川bl因子表对右端项变换*7a4:for(i

47、=l ;i<=n-2;i+)t=0;for(j=l;j<=i;j+)t=t+bii+iu*ppu；ppi+l=ppi+l+t;nl=n-l;/*回代求解*/ppnl=ppnl/b1nlnl;fbr(k=l;k<=nl-l;k+)t=0;for(j=nl-k+1 ;j<=nl;j+)t=t+blnl-kj*ppu;ppnl-k=(ppnl-k-t)/b 1 nl-k nl-k; 修正电压角*/fbr(i=l;i<=nl;i+)wi=wi-ppi;/*利用b2因子表对右端项变换*/a2:fbr(i=1;i<=m-l;i+)t=0;fbr(j=l;j<=i;

48、j-h-)t=t+b2i+ld*qq|j;qqi+l 戶 qqi+l+t;qqm=qqm/b2m m;/*回代求解*/ fdr(k-l;k<-m-l;k+)t=0；fbr(j=n-k+1 ;j<=m;j+)t=t+qqj*b2m-kj;qqm-k=(qqm-k-t)/b2m-k m-k;/*修正电压模*/fbr(i=l;i<-m;i+)vi=vi-qqi; kk=kk+; goto al;a3:for(i=l ;i<=n;i-h-)ei=vi *cos(wi);fli=vi*sin(wi);wwi 卜 wi*180/3.1415926;计算支路功率for(k= 1 ;k

49、<=ls;k+)i=zkl;j=zk2;kb-zk5; ti=-ygij*(ei-eu)+ybilj*(tti-ftj);t2=-ybiu*(ei.eu)-ygid*(fi-fj);pl=ei*tl+f!i*t2;/*ij 节点间有功的部分 vi(uiuj)yij*/ ql=fi*tl-ei*t2;p2=ej*(-tl)+fd*(-t2); q2=fj*(-tl)-ej*(-t2);if(kbx)ggkl =p 1 +(kb 1 )*vi*vi (-ygij)； ttkl =q i -(kb-1 )*vi*vi*(-ybilj)； ggk2=p2+(l-kb)*vlj*vj*(-ygij

50、)/kb;ttk2=q2-(l-kb)*vj*vu*(-ygiuykb;elseggk l=pl ;ggk 2=p2;/* 有功 */ ttkl=ql+kb*vi*vi；/* 无功 */ ttk2=q2+kb*vj*vj;/*计算支路功率损耗*/for(k=l;k<=ls;k-h-)pwk=ggkl+ggk2; qwk=ttkl+ttk2; tp=o;tq=o;lor(k=l ;k<=ls;k+)tp=tp+pwk;tq=tq+qwk;电压角（度）nh,kk）;printf(“迭代(1次电压实部e电压虚部f电压模vfor(i=l;i<=n;i+)printff 节 %6d%14.7f%14.7f%14.7f%14.7fn,i,ei,fi,vi,wwi)； printfn”)；print”支路功率：有功无功n");fbr(k-l;k<=ls;k-h-)i=zkl;j=zk2;printf (”支路 %d-%d%16.7f% 18.7fn”,i j,ggk