A method and system for writing state equations of power simulation system

Abstract

The present invention relates to a method and system for writing state equations for an electric power simulation system, comprising: obtaining a transfer function of a control system in the electric power simulation system; determining a state equation of the control system in the electric power simulation system based on the transfer function; determining a correlation matrix between the electrical system and the control system in the electric power simulation system based on the correlation relationship between the electrical system and the control system in the electric power simulation system; and writing the state equation of the electric power simulation system using the state equation of the electrical system in the electric power simulation system, the state equation of the control system, and the correlation matrix between the electrical system and the control system in the electric power simulation system. When applied to electromagnetic transient simulation, the technical solution provided by the present invention obtains simulation results by solving the state equation of the electric power simulation system, thereby eliminating the time delay caused by the electromagnetic transient simulation gradually solving the electrical system and the control system, and improving the numerical stability and calculation accuracy of the electromagnetic transient simulation.

Description

Technical Field

[0001] The present invention relates to the field of electromagnetic transient simulation, and in particular to a method and system for writing state equations of an electric power simulation system. Background Art

[0002] like Figure 1 As shown in the figure, when performing electromagnetic transient simulation of power system, the electrical system and control system are usually modeled separately according to the physical properties of the components; then the electrical system and control system models are solved separately to obtain simulation results; this electromagnetic transient simulation method can simplify the solution method of electromagnetic transient simulation of power system.

[0003] However, the alternating iterative solution of the electrical system and the control system will produce a time step delay. This delay error will lead to problems such as inaccurate calculation accuracy and poor convergence performance of the electromagnetic simulation system. Summary of the Invention

[0004] In view of the shortcomings of the existing technology, the purpose of the present invention is to provide a method for writing state equations of an electric power simulation system. When applied to electromagnetic transient simulation, this method obtains simulation results by solving the overall state equations of the electric power simulation system, thereby eliminating the time delay caused by the gradual solution of the electrical system and control system in the electromagnetic transient simulation, and improving the numerical stability and calculation accuracy of the electromagnetic transient simulation.

[0005] The purpose of the present invention is achieved by adopting the following technical solutions:

[0006] The present invention provides a method for writing state equations of a power simulation system, wherein the method comprises:

[0007] Obtain the transfer function of the control system in the power simulation system;

[0008] determining a state equation of a control system in a power simulation system according to the transfer function;

[0009] Determine the correlation matrix between the electrical system and the control system in the power simulation system according to the correlation relationship between the electrical system and the control system in the power simulation system;

[0010] The state equation of the electric power simulation system is written by using the state equation of the electric system in the electric power simulation system, the state equation of the control system and the correlation matrix between the electric system and the control system in the electric power simulation system.

[0011] Preferably, the step of obtaining a transfer function of a control system in the power simulation system includes:

[0012] The transfer function G(s) of the control system in the power simulation system is determined as follows:

[0013]

[0014] Where, β i =b i -a i b m , b i is the i-order coefficient of the output of the control system in the power simulation system, i∈(1~m-1); a i is the i-order coefficient of the input quantity of the control system in the power simulation system; b m is the m-order coefficient of the output of the control system in the power simulation system; s is the Laplace operator, and m is the number of state variables of the control system.

[0015] Furthermore, determining the state equation of the control system in the power simulation system according to the transfer function includes:

[0016] Decompose the transfer function of the control system in the power simulation system and obtain the intermediate quantity of the transfer function of the control system in the power simulation system

[0017] Obtaining the Laplace inverse transform z(t) of the intermediate quantity z(s);

[0018] The state variable matrix ω=[ω1…ω i …ω m ] T , where ω i is the i-th element in the state variable matrix of the control system’s transfer function, ω i =[z(t)] (i-1) ;

[0019] The state equation of the control system in the power simulation system is determined as follows:

[0020]

[0021] Wherein, ω' is the first-order derivative matrix of the state variable matrix of the transfer function of the control system in the power simulation system; A is the state variable coefficient matrix of the transfer function of the control system in the power simulation system; B is the input coefficient matrix of the control system in the power simulation system; C' is the correlation coefficient matrix of the state variables and output quantities of the control system in the power simulation system; D is the correlation coefficient between the input quantity and output quantity of the control system in the power simulation system; p is the input quantity of the control system in the power simulation system; q is the output quantity of the control system in the power simulation system;

[0022] The state variable coefficient matrix A of the transfer function of the control system in the power simulation system is determined as follows:

[0023]

[0024] Where a0 is the constant term coefficient of the control system input in the power simulation system; a1 is the first-order term coefficient of the control system input in the power simulation system; a2 is the second-order term coefficient of the control system input in the power simulation system; a m-1 is the m-1-order coefficient of the input quantity of the control system in the power simulation system;

[0025] The input coefficient matrix B of the control system in the power simulation system is determined as follows:

[0026] B=[0 0 0 … 0 1] T

[0027] Where B is an m-order column vector;

[0028] The correlation coefficient matrix C' between the state variables and output quantities of the control system in the power simulation system is determined as follows:

[0029] C'=[β0 … β i … β m-1 ]

[0030] The correlation coefficient D between the input and output of the control system in the power simulation system is determined as follows:

[0031] D=b m

[0032] Where b m is the m-th order coefficient of the output of the control system in the power simulation system.

[0033] Preferably, the determining of the association matrix between the electrical system and the control system in the power simulation system according to the association relationship between the electrical system and the control system in the power simulation system includes:

[0034] Initialize the electrical system and control system correlation matrix P in the electrical system capacitance and current to control system correlation matrix in the power simulation system C , the electrical system inductance current to control system correlation matrix P L , the electrical system resistance and current control system correlation matrix P R are all m×n order zero element matrices, the control system control state quantity to electrical system correlation matrix Q is n×m order zero element matrix, the control system to electrical system capacitance current correlation matrix C C is an n×n order zero element matrix, the control system's inductance current correlation matrix C L is an n×n order zero element matrix and the control system's resistance and current correlation matrix C of the electrical system Ris an n×n order zero element matrix, where m is the number of state variables of the control system; n is the number of circuit nodes of the electrical system;

[0035] According to the signals of each branch flow control system in the electrical system, P is updated respectively. C 、P R and P L The mth row element of ;

[0036] Update Q and C respectively according to the signal flowing from the control system to the electrical system C 、C L and C R .

[0037] Furthermore, the signals of the control system of each branch in the electrical system are updated respectively according to the flow direction of P C 、P R and P L The elements of the mth row include:

[0038] If the signal flowing to the control system from branch ij in the electrical system is a capacitor current signal, then update P C The element P in the mth row and ith column C (m,i)=C ij , element P in row m and column j C (m,j)=-C ij ;

[0039] If the signal flowing from branch ij in the electrical system to the control system is the inductor current signal, then update P L The element in row m and column i The element in row m and column j

[0040] If the signal flowing to the control system from branch ij in the electrical system is a resistance current signal, then update P R The element in row m and column i The element in row m and column j

[0041] If the signal flowing from branch ij to the control system in the electrical system is the node voltage difference of branch ij, then update P R The element P in the mth row and ith column R (m,i)=1, element P in row m and column j R (m,j)=-1;

[0042] Among them, R ij is the resistance value of the resistance branch ij, L ij is the inductance value of the inductor branch ij, C ij is the capacitance value of the capacitor branch ij.

[0043] Furthermore, the signals flowing from the control system to the electrical system are updated respectively as Q and C C 、C L and C R ,include:

[0044] If the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update the i-th row element Q(i)=const·C' and the j-th row element Q(j)=-const·C' in Q;

[0045] If the signal flowing from the electrical system to the control system is a capacitor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C C The element C in row i and column i C (i,i)=const·D·C, element C in row i and column j C (i, j) = -const·D·C, the element C in the jth row and ith column C (j,i)=-const·D·C, element C in row j and column j C (j,j)=-const·D·C;

[0046] If the signal flowing from the electrical system to the control system is an inductor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0047] If the signal flowing from the electrical system to the control system is a resistance current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0048] If the signal flowing from the electrical system to the control system is the node voltage difference of the branch, the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update C R The element C in row i and column i R (i,i)=const·D,element C in row i and column j R (i, j) = -const·D, element C in row j and column i R(j,i)=-const·D, element C in row j and column j R (j,j)=const·D;

[0049] If the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update the i-th row element Q(i)=C' and the j-th row element Q(j)=-C' in Q;

[0050] If the signal flowing from the electrical system to the control system is a capacitor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C C The element C in row i and column i C (i,i)=D·C, element C in row i and column j C (i, j) = -D·C, the element C in the jth row and ith column C (j,i)=-D·C, element C in row j and column j C (j,j)=D·C;

[0051] If the signal flowing from the electrical system to the control system is an inductor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0052] If the signal flowing from the electrical system to the control system is a resistance current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0053] If the signal flowing from the electrical system to the control system is the node voltage difference of the branch, the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update C R The element C in row i and column i R (i,i)=D, element C in row i and column j R (i,j)=-D, element C in row j and column i R (j,i)=-D, element C in row j and column j R (j,j)=D;

[0054] Where const is the derivative of the internal resistance of the controlled voltage source; C' is the correlation coefficient matrix of the control system state variables and output quantities in the power simulation system; D is the correlation coefficient of the input and output quantities of the control system in the power simulation system; C is the capacitance value of the branch to which the capacitive current signal flowing from the electrical system to the control system belongs; L is the inductance value of the branch to which the inductive current signal flowing from the electrical system to the control system belongs; R is the resistance value of the branch to which the resistive current signal flowing from the electrical system to the control system belongs.

[0055] Preferably, the state equation of the power simulation system is written by using the state equation of the electrical system in the power simulation system, the state equation of the control system and the correlation matrix between the electrical system and the control system in the power simulation system, including:

[0056] Determine the state equation of the power simulation system as follows:

[0057]

[0058] Where K C K is the capacitance coefficient matrix in the electrical system state equation of the power simulation system; R K is the resistance coefficient matrix in the electrical system state equation of the power simulation system; L is the inductance matrix in the electrical system state equation of the power simulation system; C C C is the capacitance and current correlation matrix of the control system to the electrical system; R C is the control system's resistance and current correlation matrix for the electrical system; L P is the correlation matrix of the control system to the inductor current of the electrical system; C P is the correlation matrix of the electrical system capacitance and current to the control system; R P is the electrical system resistance and current control system correlation matrix; L is the correlation matrix of the inductance current of the electrical system to the control system; Q is the correlation matrix of the control state quantity of the control system to the electrical system; E is the unit diagonal matrix; is the potential of node n; w is the state variable matrix of the control system in the power simulation system; A is the state variable coefficient matrix of the transfer function of the control system in the power simulation system; I s is the input vector of the circuit in the electrical system state equation in the power simulation system; B is the input coefficient matrix of the control system in the power simulation system; p is the input quantity of the control system in the power simulation system.

[0059] Furthermore, the electrical system state equation of the power simulation system is determined as follows:

[0060]

[0061] Where K C K is the capacitance coefficient matrix in the electrical system state equation of the power simulation system; R K is the resistance coefficient matrix in the electrical system state equation of the power simulation system; L I is the inductance matrix in the electrical system state equation of the power simulation system; s is the input vector in the circuit of the electrical system state equation of the power simulation system; E is the unit diagonal matrix; Each element value represents the potential of the corresponding node;

[0062] Furthermore, the resistance coefficient matrix K in the electrical system state equation of the power simulation system R , inductance matrix K L , capacitance coefficient matrix K C and the input vector I in the circuit S The acquisition process includes:

[0063] Step a. Initialize the resistance coefficient matrix K in the electrical system state equation of the power simulation system R , inductance matrix K L and the capacitance coefficient matrix K C are all n+1 order zero element matrices; the input vector I in the initialization circuit S is an n+1-order zero-element vector;

[0064] Step b. Update K according to the type of branch in the electrical system R , K L , K C and I S ;

[0065] Step c. Eliminate the resistivity matrix K R , inductance matrix K L and the capacitance coefficient matrix K C The 0th row and column in the elimination circuit; the input vector I S Element 0 in .

[0066] Furthermore, the step b includes:

[0067] If the branch ij in the electrical system is a capacitor branch, then update K C The element K in row i and column i C (i,i)=K C (i,i)+C, element K in row i and column j C (i,j)=K C (i,j)-C, element K in row j and column i C (j,i)=K C(j,i)-C, element K in row j and column j C (j,j)=K C (j,j)+C;

[0068] If the branch ij in the electrical system is an inductive branch, then update K L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0069] If the branch ij in the electrical system is a resistance branch, then update K R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0070] If the branch ij in the electrical system is a voltage source branch, then update K R The element K in row i and column i R (i,i)=K R (i,i)+Const, the element K in row i and column j R (i,j)=K R (i,j)-Const, the element K in the jth row and ith column R (j,i)=K R (j,i)-Const, the element K in the jth row and jth column R (j,j)=K R (j,j)+Const;

[0071] And update the input vector I in the circuit S The i-th vector I s (i)=I s (i)+Const·U S , the jth vector I s (j) = I s (j)-Const·U S ;

[0072] If the branch ij in the electrical system is a current source branch, then the input vector I in the circuit is updated S The i-th vector I s (i)=I s (i)+I' S , the jth vector I s (j) = I s (j)-I' S ;

[0073] Where const is the derivative of the internal resistance of the controlled voltage source in the voltage source branch; R is the resistance value of the resistance branch, L is the inductance value of the inductance branch, and C is the capacitance value of the capacitance branch; U S is the output voltage value of the voltage source in the voltage source branch; I' S is the output current value of the current source in the current source branch.

[0074] The present invention provides a state equation writing system for a power simulation system, wherein the system comprises:

[0075] An acquisition module, used for acquiring a transfer function of a control system in a power simulation system;

[0076] A first determining module is configured to determine a state equation of a control system in a power simulation system according to the transfer function;

[0077] A second determining module is used to determine a correlation matrix between the electrical system and the control system in the power simulation system according to a correlation relationship between the electrical system and the control system in the power simulation system;

[0078] The writing module is used to write the state equation of the power simulation system by using the state equation of the electrical system in the power simulation system, the state equation of the control system and the correlation matrix between the electrical system and the control system in the power simulation system.

[0079] Preferably, the acquisition module is used to:

[0080] The transfer function G(s) of the control system in the power simulation system is determined as follows:

[0081]

[0082] Where, β i =b i -a i b m , b i is the i-order coefficient of the output of the control system in the power simulation system, i∈(1~m-1); a i is the i-order coefficient of the input quantity of the control system in the power simulation system; b m is the m-order coefficient of the output of the control system in the power simulation system; s is the Laplace operator, and m is the number of state variables of the control system.

[0083] Furthermore, the first determining module is configured to:

[0084] Decompose the transfer function of the control system in the power simulation system and obtain the intermediate quantity of the transfer function of the control system in the power simulation system

[0085] Obtaining the Laplace inverse transform z(t) of the intermediate quantity z(s);

[0086] The state variable vector ω=[ω1…ω i …ω m ] T , where ω i is the i-th element in the state variable vector of the control system’s transfer function, ω i =[z(t)] (i-1) ;

[0087] The state equation of the control system in the power simulation system is determined as follows:

[0088]

[0089] Wherein, ω' is the first-order derivative matrix of the state variable matrix of the transfer function of the control system in the power simulation system; A is the state variable coefficient matrix of the transfer function of the control system in the power simulation system; B is the input coefficient matrix of the control system in the power simulation system; C' is the correlation coefficient matrix of the state variables and output quantities of the control system in the power simulation system; D is the correlation coefficient between the input quantity and output quantity of the control system in the power simulation system; p is the input quantity of the control system in the power simulation system; q is the output quantity of the control system in the power simulation system;

[0090] The state variable coefficient matrix A of the transfer function of the control system in the power simulation system is determined as follows:

[0091]

[0092] Where a0 is the constant term coefficient of the control system input in the power simulation system; a1 is the first-order term coefficient of the control system input in the power simulation system; a2 is the second-order term coefficient of the control system input in the power simulation system; a m-1 is the m-1-order coefficient of the input quantity of the control system in the power simulation system;

[0093] The input coefficient matrix B of the control system in the power simulation system is determined as follows:

[0094] B=[0 0 0 … 0 1] T

[0095] Where B is an m-order column vector;

[0096] The correlation coefficient matrix C' between the state variables and output quantities of the control system in the power simulation system is determined as follows:

[0097] C'=[β0 … β i … β m-1 ]

[0098] The correlation coefficient D between the input and output of the control system in the power simulation system is determined as follows:

[0099] D=b m

[0100] Where b m is the m-th order coefficient of the output of the control system in the power simulation system.

[0101] Preferably, the second determining module includes:

[0102] Initialization unit, used to initialize the electrical system capacitance and current to control system correlation matrix P in the electrical system and control system correlation matrix in the power simulation system C , the electrical system inductance current to control system correlation matrix P L , the electrical system resistance and current control system correlation matrix P R are all m×n order zero element matrices, the control system control state quantity to electrical system correlation matrix Q is n×m order zero element matrix, the control system to electrical system capacitance current correlation matrix C C is an n×n order zero element matrix, the control system's inductance current correlation matrix C L is an n×n order zero element matrix and the control system's resistance and current correlation matrix C of the electrical system R is an n×n order zero element matrix, where m is the number of state variables of the control system; n is the number of circuit nodes of the electrical system;

[0103] The first updating unit is used to update P according to the signal of each branch flow control system in the electrical system. C 、P R and P L The mth row element of ;

