CN102880188B - Aircraft modeling method based on maximum information-reliability online identification criterion - Google Patents
Aircraft modeling method based on maximum information-reliability online identification criterion Download PDFInfo
- Publication number
- CN102880188B CN102880188B CN201210382610.3A CN201210382610A CN102880188B CN 102880188 B CN102880188 B CN 102880188B CN 201210382610 A CN201210382610 A CN 201210382610A CN 102880188 B CN102880188 B CN 102880188B
- Authority
- CN
- China
- Prior art keywords
- omega
- overbar
- eta
- sigma
- partiald
- 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
Landscapes
- Feedback Control In General (AREA)
Abstract
The invention discloses an aircraft modeling method based on a maximum information-reliability online identification criterion. The technical problems that the proving correctness of aerodynamic models and parameters verified by a flight test is poor caused by a simplified model identification criterion are solved. The technical scheme includes that data lengths are considered and reliable parameters are introduced in a maximum information criterion, an online modeling criterion is corrected, the division of U-D of Rj and Rj+1 is evaluated for the measurement variance, and a scalar online model selection and verified discriminant is obtained. The method has the advantages that aerodynamic models and moment models of an aircraft can be established on line directly in accordance with flight test data, and the technical problems of aircraft high attack angle model establishment and verification according to flight tests caused by a simplified model determination criterion are solved.
Description
Technical field
The present invention relates to a kind of Modeling of Vehicle method, particularly relate to a kind of Modeling of Vehicle method based on maximum fault information-confidence level on-line identification criterion.
Background technology
According to aircraft aerodynamic model and parameter, not only can determine the control stability of aircraft, also can be ground and provide correct mathematical model with aerial emulator; Wind tunnel experiment and the calculated results of checking aircraft aerodynamic parameter; For design and the improvement of aircraft control system provides master data; Identify the flight quality of sizing aircraft; The flight quality of research high performance airplane; Carry out crash analysis of aviation accident etc.; Set up exactly aircraft mathematical model problem with completely different by the theoretical method of the modellings by mechanism such as philosophy, theorem, mainly according to the input and output data of experiment gained, set up model, its main theoretical basis of the reform of Chinese economic structure is Nonlinear Systems Identification and nonlinear flight dynamics; When aircraft is done low incidence microvariations flight, flight force and moment mould can be got once item, i.e. Bryan model representation by platform labor series expansion.When Mach number, height one timing, this model is Linear Time Invariant model, and this model, because form simply and is always used till today, becomes the foundation stone of pneumatic mathematical model; Adopt this model, aerocraft system identification has been estimated the systematic parameter of known mathematical model with regard to having become; Modern combat aircraft, tactical missile doing that wartime requirement is more motor-driven, fault speed tailspin even, its angle of attack can be from tens degree, tens degree until spend 100 more, can not adopt linear model; The caused non-permanent downwash flow field in body-shedding vortex, separated whirlpool that Aircraft at High Angle of Attack forms, steady model can not be suitable for again.Research non-permanent, nonlinear aerodynamic model of aircraft under At High Angle of Attack become current aircraft development in the urgent need to problem.Yet, the identification complex of Nonlinear Aerodynamic, it is general Nonlinear Systems Identification problem, the funtcional relationship between input quantity and state is difficult to determine, need to carries out identification to model; The key of Model Distinguish is modeling criterion and optimization algorithm, and for given version, application modeling criterion is determined the optimum exponent number of model and select optimization model from candidate model; Because measured data contains noise, modeling criterion can not only be investigated the error of fitting size to available data, and considers other factors, otherwise will make model incorrect; Conventionally, modeling criterion should be able to make the model optimizing have following characteristics: 1. the model existing flying quality of matching well; 2. model is every an obvious physical significance; 3. model can be predicted the measured data under conditions of similarity; 4. under the suitable condition of performance, order is minimum; The most frequently used identification Method is method of gradual regression, and its principle is item by item the predictor that affects conspicuousness to be selected into, and the little factor of impact is rejected, and sets up the method for regression equation; This method is calculated simple, practical; But this method has two obvious shortcomings: the one, and choice criteria is determined by people, and does not provide the confidence level of result; The 2nd, the accumulation of error is large, easily leaks choosing and falsely drops; For this reason, people usually adopt maximum fault information criterion AIC method to the dummy vehicle identification problem of having relatively high expectations, but the method processing speed is slow, and signal to noise ratio (S/N ratio) is hour Model Distinguish low precision; Due under nonlinear situation; can only carry out numerical integration to aircraft nonlinear equation; carry out sensitivity matrix calculating and iterative computation; thereby make complexity and the calculated amount of calculating more much higher than Linear Estimation; also make the matching between model output and experimental data become more difficult simultaneously; for online Approach For Identification of Model Structure problem, existing AIC criterion does not directly provide result, and aerodynamic model and Verification that other short-cut method usually can cause flight test to provide are incorrect.
Summary of the invention
The aerodynamic model and the poor deficiency of Verification correctness that in order to overcome existing simplified model identification criterion, cause flight test to provide, the invention provides a kind of Modeling of Vehicle method based on maximum fault information-confidence level on-line identification criterion.The method is by analyzing the impact of data length and introducing confidence level parameter, maximum fault information criterion is revised, obtained new model on-line identification criterion, by new criterion, set up index modeling, directly can, for flight test line modeling and the modelling verification of aircraft, can avoid the technical matters according to flight test is set up and checking aircraft at high angle of attack model exists.
The technical solution adopted for the present invention to solve the technical problems is: a kind of Modeling of Vehicle method based on maximum fault information-confidence level on-line identification criterion, is characterized in comprising the following steps:
The state equation of the aircraft candidate family that step 1, flight test are to be determined is
Observation equation is
(1), in (2) formula, x (t) is n dimension state vector; Y (t) is m dimension observation vector; f{f
0[x (t), Ω
0], f
1[x (t), θ
1] ..., f
q[θ
q, x (t)], t}, g{g
0[x (t), Ω
0], g
1[x (t), θ
1] ..., g
q[θ
q, x (t)], t} is the known model structure function to be determined of expression formula, f
0[x (t), Ω
0], g
0[x (t), Ω
0] be the model that must be selected into according to physical concept, f
i[x (t), θ
i], g
i[x (t), θ
i] (i=1,2 ..., q) be candidate family, z (t
k) be at t
kconstantly to y (t
k) measured value; Ω is the parameter vector of unknown dimension, Ω
0parameter vector for known dimension; V (k), for measuring noise, supposes that variance is R
kzero-mean white Gaussian noise; f
i[x (t), θ
i], g
i[x (t), θ
i] (i=1,2 ..., q) whether in model, occur and Ω
0, θ
i(i=1,2 ..., value q) needs identification, and q is known candidate family number;
Conventionally higher to the model structure accuracy requirement of aircraft, the present invention provides following maximum fault information-confidence level criterion:
In formula, L
1[x (t), Ω
1, t], L
2[x (t), Ω
2, t] and for getting different parameter vector Ω
1and Ω
2maximum likelihood function, δ is given arithmetic number, represent confidence level parameter,
Step 2, according to supposition f
0[x (t), Ω
0], g
0[x (t), Ω
0], Ω
0=Ω
0by optimization algorithm, be selected into model, and calculated by following algorithm iteration:
Make j=0,1,2 ..., q, supposes f
j[x (t), θ
j], g
j[x (t), θ
j], Ω
jbe selected into model, selected in such a way other candidate family:
Ask (4) formula maximum value, iterative computation:
And
(5), in (6) formula:
If
And N+1 secondary data is come, θ
j+1be selected into or rejecting and verification condition are: when
Time, θ
j+1, f
j+1[x (t), θ
j+1], g
j+1[x (t), θ
j+1] be selected into model, and
Otherwise reject f
j+1[x (t), θ
j+1], g
j+1[x (t), θ
j+1] candidate item, and Ω
j+1=Ω
j;
(7) in formula: ε is given positive number,
Step 3, because the dimension m of aircraft measurement vector y is larger, adopt Gram-Schmidt orthogonalization method to R
nj, R
(N+1) j, R
n (j+1)and R
(N+1) (j+1)carry out U-D decomposition, R
nj, R
(N+1) j, R
n (j+1)and R
(N+1) (j+1)u-D be respectively:
In formula, U
rNj, U
r (N+1) j, U
rN (j+1), U
r (N+1) (j+1)for unit upper triangular matrix; D
rNj=diag[d
rNj(1), d
rNj(2) ..., d
rNj(m)], D
rN (j+1)=diag[d
rN (j+1)(1), d
rN (j+1)(2) ..., d
rN (j+1)(m)], D
r (N+1) j=diag[d
r (N+1) j(1), d
r (N+1) j(2) ..., d
r (N+1) j(m)], D
r (N+1) (j+1)=diag[d
r (N+1) (j+1)(1), d
r (N+1) (j+1)(2) ..., d
r (N+1) (j+1)(m)]; Diag is diagonal angle symbol;
The maximum fault information criterion of modelling verification is write as: when
During establishment, θ
i+1, f
j+1[x (t), θ
j+1], g
j+1[x (t), θ
j+1] to be selected into model correct, and
Otherwise reject f
j+1[x (t), θ
j+1], g
j+1[x (t), θ
j+1] candidate item, and Ω
j+1=Ω
j.
The invention has the beneficial effects as follows: due to by considering data length in maximum fault information criterion and introducing confidence level parameter, revised line modeling criterion; To measuring variance, estimate R
jand R
j+1u-D decompose, the online Model Selection of scalar and checking discriminant have been obtained, be convenient to directly according to test flight data, set up online aircraft aerodynamic force, moment model, avoid simplified model criterion to cause, avoided the technical matters according to flight test is set up and checking aircraft at high angle of attack model exists.
Below in conjunction with embodiment, the present invention is elaborated.
Embodiment
The Modeling of Vehicle method concrete steps that the present invention is based on maximum fault information-confidence level on-line identification criterion are as follows:
1, many aircraft are less than 60 conventional candidate family forms while spending at the angle of attack and are:
(1), in (2) formula,, x (t) is n dimension state vector; Y (t) is m dimension observation vector, Φ (Ω
0) f
0[x (t)], Ψ (Ω
0) g
0[x (t)] model for being selected into according to physical concept, θ
if
i[x (t)], θ
ig
i[x (t)] (i=1,2 ..., q) be candidate family, z (t
k) be at t
kconstantly to y (t
k) measured value; Ω is the parameter vector of unknown dimension, Ω
0parameter vector for known dimension; V (k), for measuring noise, supposes that variance is R
kzero-mean white Gaussian noise; θ
if
i[x (t)], θ
ig
i[x (t)] (i=1,2 ..., q) whether in model, occur and Ω
0, θ
i(i=1,2 ..., value q) needs identification., q is known candidate family number;
Conventionally higher to the model structure accuracy requirement of aircraft, the present invention provides following maximum fault information-confidence level criterion:
In formula, L
1[x (t), Ω
1, t], L
2[x (t), Ω
2, t] and for getting different parameter vector Ω
1and Ω
2maximum likelihood function, δ is given arithmetic number, represent confidence level parameter,
2, according to supposition f
0[x (t), Ω
0], g
0[x (t), Ω
0], Ω
0=Ω
0by optimization algorithm, be selected into model, and calculated by following algorithm iteration:
Make j=0,1,2 ..., q, supposes f
j[x (t), θ
j], g
j[x (t), θ
j], Ω
jbe selected into model, selected in such a way other candidate family:
Ask (4) formula maximum value, iterative computation:
And
(5), in (6) formula:
If
And N+1 secondary data is come, θ
j+1be selected into or rejecting and verification condition are: when
Time, θ
i+1, θ
jf
j[x (t)], θ
jg
j[x (t)] is selected into model, and
Otherwise reject θ
jf
j[x (t)], θ
jg
j[x (t)] candidate item, and Ω
j+1=Ω
j;
(7) in formula: ε is given positive number,
3, the dimension m of aircraft measurement vector y is larger conventionally, adopts Gram-Schmidt orthogonalization method to R
nj, R
(N+1) j, R
n (j+1)and R
(N+1) (j+1)carry out U-D decomposition, R
nj, R
(N+1) j, R
n (j+1)and R
(N+1) (j+1)u-D be respectively:
In formula, U
rNj, U
r (N+1) j, U
rN (j+1), U
r (N+1) (j+1)for unit upper triangular matrix; D
rNj=diag[d
rNj(1), d
rNj(2) ..., d
rNj(m)], D
rN (j+1)=diag[d
rN (j+1)(1), d
rN (j+1)(2) ..., d
rN (j+1)(m)], D
r (N+1) j=diag[d
r (N+1) j(1), d
r (N+1) j(2) ..., d
r (N+1) j(m)], D
r (N+1) (j+1)=diag[d
r (N+1) (j+1)(1), d
r (N+1) (j+1)(2) ..., d
r (N+1) (j+1)(m)]; Diag is diagonal angle symbol;
The maximum fault information criterion of modelling verification can be write as: when
During establishment, θ
j+1, θ
jf
j[x (t)], θ
jg
jit is correct that [x (t)] is selected into model, and
Otherwise reject θ
jf
j[x (t)], θ
jg
j[x (t)] candidate item, and Ω
j+1=Ω
j.
Claims (1)
1. the Modeling of Vehicle method based on maximum fault information-confidence level on-line identification criterion, is characterized in that comprising the following steps:
The state equation of the aircraft candidate family that step 1, flight test are to be determined is
Observation equation is
(1), in (2) formula, x (t) is n dimension state vector; Y (t) is m dimension observation vector; f{f
0[x (t), Ω
0], f
1[x (t), θ
1] ..., f
q[θ
q, x (t)], t}, g{g
0[x (t), Ω
0], g
1[x (t), θ
1] ..., g
q[θ
q, x (t)], t} is the known model structure function to be determined of expression formula, f
0[x (t), Ω
0], g
0[x (t), Ω
0] be the model that must be selected into according to physical concept, f
i[x (t), θ
i], g
i[x (t), θ
i], i=1,2 ..., q is candidate family, z (t
k) be at t
kconstantly to y (t
k) measured value; Ω is the parameter vector of unknown dimension, Ω
0parameter vector for known dimension; V (k), for measuring noise, supposes that variance is R
kzero-mean white Gaussian noise; f
i[x (t), θ
i], g
i[x (t), θ
i], i=1,2 ..., whether q occurs and Ω in model
0, θ
i, i=1,2 ..., the value of q needs identification, and q is known candidate family number;
Conventionally higher to the model structure accuracy requirement of aircraft, provide following maximum fault information-confidence level criterion:
In formula, L
1[x (t), Ω
1, t], L
2[x (t), Ω
2, t] and for getting different parameter vector Ω
1and Ω
2maximum likelihood function, δ is given arithmetic number, represent confidence level parameter,
Step 2, according to supposition f
0[x (t), Ω
0], g
0[x (t), Ω
0], Ω
0=Ω
0by optimization algorithm, be selected into model, and calculated by following algorithm iteration:
Make j=0,1,2 ..., q, supposes f
j[x (t), θ
j], g
j[x (t), θ
j], Ω
jbe selected into model, selected in such a way other candidate family:
Ask (4) formula maximum value, iterative computation:
And
(5), in (6) formula:
If
And N+1 secondary data is come, θ
j+1be selected into or rejecting and verification condition are: when
Time, θ
j+1, f
j+1[x (t), θ
j+1], g
j+1[x (t), θ
j+1] be selected into model, and
Otherwise reject f
j+1[x (t), θ
j+1], g
j+1[x (t), θ
j+1] candidate item, and Ω
j+1=Ω
j;
(7) in formula: ε is given positive number,
Step 3, because the dimension m of aircraft measurement vector y is larger, adopt Gram-Schmidt orthogonalization method to R
nj, R
(N+1) j, R
n (j+1)and R
(N+1) (j+1)carry out U-D decomposition, R
nj, R
(N+1) j, R
n (j+1)and R
(N+1) (j+1)u-D be respectively:
In formula, U
rNj, U
r (N+1) j, U
rN (j+1), U
r (N+1) (j+1)for unit upper triangular matrix;
D
RNj=diag[d
RNj(1),d
RNj(2),…,d
RNj(m)],D
RN(j+1)=diag[d
RN(j+1)(1),d
RN(j+1)(2),…,d
RN(j+1)(m)],
D
R(N+1)j=diag[d
R(N+1)j(1),d
R(N+1)j(2),…,d
R(N+1)j(m)],
D
r (N+1) (j+1)=diag[d
r (N+1) (j+1)(1), d
r (N+1) (j+1)(2) ..., d
r (N+1) (j+1)(m)]; Diag is diagonal angle symbol;
The maximum fault information criterion of modelling verification is write as: when
During establishment, θ
j+1, f
j+1[x (t), θ
j+1], g
j+1[x (t), θ
j+1] to be selected into model correct, and
Otherwise reject f
j+1[x (t), θ
j+1], g
j+1[x (t), θ
j+1] candidate item, and Ω
j+1=Ω
j.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210382610.3A CN102880188B (en) | 2012-10-11 | 2012-10-11 | Aircraft modeling method based on maximum information-reliability online identification criterion |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210382610.3A CN102880188B (en) | 2012-10-11 | 2012-10-11 | Aircraft modeling method based on maximum information-reliability online identification criterion |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102880188A CN102880188A (en) | 2013-01-16 |
CN102880188B true CN102880188B (en) | 2014-11-26 |
Family
ID=47481557
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210382610.3A Expired - Fee Related CN102880188B (en) | 2012-10-11 | 2012-10-11 | Aircraft modeling method based on maximum information-reliability online identification criterion |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102880188B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108255063B (en) * | 2018-01-26 | 2021-03-19 | 深圳禾苗通信科技有限公司 | Small rotor unmanned aerial vehicle system modeling method based on closed-loop subspace identification |
CN111611648B (en) * | 2020-05-09 | 2022-11-08 | 中国人民解放军96901部队 | Method for evaluating maximum flight distance capability of boosting-gliding aircraft |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298332A (en) * | 2010-12-30 | 2011-12-28 | 清华大学 | Compound identification method of unmanned helicopter flight dynamical model |
-
2012
- 2012-10-11 CN CN201210382610.3A patent/CN102880188B/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298332A (en) * | 2010-12-30 | 2011-12-28 | 清华大学 | Compound identification method of unmanned helicopter flight dynamical model |
Non-Patent Citations (10)
Title |
---|
"Flight Flutter Modal Parameters Identification with Atmospheric Turbulence Excitation Based on Wavelet Transformation";Zhang Bo 等;《Chinese Journal of Aeronautics》;20071231;第394-401页 * |
"Modeling and Combined Control of Aerial Vehicle Based on Aerodynamic Force and Lateral Thrust";WANG Zheng-jie 等;《Proceedings of the 2010 International Conference on Modelling, Identification and Control》;20101231;第106-111页 * |
"飞机大攻角空间机动气动力建模研究";汪清 等;《航空学报》;20040930;第25卷(第5期);第447-450页 * |
"飞行器气动参数辨识进展";蔡金狮;《力学进展》;19871125;第17卷(第4期);第467-478页 * |
WANG Zheng-jie 等."Modeling and Combined Control of Aerial Vehicle Based on Aerodynamic Force and Lateral Thrust".《Proceedings of the 2010 International Conference on Modelling, Identification and Control》.2010,第106-111页. * |
Zhang Bo 等."Flight Flutter Modal Parameters Identification with Atmospheric Turbulence Excitation Based on Wavelet Transformation".《Chinese Journal of Aeronautics》.2007,第394-401页. * |
李乃宏 等."有控飞行器非线性气动模型结构的确定(Ⅰ)".《飞行力学》.1993,第11卷(第2期),第26-33页. * |
李乃宏 等."有控飞行器非线性气动模型结构的确定(Ⅰ)".《飞行力学》.1993,第11卷(第2期),第26-33页. * |
汪清 等."飞机大攻角空间机动气动力建模研究".《航空学报》.2004,第25卷(第5期),第447-450页. * |
蔡金狮."飞行器气动参数辨识进展".《力学进展》.1987,第17卷(第4期),第467-478页. * |
Also Published As
Publication number | Publication date |
---|---|
CN102880188A (en) | 2013-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105115692B (en) | The pneumatic Forecasting Methodology that a kind of CFD numerical simulations are combined with wind tunnel test | |
CN110851957A (en) | Atmospheric data sensing system resolving method based on deep learning | |
CN105843073A (en) | Method for analyzing wing structure aero-elasticity stability based on aerodynamic force uncertain order reduction | |
CN103020456B (en) | A kind of To Aerodynamic Coefficients Identification of Vehicle error burst evaluation method | |
Cenko | Store separation lessons learned during the last 30 years | |
CN108009378B (en) | The structure changing damage appraisal procedure of guided wave HMM based on equality initialization GMM | |
CN109596145A (en) | A kind of dynamic testing method and system of vehicle mounted guidance | |
CN110516394A (en) | Aero-engine steady-state model modeling method based on deep neural network | |
CN103955714A (en) | Navy detection model construction method and system and navy detection method | |
CN102999667A (en) | Fuzzy regression based airport noise predicating method and system | |
CN104899448A (en) | Adaptive compensation method for static localization scheme of ensemble Kalman filter | |
CN102880188B (en) | Aircraft modeling method based on maximum information-reliability online identification criterion | |
CN104571087A (en) | Diagnostic determination method for spacecraft control system under influence of noise | |
CN102880057B (en) | Aircraft modeling method based on variable data length maximum information criterion | |
CN102867130B (en) | Aircraft modeling method based on variable measurement number maximum information criterion | |
CN115204041A (en) | Large aircraft approach and approach stage oil consumption prediction method | |
CN104462022A (en) | Aircraft dynamics system parameter identifiability analitical method | |
CN116992577B (en) | Simulation method, system, equipment and storage medium of cross-medium aircraft | |
CN102867102B (en) | Modeling method of aircraft based on maximum information quantity-reliability rule | |
CN102890452B (en) | Aircraft modeling method based on maximum information content-credibility criterion of variable metrical data | |
CN102867129B (en) | Aircraft modeling method based on variable data length maximum information amount-credibility criterion | |
CN102938002B (en) | Aircraft modeling method based on adjustable parameter maximum information criterion | |
CN105118332A (en) | Air traffic control analog simulation abnormality detection method and device based on clustering analysis method | |
CN102866630B (en) | Aircraft modeling method based on maximum information amount on-line identification criterion | |
CN106228031A (en) | A kind of flight dynamics model frequency range Adaptive Identification method based on test flight data |
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 |
Granted publication date: 20141126 Termination date: 20211011 |
|
CF01 | Termination of patent right due to non-payment of annual fee |