CN105242683A - Airship neural network terminal sliding mode track control method - Google Patents
Airship neural network terminal sliding mode track control method Download PDFInfo
- Publication number
- CN105242683A CN105242683A CN201510744363.0A CN201510744363A CN105242683A CN 105242683 A CN105242683 A CN 105242683A CN 201510744363 A CN201510744363 A CN 201510744363A CN 105242683 A CN105242683 A CN 105242683A
- Authority
- CN
- China
- Prior art keywords
- eta
- centerdot
- sin
- cos
- phi
- 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.)
- Granted
Links
Landscapes
- Feedback Control In General (AREA)
Abstract
The invention provides an airship neural network terminal sliding mode track control method, and aims at the problem of airship track control under the condition of model uncertainty. Firstly error amount is calculated by the given instruction tracks and actual tracks, then the track control law is designed by adopting a terminal sliding mode control method through selecting a terminal sliding mode function, and a neural network is adopted for approximation to the uncertain model of an airship. A system controlled by the method can precisely track the instruction tracks under the condition of model uncertainty so that adaptability is wider and robustness is high, and an effective scheme is provided for engineering realization of airship track control.
Description
Technical field
The present invention relates to a kind of flight control method of aviation field, it provides a kind of neural network TSM control method for dirigible, belongs to automatic control technology field.
Background technology
Dirigible fills buoyance lift gas (as helium, hydrogen etc.) in a kind of dependence to provide uplift, propulsion system and flight control system is relied on to realize handling the floating class aircraft of flight, there is the advantages such as airborne period is long, energy consumption is low, reusable, be particularly suitable as carrying platform, novel electronic information equipment is developed into by carrying multiple useful load, the fields such as land mapping, disaster monitoring, regional early warning, reconnaissance and surveillance can be widely used in, there is significant application value and wide application prospect, the current study hotspot having become aviation field.
Flight tracking control refers to the control to dirigible center of mass motion, usually requires that dirigible flies according to prebriefed pattern, to complete specific aerial mission.The spatial movement of dirigible has non-linear, passage coupling, the feature such as uncertain, and therefore, flight tracking control becomes the difficulties that airship flight controls.Existing document mostly based on inearized model, does not consider non-linear factor and the coupling longitudinally and between horizontal sideway movement to the research of dirigible flight tracking control, only effective near equilibrium state.Sliding-mode control has robustness to Parameter Perturbation and external interference, for dirigible flight tracking control provides a kind of effective means.But sliding formwork controls usually to adopt linear sliding mode, and after system arrives sliding-mode surface, status tracking error asymptotic convergence to zero, cannot at Finite-time convergence.For this problem, patent of invention " a kind of dirigible non-singular terminal sliding formwork flight tracking control method " (patent No.: ZL201410623630.4), propose a kind of non-singular terminal sliding formwork flight tracking control method, making flight tracking control error at Finite-time convergence to zero by choosing terminal sliding mode function, improve response speed and the control accuracy of flight tracking control.But the method does not consider dirigible model uncertain problem, therefore its adaptability and robustness need to be improved further.
Summary of the invention
For solving the problem, the invention provides a kind of dirigible neural network terminal sliding mode flight tracking control method.Neural network has the ability of approaching any complex nonlinear function, can be applied to the modeling problem of uncertain system, can significantly improve the control performance of system.The present invention, on the basis of non-singular terminal sliding formwork flight tracking control method, for dirigible model uncertain problem, adopts neural network to approach the ambiguous model of dirigible.
Flight path control system structured flowchart proposed by the invention as shown in Figure 1.As shown in Figure 1, first according to instruction flight path η
dwith actual flight path η error of calculation amount e; Then choose terminal sliding mode face according to margin of error e, adopt TSM control method design flight tracking control rule; For solving model uncertain problem, according to
constructing neural network approaches device, the ambiguous model of On-line Estimation dirigible, thus obtains neural network terminal sliding mode flight tracking control rule.The system controlled by the method can under model condition of uncertainty high precision tracking instruction flight path, compared to patent of invention " a kind of dirigible non-singular terminal sliding formwork flight tracking control method ", there is wider adaptability and strong robustness, for the Project Realization of dirigible flight tracking control provides effective scheme.
A kind of dirigible neural network of the present invention terminal sliding mode flight tracking control method, first by given instruction flight path and actual flight path error of calculation amount, then by choosing terminal sliding mode function, adopt TSM control method design flight tracking control rule, and adopt neural network to approach the ambiguous model of dirigible.In practical application, dirigible flight path is obtained by integrated navigation system measurement, the controlled quentity controlled variable calculated is transferred to topworks can realize flight tracking control function by the method.
A kind of dirigible neural network terminal sliding mode flight tracking control method, its concrete steps are as follows, as shown in Figure 2:
Step one: the given instruction flight path represented with generalized coordinate: η
d=[x
d, y
d, z
d, θ
d, ψ
d, φ
d]
t;
Described instruction flight path is generalized coordinate η
d=[x
d, y
d, z
d, θ
d, ψ
d, φ
d]
t, wherein: x
d, y
d, z
d, θ
d, ψ
dand φ
dbe respectively instruction x coordinate, instruction y coordinate, instruction z coordinate, the instruction angle of pitch, instruction crab angle and instruction roll angle, subscript T represents vector or transpose of a matrix.
Step 2: the margin of error calculates: the margin of error between computations flight path and actual flight path;
The margin of error between instruction flight path and actual flight path, its computing method are:
e=η-η
d=[x-x
d,y-y
d,z-z
d,θ-θ
d,ψ-ψ
d,φ-φ
d]
T(1)
η=[x, y, z, θ, ψ, φ]
tfor actual flight path, x, y, z, θ, ψ, φ are respectively the x coordinate of actual flight path, y coordinate, z coordinate, the angle of pitch, crab angle and roll angle.
Step 3: TSM control rule design: choose terminal sliding mode function, adopts TSM control method design flight tracking control rule, and adopts neural network to approach the ambiguous model of dirigible, calculate flight tracking control amount;
1) mathematical model of dirigible spatial movement is set up
For ease of describing, coordinate system and the kinematic parameter of dirigible spatial movement are defined as follows.As shown in Figure 3, earth axes O is adopted
ex
ey
ez
ewith body coordinate system o
bx
by
bz
bbe described the spatial movement of dirigible, CV is centre of buoyancy, and CG is center of gravity, and centre of buoyancy is r to the vector of center of gravity
g=[x
g, y
g, z
g]
t, x
gfor centre of buoyancy is to the axial component of center of gravity, y
gfor centre of buoyancy is to the cross component of center of gravity, z
gfor centre of buoyancy is to the vertical component of center of gravity; Kinematic parameter defines: position P=[x, y, z]
t, x, y, z is respectively the displacement of axis, side direction and vertical direction; Attitude angle Ω=[θ, ψ, φ]
t, θ, ψ, φ are respectively the angle of pitch, crab angle and roll angle; Speed v=[u, v, w]
t, u, v, w are respectively the speed of axis, side direction and vertical direction in body coordinate system; Angular velocity omega=[p, q, r]
t, p, q, r are respectively rolling, pitching and yaw rate; Note generalized coordinate η=[x, y, z, θ, ψ, φ]
t, generalized velocity is V=[u, v, w, p, q, r]
t.
The mathematical model of dirigible spatial movement is described below:
In formula, 0
3 × 3represent the null matrix on 3 × 3 rank,
represent the first differential of η,
represent the first differential of V,
Wherein, G represents the gravity suffered by dirigible, and B represents the buoyancy suffered by dirigible, and F represents the thrust suffered by dirigible,
N
u=(m+m
22)vr-(m+m
33)wq+m[x
G(p
2+r
2)-y
Gpq-z
Gpr]
(10)
+QV
2/3(-C
Xcosαcosβ+C
Ycosαsinβ+C
Zsinα)
N
v=(m+m
33)wp-(m+m
11)ur-m[x
Gpq-y
G(p
2+r
2)+z
Gqr]
(11)
+QV
2/3(C
Xsinβ+C
Ycosβ)
N
w=(m+m
22)vp-(m+m
11)uq-m[x
Gpr+y
Gqr-z
G(p
2+q
2)]
(12)
+QV
2/3(-C
Xsinαsinβ+C
Ysinαcosβ-C
Zcosα)
N
p=[(I
y+I
22)-(I
z+I
33)]qr+I
xzpq-I
xypr-I
yz(r
2-q
2)+
(13)
[mz
G(ur-wp)+y
G(uq-vp)]+QVC
l
N
q=[(I
z+I
33)-(I
x+I
11)]pr+I
xyqr-I
yzpq-I
xz(p
2-r
2)
(14)
+m[x
G(vp-uq)-z
G(wp-vr)]+QVC
m
N
r=[(I
y+I
22)-(I
x+I
11)]pq-I
xzqr-I
xy(q
2-p
2)+I
yzpr
(15)
+m[y
G(wq-vr)-x
G(ur-wp)]+QVC
n
In formula, m is dirigible quality, m
11, m
22, m
33for additional mass, I
11, I
22, I
33for additional inertial; Q is dynamic pressure, and α is the angle of attack, and β is yaw angle, C
x, C
y, C
z, C
l, C
m, C
nfor Aerodynamic Coefficient; I
x, I
y, I
zbe respectively around o
bx
b, o
by
b, o
bz
bprincipal moments; I
xy, I
xz, I
yzbe respectively about plane o
bx
by
b, o
bx
bz
b, o
by
bz
bproduct of inertia; T is thrust size, and μ is thrust vectoring and o
bx
bz
bangle between face, specifies that it is at o
bx
bz
bthe left side in face is just, υ is that thrust vectoring is at o
bx
bz
bthe projection in face and o
bx
bangle between axle, specifies that it is projected in o
bx
bjust be under axle; l
x, l
y, l
zrepresent that thrust point is apart from initial point o
bdistance.
Formula (3) is the expression formula about generalized velocity V, needs to be transformed to the expression formula about generalized coordinate η.
Can be obtained by formula (1):
In formula, J
-1(η) be the inverse matrix of J (η),
To formula (16) differential, can obtain
In formula
Formula (19) premultiplication
can obtain
Composite type (3), formula (19) and formula (21) can obtain:
In formula
M
η=R
TMR(23)
In practical flight process, M
η, N
ηand G
ηall there is indeterminate, be respectively Δ M
η, Δ N
ηwith Δ G
η; Formula (22) can be written as thus:
Note
Then formula (27) can be written as:
With the mathematical model described by formula (28) for controlled device, adopt neural network TSM control method design flight tracking control rule.
2) sliding-mode surface design
Design terminal sliding-mode surface is:
Wherein, e
1, e
2, e
3, e
4, e
5, e
6be respectively the 1st to the 6th element of vectorial e, s=[s
1, s
2, s
3, s
4, s
5, s
6]
t, s
1, s
2, s
3, s
4, s
5, s
6be respectively the 1st to the 6th element of vectorial s, λ=diag (λ
1, λ
2, λ
3, λ
4, λ
5, λ
6), λ
1, λ
2, λ
3, λ
4, λ
5, λ
6be respectively the 1st to the 6th element of vectorial λ, diag () represents diagonal matrix, and λ is positive definite matrix, and κ is arithmetic number and meets 1 < κ < 2.
3) design terminal sliding formwork control law:
In formula, λ
-1represent the inverse matrix of λ,
represent M
ηinverse matrix, || || represent euclideam norm.
The stability analysis of sliding-mode surface is as follows.
Be defined as follows Lyapunov function
Formula (29) utilizes to formula (31) differential, can obtain:
In formula,
represent
first differential, s
trepresent the transposition of s.
Second-order differential asked to formula (1) and utilize formula (22) and formula (29), can obtain:
Formula (33) is substituted into formula (32), can obtain:
Formula (34) namely demonstrate,proves the stability of sliding-mode surface.
Model indeterminate Δ f in formula (30) is in fact not known, and therefore, through type (30) can not accurately provide flight tracking control amount.For this problem, the present invention adopts neural network approximate model indeterminate Δ f.
4) constructing neural network approaches device, design neural network TSM control rule:
Neural network approximator comprises input layer, hidden layer and output layer, as shown in Figure 4.
Input layer: the input variable choosing network is
E
tfor margin of error e transposition,
for the first differential of margin of error e transposition,
for transposition, the η of the second-order differential of margin of error e
tfor η transposition,
for the first differential of η transposition,
for the transposition of the second-order differential of η.
Hidden layer: choose the basis function of Gaussian function as hidden node
Wherein, c
j=[c
j1, c
j2..., c
jn] be the intermediate value of a jth Gaussian function, c
j1, c
j2..., c
jnbe respectively c
jthe 1st, the 2nd ..., the n-th component, n is the nodes of neural network; b
jfor the standard deviation of a jth Gaussian function, j=1,2 ..., n, n are the nodes of neural network, || || represent euclideam norm.
Output layer: the output of neural network approximator is
In formula,
for || the estimated value of Δ f||,
for optimal weights coefficient vector, h (x)=[h
1(x), h
2(x) ..., h
n(x)]
t, h
1(x), h
2(x) ..., h
n(x) be respectively the 1st of vector function, the 2nd ..., the n-th component, n is the nodes of neural network.
Thus, neural network TSM control rule is:
(3) advantage and effect:
Compared with prior art, its advantage is in the present invention:
1) the method directly designs based on the non-linear dynamic model of dirigible spatial movement, considers every non-linear factor and the coupling longitudinally and between horizontal sideway movement, overcomes the limitation that inearized model is only suitable for equilibrium state.
2) TSM control makes attitude control error at Finite-time convergence to zero by choosing terminal sliding mode function, overcome the asymptotic convergence problem that traditional sliding formwork controls, there is rapid dynamic response speed, finite time convergence control, steady-state tracking precision advantages of higher.
3) the method adopts neural network to approach the ambiguous model of dirigible, overcomes the flight model uncertain problem in practical flight process, has wider adaptability and strong robustness.
Control engineering teacher can according to the given arbitrary instruction flight path of actual dirigible in application process, and the controlled quentity controlled variable obtained by the method is transferred to topworks realizes flight tracking control function.
Accompanying drawing explanation
Fig. 1 is dirigible flight path control system structural drawing of the present invention
Fig. 2 is dirigible flight tracking control method step process flow diagram of the present invention
Fig. 3 is dirigible coordinate system of the present invention and kinematic parameter definition
Fig. 4 is neural network approximator structural drawing of the present invention
Fig. 5 is dirigible flight tracking control result of the present invention
Fig. 6 is dirigible flight tracking control error of the present invention
Fig. 7 is neural network Approaching Results of the present invention
In figure, symbol description is as follows:
η η=[x, y, z, θ, ψ, φ]
tfor dirigible flight path, wherein x, y, z, θ, ψ, φ are respectively the x coordinate of actual flight path, y coordinate, z coordinate, the angle of pitch, crab angle and roll angle;
for the first differential of η;
for the second-order differential of η;
η
dη
d=[x
d, y
d, z
d, θ
d, ψ
d, φ
d]
tfor instruction flight path, wherein x
d, y
d, z
d, θ
d, ψ
dand φ
dbe respectively instruction x coordinate, instruction y coordinate, instruction z coordinate, the instruction angle of pitch, instruction crab angle and instruction roll angle;
O
ex
ey
ez
eo
ex
ey
ez
erepresent earth axes;
O
bx
by
bz
bo
bx
by
bz
brepresent dirigible body coordinate system;
Ee=[x
e, y
e, z
e, θ
e, ψ
e, φ
e]
tfor flight tracking control error, x
e, y
e, z
e, θ
e, ψ
eand φ
ebe respectively the x error of coordinate of flight tracking control, y error of coordinate, z coordinate error, angle of pitch error, crab angle error and roll angle error;
for the first differential of margin of error e;
for the second-order differential of margin of error e; Ss is sliding-mode surface;
U
nu
nfor neural network terminal sliding mode flight tracking control amount;
represent first differential computing.
Embodiment
In order to make object of the present invention, technical scheme and beneficial effect clearly understand, below in conjunction with drawings and Examples, the present invention is further elaborated.It should be noted that specific embodiment described herein only in order to explain the present invention, be not intended to limit the present invention.
Step one: given instruction flight path
Given instruction flight path is:
η
d=[x
d, y
d, z
d, θ
d, ψ
d, φ
d]
t=[(1.5t) m, 200sin (0.005t) m, 10m, 0rad, 0.02rad, 0rad]
t, x
d, y
d, z
d, θ
d, ψ
dand φ
dbe respectively instruction x coordinate, instruction y coordinate, instruction z coordinate, the instruction angle of pitch, instruction crab angle and instruction roll angle;
Step 2: the margin of error calculates
The margin of error between computations flight path and actual flight path:
e=η-η
d=[x-x
d,y-y
d,z-z
d,θ-θ
d,ψ-ψ
d,φ-φ
d]
T,
Wherein, η=[x, y, z, θ, ψ, φ]
tfor actual flight path, x, y, z, θ, ψ, φ are respectively the x coordinate of actual flight path, y coordinate, z coordinate, the angle of pitch, crab angle and roll angle, are consecutive variations value.
Initial flight path is:
η
0=[x
0,y
0,z
0,θ
0,ψ
0,φ
0]
T=[50m,-100m,8m,0.01rad,0.01rad,0.01rad]
T。
Initial velocity:
V
0=[u
0,v
0,w
0,p
0,q
0,r
0]
T=[5m/s,2.5m/s,0m/s,0rad/s,0rad/s,0rad/s]
T
Step 3: design flight tracking control rule:
1) mathematical model of dirigible spatial movement is set up
The mathematical model of dirigible spatial movement can be expressed as:
In formula
Wherein
N
u=(m+m
22)vr-(m+m
33)wq+m[x
G(p
2+r
2)-y
Gpq-z
Gpr]
(46)
+QV
2/3(-C
Xcosαcosβ+C
Ycosαsinβ+C
Zsinα)
N
v=(m+m
33)wp-(m+m
11)ur-m[x
Gpq-y
G(p
2+r
2)+z
Gqr]
(47)
+QV
2/3(C
Xsinβ+C
Ycosβ)
N
w=(m+m
22)vp-(m+m
11)uq-m[x
Gpr+y
Gqr-z
G(p
2+q
2)]
(48)
+QV
2/3(-C
Xsinαsinβ+C
Ysinαcosβ-C
Zcosα)
N
p=[(I
y+I
22)-(I
z+I
33)]qr+I
xzpq-I
xypr-I
yz(r
2-q
2)+
(49)
[mz
G(ur-wp)+y
G(uq-vp)]+QVC
l
N
q=[(I
z+I
33)-(I
x+I
11)]pr+I
xyqr-I
yzpq-I
xz(p
2-r
2)
(50)
+m[x
G(vp-uq)-z
G(wp-vr)]+QVC
m
N
r=[(I
y+I
22)-(I
x+I
11)]pq-I
xzqr-I
xy(q
2-p
2)+I
yzpr
(51)
+m[y
G(wq-vr)-x
G(ur-wp)]+QVC
n
In formula, m is dirigible quality, m
11, m
22, m
33for additional mass, I
11, I
22, I
33for additional inertial; Q is dynamic pressure, and α is the angle of attack, and β is yaw angle, C
x, C
y, C
z, C
l, C
m, C
nfor Aerodynamic Coefficient; I
x, I
y, I
zbe respectively around o
bx
b, o
by
b, o
bz
bprincipal moments; I
xy, I
xz, I
yzbe respectively about plane o
bx
by
b, o
bx
bz
b, o
by
bz
bproduct of inertia; T is thrust size, and μ is thrust vectoring and o
bx
bz
bangle between face, specifies that it is at o
bx
bz
bthe left side in face is just, υ is that thrust vectoring is at o
bx
bz
bthe projection in face and o
bx
bangle between axle, specifies that it is projected in o
bx
bjust be under axle; l
x, l
y, l
zrepresent that thrust point is apart from initial point o
bdistance.
Formula (39) is transformed to the expression formula about generalized coordinate η.
Can be obtained by formula (38):
In formula, J
-1(η) be the inverse matrix of J (η),
To formula (52) differential, can obtain
In formula
Formula (55) premultiplication
can obtain
Composite type (39), formula (55) and formula (57) can obtain:
In formula
M
η=R
TMR(59)
Consider model indeterminate, formula (58) can be written as:
In formula,
wherein, Δ M
η, Δ N
ηwith Δ G
ηbe respectively M
η, N
ηand G
ηunknown term and indeterminate.
Dirigible parameter in the present embodiment is in table 2.
Table 2 dirigible parameter
2) sliding-mode surface design
Design terminal sliding-mode surface is:
Wherein, diag (2,2,2,2,2,2), diag () represents diagonal matrix, κ=5/3.
3) design terminal sliding formwork control law:
4) constructing neural network approaches device, design neural network TSM control rule:
Neural network approximator comprises input layer, hidden layer and output layer.
Input layer: the input variable choosing network is
Hidden layer: choose the basis function of Gaussian function as hidden node, nodes n=7,
Wherein,
b
j=0.12。
Output layer: the output of neural network approximator is
Wherein,
be taken as that element is 0.01, dimension is the vector of 42.
Neural network TSM control rule is:
Dirigible Three-dimensional Track tracking results in embodiment as shown in Figure 5-Figure 7.Fig. 5 gives dirigible flight tracking control result, can be obtained by Fig. 5: dirigible is by initial position (50m,-100m, 8m) set out and to reach home position (450m through 300s, 200m, 10m), can accurate trace command flight path, demonstrate the validity of flight tracking control method proposed by the invention; Fig. 6 gives flight tracking control error, can be obtained by Fig. 6: dirigible with zero steady-state error trace command flight path, can have higher control accuracy.Fig. 7 gives neural network Approaching Results, can be obtained by Fig. 7, and neural network approximator proposed by the invention can approach the ambiguous model of dirigible preferably.
Claims (1)
1. a dirigible neural network terminal sliding mode flight tracking control method, is characterized in that, comprise the following steps:
Step one: the given instruction flight path represented with generalized coordinate: η
d=[x
d, y
d, z
d, θ
d, ψ
d, φ
d]
t;
Described instruction flight path is generalized coordinate η
d=[x
d, y
d, z
d, θ
d, ψ
d, φ
d]
t, wherein: x
d, y
d, z
d, θ
d, ψ
dand φ
dbe respectively instruction x coordinate, instruction y coordinate, instruction z coordinate, the instruction angle of pitch, instruction crab angle and instruction roll angle, subscript T represents vector or transpose of a matrix;
Step 2: the margin of error calculates: the margin of error between computations flight path and actual flight path;
The margin of error between instruction flight path and actual flight path, its computing method are:
e=η-η
d=[x-x
d,y-y
d,z-z
d,θ-θ
d,ψ-ψ
d,φ-φ
d]
T(1)
η=[x, y, z, θ, ψ, φ]
tfor actual flight path, x, y, z, θ, ψ, φ are respectively the x coordinate of actual flight path, y coordinate, z coordinate, the angle of pitch, crab angle and roll angle;
Step 3: TSM control rule design: choose terminal sliding mode function, adopts TSM control method design flight tracking control rule, and adopts neural network to approach the ambiguous model of dirigible, calculate flight tracking control amount;
1) mathematical model of dirigible spatial movement is set up
Coordinate system and the kinematic parameter of dirigible spatial movement are defined as follows: adopt earth axes O
ex
ey
ez
ewith body coordinate system o
bx
by
bz
bbe described the spatial movement of dirigible, CV is centre of buoyancy, and CG is center of gravity, and centre of buoyancy is r to the vector of center of gravity
g=[x
g, y
g, z
g]
t, x
gfor centre of buoyancy is to the axial component of center of gravity, y
gfor centre of buoyancy is to the cross component of center of gravity, z
gfor centre of buoyancy is to the vertical component of center of gravity; Kinematic parameter defines: position P=[x, y, z]
t, x, y, z is respectively the displacement of axis, side direction and vertical direction; Attitude angle Ω=[θ, ψ, φ]
t, θ, ψ, φ are respectively the angle of pitch, crab angle and roll angle; Speed v=[u, v, w]
t, u, v, w are respectively the speed of axis, side direction and vertical direction in body coordinate system; Angular velocity omega=[p, q, r]
t, p, q, r are respectively rolling, pitching and yaw rate; Note generalized coordinate η=[x, y, z, θ, ψ, φ]
t, generalized velocity is V=[u, v, w, p, q, r]
t;
The mathematical model of dirigible spatial movement is described below:
In formula, 0
3 × 3represent the null matrix on 3 × 3 rank,
represent the first differential of η,
represent the first differential of V,
Wherein, G represents the gravity suffered by dirigible, and B represents the buoyancy suffered by dirigible, and F represents the thrust suffered by dirigible,
In formula, m is dirigible quality, m
11, m
22, m
33for additional mass, I
11, I
22, I
33for additional inertial; Q is dynamic pressure, and α is the angle of attack, and β is yaw angle, C
x, C
y, C
z, C
l, C
m, C
nfor Aerodynamic Coefficient; I
x, I
y, I
zbe respectively around o
bx
b, o
by
b, o
bz
bprincipal moments; I
xy, I
xz, I
yzbe respectively about plane o
bx
by
b, o
bx
bz
b, o
by
bz
bproduct of inertia; T is thrust size, and μ is thrust vectoring and o
bx
bz
bangle between face, specifies that it is at o
bx
bz
bthe left side in face is just, υ is that thrust vectoring is at o
bx
bz
bthe projection in face and o
bx
bangle between axle, specifies that it is projected in o
bx
bjust be under axle; l
x, l
y, l
zrepresent that thrust point is apart from initial point o
bdistance;
Formula (3) is the expression formula about generalized velocity V, needs to be transformed to the expression formula about generalized coordinate η;
Can be obtained by formula (1):
J in formula
-1(η) be the inverse matrix of J (η),
To formula (16) differential, can obtain
In formula
Formula (19) premultiplication
can obtain
Composite type (3), formula (19) and formula (21) can obtain:
In formula
M
η=R
TMR(23)
In practical flight process, M
η, N
ηand G
ηall there is indeterminate, be respectively Δ M
η, Δ N
ηwith Δ G
η; Formula (22) can be written as thus:
Note
Then formula (27) can be written as:
With the mathematical model described by formula (28) for controlled device, adopt neural network TSM control method design flight tracking control rule;
2) sliding-mode surface design
Design terminal sliding-mode surface is:
Wherein, e=[e
1, e
2, e
3, e
4, e
5, e
6]
t, e
1, e
2, e
3, e
4, e
5, e
6be respectively the 1st to the 6th element of vectorial e, s=[s
1, s
2, s
3, s
4, s
5, s
6]
t, s
1, s
2, s
3, s
4, s
5, s
6be respectively the 1st to the 6th element of vectorial s, λ=diag (λ
1, λ
2, λ
3, λ
4, λ
5, λ
6), λ
1, λ
2, λ
3, λ
4, λ
5, λ
6be respectively the 1st to the 6th element of vectorial λ, diag () represents diagonal matrix, and λ is positive definite matrix, and κ is arithmetic number and meets 1 < κ < 2;
3) design terminal sliding formwork control law:
In formula, λ
-1represent the inverse matrix of λ,
represent M
ηinverse matrix, || || represent euclideam norm;
The stability analysis of sliding-mode surface is as follows;
Be defined as follows Lyapunov function
Formula (29) utilizes to formula (31) differential, can obtain:
In formula,
represent
first differential, s
trepresent the transposition of s;
Second-order differential asked to formula (1) and utilize formula (22) and formula (29), can obtain:
Formula (33) is substituted into formula (32), can obtain:
Formula (34) namely demonstrate,proves the stability of sliding-mode surface;
For the model indeterminate Δ f in formula (30), adopt neural network approximate model indeterminate Δ f;
4) constructing neural network approaches device, design neural network TSM control rule:
Neural network approximator comprises input layer, hidden layer and output layer;
Input layer: the input variable choosing network is
e
tfor margin of error e transposition,
for the first differential of margin of error e transposition,
for transposition, the η of the second-order differential of margin of error e
tfor η transposition,
for the first differential of η transposition,
for the transposition of the second-order differential of η;
Hidden layer: choose the basis function of Gaussian function as hidden node
Wherein, c
j=[c
j1, c
j2..., c
jn] be the intermediate value of a jth Gaussian function, c
j1, c
j2..., c
jnbe respectively c
jthe 1st, the 2nd ..., the n-th component, n is the nodes of neural network; b
jfor the standard deviation of a jth Gaussian function, j=1,2 ..., n, n are the nodes of neural network, || || represent euclideam norm;
Output layer: the output of neural network approximator is
In formula,
for || Δ f| estimated value,
for optimal weights coefficient vector, h (x)=[h
1(x), h
2(x) ..., h
n(x)]
t, h
1(x), h
2(x) ..., h
n(x) be respectively the 1st of vector function, the 2nd ..., the n-th component, n is the nodes of neural network;
Thus, neural network TSM control rule is:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510744363.0A CN105242683B (en) | 2015-11-04 | 2015-11-04 | A kind of dirigible neutral net terminal sliding mode flight tracking control method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510744363.0A CN105242683B (en) | 2015-11-04 | 2015-11-04 | A kind of dirigible neutral net terminal sliding mode flight tracking control method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105242683A true CN105242683A (en) | 2016-01-13 |
CN105242683B CN105242683B (en) | 2016-07-06 |
Family
ID=55040359
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510744363.0A Active CN105242683B (en) | 2015-11-04 | 2015-11-04 | A kind of dirigible neutral net terminal sliding mode flight tracking control method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105242683B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106372307A (en) * | 2016-08-30 | 2017-02-01 | 中国航空工业集团公司西安飞行自动控制研究所 | Civil aircraft air flow angle estimation method based on pneumatic model |
CN107065847A (en) * | 2016-03-31 | 2017-08-18 | 北京航空航天大学 | A kind of surface vessel Trajectory Tracking Control method of the asymmetric saturation of actuator |
CN107450319A (en) * | 2017-08-23 | 2017-12-08 | 中国人民解放军国防科技大学 | Designated time nonsingular terminal sliding mode control method for airship track tracking |
CN107703967A (en) * | 2017-10-24 | 2018-02-16 | 中国人民解放军国防科技大学 | Control method for controlling track of limited airship |
CN107817818A (en) * | 2017-12-07 | 2018-03-20 | 中国人民解放军国防科技大学 | Finite time control method for flight path tracking of uncertain model airship |
CN108254741A (en) * | 2018-01-16 | 2018-07-06 | 中国人民解放军海军航空大学 | Targetpath Forecasting Methodology based on Recognition with Recurrent Neural Network |
CN118151678A (en) * | 2024-05-13 | 2024-06-07 | 山东科技大学 | Stratospheric airship fixed time track tracking control method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6317658B1 (en) * | 1997-10-15 | 2001-11-13 | The Boeing Company | Neurocomputing control distribution system |
CN102759928A (en) * | 2012-04-28 | 2012-10-31 | 中国人民解放军国防科学技术大学 | Control method for flight path of airship on stratosphere |
CN103777635A (en) * | 2014-01-13 | 2014-05-07 | 哈尔滨工程大学 | Robust self-adaptive track tracking control system for dynamic positioning vessel |
CN104281155A (en) * | 2014-11-07 | 2015-01-14 | 中国人民解放军国防科学技术大学 | Three-dimensional flight path tracking method for unmanned airship |
CN104360686A (en) * | 2014-11-07 | 2015-02-18 | 中国人民解放军国防科学技术大学 | Nonsingular terminal sliding mode flight path control method for airships |
-
2015
- 2015-11-04 CN CN201510744363.0A patent/CN105242683B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6317658B1 (en) * | 1997-10-15 | 2001-11-13 | The Boeing Company | Neurocomputing control distribution system |
CN102759928A (en) * | 2012-04-28 | 2012-10-31 | 中国人民解放军国防科学技术大学 | Control method for flight path of airship on stratosphere |
CN103777635A (en) * | 2014-01-13 | 2014-05-07 | 哈尔滨工程大学 | Robust self-adaptive track tracking control system for dynamic positioning vessel |
CN104281155A (en) * | 2014-11-07 | 2015-01-14 | 中国人民解放军国防科学技术大学 | Three-dimensional flight path tracking method for unmanned airship |
CN104360686A (en) * | 2014-11-07 | 2015-02-18 | 中国人民解放军国防科学技术大学 | Nonsingular terminal sliding mode flight path control method for airships |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107065847A (en) * | 2016-03-31 | 2017-08-18 | 北京航空航天大学 | A kind of surface vessel Trajectory Tracking Control method of the asymmetric saturation of actuator |
CN106372307A (en) * | 2016-08-30 | 2017-02-01 | 中国航空工业集团公司西安飞行自动控制研究所 | Civil aircraft air flow angle estimation method based on pneumatic model |
CN106372307B (en) * | 2016-08-30 | 2020-04-07 | 中国航空工业集团公司西安飞行自动控制研究所 | Civil aircraft airflow angle estimation method based on pneumatic model |
CN107450319A (en) * | 2017-08-23 | 2017-12-08 | 中国人民解放军国防科技大学 | Designated time nonsingular terminal sliding mode control method for airship track tracking |
CN107450319B (en) * | 2017-08-23 | 2020-10-30 | 中国人民解放军国防科技大学 | Designated time nonsingular terminal sliding mode control method for airship track tracking |
CN107703967A (en) * | 2017-10-24 | 2018-02-16 | 中国人民解放军国防科技大学 | Control method for controlling track of limited airship |
CN107817818A (en) * | 2017-12-07 | 2018-03-20 | 中国人民解放军国防科技大学 | Finite time control method for flight path tracking of uncertain model airship |
CN108254741A (en) * | 2018-01-16 | 2018-07-06 | 中国人民解放军海军航空大学 | Targetpath Forecasting Methodology based on Recognition with Recurrent Neural Network |
CN108254741B (en) * | 2018-01-16 | 2021-02-09 | 中国人民解放军海军航空大学 | Target track prediction method based on cyclic neural network |
CN118151678A (en) * | 2024-05-13 | 2024-06-07 | 山东科技大学 | Stratospheric airship fixed time track tracking control method |
Also Published As
Publication number | Publication date |
---|---|
CN105242683B (en) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104281155B (en) | Three-dimensional flight path tracking method for unmanned airship | |
CN105242683A (en) | Airship neural network terminal sliding mode track control method | |
CN104360686B (en) | Nonsingular terminal sliding mode flight path control method for airships | |
Jehle et al. | Applied tracking control for kite power systems | |
CN104793629B (en) | A kind of contragradience neural network control method of dirigible Three-dimensional Track tracking | |
CN102495633B (en) | Posture control method for airship staying on stratosphere | |
CN103425135A (en) | Near space vehicle robust control method with input saturation | |
CN102759928A (en) | Control method for flight path of airship on stratosphere | |
CN111142562B (en) | Formation transformation control method under hybrid condition constraint based on stress matrix | |
CN103853157A (en) | Aircraft attitude control method based on self-adaptive sliding mode | |
CN104950901A (en) | Nonlinear robust control method with finite-time convergence capacity for unmanned helicopter attitude error | |
CN108181806B (en) | Space robot position and attitude active disturbance rejection control method based on sampling output | |
Liu et al. | Underactuated stratospheric airship trajectory control using an adaptive integral backstepping approach | |
Yuan et al. | Finite-time trajectory tracking control for a stratospheric airship with full-state constraint and disturbances | |
CN104793498A (en) | Unmanned dirigible airship fixed-point resident control method | |
Li et al. | Anti-disturbance control for attitude and altitude systems of the helicopter under random disturbances | |
CN107817818A (en) | Finite time control method for flight path tracking of uncertain model airship | |
Baldi et al. | Towards adaptive autopilots for fixed-wing unmanned aerial vehicles | |
Liang et al. | Active disturbance rejection attitude control for a bird-like flapping wing micro air vehicle during automatic landing | |
Zhou et al. | Spatial path following control of an autonomous underactuated airship | |
CN105116905A (en) | Aircraft attitude control method | |
Xu et al. | Multivariable adaptive distributed leader-follower flight control for multiple UAVs formation | |
Rosa et al. | Autolanding controller for a fixed wing unmanned air vehicle | |
Cao et al. | Passive fault tolerant control approach for hypersonic vehicle with actuator loss of effectiveness faults | |
CN107450319A (en) | Designated time nonsingular terminal sliding mode control method for airship track tracking |
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 |