CN108733910A - A method of realizing POWER SYSTEM STATE equation automatic deduction - Google Patents

A method of realizing POWER SYSTEM STATE equation automatic deduction Download PDF

Info

Publication number
CN108733910A
CN108733910A CN201810470024.1A CN201810470024A CN108733910A CN 108733910 A CN108733910 A CN 108733910A CN 201810470024 A CN201810470024 A CN 201810470024A CN 108733910 A CN108733910 A CN 108733910A
Authority
CN
China
Prior art keywords
single order
universal model
order link
output
state equation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810470024.1A
Other languages
Chinese (zh)
Other versions
CN108733910B (en
Inventor
孙志媛
刘光时
孙艳
周柯
刘默斯
李明珀
窦骞
梁水莹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Electric Power Research Institute of Guangxi Power Grid Co Ltd
Original Assignee
Electric Power Research Institute of Guangxi Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Electric Power Research Institute of Guangxi Power Grid Co Ltd filed Critical Electric Power Research Institute of Guangxi Power Grid Co Ltd
Priority to CN201810470024.1A priority Critical patent/CN108733910B/en
Publication of CN108733910A publication Critical patent/CN108733910A/en
Application granted granted Critical
Publication of CN108733910B publication Critical patent/CN108733910B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Abstract

The invention discloses a kind of methods for realizing POWER SYSTEM STATE equation automatic deduction, are related to electric system Small signal stability analysis technical field.The method, using that can identify that the join algorithm of each single order link universal model connects each single order link universal model, simplifies the derivation of state equation and output equation by the way that high-order link is converted into multiple single order link universal models;Simultaneously, single order link universal model, high-order link universal model and join algorithm can be write and be stored by computer program, once completing the programming of join algorithm and common single order link universal model, these join algorithm and common single order link universal model can be reused, the automatic deduction of distinct device state equation is carried out;When there is new equipment in system, without deriving again manually, existing single order link universal model, high-order link universal model and join algorithm can equally be utilized to complete the automatic deduction of its state equation.

Description

A method of realizing POWER SYSTEM STATE equation automatic deduction
Technical field
The invention belongs to the small interference of electric system Small signal stability analysis technical field more particularly to a kind of electric system is steady The automatic deduction method of each equipment state equation in setting analysis.
Background technology
It is continuously improved with the complexity of electric system, the access of especially extensive novel energy, power electronic element Hypersynchronous and extensive trans-regional long distance powedr transmission increase, lead to Study of Power System Small Disturbance problem, especially thus spread out The low-frequency oscillation of electric power system that bears, sub-synchronous oscillation problem, are paid attention to by electric system researcher further.
Study of Power System Small Disturbance analysis theories are established on the basis of dynamical system Method of Stability Analysis, its basic thought It is in the way of each element model in electric system and basic parameter, to construct the differential equation and the generation of whole system Number equation.According to Liapunov First Law, linearize the differential equation and algebraic equation to obtain the state side of whole system Journey is finally analyzed according to the characteristic root of state matrix in state equation to judge the small signal stability of system.
Dynamic device for the small interference analysis modeling of electric system generally has:Generator, magnetizing exciter, power system stability Device (Power System Stabilizer, PSS), prime mover, governor etc..Constitute the dynamic variable of their state equations Referred to as state variable, remaining variables are known as algebraic variable, and general device has its fixed state variable and its state equation.So And when the equipment that electric system is related to is more, the model of equipment is wide, and when equipment operates in different conditions because parameters default leads to it Link increases and decreases, and causes state variable quantity uncertain, and state equation also changes therewith.Therefore, it is necessary to be directed to distinct device, same The different model of one equipment, the different operating conditions of same equipment can obtain best suiting actual conditions to derive state equation POWER SYSTEM STATE equation.
All the time, the method for obtaining equipment state equation in electric system is directly pushed away by hand according to the block diagram of equipment It leads, the slave system block diagram that basic skills refers to Zhu Xiyan propositions directly seeks state equation;In the publication, first, biography Delivery function turns to each Dynamic mode only form containing integral element, and the coefficient of integral element is turned to 1;Then, exist The exit of each integral element sets a state variable, from the block diagram of equipment each state variable of the equipment and The relationship of other link state variables and algebraic variable, writes out the equation between them, and abbreviation arranges, you can writes out the equipment State equation.When there are when the integral element of high-order, then need to be broken down into multiple independent integral element connection in series-parallel for system The case where, in the block diagram before being included in, then with the method for observation their state equation is write out according to connection type.Hand The method of work is simply direct, for only having the equipment block scheme of several state variables that can acquire its state equation easily, still For existing compared with the large-scale block diagram for connecting complexity between multi-state variable, high step integration link, links, manual derivation method It is complicated for operation, and easily make mistakes.When being related to multiple models, multiple models and multiple operation conditions, this method needs Machinery repeats, and is repeatedly derived according to actual conditions, inefficiency.
The novel device being continuously added in face of electric system, the model of existing equipment and its simple state equation can not expire The needs of the small interference analysis of foot nowadays electric system.And carry out high-order state equation derive, especially face high-volume equipment, In the case of polytypic, more operation conditions, the time-consuming and laborious and error rate of original derivation is high, does not adapt to current electric system Development.Therefore, in today of computerization, become using the state equation of computer derivation power equipment with a set of algorithm One very significant thing.
Invention content
In view of the deficiencies of the prior art, the present invention provides a kind of method for realizing POWER SYSTEM STATE equation automatic deduction, This method proposes to indicate the general Link Model of all link state equations and output equation, and common single order link is converted For universal model;Based on the universal model, then a set of algorithm for making general link be connected is write, which can identify general Model connects inverted each link universal model of good equipment, just according to the connection of equipment block scheme using the algorithm The state equation and output equation of the equipment can be calculated.
The present invention is to solve above-mentioned technical problem by the following technical solutions:A kind of realization POWER SYSTEM STATE side The method of journey automatic deduction, including following steps:
(1) equipment is often converted into single order link universal model with single order link, calculates the single order link universal model State equation and output equation, and its state equation and output equation are programmed, are stored;
(2) derive and write the join algorithm of single order link universal model;
(3) when there are high-order link, the single order link universal model derived using the step (1) constructs high-order Link universal model, calculates the state equation and output equation of the high-order link universal model, and to its state equation and defeated Go out equation to be programmed;To common high-order link, the corresponding state equation of high-order link universal model and output equation are compiled Journey, storage;
(4) it according to the abbreviation of equipment transmission function block diagram sequence, is sequentially connected by the join algorithm of the step (2) Each block diagram, automatic deduction go out the state equation of the equipment;
(5) it when having novel device or other model block diagrams of equipment, repeats step (3) (4) and derives novel device or be somebody's turn to do The state equation of other models of equipment.
Further, the state equation of step (1) the single order link universal model and output equation use following form It indicates:
Px=Ax+Bu
Y=Cx+Du
Wherein, p indicates that differential operator, x indicate state variable;U indicates that input variable, y indicate that output variable, A indicate shape State matrix, B indicate that input matrix, C indicate that output matrix, D indicate feedforward matrix.The state equation and output equation can indicate The state equation and output equation of all links.
Further, the computational methods of state matrix A, input matrix B, output matrix C and feedforward matrix D are:It will be normal The transmission function form of corresponding single order link universal model is converted to the transmission function of single order link, it is logical further according to single order link With model state equation and output equation, A, B, C, D are determined with the method for undetermined coefficients.
Further, the transmission function form of the single order link universal model is indicated using following form:
Or
Wherein, G (s) indicates transmission function, A, B, C, D be respectively state matrix in original state equation and output equation, Input matrix, output matrix and feedforward matrix.
Further, the step (2) is to carry out single order link universal model after determining condition of contact according to connection type The derivation of join algorithm, if single order link universal model at least there are two, respectively single order link universal model one and single order ring Universal model two is saved, the definition of connection type and its condition of contact are:
(2.1) it connects:The output end of single order link universal model one is connected to the input of single order link universal model two End;Its condition of contact is the output variable y of single order link universal model one1Equal to the input variable of single order link universal model two u2
(2.2) phase opposite sign is taken:The output variable y of single order link universal model one1It is inputted for single order link universal model two Variable u1Opposite number, condition of contact y1=-ku1
(2.3) addition input:Single order link universal model one increases an input signal, and it is u to enable former input variable1, add It is u to add input variable2, condition of contact is new input variable u1'=u1+u2
(2.4) feedback link:By two single order link universal models by two amplification factors be p12iAnd p21iRatio is put Big link joins end to end, if the input variable and output variable of single order link universal model one are respectively u before connection1And y1, single order The input variable and output variable of link universal model two are respectively u2And y2;External signal ur1It is p by amplification factor1iRatio Example amplifying element input signal is to single order link universal model one, external signal ur2It is p by amplification factor2iRatio enlargement Link input signal is to single order link universal model two;The output signal of single order link universal model one is p by amplification factor1o Ratio enlargement link output be yr1, the output signal of single order link universal model two is p by amplification factor2oRatio put Big link exports yr2;Its condition of contact is the new input variable u of single order link universal model one1'=p1iur1+p21iy2, single order ring Save the new input variable u of universal model two2'=p2iur2+p12iy1;The output variable y of entire reponse systemr1=P1oy1, yr2= P2oy2
Above four connection types are the common connection type of equipment block scheme, other connection types and condition of contact can be with It similarly derives, computer programming is used in combination to realize.
Further, it is to split into high-order link more to the building method of the step (3) high-order link universal model The form that a single order link is connected utilizes the connection of the step (2) then according to the split single order link order of connection Multiple single order link universal models are connected and constitute high-order link Universal Die by algorithm and corresponding single order link universal model Type.If high-order link is more complicated, the transmission function in its block diagram can be utilized to derive, but the link in block diagram is general It is relatively low, be not in the link of higher-order.
Further, for the step (4) abbreviation sequence, if presence signal is intersected in block diagram, first movable signal Summing point, leading point cancellation signal are intersected, then utilize the single order link universal model that the step (2) and (3) have derived with High-order link universal model, according to from inner looping to the sequence of external loop piecemeal gradually abbreviation, using join algorithm by each ring Section links together, and ultimately generates the state equation of the equipment.
With in the prior art, each state variable relationship is listed by observation device block diagram manually, then derive and do well The manual derivation method of equation, a kind of method for realizing POWER SYSTEM STATE equation automatic deduction proposed by the present invention, pass through by High-order link is converted into multiple single order link universal models, recycles the join algorithm that can identify each single order link universal model Each single order link universal model is connected, the derivation of state equation and output equation is simplified;Meanwhile single order link is general Model, high-order link universal model and join algorithm can be write and be stored by computer program, once the company of completion Connect the programming of algorithm and common single order link universal model, it will be able to reuse these join algorithm and common single order link is logical With model, the automatic deduction of distinct device state equation is carried out;When there is new equipment in system, without deriving again manually, together Sample can complete its state equation using existing single order link universal model, high-order link universal model and join algorithm Automatic deduction;
The method of the present invention has the advantages that efficient, highly reliable, favorable expandability, facilitates large-scale calculations.When in face of working as The case where modern Electrical Power System Dynamic device category is more, operating condition is complicated and novel electric power electric equipment is largely put into, calculate Machine program is the effective tool for the state equation automatic deduction for realizing large-scale electrical power system, and it is small dry that it to write electric system Stability Calculation is disturbed to greatly improve with the efficiency and accuracy for analyzing program.
Description of the drawings
It, below will be to attached drawing needed in embodiment description in order to illustrate more clearly of technical scheme of the present invention It is briefly described, it should be apparent that, the accompanying drawings in the following description is only one embodiment of the present of invention, general for this field For logical technical staff, without creative efforts, other drawings may also be obtained based on these drawings.
Fig. 1 is the present embodiment single order link universal model block diagram;
Fig. 2 is that two link universal models of the present embodiment are connected in series with block diagram;
Fig. 3 is the present embodiment single order link universal model addition input block diagram;
Fig. 4 is the present embodiment universal model feedback link block diagram;
Fig. 5 is the present embodiment IEEE DC1A direct-current exciter exemplary block diagrams.
Specific implementation mode
With reference to the attached drawing in the embodiment of the present invention, the technical solution in the present invention is clearly and completely described, Obviously, described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Based in the present invention Embodiment, the every other embodiment that those of ordinary skill in the art are obtained under the premise of not making creative work, It shall fall within the protection scope of the present invention.
A kind of method for realizing POWER SYSTEM STATE equation automatic deduction provided by the present invention, including following step Suddenly:
(1) equipment is often converted into single order link universal model with single order link, calculates the single order link universal model State equation and output equation, and its state equation and output equation are programmed, are stored.
The state equation and output equation of step (1) single order link universal model are indicated using following form:
Wherein, p indicates that differential operator, x indicate state variable;U indicates that input variable, y indicate that output variable, A indicate shape State matrix, B indicate that input matrix, C indicate that output matrix, D indicate feedforward matrix;And with following symbolic indication variable number:
nst:State variable number;nin:Input variable number;nout:Output variable number.
From the foregoing, it will be observed that matrix A is nst×nstThe square formation of dimension, matrix B nst×ninThe matrix of dimension, Matrix C nout ×nstThe matrix of dimension, matrix D nout×ninThe matrix of dimension.
State matrix A, input matrix B, output matrix C and feedforward matrix D computational methods be:By common single order link Transmission function be converted to the block diagram and transmission function form (as shown in Figure 1) of corresponding single order link universal model, further according to Single order link universal model state equation and output equation (formula (1)), A, B, C, D related coefficient square are determined with the method for undetermined coefficients Battle array.
The transmission function of block diagram according to figure 1, single order link universal model is indicated using following form:
Other single order links can be converted to single order link universal model form using formula (2):By all common single orders The transmission function of link turns to the transmission function form of single order link universal model, its state equation is found out with the method for undetermined coefficients Coefficient expressions, and in the state equation program storage to computer that coefficient is determined.
For example, for delay component, transmission function is:
The transmission function form for being turned to universal model is:
Formula (4) and formula (2) are compared, the coefficient in corresponding universal model is obtainedD=0 can be enabledC=K obtains delay component generic state equation and output equation form:
It is corresponding:nst=1, nin=1, nout=1.
Alternatively, it is also possible to enableC=1 obtains another state equation and output equation shape of the delay component Formula:
This variation has no effect on end-state formulation character value result.
When block diagram ignores the delay component that original possesses, the parameter T of the link is zero or default, at this time its state Equation is:
Y=Ku (7)
Then corresponding A, B, C are empty matrix, D=K, corresponding nst=0, nin=1, nout=1.
For lead-lag link, transmission function is:
Transformation obtains the transmission function form of its universal model:
Through-beam Series (2) have:It can then enable: Obtain the universal model state equation and output equation form of lead-lag link:
It is corresponding:nst=1, nin=1, nout=1.
When block diagram ignores the lead-lag link that original possesses, the parameter T of the linknAnd TdIt is simultaneously zero, this When its state equation and output equation be:
Y=Ku (11)
Then corresponding A, B, C are empty matrix, D=K, corresponding nst=0, nin=1, nout=1.
To unit integral element, transmission function is:
It is transformed to universal model transmission function form:
Contrast (2) have:C=1, D=0 obtain its universal model state equation and output Equation form:
It is corresponding:nst=1, nin=1, nout=1.
Because link cannot lack thus, TeCannot be zero.
For feedback element, transmission function is:
It is transformed to the transmission function form of universal model:
Contrast (2) have:
It can then enable:C=-1 obtains its universal model state equation and output equation form:
It is corresponding:nst=1, nin=1, nout=1.
This link cannot lack, therefore TfCannot be zero.
It is the expression-form of common single order link universal model state equation and output equation, other single order links above Its universal model expression formula can be derived with the method, and writes program storage.For common single order link, universal model State equation and output equation expression formula can the link occur place recycling.
(2) derive and write the join algorithm of single order link universal model.
The derivation that progress single order link universal model join algorithm after condition of contact is determined according to connection type, if single order ring Section universal model has multiple, respectively single order link universal model one, single order link universal model two, single order link universal model Three, single order link universal model four ..., the definition of connection type and its condition of contact are:
(2.1) it connects:The output end of single order link universal model one is connected to the input of single order link universal model two End;Its condition of contact is the output variable y of single order link universal model one1Equal to the input variable of single order link universal model two u2;Its connection as shown in Fig. 2, the universal model of two links of simultaneous state equation and output equation, obtain:
It derives:
Then new universal model state equation and output equation are:
Then corresponding to variable number is:
nst=nst1+nst2
nin=nin1
nout=nout2
Consider that universal model therein has the case where default link:
Work as nst1、nst2When being not zero, have:
Work as nst1、nst2When one of which is zero, then matrix, which removes, lacks the corresponding part of dimension, such as:
If a):nst1=1, nst2=0, nst=nst1+nst2=nst1.N should be removed at this timest2Corresponding matrix part, i.e.,:
If b):nst1=0, nst2=1, nst=nst1+nst2=nst2.N should be removed at this timest1Corresponding matrix part, i.e.,:
If c):nst1=0, nst2=0, n at this timest=nst1+nst2=0.Then have: ([] indicates empty matrix).
(2.2) phase opposite sign is taken:The output variable y of single order link universal model one1It is inputted for single order link universal model two Variable u1Opposite number, condition of contact y1=-ku1
The output variable y of universal model1' it is former output variable y1Opposite number, i.e., its condition of contact be y1'=- y1, have:
It derives:
Therefore have:
Relevant variable number is constant.
(2.3) addition input:Single order link universal model one increases an input signal, and it is u to enable former input variable1, add It is u to add input variable2, condition of contact is new input variable u1'=u1+u2, as shown in figure 3, then:
Therefore have:
Relevant variable number is:
nst=nst
nin=nin+1
nout=nout
(2.4) feedback link:By two single order link universal models by two amplification factors be p12iAnd p21iRatio is put Big link joins end to end, if the input variable and output variable of single order link universal model one are respectively u before connection1And y1, single order The input variable and output variable of link universal model two are respectively u2And y2;External signal ur1It is p by amplification factor1iRatio Example amplifying element input signal is to single order link universal model one, external signal ur2It is p by amplification factor2iRatio enlargement Link input signal is to single order link universal model two;The output signal of single order link universal model one is p by amplification factor1o Ratio enlargement link output be yr1, the output signal of single order link universal model two is p by amplification factor2oRatio put Big link exports yr2;Its condition of contact is the new input variable u of single order link universal model one1'=p1iur1+p21iy2, single order ring Save the new input variable u of universal model two2'=p2iur2+p12iy1;The output variable y of entire reponse systemr1=P1oy1, yr2= P2oy2, connection is as shown in Figure 4.
If two of which universal model S1、S2Equation be:
According to condition of contact, the equation of each tie point is:
It can be released by formula (27) and formula (28):
Wherein:I is unit diagonal matrix.
New state equation is:
It derives:
Above four connection types are the common connection type of equipment block diagram, other connection types can similarly derive and It is realized with computer programming.
(3) when there are high-order link, the single order link universal model derived using the step (1) constructs high-order Link universal model, calculates the state equation and output equation of the high-order link universal model, and to its state equation and defeated Go out equation to be programmed;To common high-order link, the corresponding state equation of high-order link universal model and output equation are compiled Journey, storage.
The building method of high-order link universal model is that high-order link is split into the shape that multiple single order links are connected Formula utilizes the join algorithm and its single order link Universal Die of step (2) then according to the split single order link order of connection Multiple single order link universal models are connected and composed high-order link universal model by type.
Specific implementation is as follows:
If certain link exists such asTransmission function can be then split asWithTwo concatenated forms of delay component indicate the two delay components with universal model, you can are connected with series connection It connects algorithm and calculates the universal model after the connection of the two delay component universal models, i.e. the high-order link universal model of the link.
If there are still the link of higher-order, need to use the method in control theory by transmission function row write state equation (such as:The quick Theories of Automatic Control of uncle Zou:3rd edition [M] China Machine Press, 2007:386-395.), the shape of the link is found out State equation and output equation.The link of the general rare high-order of equipment block scheme, the method will not generally be used.
(4) it according to the abbreviation of equipment transmission function block diagram sequence, is sequentially connected by the join algorithm of the step (2) Each block diagram, automatic deduction go out the state equation of the equipment.It is as follows:
For the step (4) abbreviation sequence, if presence signal is intersected in block diagram, first movable signal summing point draws Go out a cancellation signal to intersect, the single order link universal model and high-order link for then utilizing the step (2) and (3) to derive are logical With model, to the sequence of external loop piecemeal gradually abbreviation, links one is connected to using join algorithm according to from inner looping It rises, ultimately generates the state equation of the equipment.
IEEE DC1A direct-current exciter models as shown in Figure 5 then utilize computer to realize the derivation step of its state equation Suddenly it is:
(4.1) it generates in first feedback loop on forward passThe Universal Die of link Type, and they are cascaded using formula (20), generate the concatenated universal model S of these three links;
(4.2) feedback element formula (23) will be had been converted into universal modelIt takes negative, obtains's Universal model;
(4.3) P in formula (32) is enabled2i、P2oIt is 0, P12i、P21i、P1i、P1oIt is 1, (32) is then utilized will to have turned before Turn to universal modelThe universal model that feedback link generates S' is carried out with the universal model of S;
(4.4) utilize formula (26) to S' addition input variables VrefGenerate new universal model S ";
(4.5) formula (20) is utilized, it willUniversal model with before gained S " universal model connect, you can generation The universal model of DC1A type magnetizing exciters is to get to its corresponding state equation and output equation.
(5) it when having novel device or other model block diagrams of equipment, repeats step (3) (4) and derives novel device or be somebody's turn to do The state equation of other models of equipment.
It can be seen that the join algorithm of the present invention and the universal model of most of links once derives programming can be repeatedly sharp With effectively preventing deriving the problem of efficiency is low, poor reliability of equipment state equation by hand.Utilize the equipment of computer derivation State equation, flexibility is high, autgmentability is strong, easily processing equipment can ignore link, newly-increased ring under different operating conditions Situations such as section, is effectively improved accuracy and efficiency that electric system small interference stability is calculated and analyzed.
Above disclosed is only the specific implementation mode of the present invention, but scope of protection of the present invention is not limited thereto, Any one skilled in the art in the technical scope disclosed by the present invention, can readily occur in variation or modification, It is covered by the protection scope of the present invention.