[0104] The second updating unit is used to update Q and C respectively according to the signal flowing from the control system to the electrical system. C 、C L and C R .

[0105] Furthermore, the first updating unit is configured to:

[0106] If the signal flowing to the control system from branch ij in the electrical system is a capacitor current signal, then update P C The element P in the mth row and ith column C (m,i)=C ij , element P in row m and column j C (m,j)=-Cij ;

[0107] If the signal flowing from branch ij in the electrical system to the control system is the inductor current signal, then update P L The element in row m and column i The element in row m and column j

[0108] If the signal flowing to the control system from branch ij in the electrical system is a resistance current signal, then update P R The element in row m and column i The element in row m and column j

[0109] If the signal flowing from branch ij to the control system in the electrical system is the node voltage difference of branch ij, then update P R The element P in the mth row and ith column R (m,i)=1, element P in row m and column j R (m,j)=-1;

[0110] Among them, R ij is the resistance value of the resistance branch ij, L ij is the inductance value of the inductor branch ij, C ij is the capacitance value of the capacitor branch ij.

[0111] Furthermore, the second updating unit is configured to:

[0112] If the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update the i-th row element Q(i)=const·C' and the j-th row element Q(j)=-const·C' in Q;

[0113] If the signal flowing from the electrical system to the control system is a capacitor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C C The element C in row i and column i C (i,i)=const·D·C, element C in row i and column j C (i, j) = -const·D·C, the element C in the jth row and ith column C (j,i)=-const·D·C, element C in row j and column j C (j,j)=-const·D·C;

[0114] If the signal flowing from the electrical system to the control system is an inductor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0115] If the signal flowing from the electrical system to the control system is a resistance current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0116] If the signal flowing from the electrical system to the control system is the node voltage difference of the branch, the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update C R The element C in row i and column i R (i,i)=const·D,element C in row i and column j R (i, j) = -const·D, element C in row j and column i R (j,i)=-const·D, element C in row j and column j R (j,j)=const·D;

[0117] If the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update the i-th row element Q(i)=C' and the j-th row element Q(j)=-C' in Q;

[0118] If the signal flowing from the electrical system to the control system is a capacitor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C C The element C in row i and column i C (i,i)=D·C, element C in row i and column j C (i, j) = -D·C, the element C in the jth row and ith column C (j,i)=-D·C, element C in row j and column j C (j,j)=D·C;

[0119] If the signal flowing from the electrical system to the control system is an inductor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0120] If the signal flowing from the electrical system to the control system is a resistance current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0121] If the signal flowing from the electrical system to the control system is the node voltage difference of the branch, the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update C R The element C in row i and column i R (i,i)=D, element C in row i and column j R (i,j)=-D, element C in row j and column i R (j,i)=-D, element C in row j and column j R (j,j)=D;

[0122] Where const is the derivative of the internal resistance of the controlled voltage source; C' is the correlation coefficient matrix of the control system state variables and output quantities in the power simulation system; D is the correlation coefficient of the input and output quantities of the control system in the power simulation system; C is the capacitance value of the branch to which the capacitance signal from the electrical system flows to the control system; L is the inductance value of the branch to which the inductance signal from the electrical system flows to the control system; and R is the resistance value of the branch to which the resistance signal from the electrical system flows to the control system.

[0123] Preferably, the column writing module is used to:

[0124] Determine the state equation of the power simulation system as follows:

[0125]

[0126] Where K C K is the capacitance coefficient matrix in the electrical system state equation of the power simulation system; R K is the resistance coefficient matrix in the electrical system state equation of the power simulation system; L is the inductance matrix in the electrical system state equation of the power simulation system; C C C is the capacitance and current correlation matrix of the control system to the electrical system; R C is the control system's resistance and current correlation matrix for the electrical system; L P is the correlation matrix of the control system to the inductor current of the electrical system; C P is the correlation matrix of the electrical system capacitance and current to the control system; RP is the electrical system resistance and current control system correlation matrix; L is the correlation matrix of the inductance current of the electrical system to the control system; Q is the correlation matrix of the control state quantity of the control system to the electrical system; E is the unit diagonal matrix; is the potential of node n; w is the state variable matrix of the control system in the power simulation system; A is the state variable coefficient matrix of the transfer function of the control system in the power simulation system; I s is the input vector of the circuit in the electrical system state equation in the power simulation system; B is the input coefficient matrix of the control system in the power simulation system; p is the input quantity of the control system in the power simulation system.

[0127] Furthermore, the electrical system state equation of the power simulation system is determined as follows:

[0128]

[0129] Where K C K is the capacitance coefficient matrix in the electrical system state equation of the power simulation system; R K is the resistance coefficient matrix in the electrical system state equation of the power simulation system; L I is the inductance matrix in the electrical system state equation of the power simulation system; s is the input vector in the circuit of the electrical system state equation of the power simulation system; E is the unit diagonal matrix; Each element value represents the potential of the corresponding node;

[0130] Furthermore, the resistance coefficient matrix K in the electrical system state equation of the power simulation system R , inductance matrix K L , capacitance coefficient matrix K C and the input vector I in the circuit S The acquisition process includes:

[0131] Step a. Initialize the resistance coefficient matrix K in the electrical system state equation of the power simulation system R , inductance matrix K L and the capacitance coefficient matrix K C are all n+1 order zero element matrices; the input vector I in the initialization circuit S is an n+1-order zero-element vector;

[0132] Step b. Update K according to the type of branch in the electrical system R , K L , K C and I S ;

[0133] Step c. Eliminate the resistivity matrix K R , inductance matrix K L and the capacitance coefficient matrix K C The 0th row and column in the elimination circuit; the input vector I S Element 0 in .

[0134] Furthermore, the step b includes:

[0135] If the branch ij in the electrical system is a capacitor branch, then update K C The element K in row i and column i C (i,i)=K C (i,i)+C, element K in row i and column j C (i,j)=K C (i,j)-C, element K in row j and column i C (j,i)=K C (j,i)-C, element K in row j and column j C (j,j)=K C (j,j)+C;

[0136] If the branch ij in the electrical system is an inductive branch, then update K L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0137] If the branch ij in the electrical system is a resistance branch, then update K R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0138] If the branch ij in the electrical system is a voltage source branch, then update K R The element K in row i and column i R (i,i)=K R (i,i)+Const, the element K in row i and column j R (i,j)=K R (i,j)-Const, the element K in the jth row and ith column R (j,i)=K R (j,i)-Const, the element K in the jth row and jth column R (j,j)=K R (j,j)+Const;

[0139] And update the input vector I in the circuit S The i-th vector I s (i)=I s (i)+Const·U S , the jth vector I s (j) = I s (j)-Const·U S ;

[0140] If the branch ij in the electrical system is a current source branch, then the input vector I in the circuit is updated S The i-th vector I s (i)=I s (i)+I' S , the jth vector I s (j) = I s (j)-I' S ;

[0141] Where const is the derivative of the internal resistance of the controlled voltage source in the voltage source branch; R is the resistance value of the resistance branch, L is the inductance value of the inductance branch, and C is the capacitance value of the capacitance branch; U S is the output voltage value of the voltage source in the voltage source branch; I' S is the output current value of the current source in the current source branch.

[0142] Compared with the closest prior art, the present invention has the following beneficial effects:

[0143] The technical solution provided by the present invention obtains the transfer function of the control system in the power simulation system; determines the state equation of the control system in the power simulation system based on the transfer function; determines the correlation matrix between the electrical system and the control system in the power simulation system based on the correlation relationship between the electrical system and the control system in the power simulation system; and writes the state equation of the power simulation system using the circuit network state equation of the electrical system in the power simulation system, the state equation of the control system, and the correlation matrix between the electrical system and the control system in the power simulation system. The state equation writing method of the technical solution provided by the present invention is simple and suitable for computer solution. When applied to electromagnetic transient simulation, the simulation results are obtained by solving the state equation of the power simulation system, thereby eliminating the time delay caused by the electromagnetic transient simulation gradually solving the electrical system and the control system, thereby improving the numerical stability and calculation accuracy of the electromagnetic transient simulation. BRIEF DESCRIPTION OF THE DRAWINGS

[0144] Figure 1 It is a schematic diagram of the solution structure of traditional electrical system and control system;

[0145] Figure 2 It is a flow chart of a method for writing state equations of an electric power simulation system;

[0146] Figure 3 Schematic diagram of branch ij in an embodiment of the present invention;

[0147] Figure 4 It is a system framework diagram for writing state equations of an electric power simulation system. DETAILED DESCRIPTION

[0148] The specific embodiments of the present invention will be further described in detail below with reference to the accompanying drawings.

[0149] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts shall fall within the scope of protection of the present invention.

[0150] The present invention provides a method for writing state equations of a power simulation system, such as Figure 2 As shown, the method includes:

[0151] Step 101. Obtain the transfer function of the control system in the power simulation system;

[0152] Step 102: Determine the state equation of the control system in the power simulation system according to the transfer function;

[0153] Step 103. Determine a correlation matrix between the electrical system and the control system in the power simulation system according to the correlation relationship between the electrical system and the control system in the power simulation system;

[0154] Step 104 . Write the state equation of the power simulation system using the state equation of the electrical system in the power simulation system, the state equation of the control system, and the correlation matrix between the electrical system and the control system in the power simulation system.

[0155] In the most preferred embodiment of the present invention, Figure 3 As shown, the branch ij may be any one of a capacitance branch, an inductance branch, a resistance branch, a voltage source branch and a current source branch in the electrical system;

[0156] Specifically, step 101 includes:

[0157] The transfer function G(s) of the control system in the power simulation system is determined as follows:

[0158]

[0159] Where, β i =b i-a i b m , b i is the i-order coefficient of the output of the control system in the power simulation system, i∈(1~m-1); a i is the i-order coefficient of the input quantity of the control system in the power simulation system; b m is the m-order coefficient of the output of the control system in the power simulation system; s is the Laplace operator, and m is the number of state variables of the control system.

[0160] Furthermore, the step 102 includes:

[0161] Decompose the transfer function of the control system in the power simulation system and obtain the intermediate quantity of the transfer function of the control system in the power simulation system

[0162] Obtaining the Laplace inverse transform z(t) of the intermediate quantity z(s);

[0163] The state variable matrix ω=[ω1…ω i …ω m ] T , where ω i is the i-th element in the state variable matrix of the control system’s transfer function, ω i =[z(t)] (i-1) ;

[0164] The state equation of the control system in the power simulation system is determined as follows:

[0165]

[0166] Wherein, ω' is the first-order derivative matrix of the state variable matrix of the transfer function of the control system in the power simulation system; A is the state variable coefficient matrix of the transfer function of the control system in the power simulation system; B is the input coefficient matrix of the control system in the power simulation system; C' is the correlation coefficient matrix of the state variables and output quantities of the control system in the power simulation system; D is the correlation coefficient between the input quantity and output quantity of the control system in the power simulation system; p is the input quantity of the control system in the power simulation system; q is the output quantity of the control system in the power simulation system;

[0167] The state variable coefficient matrix A of the transfer function of the control system in the power simulation system is determined as follows:

[0168]

[0169] Where a0 is the constant term coefficient of the control system input in the power simulation system; a1 is the first-order term coefficient of the control system input in the power simulation system; a2 is the second-order term coefficient of the control system input in the power simulation system; a m-1 is the m-1-order coefficient of the input quantity of the control system in the power simulation system;

[0170] The input coefficient matrix B of the control system in the power simulation system is determined as follows:

[0171] B=[0 0 0 … 0 1] T

[0172] Where B is an m-order column vector;

[0173] The correlation coefficient matrix C' between the state variables and output quantities of the control system in the power simulation system is determined as follows:

[0174] C'=[β0 … β i … β m-1 ]

[0175] The correlation coefficient D between the input and output of the control system in the power simulation system is determined as follows:

[0176] D=b m

[0177] Where b m is the m-th order coefficient of the output of the control system in the power simulation system.

[0178] Specifically, step 103 includes:

[0179] Step a. Initialize the electrical system capacitance and current to control system correlation matrix P in the electrical system and control system correlation matrix in the power simulation system. C , the electrical system inductance current to control system correlation matrix P L , the electrical system resistance and current control system correlation matrix P R are all m×n order zero element matrices, the control system control state quantity to electrical system correlation matrix Q is n×m order zero element matrix, the control system to electrical system capacitance current correlation matrix C C is an n×n order zero element matrix, the control system's inductance current correlation matrix C L is an n×n order zero element matrix and the control system's resistance and current correlation matrix C of the electrical system R is an n×n order zero element matrix, where m is the number of state variables of the control system; n is the number of circuit nodes of the electrical system;

[0180] Step b. Update P according to the signal of each branch flow control system in the electrical systemC 、P R and P L The mth row element of ;

[0181] Step c. Update Q and C respectively according to the signal flowing from the control system to the electrical system C 、C L and C R .

[0182] Furthermore, the step b includes:

[0183] If the signal flowing to the control system from branch ij in the electrical system is a capacitor current signal, then update P C The element P in the mth row and ith column C (m,i)=C ij , element P in row m and column j C (m,j)=-C ij ;

[0184] If the signal flowing from branch ij in the electrical system to the control system is the inductor current signal, then update P L The element in row m and column i The element in row m and column j

[0185] If the signal flowing to the control system from branch ij in the electrical system is a resistance current signal, then update P R The element in row m and column i The element in row m and column j

[0186] If the signal flowing from branch ij to the control system in the electrical system is the node voltage difference of branch ij, then update P R The element P in the mth row and ith column R (m,i)=1, element P in row m and column j R (m,j)=-1;

[0187] Among them, R ij is the resistance value of the resistance branch ij, L ij is the inductance value of the inductor branch ij, C ij is the capacitance value of the capacitor branch ij.

[0188] Furthermore, the step c comprises:

[0189] If the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update the i-th row element Q(i)=const·C' and the j-th row element Q(j)=-const·C' in Q;

[0190] If the signal flowing from the electrical system to the control system is a capacitor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C C The element C in row i and column i C (i,i)=const·D·C, element C in row i and column j C (i, j) = -const·D·C, the element C in the jth row and ith column C (j,i)=-const·D·C, element C in row j and column j C (j,j)=-const·D·C;

[0191] If the signal flowing from the electrical system to the control system is an inductor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0192] If the signal flowing from the electrical system to the control system is a resistance current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0193] If the signal flowing from the electrical system to the control system is the node voltage difference of the branch, the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update C R The element C in row i and column i R (i,i)=const·D,element C in row i and column j R (i, j) = -const·D, element C in row j and column i R (j,i)=-const·D, element C in row j and column j R (j,j)=const·D;

[0194] If the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update the i-th row element Q(i)=C' and the j-th row element Q(j)=-C' in Q;

[0195] If the signal flowing from the electrical system to the control system is a capacitor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C C The element C in row i and column i C (i,i)=D·C, element C in row i and column j C (i, j) = -D·C, the element C in the jth row and ith column C (j,i)=-D·C, element C in row j and column j C (j,j)=D·C;

[0196] If the signal flowing from the electrical system to the control system is an inductor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0197] If the signal flowing from the electrical system to the control system is a resistance current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0198] If the signal flowing from the electrical system to the control system is the node voltage difference of the branch, the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update C R The element C in row i and column i R (i,i)=D, element C in row i and column j R (i,j)=-D, element C in row j and column i R (j,i)=-D, element C in row j and column j R (j,j)=D;

[0199] Where const is the derivative of the internal resistance of the controlled voltage source; C' is the correlation coefficient matrix of the control system state variables and output quantities in the power simulation system; D is the correlation coefficient of the input and output quantities of the control system in the power simulation system; C is the capacitance value of the branch to which the capacitance signal from the electrical system flows to the control system; L is the inductance value of the branch to which the inductance signal from the electrical system flows to the control system; and R is the resistance value of the branch to which the resistance signal from the electrical system flows to the control system.

[0200] Specifically, step 104 includes:

[0201] Determine the state equation of the power simulation system as follows:

[0202]

[0203] Where K C K is the capacitance coefficient matrix in the electrical system state equation of the power simulation system; R K is the resistance coefficient matrix in the electrical system state equation of the power simulation system; L is the inductance matrix in the electrical system state equation of the power simulation system; C C C is the capacitance and current correlation matrix of the control system to the electrical system; R C is the control system's resistance and current correlation matrix for the electrical system; L P is the correlation matrix of the control system to the inductor current of the electrical system; C P is the correlation matrix of the electrical system capacitance and current to the control system; R P is the electrical system resistance and current control system correlation matrix; L is the correlation matrix of the inductance current of the electrical system to the control system; Q is the correlation matrix of the control state quantity of the control system to the electrical system; E is the unit diagonal matrix; is the potential of node n; w is the state variable matrix of the control system in the power simulation system; A is the state variable coefficient matrix of the transfer function of the control system in the power simulation system; I s is the input vector of the circuit in the electrical system state equation in the power simulation system; B is the input coefficient matrix of the control system in the power simulation system; p is the input quantity of the control system in the power simulation system.

[0204] Furthermore, the electrical system state equation of the power simulation system is determined as follows:

[0205]

[0206] Where K C K is the capacitance coefficient matrix in the electrical system state equation of the power simulation system; R K is the resistance coefficient matrix in the electrical system state equation of the power simulation system; L I is the inductance matrix in the electrical system state equation of the power simulation system; s is the input vector in the circuit of the electrical system state equation of the power simulation system; E is the unit diagonal matrix; Each element value represents the potential of the corresponding node;

