CN109032161B - Small-inertia spacecraft attitude jitter determination method based on fourth-order Runge Kutta method - Google Patents
Small-inertia spacecraft attitude jitter determination method based on fourth-order Runge Kutta method Download PDFInfo
- Publication number
- CN109032161B CN109032161B CN201810870904.8A CN201810870904A CN109032161B CN 109032161 B CN109032161 B CN 109032161B CN 201810870904 A CN201810870904 A CN 201810870904A CN 109032161 B CN109032161 B CN 109032161B
- Authority
- CN
- China
- Prior art keywords
- attitude
- spacecraft
- sensor
- equation
- time
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 239000011159 matrix material Substances 0.000 claims abstract description 45
- 230000004927 fusion Effects 0.000 claims abstract description 17
- 238000012937 correction Methods 0.000 claims abstract description 14
- 238000005259 measurement Methods 0.000 claims abstract description 10
- 238000006073 displacement reaction Methods 0.000 claims description 27
- 239000013598 vector Substances 0.000 claims description 22
- 238000005070 sampling Methods 0.000 claims description 13
- 238000005096 rolling process Methods 0.000 claims description 2
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 101100269157 Caenorhabditis elegans ads-1 gene Proteins 0.000 description 4
- 101100119135 Mus musculus Esrrb gene Proteins 0.000 description 4
- 101100080600 Schizosaccharomyces pombe (strain 972 / ATCC 24843) nse6 gene Proteins 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005094 computer simulation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Navigation (AREA)
Abstract
The invention provides a small inertia spacecraft attitude jitter determination method based on a fourth-order Runge Kutta method, which comprises the following steps of: s1, establishing an attitude kinematics equation of the spacecraft based on the direction cosine matrix; s2, dispersing the attitude kinematics equation by using a fourth-order Runge Kutta method to obtain a recursion equation of a direction cosine matrix; and S3, correcting the recursion equation by using two different types of sensors, wherein the measurement data of one type of sensor is used as a dominant term, and the sensor of the other type is used as a correction term through a data fusion algorithm. The invention has the beneficial effects that: the measurement bandwidth of the attitude jitter is improved, and the attitude jitter of the small-inertia spacecraft under a small angle can be accurately measured in a relatively wide bandwidth range.
Description
Technical Field
The invention relates to a spacecraft, in particular to a small-inertia spacecraft attitude jitter determination method based on a four-order Runge Kutta method.
Background
Small spacecraft are beginning to become more and more powerful for a variety of reasons. Such spacecraft possess less inertia than conventional spacecraft, which makes them highly susceptible to disturbances that produce attitude jitter that is more severe than large spacecraft. This level of jitter is unacceptable to scientific instruments and must be compensated for. The attitude jitter must be measured over a wide bandwidth before compensation.
The existing attitude sensors have the problem of bandwidth limitation, and can only measure within respective bandwidth ranges.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides a small inertia spacecraft attitude jitter determination method based on a four-order Runge Kutta method, which is used for solving the problem of expanding the measurement bandwidth of an attitude sensor.
The invention provides a small inertia spacecraft attitude jitter determination method based on a fourth-order Runge Kutta method, which comprises the following steps of:
s1, establishing an attitude kinematics equation of the spacecraft based on the direction cosine matrix;
s2, dispersing the attitude kinematics equation by using a fourth-order Runge Kutta method to obtain a recursion equation of a direction cosine matrix;
and S3, correcting the recursion equation by using two different types of sensors, wherein the measurement data of one type of sensor is used as a dominant term, and the sensor of the other type is used as a correction term through a data fusion algorithm.
As a further improvement of the present invention, step S1 includes:
establishing an attitude kinematics equation of the spacecraft based on the direction cosine matrix as follows:
wherein C is a direction cosine matrix of the spacecraft,is the cosine of the directionFirst derivative of the matrix with respect to time, ω×Is the angular velocity vector ω ═ ωx ωy ωz]TCross-product matrix of (a).
As a further improvement of the present invention, step S2 includes:
a four-step Runge Kutta method is used to perform discrete approximation on the formula (1):
wherein, K1Is the slope at time n; k2Introducing K by an Euler method for the slope of the midpoint between the n moment and the n +1 moment1To obtain the value of (1); k3Also is the slope of the midpoint between the n time and the n +1 time, and introduces K by an Euler method2To obtain the value of (1); k4The slope at the time n + 1; and T is the sampling time of the sensor, and because two different types of sensors are involved in the data fusion algorithm, if the sampling time of the two different types of sensors is different, the T is taken as the least common multiple of the sampling time of the two different types of sensors.
As a further improvement of the present invention, step S2 includes:
the recursive equation of the direction cosine matrix obtained by the fourth-order Runge Kutta method is as follows:
wherein, Cn+1A direction cosine matrix of the spacecraft at time n +1, CnIs a spacecraft direction cosine matrix at the time of n,is a cross-product matrix of the attitude angular velocity vector at the moment n,is a cross-product matrix of the attitude angular velocity vector at the midpoint of the time n and the time n +1,a cross-product matrix of the attitude angular velocity vector at the moment n + 1;andall the attitude angle displacement measured by the sensor is obtained by difference; since the sampling time of the sensor is T, the attitude angular velocity at the midpoint between n and n +1 timesThe sensor can not be obtained and can only be usedApproximately, so there are:
wherein,is a cross-multiplication matrix of the spacecraft attitude angular displacement vector measured by the sensor at the moment of n +1,is a cross-product matrix of the attitude angular displacement vector of the spacecraft measured by the sensor at the moment n,a cross-multiplication matrix of the spacecraft attitude angular displacement vector measured by the sensor at the n-1 moment;
the formula (4) can be substituted for the formula (3):
as a further improvement of the present invention, step S3 includes:
when the angular displacement sensor and the data measured by the inertial reference unit are used for data fusion, the data measured by the angular displacement sensor is used as a leading term, the data measured by the inertial reference unit is used as a low-frequency correction term, and the following steps are performed:
wherein,are cross multiplication matrixes of angular displacement vectors of the spacecraft measured by the angular displacement sensor at the moments of n-1, n and n +1 respectively,the correction terms of the angular displacement are respectively obtained after data measured by the inertia reference unit at the moments of n-1, n and n +1 pass through the closed-loop controller;
the specific expression is as follows:
whereinADSθx、ADSθy、ADSθzThe angular displacement of a rolling shaft, a pitching shaft and a yawing shaft of the spacecraft is measured by an angular displacement sensor; whilecorθx,corθy,corθzIt is the low frequency correction term produced by the closed loop controller.
By substituting equation (6) for equation (5), a recurrence equation after data fusion can be obtained:
the invention has the beneficial effects that: the measurement bandwidth of the attitude jitter is improved, and the attitude jitter of the small-inertia spacecraft under a small angle can be accurately measured in a relatively wide bandwidth range.
Drawings
Fig. 1 is a schematic block diagram of a data fusion algorithm of a small inertia spacecraft attitude jitter determination method based on a fourth-order longstota method according to the present invention.
FIG. 2 is a simulink simulation block diagram of a data fusion algorithm of a small inertia spacecraft attitude jitter determination method based on a fourth-order Runge Kutta method.
Detailed Description
The invention is further described with reference to the following description and embodiments in conjunction with the accompanying drawings.
As shown in fig. 1 to 2, a method for determining small inertia spacecraft attitude jitter based on a fourth-order longge stoke method includes the following steps:
firstly, establishing an attitude kinematics equation of the spacecraft based on a direction cosine matrix as follows:
wherein C is a direction cosine matrix of the spacecraft,is the first derivative of the directional cosine matrix with respect to time, ω×Is the angular velocity vector ω ═ ωx ωy ωz]TCross-product matrix of (a). And because the attitude sensors have sampling time, carrying out discretization processing on the continuous equation by using a four-order Runge Kutta method. According to the fourth-order longgedusta method, there are:
wherein, K1Is the slope at time n; k2Introducing K by an Euler method for the slope of the midpoint between the n moment and the n +1 moment1To obtain the value of (1); k3Also is the slope of the midpoint between the n time and the n +1 time, and introduces K by an Euler method2To obtain the value of (1); k4The slope at time n + 1. T is the sampling time of the sensor, and because two different types of sensors are involved in the data fusion algorithm provided by the patent, the sampling time of the two sensors may be different, and the T can be taken as the least common multiple of the sampling time of the two sensors. The discrete form of the spacecraft attitude kinematics equation obtained by the fourth-order longge stoke method is:
wherein, Cn+1A direction cosine matrix of the spacecraft at time n +1, CnIs a spacecraft direction cosine matrix at the time of n,is a cross-product matrix of the attitude angular velocity vector at the moment n,is a cross-product matrix of the attitude angular velocity vector at the midpoint time of n time and n +1 time,is a cross-product matrix of the attitude angular velocity vector at the moment n + 1.Andcan be obtained by differentiating the angular displacement of the attitude measured by the sensor. Since attitude sensor sampling time T is constant, the time is the midpoint between n and n +1Attitude angular velocityThe sensors are not directly available, as used hereinApproximately, so there are:
wherein,and the cross multiplication matrixes are respectively expressed by the three-axis attitude angular displacement vectors of the spacecraft measured by the attitude sensor at the moments of n-1, n and n + 1.
By substituting equation (13) into equation (12), a recursive form of the spacecraft attitude kinematics equation can be obtained:
in order to expand the measurement bandwidth of the attitude sensor, two sensors with different bandwidths can be used for data fusion, data of one sensor is used as a leading term, data obtained by the other sensor through a closed-loop controller is used as a correction term, a schematic block diagram of the attitude sensor is shown in fig. 1, a sensor 1 is used as the leading term, and a sensor 2 is used as the correction term.
For example, when the Angular Displacement Sensor (ADS) is selected to be fused with the Inertial Reference Unit (IRU) data, the ADS data is used as a dominant term, and the IRU data is used as a correction term. In the course of updating directional cosine matrix, order
Wherein,are cross multiplication matrixes of spacecraft angular displacement vectors measured by the ADS sensor at the moments of n-1, n and n +1 respectively,the correction terms of the angular displacement are respectively obtained after data measured by the IRU sensor at the time of n-1, n and n +1 passes through the closed-loop controller. The specific expression is as follows:
by substituting equation (15) into equation (14), a recursive algorithm of the attitude kinematics equation after data fusion can be obtained:
when the attitude angle of the spacecraft is calculated by the direction cosine matrix obtained after data fusion, because the problem of attitude jitter is concerned, the attitude jitter is very small, and the method can be suitable for small-angle approximation of the attitude angle. The attitude angle of the spacecraft can be directly extracted from the directional cosine matrix:
θ=f(C)=[θx θy θz]T (19)
wherein,
fig. 2 shows a block diagram of a computer simulation program constructed by simulink, in which the M code of the fcn module:
function[x,y,z,c]=fcn(ads0,ads1,ads2,cor0,cor1,cor2,c0)
adsx0=[0 -ads0(3) ads0(2);ads0(3) 0 -ads0(1);-ads0(2) ads0(1) 0];
adsx1=[0 -ads1(3) ads1(2);ads1(3) 0 -ads1(1);-ads1(2) ads1(1) 0];
adsx2=[0 -ads2(3) ads2(2);ads2(3) 0 -ads2(1);-ads2(2) ads2(1) 0];
corx0=[0 -cor0(3) cor0(2);cor0(3) 0 -cor0(1);-cor0(2) cor0(1) 0];
corx1=[0 -cor1(3) cor1(2);cor1(3) 0 -cor1(1);-cor1(2) cor1(1) 0];
corx2=[0 -cor2(3) cor2(2);cor2(3) 0 -cor2(1);-cor2(2) cor2(1) 0];
cx0=[c0(1) c0(2) c0(3);c0(4) c0(5) c0(6);c0(7) c0(8) c0(9)];
fused0=adsx0+corx0;
fused1=adsx1+corx1;
fused2=adsx2+corx2;
err1=fused1-fused0;
err2=fused2-fused1;
cx1=(eye(3)-err1+((err1^2)/3)-((err1^3)/12)+(err2*err1/6)-(err2*(err1^2)/12)+(err2*(err1^3)/24))*cx0;
x=0.5*(cx1(2,3)-cx1(3,2));
y=0.5*(cx1(3,1)-cx1(1,3));
z=0.5*(cx1(1,2)-cx1(2,1));
c=[cx1(1,1);cx1(1,2);cx1(1,3);cx1(2,1);cx1(2,2);cx1(2,3);cx1(3,1);cx1(3,2);cx1(3,3)]。
the invention provides a small inertia spacecraft attitude jitter determination method based on a fourth-order Runge Kutta method, aiming at the problem that a single sensor has measurement bandwidth limitation, a data fusion algorithm of two sensors is designed based on the fourth-order Runge Kutta method. The invention aims to expand the measurement bandwidth of the attitude jitter, thereby laying a foundation for inhibiting the attitude jitter. The method adopts a spacecraft attitude kinematics equation expressed by a directional cosine matrix, uses a four-order Runge Kutta method to carry out discretization, uses the measurement data of one sensor as a leading item and the data of the other sensor as a correction item, processes the discretized attitude kinematics equation, and obtains an attitude kinematics recurrence equation of a data fusion algorithm. The ADS sensor and the IRU sensor are taken as examples, and specific expression forms of the recurrence equation are given. And finally, verifying the effectiveness of the data fusion algorithm through simulink computer simulation design.
The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and it is not intended that the invention be limited to these specific details. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.
Claims (1)
1. A small inertia spacecraft attitude jitter determination method based on a fourth-order Runge Kutta method is characterized by comprising the following steps:
s1, establishing an attitude kinematics equation of the spacecraft based on the direction cosine matrix;
s2, dispersing the attitude kinematics equation by using a fourth-order Runge Kutta method to obtain a recursion equation of a direction cosine matrix;
s3, correcting the recursion equation by using two sensors of different types, wherein the measurement data of one sensor is used as a leading term, and the sensor of the other sensor is used as a correction term;
wherein,
step S1 includes:
establishing an attitude kinematics equation of the spacecraft based on the direction cosine matrix as follows:
wherein C is a direction cosine matrix of the spacecraft,is the first derivative of the directional cosine matrix with respect to time, ω×Is the angular velocity vector ω ═ ωx ωy ωz]TA cross-product matrix of;
step S2 includes:
a four-step Runge Kutta method is used to perform discrete approximation on the formula (1):
wherein, K1Is the slope at time n; k2Introducing K by an Euler method for the slope of the midpoint between the n moment and the n +1 moment1To obtain the value of (1); k3Also is the slope of the midpoint between the n time and the n +1 time, and introduces K by an Euler method2To obtain the value of (1); k4The slope at the time n + 1; t is sampling time of the sensors, and because two different types of sensors are involved in the data fusion algorithm, if the sampling time of the two different types of sensors is different, the T is taken as the least common multiple of the sampling time of the two different types of sensors;
step S2 includes:
the recursive equation of the direction cosine matrix obtained by the fourth-order Runge Kutta method is as follows:
wherein, Cn+1A direction cosine matrix of the spacecraft at time n +1, CnIs a spacecraft direction cosine matrix at the time of n,is a cross-product matrix of the attitude angular velocity vector at the moment n,is a cross-product matrix of the attitude angular velocity vector at the midpoint of the time n and the time n +1,a cross-product matrix of the attitude angular velocity vector at the moment n + 1;andall the attitude angle displacement measured by the sensor is obtained by difference; since the sampling time of the sensor is T, the attitude angular velocity at the midpoint between n and n +1 timesThe sensor can not be obtained and can only be usedApproximately, so there are:
wherein,is a cross-multiplication matrix of the spacecraft attitude angular displacement vector measured by the sensor at the moment of n +1,fork for spacecraft attitude angular displacement vector measured by sensor at n momentThe multiplication matrix is used to multiply the data,a cross-multiplication matrix of the spacecraft attitude angular displacement vector measured by the sensor at the n-1 moment;
the formula (4) can be substituted for the formula (3):
step S3 includes:
when the angular displacement sensor and the data measured by the inertial reference unit are used for data fusion, the data measured by the angular displacement sensor is used as a leading term, the data measured by the inertial reference unit is used as a low-frequency correction term, and the following steps are performed:
wherein,are cross multiplication matrixes of angular displacement vectors of the spacecraft measured by the angular displacement sensor at the moments of n-1, n and n +1 respectively,the correction terms of the angular displacement are respectively obtained after data measured by the inertia reference unit at the moments of n-1, n and n +1 pass through the closed-loop controller;
the specific expression is as follows:
whereinADSθx、ADSθy、ADSθzThe angular displacement of a rolling shaft, a pitching shaft and a yawing shaft of the spacecraft is measured by an angular displacement sensor; whilecorθx,corθy,corθzThen it is the low frequency correction term produced by the closed loop controller;
by substituting equation (6) for equation (5), a recurrence equation after data fusion can be obtained:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810870904.8A CN109032161B (en) | 2018-08-02 | 2018-08-02 | Small-inertia spacecraft attitude jitter determination method based on fourth-order Runge Kutta method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810870904.8A CN109032161B (en) | 2018-08-02 | 2018-08-02 | Small-inertia spacecraft attitude jitter determination method based on fourth-order Runge Kutta method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109032161A CN109032161A (en) | 2018-12-18 |
CN109032161B true CN109032161B (en) | 2021-05-07 |
Family
ID=64647656
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810870904.8A Expired - Fee Related CN109032161B (en) | 2018-08-02 | 2018-08-02 | Small-inertia spacecraft attitude jitter determination method based on fourth-order Runge Kutta method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109032161B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112550766B (en) * | 2020-11-27 | 2022-04-08 | 上海航天控制技术研究所 | Method for improving satellite attitude control precision in thruster dead zone |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103512584A (en) * | 2012-06-26 | 2014-01-15 | 北京赛佰特科技有限公司 | Navigation attitude information output method, device and strapdown navigation attitude reference system |
CN105629732A (en) * | 2016-01-29 | 2016-06-01 | 北京航空航天大学 | Spacecraft attitude output feedback tracking control method giving consideration to limited control |
CN106324643A (en) * | 2016-10-19 | 2017-01-11 | 山东科技大学 | Method for estimating airspeed of unmanned aerial vehicle and detecting failures of airspeed head |
CN106643737A (en) * | 2017-02-07 | 2017-05-10 | 大连大学 | Four-rotor aircraft attitude calculation method in wind power interference environments |
CN107831775A (en) * | 2017-11-14 | 2018-03-23 | 哈尔滨工业大学深圳研究生院 | The attitude control method without angular velocity measurement based on Spacecraft |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104850746B (en) * | 2015-05-22 | 2017-08-18 | 国家电网公司 | A kind of equivalent salt density Forecasting Methodology based on fourth order Runge-Kutta and simulated annealing |
CN107655494A (en) * | 2017-09-15 | 2018-02-02 | 哈尔滨工程大学 | Inertial navigation system coarse alignment method under the conditions of a kind of swaying base |
-
2018
- 2018-08-02 CN CN201810870904.8A patent/CN109032161B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103512584A (en) * | 2012-06-26 | 2014-01-15 | 北京赛佰特科技有限公司 | Navigation attitude information output method, device and strapdown navigation attitude reference system |
CN105629732A (en) * | 2016-01-29 | 2016-06-01 | 北京航空航天大学 | Spacecraft attitude output feedback tracking control method giving consideration to limited control |
CN106324643A (en) * | 2016-10-19 | 2017-01-11 | 山东科技大学 | Method for estimating airspeed of unmanned aerial vehicle and detecting failures of airspeed head |
CN106643737A (en) * | 2017-02-07 | 2017-05-10 | 大连大学 | Four-rotor aircraft attitude calculation method in wind power interference environments |
CN107831775A (en) * | 2017-11-14 | 2018-03-23 | 哈尔滨工业大学深圳研究生院 | The attitude control method without angular velocity measurement based on Spacecraft |
Non-Patent Citations (5)
Title |
---|
Attitude Calculating for a New Strap-down Inertial Navigation System Based on Angular Accelerometers;Meiling Wang等;《International Conference on Information and Automation》;20161231;第1467-1472页 * |
Detection of RF signal real-time phase based on spectralhole-burning material;Xiu-Rong Ma等;《Optik》;20180129;第68-76页 * |
四旋翼无人飞行器控制***设计与实现;王福超;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;中国学术期刊(光盘版)电子杂志社;20140615(第06期);全文 * |
基于微惯性组合的旋翼飞行器姿态检测及控制算法研究;高洪涛;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;中国学术期刊(光盘版)电子杂志社;20170215(第02期);第10-19、35-40页 * |
基于虚拟仪器的捷联惯性导航算法测试平台;冯为荣;《中国优秀硕士论文电子期刊网》;20101115(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109032161A (en) | 2018-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108828955B (en) | Precise track tracking control method based on finite time extended state observer | |
CN103116275B (en) | Based on the gyroscope Robust Neural Network Control system and method that sliding formwork compensates | |
CN105159084B (en) | A kind of manipulator nerve network control system and method with interference observer | |
CN108646556B (en) | Input saturation spacecraft is without unwinding Attitude tracking control method | |
CN104811588B (en) | A kind of boat-carrying based on gyroscope is surely as control method | |
US20100138180A1 (en) | Moving Body Posture Angle Detecting Apparatus | |
CN109625335B (en) | Sun capturing method based on angular velocity estimation information and sun sensor | |
TW201411096A (en) | Method and apparatus for data fusion of a three axis magnetometer and three axis accelerometer | |
CN107992084B (en) | Unmanned aerial vehicle robust attitude control method and device independent of angular velocity feedback | |
CN106772351B (en) | Kalman filter method based on the memory of limited step | |
CN110926460A (en) | Uwb positioning abnormal value processing method based on IMU | |
CN106813679B (en) | Method and device for estimating attitude of moving object | |
CN109164819B (en) | Backstepping self-adaptive sliding mode large-angle attitude maneuver control method of rigid-body spacecraft | |
CN110986928A (en) | Real-time drift correction method for triaxial gyroscope of photoelectric pod | |
CN111189474A (en) | Autonomous calibration method of MARG sensor based on MEMS | |
CN109032161B (en) | Small-inertia spacecraft attitude jitter determination method based on fourth-order Runge Kutta method | |
CN112859600A (en) | Mechanical system finite time control method based on extended state observer | |
Zhang et al. | Integral barrier Lyapunov function-based three-dimensional low-order integrated guidance and control design with seeker's field-of-view constraint | |
Parikh et al. | Comparison of attitude estimation algorithms with IMU under external acceleration | |
CN110345942B (en) | Interpolation three subsampled cone error compensation algorithm based on angular rate input | |
CN110375773B (en) | Attitude initialization method for MEMS inertial navigation system | |
CN108871312B (en) | Combined attitude determination method for gravity gradiometer and star sensor | |
Oliva-Palomo et al. | Robust global observer position-yaw control based on ellipsoid method for quadrotors | |
Faria | Sensor fusion and rotational motion reconstruction via nonlinear state-observers | |
CN105180946A (en) | Wideband measurement-based satellite high-precision attitude determination method and system thereof |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210507 |
|
CF01 | Termination of patent right due to non-payment of annual fee |