CN106997208B - A kind of control method towards the hypersonic aircraft under condition of uncertainty - Google Patents
A kind of control method towards the hypersonic aircraft under condition of uncertainty Download PDFInfo
- Publication number
- CN106997208B CN106997208B CN201710325189.5A CN201710325189A CN106997208B CN 106997208 B CN106997208 B CN 106997208B CN 201710325189 A CN201710325189 A CN 201710325189A CN 106997208 B CN106997208 B CN 106997208B
- Authority
- CN
- China
- Prior art keywords
- hypersonic aircraft
- follows
- parameter
- virtual controlling
- formula
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000013461 design Methods 0.000 claims abstract description 36
- 238000009415 formwork Methods 0.000 claims abstract description 20
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 230000008569 process Effects 0.000 claims abstract description 6
- 238000012360 testing method Methods 0.000 claims abstract description 4
- 238000009795 derivation Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000012938 design process Methods 0.000 claims description 2
- 235000021170 buffet Nutrition 0.000 abstract description 6
- 210000001258 synovial membrane Anatomy 0.000 abstract description 3
- 230000004044 response Effects 0.000 description 22
- 230000001276 controlling effect Effects 0.000 description 17
- 238000013459 approach Methods 0.000 description 5
- 230000008878 coupling Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 244000145845 chattering Species 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses a kind of control methods towards the hypersonic aircraft under condition of uncertainty, the following steps are included: step S1, according to hypersonic aircraft longitudinal dynamics equation, overall-finished housing linearization process is used to flying speed V and flying height h, obtains its corresponding state equation；Step S2 designs contragradience sliding mode controller according to this state equation, selects double power Reaching Laws in synovial membrane face are as follows:In formula, k1> 0, k21,1 > λ > 0 of > 0, η > is parameter；Step S3 builds test platform based on this contragradience sliding mode controller, carries out performance simulation.The method of the present invention uses double power sliding mode controllers based on Backstepping, and this method, which can be effectively reduced handoff gain and eliminate, buffets problem present in sliding formwork control.
Description
Technical field
The present invention relates to hypersonic aircraft technical fields, and in particular to a kind of high ultrasound towards under condition of uncertainty
The flight control method of fast aircraft.
Background technique
Hypersonic aircraft mainly executes aerial mission, near space in the near space region apart from 20~100km of ground
The mankind space that develops and uses very well not yet so far, in terms of great potential, it is hypersonic to fly
Concealment of the row device near space region is very high and aircraft communication signal in the field is very strong, hypersonic to fly
Row device, which determines it not only near space flight, has the advantage of aeronautical technology, but also has the advantages of spacecraft can not possess.In army
For thing using upper, it is very fast that hypersonic aircraft executes speed, can first time with reaching target, will cause quickly accurate
Hit the arriving in weapon epoch.Especially information acquisition, hasty breaching, communication support, long-range strike, in terms of it is great
Development potentiality.Just because of these huge military values and potential economic value, more and more military powers are all competitively
Develop hypersonic aircraft.
Due to changeable, complicated flight environment of vehicle and the extraneous strong jamming of the self structure of hypersonic aircraft, it is necessary to
Design has good stability, fast response time and control controller with high accuracy, this just proposes advanced control method urgent
The demand cut.Expert applies to many advanced control technologies in hypersonic aircraft both at home and abroad at present, and achieves very
Good control effect.Wherein several important and successful control method is introduced below.
The main thought of gain preset control method is by complicated nonlinear Control PROBLEM DECOMPOSITION into multiple linear models
And the design problem of multiple linear controllers.The method is widely used in Design of Flight Control, and U.S. X-43A flies
Row device introduces gain preset method in Control System Design, but the flight time of aircraft is too short, under high maneuver state,
The validity of this method can not be proved to.Under High Angle of Attack and strong maneuvering condition, the flight of hypersonic aircraft has strong
Non-linear and high coupling is much unable to satisfy the requirement of performance indicator just with the controller that gain preset method designs.
It is also one of most widely used method that control theory of feedback linearization, which is most important in nonlinear control method, real
Two existing effective ways are as follows: Differential Geometry method and dynamic inversion.They are different from the part of traditional Taylor expansion
Linearisation, but nonlinear system is subjected to exact linearization method, wherein the higher order term of Taylor expansion is contained, but because this
Kind method needs to establish on the basis of accurate model, and very sensitive to error existing for model, to reduce entire non-
The robustness of linear system.
Sliding mode variable structure control method (SMVSC) is an important method for handling nonlinear system, so-called structure changes sheet
Refer to the discontinuous non-linear shear that the feedback controller structure of internal system is occurred in matter.When system mode passes through not same district
When domain, the structure of feedback control requires the switch logic formulated to become by designer according to a set of according to system performance index
Change, enables control system that there is certain adaptation to factors such as inherent Parameters variation and the external environment disturbances of controlled device
Power guarantees that system performance reaches desired performance indicator requirement.But system mode easily causes to buffet in motion switch.For
Overcome these defects, many domestic and foreign scholars are proposed some relatively effective methods, such as saturation function method, Reaching Law
Method, boundary layer method, High-Order Sliding Mode method etc..
Since hypersonic aircraft using body coupled structure and flies under conditions of High aititude and big Mach number,
Cause it quite sensitive to the variation of gas condition, and hypersonic aircraft also receives structure dynamics in flight course,
The influence for promoting dynamic and coupling between them, pneumatic and propulsion characteristic is uncertain, or even is difficult to estimate.These because
The influence of element so that the model of hypersonic aircraft be it is uncertain, it is changeable, it is unstable, and there are input and output
Between close coupling, therefore design the controller with nonlinearity and strong robustness for hypersonic aircraft and become especially
It is important.Meanwhile there are convergence rates slowly and the deficiencies of buffet for Reaching Law traditional in sliding mode controller design in the prior art,
The present invention provides a kind of new resolving ideas.
Summary of the invention
It is an object of the invention to overcome deficiency in the prior art, provide a kind of towards superb under condition of uncertainty
The control method of velocity of sound aircraft, using double power sliding mode controllers based on Backstepping, this method can be effectively reduced and cut
It changes gain and eliminates and buffet problem present in sliding formwork control.
In order to solve the above technical problems, the present invention provides a kind of hypersonic aircrafts towards under condition of uncertainty
Control method, characterized in that the following steps are included:
Step S1, according to hypersonic aircraft longitudinal dynamics equation, to flying speed V and flying height h using complete
State feedback linearization processing, obtains its corresponding state equation；
Step S2 designs contragradience sliding mode controller according to this state equation,
Define sliding-mode surface:
Wherein, ci> 0 is parameter, eiAnd e3For system tracking error.
Select double power Reaching Laws in synovial membrane face are as follows:
In formula, k1> 0, k21,1 > λ > 0 of > 0, η > is parameter；
Contragradience sliding formwork control ratio are as follows:
Step S3 builds test platform based on this contragradience sliding mode controller, carries out performance simulation.
Further, in step S1, the dynamical equation of hypersonic aircraft longitudinal direction model can be converted into following state
Equation form:
WhereinU=[βc
δe]T, V is flying speed, h is height.
Further, in step s 2, the design process of contragradience sliding mode controller are as follows: be decomposed into speed different three and return
Road: fast circuit, slower circuit, slow circuit；Virtual controlling rule is calculated since the slow circuit subsystem farthest from control input, it is first
First define the error of subsystems:
WhereinFor the desired command signal of tracking, xidIt is restrained for the virtual controlling of subsystem；
First subsystem, that is, slow circuit subsystem tracking error is e1=x1-xd, after both members difference derivationAgain because of e2=x2-x2dIt substitutes into:
x2dIt is the virtual controlling rule of slower circuit subsystem for second subsystem, design virtual controlling restrains x2dAre as follows:Wherein, kx1> 0 is that virtual controlling restrains design parameter to be asked, and is obtained after being updated to formula (11):
Second sub- system tracking error e2=x2-x2dBoth members obtain after distinguishing derivationThen again because of e3
=x3-x3d, after substitution:
x3dIt is to return to the virtual controlling rule of subsystems for three subsystems, by formula (13) design virtual controlling rule
x3dAre as follows:Wherein, kx2> 0 is that virtual controlling restrains design parameter to be asked.After being updated to formula (13)
?
Three subsystems return to subsystems tracking error e3=x3-x3dBoth members obtain after distinguishing derivation:
Design control input:
kx3> 0 is controller design parameter, is obtained:
Further, the calculating process of the stability of contragradience sliding mode controller are as follows:
Define Liapunov function are as follows:
The derivative of the function against time has:
Formula (24) and (25) are substituted into (39) to obtain:
Because of coefficient kxi> 0, k1> 0, k2> 0, obtains:
In conclusion system mode can reach diverter surface within the limited time, meet system stable condition.
Further, aircraft altitude is chosen when performance simulation and the initial equilibrium conditions of flying speed are respectively
The simulation parameter of nominal parameters under 33528m and 4590.3m/s model, double power sliding mode controllers is chosen are as follows:
Control parameter: kx1=3, kx2=3, c1=2, c2=2, k1=0.8, k2=3, η=0.5, λ=1.5.
Further, also this maximum Parameter Perturbation is introduced to dummy vehicle when emulation to verify.
Compared with prior art, the beneficial effects obtained by the present invention are as follows being: using double power sliding formwork controls based on Backstepping
Device processed, which, which can be effectively reduced handoff gain and eliminate, buffets problem present in sliding formwork control.Based on contragradience
Double power sliding mode controllers of method, which combine the uncertain problem that sliding formwork control meets system under matching condition, to be had relatively by force
Robustness and Backstepping to the advantage of processing system mismatched uncertainties, be capable of rapidity and the Shandong of Guarantee control system
Stick.
Detailed description of the invention
Fig. 1 is the model framework chart that hypersonic aircraft longitudinal direction model feedback linearizes；
Fig. 2 is the simulation result using contragradience sliding formwork control ratio of the present invention: wherein (a) is the tracking response curve pair of speed
Than figure, (b) be height tracking response curve comparison figure；
Fig. 3 is the comparison result using contragradience sliding formwork control ratio of the present invention under different parameters perturbation: wherein (a) is flight
The tracking response Dependence Results of speed；(b) be flying height tracking response Dependence Results；(c) be flight angular speed tracking
Response curve result；It (d) is the tracking response Dependence Results of flying angle；It (e) is the tracking response curve knot of flight angle of rudder reflection
Fruit；(f) be engine regulating valve tracking response Dependence Results.
Specific embodiment
The invention will be further described below in conjunction with the accompanying drawings.Following embodiment is only used for clearly illustrating the present invention
Technical solution, and not intended to limit the protection scope of the present invention.
When hypersonic flight cruising flight, it is assumed that roll angle and sideslip are all zero, then approximation ignore it is longitudinal and horizontal
Lateral coupling carries out decoupling processing to hypersonic aircraft longitudinal direction model, and it is vertical to obtain hypersonic aircraft rigid body
To kinetic model.The control input of aircraft is engine throttle opening degree instruction βcWith elevator angle δe, take aircraft
Quantity of state have: flying speed V, track inclination angle γ, angle of attack α, pitch rate q and flying height h, then hypersonic flight
Device Longitudinal Dynamic Model are as follows:
Wherein, L, D, T respectively indicate aircraft lift, resistance and motor power, and M indicates pitching moment, IyyIndicate vertical
To rotary inertia, m is quality, and R indicates that earth radius, r indicate the earth's core of aircraft away from r=R+h, g are terrestrial gravitation acceleration
Constant.
The purpose of hypersonic aircraft control is selection engine's throttling valve opening control amount βcWith lifting angle of rudder reflection δe,
Guarantee that the flying speed V and flying height h of aircraft be capable of fast accurate traces into designated value VdAnd hd.Formula (2)~(3) are right
It exports flying speed V and flying height h and uses overall-finished housing linearization process, i.e., flying speed V and flying height h is carried out
Differential process, until control inputs βcOr δeIt appears in differential formula.Thus according to kinetics equation to output flying speed V
With flying height h using the expression formula after overall-finished housing linearization process are as follows:
In formula:
b11=(ρ V2swcω2/2m)cosα (5)
b21=(ρ V2swcω2/2m)sin(α+γ)(7)
In formula: x0Indicate the initial value of quantity of state, DαIt indicates because of resistance caused by the angle of attack, LαCaused by indicating because of the angle of attack
Lift, TαThrust caused by indicating because of the angle of attack, ρ represent atmospheric density, SwAircraft area of reference is represented, c represents Average aerodynamic string
It is long, ceFor elevator coefficient, it is 33528m that other coefficient values, which are in nominal height h, and datum speed v is 4590.3m/s flight shape
It is obtained under state.It can be obtained by above formulah(4)Expression formula contain control input βcAnd δe。
A kind of control method towards the hypersonic aircraft under condition of uncertainty of the invention, including following procedure:
According to above-mentioned coordinate transform, the dynamical equation (2) of hypersonic aircraft longitudinal direction model can be converted into following shape
State equation form:
WhereinU=[βc
δe]T。
Fig. 1 provides the longitudinal model framework chart obtained after the linearisation of hypersonic aircraft longitudinal direction model feedback, according to
The characteristics of state variable is broken down into three different circuits of speed: fast circuit, slower circuit, slow circuit respectively correspond three
A subsystem returns to subsystems, slower circuit subsystem, slow circuit subsystem.Only consider vertical passage, flying speed V,
Height h is slow state.
With reference to the contragradience sliding formwork control of hypersonic aircraft track following " control ", according to the design philosophy of Backstepping,
The design that controller is carried out since the slow circuit subsystem farthest from control input u is design that virtual controlling is restrained, gradually after
It moves back.The error of subsystems is defined first:
WhereinFor the desired command signal of tracking, xidIt is restrained for the virtual controlling of subsystem.
First subsystem, that is, slow circuit subsystem tracking error is e1=x1-xd, control target is e1→ 0, both members
Respectively after derivationAgain because of e2=x2-x2dIt substitutes into:
x2dIt is the virtual controlling rule of slower circuit subsystem for second subsystem, design virtual controlling restrains x2dAre as follows:Wherein, kx1> 0 is that virtual controlling restrains design parameter to be asked, and is obtained after being updated to formula (11):
Second sub- system tracking error e2=x2-x2d, in order to make error e2Minimum, both members obtain after distinguishing derivationThen again because of e3=x3-x3d, after substitution:
x3dIt is to return to the virtual controlling rule of subsystems for three subsystems, by formula (13) design virtual controlling rule
x3dAre as follows:Wherein, kx2> 0 is that virtual controlling restrains design parameter to be asked.After being updated to formula (13)
?
Three subsystems return to subsystems tracking error e3=x3-x3dBoth members obtain after distinguishing derivation:
Design control input:
kx3> 0 is controller design parameter, is obtained:
What sliding formwork control solved is matching uncertain problem, and Backstepping solves the problems, such as mismatched uncertainties, so
Sliding mode design, final design contragradience sliding formwork control ratio have been carried out after design Backstepping.
Define sliding-mode surface:
Wherein ci> 0, for design parameter to be asked.
For Reaching Law traditional in sliding formwork control, there are convergence rates slowly and the deficiencies of buffet, in order to realize that finite time arrives
Up to sliding-mode surface, and weaken chattering phenomenon, select double power Reaching Laws in synovial membrane face are as follows:
In formula, k1> 0, k21,1 > λ > 0 of > 0, η > is parameter to be asked.
When | S | > 1 indicates that system mode far from sliding mode, that is, only has the first item in (19) to play a leading role；When | s |
< 1 indicates system mode close to sliding formwork state, and only Section 2 plays a major role in formula (19), sufficiently combines this two excellent
Gesture, so that system mode has better motion qualities.
In conjunction with reachable condition, the sliding moding structure flight control system based on double power Reaching Laws converges to for zero time
Are as follows:
Wherein t1,t2Respectively indicate convergence time to be asked.
It proves: when system mode is far from sliding-mode surface, because of 0 < λ < 1, η > 1, so velocity of approach is mainly by first item
It determines, does not consider the influence of Section 2 at this time, formula (19) can write a Chinese character in simplified form are as follows:
Above formula both sides integral can be obtained:
S1-η=-(1- η) k1t+S(0)1-η (22)
Therefore the time required to sliding-mode surface S=0 → S=1 can be obtained are as follows:
When system mode moves closer to sliding-mode surface, because of 0 < λ < 1, η > 1, velocity of approach is mainly determined by Section 2,
Therefore the influence of first item and Parameter Perturbation is not considered, and formula (19) can write a Chinese character in simplified form are as follows:
Formula (32) both sides are integrated:
S1-λ=-(1- λ) k2t+1 (25)
Therefore the time required to sliding-mode surface S=1 → S=0 can be obtained are as follows:
Therefore, convergence time is sum of the two, it may be assumed that
From the above analysis, work as S=0, because speed is gradually decreased as when system mode reaches sliding mode
Zero, with sliding mode realize it is smooth excessively, largely reduced system chatter.As long as suitably increasing k1It can add with the value of η
Fast velocity of approach of the quantity of state far from Fault slip rate, suitably increase k2It can accelerate quantity of state close to Fault slip rate with λ
Velocity of approach.Theoretical analysis shows that: double power sliding formworks can effectively eliminate buffeting, and when far from and close to sliding mode
All there is cracking velocity of approach, with the Second Order Sliding Mode in High-Order Sliding Mode, there is similar convergence property, i.e. S=0.
Derivation is carried out to (18) to obtain:
Formula (12), (14) and (15) are substituted into (27), and (19) and (27) are combined to obtained contragradience sliding formwork control ratio
Are as follows:
Stability analysis
Define Liapunov function are as follows:
The derivative of the function against time has:
Formula (24) and (25) are substituted into (39) to obtain:
Because of coefficient kxi> 0, k1> 0, k2> 0, obtains:
In conclusion system mode can reach diverter surface within the limited time, meet system stable condition.
The present invention solves chattering phenomenon generally existing in sliding formwork control using the control method of double power sliding formworks first,
Greatly enhance the accuracy of the control method.
Embodiment
In order to verify the control effect of designed contragradience sliding mode controller, emulation point is carried out to hypersonic aircraft
Analysis gives instruction trace signal, while adjusted design controller parameter, obtains corresponding instruction trace effect, while being verifying
The robustness of contragradience sliding-mode control introduces this maximum Parameter Perturbation to dummy vehicle and verifies, shows the controlling party
Method greatly improves the stability and accuracy of system.
Simulating, verifying is carried out to the hypersonic vehicle having built up from the angle of emulation below, it is imitative in Matlab
In true test, referring to simulating, verifying condition in " the contragradience sliding formwork control of hypersonic aircraft track following " document, chooses and fly
The initial equilibrium conditions of row device flying height and flying speed are respectively the nominal parameters under 33528m and 4590.3m/s model,
The simulation parameter of double power sliding mode controllers is chosen are as follows:
Control parameter: kx1=3, kx2=3, c1=2, c2=2, k1=0.8, k2=3, η=0.5, λ=1.5.
Fig. 2 is the response curve of the speed and height under the conditions of contragradience sliding formwork control ratio, wherein (a) is the tracking of speed
Response curve comparison diagram, (b) be height tracking response curve comparison figure.It can be seen from Fig. 2 that contragradience sliding-mode control is a kind of
Very effective nonlinear control method, speed V and height the h output valve of hypersonic aircraft can preferably trace commands
Signal has lesser overshoot, and tenacious tracking may be implemented in 15s, and system has preferable tracking performance.
Maximum Parameter Perturbation amount, the maximum perturbation parameter that Selection Model allows are introduced in the model parameter nominal value of foundation
It is worth as follows: | Δ m |/m0=0.03, | Δ ce|/ce0=0.02, | Δ sw|/sw0=0.03, | Δ ρ |/ρ0=0.03, wherein subscript
The nominal value of " 0 " expression relevant parameter.
Response curve with Parameter Perturbation is compared with the model response curve that Parameter Perturbation is not added, i.e. Fig. 3
The response curve of each state under Parameter Perturbation, wherein (a) is the tracking response Dependence Results of flying speed；It (b) is flying height
Tracking response Dependence Results；(c) be flight angular speed tracking response Dependence Results；It (d) is the tracking response of flying angle
Dependence Results；It (e) is the tracking response Dependence Results of flight angle of rudder reflection；(f) be engine regulating valve tracking response curve knot
Fruit.Positive and negative variation is gone respectively to each variable parameter in Fig. 3, when to take positive Parameter Perturbation be Parameter Perturbation+100%, negative when taking
Parameter Perturbation be -100%.
From figure 3, it can be seen that when Parameter Perturbation reaches maximum, there is lesser overshoot in speed tracing, and height with
There is subtle steady-state error in track and angle of attack response, and pitch rate and change in angle of attack are small, and whole control effect is good, instead
Step sliding-mode control has good compensating action for Parameter Perturbation and external interference, and whole system has preferable tracing property
It can be with stronger robust performance.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, without departing from the technical principles of the invention, several improvements and modifications, these improvements and modifications can also be made
Also it should be regarded as protection scope of the present invention.
Claims (6)
1. a kind of control method towards the hypersonic aircraft under condition of uncertainty, characterized in that the following steps are included:
Step S1 uses total state to flying speed V and flying height h according to hypersonic aircraft longitudinal dynamics equation
Feedback linearization processing, obtains its corresponding state equation；
Step S2 designs contragradience sliding mode controller according to this state equation,
Define sliding-mode surface:
Wherein, ci> 0 is parameter, eiAnd e3For system tracking error,
Select double power Reaching Laws of sliding-mode surface are as follows:
In formula, k1> 0, k21,1 > λ > 0 of > 0, η > is parameter；
Contragradience sliding formwork control ratio are as follows:
In formula,x3dIt is to return to the virtual controlling rule of subsystems for three subsystems；
Step S3 builds test platform based on this contragradience sliding mode controller, carries out performance simulation.
2. a kind of control method towards the hypersonic aircraft under condition of uncertainty according to claim 1, special
Sign is, in step S1, the dynamical equation of hypersonic aircraft longitudinal direction model is converted into following state equation form:
WhereinU=[βc δe]T, V
It is height for flying speed, h.
3. a kind of control method towards the hypersonic aircraft under condition of uncertainty according to claim 1, special
Sign is, in step s 2, the design process of contragradience sliding mode controller are as follows: be decomposed into three different circuits of speed: fast circuit, compared with
Slow circuit, slow circuit；Virtual controlling rule is calculated since the slow circuit subsystem farthest from control input, defines each height first
The error of system:
WhereinFor the desired command signal of tracking, xidIt is restrained for the virtual controlling of subsystem；
First subsystem, that is, slow circuit subsystem tracking error is e1=x1-xd, after both members difference derivationAgain because of e2=x2-x2dIt substitutes into:
x2dIt is the virtual controlling rule of slower circuit subsystem for second subsystem, design virtual controlling restrains x2dAre as follows:Wherein, kx1> 0 is that virtual controlling restrains design parameter to be asked, and is obtained after being updated to formula (11):
Second sub- system tracking error e2=x2-x2dBoth members obtain after distinguishing derivationThen again because of e3=x3-
x3d, after substitution:
x3dIt is to return to the virtual controlling rule of subsystems for three subsystems, x is restrained by formula (13) design virtual controlling3dAre as follows:Wherein, kx2> 0 is that virtual controlling restrains design parameter to be asked, and is obtained after being updated to formula (13)
Three subsystems return to subsystems tracking error e3=x3-x3dBoth members obtain after distinguishing derivation:
Design control input:
kx3> 0 is controller design parameter, is obtained:
4. a kind of control method towards the hypersonic aircraft under condition of uncertainty according to claim 1, special
Sign is the calculating process of the stability of contragradience sliding mode controller are as follows:
Define Liapunov function are as follows:
The derivative of the function against time has:
Formula (24) and (25) are substituted into (39) to obtain:
Because of coefficient kxi> 0, k1> 0, k2> 0, obtains:
In conclusion system mode can reach diverter surface within the limited time, meet system stable condition.
5. a kind of control method towards the hypersonic aircraft under condition of uncertainty according to claim 1, special
Sign is, when performance simulation choose aircraft altitude and flying speed initial equilibrium conditions be respectively 33528m and
The simulation parameter of nominal parameters under 4590.3m/s model, double power sliding mode controllers is chosen are as follows: kx1=3, kx2=3, c1=2,
c2=2, k1=0.8, k2=3, η=0.5, λ=1.5.
6. a kind of control method towards the hypersonic aircraft under condition of uncertainty according to claim 1, special
Sign is, when emulation also introduces this maximum Parameter Perturbation to dummy vehicle and verifies.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710325189.5A CN106997208B (en) | 2017-05-10 | 2017-05-10 | A kind of control method towards the hypersonic aircraft under condition of uncertainty |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710325189.5A CN106997208B (en) | 2017-05-10 | 2017-05-10 | A kind of control method towards the hypersonic aircraft under condition of uncertainty |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106997208A CN106997208A (en) | 2017-08-01 |
| CN106997208B true CN106997208B (en) | 2019-09-27 |
Family
ID=59435717
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710325189.5A Active CN106997208B (en) | 2017-05-10 | 2017-05-10 | A kind of control method towards the hypersonic aircraft under condition of uncertainty |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106997208B (en) |
Families Citing this family (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107608212B (en) * | 2017-08-30 | 2020-01-24 | 湖北工业大学 | Integral sliding mode control maglev train system control method based on improved double-power approximation law |
| CN107831653B (en) * | 2017-10-16 | 2020-07-28 | 南京航空航天大学 | Hypersonic aircraft instruction tracking control method for inhibiting parameter perturbation |
| CN107992082B (en) * | 2017-12-26 | 2020-05-08 | 电子科技大学 | Four-rotor unmanned aerial vehicle flight control method based on fractional order power switching law |
| CN107943097B (en) * | 2017-12-29 | 2020-07-17 | 北京航空航天大学 | Aircraft control method and device and aircraft |
| CN108181920B (en) * | 2018-01-31 | 2021-08-31 | 天津大学 | High-precision attitude tracking control method for quad-rotor unmanned aerial vehicle based on given time |
| CN108549398B (en) * | 2018-04-24 | 2020-05-08 | 电子科技大学 | Four-rotor flight control method based on fractional order saturation function power switching law |
| CN108490786A (en) * | 2018-04-27 | 2018-09-04 | 哈尔滨工业大学 | A kind of hypersonic aircraft Robust Tracking Control based on terminal sliding mode |
| CN108829120B (en) * | 2018-05-28 | 2021-08-03 | 浙江工业大学 | Self-adaptive control method of four-rotor aircraft based on arc tangent enhanced power approach law and fast terminal sliding mode surface |
| CN108829117B (en) * | 2018-05-28 | 2021-08-03 | 浙江工业大学 | Self-adaptive control method of four-rotor aircraft based on logarithm enhancement type power approach law and fast terminal sliding mode surface |
| CN109947123B (en) * | 2019-02-27 | 2021-06-22 | 南京航空航天大学 | Unmanned aerial vehicle path tracking and autonomous obstacle avoidance method based on sight guidance law |
| CN110244768B (en) * | 2019-07-19 | 2021-11-30 | 哈尔滨工业大学 | Hypersonic aircraft modeling and anti-saturation control method based on switching system |
| CN110673544B (en) * | 2019-09-27 | 2022-11-18 | 上海大学 | Upper limb rehabilitation robot control method based on adaptive online learning |
| CN112000127B (en) * | 2020-08-13 | 2022-01-14 | 上海航天控制技术研究所 | Reverse-step-method-based aircraft lateral combined control method |
| CN112394645B (en) * | 2021-01-20 | 2021-04-16 | 中国人民解放军国防科技大学 | Neural network backstepping sliding mode control method and system for spacecraft attitude tracking |
| CN113031449B (en) * | 2021-03-29 | 2022-08-09 | 西北工业大学 | Pneumatic unknown-oriented intelligent backstepping switching control method for wide-area flight |
| CN113110028B (en) * | 2021-04-14 | 2022-06-21 | 福州大学 | Environment signal processing method of flight environment simulation control system |
| CN115657729B (en) * | 2022-12-27 | 2023-03-10 | 北京航空航天大学 | Unmanned aerial vehicle boundary protection control method considering detection task constraint |
| CN117289598B (en) * | 2023-08-01 | 2024-06-11 | 北京理工大学重庆创新中心 | Method and system for controlling backstepping sliding mode of aircraft |
| CN117394421B (en) * | 2023-09-28 | 2024-05-07 | 陕西理工大学 | Improved active disturbance rejection control method of energy storage converter based on supercoiled sliding mode observer |
| CN117519257B (en) * | 2024-01-04 | 2024-03-29 | 中国人民解放军国防科技大学 | Supersonic speed cruising altitude control method based on back-stepping method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5371669A (en) * | 1992-06-18 | 1994-12-06 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Sliding mode control method having terminal convergence in finite time |
| CN101937233A (en) * | 2010-08-10 | 2011-01-05 | 南京航空航天大学 | Nonlinear self-adaption control method of near-space hypersonic vehicle |
| CN102880053A (en) * | 2012-09-29 | 2013-01-16 | 西北工业大学 | Prediction model based hypersonic aircraft sliding-mode control method |
| CN103425135A (en) * | 2013-07-30 | 2013-12-04 | 南京航空航天大学 | Near space vehicle robust control method with input saturation |
| CN105159305A (en) * | 2015-08-03 | 2015-12-16 | 南京航空航天大学 | Four-rotor flight control method based on sliding mode variable structure |
-
2017
- 2017-05-10 CN CN201710325189.5A patent/CN106997208B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5371669A (en) * | 1992-06-18 | 1994-12-06 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Sliding mode control method having terminal convergence in finite time |
| CN101937233A (en) * | 2010-08-10 | 2011-01-05 | 南京航空航天大学 | Nonlinear self-adaption control method of near-space hypersonic vehicle |
| CN102880053A (en) * | 2012-09-29 | 2013-01-16 | 西北工业大学 | Prediction model based hypersonic aircraft sliding-mode control method |
| CN103425135A (en) * | 2013-07-30 | 2013-12-04 | 南京航空航天大学 | Near space vehicle robust control method with input saturation |
| CN105159305A (en) * | 2015-08-03 | 2015-12-16 | 南京航空航天大学 | Four-rotor flight control method based on sliding mode variable structure |
Non-Patent Citations (2)
| Title |
|---|
| 吸气式高超声速飞行器鲁棒非奇异Terminal滑模反步控制;王肖 等;《航空学报》;20170325;第38卷(第3期);第320287-1-13页 * |
| 四旋翼飞行器准滑模反步法控制设计;王婵 等;《测控技术》;20161231;第35卷(第1期);第77-80页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106997208A (en) | 2017-08-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106997208B (en) | A kind of control method towards the hypersonic aircraft under condition of uncertainty | |
| CN102759928B (en) | Control method for flight path of airship on stratosphere | |
| CN103838145B (en) | VTOL aircraft Robust Fault-Tolerant Control Systems based on cascade observer and method | |
| CN103090728B (en) | Tail angle restraining guidance method based on sliding mode control | |
| CN109426146A (en) | The high-order Non-Singular Terminal Sliding Mode control method of hypersonic aircraft | |
| Wang et al. | Nonlinear hierarchy-structured predictive control design for a generic hypersonic vehicle | |
| CN107942651A (en) | A kind of Near Space Flying Vehicles control system | |
| CN109062055A (en) | A kind of Near Space Flying Vehicles control system based on Back-stepping robust adaptive dynamic surface | |
| CN110109361B (en) | Rapid smooth self-adaptive sliding mode control method for hypersonic aircraft | |
| CN111290278B (en) | Hypersonic aircraft robust attitude control method based on prediction sliding mode | |
| CN105138012A (en) | GPS-guided unmanned aerial vehicle automatic carrier-landing adaptive control system and method | |
| CN109858106A (en) | Aircraft winglet stroke optimization method based on Gauss puppet spectrometry | |
| CN109358634A (en) | A kind of hypersonic aircraft Robust Adaptive Control method | |
| CN105759609B (en) | A kind of carrier-borne aircraft autonomous landing on the ship method based on explicit non-linear Model Predictive Control | |
| CN109164708B (en) | Neural network self-adaptive fault-tolerant control method for hypersonic aircraft | |
| CN107390531A (en) | The hypersonic aircraft control method of parameter learning finite time convergence control | |
| CN110320927A (en) | Flight control method and system of intelligent deformable aircraft | |
| CN103592847A (en) | Hypersonic aerocraft nonlinear control method based on high-gain observer | |
| CN111007720A (en) | Nonlinear self-healing control method for hypersonic aircraft multi-sensor fault | |
| CN103809446B (en) | Aircraft multiloop model bunch Flutter Suppression combination frequency robust Controller Design method | |
| CN116360255A (en) | Self-adaptive adjusting control method for nonlinear parameterized hypersonic aircraft | |
| CN114815625A (en) | Intelligent learning control method for composite control aircraft | |
| CN103777523B (en) | Aircraft multiloop model bunch Composite PID robust Controller Design method | |
| CN103809442B (en) | Aircraft multiloop model bunch combination frequency robust Controller Design method | |
| CN107489531B (en) | Hypersonic jets fuel supply rate curve design method based on semi-integral and gain-adaptive |
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 |