CN108536014A - A kind of model predictive control method for considering the spacecraft attitude of flywheel dynamic characteristic and evading - Google Patents
A kind of model predictive control method for considering the spacecraft attitude of flywheel dynamic characteristic and evading Download PDFInfo
- Publication number
- CN108536014A CN108536014A CN201810298997.1A CN201810298997A CN108536014A CN 108536014 A CN108536014 A CN 108536014A CN 201810298997 A CN201810298997 A CN 201810298997A CN 108536014 A CN108536014 A CN 108536014A
- Authority
- CN
- China
- Prior art keywords
- flywheel
- spacecraft
- indicate
- model
- attitude
- 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
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000005457 optimization Methods 0.000 claims abstract description 33
- 238000013178 mathematical model Methods 0.000 claims abstract description 8
- 238000013461 design Methods 0.000 claims abstract description 5
- 238000012545 processing Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 23
- 230000014509 gene expression Effects 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000002360 preparation method Methods 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 abstract description 6
- 238000003786 synthesis reaction Methods 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000026676 system process Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
Landscapes
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Artificial Intelligence (AREA)
- Aviation & Aerospace Engineering (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
A kind of model predictive control method for considering the spacecraft attitude of flywheel dynamic characteristic and evading, includes the following steps:The dynamics of attitude dynamics and flywheel based on spacecraft establishes prediction model;Secondly the mathematical model of constraint is established according to the dynamic characteristic and the instrument angle of sight of flywheel;Then the performance index function towards different task demand is designed, control problem is converted into seek the extreme-value problem of object function under the conditions of equation and inequality constraints;Finally, pass through the optimization method based on real-time iterative, the rapid solving above problem, this method can be good at handling posture evasion of the spacecraft under executing agency's restraint condition, it is optimal by reaching the synthesis of energy and time to the design of object function, and by the processing of real-time iterative and hot start method, and the calculation amount of solving-optimizing problem can be reduced.
Description
Technical field
The present invention relates to technical field of spacecraft control is belonged to, it is mainly used in the spacecraft controlled using counteraction flyback
Attitude maneuver evades control, and in particular to a kind of Model Predictive Control for considering the spacecraft attitude of flywheel dynamic characteristic and evading
Method.
Background technology
The development of space technology is at full speed in recent years, and the mission requirements of in-orbit spacecraft are also more and more, so in spacecraft
On all equipped with various optical instruments, such as CCD camera, infrared interferometer, be required for making during these instruments work
Its sight avoids directly opposite strong light, to protect in instrument to illumination and the more sensitive component of temperature ratio.So this requires
Spacecraft makes the direction of these instruments that should get around the direction of strong light during attitude maneuver.Meanwhile most of spacecraft appearances
Executing agency in state control system combines for counteraction flyback, and there is counteraction flyback control accuracy height, output torque only to need
The advantages that consuming electric energy.But counteraction flyback is easy to reach full in use there is also the moment of reaction is smaller
And the problems such as, to influence the gesture stability of spacecraft.So research spacecraft is under the constraints of executing agency's performance
Posture evasion meaning just it is particularly significant.
The method evaded about posture can have the potential function method mentioned in patent 201710521561.X, potential function method can
Effectively to evade constraint, but using the method for Model Predictive Control, it can not only consider that posture evades constraint, and can
Consider the performance constraints of executing agency, can also realize the optimization of energy and time according to performance function, be more suitable for holding in addition
Attitude maneuver control under the constraint of row mechanism performance.
Invention content
It is an object of the invention to overcome the deficiencies of the prior art and provide a kind of spacecraft appearances considering flywheel dynamic characteristic
The model predictive control method that state is evaded, existing executing agency's performance constraints and appearance during the attitude maneuver of in-orbit spacecraft
The problem of modal constraint, the present invention provide a kind of Model Predictive Control side for considering the spacecraft attitude of flywheel dynamic characteristic and evading
Method.This method be it is a kind of can handle the performance constraints of executing agency simultaneously and posture evades model predictive control method, it
Real-time iterative and the method for thermal starting have been used in optimization process to improve the solution efficiency for solving optimization problem in MPC, and have been set
Performance function has been counted to realize that energy consumption and the synthesis of time are optimal during attitude maneuver.
The present invention provides a kind of model predictive control method for considering the spacecraft attitude of flywheel dynamic characteristic and evading, packets
Include following steps:
(1) include according to the attitude kinematics of spacecraft and the foundation of the rotational power of kinetic characteristics and flywheel characteristic
Prediction model of the spacecraft attitude model of executing agency as MPC;
(2) it carries the angle of sight of instrument according to spacecraft and evades vector and establish the mathematical model that posture evades constraint,
And the mathematical model of executing agency's constraint is established according to the angular momentum of flywheel saturation and torque saturated characteristic;
(3) corresponding optimization object function is designed according to mission requirements, including the quadratic form of executing agency's input and being
The quadratic form of system state error executes time and energy expenditure to consider;
(4) control problem is converted into the equation and state of system dynamics equation and the inequality constraints of input-bound
Under the conditions of, then the problem of seeking object function extreme value, uses the optimized treatment method rapid solving of real-time iterative, optimization is obtained
It solves and is exported as the controlled quentity controlled variable of system.
It is by the way that the kinetics equation of flywheel is attached to boat in the spacecraft attitude dynamics model described in step (1)
In its device attitude dynamic equations, and sliding-model control arrives, and representation is as follows:
Wherein, ω=[ω1,ω2,ω3]TIndicate the attitude angle of spacecraft relative inertness coordinate system under body coordinate system
Velocity vector, ω1,ω2,ω3Respectively spacecraft is about the angular speed on roll axis, yaw axis and the pitch axis in this system
Component;Indicate derivatives of the ω to the time;J is the total inertia matrix of spacecraft, is simplified shown as diagonal matrix J=diag (J1,J2,
J3), J1,J2,J3For around the rotary inertia of the principal axis of inertia;S (ω) is skew symmetric matrix, and form isτ indicates the output torque of executing agency;Q=[q0,q1,q2,q3]TIndicate the posture list of spacecraft
Position quaternary number,Indicate that the scalar component in the posture unit quaternion of spacecraft, θ expressions are crossed around Euler's shaft rotation
One angle,ex,ey,ezThe rotary shaft on three directions of Euler's axis is represented,
And meet Indicate derivatives of the q to the time;Ω (ω) is skew symmetric matrix, and form is as follows:
The model of counteraction flyback combination is as follows:
Wherein HrwFor the angular momentum of flywheel,For the derivative of the angular momentum relative time of flywheel, by four
The relationship of its angular momentum and rotating speed is as follows in the flywheel of flywheel composition:
Hrw=CJrwN
Wherein C is that 3 × 4 flywheels install matrix, N=[n1,n2,n3,n4]TFor the angular velocity vector of flywheel, n1,n2,n3,n4
The angular speed of each flywheel is indicated respectively;JrwIndicate the moment of inertia matrix of flywheel, form Jrw=JαI4×4, JαTable
Show the rotary inertia of single flywheel, I4×4For 4 rank unit matrixs;
Since executing agency is to control posture by way of exchanging angular momentum with spacecraft, so the angular motion that system is total
Measure conservation:
H=Hrw+Jω
H is the total angular momentum of system, is a constant in no moment of face.
The model integration of the model of executing agency and spacecraft attitude dynamics is got up to obtain:
By above-mentioned model discretization, the setting sampling interval is that Δ t is obtained, and is carved in kth:
ωk+1=J-1ΔtS(ωk)H+J-1JαCΔNk+ωk
Wherein subscript k indicates the value to dependent variable at the kth moment, Δ Nk=Nk-Nk-1, I3×3For three rank unit matrixs.
In step (2):
(a) counteraction flyback maximum output torque constrains:
The torque output of counteraction flyback is to change the angular momentum realization of flywheel, so there is following form to state:
Wherein TmaxFor maximum output torque vector.It will be obtained after above-mentioned formula discretization:
It is obtained after arrangement:
(b) counteraction flyback maximum angular momentum constrains:
The angular momentum saturation of flywheel is presented as that the rotating speed of flywheel rotor reaches the upper limit, so the full constraint of angular momentum can be used and fly
Angular speed constraint is taken turns to indicate:
-Nmax≤ΔNk+Nk-1≤Nmax
Wherein NmaxFor flywheel maximum angular rate vector.
It is obtained after arrangement:
(c) angle of sight constraint that spacecraft attitude is directed toward:
Consider that the direction of spacecraft will evade certain conical lines-of-sight areas, the following form of posture restraint of design:
Wherein α indicates that the unit vector that spacecraft is directed toward under body coordinate system, β indicate spacecraft under body coordinate system
The unit vector in direction need to be evaded, θ indicates the size for evading the angle of sight in region.
The optimization object function function V (x in step (3)k,uk) be expressed as:
Wherein xk=[qk,ωk]TThe quantity of state of expression system, physical significance are the posture and angular speed of spacecraft;With
Flywheel rotation speed change amount is as input, i.e. uk=Δ Nk;NPFor the estimation range of MPC;Q and P is that state variable and input become
The weight matrix of amount illustrates that optimization aim is more concerned with the stable time, if Q is larger relative to P if Q relative to P if larger
Illustrate that optimization aim is more concerned with the consumption of energy;
It is organized into compact form:
Wherein
In step (4), control problem is converted to the extreme-value problem for seeking object function under constraints, is used in combination in real time
The method of iteration solves controlled quentity controlled variable.The problem of control input will be solved, is converted to following mathematical problem:
xi,k+1=F (xi,k,ui,k)
L(xi,k,ui,k)≤017×1K=0,1 ..., NP-1
WhereinIndicate the optimization problem in moment i,Indicate the feedback of status in current time i system, xi,k
It indicates by current time state xi,0The system mode at the kth moment of reckoning;ui,kThe estimated input at expression kth moment, 017×1Table
Show 17 × 1 null matrix;F(xi,k,ui,k) it is system state equation, it is embodied as:
L(xi,k,ui,k) it is system restriction, it is embodied as:
Processing procedure using the method for real-time iterative is:Known previous moment optimization problemIt solves
It finds outWithAnd then solving-optimizing problemIt is specifically divided into
Two stages:Preparation stage and response phase:
Preparation phase procedures are:By xi-1And ui-1The displacement for doing a sampling instant keeps the last one element constant, obtains
It arrivesWithAccording toWithCalculate sensitive matrix Ai,k,
Bi,k,Ci,k,Di,kWith error li,k,ri,k, expression is:
By linearisation, the optimization problem of a nonlinear restriction is converted to the quadratic programming problem of linear restriction:
Δui,k=ui,k-ui-1,k
Δxi,k+1=Ai,kΔxi,k+Bi,kΔui,k+ri,k
Ci,kΔxi,k+Di,kΔui,k+li,k≤ 0 k=0,1 ..., NP-1
The process of response phase is:The feedback of status of acquisition systemBring Solve problems intoObtain Δ xi,kAnd Δ
ui,kPass through formula:
Obtain problemSolution (xi,k,ui,k), finally again by ui=ui,0Control as system inputs.
The model predictive control method that the spacecraft attitude of the consideration flywheel dynamic characteristic of the present invention is evaded, can make space flight
Device in the case of executing agency's limited capacity, carries out posture and evades, while can meet energy and time during attitude maneuver
Synthesis it is optimal.The optimization method provided can improve the solution efficiency of MPC, and accelerating solving-optimizing using the thought of thermal starting asks
The speed of topic, the thought that real-time iterative is divided into two stages can make output for the feedback of status at a nearest moment, carry
High system executes the timeliness of control.The present invention has following advantages in summary:
(1) by the model of spacecraft together with the model integration of executing agency, using the knots modification of the rotating speed of flywheel as
Input, control moment is distributed to the control method of each flywheel by control moment again with traditional first obtaining, controlled quentity controlled variable is directly made
For executing agency, the process of the distribution of controlled quentity controlled variable is eliminated, directly whole system process is optimized.
(2) it compared with the potential function method that existing spacecraft attitude is evaded, can will be held using model predictive control method
Posture directing constraint in the performance constraints and attitude of satellite mobile process of row mechanism considers simultaneously, in addition can also be according to task
Corresponding performance function is arranged in demand so that the process of attitude maneuver realizes that energy and the synthesis of time are optimal.
(3) compared with general posture evades optimum control, the method using Model Predictive Control is closed-loop control,
Track optimizing again can be made according to feedback of status after making execution action each time so that control has robustness.
(4) the optimization thought for the real-time iterative and thermal starting that Model Predictive Control of the invention uses, compared to common
The method of Nonlinear Model Predictive Control can reduce operand and optimization efficiency, can be directed to the system shape at a nearest moment
State responds, and the output of control is more time-efficient.
Description of the drawings
Fig. 1 is the flow diagram for the model predictive control method for considering that the spacecraft attitude of flywheel dynamic characteristic is evaded;
Fig. 2 procedural block diagrams that system solution control inputs in order to control;
Specific implementation mode
The following detailed description of the specific implementation of the present invention, it is necessary to it is indicated herein to be, implement to be only intended to this hair below
Bright further explanation, should not be understood as limiting the scope of the invention, and field technology skilled person is according to above-mentioned
Some nonessential modifications and adaptations that invention content makes the present invention, still fall within protection scope of the present invention.
The present invention provides a kind of model predictive control methods for considering the spacecraft attitude of flywheel dynamic characteristic and evading, such as
Shown in Fig. 1, the specific steps are:First according to the rotational power of the attitude kinematics and kinetic characteristics and flywheel of spacecraft
Characteristic establishes prediction model of the spacecraft attitude model comprising executing agency as MPC;Secondly instrument is carried according to spacecraft
The angle of sight and evade vector and establish the mathematical model that posture evades constraint, and according to the angular momentum of flywheel saturation and torque
Saturated characteristic establishes the mathematical model of executing agency's constraint.Then corresponding optimization object function, packet are designed according to mission requirements
The quadratic form of the input containing executing agency and the quadratic form of system state amount error, time and energy expenditure are executed to consider
Problem.Control problem is finally converted into the inequality constraints item of the equation and state and input-bound in system dynamics equation
Under part, the mathematical problem of the extreme value of object function is sought;With the optimized treatment method rapid solving of real-time iterative, optimization is obtained
It solves and is exported as the controlled quentity controlled variable of system.The block diagram that control system solves control input process is as shown in Figure 2.
Specific implementation step is as follows:
The first step establishes the mathematical model of attitude control system.The kinetics equation of flywheel is attached to spacecraft attitude
In kinetics equation:
Wherein, ω=[ω1,ω2,ω3]TIndicate the attitude angle of spacecraft relative inertness coordinate system under body coordinate system
Velocity vector, ω1,ω2,ω3Respectively spacecraft is about the angular speed on roll axis, yaw axis and the pitch axis in this system
Component;Indicate the derivative of ω relative times;J is the total inertia matrix of spacecraft, is simplified shown as diagonal matrix J=diag (J1,
J2,J3), J1=30kg/m2,J2=50kg/m2,J3=40kg/m2For around the rotary inertia of the principal axis of inertia;S (ω) is skew symmetry square
Battle array, form areτ indicates the output torque of executing agency;Q=[q0,q1,q2,q3]TIndicate space flight
The posture unit quaternion of device,The scalar component in the posture unit quaternion of spacecraft is indicated, and around Euler's axis
The angle of rotation is related, and θ indicates the angle crossed around Euler's shaft rotation,
ex,ey,ezThe rotary shaft on three directions of Euler's axis is represented, and is met Indicate that q leads the time
Number;Ω (ω) is skew symmetric matrix, and form is
The model of counteraction flyback combination is as follows:
Wherein HrwFor the angular momentum of flywheel,For HrwTo the derivative of time, in the flywheel being made of four flywheels
The relationship of its angular momentum and rotating speed is as follows in combination:
Hrw=CJrwN
WhereinMatrix, N=[n are installed for flywheel1,n2,n3,n4]TFor flywheel angular speed to
Amount, n1,n2,n3,n4The angular speed of each flywheel is indicated respectively;JrwIndicate the moment of inertia matrix of flywheel, form Jrw
=JαI4×4, Jα=0.01608kg/m2Indicate the rotary inertia of single flywheel, I4×4For 4 rank unit matrixs.
Since executing agency is to control posture by way of exchanging angular momentum with spacecraft, so the angular motion that system is total
Measure conservation:
H=Hrw+Jω
H is the total angular momentum of system, is a constant when no moment of face interferes.
The model integration of the model of executing agency and spacecraft attitude dynamics is got up to obtain:
By above-mentioned model discretization, the setting sampling interval is that Δ t=0.2s is obtained, and is carved in kth:
ωk+1=J-1ΔtS(ωk)H+J-1JαCΔNk+ωk
Wherein subscript k indicates the value to dependent variable at the kth moment, Δ Nk=Nk-Nk-1, I3×3For 3 rank unit matrixs.
The initial angular velocity of the initial attitude of spacecraft, initial angular velocity and flywheel is respectively qintial=[-
0.9524,-0.3048,0,0]Tωintial=[0,0,0]TAnd Nintial=[0,0,0,0]T, control desired posture and angle speed
It spends to be not qd=[1,0,0,0]Tωd=[0,0,0]T.Therefore the initial angular momentum H=0 that system is total.
Second step, establish posture of spacecraft during attitude maneuver evade constraint and executing agency performance about
Beam:
(1) counteraction flyback maximum output torque constrains:
The torque output of counteraction flyback is to change the angular momentum realization of flywheel, so there is following form to state:
Wherein Tmax=[1 11 1]TNm is maximum output torque vector.It will be obtained after above-mentioned formula discretization:
It is obtained after arrangement:
(2) counteraction flyback maximum angular momentum constrains:
The angular momentum saturation of flywheel is presented as that the rotating speed of flywheel rotor reaches the upper limit, so the full constraint of angular momentum can be used and fly
Angular speed constraint is taken turns to indicate:
-Nmax≤ΔNk+Nk-1≤Nmax
Wherein Nmax=[200 π, 200 π, 200 π, 200 π]TRad/s is maximum angular rate vector.
It is obtained after arrangement:
(3) angle of sight constraint that spacecraft attitude is directed toward:
Consider that the direction of spacecraft will evade certain conical lines-of-sight areas, the following form of posture restraint of design:
Wherein α=[0 0 1]TIndicate that spacecraft is directed toward unit vector under body coordinate system,It indicates
Spacecraft needs the unit vector for evading direction under body coordinate system,Indicate the size for evading the angle of sight in region.
Optimization object function of the third step design based on model predictive controller:
Wherein xk=[qk,ωk]TThe quantity of state of expression system, physical significance are the posture and angular speed of spacecraft;With
Flywheel rotation speed change amount is as input, i.e. uk=Δ Nk;NP=5 be the estimation range of MPC;Q=I7×7And P=I4×4For shape
The weight matrix of state variable and input variable, I7×7For 7 rank unit matrixs, illustrate that optimization aim is more noted if Q is larger relative to P
The stable time is overweighted, illustrates that optimization aim is more concerned with the consumption of energy if larger relative to P if Q.It is organized into compact form:
Wherein
Control problem is converted to the extreme-value problem for seeking object function under constraints, real-time iterative is used in combination by the 4th step
Method solve controlled quentity controlled variable, by restricted model in the system model and second step in the first step and third step in optimization aim
Function Synthesis is got up, and will be solved the problem of control inputs and is converted to following mathematical problem:
xi,k+1=F (xi,k,ui,k)
L(xi,k,ui,k)≤017×1K=0,1 ..., NP-1
WhereinIndicate the optimization problem in moment i,Indicate the feedback of status in current time i system, xi,k
It indicates by current time state xi,0Speculate the system mode at the kth moment;ui,kIndicate the estimated input at the kth moment;F
(xi,k,ui,k) be embodied as the system state equation in the first step:
L(xi,k,ui,k) be second step in system restriction, be embodied as:
Processing procedure using the method for real-time iterative is:Known previous optimization problemSolve ask
OutWithAnd then solving-optimizing problemProcess, tool
Body is divided into two stages:Preparation stage and response phase
Preparation phase procedures are:By xi-1And ui-1The displacement of a sampling instant is done, the last one position is constant, obtainsWithAccording toWithCalculate sensitive square
Battle array Ai,k,Bi,k,Ci,k,Di,kWith error li,k,ri,k, expression is:
By linearisation as above, the optimization problem of a nonlinear restriction is thus converted to the secondary of linear restriction
Planning problem:
Δui,k=ui,k-ui-1,k
Δxi,k+1=Ai,kΔxi,k+Bi,kΔui,k+ri,k
Ci,kΔxi,k+Di,kΔui,k+li,k≤ 0 k=0,1 ..., NP-1
Response phase process is:The feedback of status of acquisition systemBring Solve problems intoObtain Δ xi,kWith Δ ui,k
Pass through formula:
Obtain problemSolution (xi,k,ui,k).Finally again by ui=ui,0Control as system inputs.
The case where control method more than utilization can make spacecraft during attitude maneuver, executing agency's limited capacity
Under, it carries out posture and evades, while the synthesis that can meet energy and time is optimal.The optimization method provided can improve the solution of MPC
Efficiency accelerates the speed of solving-optimizing problem using the thought of thermal starting, and the thought that real-time iterative is divided into two stages can
Control output is made for the feedback of status at a nearest moment, improves the timeliness that system executes control, the flow of optimization
As shown in Figure 2.
Although for illustrative purposes, it has been described that exemplary embodiments of the present invention, those skilled in the art
Member it will be understood that, can be in form and details in the case of not departing from the scope and spirit invented disclosed in appended claims
On the change that carry out various modifications, add and replace etc., and all these changes should all belong to appended claims of the present invention
Protection domain, and each step in the claimed each department of product and method, can be in any combination
Form is combined.Therefore, to disclosed in this invention the description of embodiment be not intended to limit the scope of the invention,
But for describing the present invention.Correspondingly, the scope of the present invention is not limited by embodiment of above, but by claim or
Its equivalent is defined.
Claims (7)
1. a kind of model predictive control method for considering the spacecraft attitude of flywheel dynamic characteristic and evading, which is characterized in that including
Following steps:
(1) it is established comprising execution according to the attitude kinematics of spacecraft and the rotational power of kinetic characteristics and flywheel characteristic
Prediction model of the spacecraft attitude model of mechanism as MPC;
(2) it carries the angle of sight of instrument according to spacecraft and evades vector and establish the mathematical model that posture evades constraint, and root
The mathematical model of executing agency's constraint is established according to the angular momentum saturation and torque saturated characteristic of flywheel;
(3) corresponding optimization object function is designed according to mission requirements, includes the quadratic form and system shape of executing agency's input
The quadratic form of state error executes time and energy expenditure to consider;
(4) control problem is converted into the inequality constraints condition of the equation and state and input-bound in system dynamics equation
Under, the extreme value of object function is sought, with the optimized treatment method rapid solving of real-time iterative, the solution that optimization is obtained is as system
Controlled quentity controlled variable exports.
2. according to the method described in claim 1, it is characterized in that:The kinetics equation of flywheel is attached to spacecraft attitude to move
In mechanical equation, and sliding-model control, establish the prediction model of MPC.
3. according to the method described in claim 2, it is characterized in that:Spacecraft attitude dynamics model in step (1) is as follows
It indicates:
Wherein, ω=[ω1,ω2,ω3]TIndicate spacecraft relative inertness coordinate system under body coordinate system attitude angular velocity to
Amount, ω1,ω2,ω3Respectively spacecraft is about the angular velocity component on roll axis, yaw axis and the pitch axis in this system;
Indicate derivatives of the ω to the time;J is the total inertia matrix of spacecraft, is simplified shown as diagonal matrix J=diag (J1,J2,J3), J1,
J2,J3For around the rotary inertia of the principal axis of inertia;S (ω) is skew symmetric matrix, and form isτ
Indicate the output torque of executing agency;Q=[q0,q1,q2,q3]TIndicate the posture unit quaternion of spacecraft,Table
Show that the scalar component in the posture unit quaternion of spacecraft, θ indicate the angle crossed around Euler's shaft rotation,ex,ey,ezThe rotary shaft on three directions of Euler's axis is represented, and full
Foot Indicate derivatives of the q to the time;Ω (ω) is skew symmetric matrix, and form is as follows:
The model of counteraction flyback combination is as follows:
Wherein HrwFor the angular momentum of flywheel,For the derivative of the angular momentum relative time of flywheel, by four flywheels
The relationship of its angular momentum and rotating speed is as follows in the flywheel of composition:
Hrw=CJrwN
Wherein C is that 3 × 4 flywheels install matrix, N=[n1,n2,n3,n4]TFor the angular velocity vector of flywheel, n1,n2,n3,n4Respectively
Indicate the angular speed of each flywheel;JrwIndicate the moment of inertia matrix of flywheel, form Jrw=JαI4×4, JαIndicate single
The rotary inertia of a flywheel, I4×4For 4 rank unit matrixs;
Since executing agency is to control posture by way of exchanging angular momentum with spacecraft, so the total angular momentum of system is kept
It is permanent:
H=Hrw+Jω
H is the total angular momentum of system, is a constant in no moment of face;
The model integration of the model of executing agency and spacecraft attitude dynamics is got up to obtain:
By above-mentioned model discretization, the setting sampling interval is that Δ t is obtained, and is carved in kth:
ωk+1=J-1ΔtS(ωk)H+J-1JαCΔNk+ωk
Wherein subscript k indicates the value to dependent variable at the kth moment, Δ Nk=Nk-Nk-1, I3×3For three rank unit matrixs.
4. according to the method described in claim 3, it is characterized in that:In step (2):
(a) counteraction flyback maximum output torque constrains:
The torque output of counteraction flyback is to change the angular momentum realization of flywheel, so there is following form to state:
Wherein TmaxFor maximum output torque vector;
It will be obtained after above-mentioned formula discretization:
It is obtained after arrangement:
(b) counteraction flyback maximum angular momentum constrains:
The angular momentum saturation of flywheel is presented as that the rotating speed of flywheel rotor reaches the upper limit, so angular momentum constraint of saturation can use flywheel
Angular speed constraint indicates:
-Nmax≤ΔNk+Nk-1≤Nmax
Wherein NmaxFor flywheel maximum angular rate vector;
It is obtained after arrangement:
(c) angle of sight constraint that spacecraft attitude is directed toward:
Consider that the direction of spacecraft will evade certain conical lines-of-sight areas, the following form of posture restraint of design:
Wherein α indicates that the unit vector that spacecraft is directed toward under body coordinate system, β indicate that spacecraft needs to advise under body coordinate system
The unit vector in direction is kept away, θ indicates the size for evading the angle of sight in region.
5. according to the method described in claim 4, it is characterized in that:Optimization object function V (x in step (3)k,uk) be expressed as:
Wherein xk=[qk,ωk]TThe quantity of state of expression system, physical significance are the posture and angular speed of spacecraft;With flywheel group
Rotation speed change amount is closed as input, i.e. uk=Δ Nk;NPFor the estimation range of MPC;Q and P is the power of state variable and input variable
Weight matrix, illustrate that optimization aim is more concerned with the stable time if larger relative to P if Q, if Q relative to P it is larger illustrate it is excellent
Change the consumption that target is more concerned with energy;
It is organized into compact form:
Wherein
6. according to the method described in claim 5, it is characterized in that:Control problem is converted to and seeks target letter under constraints
Several extreme-value problem, the method that real-time iterative is used in combination solve controlled quentity controlled variable.
7. according to the method described in claim 6, it is characterized in that:The problem of control input will be solved, is converted to following mathematics and asks
Topic:
xi,k+1=F (xi,k,ui,k)
L(xi,k,ui,k)≤017×1K=0,1 ..., NP-1
WhereinIndicate the optimization problem in moment i,Indicate the feedback of status in current time i system, xi,kTable
Show by current time state xi,0The system mode at the kth moment of reckoning;ui,kThe estimated input at expression kth moment, 017×1It indicates
17 × 1 null matrix;F(xi,k,ui,k) it is system state equation, it is embodied as:
L(xi,k,ui,k) it is system restriction, it is embodied as:
Processing procedure using the method for real-time iterative is:Known previous moment optimization problemSolve find out
ComeWithAnd then solving-optimizing problemSpecific point
For two stages:Preparation stage and response phase:
Preparation phase procedures are:By xi-1And ui-1The displacement for doing a sampling instant keeps the last one element constant, obtainsWithAccording toWithCalculate sensitive square
Battle array Ai,k,Bi,k,Ci,k,Di,kWith error li,k,ri,k, expression is:
By linearisation, the optimization problem of a nonlinear restriction is converted to the quadratic programming problem of linear restriction:
Δui,k=ui,k-ui-1,k
Δxi,k+1=Ai,kΔxi,k+Bi,kΔui,k+ri,k
Ci,kΔxi,k+Di,kΔui,k+li,k≤ 0k=0,1 ..., NP-1
The process of response phase is:The feedback of status of acquisition systemBring Solve problems intoObtain Δ xi,kWith Δ ui,kIt is logical
Cross formula:
Obtain problemSolution (xi,k,ui,k), finally again by ui=ui,0Control as system inputs.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810298997.1A CN108536014B (en) | 2018-04-04 | 2018-04-04 | Model prediction control method for spacecraft attitude avoidance considering dynamic characteristics of flywheel |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810298997.1A CN108536014B (en) | 2018-04-04 | 2018-04-04 | Model prediction control method for spacecraft attitude avoidance considering dynamic characteristics of flywheel |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108536014A true CN108536014A (en) | 2018-09-14 |
CN108536014B CN108536014B (en) | 2020-02-21 |
Family
ID=63483072
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810298997.1A Active CN108536014B (en) | 2018-04-04 | 2018-04-04 | Model prediction control method for spacecraft attitude avoidance considering dynamic characteristics of flywheel |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108536014B (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109782787A (en) * | 2019-03-08 | 2019-05-21 | 北京航空航天大学 | A kind of solar light pressure assists the bimodulus MPC control method of lower drive lacking spacecraft attitude |
CN110329548A (en) * | 2019-05-24 | 2019-10-15 | 中国人民解放军63789部队 | In-orbit turn of biasing of spacecraft controls lower fly wheel system reconstructing method |
CN110502025A (en) * | 2019-08-01 | 2019-11-26 | 北京航空航天大学 | A kind of Spacecraft Attitude Control method considering counteraction flyback characteristic and power limited |
CN111240194A (en) * | 2018-11-28 | 2020-06-05 | 罗伯特·博世有限公司 | Model predictive regulation better considering constraints |
CN111897355A (en) * | 2020-08-06 | 2020-11-06 | 中国科学院微小卫星创新研究院 | Satellite attitude maneuver trajectory planning method |
CN113485397A (en) * | 2021-07-08 | 2021-10-08 | 北京航空航天大学 | Spacecraft attitude maneuver path planning method based on polynomial planning |
CN114115308A (en) * | 2021-11-12 | 2022-03-01 | 北京理工大学 | Gravitational wave detection spacecraft attitude autonomous planning method |
CN114714351A (en) * | 2022-04-06 | 2022-07-08 | 上海工程技术大学 | Anti-saturation target tracking control method and control system for mobile mechanical arm |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103092208A (en) * | 2013-01-09 | 2013-05-08 | 哈尔滨工业大学 | Spacecraft high-accuracy speediness attitude maneuver method based on single gimbal control moment gyro (SGCMG) and reaction wheel (RW) |
CN103235597A (en) * | 2013-04-09 | 2013-08-07 | 北京理工大学 | Rapid stable joint control method for attitude maneuver of spacecraft |
CN104866692A (en) * | 2015-06-18 | 2015-08-26 | 北京理工大学 | Aircraft multi-objective optimization method based on self-adaptive agent model |
CN105629732A (en) * | 2016-01-29 | 2016-06-01 | 北京航空航天大学 | Spacecraft attitude output feedback tracking control method giving consideration to limited control |
US20160194095A1 (en) * | 2015-01-07 | 2016-07-07 | Mitsubishi Electric Research Laboratories, Inc. | Model Predictive Control of Spacecraft |
CN105843244A (en) * | 2016-06-02 | 2016-08-10 | 北京航空航天大学 | Output feedback-based flexible spacecraft precise attitude control method |
CN105867395A (en) * | 2015-12-28 | 2016-08-17 | 北京理工大学 | Method for planning attitude maneuvering in presence of constraints of deep space probe based on sliding film control |
CN105867401A (en) * | 2016-04-28 | 2016-08-17 | 北京航空航天大学 | Spacecraft posture fault tolerance control method of single gimbal control moment gyroscope groups |
CN107065902A (en) * | 2017-01-18 | 2017-08-18 | 中南大学 | UAV Attitude fuzzy adaptive predictive control method and system based on nonlinear model |
CN107168357A (en) * | 2017-06-30 | 2017-09-15 | 北京航空航天大学 | It is a kind of to consider posture restraint and the spacecraft attitude maneuver control method of anti-unwinding |
CN107491082A (en) * | 2016-12-31 | 2017-12-19 | 南京航空航天大学 | Spacecraft Attitude Control mixing executing agency optimal control method |
-
2018
- 2018-04-04 CN CN201810298997.1A patent/CN108536014B/en active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103092208A (en) * | 2013-01-09 | 2013-05-08 | 哈尔滨工业大学 | Spacecraft high-accuracy speediness attitude maneuver method based on single gimbal control moment gyro (SGCMG) and reaction wheel (RW) |
CN103235597A (en) * | 2013-04-09 | 2013-08-07 | 北京理工大学 | Rapid stable joint control method for attitude maneuver of spacecraft |
US20160194095A1 (en) * | 2015-01-07 | 2016-07-07 | Mitsubishi Electric Research Laboratories, Inc. | Model Predictive Control of Spacecraft |
CN104866692A (en) * | 2015-06-18 | 2015-08-26 | 北京理工大学 | Aircraft multi-objective optimization method based on self-adaptive agent model |
CN105867395A (en) * | 2015-12-28 | 2016-08-17 | 北京理工大学 | Method for planning attitude maneuvering in presence of constraints of deep space probe based on sliding film control |
CN105629732A (en) * | 2016-01-29 | 2016-06-01 | 北京航空航天大学 | Spacecraft attitude output feedback tracking control method giving consideration to limited control |
CN105867401A (en) * | 2016-04-28 | 2016-08-17 | 北京航空航天大学 | Spacecraft posture fault tolerance control method of single gimbal control moment gyroscope groups |
CN105843244A (en) * | 2016-06-02 | 2016-08-10 | 北京航空航天大学 | Output feedback-based flexible spacecraft precise attitude control method |
CN107491082A (en) * | 2016-12-31 | 2017-12-19 | 南京航空航天大学 | Spacecraft Attitude Control mixing executing agency optimal control method |
CN107065902A (en) * | 2017-01-18 | 2017-08-18 | 中南大学 | UAV Attitude fuzzy adaptive predictive control method and system based on nonlinear model |
CN107168357A (en) * | 2017-06-30 | 2017-09-15 | 北京航空航天大学 | It is a kind of to consider posture restraint and the spacecraft attitude maneuver control method of anti-unwinding |
Non-Patent Citations (1)
Title |
---|
缪远明等: "规避姿态禁区的航天器姿态机动路径规划", 《航天器工程》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111240194A (en) * | 2018-11-28 | 2020-06-05 | 罗伯特·博世有限公司 | Model predictive regulation better considering constraints |
CN111240194B (en) * | 2018-11-28 | 2024-05-24 | 罗伯特·博世有限公司 | Model predictive regulation with better consideration of constraints |
CN109782787A (en) * | 2019-03-08 | 2019-05-21 | 北京航空航天大学 | A kind of solar light pressure assists the bimodulus MPC control method of lower drive lacking spacecraft attitude |
CN110329548A (en) * | 2019-05-24 | 2019-10-15 | 中国人民解放军63789部队 | In-orbit turn of biasing of spacecraft controls lower fly wheel system reconstructing method |
CN110502025A (en) * | 2019-08-01 | 2019-11-26 | 北京航空航天大学 | A kind of Spacecraft Attitude Control method considering counteraction flyback characteristic and power limited |
CN111897355A (en) * | 2020-08-06 | 2020-11-06 | 中国科学院微小卫星创新研究院 | Satellite attitude maneuver trajectory planning method |
CN111897355B (en) * | 2020-08-06 | 2022-09-13 | 中国科学院微小卫星创新研究院 | Satellite attitude maneuver trajectory planning method |
CN113485397B (en) * | 2021-07-08 | 2024-02-02 | 北京航空航天大学 | Spacecraft attitude maneuver path planning method based on polynomial planning |
CN113485397A (en) * | 2021-07-08 | 2021-10-08 | 北京航空航天大学 | Spacecraft attitude maneuver path planning method based on polynomial planning |
CN114115308A (en) * | 2021-11-12 | 2022-03-01 | 北京理工大学 | Gravitational wave detection spacecraft attitude autonomous planning method |
CN114115308B (en) * | 2021-11-12 | 2023-02-17 | 北京理工大学 | Gravitational wave detection spacecraft attitude autonomous planning method |
CN114714351B (en) * | 2022-04-06 | 2023-06-23 | 上海工程技术大学 | Anti-saturation target tracking control method and control system for mobile mechanical arm |
CN114714351A (en) * | 2022-04-06 | 2022-07-08 | 上海工程技术大学 | Anti-saturation target tracking control method and control system for mobile mechanical arm |
Also Published As
Publication number | Publication date |
---|---|
CN108536014B (en) | 2020-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108536014A (en) | A kind of model predictive control method for considering the spacecraft attitude of flywheel dynamic characteristic and evading | |
Diehl et al. | A real-time iteration scheme for nonlinear optimization in optimal feedback control | |
Lai et al. | Stable control strategy for planar three-link underactuated mechanical system | |
Bardhan et al. | Nonlinear differential games-based impact-angle-constrained guidance law | |
CN108319148B (en) | A kind of control-moment gyro Gimbal servo system slow-speed of revolution high-accuracy control method | |
CN109283934B (en) | Spacecraft multi-constraint attitude maneuver optimization method based on rotating path quality | |
Ritz et al. | Quadrocopter performance benchmarking using optimal control | |
CN102343985B (en) | Satellite time optimal posture maneuvering method with reaction flywheel | |
CN104020778B (en) | Flexible Satellite Attitude maneuver autopilot method based on tracking time energy consumption optimal control orbit | |
CN109283841A (en) | Flexibility of joint and armed lever space manipulator control method flexible | |
Lücking et al. | Electrochromic orbit control for smart-dust devices | |
CN112660423B (en) | Method and system for controlling gaze tracking of video satellite on moving target | |
Jiao et al. | Analysis and design the controller for quadrotors based on PID control method | |
CN108663936B (en) | Model does not know spacecraft without unwinding Attitude Tracking finite-time control method | |
CN111338368B (en) | Self-adaptive robust control method for spacecraft rapid maneuver attitude tracking | |
CN108438256A (en) | It is a kind of that satellite attitude control method is stared based on permanent magnetism momentum-exchange ball over the ground | |
Yu | H∞ Tracking-based adaptive fuzzy-neural control for MIMO uncertain robotic systems with time delays | |
CN109164819A (en) | The contragradience adaptive sliding mode large angle attitude control method of rigid body spacecraft | |
CN110377044A (en) | A kind of the finite time height and Attitude tracking control method of unmanned helicopter | |
Williams et al. | Flexible tethered kite with moveable attachment points, part I: Dynamics and control | |
CN106339002B (en) | A kind of control of Solar sail spacecraft three-axis attitude and implementation method | |
CN106096206B (en) | Optimal configuration method for distributed execution mechanism and sensor of flexible spacecraft | |
Pan et al. | Quadrotor control based on self-tuning LQR | |
Guo et al. | PID parameter tuning algorithm of rotor UAV Based on Improved Particle Swarm Optimization | |
CN110134135A (en) | Based on the quadrotor control method for improving MPC-PID |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |