CN108995829B - Platform on-orbit calibration method - Google Patents

Platform on-orbit calibration method Download PDF

Info

Publication number
CN108995829B
CN108995829B CN201810716502.2A CN201810716502A CN108995829B CN 108995829 B CN108995829 B CN 108995829B CN 201810716502 A CN201810716502 A CN 201810716502A CN 108995829 B CN108995829 B CN 108995829B
Authority
CN
China
Prior art keywords
load
axis
real
mass center
force
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
Application number
CN201810716502.2A
Other languages
Chinese (zh)
Other versions
CN108995829A (en
Inventor
张科备
关新
汤亮
王有懿
齐田雨
郝仁剑
朱琦
邢林峰
张聪
张强
于国庆
周元子
李勇
马官营
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Control Engineering
Original Assignee
Beijing Institute of Control Engineering
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Control Engineering filed Critical Beijing Institute of Control Engineering
Priority to CN201810716502.2A priority Critical patent/CN108995829B/en
Publication of CN108995829A publication Critical patent/CN108995829A/en
Application granted granted Critical
Publication of CN108995829B publication Critical patent/CN108995829B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/244Spacecraft control systems
    • B64G1/245Attitude control algorithms for spacecraft attitude control

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)

Abstract

A platform on-orbit calibration method, in particular to a six-degree-of-freedom Gough-Stewart platform on-orbit calibration method, reduces the load triaxial attitude coupling coefficient of an actuator of an active pointing hyperstatic platform by compensating the actuator force coefficient of the active pointing hyperstatic platform. The method comprises the following steps: according to the three-axis attitude and the three-axis main inertia of the load whole subjected to normalization processing, determining a resultant moment normalization result actually suffered by the three axes of the mass center of the load whole; and according to the translation displacement of the load mass center subjected to normalization processing, normalizing the three-axis actual acting force of the load mass center. And performing iterative computation to give an optimal solution of the force coefficient of the actuator according to the translational displacement of the load mass center subjected to normalization processing and the three-axis actual acting force of the load mass center subjected to normalization processing. Compensating driving current of the active directional hyperstatic actuator according to the identified actuator force coefficient, realizing load triaxial attitude decoupling of the active directional hyperstatic platform, and reducing triaxial attitude coupling coefficient.

Description

Platform on-orbit calibration method
Technical Field
The invention belongs to the field of spacecraft attitude control, and relates to a platform on-orbit calibration method.
Background
At present, the optical load of the spacecraft points at high precision and the kinematics parameters of which the control can not be accurate are high in stability. And because the counter electromotive force coefficients of the motors in the actuators are different, the output forces of the motors are different under the same input condition. This allows significant coupling of the three-axis attitude of the load. Therefore, the force coefficient of an actuator for decoupling the load attitude is required to be calibrated, and the load triaxial attitude coupling coefficient is reduced.
The single actuator force coefficient calibration method has the following defects:
1. in the existing hyperstatic platform load high-precision control, the force coefficient calibration of a single actuator is usually only carried out, and the centralized calibration of the force coefficients of a plurality of actuators cannot be realized;
2. the existing calibration technology needs to obtain the load mass center attitude and the mass center translational displacement simultaneously in the calibration process of a plurality of actuators. When the translational displacement of the load mass center cannot be obtained, the existing calibration method is difficult to realize the accurate calibration of the force coefficients of a plurality of actuators.
3. The active directional hyperstatic platform formed by a plurality of actuators has a platform configuration error which directly influences the three-axis attitude error of the load. Single actuator force coefficient calibration method in prior art does not calibrate load centroid Jacobian matrix JpThe load attitude error caused by configuration deviation, namely the Jacobian matrix J of the centroid of the load can not be solvedpThe difference between the actual value and the theoretical design causes the load attitude deviation problem.
Disclosure of Invention
The technical problem solved by the invention is as follows: the method overcomes the defects of the prior art, can greatly reduce the load triaxial attitude coupling coefficient, and provides a technical basis for high-precision pointing, high-stability control and rapid and stable control of the optical load of the spacecraft.
The technical solution of the invention is as follows:
a platform on-orbit calibration method comprises the following steps:
1) taking all the loads carried by the platform as a whole, and solving the mass center of the whole asDetermining the resultant moment normalization result tau actually suffered by the three shafts of the load equivalent mass center according to the three-shaft attitude and the three-shaft main inertia of the load equivalent mass centerreal,norm(ii) a According to the translational displacement r of the equivalent center of mass of the loadpDetermining the load equivalent mass center three-axis actual acting force normalization result Freal,norm
2) Iteratively correcting the force normalization result according to the resultant moment normalization result and the resultant moment expected value taur=[τ1τ2τ3]Finding the optimal solution, tau, with the force normalization result approaching 01,τ2,τ3Respectively obtaining expected values of moments of a rolling shaft, a pitching shaft and a yawing shaft of the load equivalent mass center;
3) and determining the force coefficient of the actuator according to the optimal solution of the acting force normalization result approaching to 0 and the resultant moment normalization result, and finishing calibration.
Normalizing the result according to the resultant moment and the resultant moment expected value taur=[τ1τ2τ3]The method for finding the optimal solution of the acting force normalization result approaching 0 specifically comprises the following steps:
21) determining an index function J according to a resultant moment normalization result and an acting force normalization resultreal;
22) Judgment index function JrealWhether the absolute value is less than the critical value, if so, the index function JrealDetermining the normalized result of the acting force as the optimal solution; if not, entering step 23); the critical value is less than 10-3
23) Iteratively correcting the action force normalization result, and repeating the step 22) until an optimal solution of the action force normalization result is obtained.
The index function JrealThe method specifically comprises the following steps:
Figure BDA0001717692310000021
Freal,norm(n)=[FrealxFrealyFrealz],
τreal,norm=[τrealxτrealyτrealz],
τr=[τ1τ2τ3],
Freal,norm(0)=Freal,norm
wherein, the lower foot marks realx, realy and realz respectively represent the components on the rolling axis, the pitching axis and the yawing axis of the equivalent center of mass of the load, and tau1,τ2,τ3Roll axis, pitch axis and yaw axis torque expectation values, n being 0,1, 2.
The method for iteratively correcting the force normalization result specifically comprises the following steps:
Figure BDA0001717692310000031
wherein k isrealAnd updating the coefficient for the acting force normalization result, wherein the value range is (-1, 1), and t is time.
The platform is a six-degree-of-freedom Gough-Stewart platform;
the step 3) determines that the actuator force coefficient α is [ α ]1… α6]The specific method comprises the following steps:
Figure BDA0001717692310000032
wherein, βijElement of the ith row and jth column of the β matrix, β ═ J (J)p T)-1,JpIs a Jacobian matrix of equivalent center of mass of the load, when βi,3+kWhen equal to 0, αi=1,i=1,2,...,6,j=1,2,...,6,k=1,2,3。
Step 1) determining the resultant moment normalization result tau actually suffered by the three axes of the equivalent center of mass of the loadreal,normThe method specifically comprises the following steps:
τreal,norm=k2·θpnorm·Ipnorm
Figure BDA0001717692310000033
Figure BDA0001717692310000034
wherein, thetap=[θpxpypz]For the load equivalent centroid three-axis coupling attitude measurement, Ipx,Ipy,IpzRespectively are main inertia measured values of a rolling axis, a pitching axis and a yawing axis of the equivalent center of mass of the load, k2The load equivalent mass center rotation stiffness coefficient is obtained according to the stiffness of the platform actuator.
Step 1) determining an acting force normalization result Freal,normThe method specifically comprises the following steps:
Freal,norm=(k1·rpnorm)·m,
Figure BDA0001717692310000041
wherein r ispx,rpy,rpzDesigned values of translational displacement r of load equivalent centroid roll axis, pitch axis and yaw axispx,rpy,rpzThe value ranges of (1) and (0.1) are all (-0.1); m is the load mass, k1The rigidity coefficient of the load equivalent mass center three-axis translational displacement is obtained according to the rigidity of the platform actuator.
Compared with the prior art, the invention has the beneficial effects that:
1) the calibration model of the 6 actuator force coefficients α is established, the optimal solution of the load centroid acting force is solved by establishing an index function, the concentrated calibration of the plurality of actuator force coefficients is realized, the inconsistent unification of the plurality of actuator force coefficients is realized, and the identification calibration of the plurality of actuator force coefficients is realized.
2) According to the method, the translation displacement of the load mass center is not required to be measured, the optimal solution of the three-axis load mass center acting force is searched in an iteration mode by setting the initial value of the translation displacement of the load mass center, and the error between the actual value and the expected value of the three-axis load mass center acting force is minimized. On the basis, the optimal solution of the force coefficient of the actuator is obtained, and the calibration precision of the force coefficient of the actuator is improved.
3) The three-axis attitude of the load comprises a Jacobian matrix J due to the centroid of the loadpAnd load attitude deviation caused by configuration deviation. According to the invention, the Jacobian matrix J of the centroid of the load can be compensated by measuring the three-axis attitude of the loadpAnd load attitude errors caused by configuration deviations.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
The invention adopts the flow shown in FIG. 1 to complete the calibration method of the force coefficient of the load attitude decoupling actuator of the active pointing hyperstatic platform.
1) Taking all loads carried by the platform as a whole, solving the mass center of the whole as a load equivalent mass center, and determining a resultant moment normalization result tau actually received by three shafts of the load equivalent mass center according to the three-shaft attitude and the three-shaft main inertia of the load equivalent mass centerreal,norm(ii) a According to the translational displacement r of the equivalent center of mass of the loadpDetermining the load equivalent mass center three-axis actual acting force normalization result Freal,norm
2) Iteratively correcting the force normalization result according to the resultant moment normalization result and the resultant moment expected value taur=[τ1τ2τ3]Finding the optimal solution, tau, with the force normalization result approaching 01,τ2,τ3Respectively obtaining expected values of moments of a rolling shaft, a pitching shaft and a yawing shaft of the load equivalent mass center;
3) and determining the force coefficient of the actuator according to the optimal solution of the acting force normalization result approaching to 0 and the resultant moment normalization result, and finishing calibration.
Normalizing the result according to the resultant moment and the resultant moment expected value taur=[τ1τ2τ3]The method for finding the optimal solution of the acting force normalization result approaching 0 specifically comprises the following steps:
21) normalizing the result and the acting force according to the resultant momentDetermining an index function J from the resultsreal
22) Judgment index function JrealWhether the absolute value is less than the critical value, if so, the index function JrealDetermining the normalized result of the acting force as the optimal solution; if not, entering step 23);
23) iteratively correcting the action force normalization result, and repeating the step 22) until an optimal solution of the action force normalization result is obtained.
The specific method comprises the following steps:
(1) establishing a kinematics model of the translational displacement of the actuator and the load mass center generalized displacement, which specifically comprises the following steps:
δL=Jp[rpθp]T
where δ L is the displacement of the actuator, rpIs the translational displacement of the center of mass of the load, thetap=[θpxpypz]Three-axis coupled attitude measurements of the load. J. the design is a squarepIs a load centroid Jacobian matrix.
(2) The establishment of actuator output force, actuator force coefficient α and load mass center generalized force dynamics model specifically comprises the following steps:
Figure BDA0001717692310000051
wherein f isr=[fr1,…,fr6,]TThe expected output force of the six actuators; frealIs the three-axis actual acting force of the load centroid, taurealThe resultant moment actually suffered by the three axes of the load mass center.
(3) Expected output force and actual output force model of actuator
Constructing a force coefficient calibration model of the active pointing hyperstatic actuator by utilizing a kinematic equation of the displacement of the active pointing hyperstatic actuator and the load mass center, the output force of the actuator, the resultant force borne by the load mass center and a moment equation; wherein, the load coupling attitude measurement is used as the first input quantity of the calibration model, and the initial value F of the load translation control force is usedrealFor calibrating second input quantity of model(ii) a The force coefficient of the actuator is the output quantity of the calibration model; the actuator kinematic relationship reflects the relationship between the force, moment of the load and the actuator force, expressed as:
Figure BDA0001717692310000061
wherein, FrThe expected force of the load mass center is obtained, and F is obtained because the hyperstatic platform only carries out load attitude controlr≡0;τr=[τ1τ2τ3],τ1,τ2,τ3Roll axis, pitch axis and yaw axis moment expectation values of the load centroid, respectively. f. ofr=[fr1,…,fr6,]TThe expected output force of the six actuators; j. the design is a squarepIs a load centroid Jacobian matrix.
The model of the expected value and the actual output force of the actuator is represented as:
freal=αfr
wherein f isrealFor the actual output force of each actuator, α ═ α1… α6]Is the ratio coefficient of the actual output to the desired output force of the actuator.
(4) Load attitude testing and coupling coefficient calculation
Carrying out load three-axis attitude measurement of the active pointing hyperstatic platform to obtain a load three-axis attitude coupling measurement value thetap=[θpxpypz]And calculating the three-axis attitude coupling coefficient of the load as
γij=θji
In the above formula, γijRepresenting the attitude coupling coefficient of the ith axis attitude maneuver to the jth axis; thetaiRepresenting the measured attitude of the ith shaft during maneuvering; thetajRepresenting the attitude component of the j-th axis coupling. The results of the three-axis attitude coupling coefficient test of the load are shown in the following table.
Figure BDA0001717692310000062
Figure BDA0001717692310000071
(5) Load attitude and load three-axis inertia normalization
The three-axis attitude and inertia normalization of the load is carried out according to the following formula, and the result is
Figure BDA0001717692310000072
Wherein, Ip=[Ipx,Ipy,Ipz]Is the three-axis principal inertia of the load; the results of the normalized test are shown in the table below.
Figure BDA0001717692310000073
Three-axis attitude normalization array theta through load equivalent centroidpnormLoad-equivalent triaxial main inertia normalization array IpnormDetermining the resultant moment normalization result tau actually suffered by the three axes of the equivalent center of mass of the loadreal,norm
τreal,norm=k2·θpnorm·Ipnorm
Figure BDA0001717692310000074
Figure BDA0001717692310000075
Wherein, thetap=[θpxpypz]For the load equivalent centroid three-axis coupling attitude measurement, Ipx,Ipy,IpzRespectively are main inertia measured values of a rolling axis, a pitching axis and a yawing axis of the equivalent center of mass of the load, k2The load equivalent mass center rotation stiffness coefficient is obtained according to the stiffness of the platform actuator.
Determining a force normalization result Freal,normThe method specifically comprises the following steps:
Freal,norm=(k1·rpnorm)·m,
Figure BDA0001717692310000076
wherein r ispx,rpy,rpzDesigned values of translational displacement r of load equivalent centroid roll axis, pitch axis and yaw axispx,rpy,rpzThe value ranges of (1) and (0.1) are all (-0.1); m is the load mass, k1The rigidity coefficient of the load equivalent mass center three-axis translational displacement is obtained according to the rigidity of the platform actuator.
(6) Taking the maneuvering of a certain axis of the load as an example, a calibration model is constructed as follows:
Figure BDA0001717692310000081
wherein, β ═ J (J)p T)-1,βijElement of i row and j column of β matrix when βi4When equal to 0, αiWhen β is equal to 1i4When not equal to 0, then have
Figure BDA0001717692310000082
Freal,norm=[FrealxFrealyFrealz],
τreal,norm=[τrealxτrealyτrealz],
τr=[τ1τ2τ3]And k is 1,2 and 3 are respectively expected values of the moment of the roll axis, the moment of the pitch axis and the moment of the yaw axis of the load centroid.
(7) Taking x-axis maneuvering as an example, the resultant force of the loads in the x and y directions is 0, i.e., Frealx=0,FrealyConstructing a load z-axis resultant force F as 0realzThe model for the actuator force coefficient α is:
λi=(βi4τrealxi5τrealyi6τrealz)/βi4
γi=βi3i4
the actuator force coefficient α is then denoted as FrealzThe function is, in order,
Figure BDA0001717692310000091
(8) calculating an index function Jreal(ii) a Let initial FrealAnd (5) calculating the Euclidean norm of the composite force and moment output by the motor of the actuator and the expected force and moment of the load.
Figure BDA0001717692310000092
Freal,norm(n)=[FrealxFrealyFrealz],
τreal,norm=[τrealxτrealyτrealz],
τr=[τ1τ2τ3],
Freal,norm(0)=Freal,norm
Wherein, the lower foot marks realx, realy and realz respectively represent the components on the rolling axis, the pitching axis and the yawing axis of the equivalent center of mass of the load, and tau1,τ2,τ3Roll axis, pitch axis and yaw axis torque expectation values, n being 0,1, 2.
(9) For the index function JrealDerivation and judgment of the index function JrealThe trend of the derivative of (c). If the absolute value of the derivative meets the threshold, i.e. the indicator function JrealAbsolute value of derivative < 10-3Judging that the actuator force coefficient identification result meets the requirement and quitting; if not, entering the next step;
(10) resultant force F for updating three-axis actual acting force of load mass centerreal,norm
Figure BDA0001717692310000093
Freal,norm(0)=Freal,norm
Wherein k isrealFor load centroid triaxial force normalization result Freal,normUpdating the coefficient, wherein the optimal value range is (-1, 1); k is a positive integer and t is time.
(11) Taking x-axis maneuver as an example for illustration, F is calculated by iterationrealz-0.1590. The power coefficient identification results are shown in the following table:
item(s) α1 α2 α3 α4 α5 α6
Value of 1.0385 1.0391 0.7823 0.7823 1.0461 1.0450
(12) And (5) solidifying the actuator force compensation coefficient obtained in the step (11) in a load DSP controller, and performing actuator compensation control to realize load triaxial attitude decoupling control. The actuator force compensation coefficient is calculated as:
kfi=1/αi
the results of the actuator force compensation coefficient test are shown in the following table.
Item(s) kf1 kf2 kf3 kf4 kf5 kf6
Value of 0.9629 0.9624 1.2782 1.2800 0.9559 0.9569
(13) And (5) carrying out a load triaxial attitude test on the actuator force compensation coefficient obtained in the step (12), wherein the test result of the load triaxial attitude coupling coefficient is as shown in the following table.
Item Before calibration After calibration
Coupling coefficient of load X to Y axis 8.22% 0.50%
Coupling coefficient of load X to Z axis 0.48% 0.28%
Coupling coefficient of load Y to X axis 7.65% 0.29%
Coupling coefficient of load Y to Z axis 0.85% 0.84%
Z-to-X axis coupling coefficient of load 1.0% 0.73%
Load Z-to-Y axis coupling coefficient 0.5% 0.35%
Those skilled in the art will appreciate that those matters not described in detail in the present specification are well known in the art.

