c语言牛拉法潮流（C language cow rafa tide）.docVIP

c语言牛拉法潮流（C language cow rafa tide） I would like to recommend to you a BBS, electric power research network. I am sure that you will be familiar with it. Please click on it. Thank you /? From user = jcxx01, there is no other purpose, because its hard to earn points in the research network, I just want to earn points # include stdio, h # include stdlib. H H # include math.h # define PI (3.14159) / / node parameter structure Struct NodeType { Int N; / / node number The int Type. / / node type Double e; / / voltage amplitude Double f; / / voltage phase Angle Double Pd; / / load work Double Qd. / / load no work Double Ps; / / work is done Double Qs. / / there is no work Double Bc; / / the reactance value of the shunt capacitance }; / / branch parameter structure Struct BranchType { Int Nbr; / / branch no. Int Nl. / / the first node Int Nr. / / small points Double R; / / branch resistance Double X; / / branch electric reactance Double Bn; / / to the earth reactance Double Kt. / / the branching ratio }; / / * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * global variable declaration * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Int n; / / the number of nodes Int nPQ; / / PQ node number Int nPV; / / number of PV nodes Int NBR; / / the way Int ng. / / generator table Int Mark = 0; Whether the / / flag branch parameter has been converted Double * * G; / / the G part of the admittance matrix Double * * B; / / the admittance matrix B Double * dS; / / power balance Mid1, double * * mid2; / / an intermediate variable in the measurement of power inequality Double * Us; / / voltage initial value Double error = 1; / / error values Double iteration = 0.000001; / / error accuracy Double * * Jacob; / / the jacobian Double * * invJac; The inverse of the jacobian matrix Double * dfe; / / node voltage correction value Struct NodeType * Node; / / the node parameter structure of the input time Struct BranchType * Branc