[0207] Furthermore, the resistance coefficient matrix K in the electrical system state equation of the power simulation system R , inductance matrix K L , capacitance coefficient matrix K C and the input vector I in the circuit S The acquisition process includes:

[0208] Step a. Initialize the resistance coefficient matrix K in the electrical system state equation of the power simulation system R , inductance matrix K L and the capacitance coefficient matrix K C are all n+1 order zero element matrices; the input vector I in the initialization circuit S is an n+1-order zero-element vector;

[0209] Step b. Update K according to the type of branch in the electrical system R , K L , K C and I S ;

[0210] Step c. Eliminate the resistivity matrix K R , inductance matrix K L and the capacitance coefficient matrix K C The 0th row and column in the elimination circuit; the input vector I S Element 0 in .

[0211] Furthermore, the step b includes:

[0212] If the branch ij in the electrical system is a capacitor branch, then update K C The element K in row i and column i C (i,i)=K C (i,i)+C, element K in row i and column j C (i,j)=K C (i,j)-C, element K in row j and column i C (j,i)=K C (j,i)-C, element K in row j and column j C (j,j)=K C (j,j)+C;

[0213] If the branch ij in the electrical system is an inductive branch, then update K L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0214] If the branch ij in the electrical system is a resistance branch, then update K R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0215] If the branch ij in the electrical system is a voltage source branch, then update K R The element K in row i and column i R (i,i)=K R (i,i)+Const, the element K in row i and column j R (i,j)=K R (i,j)-Const, the element K in the jth row and ith column R (j,i)=K R (j,i)-Const, the element K in the jth row and jth column R (j,j)=K R (j,j)+Const;

[0216] And update the input vector I in the circuit S The i-th vector I s (i)=I s (i)+Const·U S , the jth vector I s (j) = I s (j)-Const·U S ;

[0217] If the branch ij in the electrical system is a current source branch, then the input vector I in the circuit is updated S The i-th vector I s (i)=I s (i)+I' S , the jth vector I s (j) = I s (j)-I' S ;

[0218] Where const is the derivative of the internal resistance of the controlled voltage source in the voltage source branch; R is the resistance value of the resistance branch, L is the inductance value of the inductance branch, and C is the capacitance value of the capacitance branch; U S is the output voltage value of the voltage source in the voltage source branch; I' S is the output current value of the current source in the current source branch.

[0219] The present invention provides a state equation writing system for a power simulation system, such as Figure 4 As shown, the system includes:

[0220] An acquisition module, used for acquiring a transfer function of a control system in a power simulation system;

[0221] A first determining module is configured to determine a state equation of a control system in a power simulation system according to the transfer function;

[0222] A second determining module is used to determine a correlation matrix between the electrical system and the control system in the power simulation system according to a correlation relationship between the electrical system and the control system in the power simulation system;

[0223] The writing module is used to write the state equation of the power simulation system by using the state equation of the electrical system in the power simulation system, the state equation of the control system and the correlation matrix between the electrical system and the control system in the power simulation system.

[0224] Specifically, the acquisition module is used to:

[0225] The transfer function G(s) of the control system in the power simulation system is determined as follows:

[0226]

[0227] Where, β i =b i -a i b m , b i is the i-order coefficient of the output of the control system in the power simulation system, i∈(1~m-1); a i is the i-order coefficient of the input quantity of the control system in the power simulation system; b m is the m-order coefficient of the output of the control system in the power simulation system; s is the Laplace operator, and m is the number of state variables of the control system.

[0228] Specifically, the first determining module is used to:

[0229] Decompose the transfer function of the control system in the power simulation system and obtain the intermediate quantity of the transfer function of the control system in the power simulation system

[0230] Obtaining the Laplace inverse transform z(t) of the intermediate quantity z(s);

[0231] The state variable matrix ω=[ω1…ω i …ω m ] T , where ω i is the i-th element in the state variable matrix of the control system’s transfer function, ω i =[z(t)] (i-1) ;

[0232] The circuit network state equation of the control system in the power simulation system is determined as follows:

[0233]

[0234] Wherein, ω' is the first-order derivative matrix of the state variable matrix of the transfer function of the control system in the power simulation system; A is the state variable coefficient matrix of the transfer function of the control system in the power simulation system; B is the input coefficient matrix of the control system in the power simulation system; C' is the correlation coefficient matrix of the state variables and output quantities of the control system in the power simulation system; D is the correlation coefficient between the input quantity and output quantity of the control system in the power simulation system; p is the input quantity of the control system in the power simulation system; q is the output quantity of the control system in the power simulation system;

[0235] The state variable coefficient matrix A of the transfer function of the control system in the power simulation system is determined as follows:

[0236]

[0237] Where a0 is the constant term coefficient of the control system input in the power simulation system; a1 is the first-order term coefficient of the control system input in the power simulation system; a2 is the second-order term coefficient of the control system input in the power simulation system; a m-1 is the m-1-order coefficient of the input quantity of the control system in the power simulation system;

[0238] The input coefficient matrix B of the control system in the power simulation system is determined as follows:

[0239] B=[0 0 0 … 0 1] T

[0240] Where B is an m-order column vector;

[0241] The correlation coefficient matrix C' between the state variables and output quantities of the control system in the power simulation system is determined as follows:

[0242] C'=[β0 … β i … β m-1 ]

[0243] The correlation coefficient D between the input and output of the control system in the power simulation system is determined as follows:

[0244] D=b m

[0245] Where b m is the m-th order coefficient of the output of the control system in the power simulation system.

[0246] Specifically, the second determining module includes:

[0247] Initialization unit, used to initialize the electrical system capacitance and current to control system correlation matrix P in the electrical system and control system correlation matrix in the power simulation system C , the electrical system inductance current to control system correlation matrix P L , the electrical system resistance and current control system correlation matrix P R are all m×n order zero element matrices, the control system control state quantity to electrical system correlation matrix Q is n×m order zero element matrix, the control system to electrical system capacitance current correlation matrix C C is an n×n order zero element matrix, the control system's inductance current correlation matrix C L is an n×n order zero element matrix and the control system's resistance and current correlation matrix C of the electrical system R is an n×n order zero element matrix, where m is the number of state variables of the control system; n is the number of circuit nodes of the electrical system;

[0248] The first updating unit is used to update P according to the signal of each branch flow control system in the electrical system. C 、P R and P L The mth row element of ;

[0249] The second updating unit is used to update Q and C respectively according to the signal flowing from the control system to the electrical system. C 、C L and C R .

[0250] Specifically, the first updating unit is used to:

[0251] If the signal flowing to the control system from branch ij in the electrical system is a capacitor current signal, then update P C The element P in the mth row and ith column C (m,i)=C ij , element P in row m and column j C (m,j)=-C ij ;

[0252] If the signal flowing from branch ij in the electrical system to the control system is the inductor current signal, then update P L The element in row m and column i The element in row m and column j

[0253] If the signal flowing to the control system from branch ij in the electrical system is a resistance current signal, then update P R The element in row m and column i The element in row m and column j

[0254] If the signal flowing from branch ij to the control system in the electrical system is the node voltage difference of branch ij, then update P R The element P in the mth row and ith column R (m,i)=1, element P in row m and column j R (m,j)=-1;

[0255] Among them, R ij is the resistance value of the resistance branch ij, L ij is the inductance value of the inductor branch ij, C ij is the capacitance value of the capacitor branch ij.

[0256] Specifically, the second updating unit is used to:

[0257] If the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update the i-th row element Q(i)=const·C' and the j-th row element Q(j)=-const·C' in Q;

[0258] If the signal flowing from the electrical system to the control system is a capacitor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C C The element C in row i and column i C (i,i)=const·D·C, element C in row i and column j C (i, j) = -const·D·C, the element C in the jth row and ith column C (j,i)=-const·D·C, element C in row j and column j C (j,j)=-const·D·C;

[0259] If the signal flowing from the electrical system to the control system is an inductor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0260] If the signal flowing from the electrical system to the control system is a resistance current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0261] If the signal flowing from the electrical system to the control system is the node voltage difference of the branch, the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update C R The element C in row i and column i R (i,i)=const·D,element C in row i and column j R (i, j) = -const·D, element C in row j and column i R (j,i)=-const·D, element C in row j and column j R (j,j)=const·D;

[0262] If the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update the i-th row element Q(i)=C' and the j-th row element Q(j)=-C' in Q;

[0263] If the signal flowing from the electrical system to the control system is a capacitor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C C The element C in row i and column i C (i,i)=D·C, element C in row i and column j C (i, j) = -D·C, the element C in the jth row and ith column C (j,i)=-D·C, element C in row j and column j C (j,j)=D·C;

