CN102507131B - Method for acquiring unknown pneumatic parameters in UAV (Unmanned Aerial Vehicle) mathematical model - Google Patents
Method for acquiring unknown pneumatic parameters in UAV (Unmanned Aerial Vehicle) mathematical model Download PDFInfo
- Publication number
- CN102507131B CN102507131B CN2011102944451A CN201110294445A CN102507131B CN 102507131 B CN102507131 B CN 102507131B CN 2011102944451 A CN2011102944451 A CN 2011102944451A CN 201110294445 A CN201110294445 A CN 201110294445A CN 102507131 B CN102507131 B CN 102507131B
- Authority
- CN
- China
- Prior art keywords
- unmanned plane
- control command
- model
- test
- aerodynamic parameter
- 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
Landscapes
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
The invention relates to an improved method for acquiring unknown pneumatic parameters in a UAV (Unmanned Aerial Vehicle) mathematical model, belonging to a UAV model identification technology. The improved method comprises the following steps of: determining a model structure of a UAV according to theoretical derivation; acquiring test flight data by a test flight step; and acquiring the unknown parameters in a parameterized model by utilizing a method combining a predication error algorithm with an improved particle swarm optimization algorithm so as to realize identification work of the model. The improved method disclosed by the invention has the advantages of no need of a wind tunnel test, convenience for a test flight test, simple identification process and high model precision. According to the improved method disclosed by the invention, dependence on test equipment and test conditions is greatly reduced, and convenience and generality of model parameter acquisition are improved.
Description
Technical field
The invention belongs to unmanned plane Model Distinguish technology, relate to the improvement to unknown aerodynamic parameter acquisition methods in the unmanned plane mathematical model.
Background technology
At present, for the problem of obtaining of unknown aerodynamic parameter in the unmanned plane mathematical model, a kind of method is to adopt wind tunnel technique, obtains the pneumatic mathematical model of unmanned plane by blowing, but its experimental cost is higher, and be difficult to guarantee its model accuracy for rotor class unmanned plane; Another kind is to utilize the CIFER software based on the frequency domain identification technology to obtain unknown parameter, but it is had relatively high expectations to experiment in flight test, realize difficulty, and identification process is very complicated, versatility is poor, referring to " model aircraft frequency domain identification method-CIFER algorithm research ", Zou Yu, Pei Hailong, Liu Xin, Zhou Hongbo etc., electric light with control in May, 2010, the 17th the 5th phase of volume.
Summary of the invention
The objective of the invention is: propose a kind of without wind tunnel experiment, be convenient to carry out experiment in flight test, identification process is simple, model accuracy the is high unknown aerodynamic parameter acquisition methods of unmanned plane mathematical model.
Technical scheme of the present invention is: a kind of method of obtaining unknown aerodynamic parameter in the unmanned plane mathematical model, it is characterized in that, and the step of obtaining unknown aerodynamic parameter is as follows:
1, determine the structure of unmanned plane mathematical model: determine the linear parameterization model of this unmanned plane according to the structure type of unmanned plane, comprised unknown aerodynamic parameter in model, structure type is single rotor unmanned helicopter; The linear parameterization model of determining single rotor unmanned helicopter is:
u=[δ
lat?δ
lon?δ
col?δ
ped]
T
Wherein:
X=[u v p q φ θ a
1sb
1sw r r
fb]
t, and [u v w] be respectively the axis speed component, and [φ θ ψ] is respectively rolling, pitching and crab angle, and [p q r] is respectively rolling, pitching and yawrate, [a
1sb
1s] be flapping angle, be respectively rear chamfering and side chamfering, r
fbthe gyrosystem state of feedback, [δ
latδ
lonδ
colδ
ped]
tbe successively the horizontal cyclic pitch control amount of main oar, the vertical cyclic pitch control amount of main oar, main oar always apart from manipulated variable and tailrotorpiston manipulated variable, comprised unknown aerodynamic parameter in matrix A, B;
2, obtain test flight data:
2.1, unmanned plane is taken a flight test, unmanned vehicle is operated under the remote manual control operating pattern, by the aircraft trim;
2.2, obtain the pitch channel test flight data: send the pitch channel control command to unmanned plane, control command is sinusoidal constant amplitude swept-frequency signal, the amplitude of sinusoidal constant amplitude swept-frequency signal represents manipulated variable, the frequency range that frequency sweep covers is 1/3 ω c~3 ω c, the cutoff frequency that ω c is the unmanned plane pitch channel, ω c is provided by unmanned plane designing unit; Unmanned plane is responded the control command of receiving, the Airborne Flight Parameter logging modle records control command and corresponding aircraft response is stored as the pitch channel test flight data;
2.3, obtain course passage test flight data: send course passage control command to unmanned plane, control command is sinusoidal constant amplitude swept-frequency signal, the amplitude of sinusoidal constant amplitude swept-frequency signal represents manipulated variable, the frequency range that frequency sweep covers is 1/3 ω c~3 ω c, the cutoff frequency that ω c is unmanned plane course passage, ω c is provided by unmanned plane designing unit; Unmanned plane is responded the control command of receiving, the Airborne Flight Parameter logging modle records control command and corresponding aircraft response is stored as course passage test flight data;
2.4, obtain the roll channel test flight data: send the roll channel control command to unmanned plane, control command is sinusoidal constant amplitude swept-frequency signal, the amplitude of sinusoidal constant amplitude swept-frequency signal represents manipulated variable, the frequency range that frequency sweep covers is 1/3 ω c~3 ω c, the cutoff frequency that ω c is the unmanned machine rolling passage, ω c is provided by unmanned plane designing unit; Unmanned plane is responded the control command of receiving, the Airborne Flight Parameter logging modle records control command and corresponding aircraft response is stored as the roll channel test flight data;
3, carry out the aerodynamic parameter estimation: adopt recursive prediction error algorithms to carry out according to a preliminary estimate the unknown aerodynamic parameter in the unmanned plane mathematical model, obtain the thick value of aerodynamic parameter;
4, obtain aerodynamic parameter: using the thick value of aerodynamic parameter as the initial value that improves the particle cluster algorithm operation, adopt the improvement particle cluster algorithm to be calculated, obtain final aerodynamic parameter.
Advantage of the present invention is: without wind tunnel experiment, be convenient to carry out experiment in flight test, identification process is simple, and model accuracy is high.The present invention has greatly reduced the dependence to experimental facilities and experiment condition, has improved convenience and versatility that model parameter is obtained.
Embodiment
Below the present invention is described in further details.A kind of method of obtaining unknown aerodynamic parameter in the unmanned plane mathematical model, is characterized in that, the step of obtaining unknown aerodynamic parameter is as follows:
1, determine the structure of unmanned plane mathematical model: the linear parameterization model of determining this unmanned plane according to the structure type of unmanned plane, comprised unknown aerodynamic parameter in model, structure type is single rotor unmanned helicopter, for single rotor unmanned helicopter, adopt ten single order parameterized models, referring to " the layer rank flight control system design of depopulated helicopter " (" Hierarchical Flight Control System Synthesis for Rotorcraft-based Unmanned Aerial Vehicles "), Hyunchul Shim, PhD Thesis, University of California, Berkeley, 2000, for the fixed-wing unmanned plane, adopt nine rank parameterized models, referring to " aerodynamics and flight mechanics ", colleague Liu, publishing house of Beijing Aeronaution College, 1987.
2, obtain test flight data:
2.1, unmanned plane is taken a flight test, unmanned vehicle is operated under the remote manual control operating pattern, by the aircraft trim;
2.2, obtain the pitch channel test flight data: send the pitch channel control command to unmanned plane, control command is sinusoidal constant amplitude swept-frequency signal, the amplitude of sinusoidal constant amplitude swept-frequency signal represents manipulated variable, the frequency range that frequency sweep covers is 1/3 ω c~3 ω c, the cutoff frequency that ω c is the unmanned plane pitch channel, ω c is provided by unmanned plane designing unit; Unmanned plane is responded the control command of receiving, the Airborne Flight Parameter logging modle records control command and corresponding aircraft response is stored as the pitch channel test flight data;
2.3, obtain course passage test flight data: send course passage control command to unmanned plane, control command is sinusoidal constant amplitude swept-frequency signal, the amplitude of sinusoidal constant amplitude swept-frequency signal represents manipulated variable, the frequency range that frequency sweep covers is 1/3 ω c~3 ω c, the cutoff frequency that ω c is unmanned plane course passage, ω c is provided by unmanned plane designing unit; Unmanned plane is responded the control command of receiving, the Airborne Flight Parameter logging modle records control command and corresponding aircraft response is stored as course passage test flight data;
2.4, obtain the roll channel test flight data: send the roll channel control command to unmanned plane, control command is sinusoidal constant amplitude swept-frequency signal, the amplitude of sinusoidal constant amplitude swept-frequency signal represents manipulated variable, the frequency range that frequency sweep covers is 1/3 ω c~3 ω c, the cutoff frequency that ω c is the unmanned machine rolling passage, ω c is provided by unmanned plane designing unit; Unmanned plane is responded the control command of receiving, the Airborne Flight Parameter logging modle records control command and corresponding aircraft response is stored as the roll channel test flight data;
3, carry out the aerodynamic parameter estimation: adopt recursive prediction error algorithms to carry out according to a preliminary estimate the unknown aerodynamic parameter in the unmanned plane mathematical model, obtain the thick value of aerodynamic parameter, referring to " research of small-sized depopulated helicopter physical parameters identification problem ", Yan Chao, South China Science & Engineering University's master thesis, 2004;
4, obtain aerodynamic parameter: using the thick value of aerodynamic parameter as the initial value that improves the particle cluster algorithm operation, adopt the improvement particle cluster algorithm to be calculated, obtain final aerodynamic parameter, referring to " parameter identification based on the APSO algorithm and optimization ", Wu Yanxiang, Li Xiaobin, Sun Haiyan etc., science and technology and engineering in July, 2008, the 8th the 14th phase of volume.
Principle of work of the present invention is: at first according to theory, derive and determine the model structure without mankind's aircraft, then gather test flight data by the step of taking a flight test, the method of utilizing afterwards recursive prediction error algorithms to combine with the improvement particle cluster algorithm, unknown parameter in getting parms model, realize the identification work to model.
Embodiment 1
Single rotor unmanned helicopter has carried the flight data recording unit, adopts said method to complete obtaining of unknown aerodynamic parameter in the depopulated helicopter mathematical model.
1, the linear parameterization model of determining single rotor unmanned helicopter is:
u=[δ
lat?δ
lon?δ
col?δ
ped]
T
Wherein:
X=[u v p q φ θ a
1sb
1sw r r
fb]
t, and [u v w] be respectively the axis speed component, and [φ θ ψ] is respectively rolling, pitching and crab angle, and [p q r] is respectively rolling, pitching and yawrate, [a
1sb
1s] be flapping angle, be respectively rear chamfering and side chamfering, r
fbthe gyrosystem state of feedback, [δ
latδ
lonδ
colδ
ped]
tbe successively the horizontal cyclic pitch control amount of main oar, the vertical cyclic pitch control amount of main oar, main oar always apart from manipulated variable and tailrotorpiston manipulated variable, comprised unknown aerodynamic parameter in matrix A, B;
2, obtain test flight data:
2.1, single rotor unmanned helicopter is taken a flight test, it is operated under the remote manual control operating pattern, and by the aircraft trim;
2.2, the cutoff frequency of single rotor unmanned helicopter pitching, course and roll channel is 3~4rad/s, respectively depopulated helicopter pitching, course and roll channel are sent to the roll channel control command, control command is sinusoidal constant amplitude swept-frequency signal, and the frequency range that frequency sweep covers is 1~12rad/s;
3, adopt recursive prediction error algorithms to carry out according to a preliminary estimate the unknown aerodynamic parameter in the unmanned plane mathematical model, obtain the thick value of aerodynamic parameter;
4, using the thick value of aerodynamic parameter as the initial value that improves the particle cluster algorithm operation, adopt the improvement particle cluster algorithm to be calculated, obtain final aerodynamic parameter:
x=[u(m/s)?v(m/s)?p(rad/s)?q(rad/s)?φ(rad)?θ(rad)?a
1s(rad)?b
1s(rad)?w(m/s)?r(rad/s)?r
fb]
T
u=[δ
lat(V)?δ
lon(V)?δ
col(V)?δ
ped(V)]
T
With regard to the unknown parameter that has completed single rotor unmanned helicopter aerodynamic model, obtain like this.
Claims (1)
1. a method of obtaining unknown aerodynamic parameter in the unmanned plane mathematical model, is characterized in that, the step of obtaining unknown aerodynamic parameter is as follows:
1.1, determine the structure of unmanned plane mathematical model: determine the linear parameterization model of this unmanned plane according to the structure type of unmanned plane, comprised unknown aerodynamic parameter in model, structure type is single rotor unmanned helicopter; The linear parameterization model of determining single rotor unmanned helicopter is:
u=[δ
lat?δ
lon?δ
col?δ
ped]
T
Wherein:
X=[u v p q φ θ a
1sb
1sw r r
fb]
t, and [u v w] be respectively the axis speed component, and [φ θ ψ] is respectively rolling, pitching and crab angle, and [p q r] is respectively rolling, pitching and yawrate, [a
1sb
1s] be flapping angle, be respectively rear chamfering and side chamfering, r
fbthe gyrosystem state of feedback, [δ
latδ
lonδ
colδ
ped]
tbe successively the horizontal cyclic pitch control amount of main oar, the vertical cyclic pitch control amount of main oar, main oar always apart from manipulated variable and tailrotorpiston manipulated variable, comprised unknown aerodynamic parameter in matrix A, B;
1.2, obtain test flight data:
1.2.1, unmanned plane is taken a flight test, unmanned vehicle is operated under the remote manual control operating pattern, by the aircraft trim;
1.2.2, obtain the pitch channel test flight data: send the pitch channel control command to unmanned plane, control command is sinusoidal constant amplitude swept-frequency signal, the amplitude of sinusoidal constant amplitude swept-frequency signal represents manipulated variable, the frequency range that frequency sweep covers is 1/3 ω c~3 ω c, the cutoff frequency that ω c is the unmanned plane pitch channel, ω c is provided by unmanned plane designing unit; Unmanned plane is responded the control command of receiving, the Airborne Flight Parameter logging modle records control command and corresponding aircraft response is stored as the pitch channel test flight data;
1.2.3, obtain course passage test flight data: send course passage control command to unmanned plane, control command is sinusoidal constant amplitude swept-frequency signal, the amplitude of sinusoidal constant amplitude swept-frequency signal represents manipulated variable, the frequency range that frequency sweep covers is 1/3 ω c~3 ω c, the cutoff frequency that ω c is unmanned plane course passage, ω c is provided by unmanned plane designing unit; Unmanned plane is responded the control command of receiving, the Airborne Flight Parameter logging modle records control command and corresponding aircraft response is stored as course passage test flight data;
1.2.4, obtain the roll channel test flight data: send the roll channel control command to unmanned plane, control command is sinusoidal constant amplitude swept-frequency signal, the amplitude of sinusoidal constant amplitude swept-frequency signal represents manipulated variable, the frequency range that frequency sweep covers is 1/3 ω c~3 ω c, the cutoff frequency that ω c is the unmanned machine rolling passage, ω c is provided by unmanned plane designing unit; Unmanned plane is responded the control command of receiving, the Airborne Flight Parameter logging modle records control command and corresponding aircraft response is stored as the roll channel test flight data;
1.3, carry out the aerodynamic parameter estimation: adopt recursive prediction error algorithms to carry out according to a preliminary estimate the unknown aerodynamic parameter in the unmanned plane mathematical model, obtain the thick value of aerodynamic parameter;
1.4, obtain aerodynamic parameter: using the thick value of aerodynamic parameter as the initial value that improves the particle cluster algorithm operation, adopt and improve particle cluster algorithm and calculated, obtain final aerodynamic parameter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102944451A CN102507131B (en) | 2011-09-26 | 2011-09-26 | Method for acquiring unknown pneumatic parameters in UAV (Unmanned Aerial Vehicle) mathematical model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102944451A CN102507131B (en) | 2011-09-26 | 2011-09-26 | Method for acquiring unknown pneumatic parameters in UAV (Unmanned Aerial Vehicle) mathematical model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102507131A CN102507131A (en) | 2012-06-20 |
CN102507131B true CN102507131B (en) | 2013-12-11 |
Family
ID=46219238
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011102944451A Active CN102507131B (en) | 2011-09-26 | 2011-09-26 | Method for acquiring unknown pneumatic parameters in UAV (Unmanned Aerial Vehicle) mathematical model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102507131B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9652904B2 (en) | 2014-09-30 | 2017-05-16 | SZ DJI Technology Co., Ltd. | System and method for data recording and analysis |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102799187B (en) * | 2012-06-26 | 2014-09-17 | 中国航空工业第六一八研究所 | Spraying breakpoint continuing method for unmanned helicopter |
CN102789527A (en) * | 2012-07-11 | 2012-11-21 | 南京航空航天大学 | Particle swarm optimization method for airplane trim |
CN104102127B (en) * | 2014-07-17 | 2016-06-01 | 厦门大学 | A kind of machine load Aerodynamic Parameter Identification system |
CN106650095B (en) * | 2016-12-21 | 2020-07-14 | 中国航天空气动力技术研究院 | Correction method of unmanned aerial vehicle control matrix based on wind tunnel test data and CFD calculation |
CN110704948A (en) * | 2019-09-24 | 2020-01-17 | 江西慧识智能科技有限公司 | Design method of intelligent controller of unmanned aerial vehicle |
CN111538237B (en) * | 2020-03-20 | 2021-10-29 | 北京航空航天大学 | Method for identifying and correcting non-linear light gray model of tilt rotor unmanned aerial vehicle |
CN113029512B (en) * | 2021-04-01 | 2022-08-05 | 中国空气动力研究与发展中心高速空气动力研究所 | Air flow field model correction control method for temporary transonic wind tunnel |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101286071A (en) * | 2008-04-24 | 2008-10-15 | 北京航空航天大学 | Multiple no-manned plane three-dimensional formation reconfiguration method based on particle swarm optimization and genetic algorithm |
CN101321667A (en) * | 2005-12-06 | 2008-12-10 | 空中客车德国有限公司 | Method for reconstructing gusts and structural loads at aircraft, in particular passenger aircraft |
US20090157364A1 (en) * | 2007-12-18 | 2009-06-18 | Airbus Espana, S.L.. | Method and system for a quick calculation of aerodynamic forces on an aircraft |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3691023B2 (en) * | 2002-03-29 | 2005-08-31 | 三菱重工業株式会社 | Model support device control method |
-
2011
- 2011-09-26 CN CN2011102944451A patent/CN102507131B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101321667A (en) * | 2005-12-06 | 2008-12-10 | 空中客车德国有限公司 | Method for reconstructing gusts and structural loads at aircraft, in particular passenger aircraft |
US20090157364A1 (en) * | 2007-12-18 | 2009-06-18 | Airbus Espana, S.L.. | Method and system for a quick calculation of aerodynamic forces on an aircraft |
CN101286071A (en) * | 2008-04-24 | 2008-10-15 | 北京航空航天大学 | Multiple no-manned plane three-dimensional formation reconfiguration method based on particle swarm optimization and genetic algorithm |
Non-Patent Citations (5)
Title |
---|
《小型无人直升机物理参数辨别算法研究》;李德根等;《福建电脑》;20070430(第4期);第96、106页 * |
《某型无人直升机飞控系统仿真研究》;陈巍等;《飞行力学》;20040331;第22卷(第1期);第41-44页 * |
JP特开2003-294575A 2003.10.15 |
李德根等.《小型无人直升机物理参数辨别算法研究》.《福建电脑》.2007,(第4期), |
陈巍等.《某型无人直升机飞控系统仿真研究》.《飞行力学》.2004,第22卷(第1期), |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9652904B2 (en) | 2014-09-30 | 2017-05-16 | SZ DJI Technology Co., Ltd. | System and method for data recording and analysis |
US9905060B2 (en) | 2014-09-30 | 2018-02-27 | SZ DJI Technology Co., Ltd. | System and method for data recording and analysis |
Also Published As
Publication number | Publication date |
---|---|
CN102507131A (en) | 2012-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102507131B (en) | Method for acquiring unknown pneumatic parameters in UAV (Unmanned Aerial Vehicle) mathematical model | |
Madani et al. | Control of a quadrotor mini-helicopter via full state backstepping technique | |
CN111873804B (en) | Electric automobile working condition data generation method and energy consumption evaluation working condition generation method | |
CN106840097B (en) | Road slope estimation method based on adaptive extended Kalman filtering | |
Madani et al. | Sliding mode observer and backstepping control for a quadrotor unmanned aerial vehicles | |
CN106844887B (en) | Dynamics modeling method and device for rotor unmanned aerial vehicle | |
Zhang et al. | Real-time estimation of vehicle mass and road grade based on multi-sensor data fusion | |
CN103246203B (en) | A kind of microminiature four rotor wing unmanned aerial vehicle speed state Forecasting Methodology based on GPS | |
CN104698845B (en) | Radar antenna servo tracking method and system based on automatic disturbance rejection controller | |
González et al. | Real-time attitude stabilization of a mini-uav quad-rotor using motor speed feedback | |
CN103217981A (en) | Four-rotor aircraft speed control method based on integral variable structure control | |
CN104391497A (en) | On-line joint servo system parameter identification and controller parameter optimization method | |
CN105680736A (en) | Yaw angle-based double-motor speed synchronization and balance control method | |
CN110007599B (en) | Self-adaptive integral sliding mode control method and system of laser radar scanning mechanism | |
CN112068594A (en) | JAYA algorithm optimization-based course control method for small unmanned helicopter | |
CN109871044B (en) | Rotating speed tracking method and device | |
CN109635376B (en) | Modeling method and system based on unmanned aerial vehicle | |
US20220381649A1 (en) | Vehicle data communication and interaction system and construction method thereof | |
Gelbal et al. | A unified architecture for scalable and replicable autonomous shuttles in a smart city | |
CN111879287B (en) | Forward terrain three-dimensional construction method of low-speed vehicle based on multiple sensors | |
CN110031646B (en) | Method for correcting driving speed by using driving direction vibration acceleration | |
CN105631179A (en) | Carrier transient response indication method based on star satellite load couple analysis result | |
CN114608841A (en) | Intelligent networking automobile test platform | |
Luo et al. | High precision positioning system for autopilot based on multi-sensor fusion | |
Nong et al. | System identification of a small unmanned aerial vehicle based on time and frequency domain technologies |
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 |