Claims (3)

1. A platform on-orbit calibration method is characterized by comprising the following steps:
1) taking all loads carried by the platform as a whole, solving the mass center of the whole as a load equivalent mass center, and determining a resultant moment normalization result tau actually received by three shafts of the load equivalent mass center according to the three-shaft attitude and the three-shaft main inertia of the load equivalent mass centerreal,norm(ii) a According to the translational displacement r of the equivalent center of mass of the loadpDetermining the load equivalent mass center three-axis actual acting force normalization result Freal,norm(ii) a The three axes are a rolling axis, a pitching axis and a yawing axis of the load equivalent mass center;
2) iteratively correcting the force normalization result according to the resultant moment normalization result and the resultant moment expected value taur=[τ1τ2τ3]Finding the optimal solution, tau, with the force normalization result approaching 01,τ2,τ3Respectively obtaining expected values of moments of a rolling shaft, a pitching shaft and a yawing shaft of the load equivalent mass center;
3) determining the force coefficient of the actuator according to the optimal solution of the force normalization result approaching 0 and the resultant moment normalization result, and completing calibration;
normalizing the result according to the resultant moment and the resultant moment expected value taur=[τ1τ2τ3]The method for finding the optimal solution of the acting force normalization result approaching 0 specifically comprises the following steps:
21) determining an index function J according to a resultant moment normalization result and an acting force normalization resultreal
22) Judgment index function JrealWhether the absolute value is less than the critical value, if so, the index function JrealDetermining the normalized result of the acting force as the optimal solution; if not, entering step 23);
23) iteratively correcting the acting force normalization result, and repeating the step 22) until an optimal solution of the acting force normalization result is obtained;
the index function JrealThe method specifically comprises the following steps:
Figure FDA0002309459760000011
Freal,norm(n)=[FrealxFrealyFrealz],
τreal,norm=[τrealxτrealyτrealz],
τr=[τ1τ2τ3],
Freal,norm(0)=Freal,norm
wherein, the lower foot marks realx, realy and realz respectively represent the components on the rolling axis, the pitching axis and the yawing axis of the equivalent center of mass of the load, and tau1,τ2,τ3The method comprises the following steps of (1) respectively obtaining expected values of moments of a rolling shaft, a pitching shaft and a yawing shaft of a load equivalent mass center, wherein n is 0,1 and 2;
step 1) determining the resultant moment normalization result tau actually suffered by the three axes of the equivalent center of mass of the loadreal,normThe method specifically comprises the following steps:
τreal,norm=k2·θpnorm·Ιpnorm,
Figure FDA0002309459760000021
Figure FDA0002309459760000022
wherein, thetap=[θpxpypz]For the load equivalent centroid three-axis coupling attitude measurement, Ipx,Ipy,IpzRespectively are main inertia measured values of a rolling axis, a pitching axis and a yawing axis of the equivalent center of mass of the load, k2The load equivalent mass center rotation stiffness coefficient is obtained according to the stiffness of the platform actuator;
said step 1) determining the forceNormalized result Freal,normThe method specifically comprises the following steps:
Freal,norm=(k1·rpnorm)·m,
Figure FDA0002309459760000023
wherein r ispx,rpy,rpzDesigned values of translational displacement r of load equivalent centroid roll axis, pitch axis and yaw axispx,rpy,rpzThe value ranges of (1) and (0.1) are all (-0.1); m is the load mass, k1The rigidity coefficient of the translational displacement of the three axes of the load equivalent mass center is obtained according to the rigidity of the platform actuator;
the critical value is less than 10-3
2. A calibration method according to claim 1, wherein the method for iteratively correcting the force normalization result specifically comprises:
Figure FDA0002309459760000031
wherein k isrealAnd updating the coefficient for the acting force normalization result, wherein the value range is (-1, 1), and t is time.
3. A calibration method according to claim 2, characterized in that:
the platform is a six-degree-of-freedom Gough-Stewart platform;
the step 3) determines that the actuator force coefficient α is [ α ]1… α6]The specific method comprises the following steps:
Figure FDA0002309459760000032
wherein, βijElement of the ith row and jth column of the β matrix, β ═ J (J)p T)-1,JpIs a Jacobian matrix of equivalent center of mass of the load, when βi,3+kWhen equal to 0, αi=1,i=1,2,...,6,j=1,2,...,6,k=1,2,3。
CN201810716502.2A 2018-06-29 2018-06-29 Platform on-orbit calibration method Active CN108995829B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810716502.2A CN108995829B (en) 2018-06-29 2018-06-29 Platform on-orbit calibration method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810716502.2A CN108995829B (en) 2018-06-29 2018-06-29 Platform on-orbit calibration method

