CN102384746A - Chebyshev output method for space motion state of rigid body - Google Patents
Chebyshev output method for space motion state of rigid body Download PDFInfo
- Publication number
- CN102384746A CN102384746A CN2011102805995A CN201110280599A CN102384746A CN 102384746 A CN102384746 A CN 102384746A CN 2011102805995 A CN2011102805995 A CN 2011102805995A CN 201110280599 A CN201110280599 A CN 201110280599A CN 102384746 A CN102384746 A CN 102384746A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- chebyshev
- phi
- rigid body
- output
- 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
Abstract
The invention discloses a Chebyshev output method for the space motion state of a rigid body. The method comprises the steps of: defining a ternary number so that three speed components of the set of axes of the body and the ternary number form a system of linear differential equations; and performing approximation description on rolling, pitching and yawing angular speeds p, q and r by using Shifted Chebyshev orthogonal polynomials, wherein the state transfer matrix of the system can be solved in the manner of any order holder, thus obtaining the expression of the discrete state equation of the motion of the rigid body; and the problem of singularity of the attitude equation is avoided; therefore, the main motion state of the rigid body is obtained. The Chebyshev output method is characterized in that the ternary number is introduced so that the state transfer matrix is in the form of portioned upper triangles and thereby can be solved through deflation; therefore, the complexity of calculation is greatly simplified; and the Chebyshev output method is convenient for engineering use.
Description
Technical field
The present invention relates to the spatial movement rigid model, particularly the big maneuvering flight state output of aircraft problem.
Background technology
Axis is that the rigid motion differential equation is a fundamental equation of describing spatial movements such as aircraft, torpedo, spacecraft.Usually, in data processing etc. was used, the state variable of axon system mainly comprised the X of 3 speed components, three Eulerian angle and earth axes
E, Y
E, Z
EDeng, because Z
EBe defined as vertical ground and point to ground ball center, so Z
EActual flying height for bearing; X
E, Y
EUsually main dependence GPS, GNSS, the Big Dipper etc. directly provide; Eulerian angle are represented the rigid space athletic posture, and the differential equation of portrayal rigid body attitude is a core wherein, is that pitching, lift-over and crab angle are described with three Eulerian angle usually.When the angle of pitch of rigid body was ± 90 °, roll angle and crab angle can't definite values, and it is excessive that error is found the solution in the zone of closing on this singular point simultaneously, caused intolerable error on the engineering and can not use; For fear of this problem, people at first adopt the method for restriction angle of pitch span, and this makes equation degenerate, attitude work entirely, thereby be difficult to be widely used in engineering practice.Along with the research to the aircraft extreme flight, people have adopted direction cosine method, equivalent gyration vector method, hypercomplex number method etc. to calculate the rigid motion attitude again in succession.
Direction cosine method has been avoided " unusual " phenomenon of Eulerian angle describing methods, and calculating attitude matrix with direction cosine method does not have the equation degenerate problem, attitude work entirely; But need find the solution 9 differential equations; Calculated amount is bigger, and real-time is relatively poor, can't satisfy the engineering practice requirement.Equivalence gyration vector method such as list appearance recursion, Shuangzi appearance gyration vector, three increment gyration vectors and four increment rotating vector methods and various correction algorithms on this basis and recursive algorithm etc.When studying rotating vector in the document, all be based on the algorithm that rate gyro is output as angle increment.Yet in actual engineering, the output of some gyros is angle rate signals, like optical fibre gyro, dynamic tuned gyroscope etc.When rate gyro was output as angle rate signal, the Algorithm Error of rotating vector method obviously increased.The hypercomplex number method is that the function of 4 Eulerian angle of definition calculates the boat appearance; Can effectively remedy the singularity of Eulerian angle describing method; As long as separate 4 differential equation of first order formula groups; Analogy has tangible minimizing to cosine attitude matrix differential equation calculated amount, can satisfy in the engineering practice requirement to real-time.Its The common calculation methods has the card of finishing approximatioss, second order, fourth-order Runge-Kutta method and three rank Taylor expansion methods etc.Finishing card approximatioss essence is list appearance algorithm, can not compensate by exchange error what limited rotation caused, and the algorithm drift under high current intelligence in the attitude algorithm can be very serious.When adopting fourth-order Runge-Kutta method to find the solution the hypercomplex number differential equation,, the trigonometric function value can occur to exceed ± 1 phenomenon, disperse thereby cause calculating along with the continuous accumulation of integral error; The Taylor expansion method also is restricted because of the deficiency of computational accuracy.When rigid body is big when motor-driven, angular speed causes the error of said method bigger more greatly; Moreover, the error of attitude estimation usually can cause the error of 4 components of speed, highly output sharply to increase.
Summary of the invention
In order to overcome the existing big problem of rigid motion model output error, the present invention provides a kind of Chebyshev's output intent of rigid space motion state, and this method is through definition ternary number; Make that axis is that three speed components and ternary number constitute the linear differential equation group; And adopting Shifted Chebyshev (change Chebyshev) orthogonal polynomial to lift-over, pitching, yaw rate p, q, r carry out close approximation to be described; Can be according to the state-transition matrix of the mode solving system of any rank retainer; And then obtain the expression formula of rigid motion discrete state equations, avoid attitude equation singular problem, thereby obtained the main motion state of rigid body.
The present invention solves the technical scheme that its technical matters adopts, a kind of Chebyshev's output intent of rigid space motion state, and its characteristic may further comprise the steps:
1, axis is that three speed components are output as:
Wherein: u, it is x that v, w are respectively along the rigid body axis, y, the speed component of z axle, n
x, n
y, n
zBe respectively along x, y, the overload of z axle, g is an acceleration of gravity, s
1, s
2, s
3Be the ternary number of definition, and
P, q, r are respectively lift-over, pitching, yaw rate, and T is the sampling period; Parameter-definition is identical in full;
Be the recursive form of Chebyshev (Chebyshev) orthogonal polynomial, b=NT, lift-over, pitching, yaw rate p, q, the expansion of r is respectively
p(t)=[p
0?p
1…p
n-1?p
n][ξ
0(t)ξ
1(t)…ξ
n-1(t)ξ
n(t)]
T
q(t)=[q
0?q
1…q
n-1?q
n][ξ
0(t)ξ
1(t)…ξ
n-1(t)ξ
n(t)]
T
r(t)=[r
0?r
1…r
n-1?r
n][ξ
0(t)ξ
1(t)…ξ
n-1(t)ξ
n(t)]
T
Work as p, q, when the high-order term n of the expansion of r is odd number, m=4,6 ..., n+1, m=5 when high-order term n is even number, 7 ..., n+1, H
i(i=1,2 ..., n) be H row vector accordingly,
2, highly be output as:
Wherein: h is a height;
3, attitude angle is output as:
The invention has the beneficial effects as follows: make that through introducing the ternary number state-transition matrix is a triangular form on the piecemeal, can find the solution state-transition matrix by depression of order, simplified computation complexity greatly, be convenient to engineering and use.
Below in conjunction with embodiment the present invention is elaborated.
Embodiment
1, axis is that three speed components are output as:
Wherein: u, it is x that v, w are respectively along the rigid body axis, y, the speed component of z axle, n
x, n
y, n
zBe respectively along x, y, the overload of z axle, g is an acceleration of gravity, s
1, s
2, s
3Be the ternary number of definition, and
P, q, r are respectively lift-over, pitching, yaw rate, and T is the sampling period; Parameter-definition is identical in full;
p(t)=[p
0?p
1…p
n-1?p
n][ξ
0(t)ξ
1(t)…ξ
n-1(t)ξ
n(t)]
T
q(t)=[q
0?q
1…q
n-1?q
n][ξ
0(t)ξ
1(t)…ξ
n-1(t)ξ
n(t)]
T
r(t)=[r
0?r
1…r
n-1?r
n][ξ
0(t)ξ
1(t)…ξ
n-1(t)ξ
n(t)]
T
Work as p, q, when the high-order term n of the expansion of r is odd number, m=4,6 ..., n+1, m=5 when high-order term n is even number, 7 ..., n+1, H
i(i=1,2 ..., n) be H row vector accordingly,
2, highly be output as:
Wherein: h is a height;
3, attitude angle is output as:.
Claims (1)
1. Chebyshev's output intent of a rigid space motion state, its characteristic may further comprise the steps:
Axis is that three speed components are output as:
Wherein: u, it is x that v, w are respectively along the rigid body axis, y, the speed component of z axle, n
x, n
y, n
zBe respectively along x, y, the overload of z axle, g is an acceleration of gravity, s
1, s
2, s
3Be the ternary number of definition, and
P, q, r are respectively lift-over, pitching, yaw rate, and T is the sampling period; Parameter-definition is identical in full;
Be the recursive form of Chebyshev (Chebyshev) orthogonal polynomial, b=NT, lift-over, pitching, yaw rate p, q, the expansion of r is respectively
p(t)=[p
0?p
1…p
n-1?p
n][ξ
0(t)ξ
1(t)…ξ
n-1(t)ξ
n(t)]
T
q(t)=[q
0?q
1…q
n-1?q
n][ξ
0(t)ξ
1(t)…ξ
n-1(t)ξ
n(t)]
T
r(t)=[r
0?r
1…r
n-1?r
n][ξ
0(t)ξ
1(t)…ξ
n-1(t)ξ
n(t)]
T
Work as p, q, when the high-order term n of the expansion of r is odd number, m=4,6 ..., n+1, m=5 when high-order term n is even number, 7 ..., n+1, H
i(i=1,2 ..., n) be H row vector accordingly,
Highly be output as:
Wherein: h is a height;
Attitude angle is output as:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110280599.5A CN102384746B (en) | 2011-09-20 | 2011-09-20 | Chebyshev output method for space motion state of rigid body |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110280599.5A CN102384746B (en) | 2011-09-20 | 2011-09-20 | Chebyshev output method for space motion state of rigid body |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102384746A true CN102384746A (en) | 2012-03-21 |
CN102384746B CN102384746B (en) | 2014-05-07 |
Family
ID=45824389
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110280599.5A Active CN102384746B (en) | 2011-09-20 | 2011-09-20 | Chebyshev output method for space motion state of rigid body |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102384746B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109724597A (en) * | 2018-12-19 | 2019-05-07 | 上海交通大学 | A kind of inertial navigation resolving method and system based on function iterated integral |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1422380A (en) * | 2000-02-03 | 2003-06-04 | 独立技术有限责任公司 | Attitude estimation in tiltable body using modified quaternion data representation |
US20070213889A1 (en) * | 2004-02-27 | 2007-09-13 | Instituto Nacional De Tecnica Aeroespacial "Esteba Terradas" | Sensor Fusion System and Method for Estimating Position, Speed and Orientation of a Vehicle, in Particular an Aircraft |
CN102072729A (en) * | 2009-11-23 | 2011-05-25 | 北京信息科技大学 | Measurement device and measurement method for measuring posture of rotary aircraft |
-
2011
- 2011-09-20 CN CN201110280599.5A patent/CN102384746B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1422380A (en) * | 2000-02-03 | 2003-06-04 | 独立技术有限责任公司 | Attitude estimation in tiltable body using modified quaternion data representation |
US20070213889A1 (en) * | 2004-02-27 | 2007-09-13 | Instituto Nacional De Tecnica Aeroespacial "Esteba Terradas" | Sensor Fusion System and Method for Estimating Position, Speed and Orientation of a Vehicle, in Particular an Aircraft |
CN102072729A (en) * | 2009-11-23 | 2011-05-25 | 北京信息科技大学 | Measurement device and measurement method for measuring posture of rotary aircraft |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109724597A (en) * | 2018-12-19 | 2019-05-07 | 上海交通大学 | A kind of inertial navigation resolving method and system based on function iterated integral |
Also Published As
Publication number | Publication date |
---|---|
CN102384746B (en) | 2014-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105222780B (en) | A kind of ellipsoid set-membership filtering method approached based on Stirling interpolation polynomial | |
CN102589553B (en) | Switching method for building aircraft motion model | |
CN102436437B (en) | Quaternion Fourier approximate output method in extreme flight of aircraft based on angular speed | |
CN102445202B (en) | Laguerre output method for rigid body space motion state | |
CN102384747A (en) | Hartley output method of rigid body space motion states | |
CN102384746B (en) | Chebyshev output method for space motion state of rigid body | |
CN102359789B (en) | Arbitrary order output method for rigid body space motion state | |
CN102346729B (en) | Legendre output method for space motion space of rigid body | |
CN102323990B (en) | Method for modeling pneumatic model for rigid body space motion | |
CN102323992B (en) | Polynomial type output method for spatial motion state of rigid body | |
CN102508818A (en) | Arbitrary-step orthogonal series output method of space motion state of rigid body | |
CN102508819B (en) | Angular-speed-based quaternion Legendre approximate output method during extreme flying of aircraft | |
CN102445203B (en) | Emmett output method for rigid body space motion state | |
CN102495829B (en) | Quaternion Walsh approximate output method based on angular velocities for aircraft during extreme flight | |
CN102495830B (en) | Quaternion Hartley approximate output method based on angular velocities for aircraft during extreme flight | |
CN102495831B (en) | Quaternion Hermitian approximate output method based on angular velocities for aircraft during extreme flight | |
CN102323991B (en) | Walsh output method for spatial motion state of rigid body | |
CN102359790B (en) | Fourier outputting method for spatial movement state of rigid body | |
CN102508821B (en) | State output method for space motion of rigid body | |
CN102506866B (en) | Angle speed-based Chebyshev approximate output method of quaternion numbers in ultimate flight of aircraft | |
CN102346727B (en) | Body axis model for rigid space motion | |
CN102506865B (en) | Four-ary number polynomial approximate output method during extreme aerobat flight based on angular velocity | |
CN102506864B (en) | Method for approximately outputting quaternion numbers with arbitrary step size in orthogonal series during extreme flight of aircraft | |
CN102494688B (en) | Quaternion Laguerre approximate output method based on angular speed used during extreme flight of flying vehicle | |
CN102375803A (en) | Air stream axle system model for rigid space motion |
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 |