电力系统潮流及短路电流计算程序

班级：姓名：学号：一、作业要求编写程序计算图 1 所示算例系统的潮流及三相短路电流。潮流计算：方法不限，计算系统的节点电压和相角。短路电流：4 号母线发生金属性三相短路时(zf=0）,分别根据精确算法和近似算法计算短路电流、系统中各节点电压以及网络中各支路的电流分布，并对两种情况下的计算结果进行比较.二、电路图及参数图 1 3 机 9 节点系统表 1 9 节点系统支路参数支路R（p.u。）X（p。u。)B/2（TK） （p.u。）1~400.05761.02~700.06251。03～900.05861.04～50.010.0850。0884～60。0170.0920.0795～70.0320.1610.1536~90。0390。170.1797~80.00850.0720.07458~90。01190。10080。1045表 2 9 节点系统发电机参数发电机编号节点类型PG（p.u。)VG（p。u。)（p。u.）（p.u.)1平衡1.040.31.1372PV1.631.0250.31。2113PV0。851。0250。31。043表 3 9 节点系统负荷参数节点编号节点类型Pi(p。u。）Qi(p.u.）4PQ005PQ1。250。56PQ0.90.37PQ008PQ10.359PQ00三、计算步骤(1)进行系统正常运行状态的潮流计算，求得(2)形成不含发电机和负荷的节点导纳矩阵 YN；(3)将发电机表示为电流源（)和导纳(）的并联组合;节点负荷用恒阻抗的接地支路表示;形成包括所有发电机支路和负荷支路的节点导纳矩阵 Y,即在 YN中的发电机节点和负荷节点的自导纳上分别增加发电机导纳和负荷导纳（）；(4)利用，计算节点阻抗矩阵,从而得到阻抗矩阵中的第 f 列；(5)利用公式(6-7）或(6—10)计算短路电流；(6)利用公式（6-8)或(6—11）计算系统中各节点电压；(7)利用公式（6—9）计算变压器支路的电流;对输电线路利用 П 型等值电路计算支路电流.四、计算结果节点导纳矩阵 Yn： Columns 1 through 5 0 —17.3611i 0 0 0 +17.3611i 0 0 0 -16。0000i 0 0 0 0 0 0 —17.0648i 0 0 0 +17。3611i 0 0 3。3074 -39。3089i —1.3652 +11.6041i 0 0 0 -1。3652 +11.6041i 2.5528 -17。3382i 0 0 0 -1。9422 +10.5107i 0 0 0 +16。0000i 0 0 -1.1876 + 5.9751i 0 0 0 0 0 0 0 0 +17.0648i 0 0 Columns 6 through 9 0 0 0 0 0 0 +16.0000i 0 0 0 0 0 0 +17.0648i -1.9422 +10.5107i 0 0 0 0 -1.1876 + 5。9751i 0 0 3.2242 —15.8409i 0 0 —1.2820 + 5。5882i 0 2.8047 —35.4456i -1.6171 +13。6980i 0 0 -1.6171 +13.6980i 2。7722 —23...