Claims (7)

1. a kind of method for realizing POWER SYSTEM STATE equation automatic deduction, which is characterized in that including following steps:
(1) equipment is often converted into single order link universal model with single order link, calculates the shape of the single order link universal model State equation and output equation, and its state equation and output equation are programmed, it stores;
(2) derive and write the join algorithm of single order link universal model;
(3) when there are high-order link, the single order link universal model derived using the step (1) constructs high-order link Universal model calculates the state equation and output equation of the high-order link universal model, and to its state equation and output side Journey is programmed;To common high-order link, by the corresponding state equation of high-order link universal model and output equation programming, deposit Storage;
(4) it according to the abbreviation of equipment transmission function block diagram sequence, is sequentially connected by the join algorithm of the step (2) each Block diagram, automatic deduction go out the state equation of the equipment;
(5) it when having novel device or other model block diagrams of equipment, repeats step (3) (4) and derives novel device or the equipment The state equation of other models.
2. a kind of method for realizing POWER SYSTEM STATE equation automatic deduction according to claim 1, which is characterized in that institute The state equation and output equation for stating step (1) single order link universal model are indicated using following form:
Px=Ax+Bu
Y=Cx+Du
Wherein, p indicates that differential operator, x indicate state variable;U indicates that input variable, y indicate that output variable, A indicate state square Battle array, B indicate that input matrix, C indicate that output matrix, D indicate feedforward matrix.
3. a kind of method for realizing POWER SYSTEM STATE equation automatic deduction according to claim 2, which is characterized in that shape State matrix A, input matrix B, output matrix C and feedforward matrix D computational methods be:By the transmission function of common single order link It is converted to the transmission function form of corresponding single order link universal model, further according to single order link universal model state equation and output Equation determines A, B, C, D with the method for undetermined coefficients.
4. a kind of method for realizing POWER SYSTEM STATE equation automatic deduction according to claim 3, which is characterized in that institute The transmission function form for stating single order link universal model is indicated using following form:
Or
Wherein, A, B, C, D are respectively state matrix in original state equation and output equation, input matrix, output matrix and preceding Present matrix.
5. a kind of method for realizing POWER SYSTEM STATE equation automatic deduction according to claim 1, which is characterized in that institute State step (2) be the derivation that single order link universal model join algorithm is carried out after condition of contact is determined according to connection type, if one Rank link universal model at least there are two, respectively single order link universal model one and single order link universal model two, connection side The definition of formula and its condition of contact are:
(2.1) it connects:The output end of single order link universal model one is connected to the input terminal of single order link universal model two;Its Condition of contact is the output variable y of single order link universal model one1Equal to the input variable u of single order link universal model two2
(2.2) phase opposite sign is taken:The output variable y of single order link universal model one1For two input variable u of single order link universal model1 Opposite number, condition of contact y1=-ku1
(2.3) addition input:Single order link universal model one increases an input signal, and it is u to enable former input variable1, addition input Variable is u2, condition of contact is new input variable u1'=u1+u2
(2.4) feedback link:By two single order link universal models by two amplification factors be p12iAnd p21iRatio enlargement ring Section joins end to end, if the input variable and output variable of single order link universal model one are respectively u before connection1And y1, single order link The input variable and output variable of universal model two are respectively u2And y2;External signal ur1It is p by amplification factor1iRatio put Big link input signal is to single order link universal model one, external signal ur2It is p by amplification factor2iRatio enlargement link Input signal is to single order link universal model two;The output signal of single order link universal model one is p by amplification factor1oRatio Example amplifying element output is yr1, the output signal of single order link universal model two is p by amplification factor2oRatio enlargement ring Section output yr2;Its condition of contact is the new input variable u of single order link universal model one1'=p1iur1+p21iy2, single order link is logical With the new input variable u of model two2'=p2iur2+p12iy1;The output variable y of entire reponse systemr1=P1oy1, yr2=P2oy2
6. a kind of method for realizing POWER SYSTEM STATE equation automatic deduction according to claim 1, which is characterized in that right The building method of high-order link universal model is in the step (3), and high-order link, which is split into multiple single order links, to be connected Form, then according to the split single order link order of connection, the join algorithm using the step (2) and corresponding single order Multiple single order link universal models are connected and constitute high-order link universal model by link universal model.
7. a kind of method for realizing POWER SYSTEM STATE equation automatic deduction according to claim 1, which is characterized in that right In the step (4) abbreviation sequence, if presence signal is intersected in block diagram, first movable signal summing point, leading point eliminate letter Number intersect, the single order link universal model and high-order link universal model for then utilizing the step (2) and (3) to derive are pressed , to the sequence of external loop piecemeal gradually abbreviation, links are linked together using join algorithm, most throughout one's life according to from inner looping At the state equation of the equipment.
CN201810470024.1A 2018-05-16 2018-05-16 Method for realizing automatic derivation of state equation of power system Active CN108733910B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810470024.1A CN108733910B (en) 2018-05-16 2018-05-16 Method for realizing automatic derivation of state equation of power system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810470024.1A CN108733910B (en) 2018-05-16 2018-05-16 Method for realizing automatic derivation of state equation of power system

