南昌大学电力系统分析实验报告

1、南昌大学实验报告实验类型：口验证综合口设计口创新实验日期：实验成绩：一、实验项目名称电力网数学模型模拟实验二、实验目的与要求：本实验通过对电力网数学模型形成的计算机程序的编制与调试，获得形成电力网数学模型：节点导纳矩阵的计算机程序，使数学模型能够由计算机自行形成，即根据已知的电力网的接线图及各支路参数由计算程序运行形成该电力网的节点导纳矩阵。通过实验教学加深学生对电力网数学模型概念的理解，学会运用数学知识建立电力系统的数学模型，掌握数学模型的形成过程及其特点，熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来

2、。三、主要仪器设备及耗材计算机、软件（已安装，包括各类编程软件C语言、C+、VB、VC等、应用软件MATLAB等）、移动存储设备（学生自备，软盘、U盘等）四、实验步骤1、将事先编制好的形成电力网数学模型的计算程序原代码由自备移动存储设备导入计算机。（1）编程思想I无变压器支路的节点导纳矩阵计算方法。以下语句用于输入题目已知节点之间的导纳和阻抗值：z=input（请输入由节点号对应的阻抗形成的矩阵：z=）;y=input（请输入由节点号对应导纳形成的矩阵：y=）;%其中,.：即为i节点对地导纳，两节点之间无直接相连的通路则/输入为0，输入为inf（即无穷大），也输入为inf。以下语句用于计算无变

