CN107203139B - Stabilization control method of multi-input multi-output nonlinear differential algebra subsystem - Google Patents
Stabilization control method of multi-input multi-output nonlinear differential algebra subsystem Download PDFInfo
- Publication number
- CN107203139B CN107203139B CN201710549039.2A CN201710549039A CN107203139B CN 107203139 B CN107203139 B CN 107203139B CN 201710549039 A CN201710549039 A CN 201710549039A CN 107203139 B CN107203139 B CN 107203139B
- Authority
- CN
- China
- Prior art keywords
- subsystem
- differential
- input
- generator
- output nonlinear
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
A stabilization control method of a multi-input multi-output nonlinear differential algebraic subsystem comprises the following steps: establishing a model of a multi-input multi-output nonlinear differential algebra subsystem; equivalently converting the model into a multi-input multi-output nonlinear ordinary differential system through differential homomorphism and state feedback, wherein the multi-input multi-output nonlinear ordinary differential system comprises each group of ordinary differential subsystems; and obtaining subsystem controllers corresponding to the ordinary differential subsystems through a backstepping method, and forming a system stabilizing controller by each group of subsystem controllers and state feedback. The invention has wider range of the controlled object, simpler and more practical backstepping method, and good effect when the model and the method are applied to the power system, thus leading the system to be more stable.
Description
Technical Field
The invention belongs to the field of engineering control, and particularly relates to a stabilization control method of a multi-input multi-output nonlinear differential algebra subsystem based on a backstepping method.
Background
The stabilizing controller is converted by adding controlled input into the subsystem controller, and can convert the nonlinear differential algebra subsystem into a nonlinear ordinary differential system by substituting the nonlinear ordinary differential algebra subsystem, which plays a crucial role in the development of the power system, however, many existing achievements relate to a single-input single-output nonlinear differential algebra subsystem. In practical application, the controlled object is often described by a multi-input multi-output nonlinear differential algebra subsystem. The multi-input multi-output nonlinear differential algebraic subsystem has the following advantages: 1. the single-input single-output model has no universality; 2. the design of the multi-input multi-output stabilizing controller is expanded from the design of the controller of each subsystem to the design of the stabilizing controller of the whole large system, so that the large system is gradually stabilized; 3. the effect of zero dynamics is not considered. Through the search of the design of the stabilizing controller of the multi-input multi-output nonlinear differential algebra subsystem, no papers and patents related to the aspect are found.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a stabilization control method of a multi-input multi-output nonlinear differential algebra subsystem.
In order to achieve the purpose, the invention adopts the following technical scheme:
establishing a model of a multi-input multi-output nonlinear differential algebra subsystem;
equivalently converting the model into a multi-input multi-output nonlinear ordinary differential system through differential homomorphism and state feedback, wherein the multi-input multi-output nonlinear ordinary differential system comprises each group of ordinary differential subsystems;
and obtaining subsystem controllers corresponding to the ordinary differential subsystems through a backstepping method, and forming a system stabilizing controller by each group of subsystem controllers and state feedback.
In order to optimize the technical scheme, the specific measures adopted further comprise:
establishing a model of a multi-input multi-output nonlinear differential algebra subsystem:
wherein, x ∈ Rn,z∈RlX and z are differential variables and algebraic variables respectively; u. ofi∈R,yi∈ R, i ═ 1, …, m, m is a positive integer tending to infinity, i is a variable subscript, uiAnd yiRespectively a control input and a control output;for relating input variables, o is the dimension of the variable, f ∈ Rn,si∈Rn,g∈Rl,hi∈ R are smooth maps and n, l are variable dimensions.
The model is equivalently converted into a multi-input multi-output nonlinear ordinary differential system:
the vector relative order (gamma) of the multi-input multi-output nonlinear differential algebra subsystem is proposed1,γ2,…,γm),γ1+γ2+…+γmN is the system dimension, γ1,…,γmRelative rank of each group of subsystems;
obtaining a differential homoembryo:
wherein the content of the first and second substances,the expression sign of each group of variables of differential isoembryo, chi and g are algebraic equations, L represents lie derivative operation, gammaiFor the consistent relative rank of the ith subsystem,Inrepresenting an n-order identity matrix;
obtaining state feedback:
wherein u is1,…,umFor each group of control inputs, r1…rmFor each set of sub-system controllers,for non-linear terms in each set of system equations, i.e.
Therefore, the multi-input multi-output nonlinear differential algebraic subsystem is equivalently converted into a multi-input multi-output nonlinear ordinary differential system.
The step of obtaining the system stabilization controller comprises:
firstly, obtaining a subsystem controller corresponding to each ordinary differential subsystem:
defining error variables Is gammaq-1 virtual controller, selecting the Lyapunov function asObtaining a subsystem controller r of the qth subsystemqI.e. byWherein the content of the first and second substances,is a design parameter;
thereby obtaining a controller of the multi-input multi-output nonlinear differential algebra subsystem:
and combining each group of subsystem controllers with state feedback to form a system stabilizing controller:
wherein u is1,…,umFor each group control input, h1,h2,…,hmIs output for each group.
The invention has the beneficial effects that: based on the prior single-input single-output nonlinear differential algebra subsystem stabilizing controller, the stabilizing controller of the multi-input multi-output nonlinear differential algebra subsystem is provided, so that the range of a controlled object is wider, the backstepping method is simpler and more practical, the model and the method have good effects when being applied to an electric power system, and the system tends to be more stable.
Drawings
Fig. 1 is a diagram of a design process of a calm controller of the present invention.
FIG. 2 is a subsystem state diagram of the present invention.
Detailed Description
The present invention will now be described in further detail with reference to the accompanying drawings.
As shown in fig. 1, the mimo nonlinear differential algebraic subsystem is equivalently converted into the mimo nonlinear ordinary differential system through differential homomorphism and state feedback, and then a stabilizing controller of a closed-loop large system is designed by using a backstepping method, so that the large system is gradually stabilized. The main structural module comprises: the system comprises a multi-input multi-output nonlinear differential algebra subsystem, a multi-input multi-output nonlinear ordinary differential system, a system stabilizing controller, a nonlinear differential algebra subsystem module, a nonlinear ordinary differential subsystem module and a subsystem controller module. When the multi-input multi-output nonlinear differential algebraic subsystem module has vector relative orders, the multi-input multi-output nonlinear differential algebraic subsystem module can be converted into the multi-input multi-output nonlinear ordinary differential system module through one differential homomorphism and state feedback. Each group of nonlinear ordinary differential subsystem modules can be designed into a subsystem controller module of each group of subsystem modules through a backstepping method. The large-system stabilizing controller module comprises all groups of subsystem controller modules, all the groups of subsystem controller modules and state feedback form the large-system stabilizing controller, and the multiple-input multiple-output nonlinear differential algebraic subsystem module is gradually stabilized through stabilizing control for derivation of a large-system Lyapunov function.
As can be seen in connection with fig. 2, the subsystem is associated with the rest of the large system by interface variables, the subsystem controller is composed of interface variable feedback, internal variable feedback and output feedback, the interface variables are generated by the effect of the rest of the large system on the subsystem and the effect of the subsystem on the rest of the large system.
The design method of the multi-input multi-output nonlinear differential algebraic subsystem stabilizing controller specifically comprises the following working steps.
First, model establishment
Modeling of multiple-input multiple-output non-linear differential algebraic subsystems, i.e.
Wherein x ∈ Rn,z∈RlAre respectively differential variable, algebraic variable, ui∈R,yi∈ R, i is 1, …, m (m is a positive integer tending to infinity, i is a variable subscript) is a control input and a control output respectively,(o is the variable dimension) is the associated input variable reflecting the influence of the rest of the larger system on the subsystem, f ∈ Rn,si∈Rn,g∈Rl,hi∈ R (n, l are variable dimensions) are smooth mappings, while the multi-input multi-output nonlinear differential algebraic subsystem is required to be exponential 1, i.e. the algebraic equation g (-) is always full rank with respect to the Jacobian matrix for the algebraic variable z, and secondly, the input variables are correlatedAnd its sufficient order derivative are both locally bounded measurable signals.
In the specific embodiment, a dual-input dual-output nonlinear differential algebraic subsystem model of the synchronous generator is established, namely
Wherein the differential variable x ═ ω, E'q,PH)TThe brackets respectively indicate the power angle of the synchronous generator, the rotating speed of a generator rotor, q-axis transient potential and the output power of a high-pressure cylinder of the steam turbine; choosing algebraic variable z ═ (P)t,θU,Id,Iq)TIn the brackets are respectively the generator active power, the generator bus voltage phase angle and the generator endD-axis and q-axis components of the current; the associated input variable is selected asThe brackets are respectively the stator current of the generator and the reactive power of the generator; in the normal operating range of the generator, It,QtAndare all locally bounded testable; control input u ═ E (E)f,Uc)TThe brackets respectively indicate excitation electromotive force and valve opening; the control output is selected as h (·) ═ (V)tAnd brackets are the generator end voltage and the generator power angle respectively; h is mechanical moment of inertia, D is damping coefficient, omega0Is synchronous angular velocity, T'd0Is the transient time constant, x, of the field windingd,xq,x′d,x′qD-axis and q-axis synchronous reactance and transient reactance, rαIs armature resistance, CH,CML,Pm0The power distribution coefficients of the high, medium and low pressure cylinders, the total output power of the steam turbine and TH∑The equivalent time constants of the high-pressure cylinder of the steam turbine comprise the time constant of the high-pressure cylinder and the time constant of a hydraulic engine of the high-pressure cylinder.
Two, proposing vector relative order
The following symbols are defined, wherein F represents an operation, InRepresenting an identity matrix of order n, i.e.
If it is paired withAnd i, j-1, …, m, k-0, …, γi-2, having the following equation (i, j are subscripts, respectively, k is 0 to γ)iA positive integer between-2, γiFor the uniform relative order of the ith subsystem, L represents the lie derivative operation):
then the integer vector (gamma)1,γ2,…,γm) Is the vector relative order of the multi-input multi-output nonlinear differential algebra subsystem.
In a specific embodiment, the vector relative order of the synchronous generator dual-input dual-output nonlinear differential algebraic subsystem model is (1, 3).
Third, model transformation
When the multi-input multi-output nonlinear differential algebra subsystem has a vector relative order (gamma)1,…,γm) And gamma is1+γ2+…+γmN (n is the system dimension, γ)1,…,γmRelative order of each set of subsystems), then there is one differential isoembryo:
and there is one state feedback:
so that the large system can be equivalently converted into a multi-input multi-output nonlinear ordinary differential system, wherein,representing symbols of each group of variables of the differential isoembryo, wherein x and g are algebraic equations; u. of1,…,umFor each group of control inputs, r1…rmFor each set of sub-system controllers,for non-linear terms in each set of system equations, i.e.
In a specific embodiment, a synchronous generator dual-input dual-output nonlinear differential algebraic subsystem model selects a differential homomorphism:
and a state feedback:
the dual-input dual-output nonlinear differential algebraic subsystem of the synchronous generator is converted into a dual-input dual-output nonlinear ordinary differential system, namely
Design of stabilizing controller
Designing a subsystem controller of each subsystem, and for a qth subsystem (q is an integer between 0 and m):
step 4-1:
defining an error variable eq,1=ξq,1,eq,2=ξq,2-αq,1Wherein αq,1Is the first virtual controller to be designed. To eq,1The derivation is carried out to obtain the result,selecting a first Lyapunov functionThenAnd isWherein c isq,1> 0 is a design parameter. … …
Step 4-k (k is 2 to γ)q-an integer between 1):
definition ofSelecting Lyapunov functionsThen α is obtained by calculationq,kAndαq,kfor the k-th virtual controller, the virtual controller,is the kth derivative of the Lyapunov function, cq,i> 0, i ═ 1, …, k is the design parameter, integer i is the parameter subscript, eq,k,eq,k+1As a function of error, i.e.
……
Step 4-. gamma.q:
Defining error variables Is gammaq-1 virtual controller, selecting the Lyapunov function asThen the subsystem controller r of the qth subsystem is availableqI.e. by
Wherein the content of the first and second substances,in order to design the parameters of the device,for error variables, a controller of a multi-input multi-output nonlinear differential algebraic subsystem is then available, i.e.
Combining each set of subsystem controllers with state feedback constitutes a large system calm controller, i.e. a large system calm controller
Wherein u is1,…,umFor each group control input, h1,h2,…,hmFor each group of output, the multiple-input multiple-output nonlinear differential algebraic subsystem is known to be asymptotically stable by differentiating the large system Lyapunov function.
In a specific embodiment, a controller of a dual-input dual-output nonlinear differential algebra subsystem of a synchronous generator is as follows:
substituting the controller into state feedback, wherein c1,1,c2,3For a parameter to be designed greater than 0, e1,1,e2,2,e2,3For error variable, the stabilizing controller of the dual-input dual-output nonlinear differential algebraic subsystem of the synchronous generator is as follows:
by derivation of a large-system Lyapunov function, the dual-input dual-output nonlinear differential algebraic subsystem of the synchronous generator is asymptotically stable.
The above is only a preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above-mentioned embodiments, and all technical solutions belonging to the idea of the present invention belong to the protection scope of the present invention. It should be noted that modifications and embellishments within the scope of the invention may be made by those skilled in the art without departing from the principle of the invention.
Claims (1)
1. A stabilization control method of a multi-input multi-output nonlinear differential algebraic subsystem comprises the following steps:
1) establishing a dual-input dual-output nonlinear differential algebraic subsystem model of the synchronous generator, i.e.
yi=hi,i=1,2
Wherein the differential variable x ═ ω, E'q,PH)TThe brackets respectively indicate the power angle of the synchronous generator, the rotating speed of a generator rotor, q-axis transient potential and the output power of a high-pressure cylinder of the steam turbine; choosing algebraic variable z ═ (P)t,θU,Id,Iq)TThe brackets are respectively the active power of the generator, the voltage phase angle of the generator bus and the d-axis and q-axis components of the generator terminal current; the associated input variable is selected asThe brackets are respectively the stator current of the generator and the reactive power of the generator; in the normal operating range of the generator, It,QtAndare all locally bounded testable; control input u ═ E (E)f,Uc)TThe brackets respectively indicate excitation electromotive force and valve opening; the control output is selected as h (·) ═ (V)tAnd brackets are the generator end voltage and the generator power angle respectively; h is mechanical moment of inertia, D is damping coefficient, omega0Is synchronous angular velocity, T'd0Is the transient time constant, x, of the field windingd,xq,x′d,x′qD-axis and q-axis synchronous reactance and transient reactance, rαIs armature resistance, CH,CML,Pm0The power distribution coefficients of the high pressure cylinder, the medium pressure cylinder and the low pressure cylinder and the total output work of the steam turbine are respectivelyRate, THΣThe equivalent time constants of the high-pressure cylinder of the steam turbine comprise the time constant of the high-pressure cylinder and the time constant of a hydraulic engine of the high-pressure cylinder;
2) the model is equivalently converted into a dual-input dual-output nonlinear ordinary differential system through differential homomorphism and state feedback, the dual-input dual-output nonlinear ordinary differential system comprises various groups of ordinary differential subsystems, and the method specifically comprises the following steps:
selecting a differential homomorphism for a synchronous generator double-input double-output nonlinear differential algebraic subsystem model:
and a state feedback:
wherein, ξ1,1,…,ξ2,3The expression sign of each group of variables of differential isoembryo, x and g are algebraic equations, L represents lie derivative operation, u is1,u2For each group of control inputs, r1,r2For each set of subsystem controllers, β1,1,β2,3Is a non-linear term in each set of system formula;
the dual-input dual-output nonlinear differential algebraic subsystem of the synchronous generator is converted into a dual-input dual-output nonlinear ordinary differential system, namely
3) The subsystem controller corresponding to the ordinary differential subsystem is obtained through a backstepping method, and each group of subsystem controllers and state feedback are combined to obtain a system stabilizing controller, which specifically comprises the following steps:
firstly, obtaining a subsystem controller corresponding to each ordinary differential subsystem:
defining error variables Is gammaq-1 virtual controller, selecting the Lyapunov function asObtaining a subsystem controller r of the qth subsystemqI.e. byWherein the content of the first and second substances,is a design parameter;
thereby obtaining a controller of the multi-input multi-output nonlinear differential algebra subsystem:
wherein the content of the first and second substances,for a design parameter greater than 0 a and,is an error variable; (gamma. rays)1,γ2,…,γm) Vector relative order of the multi-input multi-output nonlinear differential algebra subsystem;
and combining each group of subsystem controllers with state feedback to form a system stabilizing controller:
wherein u is1,…,umFor each group control input, h1,h2,…,hmOutputting for each group;
the controller of the double-input double-output nonlinear differential algebraic subsystem of the synchronous generator obtained from the above steps is as follows:
substituting the controller into state feedback, the stabilizing controller of the double-input double-output nonlinear differential algebraic subsystem of the synchronous generator is as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710549039.2A CN107203139B (en) | 2017-07-06 | 2017-07-06 | Stabilization control method of multi-input multi-output nonlinear differential algebra subsystem |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710549039.2A CN107203139B (en) | 2017-07-06 | 2017-07-06 | Stabilization control method of multi-input multi-output nonlinear differential algebra subsystem |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107203139A CN107203139A (en) | 2017-09-26 |
CN107203139B true CN107203139B (en) | 2020-09-08 |
Family
ID=59910611
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710549039.2A Active CN107203139B (en) | 2017-07-06 | 2017-07-06 | Stabilization control method of multi-input multi-output nonlinear differential algebra subsystem |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107203139B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109459932A (en) * | 2018-12-27 | 2019-03-12 | 南京信息工程大学 | The output calibration control method of multiple-input and multiple-output non-linear differential algebraic subsystem |
CN109740110A (en) * | 2018-12-27 | 2019-05-10 | 南京信息工程大学 | The generation method of the sampling observer of multiple-input and multiple-output Nonlinear Differential Algebraic Systems |
CN110619162A (en) * | 2019-09-03 | 2019-12-27 | 中国航空工业集团公司沈阳飞机设计研究所 | Stability analysis method and equipment for combined simulation platform |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102938567A (en) * | 2012-10-30 | 2013-02-20 | 东北电力大学 | Dynamic response factor-based direct-current power modulation control method of multicomputer system |
-
2017
- 2017-07-06 CN CN201710549039.2A patent/CN107203139B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN107203139A (en) | 2017-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Improvement of boom control performance for hybrid hydraulic excavator with potential energy recovery | |
CN107203139B (en) | Stabilization control method of multi-input multi-output nonlinear differential algebra subsystem | |
CN102969968B (en) | Permanent magnet synchronous motor control method | |
CN107179682B (en) | Passive load simulator and redundant moment restraining method | |
CN103051274B (en) | Variable damping-based passive control method for two-degree-of-freedom permanent magnetic synchronous motor | |
CN105786036A (en) | Control moment gyroscope framework control system and control moment gyroscope framework control method for restraining dynamic unbalance disturbance of rotor | |
CN105048917A (en) | ESO-based control method of double-fed wind power generation system integral sliding mode controller | |
CN108983801A (en) | A kind of anti-interference attitude control method of spacecraft based on counteraction flyback dynamic characteristic | |
Xu | Permanent magnet synchronous motor with linear quadraticspeed controller | |
CN109347141B (en) | Design method of grid-side terminal sliding mode controller of double-fed wind power generation system | |
Rui et al. | Fractional‐order sliding mode control for hybrid drive wind power generation system with disturbances in the grid | |
CN103939082A (en) | System and method for controlling drill rod stick-slip vibration based on active damping method | |
CN116317794A (en) | High-precision control method for electric actuator of aero-engine | |
CN103941635A (en) | System and method for restraining stick-slip vibration of drill rod | |
Rehman et al. | Adaptive control for motion synchronization of HA/EHA system by using modified MIT rule | |
CN108429501B (en) | Method for observing load disturbance of permanent magnet synchronous motor | |
CN106655962B (en) | Electric vehicle Induction Motor-Driven system control method based on extreme learning machine | |
CN112398369B (en) | Multi-motor total amount cooperative finite time anti-saturation control method | |
CN101655690A (en) | Method for simulating electric drive control system under traction working condition of electric-wheel truck | |
Belabbes et al. | Power control of a wind energy conversion system based on a doubly fed induction generator using RST and sliding mode controllers | |
Wan et al. | Practical nonlinear excitation control for a single-machine infinite-bus power system based on a detailed model | |
CN113791543A (en) | Finite time quantization control method of static var compensator based on disturbance observer | |
Zhang et al. | Research of coordination control system between nonlinear robust excitation control and governor power system stabilizer in multi-machine power system | |
CN105929682A (en) | Steam turbine generator set integrated control system | |
Bozhko et al. | Control design for electric starter-generator based on a high-speed permanent-magnet machine fed by an active front-end rectifier |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220630 Address after: 210000 room 4748, building F8, No. 9, Kechuang Avenue, Zhongshan Science Park, Jiangbei new district, Nanjing, Jiangsu Province Patentee after: Nanjing efei Scientific Instrument Co.,Ltd. Address before: 210000 No. 219 Ningliu Road, Pukou District, Nanjing City, Jiangsu Province Patentee before: Nanjing University of Information Science and Technology |
|
TR01 | Transfer of patent right |