CN104267605A - Smooth nonsingular terminal sliding-mode control method suitable for control system with relative degree of 1 - Google Patents
Smooth nonsingular terminal sliding-mode control method suitable for control system with relative degree of 1 Download PDFInfo
- Publication number
- CN104267605A CN104267605A CN201410578260.7A CN201410578260A CN104267605A CN 104267605 A CN104267605 A CN 104267605A CN 201410578260 A CN201410578260 A CN 201410578260A CN 104267605 A CN104267605 A CN 104267605A
- Authority
- CN
- China
- Prior art keywords
- prime
- control system
- state
- terminal sliding
- represent
- 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 discloses a smooth nonsingular terminal sliding-mode control method suitable for a control system with a relative degree of 1, relates to a nonsingular terminal sliding-mode control method, and solves the problems that the existing nonsingular terminal sliding-mode control method exists buffeting to cause a controller not to output continuous smooth control signals and not be applied to the control system with the relative degree of 1. The smooth nonsingular terminal sliding-mode control method comprises the following steps: if the controlled system serves as a single-input single-output control system with the relative degree of 1, then acquiring the system state differential in real time, designing a nonsingular terminal sliding mode, introducing a virtual controlled quantity and utilizing an integral action to enable the actual output controlled quantity to be smooth and continuous; if the controlled system serves as a matched multi-input multi-output control system, then converting a scalar realization form of the control system to a matrix vector realization form; if the controlled system serves as a non-matching indeterminate multi-input multi-output control system, then conducting nonsingular state change for two times according to the controllability index r to resolve the system into r subsystems, introducing a reference model to eliminate time-varying uncertainty of input channels of the subsystems, and further designing a control law and a minor control law. The smooth nonsingular terminal sliding-mode control method disclosed by the invention can be applied to smooth nonsingular terminal sliding-mode control of the control system with the relative degree of 1.
Description
Technical field
The present invention relates to and be a kind ofly applicable to the level and smooth non-singular terminal sliding-mode control that Relative order is 1 control system.
Background technology
Existing non-singular terminal sliding-mode control is only applicable to the control system of Relative order >=2, and application has certain limitation.
(1) non-singular terminal sliding-mode control is a kind of nonlinear sliding mode control method occurred in recent years, can make system state overall situation finite time convergence control, therefore have applications well prospect at high speed, high precision control field.But, because Controller gain variations still continues to use traditional Sliding Mode Control Design Method, be only applied to the Mechatronic Systems of Relative order >=2 at present, as mechanical arm, permagnetic synchronous motor, DC converter etc.
(2) there is buffeting problem in existing non-singular terminal sliding-mode control, serious its practical engineering application of restriction.
With shape as
second order control be example, non-singular terminal sliding mode is typically designed to:
wherein, system state
known, design parameter β >0, p and q are odd number, and 1<p/q<2.Visible, the Relative order of non-singular terminal sliding mode s is 1, and namely controlled quentity controlled variable u is aobvious is contained in
in.At Lyapunov stable condition
under constraint, directly comprise in controlled quentity controlled variable u and switch control item sgn (.), to ensure that system state arrives and maintains on the sliding-mode surface that designs in advance, and to Parameter Perturbation and external disturbance, there is robustness.But in actual control system, the fast switching controlled frequency of theory unlimited cannot realize, and makes system pass through sliding-mode surface back and forth, produces buffeting problem, cause system that vibration and unstable, serious its practical application of restriction occur.
Be different from single input single output control system, the robust stability of strong coupling restriction multi-input multi-output system controls.Complicated relevance is there is between different input channel and input quantity, system state, output quantity, utilize existing decoupling algorithm that System Model Reduction can be made for canonical structure, namely certain controlled quentity controlled variable can according to the one or more controlled volume of a certain pass effect, from simply disturb different, the stable single-loop system of multiple script can be caused in a lot of situation of coupling to become instability.Especially, in high dimension actual control system, when there is Parameter Perturbation and external disturbance in multiple input channel, cannot by realize Reasonable between input/output variable to and internal state between decoupling zero, set up the uncertain mathematics interact relation to system input and output and internal stability, and then the robust stability realizing non-matching uncertain multi-input multi-output system controls.
Summary of the invention
The object of the invention is to make controller cannot export the problem of the control signal of continuously smooth to solve existing non-singular terminal sliding-mode control existence buffeting, and the problem that Relative order is the control system of 1 cannot be applied to.
Be applicable to the level and smooth non-singular terminal sliding-mode control that Relative order is 1 control system, described control method is realized by following steps:
If step one controlled system is Relative order be 1 single input single output control system A:
in formula,
represent the state differential signal of single input single output control system A, x represents the system state of single input single output control system A, and u represents the single output control system A controlled quentity controlled variable of input, and t represents the time, then perform the control method of step 2; If controlled system is Relative order be 1 multi-input multi-output control system:
in formula,
represent multi-input multi-output control system state differential signal, x
srepresent multi-input multi-output control system state, u
srepresent multi-input multi-output control system controlled quentity controlled variable, t represents the time, then perform the control method of step 3;
Step 2, described single input single output control system:
control method be specially:
Step 2 one, due to single input single output control system A:
relative order be 1, described system state differential signal
be unknown quantity at real system, then utilize High-Order Sliding Mode robust precision differential device Real-time Obtaining
Step 2 two, design non-singular terminal sliding mode s (t), if
in formula, the design parameter c>0 of non-singular terminal sliding mode s (t); Design parameter p, q of non-singular terminal sliding mode s (t) are odd number, and meet p>q>0,1<p/q<2;
Step 2 three, introduces virtual controlling amount v (4), makes described system state x be 2 relative to the Relative order of virtual controlling amount v, add the Relative order of system; Based on Lyapunov Stability Theorem design virtual controlling amount v, ensure that system state arrives and maintains on sliding-mode surface s (t)=0 designed in advance, and to Parameter Perturbation and external disturbance, there is robustness, and due to integral action u=∫ vdt (5), make actual output controlled quentity controlled variable u smoothly continuous;
Step 3, described multi-input multi-output control system:
concrete control method be:
If the controlled Relative order of step 3 one is the multi-input multi-output control system of 1:
interior indeterminate or disturbance term are matching uncertainties, and namely described multi-input multi-output control system is coupling multi-input multi-output control system B:
If controlled system is Relative order be 1 non-matching uncertain multi-input multi-output control system C:
Δ B (t)=Bd (t) (8), in formula, d (t) ∈ R
m × mtime become bounded, its scope || d (t) ||≤l
d; F (t) represents the non-matching external disturbance of non-matching uncertain multi-input multi-output control system C, is R
nsmooth limited function, then perform the control method of step 3 two;
Step 3 two, described non-matching multi-input multi-output control system C:
concrete control method be:
Step 1, according to the controllability conditions r of non-matching multi-input multi-output control system C, first time nonsingular state transformation is done to non-matching multi-input multi-output control system C:
Y=F
1x (9), in formula, y represents the state after first time nonsingular state transformation, and x represents the system state before first time nonsingular state transformation, F
1∈ R
n × nfor transformation matrix; Then non-matching multi-input multi-output control system C is transformed to block control standard form control system D:
Step 2, block control standard form control system D step 1 obtained is written as block form, and corresponding plot control standard form control system D is decomposed into internal subsystems one and IOS one,
Internal subsystems one:
IOS one:
Controlled quentity controlled variable u only appears in IOS one; Y represents the state after first time nonsingular state transformation,
For the component of y, dimension is n
i, i=1 ... r, n
1+ ... + n
r=n;
for the component of f ',
for the component of f '; A '
i,jfor the design matrix of internal subsystems one, i=2 ... r;
Step 3, the internal subsystems one of existence step 2 obtained coupling does the nonsingular state transformation of second time:
Submatrix
Submatrix
for B
i, i-1moore-Penrose inverse, N
ifor transformation matrix F
2the design matrix of middle submatrix, corresponding plot control standard form control system D is converted to control system E further:
In formula, A "=(F
2)
-1a ' F
2represent coefficient matrices A ' through second time nonsingular state transformation after matrix;
Represent the matrix of matrix of coefficients B ' after the nonsingular state transformation of second time; Δ B " (t)=B " d (t) represents the matching uncertainties in the control system E multiple input path after the nonsingular state transformation of second time; F " (t)=(F
2)
-1f ' (t) represents the non-matching external disturbance of control system E after the nonsingular state transformation of second time;
Step 4, control system E step 3 obtained is write as block form, and correspondingly control system E is decomposed into internal subsystems two and IOS two:
Internal subsystems two:
IOS two:
In formula,
Represent the component of z;
Represent f " component i=1 ..., r,
represent f " component;
for the design matrix of internal subsystems two, i=2 ... r;
Step 5, the IOS two obtained step 4 is introduced reference model one and is asked for B to avoid applying classic method in step 4
1,0the process of d (t) inverse matrix, described reference model one is:
Beneficial effect of the present invention is:
First the inventive method judges the control system that Relative order is 1, then respectively by: be the single input single output control system of 1 for Relative order, utilize High-Order Sliding Mode robust precision differential device real-time acquisition system state differential, design non-singular terminal sliding mode, introducing virtual controlling amount utilizes the method for integral action to eliminate high frequency switching control item sgn (.) completely, the non-singular terminal sliding formwork realizing continuously smooth controls, reach the effect inherently eliminating the intrinsic buffeting problem of existing non-singular terminal sliding-mode control, the control signal of direct output continuously smooth.For coupling multi-input multi-output control system B, then the scalar way of realization of single input single output control system A is transformed to matrix vector way of realization, and level and smooth non-singular terminal sliding-mode control has unchangeability because of its robustness had to matching uncertainties disturbance.In addition, to non-matching uncertain multi-input multi-output control system C, then do twice nonsingular state transformation according to controllability conditions r, control system is resolved into r subsystem, then introduce reference model one, play and avoid application classic method to ask for B
1,0the process of d (t) inverse matrix, and there is the effect of the time-varying Hurst index eliminating IOS input channel, there is short-circuit evolution computation process and eliminate the effect of the intrinsic buffeting problem of existing non-singular terminal sliding-mode control, and then the robust stability that design control law and auxiliary control law realize control system adaptively controls.Having expanded existing non-singular terminal synovial membrane control method by said method is application in the control system of 1 at Relative order, has ubiquity application.Make existing non-singular terminal sliding-mode control correctly, to be effectively applied to Relative order be in 1 control system.
Accompanying drawing explanation
Fig. 1 to be embodiment 1 Relative order be 1 non-matching uncertain multi-input multi-output control system deviation differential
component signal convergence process figure, ordinate
expression state differential vector value, horizontal ordinate t represents the time, and unit is second (s);
Fig. 2 to be embodiment 1 Relative order be 1 non-matching uncertain multi-input multi-output control system non-singular terminal sliding mode s
1component convergence process figure, ordinate s
1represent non-singular terminal sliding mode vector value, horizontal ordinate t represents the time, and unit is second (s);
The component convergence process figure of Fig. 3 to be embodiment 1 Relative order be non-matching uncertain multi-input multi-output control system control law u of 1, ordinate represents the vector value of control law u, and horizontal ordinate t represents the time, and unit is second (s);
Fig. 4 to be embodiment 1 Relative order be 1 non-matching uncertain multi-input multi-output control system bias state e component at the process schematic of finite time convergence control to zero, ordinate represents system deviation state e vector value, horizontal ordinate t represents the time, and unit is second (s);
Fig. 5 to be embodiment 1 Relative order be 1 non-matching uncertain multi-input multi-output control system reference model state differential
component convergence process figure, ordinate
represent reference model state differential vector value, horizontal ordinate t represents the time, and unit is second (s);
Fig. 6 to be embodiment 1 Relative order be 1 non-matching uncertain multi-input multi-output control system non-singular terminal sliding mode s
2component convergence process figure, ordinate s
2represent non-singular terminal sliding mode vector value, horizontal ordinate t represents the time, and unit is second (s);
Fig. 7 to be embodiment 1 Relative order be 1 non-matching uncertain multi-input multi-output control system assist the component convergence process figure of control law v, ordinate represents the vector value of auxiliary control law v, and horizontal ordinate t represents the time, and unit is second (s);
Fig. 8 is embodiment 1 Relative order is IOS state z in the non-matching uncertain multi-input multi-output control system of 1
1component is at the process schematic of finite time convergence control to zero, and ordinate represents IOS state z
1vector value, horizontal ordinate t represents the time, and unit is second (s);
Fig. 9 to be embodiment 1 Relative order be 1 non-matching uncertain multi-input multi-output control system zero dy namics subsystem state z
2component and z
3convergence process figure, ordinate represents zero dy namics subsystem state z
2vector value and z
3convergency value, horizontal ordinate t represents the time, and unit is second (s);
The component convergence process figure of Figure 10 to be embodiment 1 Relative order be non-matching uncertain multi-input multi-output control system original state x component of 1, ordinate represents the vector value of system initial state x, and horizontal ordinate t represents the time, and unit is second (s).
Embodiment
Embodiment one: present embodiment be applicable to the level and smooth non-singular terminal sliding-mode control that Relative order is 1 control system, described control method is realized by following steps:
If step one controlled system is Relative order be 1 single input single output control system A:
in formula,
represent the state differential signal of single input single output control system A, x represents the system state of single input single output control system A, and u represents the single output control system A controlled quentity controlled variable of input, and t represents the time, then perform the control method of step 2; If controlled system is Relative order be 1 multi-input multi-output control system:
in formula,
represent multi-input multi-output control system state differential signal, x
srepresent multi-input multi-output control system state, u
srepresent multi-input multi-output control system controlled quentity controlled variable, t represents the time, then perform the control method of step 3;
Step 2, described single input single output control system:
control method be specially:
Step 2 one, due to single input single output control system A:
relative order be 1, described system state differential signal
be unknown quantity at real system, then utilize High-Order Sliding Mode robust precision differential device Real-time Obtaining
Step 2 two, design non-singular terminal sliding mode s (t), if
in formula, the design parameter c>0 of non-singular terminal sliding mode s (t); Design parameter p, q of non-singular terminal sliding mode s (t) are odd number, and meet p>q>0,1<p/q<2;
Step 2 three, introduces virtual controlling amount v (4), makes described system state x be 2 relative to the Relative order of virtual controlling amount v, add the Relative order of system; Based on Lyapunov Stability Theorem design virtual controlling amount v, ensure that system state arrives and maintains on sliding-mode surface s (t)=0 designed in advance, and to Parameter Perturbation and external disturbance, there is robustness, and due to integral action u=∫ vdt (5), make actual output controlled quentity controlled variable u smoothly continuous;
Step 3, described multi-input multi-output control system:
concrete control method be:
If the controlled Relative order of step 3 one is the multi-input multi-output control system of 1:
interior indeterminate or disturbance term are matching uncertainties, and namely described multi-input multi-output control system is coupling multi-input multi-output control system B:
If controlled system is Relative order be 1 non-matching uncertain multi-input multi-output control system C:
Step 3 two, described non-matching multi-input multi-output control system C:
concrete control method be:
Step 1, according to the controllability conditions r of non-matching multi-input multi-output control system C, first time nonsingular state transformation is done to non-matching multi-input multi-output control system C:
Y=F
1x (9), in formula, y represents the state after first time nonsingular state transformation, and x represents the system state before first time nonsingular state transformation, F
1∈ R
n × nfor transformation matrix; Then non-matching multi-input multi-output control system C is transformed to block control standard form control system D:
Step 2, block control standard form control system D step 1 obtained is written as block form, and corresponding plot control standard form control system D is decomposed into internal subsystems one and IOS one,
Internal subsystems one:
IOS one:
Controlled quentity controlled variable u only appears in IOS one; Y represents the state after first time nonsingular state transformation,
For the component of y, dimension is n
i, i=1 ... r, n
1+ ... + n
r=n;
for the component of f ',
for the component of f '; A '
i,jfor the design matrix of internal subsystems one, i=2 ... r;
Step 3, the internal subsystems one of existence step 2 obtained coupling does the nonsingular state transformation of second time:
Submatrix
Submatrix
for B
i, i-1moore-Penrose inverse, N
ifor transformation matrix F
2the design matrix of middle submatrix, corresponding plot control standard form control system D is converted to control system E further:
In formula, A "=(F
2)
-1a ' F
2represent coefficient matrices A ' through second time nonsingular state transformation after matrix;
Δ B " (t)=B " d (t) represents the matching uncertainties in the control system E multiple input path after the nonsingular state transformation of second time; F " (t)=(F
2)
-1f ' (t) represents the non-matching external disturbance of control system E after the nonsingular state transformation of second time;
Step 4, control system E step 3 obtained is write as block form, and correspondingly control system E is decomposed into internal subsystems two and IOS two:
Internal subsystems two:
IOS two:
In formula,
Represent the component of z;
Represent f " component i=1 ..., r,
represent f " component;
for the design matrix of internal subsystems two, i=2 ... r;
Step 5, the IOS two obtained step 4 is introduced reference model one and is asked for B to avoid applying classic method in step 4
1,0the process of d (t) inverse matrix, described reference model one is:
In formula,
represent the state of reference model one, ε (R
1represent the design parameter of reference model one, the control law of the reference model one of u ' expression design,
represent the auxiliary control law of the reference model one of design.
Design control law u (t), makes IOS two finite time follow the tracks of upper reference model one;
Design assistant control law w (t), makes reference model one finite time convergence control to zero.
Embodiment two:
With embodiment one unlike, present embodiment is applicable to the level and smooth non-singular terminal sliding-mode control that Relative order is 1 control system, in order to make IOS two finite time follow the tracks of upper reference model one, the design process of the control law u ' of reference model one described in the step 5 in step 3 two is specially:
First, deviation variables e=z is defined
1– ξ, can obtain by IOS two and reference model one bias system that Relative order is 1:
Make ε >>l
d, then || d (t)/(ε+1) || < < 1, namely has (I+d (t)/(ε+1)) ≈ I, so Relative order be 1 bias system be approximately:
Secondly, design non-singular terminal sliding mode one:
In formula,
represent non-singular terminal sliding mode one; c
1=diag (c11 ..., c1ni) and represent the design parameter of non-singular terminal sliding mode one non-singular terminal sliding mode one and c
1i>0, i=1 ... n
1; The design parameter p of non-singular terminal sliding mode one
1, q
1for odd number, and p
1>q
1>0,1<p
1/ q
1<2;
Finally, based on Relative order and sliding formwork equivalent control measurements, the level and smooth control law of design robust:
In formula,
U
eqt () represents equivalent control term, u
nt () represents switching control item, then introduce virtual controlling item:
based on Lyapunov Stability Theorem virtual controlling item v
1t () is designed to:
Namely by introducing virtual controlling item v
1and add the relative exponent number of system, and then obtain level and smooth continuous print working control rule u ', wherein, η
1>0 is ride gain, v
1eqand v
1nfor state variable, and convergence time is:
Its stability proof procedure: choose Lyapunov function
and first order derivative is asked to V (t), then have:
By equivalent control term u
eqt () substitutes into the approximate formula of bias system:
in, then have:
Thus, make:
Again according to virtual controlling item v
1t () design switches control formula, due to c
1for diagonal matrix, therefore
also be diagonal matrix,
When system does not also reach stable, || s
1(t) || ≠ 0, to above-mentioned inequality, and if only if
and
the transient state of system, explanation
state can not keep.
Embodiment three:
With embodiment one or two unlike, present embodiment be applicable to the level and smooth non-singular terminal sliding-mode control that Relative order is 1 control system, in order to make reference model one at finite time convergence control to zero and after bias system converges to zero described in step 5, reference model one is of equal value with IOS two, be zero by described bias system approximate formula, namely
Then:
substitute into reference model one and obtain reference model two:
First, non-singular terminal sliding mode two is designed to:
In formula,
represent non-singular terminal sliding mode two, c
2=diag (c
22..., c
2ni) represent non-singular terminal sliding mode two design parameter and c
2i>0, i=1 ... n
2; Non-singular terminal sliding mode two design parameter p
2, q
2for odd number and p
2>q
2>0,1<p
2/ q
2<2;
Secondly, for make reference model two system state ξ and
at convergence time
inside converge to zero, the Controller gain variations of reference model two is: w " (t)=w
eq(t)+w
nt () (27), in formula, equivalent control term is
Introduce virtual controlling item
to increase the relative exponent number of system, and actual switching control item w
n(t)=∫ v
2(t) dt continuously smooth; Based on Lyapunov Stability Theorem virtual controlling item v
2t () is designed to:
In formula, η
2>0 is ride gain.The design matrix N that in embodiment one, step 4 relates to
iwith non-matching external disturbance f "
uit () impact on internal system stability is described:
As IOS state z
1after converging to zero, internal subsystems one becomes the zero dy namics subsystem with disturbance:
According to Lyapunov Stability Theorem, N
i, i=2 ..., the eigenwert of r should be negative entirely, to ensure that described zero dy namics subsystem is at equilibrium point z
ithe local stability of=0, designs here:
Wherein ,-λ
2<...<-λ
r<0, makes z
2..., z
rorder convergence.
For described zero dy namics subsystem, there is symmetric positive definite matrix P
i, Q
i∈ R
ni × ni, i=2 ..., r meets
then z
2..., z
rasymptotic convergence is to Ω
iin:
In formula, τ
i=2 (|| f "
ui||+τ
i-1|| B
i, i-1||)/μ
i, μ
i=λ
min(Q
i)/λ
max(P
i), λ
min(Q
i) and λ
max(P
i) represent Q respectively
iand P
iminimum and maximum eigenwert; τ
i-1=0.
Asymptotic convergence region Ω
iderivation be: because N
i, i=2 ..., r eigenwert is designed to bear, to any positive definite symmetrical matrix Q
i, necessarily there is positive definite symmetrical matrix P
imeet
choose Lyapunov function
and first order derivative is asked to it, then have
According to Lyapunov Stability Theorem, then when
in,
then expression formula is derived:
Because
and twice nonsingular state transformation y=F
1x and
for nonsingular, thus former non-matching uncertain multi-input multi-output control system C can be obtained:
state x
i, i=1 ..., the convergence range of n:
Embodiment 1:
Design shape such as Relative order is the non-matching uncertain multi-input multi-output control system C of 1:
System initial state x (0)=[10.94,19.09,27.03,45.56,31.39 ,-16.68 ,-14.09]
t; Due to system controllability index r=3, therefore can be analyzed to 3 subsystem B
1,0, N
2and N
3be respectively:
Thus, according to first time nonsingular state transformation y=F
1x and the nonsingular state transformation of second time
described formula, twice transformation matrix F
1and F
2for:
Thus obtain the IOS after decoupling zero and be:
Internal subsystems is write as block form:
Controller gain variations parameter p
1=p
2=5, q
1=q
2=3, c
1=diag (0.3,0.3,0.3), c
2=diag (0.8,1,0.7), ε=10, η
1=20, η
2=80; High-Order Sliding Mode robust precision differential device described in step one is utilized to ask deviation
with
design parameter λ
0e1=22, λ
0e2=λ
0e3=15, λ
1e1=λ
1e2=λ
1e3=50; Ask reference model state differential
with
design parameter λ
0z1=22, λ
0z2=15, λ
0z3=13, λ
1z1=λ
1z2=λ
1z3=50, simulation result as shown in figures 1-8.
Fig. 1-4 and Fig. 5-8 adopts the control law u of level and smooth non-singular terminal sliding-mode control design and the simulation process of auxiliary control law v, and wherein Fig. 1 is system state differential
, Fig. 5 is system state differential
in Fig. 1, Fig. 5, solid line represents real differential signal, and dotted line represent by High-Order Sliding Mode robust precision differential device obtain differential signal, all can extract differential signal quickly and accurately; Fig. 2 and Fig. 6 is the non-singular terminal sliding mode s that the differential signal of application fetches is formed
1and s
2all can rapidly converge to zero; As Fig. 3 and Fig. 7 is respectively control law u and auxiliary control law v, signal smoothing is continuous, and essence eliminates high-frequency buffeting signal, and makes the IOS state z of bias state e and Fig. 8 as Fig. 4 respectively
1equal finite time convergence control to zero, thus demonstrates correctness and the validity of put forward level and smooth non-singular terminal sliding-mode control.
Fig. 9 is zero dy namics subsystem state z
2and z
3convergence curve, according to expression formula
utilize LMI tool box can obtain its convergence region to be || z
2||≤0.0710, || z
3||≤0.0533.
Figure 10 is the convergence curve of original system original state x, can obtain its maximum convergence range to be further || x||≤4.4946.The correctness that Relative order is the level and smooth non-singular terminal sliding-mode control of the control system of 1 is applicable to described in above simulation results show.
Claims (3)
1. be applicable to the level and smooth non-singular terminal sliding-mode control that Relative order is 1 control system, it is characterized in that: described control method is realized by following steps:
If step one controlled system is Relative order be 1 single input single output control system A:
in formula,
represent the state differential signal of single input single output control system A, x represents the system state of single input single output control system A, and u represents the single output control system A controlled quentity controlled variable of input, and t represents the time, then perform the control method of step 2; If controlled system is Relative order be 1 multi-input multi-output control system:
in formula,
represent multi-input multi-output control system state differential signal, x
srepresent multi-input multi-output control system state, u
srepresent multi-input multi-output control system controlled quentity controlled variable, t represents the time, then perform the control method of step 3;
Step 2, described single input single output control system:
control method be specially:
Step 2 one, due to single input single output control system A:
relative order be 1, described system state differential signal
be unknown quantity at real system, then utilize High-Order Sliding Mode robust precision differential device Real-time Obtaining
In formula, λ
0for design parameter one, λ
1for design parameter two; Y, v
0and v
1for the intermediate variable of High-Order Sliding Mode robust precision differential device state;
Step 2 two, design non-singular terminal sliding mode s (t), if
in formula, the design parameter c>0 of non-singular terminal sliding mode s (t); Design parameter p, q of non-singular terminal sliding mode s (t) are odd number, and meet p>q>0,1<p/q<2;
Step 2 three, introduces virtual controlling amount v (4), makes described system state x be 2 relative to the Relative order of virtual controlling amount v, add the Relative order of system; Based on Lyapunov Stability Theorem design virtual controlling amount v, ensure that system state arrives and maintains on sliding-mode surface s (t)=0 designed in advance, and to Parameter Perturbation and external disturbance, there is robustness, and due to integral action u=∫ vdt (5), make actual output controlled quentity controlled variable u smoothly continuous;
Step 3, described multi-input multi-output control system:
concrete control method be:
If the controlled Relative order of step 3 one is the multi-input multi-output control system of 1:
interior indeterminate or disturbance term are matching uncertainties, and namely described multi-input multi-output control system is coupling multi-input multi-output control system B:
Only the scalar way of realization of single input single output control system A in step 2 one to step 2 three need be transformed to matrix vector way of realization, and level and smooth non-singular terminal sliding-mode control has unchangeability because of its robustness had to matching uncertainties disturbance;
If controlled system is Relative order be 1 non-matching uncertain multi-input multi-output control system C:
In formula, x ∈ R
nfor system state, u ∈ R
mfor controlled quentity controlled variable, and 1≤m≤n; A ∈ R
n × nknown constant matrix, B ∈ R
n × mbe known constant matrix, dimension is n
1, n
1≤ m; (A, B) is controlled, if r is the controllability conditions of constant matrices (A, B) in non-matching uncertain multi-input multi-output control system C; Δ B (t) ∈ R
n × mrepresent the matching uncertainties in multiple input path, that is:
ΔB(t)=Bd(t) (8),
In formula, d (t) ∈ R
m × mtime become bounded, its scope || d (t) ||≤l
d; F (t) represents the non-matching external disturbance of non-matching uncertain multi-input multi-output control system C, is R
nsmooth limited function, then perform the control method of step 3 two;
Step 3 two, described non-matching multi-input multi-output control system C:
concrete control method be:
Step 1, according to the controllability conditions r of non-matching multi-input multi-output control system C, first time nonsingular state transformation is done to non-matching multi-input multi-output control system C:
y=F
1x (9),
In formula, y represents the state after first time nonsingular state transformation, and x represents the system state before first time nonsingular state transformation, F
1∈ R
n × nfor transformation matrix; Then non-matching multi-input multi-output control system C is transformed to block control standard form control system D:
In formula, A '=F
1aF
1 -1represent coefficient matrices A through first time nonsingular state transformation after matrix;
Represent the matrix of coefficients B matrix after first time nonsingular state transformation,
represent the submatrix of B ', dimension is n
1; Δ B ' (t)=B ' d (t) represents the matching uncertainties in multiple input path; F ' (t) represents the non-matching external disturbance of block control standard form control system D, and there is f ' (t)=F
1f (t);
Step 2, block control standard form control system D step 1 obtained is written as block form, and corresponding plot control standard form control system D is decomposed into internal subsystems one and IOS one,
Internal subsystems one:
IOS one:
Controlled quentity controlled variable u only appears in IOS one; Y represents the state after first time nonsingular state transformation,
for the component of y, dimension is n
i, i=1 ... r, n
1+ ... + n
r=n;
for the component of f ',
for the component of f '; A '
i, jfor the design matrix of internal subsystems one, i=2 ... r;
Step 3, the internal subsystems one of existence step 2 obtained coupling does the nonsingular state transformation of second time:
In formula, y is the state after first time nonsingular state transformation, and z is the state after the nonsingular state transformation of second time, F
2for transformation matrix,
Submatrix
Submatrix
for B
i, i-1moore-Penrose inverse, N
ifor transformation matrix F
2the design matrix of middle submatrix, corresponding plot control standard form control system D is converted to control system E further:
In formula, A "=(F
2)
-1a ' F
2represent coefficient matrices A ' through second time nonsingular state transformation after matrix;
Represent the matrix of matrix of coefficients B ' after the nonsingular state transformation of second time; Δ B " (t)=B " d (t) represents the matching uncertainties in the control system E multiple input path after the nonsingular state transformation of second time; F " (t)=(F
2)
-1f ' (t) represents the non-matching external disturbance of control system E after the nonsingular state transformation of second time;
Step 4, control system E step 3 obtained is write as block form, and correspondingly control system E is decomposed into internal subsystems two and IOS two:
Internal subsystems two:
IOS two:
In formula,
represent the component of z;
represent f " component i=1 ..., r,
represent f " component;
for the design matrix of internal subsystems two, i=2 ... r;
Step 5, the IOS two obtained step 4 is introduced reference model one and is asked for B to avoid applying classic method in step 4
1,0the process of d (t) inverse matrix, described reference model one is:
In formula,
represent the state of reference model one, ε ∈ R
1represent the design parameter of reference model one, the control law of the reference model one of u ' expression design,
represent the auxiliary control law of the reference model one of design.
2. be applicable to the level and smooth non-singular terminal sliding-mode control that Relative order is 1 control system according to claim 1, it is characterized in that: in order to make IOS two finite time follow the tracks of upper reference model one, the design process of the control law u ' of reference model one described in the step 5 in step 3 two is specially:
First, deviation variables e=z is defined
1– ξ, can obtain by IOS two and reference model one bias system that Relative order is 1:
Make ε > > l
d, then || d (t)/(ε+1) || < < 1, namely has (I+d (t)/(ε+1)) ≈ I, so Relative order be 1 bias system be approximately:
Secondly, design non-singular terminal sliding mode one:
In formula,
represent non-singular terminal sliding mode one; c
1=diag (c11 ..., c1ni) and represent the design parameter of non-singular terminal sliding mode one non-singular terminal sliding mode one and c
1i>0, i=1 ... n
1; The design parameter p of non-singular terminal sliding mode one
1, q
1for odd number, and p
1>q
1>0,1<p
1/ q
1<2;
Finally, based on Relative order and sliding formwork equivalent control measurements, the level and smooth control law of design robust:
u′(t)=u
eq(t)+u
n(t) (19),
In formula,
U
eqt () represents equivalent control term, u
nt () represents switching control item, then introduce virtual controlling item:
Based on Lyapunov Stability Theorem virtual controlling item v
1t () is designed to:
Namely by introducing virtual controlling item v
1and add the relative exponent number of system, and then obtain level and smooth continuous print working control rule u ', wherein, η
1>0 is ride gain, v
1eqand v
1nfor state variable, and convergence time is:
In formula, e
i(0) be deviation variables e
ithe original state value of (t).
3. according to claim 1 or 2, be applicable to the level and smooth non-singular terminal sliding-mode control that Relative order is 1 control system, it is characterized in that: in order to make reference model one at finite time convergence control to zero and after bias system converges to zero described in step 5, reference model one is of equal value with IOS two, be zero by described bias system approximate formula, namely
substitute into reference model one and obtain reference model two:
Relative order due to reference model two is 1, with the level and smooth non-singular terminal sliding mode controller design process of the step 2 controlled quentity controlled variable u in step 3 two, direct Design assistant control law w ' (t), then described in step 5, the design process of reference model one control law w ' (t) is specially:
First, non-singular terminal sliding mode two is designed to:
In formula,
represent non-singular terminal sliding mode two, c
2=diag (c
22..., c
2ni) represent non-singular terminal sliding mode two design parameter and c
2i>0, i=1 ... n
2; Non-singular terminal sliding mode two design parameter p
2, q
2for odd number and p
2>q
2>0,1<p
2/ q
2<2;
Secondly, for make reference model two system state ξ and
at convergence time
inside converge to zero, the Controller gain variations of reference model two is: w " (t)=w
eq(t)+w
nt () (27), in formula, equivalent control term is
Introduce virtual controlling item
to increase the relative exponent number of system, and actual switching control item w
n(t)=∫ v
2(t) dt continuously smooth; Based on Lyapunov Stability Theorem virtual controlling item v
2t () is designed to:
In formula, η
2>0 is ride gain.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410578260.7A CN104267605B (en) | 2014-10-25 | 2014-10-25 | It is applicable to the smooth non-singular terminal sliding-mode control that Relative order is 1 control system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410578260.7A CN104267605B (en) | 2014-10-25 | 2014-10-25 | It is applicable to the smooth non-singular terminal sliding-mode control that Relative order is 1 control system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104267605A true CN104267605A (en) | 2015-01-07 |
CN104267605B CN104267605B (en) | 2016-10-05 |
Family
ID=52159133
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410578260.7A Active CN104267605B (en) | 2014-10-25 | 2014-10-25 | It is applicable to the smooth non-singular terminal sliding-mode control that Relative order is 1 control system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104267605B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104614995A (en) * | 2015-02-16 | 2015-05-13 | 天津大学 | General design method for second-order system limited time slip form controller |
CN105186863A (en) * | 2015-08-14 | 2015-12-23 | 哈尔滨工业大学 | Continuous nonsingular terminal sliding mode control method based on Buck converter |
CN107147283A (en) * | 2017-06-29 | 2017-09-08 | 江苏大学 | A kind of Buck converter compensating control methods based on disturbance observer and Second Order Sliding Mode |
CN108833075A (en) * | 2018-06-21 | 2018-11-16 | 东北大学 | A kind of second order Projective synchronization in chaotic method based on non-singular terminal sliding mode controller |
CN113517835A (en) * | 2021-04-22 | 2021-10-19 | 湖南工业大学 | PMSM drive system field loss fault control method and permanent magnet synchronous motor |
CN115818490A (en) * | 2022-11-22 | 2023-03-21 | 哈尔滨工程大学 | Semi-active heave compensation device for ROV release and recovery |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4358822A (en) * | 1976-08-04 | 1982-11-09 | Sanchez Juan M | Adaptive-predictive control system |
CN102403719A (en) * | 2011-11-22 | 2012-04-04 | 上海电力学院 | Method for designing sliding mode load frequency controller of non matched indefinite power system |
CN102709910A (en) * | 2012-06-01 | 2012-10-03 | 河海大学常州校区 | Method for controlling self-adpative sliding mode of single-phase SAPF (Shunt Active Power Filter) |
CN103077729A (en) * | 2013-02-05 | 2013-05-01 | 西安电子科技大学 | High-order nonsingular terminal slip mode control method for two-stage magnetic head positioning system |
CN103728882A (en) * | 2014-01-07 | 2014-04-16 | 河海大学常州校区 | Adaptive inversion nonsingular terminal sliding mode control method of micro gyroscope |
-
2014
- 2014-10-25 CN CN201410578260.7A patent/CN104267605B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4358822A (en) * | 1976-08-04 | 1982-11-09 | Sanchez Juan M | Adaptive-predictive control system |
CN102403719A (en) * | 2011-11-22 | 2012-04-04 | 上海电力学院 | Method for designing sliding mode load frequency controller of non matched indefinite power system |
CN102709910A (en) * | 2012-06-01 | 2012-10-03 | 河海大学常州校区 | Method for controlling self-adpative sliding mode of single-phase SAPF (Shunt Active Power Filter) |
CN103077729A (en) * | 2013-02-05 | 2013-05-01 | 西安电子科技大学 | High-order nonsingular terminal slip mode control method for two-stage magnetic head positioning system |
CN103728882A (en) * | 2014-01-07 | 2014-04-16 | 河海大学常州校区 | Adaptive inversion nonsingular terminal sliding mode control method of micro gyroscope |
Non-Patent Citations (3)
Title |
---|
WANG YAN-MIN等: "Trajectory Tracking Control of Uncertain Two-link Flexible Manipulators", 《***仿真学报》 * |
冯勇等: "非奇异终端滑模控制***的设计方法", 《控制与决策》 * |
郑剑飞等: "非匹配不确定多变量***高阶终端滑模控制", 《电机与控制学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104614995A (en) * | 2015-02-16 | 2015-05-13 | 天津大学 | General design method for second-order system limited time slip form controller |
CN104614995B (en) * | 2015-02-16 | 2017-03-29 | 天津大学 | A kind of general design method of second-order system finite time sliding mode controller |
CN105186863A (en) * | 2015-08-14 | 2015-12-23 | 哈尔滨工业大学 | Continuous nonsingular terminal sliding mode control method based on Buck converter |
CN107147283A (en) * | 2017-06-29 | 2017-09-08 | 江苏大学 | A kind of Buck converter compensating control methods based on disturbance observer and Second Order Sliding Mode |
CN107147283B (en) * | 2017-06-29 | 2019-10-01 | 江苏大学 | A kind of Buck converter compensating control method based on disturbance observer and Second Order Sliding Mode |
CN108833075A (en) * | 2018-06-21 | 2018-11-16 | 东北大学 | A kind of second order Projective synchronization in chaotic method based on non-singular terminal sliding mode controller |
CN108833075B (en) * | 2018-06-21 | 2021-06-11 | 东北大学 | Second-order chaotic projection synchronization method based on nonsingular terminal sliding mode controller |
CN113517835A (en) * | 2021-04-22 | 2021-10-19 | 湖南工业大学 | PMSM drive system field loss fault control method and permanent magnet synchronous motor |
CN113517835B (en) * | 2021-04-22 | 2023-06-06 | 湖南工业大学 | PMSM driving system loss-of-magnetic fault control method and permanent magnet synchronous motor |
CN115818490A (en) * | 2022-11-22 | 2023-03-21 | 哈尔滨工程大学 | Semi-active heave compensation device for ROV release and recovery |
CN115818490B (en) * | 2022-11-22 | 2023-07-21 | 哈尔滨工程大学 | Semi-active heave compensation device for ROV release recovery |
Also Published As
Publication number | Publication date |
---|---|
CN104267605B (en) | 2016-10-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104267605A (en) | Smooth nonsingular terminal sliding-mode control method suitable for control system with relative degree of 1 | |
CN106886149B (en) | A kind of spacecraft robust finite time saturation Attitude tracking control method | |
CN106843254B (en) | It is a kind of actively to reconstruct fault tolerant control method in real time | |
CN106547207B (en) | Construction method of nonlinear multi-input multi-output system hybrid observer | |
CN103728882B (en) | The self-adaptation inverting non-singular terminal sliding-mode control of gyroscope | |
CN108733031B (en) | Network control system fault estimation method based on intermediate estimator | |
CN104678765A (en) | Piezoelectric ceramic actuator hysteretic model and control method thereof | |
CN102929142A (en) | Method for designing controller of time-varying aircraft model with uncertainty | |
CN107679287B (en) | Electromagnetic transient numerical calculation method based on 3-step 4-order implicit Taylor series method | |
Guan et al. | Improved H∞ filter design for discrete-time Markovian jump systems with time-varying delay | |
CN106227964B (en) | Nonlinear systems with hysteresis parameter identification method based on extended state observer | |
Sadabadi et al. | An LMI formulation of fixed-order H∞ and H 2 controller design for discrete-time systems with polytopic uncertainty | |
CN109314390A (en) | Obtain the global linear symmetry approach of equivalent conduction compensation type of direct-current mains trend | |
CN105186863A (en) | Continuous nonsingular terminal sliding mode control method based on Buck converter | |
CN103795071A (en) | Wide-area damping control method based on identification of power system model | |
Verriest et al. | Advances in optimal control of differential systems with the state suprema | |
CN102929144A (en) | Method for designing controller of interval aircraft model | |
CN105549385A (en) | Decoupling internal die controller, control system, and control method of multivariable time-lag non-minimum-phase non-square system | |
Qi et al. | Observer design for stochastic time‐delayed Markovian jump systems with incomplete transition rates and actuator saturation | |
CN102929129B (en) | Aircraft controller designing method containing multiple time lag varying models | |
Zhao et al. | Finite-time stochastic input-to-state stability and observer-based controller design for singular nonlinear systems | |
Liu et al. | An adaptive sliding mode control of delta operator systems with input nonlinearity containing unknown slope parameters | |
Sun et al. | Robust actuator and sensor fault estimation for fuzzy delay singularly perturbed systems via proportional multiple-integral observer | |
Rafiq et al. | Structure preserving nonlinear reduced order modeling technique for power systems | |
CN102929135B (en) | Method for designing controller for aircraft time lagging and varying model |
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 |