CN104570734B - A kind of all-coefficient adaptive control method based on single order characteristic model - Google Patents
A kind of all-coefficient adaptive control method based on single order characteristic model Download PDFInfo
- Publication number
- CN104570734B CN104570734B CN201410789593.4A CN201410789593A CN104570734B CN 104570734 B CN104570734 B CN 104570734B CN 201410789593 A CN201410789593 A CN 201410789593A CN 104570734 B CN104570734 B CN 104570734B
- Authority
- CN
- China
- Prior art keywords
- theta
- current period
- alpha
- parameter
- feature model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 24
- 239000013598 vector Substances 0.000 claims abstract description 24
- 230000009466 transformation Effects 0.000 claims description 15
- 230000010354 integration Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000009191 jumping Effects 0.000 description 2
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Landscapes
- Feedback Control In General (AREA)
Abstract
The invention discloses a kind of all-coefficient adaptive control method based on single order characteristic model, according to formula <maths num=" 0001 " >
</maths> calculates the preliminary parameters estimate vector θ of current period
1(k); According to the preliminary parameters estimate vector θ of current period
1(k) and a upper cycle estimated parameter vector θ (k) by what estimate that parameter vector θ (k-1) calculates current period; Parameter is estimated according to current period
with
carry out the controlled quentity controlled variable u (k) that linear Feedback Control obtains current period.All-coefficient adaptive control method based on single order characteristic model of the present invention, parameter is few, and algorithm is simple, and greatly and easily determine, the dynamic gain tracking velocity calculated by estimates of parameters is fast, adaptable for parameter area.
Description
Technical Field
The invention relates to the field of intelligent and self-adaptive control, in particular to a full-coefficient self-adaptive control method based on a first-order feature model.
Background
The existing full-coefficient adaptive control method based on the feature model is based on the feature model with the second order or more, for example, in 1998, "manned spacecraft full-coefficient adaptive reentry lift control" published in aerospace science (vol.19, No.1, jan, 1998, p8-12), the paper proposes a full-coefficient adaptive control method based on the drop point prediction and adopting the second-order feature model. The second-order characteristic model is adopted, the upper and lower limits are required to be carried out on each model parameter, and the upper and lower limits are also required to be carried out on the sum of the model parameters, so that the stability of the closed-loop system can be ensured. For a feature model with a second order or more, the upper and lower limit values of the parameter are related to the ratio of the control period to the time constant, and because the controlled object is unknown and the time constant is unknown, the upper and lower limit values of the parameter can only be determined empirically. The method is adopted for controlling the jumping reentry lift force of certain types of spaceflight. In the process of carrying out deep research on the control of the jumping reentry lift force of a certain model, the fact that model parameters sometimes touch or stay at the boundary due to large dynamic change of a controlled object is found, and the parameter tracking and self-adaptive capacity is lost, so that the further improvement of the reentry guidance precision is limited.
Disclosure of Invention
The invention aims to solve the technical problem of providing a full-coefficient self-adaptive control method based on a first-order characteristic model, the adopted first-order characteristic model has few parameters, only limits the parameter range, does not have the limitation of the sum of the parameters, has simple algorithm and large parameter range, and the dynamic gain tracking speed calculated by the parameter estimation value is high, the adaptive capacity is strong, and the first-order characteristic model is constant and stable, and the control effect on a reentry vehicle is superior to that of the full-coefficient self-adaptive control method of the second-order characteristic model.
The technical scheme of the invention is as follows:
a full-coefficient self-adaptive control method based on a first-order feature model is a first-order variable coefficient difference equation and is represented as follows:
y(k+1)=α(k)y(k)+β0(k)u(k)
wherein y (k) is the output transformation quantity of the current period of the controlled object, y (k +1) is the output transformation quantity of the next period of the controlled object, u (k) is the control quantity of the current period, α (k), β0(k) Time-varying parameters of a current period of the first-order feature model;
the full-coefficient self-adaptive control method based on the first-order feature model in each period comprises the following steps:
according to the formula Calculating a preliminary parameter estimation vector theta for the current cycle1(k) Wherein λ is1Is a gain factor, λ2To remove the zero protection factor, 0<λ1<1,λ2>0; phi (k-1) is the regression vector of the previous cycle, phiT(k-1)=[y(k-1),u(k-1)]T(ii) a Theta (k-1) is the estimated parameter vector of the previous cycle,
estimating a vector theta from the preliminary parameters of the current period1(k) And calculating the estimated parameter vector theta (k) of the current period by using the estimated parameter vector theta (k-1) of the previous period, wherein theta (k) is Fθ1θ1(k)+(1-Fθ1) Theta (k-1) wherein Fθ1Is a filter coefficient, 0<Fθ1Less than or equal to 1; determining the estimated time-varying parameter of the current period according to the estimated parameter vector theta (k) of the current periodWherein,
judging the estimated time-varying parameter of the current periodWhether or not to satisfy αIs the minimum value of the estimated time-varying parameter, if satisfied, the estimated time-varying parameter of the current period is retainedIf it is not Then order If it is not Then order The estimated time-varying parameter of the current period is judged by the same methodWhether or not to satisfy
Estimated parameters according to the current cycleAndlinear feedback control is performed to obtain a controlled variable u (k) in the current cycle, u (k) ═ L α (k) y (k)/(β)0(k) + λ), L, λ are gain adjustment parameters, 0, respectively<L≤1,λ≥0。
The method is used for the jump-reentry lift control of the reentry vehicle.
Minimum of estimated time-varying parametersαThe value range is as follows: 0.01<α≤0.4。
αTake 0.3.
And the control quantity u (k) is subjected to output conversion to obtain a control quantity for feedback, and the control quantity for feedback is applied to a controlled object.
The control quantity u (k) is subjected to integration and output transformation to obtain a lift force control quantity for feedback, and the control quantity for feedback is applied to a controlled object; the formula for the integral is: u. of1(k)=u1(k-1) + u (k), wherein u1(k) Is the integrated control quantity of the current period, u1And (k-1) is the control quantity after integration of the previous period.
The output conversion amount y (k) is determined by the following method: and calculating the error between the output of the controlled object and the reference output, and then performing input transformation on the error to obtain an output transformation quantity y (k).
Compared with the prior art, the invention has the advantages that:
(1) the invention provides a first-order variable coefficient feature model for the first time, wherein the first-order feature model has few parameters, and the first-order feature model has only 2 parameters compared with 3 parameters of a second-order feature model.
(2) The dynamic gain calculated by the parameter estimation value has high tracking speed and strong adaptability. Compared with a second-order characteristic model and a first-order characteristic model, the dynamic gain tracking speed calculated by the parameter estimation value is high and the adaptability is strong because the parameter sum limitation does not exist.
(3) The parameter range is large and easy to determine. Compared with the parameter range of the second-order characteristic model, the parameter range of the first-order characteristic model is large, and the dynamic gain variation range calculated by the method is also large.
(4) The method is already used in the actual flight of the jump reentry of the reentry aircraft, and the high-precision control effect is achieved.
Drawings
FIG. 1 is a block diagram of an incremental first-order feature model and its full-coefficient adaptive control method;
FIG. 2 is a block diagram of a full-scale first-order feature model and its full-coefficient adaptive control method;
FIG. 3 is a graph of height (m) versus time(s) for an example;
FIG. 4 is a plot of velocity (m/s) versus time(s) of the reflector of an embodiment relative to the earth's surface;
FIG. 5 is an axial apparent velocity increment VXP (m/s) versus time(s) plot for an embodiment;
FIG. 6 is an Nx-time(s) curve of an axial overload of an embodiment.
Detailed Description
The conventional full-coefficient adaptive control based on the feature model is based on the feature model of more than or equal to the second order, but the conventional full-coefficient adaptive control based on the feature model is only applied to the second order feature model. The full-coefficient self-adaptive control based on the first-order feature model provided by the invention provides a lower-order first-order feature model only with the lower order, and the effectiveness of the method is shown for the mathematical simulation and the actual flight application of certain type of jump reentry. The full-coefficient adaptive control based on the first-order feature model is not available in the field of intelligent adaptive control in the past, is innovative and developed in the field of intelligent adaptive control, and has great significance for the development of the intelligent adaptive control. Since the first-order feature model is applied to the reentry aircraft in many design steps compared with the first-order feature model applied to other fields, the specific implementation steps are described by taking the first-order feature model as an example, and other applications can be referred to or reduced for use.
(1) Input and output of definite controlled object and input and output of controller
The input of the controlled object is recorded as uallThe output is r; the required reference output is r0The error of the output from the reference output is y1=r-r0(ii) a For y1The error result after input transformation is y, and the output u is output after the controller based on the first-order characteristic model.
If the incremental controller is selected according to the characteristics of the controlled object(see figure 1), u is the increment of the control quantity, and the control quantity u is obtained by integration1U obtained by output conversion2Is a feedback control amount;
if a full-scale controller (see FIG. 2) is selected according to the characteristics of the controlled object, u is the controlled variable, and u is obtained by output conversion2Is a feedback control amount;
feedback control quantity u2And a reference input u0The sum is the input u of the controlled objectallIt is the total control quantity of the controlled object.
Here, the reference input, the reference output, the input transformation and the output transformation are determined off-line in the task preparation process. For an incremental controller, the integrator is invariant and does not fit within the adaptive controller. Therefore, the relation between y and u is the real-time algorithm of the adaptive controller.
The invention proposes for the first time to represent the relationship between y and u with a first order feature model,
the first-order feature model is a first-order coefficient-of-variation difference equation expressed as follows:
y(k+1)=α(k)y(k)+β0(k)u(k)
wherein y (k) is the output quantity of the current period of the controlled object, u (k) is the control quantity of the current period, y (k +1) is the output quantity of the next period of the controlled object, y (k) and y (k +1) are obtained by measuring values or transforming the measuring values, α (k) and β0(k) Time-varying parameters of the current period of the first-order feature model;
the process of the full-coefficient self-adaptive control method based on the first-order feature model in each period is as follows:
defining a regression vector phiT(k-1)=[y(k-1),u(k-1)]TPhi (k-1) is the regression vector of the previous cycle;
defining estimated parameter vectorsθ (k) is the estimated parameter vector of the period;
calculating a preliminary parameter estimation vector theta of the current period according to the following formula1(k),
Wherein λ1Is a gain factor, λ2To remove the zero protection factor, 0<λ1<1,λ2>0。
θ(k)=Fθ1θ1(k)+(1-Fθ1)θ(k-1),
Wherein Fθ1Is a filter coefficient, 0<Fθ1≤1。
Determining the estimated time-varying parameter of the current period according to the estimated parameter vector theta (k) of the current period Wherein,
judging the estimated time-varying parameter of the current periodWhether or not to satisfy αIs the minimum value of the estimated time-varying parameter, if satisfied, the estimated time-varying parameter of the current period is retainedIf it is not Then order If it is not Then order Judging the estimated time-varying parameter of the current periodWhether or not to satisfy If it is not Then order If it is notThen orderWherein 0.01<αLess than or equal to 0.4, and is determined according to simulation research on an actual object.
Estimated parameters according to the current cycleAndlinear feedback control is performed to obtain a controlled variable u (k) in the current cycle, u (k) ═ L α (k) y (k)/(β)0(k) + λ), L, λ are gain adjustment parameters, 0, respectively<L≤1,λ≥0。
In the first period, the initial values of the parameters are α (k) ═ α (k-1) ═ 0.9, and β (k) ═ β (k-1) ═ 0.1; u (k-1) ═ 0.
Examples
Aiming at the jump type reentry lift control of the reentry aircraft, a full-coefficient adaptive controller based on a first-order characteristic model is designed according to a block diagram shown in figure 1, and the reference input is a nominal roll angle u2Is a lift feedback control quantity; input transformationOutput conversion toA takes the value 36 and D (tcf) is a known dynamic gain function.
Multiple iterations through design and mathematical simulation to determineα=0.3、λ1=1、λ2=4、Fθ10.1, 0.012 and 1. The actual flight has good effect, the expected aim is achieved, the parachute opening point error is 680m, and the control precision is at the international leading level.
Fig. 3-6 show the altitude curve, the speed curve, the axial apparent speed increment curve, and the axial overload curve, respectively. 3-6, the returner realizes perfect jump reentry, the primary reentry section attenuates the speed from 10.674km/s to 7.2km/s within the first cosmic speed, the secondary reentry section reduces the speed from 7.2km/s to 0.143km/s of the parachute opening point, and the speed attenuation and overload process of the primary reentry section, the leap-out section and the secondary reentry section of the whole process are clearly distinguished.
Claims (7)
1. A full-coefficient self-adaptive control method based on a first-order feature model is characterized in that the first-order feature model is a first-order variable coefficient difference equation and is expressed as follows:
y(k+1)=α(k)y(k)+β0(k)u(k)
wherein y (k) is the output transformation quantity of the current period of the controlled object, y (k +1) is the output transformation quantity of the next period of the controlled object, u (k) is the control quantity of the current period, α (k), β0(k) Time-varying parameters of the current period of the first-order feature model;
the full-coefficient self-adaptive control method based on the first-order feature model in each period comprises the following steps:
according to the formula Calculating a preliminary parameter estimation vector theta for the current cycle1(k) Wherein λ is1Is a gain factor, λ2To remove the zero protection factor, 0<λ1<1,λ2>0; phi (k-1) is the regression vector of the previous cycle, phiT(k-1)=[y(k-1),u(k-1)]T(ii) a Theta (k-1) is the estimated parameter vector of the previous cycle, estimated time-varying parameters for the previous cycle;
estimating a vector theta from the preliminary parameters of the current period1(k) And calculating the estimated parameter vector theta (k) of the current period by using the estimated parameter vector theta (k-1) of the previous period, wherein theta (k) is Fθ1θ1(k)+(1-Fθ1) Theta (k-1) wherein Fθ1Is a filter coefficient, 0<Fθ1Less than or equal to 1; determining the estimated time-varying parameter of the current period according to the estimated parameter vector theta (k) of the current period Wherein,
judging the estimated time-varying parameter of the current control periodWhether or not to satisfy αIs the minimum value of the estimated time-varying parameter, if satisfied, the estimated time-varying parameter of the current period is retainedIf it is not Then order If it is not Then order Judging the estimated time-varying parameter of the current periodWhether or not to satisfy If it is not Then order If it is not Then order
Estimated parameters according to the current control periodAndlinear feedback control is performed to obtain a controlled variable u (k) in the current cycle, u (k) ═ L α (k) y (k)/(β)0(k) + λ), L, λ are gain adjustment parameters, 0, respectively<L≤1,λ≥0。
2. The full-coefficient adaptive control method based on the first-order feature model is used for the jump-type reentry lift control of the reentry vehicle according to claim 1.
3. The full-coefficient adaptive control method based on the first-order feature model as claimed in claim 1, wherein the minimum value of the estimated time-varying parameterαThe value range is as follows: 0.01 < (R) >α≤0.4。
4. The full-coefficient adaptive control method based on the first-order feature model according to claim 1,αtake 0.3.
5. The full-coefficient adaptive control method based on the first-order feature model of claim 1, wherein the control quantity u (k) is subjected to output transformation to obtain a control quantity for feedback, and the control quantity for feedback is applied to a controlled object.
6. The full-coefficient adaptive control method based on the first-order feature model is characterized in that the control quantity u (k) is subjected to integration and output transformation to obtain a lift control quantity for feedback, and the lift control quantity for feedback is applied to a controlled object; the formula for the integral is: u. of1(k)=u1(k-1) + u (k), wherein u1(k) Is the integrated control quantity of the current period, u1And (k-1) is the control quantity after integration of the previous period.
7. The full-coefficient adaptive control method based on the first-order feature model of claim 1, wherein the determination method of the output transformation quantity y (k) is as follows: and calculating the error between the output of the controlled object and the reference output, and then performing input transformation on the error to obtain an output transformation quantity y (k).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410789593.4A CN104570734B (en) | 2014-12-18 | 2014-12-18 | A kind of all-coefficient adaptive control method based on single order characteristic model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410789593.4A CN104570734B (en) | 2014-12-18 | 2014-12-18 | A kind of all-coefficient adaptive control method based on single order characteristic model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104570734A CN104570734A (en) | 2015-04-29 |
CN104570734B true CN104570734B (en) | 2016-01-27 |
Family
ID=53087104
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410789593.4A Active CN104570734B (en) | 2014-12-18 | 2014-12-18 | A kind of all-coefficient adaptive control method based on single order characteristic model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104570734B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101364085A (en) * | 2008-09-12 | 2009-02-11 | 北京控制工程研究所 | Non-linear differential calculus golden cut adaptive control method |
CN101393458A (en) * | 2008-10-30 | 2009-03-25 | 北京控制工程研究所 | Longitudinal control method for high altitude lifting of aerospaceplane |
CN102033491A (en) * | 2010-09-29 | 2011-04-27 | 北京控制工程研究所 | Method for controlling flexible satellite based on feature model |
CN201917426U (en) * | 2010-12-27 | 2011-08-03 | 朱力 | Performance testing platform for miniature propulsion system |
CN102176160A (en) * | 2010-12-31 | 2011-09-07 | 北京控制工程研究所 | Hypersonic aircraft self-adaptive learning control method based on deviator characteristic model |
CN103267517A (en) * | 2013-05-24 | 2013-08-28 | 北京控制工程研究所 | Man-controlled rendezvous and docking measurement method based on television cameras and target |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005063073A1 (en) * | 2005-12-29 | 2007-07-12 | Airbus Deutschland Gmbh | Method of sizing and fabricating stiffened structural components, use of structural condition sensors and aircraft |
-
2014
- 2014-12-18 CN CN201410789593.4A patent/CN104570734B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101364085A (en) * | 2008-09-12 | 2009-02-11 | 北京控制工程研究所 | Non-linear differential calculus golden cut adaptive control method |
CN101393458A (en) * | 2008-10-30 | 2009-03-25 | 北京控制工程研究所 | Longitudinal control method for high altitude lifting of aerospaceplane |
CN102033491A (en) * | 2010-09-29 | 2011-04-27 | 北京控制工程研究所 | Method for controlling flexible satellite based on feature model |
CN201917426U (en) * | 2010-12-27 | 2011-08-03 | 朱力 | Performance testing platform for miniature propulsion system |
CN102176160A (en) * | 2010-12-31 | 2011-09-07 | 北京控制工程研究所 | Hypersonic aircraft self-adaptive learning control method based on deviator characteristic model |
CN103267517A (en) * | 2013-05-24 | 2013-08-28 | 北京控制工程研究所 | Man-controlled rendezvous and docking measurement method based on television cameras and target |
Non-Patent Citations (1)
Title |
---|
胡军.载人飞船全系数自适应再入升力控制.《宇航学报》.1998,(第1期), * |
Also Published As
Publication number | Publication date |
---|---|
CN104570734A (en) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110928182B (en) | Robust self-adaptive repetitive control method of hydraulic servo system based on state estimation | |
WO2021240905A1 (en) | Nonlinear optimization method for stochastic predictive control | |
Veronesi et al. | Simultaneous closed-loop automatic tuning method for cascade controllers | |
US11932262B2 (en) | Stochastic nonlinear predictive controller and method based on uncertainty propagation by Gaussian-assumed density filters | |
CN103927436A (en) | Self-adaptive high-order volume Kalman filtering method | |
CN111562794B (en) | Spacecraft attitude control method for actuator fault and input quantization | |
Abel et al. | Prescribed-time safety design for a chain of integrators | |
Tong et al. | Adaptive fuzzy observer backstepping control for a class of uncertain nonlinear systems with unknown time-delay | |
CN105093932B (en) | A kind of method of robustness determining LPV Gain-scheduling control device | |
CN104570734B (en) | A kind of all-coefficient adaptive control method based on single order characteristic model | |
KR102605904B1 (en) | Method and apparatus for attitude tracking control for rigid spacecraft with arbitrary disturbance | |
Zhu et al. | A singular perturbation approach for time-domain assessment of phase margin | |
CN107421543B (en) | Implicit function measurement model filtering method based on state dimension expansion | |
Patete et al. | Stability of self‐tuning control based on Lyapunov function | |
CN111856941B (en) | Adaptive terminal dynamic sliding mode control method based on active disturbance rejection | |
Hirose et al. | Simplified internal model control for time delay processes | |
Cooper et al. | L1 adaptive control for aircraft air management system pressure-regulating bleed valve | |
CN105955024B (en) | Ship supercharged steam generator main steam pressure forecast Control Algorithm | |
De La Cruz et al. | Switching adaptive control of mobile robots | |
Sawyer | Gain-scheduled control of a quadcopter UAV | |
KR102017522B1 (en) | Control Apparatus for Nonlinear System Based on Nonstrict-Feedback | |
Wang et al. | Filtering controller in sliding mode: From the estimation to control | |
Yang et al. | Practical non-linear model predictive control pnmpc: Algorithm implementations | |
Ulsoy | Improving stability margins via time delay control | |
Mustafa et al. | Reduced order observer design with DMPC and LQR for system with backlash nonlinearity |
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 |