电力系统导论实验报告

1、. Experiment Introduction to power systems学生姓名： 学 号： 专业班级： 实验名称： 电力系统导论（双语） 2014年 6 月 5 日 ; CONTENTS1、EXPERIMENT 11Bus Admittance Matrix1-62、EXPERIMENT 26Bus Impedance Matrix6-133、EXPERIMENT 313Gauss-Seidel and newton method13-164、PERSENAL SUMMARY16Experiment 1Bus Admittance Matrix1. Objective To wr

2、ite a simple program in MATLAB® for the algorithm of bus admittance matrix. 2. System Requirement Computer with MATLAB® 6 or above installed. 3. Procedure 1.0 Launch the MATLAB program.2.0 Go to FILE NEW M-file. 3.0 Write a function Y = The_Node_Admittance_Matrix(TopoStructureAndBranchPara

3、) for the formation of the bus admittance matrix. 4.0 TopoStructureAndBranchPara is the transmission line, cable and transformer input data and contains five columns parameters. The first two columns are the line bus numbers and the remaining columns contain the line resistance and reactance in per-

4、unit and transformer tap ratio or capacitor of transmission line. 5.0 The function should return the bus admittance matrix. 4. Exercises Use the written function, Y = The_Node_Admittance_Matrix (TopoStructureAndBranchPara) to obtain the Ybus of the following power system network: Q1. You are require

5、d to write the Ybus topological structure and parameter into a text file. (Hint: use the matlab text compiler to write down the table 1 data, using the comma to separate the parameters, and save it use the name of 4_Power_System_Data.dbf) Q2. You are required to write out the program flow figure of

6、forming a nodal admittance matrix.Hint. You are required to compile a program to form the Ybus Matrix, the following program is a reference program to you. Figure : One-line diagram of power systemFor example ,from the textbook “power system analysis” No.2 edition 3 on page 6162Table 1：Transformer a

7、nd transmissssion Line dataFrom Bus#To Bus#R(p.u)X(p.u)B(p.u)or ratio KOthers120.10.4j0.015281300.31.1140.120.5j0.01920240.080.40J0.01413 NodalAdmittanceMatrix =     1.0421 - 8.2429i   -0.5882 + 2.3529i     

8、60;  0 + 3.6667i   -0.4539 + 1.8911i   -0.5882 + 2.3529i    1.0690 - 4.7274i         0              

9、     0                   0 + 3.6667i                       &

10、#160; 0  0 - 3.3333i             0      -0.4539 + 1.8911i         0      0   0.9346 - 4.2616i5.The flow chart

11、Figure : The flow chart of Forming Nodal Admittance Matrix15The program is:%function OutPut=The_Node_Admittance_Matrix(handles)%is a subroutine of PowerSystemCalculationfunction OutPut=The_Node_Admittance_Matrix(handles)%the following program is open a data file and get the Number of % Node and Bran

12、ch data to form a nodal addmittance matrix%the following code is open a file and read the data of power system networkfname,pname = uigetfile('*.dbf','Select the network parametre data-file');TopoStructureAndBranchPara= csvread(fname);NumberOfBranch,NumberOfPara=size(TopoStructureAnd

13、BranchPara);Temporary1=max(TopoStructureAndBranchPara(:,1);Temporary2=max(TopoStructureAndBranchPara(:,2); if Temporary1 > Temporary2 NumberOfNode=Temporary1; else NumberOfNode=Temporary2; end%The following program is to form the Nodal Admittance Matrix % and the Topologic structure and Branch Pa

14、rametres are arranged% I,J,R,X,C/K, and pay attention to the inpedence of transformer is in the% side of Node J and the ratio of transformer 1:K is in the side of Node Ifor CircleNumber1=1:NumberOfBranch for CircleNumber2=1:NumberOfBranch NodalAdmittanceMatrix(CircleNumber1,CircleNumber2)=0; endendf

15、or CircleNumber=1:NumberOfBranch if TopoStructureAndBranchPara(CircleNumber,5) > 0.85 NodalAdmittanceMatrix(TopoStructureAndBranchPara(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,1)=. NodalAdmittanceMatrix(TopoStructureAndBranchPara(TopoStructureAndBranchPar

16、a(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,1)+. TopoStructureAndBranchPara(CircleNumber,5)2/. (TopoStructureAndBranchPara(CircleNumber,3)+. j*TopoStructureAndBranchPara(CircleNumber,4) ; NodalAdmittanceMatrix(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(Circle

17、Number,2)=. NodalAdmittanceMatrix(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2)+. 1/(TopoStructureAndBranchPara(CircleNumber,3)+j*TopoStructureAndBranchPara(CircleNumber,4); NodalAdmittanceMatrix(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranc

18、hPara(CircleNumber,2)=. NodalAdmittanceMatrix(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,2). -TopoStructureAndBranchPara(CircleNumber,5)/. (TopoStructureAndBranchPara(CircleNumber,3)+j*TopoStructureAndBranchPara(CircleNumber,4); NodalAdmittanceMatrix(TopoStruc

19、tureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,1)=. NodalAdmittanceMatrix(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,2); else NodalAdmittanceMatrix(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,1)

20、=. NodalAdmittanceMatrix(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,1)+. +1/(TopoStructureAndBranchPara(CircleNumber,3)+. j*TopoStructureAndBranchPara(CircleNumber,4)+j*TopoStructureAndBranchPara(CircleNumber,5); NodalAdmittanceMatrix(TopoStructureAndBranchPar

21、a(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2)=. NodalAdmittanceMatrix(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2)+. +1/(TopoStructureAndBranchPara(CircleNumber,3)+. j*TopoStructureAndBranchPara(CircleNumber,4)+j*TopoStructureAndBranchPara(Circ

22、leNumber,5) NodalAdmittanceMatrix(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara( CircleNumber,2)=. NodalAdmittanceMatrix(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara( CircleNumber,2). -1/(TopoStructureAndBranchPara(CircleNumber,3)+. j*TopoStructureAnd

23、BranchPara(CircleNumber,4); NodalAdmittanceMatrix(TopoStructureAndBranchPara( CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,1)=. NodalAdmittanceMatrix(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara( CircleNumber,2); endendThe result is： NodalAdmittanceMatrix = 1.0421

24、 - 8.2429i -0.5882 + 2.3529i 0 + 3.6667i -0.4539 + 1.8911i -0.5882 + 2.3529i 1.0690 - 4.7274i 0 0 0 + 3.6667i 0 0 - 3.3333i 0 -0.4539 + 1.8911i 0 0 0.9346 - 4.2616iExperiment 2Power GridBus Impedance Matrix 1. Objective To write a simple program in MATLAB® for the algorithm of bus impedance mat

25、rix. 2. System Requirement Computer with MATLAB® 6 or above installed. 3. Procedure 1.0 Launch the MATLAB program.2.0 Go to FILE NEW M-file. 3.0 Write a function Z = znbus (z) for the formation of the bus impedance matrix. 4.0 z is the line input and contains three columns. The first two column

26、s are the line bus numbers and the remaining columns contain the line resistance in per-unit. 5.0 The function should return the bus impedance matrix. 4. Exercises Use the written function, Z = znbus(z) to obtain the Ybus of the following power system network: Example 1Figure 3: One-line diagram of

27、power systemFor example ,from the textbook “power system analysis” No.2 edition 3 on page 6162Table 1：Transformer and transmissssion Line dataFrom Bus#To Bus#R(p.u)X(p.u)B(p.u)or ratio KOthers120.10.4j0.015281300.31.1140.120.5j0.01920240.080.40J0.01413Q2. You are required to write the Zbus into a te

28、xt file. (Hint: use the matlab text compiler)Example 2 For the system shown, form Zbus matrix using the building algorithmSolution A line listApply Kron reduction to eliminate the last rowHint. You are required to compile a program to form the Zbus Matrix.the following program is a reference program

29、 to you. The program is:%function OutPut=The_Node_impedance_Matrix(handles)%is a subroutine of PowerSystemCalculationfunction OutPut=The_Node_impedance_Matrix(handles)%the following program is open a data file and get the Number of % Node and Branch data to form a nodal impedance matrix%the followin

30、g code is open a file and read the data of power system networkfname,pname = uigetfile('*.dbf','Select the network parametre data-file');Topo_Structure_And_Branch_Para= csvread(fname);%get the electric power system the number of branch and the parametre of% elements NumberOfBranch,Nu

31、mberOfPara=size(Topo_Structure_And_Branch_Para);%Temporary1-temporary variable 1%Temporary2-temporary variable 2Temporary1=max(Topo_Structure_And_Branch_Para(:,1);Temporary2=max(Topo_Structure_And_Branch_Para(:,2); if Temporary1 > Temporary2 NumberOfNode=Temporary1; else NumberOfNode=Temporary2;

32、end% The following program is to form the Nodal impedance Matrix % and the Topologic structure and Branch Parametres are arranged% I,J,R,X,C/K, and pay attention to the inpedence of transformer is in the% side of Node J and the ratio of transformer 1:K is in the side of Node % set the initial value

33、of Nodal Admittance Matrix to zerofor CircleNumber1=1:NumberOfNode for CircleNumber2=1:NumberOfNode Nodal_impedance_Matrix(CircleNumber1,CircleNumber2)=0; endendfor CircleNumber=1:NumberOfBranchif Topo_Structure_And_Branch_Para(CircleNumber,5) > 0.85Nodal_impedance_Matrix(Topo_Structure_And_Branc

34、h_Para(Topo_Structure_And_Branch_Para(CircleNumber,1),Topo_Structure_And_Branch_Para(CircleNumber,1)=. Nodal_impedance_Matrix(Topo_Structure_And_Branch_Para(Topo_Structure_And_Branch_Para(CircleNumber,1),Topo_Structure_And_Branch_Para(CircleNumber,1)+Topo_Structure_And_Branch_Para(CircleNumber,5)2/(

35、Topo_Structure_And_Branch_Para(CircleNumber,3)+. j*Topo_Structure_And_Branch_Para(CircleNumber,4) ; Nodal_impedance_Matrix(Topo_Structure_And_Branch_Para(CircleNumber,2),Topo_Structure_And_Branch_Para(CircleNumber,2)=. Nodal_impedance_Matrix(Topo_Structure_And_Branch_Para(CircleNumber,2),Topo_Struct

36、ure_And_Branch_Para(CircleNumber,2)+. 1/(Topo_Structure_And_Branch_Para(CircleNumber,3)+j*Topo_Structure_And_Branch_Para(CircleNumber,4); Nodal_impedance_Matrix(Topo_Structure_And_Branch_Para(CircleNumber,1),Topo_Structure_And_Branch_Para(CircleNumber,2)=. Nodal_impedance_Matrix(Topo_Structure_And_B

37、ranch_Para(CircleNumber,1),Topo_Structure_And_Branch_Para(CircleNumber,2). -Topo_Structure_And_Branch_Para(CircleNumber,5)/. (Topo_Structure_And_Branch_Para(CircleNumber,3)+j*Topo_Structure_And_Branch_Para(CircleNumber,4); Nodal_impedance_Matrix(Topo_Structure_And_Branch_Para(CircleNumber,2),Topo_St

38、ructure_And_Branch_Para(CircleNumber,1)=. Nodal_impedance_Matrix(Topo_Structure_And_Branch_Para(CircleNumber,1),Topo_Structure_And_Branch_Para(CircleNumber,2); else Nodal_impedance_Matrix(Topo_Structure_And_Branch_Para(CircleNumber,1),Topo_Structure_And_Branch_Para(CircleNumber,1)=. Nodal_impedance_

39、Matrix(Topo_Structure_And_Branch_Para(CircleNumber,1),Topo_Structure_And_Branch_Para(CircleNumber,1)+. +1/(Topo_Structure_And_Branch_Para(CircleNumber,3)+. j*Topo_Structure_And_Branch_Para(CircleNumber,4)+j*Topo_Structure_And_Branch_Para(CircleNumber,5); Nodal_impedance_Matrix(Topo_Structure_And_Bra

40、nch_Para(CircleNumber,2),Topo_Structure_And_Branch_Para(CircleNumber,2)=. Nodal_impedance_Matrix(Topo_Structure_And_Branch_Para(CircleNumber,2),Topo_Structure_And_Branch_Para(CircleNumber,2)+. +1/(Topo_Structure_And_Branch_Para(CircleNumber,3)+. j*Topo_Structure_And_Branch_Para(CircleNumber,4)+j*Top

41、o_Structure_And_Branch_Para(CircleNumber,5) Nodal_impedance_Matrix(Topo_Structure_And_Branch_Para(CircleNumber,1),Topo_Structure_And_Branch_Para( CircleNumber,2)=. Nodal_impedance_Matrix(Topo_Structure_And_Branch_Para(CircleNumber,1),Topo_Structure_And_Branch_Para( CircleNumber,2). -1/(Topo_Structur

42、e_And_Branch_Para(CircleNumber,3)+. j*Topo_Structure_And_Branch_Para(CircleNumber,4); Nodal_impedance_Matrix(Topo_Structure_And_Branch_Para( CircleNumber,2),Topo_Structure_And_Branch_Para(CircleNumber,1)=. Nodal_impedance_Matrix(Topo_Structure_And_Branch_Para(CircleNumber,1),Topo_Structure_And_Branc

43、h_Para( CircleNumber,2); endendformat shortNodal_impedance_Matrix*inv(Nodal_impedance_Matrix)运行结果：Nodal_impedance_Matrix = 1.0421e+000 -8.2429e+000i -5.8824e-001 +2.3529e+000i 0 +3.6667e+000i 0 -5.8824e-001 +2.3529e+000i 5.8824e-001 -2.3377e+000i 0 0 0 +3.6667e+000i 0 0 -3.3333e+000i 0 0 0 0 4.5386e

44、-001 -1.8719e+000iNodal_impedance_Matrix = 1.0421e+000 -8.2429e+000i -5.8824e-001 +2.3529e+000i 0 +3.6667e+000i -4.5386e-001 +1.8911e+000i -5.8824e-001 +2.3529e+000i 1.0690e+000 -4.7274e+000i 0 0 0 +3.6667e+000i 0 0 -3.3333e+000i 0 -4.5386e-001 +1.8911e+000i 0 0 9.3463e-001 -4.2616e+000ians = 1.0000

45、 + 0.0000i 0.0000 + 0.0000i 0.0000 - 0.0000i 0.0000 - 0.0000i -0.0000 - 0.0000i 1.0000 - 0.0000i -0.0000 + 0.0000i -0.0000 - 0.0000i -0.0000 - 0.0000i -0.0000 - 0.0000i 1.0000 - 0.0000i -0.0000 0 - 0.0000i 0 + 0.0000i 0.0000 - 0.0000i 1.0000 + 0.0000i以上就是对阻抗矩阵的验证，其和其逆相乘为单位对角矩阵Experiment 3Gauss-Seide

46、l Method 1. Objective To write a simple program in MATLAB® for the algorithm to solution of nonlinear algebraic equations; Known as the method of successive displacements. 2. Discussion The most common methods for solving nonlinear algebraic equations are Gauss-Seidel, Newtow-Rahpson, and quasi

47、-Newton-Raphson methods. We start with one dimensional equations and then generalize to n-dimensional equations. 3. Mathmatics modelConsider the nonlinear equation .The equation is broken into two parts thus:. We assume is an initial "guess" of the solution, then "refine" the sol

48、ution using: This process is repeated thus and on the iteration we have: . If this process is convergent, then the successive solutions approach a value which is declared as the solution. Thus if at some step we have: where e is the desired "accuracy", then we claim the solution has been f

49、ound to the accuracy specified.4. System Requirement Computer with MATLAB® 6 or above installed. 5. Procedure 1.0 Launch the MATLAB program.2.0 Go to FILE NEW M-file. 3.0 Write a function program of Gauss Seidel Method. 6. Exercises Example: Using the Gauss-Seidel method to obtain the roots of

50、the equation:First the equation is expressed in a different form thusAnd the iteration can proceed. Take a good look at the shape of the iterations! Below is the program showing the process graphically (later showing how to do it iteratively).7.The flow chart of Gauss Seidel method (Omitted) 8.Reference Program and result.程序是：clear allclcx0=0.5;n=0;while (abs(x03-6*x02+9*x0-4)>0.00001)x0=-(x03-6*x02-4)/9;y=x0;n=n+1;end结果是：n=1627 y=x0=0.99818clear allclcx0=2.5;n=0;while (abs(x03-6*x02+9*x0-4)>0.00001)x0=-(x03