[0264] If the signal flowing from the electrical system to the control system is an inductor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0265] If the signal flowing from the electrical system to the control system is a resistance current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0266] If the signal flowing from the electrical system to the control system is the node voltage difference of the branch, the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update C R The element C in row i and column i R (i,i)=D, element C in row i and column j R (i,j)=-D, element C in row j and column i R (j,i)=-D, element C in row j and column j R (j,j)=D;

[0267] Where const is the derivative of the internal resistance of the controlled voltage source; C' is the correlation coefficient matrix of the control system state variables and output quantities in the power simulation system; D is the correlation coefficient of the input and output quantities of the control system in the power simulation system; C is the capacitance value of the branch to which the capacitance signal from the electrical system flows to the control system; L is the inductance value of the branch to which the inductance signal from the electrical system flows to the control system; and R is the resistance value of the branch to which the resistance signal from the electrical system flows to the control system.

[0268] Specifically, the column writing module is used to:

[0269] Determine the state equation of the power simulation system as follows:

[0270]

[0271] Where K C K is the capacitance coefficient matrix in the electrical system state equation of the power simulation system; R K is the resistance coefficient matrix in the electrical system state equation of the power simulation system; L is the inductance matrix in the electrical system state equation of the power simulation system; C C C is the capacitance and current correlation matrix of the control system to the electrical system; R C is the control system's resistance and current correlation matrix for the electrical system; L P is the correlation matrix of the control system to the inductor current of the electrical system; C P is the correlation matrix of the electrical system capacitance and current to the control system; R P is the electrical system resistance and current control system correlation matrix; L is the correlation matrix of the inductance current of the electrical system to the control system; Q is the correlation matrix of the control state quantity of the control system to the electrical system; E is the unit diagonal matrix; is the potential of node n; w is the state variable matrix of the control system in the power simulation system; A is the state variable coefficient matrix of the transfer function of the control system in the power simulation system; I sis the input vector of the circuit in the electrical system state equation in the power simulation system; B is the input coefficient matrix of the control system in the power simulation system; p is the input quantity of the control system in the power simulation system.

[0272] Specifically, the electrical system state equation of the power simulation system is determined as follows:

[0273]

[0274] Where K C K is the capacitance coefficient matrix in the electrical system state equation of the power simulation system; R K is the resistance coefficient matrix in the electrical system state equation of the power simulation system; L I is the inductance matrix in the electrical system state equation of the power simulation system; s is the input vector in the circuit of the electrical system state equation of the power simulation system; E is the unit diagonal matrix; Each element value represents the potential of the corresponding node;

[0275] Specifically, the resistance coefficient matrix K in the electrical system state equation of the power simulation system R , inductance matrix K L , capacitance coefficient matrix K C and the input vector I in the circuit S The acquisition process includes:

[0276] Step a. Initialize the resistance coefficient matrix K in the electrical system state equation of the power simulation system R , inductance matrix K L and the capacitance coefficient matrix K C are all n+1 order zero element matrices; the input vector I in the initialization circuit S is an n+1-order zero-element vector;

[0277] Step b. Update K according to the type of branch in the electrical system R , K L , K C and I S ;

[0278] Step c. Eliminate the resistivity matrix K R , inductance matrix K L and the capacitance coefficient matrix K C The 0th row and column in the elimination circuit; the input vector I S Element 0 in .

[0279] Specifically, the step b includes:

[0280] If the branch ij in the electrical system is a capacitor branch, then update K C The element K in row i and column i C (i,i)=K C (i,i)+C, element K in row i and column j C (i,j)=K C (i,j)-C, element K in row j and column i C (j,i)=K C (j,i)-C, element K in row j and column j C (j,j)=K C (j,j)+C;

[0281] If the branch ij in the electrical system is an inductive branch, then update K L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0282] If the branch ij in the electrical system is a resistance branch, then update K R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j

[0283] If the branch ij in the electrical system is a voltage source branch, then update K R The element K in row i and column i R (i,i)=K R (i,i)+Const, the element K in row i and column j R (i,j)=K R (i,j)-Const, the element K in the jth row and ith column R (j,i)=K R (j,i)-Const, the element K in the jth row and jth column R (j,j)=K R (j,j)+Const;

[0284] And update the input vector I in the circuit S The i-th vector I s (i)=I s (i)+Const·U S , the jth vector I s (j) = I s (j)-Const·U S ;

[0285] If the branch ij in the electrical system is a current source branch, then the input vector I in the circuit is updated S The i-th vector I s (i)=I s (i)+I' S , the jth vector I s (j) = I s (j)-I' S ;

[0286] Where const is the derivative of the internal resistance of the controlled voltage source in the voltage source branch; R is the resistance value of the resistance branch, L is the inductance value of the inductance branch, and C is the capacitance value of the capacitance branch; U S is the output voltage value of the voltage source in the voltage source branch; I' S is the output current value of the current source in the current source branch.

[0287] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0288] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the steps in the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0289] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0290] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0291] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the field should understand that the specific implementation methods of the present invention can still be modified or replaced by equivalents. Any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be covered by the scope of protection of the claims of the present invention.

Claims

1. A method for writing state equations of a power simulation system, characterized in that: The method comprises: Obtain the transfer function of the control system in the power simulation system; determining a state equation of a control system in a power simulation system according to the transfer function; Determine the correlation matrix between the electrical system and the control system in the power simulation system according to the correlation relationship between the electrical system and the control system in the power simulation system; Write the state equation of the power simulation system using the state equation of the electrical system in the power simulation system, the circuit network state equation of the control system, and the correlation matrix between the electrical system and the control system in the power simulation system; The determining of the association matrix between the electrical system and the control system in the power simulation system according to the association relationship between the electrical system and the control system in the power simulation system includes: Initialize the electrical system and control system correlation matrix P in the electrical system capacitance and current to control system correlation matrix in the power simulation system C , the electrical system inductance current to control system correlation matrix P L , the electrical system resistance and current control system correlation matrix P R are all m×n order zero element matrices, the control system control state quantity to electrical system correlation matrix Q is n×m order zero element matrix, the control system to electrical system capacitance current correlation matrix C C is an n×n order zero element matrix, the control system's inductance current correlation matrix C L is an n×n order zero element matrix and the control system's resistance and current correlation matrix C of the electrical system R is an n×n order zero element matrix, where m is the number of state variables of the control system; n is the number of circuit nodes of the electrical system; According to the signals of each branch flow control system in the electrical system, P is updated respectively. C 、P R and P L The mth row element of ; Update Q and C respectively according to the signal flowing from the control system to the electrical system C 、C L and C R ; The state equation of the power simulation system is written by using the state equation of the electrical system in the power simulation system, the circuit network state equation of the control system, and the correlation matrix between the electrical system and the control system in the power simulation system, including: The circuit network state equation of the power simulation system is determined as follows: Where K C K is the capacitance coefficient matrix in the electrical system state equation of the power simulation system; R K is the resistance coefficient matrix in the electrical system state equation of the power simulation system; L is the inductance matrix in the electrical system state equation of the power simulation system; C C C is the capacitance and current correlation matrix of the control system to the electrical system; R C is the control system's resistance and current correlation matrix for the electrical system; L P is the correlation matrix of the control system to the inductor current of the electrical system; C P is the correlation matrix of the electrical system capacitance and current to the control system; R P is the electrical system resistance and current control system correlation matrix; L is the correlation matrix of the inductance current of the electrical system to the control system; Q is the correlation matrix of the control state quantity of the control system to the electrical system; E is the unit diagonal matrix; is the potential of node n; ω is the state variable vector of the control system in the power simulation system; A is the state variable coefficient matrix of the transfer function of the control system in the power simulation system; I s is the input vector of the circuit in the state equation of the electrical system in the power simulation system; B is the input coefficient matrix of the control system in the power simulation system; p is the input quantity of the control system in the power simulation system; The electrical system state equation of the power simulation system is determined as follows: Where K C K is the capacitance coefficient matrix in the electrical system state equation of the power simulation system; R K is the resistance coefficient matrix in the electrical system state equation of the power simulation system; L I is the inductance matrix in the electrical system state equation of the power simulation system; s is the input vector in the circuit of the electrical system state equation of the power simulation system; E is the unit diagonal matrix; Each element value represents the potential of the corresponding node; The state equation of the control system in the power simulation system is determined as follows: Wherein, ω' is the first-order derivative matrix of the state variable vector of the transfer function of the control system in the power simulation system; A is the state variable coefficient matrix of the transfer function of the control system in the power simulation system; B is the input coefficient matrix of the control system in the power simulation system; C' is the correlation coefficient matrix of the state variables and output quantities of the control system in the power simulation system; D is the correlation coefficient between the input and output quantities of the control system in the power simulation system; p is the input quantity of the control system in the power simulation system; q is the output quantity of the control system in the power simulation system.

2. The method according to claim 1, wherein The obtaining of the transfer function of the control system in the power simulation system includes: The transfer function G(s) of the control system in the power simulation system is determined as follows: Where, β i =b i -a i b m , b i is the i-order coefficient of the output of the control system in the power simulation system, i∈(1~m-1); a i is the i-order coefficient of the input quantity of the control system in the power simulation system; b m is the m-order coefficient of the output of the control system in the power simulation system; s is the Laplace operator, and m is the number of state variables of the control system.

3. The method according to claim 2, wherein Determining the state equation of the control system in the power simulation system according to the transfer function includes: Decompose the transfer function of the control system in the power simulation system and obtain the intermediate quantity of the transfer function of the control system in the power simulation system Obtaining the Laplace inverse transform z(t) of the intermediate quantity z(s); The state variable vector ω=[ω1…ω i …ω m ] T , where ω i is the i-th element in the state variable matrix of the control system’s transfer function, ω i =[z(t)] (i-1) ; The state variable coefficient matrix A of the transfer function of the control system in the power simulation system is determined as follows: Where a0 is the constant term coefficient of the control system input in the power simulation system; a1 is the first-order term coefficient of the control system input in the power simulation system; a2 is the second-order term coefficient of the control system input in the power simulation system; a m-1 is the m-1-order coefficient of the input quantity of the control system in the power simulation system; The input coefficient matrix B of the control system in the power simulation system is determined as follows: B＝[0 0 0…0 1] T Where B is an m-order column vector; The correlation coefficient matrix C' between the state variables and output quantities of the control system in the power simulation system is determined as follows: C'＝[β0 … β i … β m-1 ] The correlation coefficient D between the input and output of the control system in the power simulation system is determined as follows: D=b m Where b m is the m-th order coefficient of the output of the control system in the power simulation system.

4. The method according to claim 1, wherein The P is updated according to the signal of each branch flow control system in the electrical system. C 、P R and P L The τth row elements of include: If the signal flowing to the control system from branch ij in the electrical system is a capacitor current signal, then update P C The element P in the τth row and the i-th column C (τ, i) = C ij , the element P in the τth row and jth column C (τ, i) = -C ij If the signal flowing from branch ij in the electrical system to the control system is the inductor current signal, then update P L The element in row τ and column i The element in row τ and column j If the signal flowing to the control system from branch ij in the electrical system is a resistance current signal, then update P R The element in row τ and column i The element in row τ and column j If the signal flowing from branch ij to the control system in the electrical system is the node voltage difference of branch ij, then update P R The element P in the τth row and the i-th column R (τ, i) = 1, the element P in the τth row and jth column R (τ, j) = -1; Among them, R ij is the resistance value of the resistance branch ij, L ij is the inductance value of the inductor branch ij, C ij is the capacitance value of the capacitor branch ij.

5. The method according to claim 1, wherein The signals flowing from the control system to the electrical system are updated respectively for Q and C C 、C L and C R ,include: If the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update the i-th row element Q(i)=const·C' and the j-th row element Q(j)=-const·C' in Q; If the signal flowing from the electrical system to the control system is a capacitor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C C The element C in row i and column i C (i,i)=const·D·C, element C in row i and column j C (i, j) = -const·D·C, the element C in the jth row and ith column C (j,i)=-const·D·C, element C in row j and column j C (j,j)=-const·D·C; If the signal flowing from the electrical system to the control system is an inductor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j If the signal flowing from the electrical system to the control system is a resistance current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j If the signal flowing from the electrical system to the control system is the node voltage difference of the branch, the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update C R The element C in row i and column i R (i,i)=const·D,element C in row i and column j R (i, j) = -const·D, element C in row j and column i R (j,i)=-const·D, element C in row j and column j R (j,j)=const·D; If the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update the i-th row element Q(i)=C' and the j-th row element Q(j)=-C' in Q; If the signal flowing from the electrical system to the control system is a capacitor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C C The element C in row i and column i C (i,i)=D·C, element C in row i and column j C (i, j) = -D·C, the element C in the jth row and ith column C (j,i)=-D·C, element C in row j and column j C (j,j)=D·C; If the signal flowing from the electrical system to the control system is an inductor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j If the signal flowing from the electrical system to the control system is a resistance current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j If the signal flowing from the electrical system to the control system is the node voltage difference of the branch, the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update C R The element C in row i and column i R (i,i)=D, element C in row i and column j R (i,j)=-D, element C in row j and column i R (j,i)=-D, element C in row j and column j R (j,j)=D; Where const is the derivative of the internal resistance of the controlled voltage source in the voltage source branch; C' is the correlation coefficient matrix between the state variables and output quantities of the control system in the power simulation system; D is the correlation coefficient between the input and output quantities of the control system in the power simulation system; C is the capacitance value of the capacitor branch; L is the inductance value of the inductor branch; and R is the resistance value of the resistor branch.

6. The method according to claim 1, wherein The resistance coefficient matrix K in the electrical system state equation of the power simulation system R , inductance matrix K L , capacitance coefficient matrix K C and the input vector I in the circuit S The acquisition process includes: Step a. Initialize the resistance coefficient matrix K in the electrical system state equation of the power simulation system R , inductance matrix K L and the capacitance coefficient matrix K C are all n+1 order zero element matrices; the input vector I in the initialization circuit S is an n+1-order zero-element vector; Step b. Update K according to the type of branch in the electrical system R , K L , K C and I S ; Step c. Eliminate the resistivity matrix K R , inductance matrix K L and the capacitance coefficient matrix K C The 0th row and column in the elimination circuit; the input vector I S Element 0 in .

7. The method according to claim 6, wherein The step b comprises: If the branch ij in the electrical system is a capacitor branch, then update K C The element K in row i and column i C (i,i)=K C (i,i)+C, element K in row i and column j C (i,j)=K C (i,j)-C, element K in row j and column i C (j,i)=K C (j,i)-C, element K in row j and column j C (j,j)=K C (j,j)+C; If the branch ij in the electrical system is an inductive branch, then update K L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j If the branch ij in the electrical system is a resistance branch, then update K R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j If the branch ij in the electrical system is a voltage source branch, then update K R The element K in row i and column i R (i,i)=K R (i,i)+Const, the element K in row i and column j R (i,j)=K R (i,j)-Const, the element K in the jth row and ith column R (j,i)=K R (j,i)-Const, the element K in the jth row and jth column R (j,j)=K R (j,j)+Const; And update the input vector I in the circuit S The i-th vector I s (i)=I s (i)+Const·U S , the jth vector I s (j) = I s (j)-Const·U S ; If the branch ij in the electrical system is a current source branch, then the input vector I in the circuit is updated S The i-th vector I s (i)=I s (i)+I' S , the jth vector I s (j) = I s (j)-I' S ; Where const is the derivative of the internal resistance of the controlled voltage source in the voltage source branch; R is the resistance value of the resistance branch, L is the inductance value of the inductance branch, and C is the capacitance value of the capacitance branch; U S is the output voltage value of the voltage source in the voltage source branch; I' S is the output current value of the current source in the current source branch.

8. A system for writing circuit network state equations for a power simulation system, characterized in that: The system comprises: An acquisition module, used for acquiring a transfer function of a control system in a power simulation system; A first determining module is used to determine a state equation of a control system in the power simulation system according to the transfer function; A second determining module is used to determine a correlation matrix between the electrical system and the control system in the power simulation system according to a correlation relationship between the electrical system and the control system in the power simulation system; A writing module is used to write the state equation of the power simulation system by using the state equation of the electrical system in the power simulation system, the circuit network state equation of the control system and the correlation matrix between the electrical system and the control system in the power simulation system; The second determining module includes: Initialization unit, used to initialize the electrical system capacitance and current to control system correlation matrix P in the electrical system and control system correlation matrix in the power simulation system C , the electrical system inductance current to control system correlation matrix P L , the electrical system resistance and current control system correlation matrix P R are all m×n order zero element matrices, the control system control state quantity to electrical system correlation matrix Q is n×m order zero element matrix, the control system to electrical system capacitance current correlation matrix C C is an n×n order zero element matrix, the control system's inductance current correlation matrix C L is an n×n order zero element matrix and the control system's resistance and current correlation matrix C of the electrical system R is an n×n order zero element matrix, where m is the number of state variables of the control system; n is the number of circuit nodes of the electrical system; The first updating unit is used to update P according to the signal of each branch flow control system in the electrical system. C 、P R and P L The mth row element of ; The second updating unit is used to update Q and C respectively according to the signal flowing from the control system to the electrical system. C 、C L and C R ; The column writing module is used to: The circuit network state equation of the power simulation system is determined as follows: Where K C K is the capacitance coefficient matrix in the electrical system state equation of the power simulation system; R K is the resistance coefficient matrix in the electrical system state equation of the power simulation system; L is the inductance matrix in the electrical system state equation of the power simulation system; C C C is the capacitance and current correlation matrix of the control system to the electrical system; R C is the control system's resistance and current correlation matrix for the electrical system; L P is the correlation matrix of the control system to the inductor current of the electrical system; C P is the correlation matrix of the electrical system capacitance and current to the control system; R P is the electrical system resistance and current control system correlation matrix; L is the correlation matrix of the inductance current of the electrical system to the control system; Q is the correlation matrix of the control state quantity of the control system to the electrical system; E is the unit diagonal matrix; is the potential of node n; ω is the state variable matrix of the control system in the power simulation system; A is the state variable coefficient matrix of the transfer function of the control system in the power simulation system; I s is the input vector of the circuit in the state equation of the electrical system in the power simulation system; B is the input coefficient matrix of the control system in the power simulation system; p is the input quantity of the control system in the power simulation system; The electrical system state equation of the power simulation system is determined as follows: Where K C K is the capacitance coefficient matrix in the electrical system state equation of the power simulation system; R K is the resistance coefficient matrix in the electrical system state equation of the power simulation system; L I is the inductance matrix in the electrical system state equation of the power simulation system; s is the input vector in the circuit of the electrical system state equation of the power simulation system; E is the unit diagonal matrix; Each element value represents the potential of the corresponding node; The circuit network state equation of the control system in the power simulation system is determined as follows: Wherein, ω' is the first-order derivative matrix of the state variable vector of the transfer function of the control system in the power simulation system; A is the state variable coefficient matrix of the transfer function of the control system in the power simulation system; B is the input coefficient matrix of the control system in the power simulation system; C' is the correlation coefficient matrix of the state variables and output quantities of the control system in the power simulation system; D is the correlation coefficient between the input and output quantities of the control system in the power simulation system; p is the input quantity of the control system in the power simulation system; q is the output quantity of the control system in the power simulation system.

9. The system according to claim 8, wherein The acquisition module is used to: The transfer function G(s) of the control system in the power simulation system is determined as follows: Where, β i =b i -a i b m , b i is the i-order coefficient of the output of the control system in the power simulation system, i∈(1~m-1); a i is the i-order coefficient of the input quantity of the control system in the power simulation system; b m is the m-order coefficient of the output of the control system in the power simulation system; s is the Laplace operator, and m is the number of state variables of the control system.

10. The system according to claim 9, wherein: The first determining module is configured to: Decompose the transfer function of the control system in the power simulation system and obtain the intermediate quantity of the transfer function of the control system in the power simulation system Obtaining the Laplace inverse transform z(t) of the intermediate quantity z(s); The state variable vector ω=[ω1…ω i …ω m ] T , where ω i is the i-th element in the state variable matrix of the control system’s transfer function, ω i =[z(t)] (i-1) ; The state variable coefficient matrix A of the transfer function of the control system in the power simulation system is determined as follows: Where a0 is the constant term coefficient of the control system input in the power simulation system; a1 is the first-order term coefficient of the control system input in the power simulation system; a2 is the second-order term coefficient of the control system input in the power simulation system; a m-1 is the m-1-order coefficient of the input quantity of the control system in the power simulation system; The input coefficient matrix B of the control system in the power simulation system is determined as follows: B＝[000…01] T Where B is an m-order column vector; The correlation coefficient matrix C' between the state variables and output quantities of the control system in the power simulation system is determined as follows: C'＝[β0 … β i … β m-1 ] The correlation coefficient D between the input and output of the control system in the power simulation system is determined as follows: D=b m Where b m is the m-th order coefficient of the output of the control system in the power simulation system.

11. The system according to claim 8, wherein The first updating unit is configured to: If the signal flowing to the control system from branch ij in the electrical system is a capacitor current signal, then update P C The element P in the τth row and the i-th column C (τ, i) = C ij , the element P in the τth row and jth column C (τ, i) = -C ij If the signal flowing from branch ij in the electrical system to the control system is the inductor current signal, then update P L The element in row τ and column i The element in row τ and column j If the signal flowing to the control system from branch ij in the electrical system is a resistance current signal, then update P R The element in row τ and column i The element in row τ and column j If the signal flowing from branch ij to the control system in the electrical system is the node voltage difference of branch ij, then update P R The element P in the τth row and the i-th column R (τ, i) = 1, the element P in the τth row and jth column R (τ, j) = -1; Among them, R ij is the resistance value of the resistance branch ij, L ij is the inductance value of the inductor branch ij, C ij is the capacitance value of the capacitor branch ij.

12. The system according to claim 8, wherein The second updating unit is configured to: If the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update the i-th row element Q(i)=const·C' and the j-th row element Q(j)=-const·C' in Q; If the signal flowing from the electrical system to the control system is a capacitor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C C The element C in row i and column i C (i,i)=const·D·C, element C in row i and column j C (i, j) = -const·D·C, the element C in the jth row and ith column C (j,i)=-const·D·C, element C in row j and column j C (j,j)=-const·D·C; If the signal flowing from the electrical system to the control system is an inductor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j If the signal flowing from the electrical system to the control system is a resistance current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j If the signal flowing from the electrical system to the control system is the node voltage difference of the branch, the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update C R The element C in row i and column i R (i,i)=const·D,element C in row i and column j R (i, j) = -const·D, element C in row j and column i R (j,i)=-const·D, element C in row j and column j R (j,j)=const·D; If the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update the i-th row element Q(i)=C' and the j-th row element Q(j)=-C' in Q; If the signal flowing from the electrical system to the control system is a capacitor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C C The element C in row i and column i C (i,i)=D·C, element C in row i and column j C (i, j) = -D·C, the element C in the jth row and ith column C (j,i)=-D·C, element C in row j and column j C (j,j)=D·C; If the signal flowing from the electrical system to the control system is an inductor current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j If the signal flowing from the electrical system to the control system is a resistance current signal, the control system signal flows to the branch ij of the electrical system, and the branch ij is a controlled voltage source branch, then update C R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j If the signal flowing from the electrical system to the control system is the node voltage difference of the branch, the control system signal flows to the branch ij of the electrical system and the branch ij is a controlled voltage source branch, then update C R The element C in row i and column i R (i,i)=D, element C in row i and column j R (i,j)=-D, element C in row j and column i R (j,i)=-D, element C in row j and column j R (j,j)=D; Where const is the derivative of the internal resistance of the controlled voltage source in the voltage source branch; C' is the correlation coefficient matrix of the control system state variables and output quantities in the power simulation system; D is the correlation coefficient of the input and output quantities of the control system in the power simulation system; C is the capacitance value of the capacitor branch; L is the inductance value of the inductor branch; and R is the resistance value of the resistor branch.

13. The system according to claim 8, wherein The resistance coefficient matrix K in the electrical system state equation of the power simulation system R , inductance matrix K L , capacitance coefficient matrix K C and the input vector I in the circuit S The acquisition process includes: Step a. Initialize the resistance coefficient matrix K in the electrical system state equation of the power simulation system R , inductance matrix K L and the capacitance coefficient matrix K C are all n+1 order zero element matrices; the input vector I in the initialization circuit S is an n+1-order zero-element vector; Step b. Update K according to the type of branch in the electrical system R , K L , K C and I S ; Step c. Eliminate the resistivity matrix K R , inductance matrix K L and the capacitance coefficient matrix K C The 0th row and column in the elimination circuit; the input vector I S Element 0 in .

14. The system according to claim 13, wherein: The step b comprises: If the branch ij in the electrical system is a capacitor branch, then update K C The element K in row i and column i C (i,i)=K C (i,i)+C, element K in row i and column j C (i,j)=K C (i,j)-C, element K in row j and column i C (j,i)=K C (j,i)-C, element K in row j and column j C (j,j)=K C (j,j)+C; If the branch ij in the electrical system is an inductive branch, then update K L The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j If the branch ij in the electrical system is a resistance branch, then update K R The element in row i and column i The element in row i and column j The element in row j and column i The element at row j and column j If the branch ij in the electrical system is a voltage source branch, then update K R The element K in row i and column i R (i,i)=K R (i,i)+Const, the element K in row i and column j R (i,j)=K R (i,j)-Const, the element K in the jth row and ith column R (j,i)=K R (j,i)-Const, the element K in the jth row and jth column R (j,j)=K R (j,j)+Const; And update the input vector I in the circuit S The i-th vector I s (i)=I s (i)+Const·U S , the jth vector I s (j) = I s (j)-Const·U S ; If the branch ij in the electrical system is a current source branch, then the input vector I in the circuit is updated S The i-th vector I s (i)=I s (i)+I' S , the jth vector I s (j) = I s (j)-I' S ; Where const is the derivative of the internal resistance of the controlled voltage source in the voltage source branch; R is the resistance value of the resistance branch, L is the inductance value of the inductance branch, and C is the capacitance value of the capacitance branch; U S is the output voltage value of the voltage source in the voltage source branch; I' S is the output current value of the current source in the current source branch.

Images