Publications (2)

Publication Number Publication Date
CN108733910A true CN108733910A (en) 2018-11-02
CN108733910B CN108733910B (en) 2022-04-29

Family

ID=63938479

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810470024.1A Active CN108733910B (en) 2018-05-16 2018-05-16 Method for realizing automatic derivation of state equation of power system

Country Status (1)

Country Link
CN (1) CN108733910B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113128072A (en) * 2021-05-13 2021-07-16 清鸾科技(成都)有限公司 High-precision transfer function simulation method and device, storage medium and electronic equipment

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103793563A (en) * 2014-01-24 2014-05-14 广西电网公司电力科学研究院 Wind power plant equivalent simulation module, AC/DC locomotive and power supply system of AC/DC locomotive
US20150212498A1 (en) * 2012-11-30 2015-07-30 Fuji Electric Co., Ltd. Design method of control device and control device
CN106503346A (en) * 2016-10-31 2017-03-15 浙江大学 The generating set electromagnetical transient emulation method that a kind of quick and smooth starts
CN107437803A (en) * 2017-06-30 2017-12-05 广西大学 Based on the power system Small signal stability analysis method that fractional order differential is theoretical

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150212498A1 (en) * 2012-11-30 2015-07-30 Fuji Electric Co., Ltd. Design method of control device and control device
CN103793563A (en) * 2014-01-24 2014-05-14 广西电网公司电力科学研究院 Wind power plant equivalent simulation module, AC/DC locomotive and power supply system of AC/DC locomotive
CN106503346A (en) * 2016-10-31 2017-03-15 浙江大学 The generating set electromagnetical transient emulation method that a kind of quick and smooth starts
CN107437803A (en) * 2017-06-30 2017-12-05 广西大学 Based on the power system Small signal stability analysis method that fractional order differential is theoretical

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
GOMES, ED 等: "State Space Decoupling Control Design Methodology for Switching Converters", 《IEEE ENERGY CONVERSION CONGRESS AND EXPOSITION》 *
SUIJIAZHUANG11: "线性变换26状态空间表达式求传递", 《豆丁》 *
董玮 等: "根据微分方程构造SIMULINK仿真模型的研究", 《电气传动》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113128072A (en) * 2021-05-13 2021-07-16 清鸾科技(成都)有限公司 High-precision transfer function simulation method and device, storage medium and electronic equipment
CN113128072B (en) * 2021-05-13 2024-01-19 清鸾科技(成都)有限公司 Transfer function high-precision simulation method and device, storage medium and electronic equipment