3、压器支路的节点导纳矩阵：for(i=1:n)for(j=1:n)Y(i,i)=sum(y(i,:),2)+sum(1./z(i,:),2);ifj=iY(i,j)=Y(i,i);elseY(i,j)=-1/z(i,j);endendend%其中，对角线元素):=.:.二+二;:=.:.二+工_:非对角线元素Yij=y(k1:kkzkz/(k-l)k*kz/(l-k)1:kkzkz/(k-l)k*kz/(l-k)当节点a，b间接有变压器支路时（见图i）,当然可以用n型等值电路，然后按照上述原则形成导纳矩阵。但在实际应用程序中，往往直接计算变压器支路对导纳矩阵的影响。根据图1-5-1可以写出节点a

4、，b的自导纳和节点间的互导纳增量分别如下：节点a的自导纳改变量式（1-1）：图1图1变压器支路的n型等值电路I-1)节点j的自导纳改变量式I-1)节点j的自导纳改变量式(1-2)：I-2)增加节点i,j间的互导纳式(1-3)：I-3)注意：在输入电路已知参数时，将k以矩阵形式输入，若i、j两节点之间无变压器，贝応.=1,G-：=1),若i、j之间有变压器则将靠近变压器的一端记为节点j，靠近阻抗的一端记为节点i,输入人员将输入为变压器变比，而厂以1输入。程序段中for(i=1:n)for(j=1:n)ifk(i,j)=1y(j,i)=(1-k(i,j)/(k(i,j)*k(i,j)*z(i,j)

5、;y(i,j)=(k(i,j)-1)/(k(i,j)*z(i,j);z(i,j)=k(i,j)*z(i,j);z(j,i)=z(i,j);endendend%判断两个节点之间是否有变压器，如果k二】则对节点之间的导纳和阻抗矩阵进行修改，得到新的汀和二源代码如下：clcclearn=input(请输入节点数：n=);z=input(请输入由节点号对应的阻抗形成的矩阵：z=);y=input(请输入由节点号对应导纳形成的矩阵：y=);k=input(请输入变比：k=);Y=zeros(n);for(i=1:n)for(j=1:n)ifk(i,j)=1y(j,i)=(1-k(i,j)/(k(i,j)

6、*k(i,j)*z(i,j);y(i,j)=(k(i,j)-1)/(k(i,j)*z(i,j);z(i,j)=k(i,j)*z(i,j);z(j,i)=z(i,j);endendendfor(i=1:n)for(j=1:n)Y(i,i)=sum(y(i,:),2)+sum(1./z(i,:),2);ifj=iY(i,j)=Y(i,i);elseY(i,j)=-1/z(i,j);endendend2、在相应的编程环境下对程序进行组织调试。3、应用计算例题验证程序的计算效果。将例题中的已知条件如下依次输入:n=5z=inf0.105iinfinfinf;0.105iinf0.024+0.065i0

7、.03+0.08iinf;inf0.024+0.065iinf0.018+0.05iinf;inf0.03+0.08i0.018+0.05iinf0.184i;infinfinf0.184iinfy=00000;000+0.0160i0+0.0200i0;00+0.0160i00+0.0130i0;00+0.0200i0+0.0130i00;00000k=1.00001.05001.00001.00001.0000;1.00001.00001.00001.00001.0000;1.00001.00001.00001.00001.0000;1.00001.00001.00001.00001.00

8、00;1.00001.00001.00000.96001.00004、对调试正确的计算程序进行存储、打印。5、完成本次实验的实验报告。五、实验数据及处理结果运行自行设计的程序，把结果与手工计算结果相比较，验证所采用方法及所编制程序运行的正确性。1手算节点导纳矩阵：详细计算过程见课本P73例4-1其结果为：0.0000-j9.52380.0000+j9.07030000.0000+j9.07039.1085-i33.1002j13.5388-4.9989j10.9589-4.106900j13.5388-4.998911.3728-j31.2151j17.7053-6.373900j10.958

9、9-4.1069j17.7053-6.373910.4835-j34.52830.0000+j5.66120000.0000+j5.66120.0000-j5.43482用matlab程序计算节点导纳矩阵实验结果：1234510.0000-9.5238i0.0000+9.0703i0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi20.0000+9.0703i9.1.085-33.1001i-4.99&04-13.538Si-4.10964-10.9589i0.00004-O.OOOOi30.0000+O.OOOOi49990+1.3.53S8i11.3729

10、-31.2153-6.3739+17.7054i0.0000+O.OOOOi40.0000+O.OOOOi-4.1096+1.0.9589i-6.3739+17.7054i104835-34.5284i0.0000+5.6612i50.0000+O.OOOOi0.00004-O.OOOOi0.00004-O.OOOOi0.00004-5.6612i0.0000-5.434Bi对比可知通过比较，发现得到的节点导纳矩阵中大部分的数据是相同的，只是少部分的结果有微小误差。这也就说明了手算能得到同样的结果，但我们发现用一般手算计算节点导纳矩阵比用计算机编程计算节点导纳矩阵繁琐了很多，而一个系统的节点数

11、不可能是几个，而是由许多的节点构成的，若在此时用一般手算计算节点导纳矩阵，这就显得非常复杂，在实际操作过程当中也没有那么多的时间去让我们去用手算计算节点导纳矩阵。所以在电力系统中我们往往都是采用计算机编程计算节点导纳矩阵，这样大大提高了工作效率。六、思考讨论题或体会或对改进实验的建议1什么是输入阻抗？什么是转移阻抗？网络化简的方法有哪些？答：输入阻抗：在节点i单独注入电流而所有其他节点的注入电流都等于0时，在节点i产生的电压同注入电流之比，即等于节点i的自阻抗，即输入阻抗。转移阻抗：节点i单独施加Ei时，该电势与其在节点j产生的短路电流Ij之比即等于节点i、j之间的转移阻抗。网络化简方法：星网

12、变换法，高斯消去法。简述节点导纳矩阵的形成的过程，节点导纳矩阵的阶数与系统的节点数有什么关系？节点导纳矩阵的互导纳yij在数值上等于什么？举例说明当网络结构发生变化时，如何修改？答：节点导纳矩阵阶数等于系统的节点数；节点导纳矩阵的互导纳y.在数值上等于联接节点i、j的支路导纳的负值；ij节点导纳矩阵的修改（1）从网络的原有节点i引出一条导纳为.:亡的支路，同时增加一个节点k，由于节点数加1,导纳矩阵将增加一行一列。新增的对角线元素込尸二上。新增的非对角线元素中，只有匚-：丄，其余的元素都为矩阵的原有部分，只有节点i的自导纳应增加二二=在网络的原有节点i、j之间增加一条导纳为汀的支路。由于只增加

13、支路，不增加节点，故导纳矩阵的阶次不变。因而只要对与节点i、j有关的元素分别增添以下的修改增量即可：二:=二匚=：.,，二.=二：=-：其余元素都不用修改。在网络的原有节点i、j之间切除一条导纳为汀的支路，这种情况可以当做是在i、j节点间增加一条导纳为-.：,的支路来处理，因此,导纳矩阵中有关元素的修正增量为二:=二=一，二=二；=举例：如将例4-1中节点4、5间变压器变比由k=0.96改为k=0.98后作节点5454导纳矩阵的修改：变比修改后，相当于先切除变比为k的变压器，再接入变比为k的变压器。与5454之相关的4、5节点的导纳值要发生变化，而网络的其他元素不发生变化。利用例4-1中结果，

14、与节点4、5有关的导纳矩阵元素的修正增量应为AY55=0,AY4410.982xj0.18410.962xj0.184=j0.2382AY=AY544510.98xj0.18410.96xj0.184=-j0.1155则变比变化后Y=10.4835-j34.5283+j0.2382=10.4835-j34.290144=Y=j5.6612-j0.1155=j5.5457=-j5.4348节点电压方程的求解方法有哪些？答：直接求解法和高斯消去法简述节点阻抗矩阵元素的物理意义及其形成。答：节点导纳矩阵中=:=.当k=i时，上述公式说明，当网络中除节点i以外，所有节点都接地时，从节点i注入网络的电流

15、同施加于节点i的电压之比，即等于节点i的自导纳二。换句话说，自导纳二是节点i以外所有节点都接地时节点i对地的总导纳。当k=i时，公式说明，当网络中除节点k以外所有接地都接地时，从节点i流入网络的电流同施加于节点k的电压之比，即等于节点k、i之间的互导纳匚.。在这种情况下，节点i的电流实际上是自网络流出并流入地中的电流，所以应等于节点k、i之间的支路导纳的负值。七、参考资料电力系统分析何仰赞华中科技大学出版社电力系统稳态分析陈珩中国电力出版社电力系统暂态分析李光琦中国电力出版社电力系统计算水利电力出版社一、实验项目名称电力网数学模型模拟实验二、实验目的与要求：本实验通过对电力网数学模型形成的计算

16、机程序的编制与调试，获得形成电力网数学模型：节点导纳矩阵的计算机程序，使数学模型能够由计算机自行形成，即根据已知的电力网的接线图及各支路参数由计算程序运行形成该电力网的节点导纳矩阵。通过实验教学加深学生对电力网数学模型概念的理解，学会运用数学知识建立电力系统的数学模型，掌握数学模型的形成过程及其特点，熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。三、主要仪器设备及耗材计算机、软件（已安装，包括各类编程软件C语言、C+、VB、VC等、应用软件MATLAB等）、移动存储设备（学生自备，软盘、U盘等）四、实验

17、步骤6、将事先编制好的形成电力网数学模型的计算程序原代码由自备移动存储设备导入计算机。7、在相应的编程环境下对程序进行组织调试。8、应用计算例题验证程序的计算效果。9、对调试正确的计算程序进行存储、打印。10、完成本次实验的实验报告。五、实验数据及处理结果运行自行设计的程序，把结果与手工计算结果相比较，验证所采用方法及所编制程序运行的正确性。实验程序：clearclc;z=inf,0.10+0.40i,0.3i,0.12+0.50i;0.10+0.40i,inf,inf,0.08+0.40j;0.3i,inf,inf,inf;0.12+0.50i,0.08+0.40i,inf,inf;（各支路

18、的阻抗）y=0,0.01528i,0,0.01920j;0.01528i,0,0,0.01413j;各支路的导纳）0,0,0,0;各支路的导纳）0.01920i,0.01413i,0,0;Y=zeros(4,4);for(i=1:4),for(j=1:4),ifi=jY(i,j)=Y(i,j)elseY(i,j)=-1.0/z(i,j)endendendfor(i=1:4),for(j=1:4),Y(i,i)=Y(i,i)+y(i,j)+1.0/z(i,j)endendY(1,1)=Y(1,1)-1.0/z(3,1)+1.1*1.1/z(3,1)(对变压器支路的两个节点进行修正)Y(1,3)=

19、Y(1,3)*1.1Y(3,1)=Y(3,1)*1.1运行结果，实验结果导纳截图：1.0421-8.2429i-0.5882+2.3529io+:3.666Ti-0.4539+1.8911i-0.5882+2.3529i1.0690-4.T2T4i0-0.4808+2.4038i0+3.6667i0o-:3.3333i0-0.4539+1.8911i-0.4808+2.4038i00.9346-4.2616irestart手工计算结果：Yyyzzz11=140+120+1/12+1/14+1.1*1.1/13=j0.01920+j0.01528+1/(0.10+j0.40)+1/(0.12+j

20、0.50)+1.1*1.1/j0.3=1.0241-j8.2429YY12=21=-1/z12=1/(0.10+j0.40)=-0.5882+j2.3529YYz13=31=-1.1/13=-1.1/j0.3=j3.667YYz14=41=-1/14=1/(0.12+j0.50)=-0.4539+j1.8911Yyyzz22=210+240+1/21+1/24=j0.01528+j0.01413+1/(0.10+j0.40)+1/(0.08+j0.40)=1.069-j4.7274YYz24=42=-1/24=-1/(0.08+j0.40)=-0.4808+j2.4038Yz33=1/13=1

21、/(j0.3)=-j3.333Yyyzz44=410+420+1/14+1/42=j0.01920+j0.01413+1/(0.12+j0.50)+1/(0.08+j0.40)=0.9346-j4.2616六、思考讨论题或体会或对改进实验的建议1什么是输入阻抗？什么是转移阻抗？网络化简的方法有哪些？答：输入阻抗就是指节点的自阻抗Z即节点i单独注入电流I时，在节点i产生iii电压V与电流I之比；转移阻抗Z就是指节点i单独注入电势，与在短路点jiiji产生的电流之比；网络化简方法有高斯消去法：也就是带有节点电流移置的星网变换；简述节点导纳矩阵的形成的过程，节点导纳矩阵的阶数与系统的节点数有什么关系

22、？节点导纳矩阵的互导纳yij在数值上等于什么？举例说明当网络结构发生变化时，如何修改？答：(1)节点导纳矩阵的形成过程：列写节点电压方程：以零电位点为计算节点电压参考点，根据基尔霍夫电流定律，写出各独立节点电流平衡方程，经整理后形成矩阵方程：YV二I。矩阵Y就是所求节点导纳矩阵。(2)节点导纳矩阵阶数等于系统的节点数(注：不包括接地点)。节点导纳矩阵的互导纳Y数值上等于节点i、间支路导纳的负值，即-y。ijij对于网络结构发生的变化主要有这几种形式：从网络的原有节点i引出一条支路；在网络的原有节点ij之间增加一条支路；在网络的原有节点之间切除一条支路。下面如将例中节点1、间变压器变比由K13=

23、l.l改为K=1后作节点导纳矩阵的修改:变比修改后，相当于先切除变比为1.1的变压器，再接入变比为1的变压器。与之相关的1、3节点的导纳值Y33、丫弘Yi3要发生变化，而网络的其他元素不发生变化。利用上述的计算结果11=Yz11=Y11-1.1*1.1/+1/13=1.0241-j8.2429-1.1*1.1/j0.3+1/j0.3=1.0241-j7.5429Y33=33=1/(j0.3)=-j3.333Yz13=31=-1/13=-j3.333节点电压方程的求解方法有哪些？答;节点电压方程的求解方法有：直接求解法和高斯消去法简述节点阻抗矩阵元素的物理意义及其形成。答：节点阻抗矩阵元素的物理

24、意义：在节点k单独注入电流，在节点k产生的电压同注入电流之比，即等于节点k的自阻抗Z;在节点i产生的电压同节点k注kk入电流之比，即等于节点k和节点i之间的互阻抗Z。ik节点阻抗矩阵的形成：将节点电压方程：YV二I改写成ZI二V,式中Z二Y-1就是阻抗矩阵七、实验体会通过本次实验，对节点导纳矩阵的求解有了更深的认识，特别是怎么求解一个网络的支路节点导纳和有变压器的支路的导纳有了更深的了解。为以后在求解类似的问题中有了一定的基础。通过实验呢更是加深了对求解的过程，应先怎么办，再怎么做有了了解。还有在本次实验中对于MATLAB的使用有了一定的基础，虽然时间很短，但在实验之后对于用MATLAB软件来

25、求解电力系统中问题总是有一定的方便性认识更深。八、参考资料电力系统分析何仰赞华中科技大学出版社电力系统稳态分析陈珩中国电力出版社7.电力系统暂态分析李光琦中国电力出版社8.电力系统计算水利电力出版社实验素材：南昌大学实验报告学生姓名：学号：专业班级：实验类型：口验证综合设计口创新实验日期：12.7实验成绩：一、实验项目名称电力系统短路计算实验二、实验目的与要求：目的：通过实验教学加深学生的基本概念，掌握电力系统的特点，使学生通过系统进行物理模拟和数学模拟，对系统进行电力系统计算和仿真实验，以达到理论联系实际的效果。通过电子计算机对电力系统短路等计算的数学模拟，分析电力系统的故障计算方法、实现工

26、程计算的功能。提高处理电力系统工程计算问题的实际能力，以及实现对电力系统仿真的过程分析。要求：1、使学生掌握对电力系统进行计算、仿真试验的方法，了解实验对电力系统分析研究的必要性和意义。2、使学生掌握使用实验设备计算机和相关计算软件、编程语言。3、应用电子计算机完成电力系统的短路计算。4、应用电子计算机及相关软件对电力系统进行仿真。三、主要仪器设备及耗材每组计算机1台、相关计算软件1套四、实验步骤将事先编制好的形成电力网数学模型的计算程序原代码由自备移动存储设备导入计算机。在相应的编程环境下对程序进行组织调试。应用计算例题验证程序的计算效果。对调试正确的计算程序进行存储、打印。完成本次实验的实

27、验报告。五、实验数据及处理结果运行自行设计的程序，把结果与例题的计算结果相比较，验证所采用的短路电流计算方法及程序运行的正确性。如果采用的是近似计算方法，还需分析由于近似所产生的误差是否在运行范围内。实验程序：clearclc;z=0.2i,inf,0.51i,inf;inf,4i,0.59i,inf;0.51i,0.59i,inf,1.43i;inf,inf,1.43i,inf;y=0,0,0,0;0,0,0,0;0,0,0,0;0,0,0,0;f=4;Y=zeros(4,4);for(i=1:4),for(j=1:4),ifi=jY(i,j)=Y(i,j)elseY(i,j)=-1.0/z

28、(i,j)endendendfor(i=1:4),for(j=1:4),Y(i,i)=Y(i,i)+y(i,j)+1.0/z(i,j)endendZ=inv(Y);If=1/Z(f,f);实验结果：If=0-0.48902i实验结果：If=0-0.48902i0-6.9608iD+1.9603i0-6.9608iD+1.9603i0-0.0-I.S449i0*1.6949i0*IJSOBi0+rem0-3550i0+0上刘3i0t0.6993i0-Q删i实验例题所给结果短路电流：If=-j0.4895，与程序运行结果在误差允许范围之内，故验证了该程序的正确性。六、思考讨论题或体会或对改进实验的

29、建议1.理解课本上讲述的同步电机突然三相短路的物理分析。答：同步电机稳态对称运行（包括稳态对称短路）时，电枢磁势大小不随时间变化，而在空间以同步速旋转，同转子没有相对运动，故不会在转子绕组中感应电流。突然短路时，定子电流发生急剧的变化，电枢反应磁通也随着变化，定转子间电流会相互影响，这是同步电机突然短路暂态过程区别于稳态短路的显著特点。我们在进行磁感分析时，对于每个绕组都遵守磁链守恒原则。对于无阻尼绕组同步电机突然三相短路时，短路后定子侧将出现：基频电流，由三相对称绕组的基频电流产生的交变磁链，用以抵消转子主磁场对定子各相绕组产生的交变磁链；直流i：对各绕组产生的不变屮、屮、屮，来apa0b0c0维持定子绕组的磁链初值不变;倍频电流i:定子各相直流产生的恒定磁势，20当转子旋转时，因转子纵横轴磁阻不同，转子每转过1800，磁阻经历一个变化周期，为适应磁阻的变化，产生倍频电流与直流共同作用，才能维持定子侧磁链初值不变。转子侧产生：附加直流分量Ai:为抵消定子电流对转子产生的强烈电枢反应fa影响，维持磁链不变，该附加直流与原直流同向，加强了励磁绕组的磁场，而且Ai产生磁通的一部分也要穿过定子绕组，故激起定子基频电流大大超过稳态短fa路电流转子基频交流：为抵消定子的直流和倍频电流产