Publications (2)

Publication Number Publication Date
CN108995829A CN108995829A (en) 2018-12-14
CN108995829B true CN108995829B (en) 2020-03-24

Family

ID=64598886

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810716502.2A Active CN108995829B (en) 2018-06-29 2018-06-29 Platform on-orbit calibration method

Country Status (1)

Country Link
CN (1) CN108995829B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111783271B (en) * 2020-05-11 2023-08-29 北京控制工程研究所 Nonlinear correction method for three-override control system of spacecraft
CN113779829B (en) * 2021-08-30 2024-03-19 江苏大学 Rapid temperature calculation method of fault-tolerant permanent magnet motor under turn-to-turn short circuit fault
CN115675919B (en) * 2022-10-31 2024-05-31 北京控制工程研究所 On-orbit calibration method for satellite active pointing ultra-static platform

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103438907B (en) * 2013-09-11 2016-01-20 哈尔滨工业大学 A kind of star sensor six degree of freedom is as the on-orbit calibration method of plane error
CN106114912A (en) * 2016-08-12 2016-11-16 上海卫星工程研究所 GEO track Rotating Platform for High Precision Star Sensor high stable in-orbit points to ensuring method
CN106564622B (en) * 2016-11-02 2019-04-09 北京控制工程研究所 A kind of remote modification method of geosynchronous satellite relative inclination based on visual line measurement

Also Published As

Publication number Publication date
CN108995829A (en) 2018-12-14

Similar Documents

Publication Publication Date Title
CN108995829B (en) Platform on-orbit calibration method
CN109822574B (en) Industrial robot end six-dimensional force sensor calibration method
CN111367236B (en) Mobile robot system calibration method and system for machining process
EP3715780A1 (en) 1d displacement sensor-based 3d measurement model and spatial calibration method
CN108759798B (en) Method for realizing precision measurement of high-precision spacecraft
WO2021128761A1 (en) Distributed multi-dimensional force measurement system and measurement method
CN113459094B (en) Industrial robot tool coordinate system and zero point self-calibration method
CN110715769A (en) Method for calibrating stress point position of weighing sensor of multi-point method centroid measuring equipment
CN110977940B (en) Geometric error modeling method and device for parallel-series robot
CN109976255B (en) Kinematics calibration method for parallel spindle head
CN109733638B (en) Load application method under large deformation condition of long strut undercarriage
CN107246883B (en) High-precision star sensor mounting matrix on-orbit real-time calibration method
CN115229796B (en) Robot tail end pose tracking compensation method for milling aviation large component
CN103921954B (en) Based on the digitalisation calibrating method of the aircraft target ship assembly deflections of three-axis numerical control steady arm
CN113211445A (en) Robot parameter calibration method, device, equipment and storage medium
CN113459104B (en) Robot Positioning Error Prediction Method and Prediction Device Based on Support Vector Regression
CN114393448A (en) Method for improving track precision of magnetorheological robot polishing equipment
CN111605736B (en) Earth-moon L2 point transfer orbit optimal error correction point selection method
CN109214131B (en) Error-optimized static test load design method and system
JPH0433006A (en) Control method for robot system
CN113733155B (en) Six-axis industrial robot calibration device and calibration method
CN109974933A (en) Consider the rotary inertia recursive least-squares discrimination method of satellite flexibility coupling
CN111442786B (en) Zero drift deviation and attitude estimation method of aircraft gyroscope
CN114061571A (en) Attitude calculation method and system of adaptive gradient descent inertial measurement unit
CN109631813B (en) Calibration method of large-size articulated arm type coordinate measuring machine

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