Also Published As

Publication number Publication date
CN108733910B (en) 2022-04-29

Similar Documents

Publication Publication Date Title
Peng et al. New Approach on Robust Delay-Dependent ${H} _\infty $ Control for Uncertain T--S Fuzzy Systems With Interval Time-Varying Delay
CN102890735B (en) Modeling method of strategy table and strategy table searching match of stable control device
US11762345B2 (en) Process optimization by grouping mixed integer nonlinear programming constraints
CN105391057B (en) A kind of GPU threaded design methods that direction of energy Jacobi battle array calculates
CN106249717A (en) A kind of control method for coordinating based on the modeling of executor's saturated multi-agent system
Sambariya et al. A robust PID controller for load frequency control of single area re-heat thermal power plant using elephant herding optimization techniques
CN109657902B (en) Doubly-fed wind turbine dynamic equivalence method considering subsynchronous oscillation characteristics
CN106155985A (en) A kind of shortage of data fill method based on adjacent data feature
CN103761420A (en) Evaluation method for stepwise regression of thermal power equipment performances
CN108710724A (en) A kind of fuzzy double-response face method calculating leaf dish vibration reliability
CN108733910A (en) A method of realizing POWER SYSTEM STATE equation automatic deduction
CN106877769A (en) A kind of method of servomotor plus of speed controller parameter self-tuning
CN113036822B (en) Wind storage coordination control method, system, equipment and storage medium
CN108258725B (en) Doubly-fed wind turbine dynamic equivalence method based on equivalent power angle coherence
Pan et al. Reliability analysis of turbine blades based on fuzzy response surface method
CN110968958B (en) Wind power plant equivalence modeling method based on single-machine equivalence and selection modal analysis
CN108736508B (en) Voltage ride-through control simulation method and system
CN113722881B (en) Generator phase advance limit calculation method, equipment and medium based on parameter identification
CN111414573B (en) Input and output data processing method and system in engineering application field
CN101630940B (en) Fuzzy neural network inverse robust controller of induction motor speed regulation system and construction method
Selmi et al. A simple method for the steady state performances of self-excited induction generators
CN106099992A (en) Generator Governor model based on the actual units limits of therrmodynamic system and modeling
Li et al. Deep Reinforcement Learning-Based Fast Prediction of Strategies for Security Control
CN103558762A (en) Method for implementing immune genetic PID controller based on graphic configuration technology
CN101825867A (en) Equipment optimization operating system, optimal operation point computing method and optimal operation point computation program

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant