CN102929142B - Method for designing controller of time-varying aircraft model with uncertainty - Google Patents
Method for designing controller of time-varying aircraft model with uncertainty Download PDFInfo
- Publication number
- CN102929142B CN102929142B CN201210381265.1A CN201210381265A CN102929142B CN 102929142 B CN102929142 B CN 102929142B CN 201210381265 A CN201210381265 A CN 201210381265A CN 102929142 B CN102929142 B CN 102929142B
- Authority
- CN
- China
- Prior art keywords
- delta
- matrix
- sigma
- design
- time
- 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.)
- Expired - Fee Related
Links
Abstract
The invention discloses a method for designing a controller of a time-varying aircraft model with uncertainty, which is used for solving the technical problem that the existing robust control theory lacks design steps, so the flight controller is hard to design directly. The technical scheme is as follows: the time-varying system segmentation robust stability and solvability conditions are given, selection of desired closed-loop poles of linear time-varying system state feedback is directly utilized and a constraint condition inequality direct design feedback matrix is given according to the characteristic that all the real parts of all the desired closed-loop poles are negative, so that the engineering technicians in the research field directly design the flight controller for the time-varying aircraft model with uncertainty obtained through wind tunnel or flight tests, thus solving the technical problem that the current researches only give the robust stability inequality but can not directly design the flight controller.
Description
Technical field
The present invention relates to a kind of controller design method, particularly relate to the controller design method of the uncertain time-varying model of a kind of aircraft.
Background technology
Aircraft robust control is one of emphasis problem of current international airline circle research, when high performance airplane Controller gain variations, must consider robust stability and kinds of robust control problems; Practical flight device model is the non-linear differential equation of very complicated Unknown Model structure, and in order to describe the non-linear of this complexity, people adopt wind-tunnel and flight test to obtain the test model described by discrete data usually; In order to reduce risks and reduce experimentation cost, usually carry out flight maneuver test according to differing heights, Mach number, like this, the discrete data describing aircraft test model is not a lot, and this model is very practical to the good aircraft of static stability.But the modern and following fighter plane all relaxes restriction to static stability to improve " agility ", and fighter plane requires to work near open loop critical temperature rise usually; So just require that flight control system can transaction module uncertain problem well; Following subject matter to be considered: test is obtained a certain approximate model of discrete data and describes by (1), there is Unmarried pregnancy in model in practical flight Control System Design; (2) wind tunnel test can not be carried out full scale model free flight, there is constraint, the input action selections of the selection of flight test discrete point, initial flight state, maneuvering flight etc. can not, by all non-linear abundant excitations, adopt System Discrimination gained model to there is various error; (3) flight environment of vehicle and experimental enviroment are had any different, and flow field change and interference etc. make actual aerodynamic force, moment model and test model have any different; (4) there is fabrication tolerance in execution unit and control element, also there is the phenomenons such as aging, wearing and tearing in system operation, not identical with the result of flight test; (5) in Practical Project problem, need controller fairly simple, reliable, usually need to simplify with being mathematics model person, remove the factor of some complexity; Therefore, when studying the control problem of present generation aircraft, just robustness problem must be considered.
After 1980, carry out the control theory research of multiple uncertain system in the world, the H-infinit particularly proposed by Canadian scholar Zames is theoretical, Zames thinks, based on the LQG method of state-space model, why robustness is bad, mainly because represent that uncertain interference is unpractical with White Noise Model; Therefore, when supposing that interference belongs to a certain known signal collection, Zames proposes by the norm of its corresponding sensitivity function as index, design object is under contingent worst interference, make the error of system be issued to minimum in this norm meaning, thus AF panel problem is converted into solve closed-loop system is stablized; From then on, lot of domestic and international scholar expands the research of H-infinit control method; At aeronautical chart, the method is in the exploratory stage always, U.S. NASA, and the states such as German aerospace research institute, Holland are all studied robust control method, achieves a lot of emulation and experimental result; Domestic aviation universities and colleges have also carried out a series of research to aircraft robust control method, as document (Shi Zhongke, Wu Fangxiang etc., " robust control theory ", National Defense Industry Press, in January, 2003; Su Hongye. " robust control basic theory ", Science Press, in October, 2010) introduce, but these results and the distance of practical application also differ very large, are difficult to directly design practical flight controller and apply; Particularly a lot of research only gives Robust Stability according to Lyapunov theorem, becomes robust Controller Design step when can not obtain specific implementation, does not have to solve the technical matters of directly design robust flight controller.
Summary of the invention
Being difficult to directly design the technical deficiency of flight controller in order to overcome existing robust control theory shortage design procedure, the invention provides the controller design method of the uncertain time-varying model of a kind of aircraft; This method provide the design conditions of the uncertain time-varying model Robust Stability Controller of real system, the closed loop of linear time varying system feedback of status is directly utilized to expect the selection of poles, and expect that the real part of limit is all the feature of negative according to all closed loops, give qualifications inequality direct design of feedback matrix, the uncertain time-varying model of aircraft that can obtain wind-tunnel or flight test directly designs flight controller, solves current research and only provides robust stability inequality and the technical matters that directly cannot design flight controller.
The technical solution adopted for the present invention to solve the technical problems is: the controller design method of the uncertain time-varying model of a kind of aircraft, is characterized in comprising the following steps:
Step one, obtained by wind-tunnel or flight test under assigned altitute, Mach number condition containing probabilistic dummy vehicle be:
In formula, x (t) ∈ R
n, u (t) ∈ R
mbe respectively state and input vector, A (t), B (t) are known matrix of coefficients, and Δ A (t), Δ B (t) are matrix of coefficients unknown portions; According to known A (t), the variation range classification of B (t), is namely expressed as A (t), B (t) in different time sections:
In formula, A
0i, B
0ifor known constant matrices, Δ A
0iΔ B
0ifor unknown matrix, t
ij, T
ijfor time constant, r, p are positive integer, and i, j are that the A (t) of subscript, different time sections, B (t) expression formula form can be identical;
In time period t
ij≤ t<t
ij+ T
ijin, flight controller is: u (t)=-K
ix (t)
In formula, K
ifor constant feedback matrix;
Bring in (1) formula, have:
Step 2, choose (A
0i-B
0ik
0i) eigenwert different and real part is negative, design of feedback matrix K
imake to satisfy condition:
This controller makes
Robust stability;
In formula, M
ifor the matrix of a linear transformation,
σ
ik, ω
ik(k=1,2 ..., n) be real number, j ω
ix(k=1,2 ..., n) represent imaginary number, diag is diagonal matrix symbol,
Δ A
0i-Δ B
0ik
iusually Δ A is assumed to be
0i-Δ B
0ik
i=H
if
iw
i, H
i, W
iall be assumed to be matrix, 0<F
i≤ I, I=diag [1,1 ..., 1] and be unit battle array.
The invention has the beneficial effects as follows: by time-varying system segmentation robust stability solution conditions provided by the invention, the closed loop of linear time varying system feedback of status is directly utilized to expect the selection of poles, and expect that the real part of limit is all the feature of negative according to all closed loops, give qualifications inequality direct design of feedback matrix, the engineering technical personnel of this research field are made directly to design flight controller to the uncertain time-varying model of aircraft that wind-tunnel or flight test obtain, solve current research and only provide robust stability inequality and the technical matters that directly cannot design flight controller.
Below in conjunction with embodiment, the present invention is elaborated.
Embodiment
The controller design method concrete steps of the uncertain time-varying model of aircraft of the present invention are as follows:
1, obtained by wind-tunnel or flight test under assigned altitute, Mach number condition containing probabilistic dummy vehicle be:
In formula, x (t) ∈ R
n, u (t) ∈ R
mbe respectively state and input vector, A (t), B (t) are known matrix of coefficients, and Δ A (t), Δ B (t) are matrix of coefficients unknown portions; According to known A (t), the variation range classification of B (t), is namely expressed as A (t), B (t) in different time sections:
In formula, A
0i, B
0ifor known constant matrices, Δ A
0iΔ B
0ifor unknown matrix, t
ij, T
ijfor time constant, r, p are positive integer, and i, j are that the A (t) of subscript, different time sections, B (t) expression formula form can be identical;
In time period t
ij≤ t<t
ij+ T
ijin, flight controller is: u (t)=-K
ix (t)
In formula, K
ifor constant feedback matrix;
Bring in (1) formula, have:
2, (A is chosen
0i-B
0ik
0i) eigenwert different and real part is negative, design of feedback matrix K
imake to satisfy condition:
This controller makes
Robust stability;
In formula, M
ifor the matrix of a linear transformation,
σ
ik, ω
ik(k=1,2 ..., n) be real number, j ω
ik(k=1,2 ..., n) represent imaginary number, diag is diagonal matrix symbol,
Δ A
0i-Δ B
0ik
iusually Δ A is assumed to be
0i-Δ B
0ik
i=H
if
iw
i, H
i, W
iall be assumed to be matrix, 0<F
i≤ I, I=diag [1,1 ..., 1] and be unit battle array;
Getting Flight Altitude Moving state variable is x=[q α θ]
t, input variable is u=δ
e, wherein q is rate of pitch, and α is the air-flow angle of attack, and θ is the angle of pitch, δ
efor elevating rudder drift angle; In time period 20≤t<100, State Equation Coefficients matrix is:
Uncertain part is:
Closed loop is selected to expect limit and A
01-B
01k
1eigenwert σ (A
01-B
01k
1)=diag [-0.5 ,-1 ,-2], can obtain:
Controller is: K
1=[-0.3794 1.5382-0.6231].
Claims (1)
1. a controller design method for the uncertain time-varying model of aircraft, is characterized in that comprising the following steps:
Step one, obtained by wind-tunnel or flight test under assigned altitute, Mach number condition containing probabilistic dummy vehicle be:
In formula, x (t) ∈ R
n, u (t) ∈ R
mbe respectively state and input vector, A (t), B (t) are known matrix of coefficients, and Δ A (t), Δ B (t) are matrix of coefficients unknown portions; According to known A (t), the variation range classification of B (t), is namely expressed as A (t), B (t) in different time sections:
In formula, A
0i, B
0ifor known constant matrices, Δ A
0iΔ B
0ifor unknown matrix, t
ij, T
ijfor time constant, r, p are positive integer, and i, j are that the A (t) of subscript, different time sections, B (t) expression formula form is identical;
In time period t
ij≤ t < t
ij+ T
ijin, flight controller is: u (t)=-K
ix (t)
In formula, K
ifor constant feedback matrix;
Bring in (1) formula, have:
Step 2, choose (A
0i-B
0ik
i) eigenwert different and real part is negative, design of feedback matrix K
imake to satisfy condition:
This controller makes
Robust stability;
In formula, M
ifor the matrix of a linear transformation,
σ
ik, ω
ik, k=1,2 ..., n is real number, j ω
ik, k=1,2 ..., n represents imaginary number, and diag is diagonal matrix symbol,
Δ A
0i-Δ B
0ik
iusually Δ A is assumed to be
0i-Δ B
0ik
i=H
if
iw
i, H
i, W
iall be assumed to be known matrix, 0 < F
i≤ I, I=diag [1,1 ..., 1] and be unit battle array;
Getting Flight Altitude Moving state variable is
Input variable is u=δ
e, wherein q is rate of pitch, and α is the air-flow angle of attack,
for the angle of pitch, δ
efor elevating rudder drift angle.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210381265.1A CN102929142B (en) | 2012-10-10 | 2012-10-10 | Method for designing controller of time-varying aircraft model with uncertainty |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210381265.1A CN102929142B (en) | 2012-10-10 | 2012-10-10 | Method for designing controller of time-varying aircraft model with uncertainty |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102929142A CN102929142A (en) | 2013-02-13 |
CN102929142B true CN102929142B (en) | 2015-04-01 |
Family
ID=47643971
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210381265.1A Expired - Fee Related CN102929142B (en) | 2012-10-10 | 2012-10-10 | Method for designing controller of time-varying aircraft model with uncertainty |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102929142B (en) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103809451B (en) * | 2014-02-28 | 2016-03-23 | 西安费斯达自动化工程有限公司 | The multistage PID controller design method of aircraft multiloop model bunch compound root locus |
CN103792848B (en) * | 2014-02-28 | 2016-03-30 | 西安费斯达自动化工程有限公司 | Longitudinal Flight model cluster man-machine loop compound root locus multistage PID robust Controller Design method |
CN103809452B (en) * | 2014-02-28 | 2016-03-23 | 西安费斯达自动化工程有限公司 | Longitudinal Flight model cluster Flutter Suppression compound root locus multistage PID robust Controller Design method |
CN103809447B (en) * | 2014-02-28 | 2016-03-02 | 西安费斯达自动化工程有限公司 | Aircraft multiloop model bunch combination frequency controller design method |
CN103809449B (en) * | 2014-02-28 | 2016-03-23 | 西安费斯达自动化工程有限公司 | Aircraft multiloop model bunch Flutter Suppression Composite PID robust Controller Design method |
CN103809450B (en) * | 2014-02-28 | 2016-04-27 | 西安费斯达自动化工程有限公司 | Multiloop dummy vehicle bunch Flutter Suppression compound root locus multistage PID robust Controller Design method |
CN103809453B (en) * | 2014-02-28 | 2016-06-29 | 西安费斯达自动化工程有限公司 | Longitudinal Flight model cluster man-machine loop's compound root locus compensates robust Controller Design method |
CN108052115B (en) * | 2017-12-06 | 2019-11-08 | 浙江工业大学 | It is a kind of based on it is asymmetric when constant obstacle liapunov function quadrotor total state be limited backstepping control method |
CN107831670B (en) * | 2017-12-06 | 2019-11-05 | 浙江工业大学 | It is a kind of based on it is asymmetric when constant obstacle liapunov function quadrotor export limited backstepping control method |
CN107831671B (en) * | 2017-12-06 | 2019-11-08 | 浙江工业大学 | A kind of limited backstepping control method of quadrotor output based on asymmetric time-varying obstacle liapunov function |
CN108037662B (en) * | 2017-12-06 | 2019-11-08 | 浙江工业大学 | A kind of limited backstepping control method of quadrotor output based on Integral Sliding Mode obstacle liapunov function |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6643569B2 (en) * | 2001-03-30 | 2003-11-04 | The Regents Of The University Of Michigan | Method and system for detecting a failure or performance degradation in a dynamic system such as a flight vehicle |
CN102591212A (en) * | 2012-03-01 | 2012-07-18 | 北京航空航天大学 | Method for observing longitudinal motion state of aircraft by using time varying measurement delay output signal |
CN102707624A (en) * | 2012-05-31 | 2012-10-03 | 西北工业大学 | Design method of longitudinal controller region based on conventional aircraft model |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05147589A (en) * | 1991-11-28 | 1993-06-15 | Mitsubishi Heavy Ind Ltd | Robust flight control device |
-
2012
- 2012-10-10 CN CN201210381265.1A patent/CN102929142B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6643569B2 (en) * | 2001-03-30 | 2003-11-04 | The Regents Of The University Of Michigan | Method and system for detecting a failure or performance degradation in a dynamic system such as a flight vehicle |
CN102591212A (en) * | 2012-03-01 | 2012-07-18 | 北京航空航天大学 | Method for observing longitudinal motion state of aircraft by using time varying measurement delay output signal |
CN102707624A (en) * | 2012-05-31 | 2012-10-03 | 西北工业大学 | Design method of longitudinal controller region based on conventional aircraft model |
Non-Patent Citations (1)
Title |
---|
"不确定时滞系统的区间时滞相关鲁棒镇定";陈杰 等;《系统工程与电子技术》;20081231;第30卷(第12期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN102929142A (en) | 2013-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102929142B (en) | Method for designing controller of time-varying aircraft model with uncertainty | |
CN102929138B (en) | Method for designing aircraft controller with nonlinearity | |
CN102929139B (en) | Method for designing controller of multi-interval aircraft model | |
CN102929141B (en) | Aircraft time lag varying model approximation and controller designing method | |
CN102929140B (en) | Method for designing approximation and controller of time lag aircraft model | |
CN102929134A (en) | Method for designing controller of multiple time lag aircraft model | |
CN102929130B (en) | Robust flight controller design method | |
CN102929144B (en) | The controller design method of aircraft interval model | |
CN102929128B (en) | Method for designing controller of aircraft model with uncertainty | |
Frink et al. | Computational aerodynamic modeling tools for aircraft loss of control | |
CN103592847A (en) | Hypersonic aerocraft nonlinear control method based on high-gain observer | |
CN105676852A (en) | Small unmanned helicopter adaptive attitude control method of dynamic-model-free structure | |
CN109840349A (en) | A kind of fixed wing aircraft gust response modeling and analysis methods | |
Blonigan | Least squares shadowing for sensitivity analysis of large chaotic systems and fluid flows | |
CN102929143B (en) | Control design method for aircraft time lag model | |
CN102929129B (en) | Aircraft controller designing method containing multiple time lag varying models | |
CN102929135B (en) | Method for designing controller for aircraft time lagging and varying model | |
CN102929145B (en) | Design method of aerocraft multi-time delay time-varying model approximation and controller | |
Moreno et al. | A modification to the enhanced correction factor technique to correlate with experimental data | |
CN102929136B (en) | Aircraft multi-time delay model approach and controller design method | |
CN102929137B (en) | Aircraft controller designing method containing distribution parameter items | |
Kopf et al. | On the effect of enforcing stability in model predictive control for gust load alleviation | |
Stanney et al. | Response of a tethered aerostat to simulated turbulence | |
CN103777523B (en) | Aircraft multiloop model bunch Composite PID robust Controller Design method | |
Tokutake et al. | Disturbance preview controller and its application to a small UAV |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150401 Termination date: 20211010 |