Specific embodiment
The technical scheme in the embodiments of the invention will be clearly and completely described below, it is clear that described implementation
Example is only a part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, this field is common
Technical staff's every other embodiment obtained without making creative work belongs to the model that the present invention protects
It encloses.
The present invention provides the controller design method of cluster satellite system, this method mainly comprise the steps that firstly,
Formation modeling is carried out to cluster satellite system described in Hill's equation using formation error analysis technique, by cluster satellite system
Formation problem is converted into Linear Higher-Order multi-agent system consistency problem;Secondly, when using specified in the theory of optimal control
Between design method to Linear Higher-Order designing Multi-Agent system substep control strategy;Then, using motion planning technology to substep
Discrete time in control strategy clicks through professional etiquette and draws, and requires to be projected into finite time interval according to specified time;Finally,
It is theoretical with series and stability analysis, to realize the control target of the specified time consistency formation of cluster satellite system.
Specifically, the controller design method process of above-mentioned cluster satellite system are as follows:
Cluster satellite system is described using Hill's equation, establishes the model of cluster satellite system, as follows:
Wherein, ri=[rxi,ryi,rzi]TIt is the three-dimensional location coordinates of i-th of subsatellite, uiFor formation control device, N is son
The quantity of satellite, I3For three-dimensional unit matrix, nr=7.273 × 10-5s-1For the intrinsic frequency of reference orbit, orbit radius R0
=4.224 × 107m。
It is how intelligent that by cluster satellite system formation control problem high-grade linear is converted by the design method of formation error
The consistency control problem of system system.It enablesCluster satellite system can be write as general linear system:
Formation control rule can be written as follow specified time consistency:
WhereinNiFor the set of other subsatellites composition neighbouring with subsatellite i, InFor n
Dimension unit matrix, 0nNull matrix is tieed up for n.Sampling time sequence are as follows:
Wherein, Ts> 0 is to stablize the time according to offline preassign of mission requirements,It is to have multinomial receipts
Hold back the infinite power series { δ of speedk, i=1,2 ... } sum.Particularly, it can choose a special example such as
According to method for optimally controlling, the specified time consistency formation control rule of subsatellite i is given below:
tk≤t<tk+1, i=1,2 ..., N
Wherein,hiFor
Given formation configuration, tk≤t<tk+1, i=1,2 ..., N.
The design of above-mentioned specified time consistency controller allows for following Hamiltonian function:
It is wherein association's state.Then, the cost function of the corresponding Hamiltonian function of formula (4) is as follows:
Wherein it is possible to be seen as initial and terminal time respectively.So, formula (1) is write as according to method for optimally controlling:
In addition, according to extremum conditions:Have:
ui(t)=BTpi(t) (7)
Therefore, the determination of optimal controller is attributed to calculating p in formula (7)i(t)。
By obtaining in formula (7) substitution formula (6):
By above-mentioned equation from tkIt is integrated to tk+1, and then have:
Due toThen
So:
Wherein
Next design [tk,tk+1] terminal condition of formula (4) in the time, it is as follows:
Formula (11) are substituted into formula (9), are obtained:
So, if Φ is reversible, followed by have:
Therefore, for time series tk, distributed director has been obtained, such as formula (2).
It summarizes above-mentioned, there is the specified time control controller of formula (2), the state x in formula (1)i(tk) can be controlled to
The average state of prediction,In time interval [tk,tk+1] in.Intuitively,
By after steps multiple enough, the state of all intelligent bodies of system is by compliance in formula (1).
The existence and feasibility of the controller obtained below to present invention design are analyzed.
1, controller Analysis of Existence: note that only when Φ can the inverse time, above controller just exists.Therefore, in the case where continuing
Before the work in face, following analysis conclusion is given.
It can be obtained by (10):
Note
Pass through reduction to absurdity, it is assumed that ∏ (t) is unusual.Accordingly, there exist at least one non-vanishing vector α ∈ Rn, so that:
αT∏ (t)=0 (12)
For any t ∈ R.K order derivative (k=1,2 ..., n-1) of modus ponens (12) both sides about t, has:
For any t ∈ R.Due to the arbitrariness of t, equation (13) is set up in t=0.Therefore, there is following equation establishment:
Simplify above-mentioned equation, there is αTAk-1BBT=0, k=1,2 ... n, therefore αTBBT=0, αTABBT=0, αTA2BBT=0,
αTA3BBT=0 ..., αTAn-1BBT=0.
Remember Q=[BBT ABBT A2BBT…An-1BBT], then αTQ=0.
By α ≠ 0 it is found that matrix Q is linearly related.Note that when (A, B) is controllable, (A, BBT) be can
Control.This condition controllable with (A, B) contradicts.Therefore, it was demonstrated that ∏ be it is nonsingular, i.e., Φ is reversible.
Similarly, it is assumed that (A, B) is uncontrollable, then at least there is a non-vanishing vector β ∈ Rn, so that βTQ=0, i.e. βTBBT=0, βTABBT=0, βTA2BBT=0, βTA3BBT=0 ..., βTAn-1BBT=0.So βT∏ (t)=0.This is with Φ can
Inverse contradiction.Therefore, (A, B) is controllable.
So this is wanting substantially for linear system control if (A, B) is controller presence that is controllable, then being proposed
It asks.
2, controller feasibility analysis: first, it was demonstrated that in sampling time sequence { tkUnder, the system mode in formula (1) exists
Consistency can be realized under the control of formula (2) controller.Formula (2) are substituted into formula (1) and are obtained:
By (15) from tkIt is integrated to tk+1, k=0,1 ..., available:
NoteThen:
WhereinThen:
There is a spanning tree in digraph G.Random matrix IN- NL has an eigenvalue λ1=1 corresponding algebraic multiplicity is equal to 1,
Other all characteristic values meet | λi| < 1, i=1,2 ..., N, then, for matrix IN- NL, there are a column vector ξ to make:
There is t additionally, due to as k → ∞k→Ts, then there is tk-1-t0It is bounded, this makes matrixIt is every
One is bounded.Then have:
NoteHave
Then
It notices againSo as k → ∞, i.e. limk→∞||xi(tk)-xj
(tk) | |=0, discrete state will reach consistency with exponential rate.
Secondly, stablizing time T for preassigning offlines, it was demonstrated that work as tk→TsWhen discrete state xi(tk) it can reach specified
Time consistency.According to the convergent sample sequence S of multinomial rate, there is limk→∞tk=Ts.ThereforeSo, when stablizing for preassigning
Between Ts, work as tk→TsWhen, discrete state xi(tk) specified time consistency will be realized with exponential rate.
Finally, will demonstrate that in t → TsWhen continuous state xi(t) it can achieve specified time consistency.By by (15) from tk
It is integrated to t, available:
tk≤t<tk+1
Note
Then
Separately have
It is furthermore noted that:
Due to the length of time intervalIt is bounded above.HaveIt is bounded.It is controllable additionally, due to (A, B), it is known that Φ is reversible,Wherein | Φ | ≠ 0.It is assumed thatWherein fs, s=0,1 ... is coefficient, at least
There are a limited constant m to make coefficient fm≠ 0, f0=f1=...=fm-1=0.So work as tk+1-tkWhen sufficiently small, | Φ | it can
To be write as again | Φ |=fm(tk+1-tk)m+o(tk+1-tk), wherein o (tk+1-tk) it is tk+1-tkHigher-order shear deformation.Due toIt is multinomial rate, hasAccording to limk→∞tk=TS,It follows thatSimilarly, availableTherefore have
In summary, under the conditions of the controller of formula (2), the linear multi-agent system of formula (1) can stablize the time preassigning
TsLower realization consistency.
By the modeling, design, conversion of above-mentioned process and analysis it can be concluded that, the specified time consistency that this patent proposes
Formation control device can solve the specified time formation control problem of cluster satellite system.
The present invention also provides the control method of cluster satellite system, this method is using the controller in formula (2) to defending
Star group system is controlled.
Although preferred embodiments of the present invention have been described, it is created once a person skilled in the art knows basic
Property concept, then additional changes and modifications may be made to these embodiments.So it includes excellent that the following claims are intended to be interpreted as
It selects embodiment and falls into all change and modification of the scope of the invention.
Obviously, various changes and modifications can be made to the invention without departing from essence of the invention by those skilled in the art
Mind and range.In this way, if these modifications and changes of the present invention belongs to the range of the claims in the present invention and its equivalent technologies
Within, then the present invention is also intended to include these modifications and variations.