CN102121828A - Method for estimating body posture angle of humanoid robot in real time - Google Patents

Method for estimating body posture angle of humanoid robot in real time Download PDF

Info

Publication number
CN102121828A
CN102121828A CN 201010597601 CN201010597601A CN102121828A CN 102121828 A CN102121828 A CN 102121828A CN 201010597601 CN201010597601 CN 201010597601 CN 201010597601 A CN201010597601 A CN 201010597601A CN 102121828 A CN102121828 A CN 102121828A
Authority
CN
China
Prior art keywords
robot
attitude
posture
angle
axis
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
Application number
CN 201010597601
Other languages
Chinese (zh)
Other versions
CN102121828B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN 201010597601 priority Critical patent/CN102121828B/en
Publication of CN102121828A publication Critical patent/CN102121828A/en
Application granted granted Critical
Publication of CN102121828B publication Critical patent/CN102121828B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Manipulator (AREA)

Abstract

The invention discloses a method for estimating a body posture angle of a humanoid robot in real time, which comprises the following steps of: first establishing a kinematics model of a robot posture; then establishing a noisy linear discrete random system model according to the kinematics model of the robot posture, estimating the posture state of the robot based on the acceleration data of a triaxial accelerometer, the angular velocity data of a three-axis gyroscope and the three-dimensional force data of a six-dimensional force/torque sensor by utilizing a corresponding Kalman filtering method; and obtaining the actual posture angle of the robot according to the estimated posture information of the robot, thereby correcting the current posture of the robot. By the method, the shortcomings of drift of an angular velocity integral and relatively narrower measurement range of the conventional linear method are overcome, and the real-time accurate posture estimation of the humanoid robot is realized.

Description

A kind of anthropomorphic robot body posture angle real-time estimation method
Technical field
The invention belongs to the robot field, relate in particular to a kind of based on sextuple power/torque sensor, accelerometer and gyrostatic anthropomorphic robot attitude method of estimation.
Technical background
Anthropomorphic robot is a kind of robot that imitates the skeleton structure, and it and other machines people's maximum difference is the foot sport mode of its apery.Bionical the analysis showed that, anthropomorphic robot have advantages such as adaptive capacity to environment is strong, flexible movements, therefore, have more wide application and development prospect.
Different with wheeled robot and industrial robot, the supporting domain of anthropomorphic robot is littler, easier falling down.Therefore need know the attitude angle of robot, comprise pitching and roll angle, so that detection machine people's state, the prevention robot is fallen down.
A kind of attitude algorithm for estimating relatively more commonly used is the high-frequency signal that earlier integration of angular velocity is obtained angle through a Hi-pass filter, again to the angle that obtains through acceleration calculation through a low-pass filter, obtain the low frequency signal of attitude angle, the stack of these two angle signals can obtain complete attitude angle (referring to " K. Loffler; M.Gienger; F.Pfeiffer:Sensor and Control Design of a Dynamically Stable Biped Robot; International Conference on Robotics and Automation; September; 2003 "), but this method need be set the parameter of suitable low-pass filter and Hi-pass filter, estimate to reach accurately, and for different walking patterns, need to select different filter parameters, so the scope of application of this method is less.
Because it is non-linear that attitude is estimated, another common method is, attitude motion is learned model carry out linearization process at equilibrium point, obtain the linear model of this kinematical equation, utilize the method for kalman filtering to estimate the state of this linear model, thereby obtain the attitude angle (referring to " An Attitude Control System for small unmanned helicopters based on gyrometer and angular rate gyrosc
Henrik Rehbinder and Xiaoming Hu have proposed the another kind of method (referring to " Drift-free attitude estimation for accelerated rigid bodies, Automatica, 2004 ") that attitude is estimated.In the method, utilize the attitude angle of acceleration of gravity correction robot.But this method needs the transverse acceleration of robot less.For the acceleration change in the anthropomorphic robot walking, because the model error that transverse acceleration is introduced can not ignore, so this method and be not suitable for the walking of anthropomorphic robot.
Chinese invention patent application 200810224579.4 has proposed a humanoid robot foot posture real-time detection and method, and this attitude detecting method utilizes the 3-axis acceleration information calculations to obtain static tilt angles information; The angular velocity information of gathering is calculated the changing value of angle; According to described static face obliquity information and described angle changing value, obtain real-time obliquity information.But since in robot motion's process with the collision on ground, the static tilt angles information of utilizing accelerometer to calculate gained is inaccurate, and this scheme also can not overcome the integral error of angular velocity.
Therefore, be necessary to design a kind of can be real-time, fast and can effectively overcome the anthropomorphic robot attitude angle estimation scheme of gyroscope integral error.
Summary of the invention
The objective of the invention is at the deficiencies in the prior art, a kind of anthropomorphic robot body posture angle real-time estimation method is provided.
The objective of the invention is to be achieved through the following technical solutions: a kind of anthropomorphic robot body posture angle real-time estimation method, the anthropomorphic robot posture estimation system comprises three-axis gyroscope, three axis accelerometer, kalman wave filter and two sextuple power/torque sensors, wherein, three-axis gyroscope, three axis accelerometer and two sextuple power/torque sensors link to each other with the kalman wave filter respectively; This method comprises following step:
(1) sextuple power/torque sensor and three-axis gyroscope and three axis accelerometer are installed; Wherein, two sextuple power/torque sensors are installed in respectively on two ankle-joints of robot, and three axis accelerometer and three-axis gyroscope are installed in the hip central point of robot.
(2) set up the kinematics model of robot pose;
(3) according to the kinematics model of robot pose, set up the linear discrete stochastic system model of band noise, utilize the correlation method of Kalman filtering, based on the acceleration information of three axis accelerometer, the attitude state of the angular velocity data of three-axis gyroscope and sextuple power/torque sensor three-dimensional force data estimation robot;
(4), obtain the attitude angle of robot reality according to the robot pose information of estimating;
(5) according to the attitude angle of the robot reality that obtains, correct the current attitude of robot
The invention has the beneficial effects as follows, anthropomorphic robot body posture angle real-time estimation method of the present invention is realized by the technology of soft measurement, kinematical equation according to robot pose, utilize the method for kalman filtering, obtain the optimal estimation of robot pose state, thereby calculate corresponding pitching and roll angle.This method has overcome the drift and the conventional linear method measurement range smaller defect of angular velocity integration, has realized the estimation accurately in real time of anthropomorphic robot attitude.
Description of drawings
Fig. 1 is the entire block diagram of anthropomorphic robot posture estimation system;
Fig. 2 is the scheme of installation of sensor;
The roll angle data plot of body posture when Fig. 3 is robot stabilized walking;
The luffing angle data plot of body posture when Fig. 4 is robot stabilized walking.
Embodiment
The principle of Kalman filtering is as follows: utilize given system's output observation sequence
Figure 946701DEST_PATH_IMAGE001
, make state
Figure 891523DEST_PATH_IMAGE002
Optimal estimation
Figure 306324DEST_PATH_IMAGE003
, make k+1 estimation index (both estimation error covariance matrixes) J constantly reach minimum:
Figure 45610DEST_PATH_IMAGE004
As shown in Figure 1, the anthropomorphic robot posture estimation system comprises three-axis gyroscope, three axis accelerometer, kalman wave filter and two sextuple power/torque sensors, wherein, three-axis gyroscope, three axis accelerometer and two sextuple power/torque sensors link to each other with the kalman wave filter respectively.Three axis linear acceleration informations when three axis accelerometer provides the robot motion, three-axis gyroscope provides three angular velocity data, and two sextuple power/torque sensors provide the three-dimensional force data.Set up system state equation based on the kinematics model that robot pose changes, according to the force data that sensor obtains, acceleration information and angular velocity data obtain the optimal estimation of robot pose by the method for kalman filtering.Obtain robot pose according to optimal estimation, calculating can get corresponding attitude angle.
Anthropomorphic robot body posture angle real-time estimation method of the present invention comprises following step: at first, sextuple power/torque sensor and three-axis gyroscope and acceierometer sensor are installed; Secondly, set up the kinematics model of robot pose; Once more, kinematics model according to robot pose, set up the linear discrete stochastic system model of band noise, utilize the correlation method of Kalman filtering, based on the acceleration information of three axis accelerometer, the attitude state of the angular velocity data of three-axis gyroscope and sextuple power/torque sensor three-dimensional force data estimation robot; According to the robot pose information after estimating, obtain the attitude angle of robot reality.At last, according to the robot pose pitching and the roll angle that obtain, correct the current attitude of robot.
It is specific as follows,
The first step is installed sextuple power/torque sensor and three-axis gyroscope and acceierometer sensor
The scheme of installation of sensor is as shown in Figure 2: two sextuple power/torque sensors are installed in respectively on two ankle-joints (connecting portion of sole and shank) of robot, and the x axle positive dirction of sensor local coordinate system is along the robot working direction, y axle positive dirction from right to left, z axle positive dirction is straight up.Three axis accelerometer and three-axis gyroscope are installed in the hip central point of robot, and the coordinate system direction of its local coordinate system direction and sextuple power/torque sensor is consistent.
Second goes on foot, and sets up the kinematics model of robot pose.
Attitude can be described with 3 * 3 rotation matrix, corresponding rotation angle around the x axle be called roll angle ( ), corresponding rotation matrix:
Figure 65065DEST_PATH_IMAGE006
The corresponding anglec of rotation around the y axle be called the angle of pitch (
Figure 334372DEST_PATH_IMAGE007
) corresponding rotation matrix is:
The corresponding anglec of rotation around the z axle be called deflection angle ( ), corresponding rotation matrix is:
Figure 784311DEST_PATH_IMAGE010
For the attitude of robot, we with Euler angle by deflection (
Figure 173704DEST_PATH_IMAGE009
), pitching ( ), roll (
Figure 729636DEST_PATH_IMAGE007
) order describe, then Dui Ying attitude matrix is:
Figure 822882DEST_PATH_IMAGE012
To this rotation matrix differentiate:
Figure 66781DEST_PATH_IMAGE013
Be written as two combinations of matrices,
Figure 484173DEST_PATH_IMAGE015
, wherein, in the time interval
Figure 112601DEST_PATH_IMAGE016
In, rotate to be around the trace of axle K
Figure 148690DEST_PATH_IMAGE017
, then following formula can be rewritten as:
Figure 837160DEST_PATH_IMAGE018
Following formula is carried out the angle differential transform gets,
Figure 21017DEST_PATH_IMAGE019
So
Figure 390818DEST_PATH_IMAGE020
Be installed in the three axis accelerometer and the three-axis gyroscope of robot hip central point, robot measurement hip central point is along the acceleration and the corresponding angular velocity of self coordinate system in real time.When robot was static, the acceleration of accelerometer measures was an acceleration of gravity.In the time of the robot motion, the acceleration that the power that the friction force on ground provides robot to advance, robot are subjected to provides (see figure 2) by bonding force.
The sole of supposing robot contacts fully with ground, and robot is a rigidity.We are the initial point of world coordinate system with ground, set up Descartes's right-handed coordinate system, and the positive dirction of x axle is the working direction of robot, and the square first vertical ground of z axle upwards.And the accelerometer of robot hip is identical with the attitude of world coordinate system with the initial attitude of attitude sensor.The bonding force f that robot foot section is subjected to, the acceleration that causes on the accelerometer of robot hip are a, and this moment, the attitude of robot was R, then:
a?=?R*(f/M)
Consider the state transition equation of attitude R, then the kinematics model of robot pose is:
Figure 609310DEST_PATH_IMAGE021
Wherein, M is the quality of robot,
Figure 406365DEST_PATH_IMAGE022
, be making a concerted effort of measuring of two sextuple power/torque sensor, the bonding force that the expression robot foot section is received is at each component of coordinate system,
Figure 77518DEST_PATH_IMAGE023
Be 3 * 1 vector, the three axis angular rate data that the expression three-axis gyroscope is measured are used respectively ,
Figure 55499DEST_PATH_IMAGE025
,
Figure 85772DEST_PATH_IMAGE026
Expression.R is 3 * 3 matrix, the attitude matrix of expression robot;
Figure 244221DEST_PATH_IMAGE027
, the accekeration that the expression three axis accelerometer is measured.
The 3rd step, according to the kinematics model of robot pose, set up the linear discrete stochastic system model of band noise, Utilize the correlation method of Kalman filtering again, based on the acceleration information of three axis accelerometer, the attitude state of the angular velocity data of three-axis gyroscope and sextuple power/torque sensor three-dimensional force data estimation robot.
Consider that R is 3 * 3 matrix, the acceleration a of measurement and bonding force f are 3 * 1 vectors, and be stretching by row to attitude matrix, supposes the noise white noise of sensor measurement simultaneously, and then the system model restructural is:
Figure 955825DEST_PATH_IMAGE028
State variable in the formula
Figure 148909DEST_PATH_IMAGE029
,
Figure 350083DEST_PATH_IMAGE030
Be plant noise, variance is
Figure 667932DEST_PATH_IMAGE031
,
Figure 245544DEST_PATH_IMAGE032
Be observation noise, variance is
Figure 293134DEST_PATH_IMAGE033
, Be 3-axis acceleration,
Figure 470355DEST_PATH_IMAGE035
The Kalman filtering equations of discretize is adopted in the estimation of this model state, then be with Ts discrete time at interval, the discretization equation of this model is:
Figure 589008DEST_PATH_IMAGE036
Wherein,
Figure 428788DEST_PATH_IMAGE037
Figure 971765DEST_PATH_IMAGE038
Figure 529785DEST_PATH_IMAGE039
Be 9 * 1 vector, the plant noise of expression system,
Figure 183620DEST_PATH_IMAGE040
Be 3 * 1 vector, the observation noise of expression system.
Figure 205803DEST_PATH_IMAGE041
Be 9 * 1 vector, the state variable of expression system.
Put forward the back predictand with state
Figure 857364DEST_PATH_IMAGE042
Be the basic variable of recursion, can write out Kalman filtering recursion formula and be:
Figure 699418DEST_PATH_IMAGE043
Figure 156945DEST_PATH_IMAGE044
It is the accekeration of k three axis accelerometer constantly
Kalman forecasts gain matrix For:
Figure 856096DEST_PATH_IMAGE046
Figure 123129DEST_PATH_IMAGE047
Covariance matrix for prediction error:
Figure 384346DEST_PATH_IMAGE048
Figure 112612DEST_PATH_IMAGE049
Be 9 * 9 matrix, the variance of representation model noise,
Figure 168293DEST_PATH_IMAGE050
Be 3 * 3 matrix, the variance of expression observation noise.
In the 4th step,, obtain the attitude angle of robot reality according to the robot pose information after estimating.
Utilize the kalman filtering algorithm, online estimation
Figure 922623DEST_PATH_IMAGE051
, according to the attitude matrix of estimation, calculate the attitude angle of pitching and rolling:
Figure 987530DEST_PATH_IMAGE052
?
In the formula
Figure 510916DEST_PATH_IMAGE053
Represent state vector respectively
Figure 737498DEST_PATH_IMAGE054
The element of the 7th, 8,9 row
In the 5th step,, correct the current attitude of robot according to the robot pose pitching and the roll angle that obtain.
In order to keep the stable operation of robot, the rolling of the setting of robot operational process and luffing angle are 0, and the joint angles of pitching and rolling by the control robot coxa joint utilizes PID to adjust mode, realizes the stable control of robot body posture.
Embodiment 1
Fig. 3 and Fig. 4 show on anthropomorphic robot use method of the present invention and existing integration method difference, the walking speed of setting robot is 0.8km/h, as seen from Figure 3, because the deficiency of robot hip rigidity, robot in the process of walking, there is cyclical variation in roll angle, and these data have reflected that the roll angle in robot motion's process changes.Because the initial position that robot is put or the deviation of zero-bit cause roll angle and non-vanishing.From with the contrast of integration method as can be seen, this method has overcome the drifting problem of angular velocity integration effectively.As seen from Figure 4, after the robot motion finished, the data of luffing angle remained stable.In sum, this method can effectively overcome the drifting problem of angular velocity integration, and the variation of pitching and roll attitude angle realizes effective measurement of robot pitching and roll attitude angle in the reflection robot motion process.

Claims (4)

1. anthropomorphic robot body posture angle real-time estimation method, the anthropomorphic robot posture estimation system comprises three-axis gyroscope, three axis accelerometer, kalman wave filter and two sextuple power/torque sensors, wherein, three-axis gyroscope, three axis accelerometer and two sextuple power/torque sensors link to each other with the kalman wave filter respectively; It is characterized in that this method comprises following step:
(1) sextuple power/torque sensor and three-axis gyroscope and three axis accelerometer are installed; Wherein, two sextuple power/torque sensors are installed in respectively on two ankle-joints of robot, and three axis accelerometer and three-axis gyroscope are installed in the hip central point of robot;
(2) set up the kinematics model of robot pose;
(3) according to the kinematics model of robot pose, set up the linear discrete stochastic system model of band noise, utilize the correlation method of Kalman filtering, based on the acceleration information of three axis accelerometer, the attitude state of the angular velocity data of three-axis gyroscope and sextuple power/torque sensor three-dimensional force data estimation robot;
(4), obtain the attitude angle of robot reality according to the robot pose information of estimating;
(5) according to the attitude angle of the robot reality that obtains, correct the current attitude of robot
2. according to the described anthropomorphic robot body posture of claim 1 angle real-time estimation method, it is characterized in that the kinematics model that described step (2) is set up is:
Figure 576370DEST_PATH_IMAGE001
Wherein, M is the quality of robot,
Figure 664411DEST_PATH_IMAGE002
, be making a concerted effort of measuring of two sextuple power/torque sensor, the bonding force that the expression robot foot section is received is at each component of coordinate system,
Figure 993762DEST_PATH_IMAGE003
Be 3 * 1 vector, the three axis angular rate data that the expression three-axis gyroscope is measured are used respectively ,
Figure 658278DEST_PATH_IMAGE005
,
Figure 979538DEST_PATH_IMAGE006
Expression; R is 3 * 3 matrix, the attitude matrix of expression robot;
Figure 530605DEST_PATH_IMAGE007
, the accekeration that the expression three axis accelerometer is measured.
3. according to the described anthropomorphic robot body posture of claim 1 angle real-time estimation method, it is characterized in that the linear discrete stochastic system model that described step (2) is set up is:
Figure 533196DEST_PATH_IMAGE008
Wherein, state variable
Figure 118898DEST_PATH_IMAGE009
,
Figure 342551DEST_PATH_IMAGE010
Be plant noise, variance is
Figure 584176DEST_PATH_IMAGE011
,
Figure 187196DEST_PATH_IMAGE012
Be observation noise, variance is
Figure 627404DEST_PATH_IMAGE013
,
Figure 228150DEST_PATH_IMAGE014
Be 3-axis acceleration:
Figure 753809DEST_PATH_IMAGE015
4. according to the described anthropomorphic robot body posture of claim 1 angle real-time estimation method, it is characterized in that the attitude angle of the robot reality that described step (4) obtains is:
Figure 347471DEST_PATH_IMAGE016
Wherein,
Figure 910695DEST_PATH_IMAGE017
Represent state vector respectively
Figure 416762DEST_PATH_IMAGE018
The element of the 7th, 8,9 row.
CN 201010597601 2010-12-21 2010-12-21 Method for estimating body posture angle of humanoid robot in real time Expired - Fee Related CN102121828B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010597601 CN102121828B (en) 2010-12-21 2010-12-21 Method for estimating body posture angle of humanoid robot in real time

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010597601 CN102121828B (en) 2010-12-21 2010-12-21 Method for estimating body posture angle of humanoid robot in real time

Publications (2)

Publication Number Publication Date
CN102121828A true CN102121828A (en) 2011-07-13
CN102121828B CN102121828B (en) 2012-12-19

Family

ID=44250435

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010597601 Expired - Fee Related CN102121828B (en) 2010-12-21 2010-12-21 Method for estimating body posture angle of humanoid robot in real time

Country Status (1)

Country Link
CN (1) CN102121828B (en)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102670217A (en) * 2012-05-04 2012-09-19 嘉兴市制衡精仪有限公司 Wearable sensor measuring device and method for lower limb joint acting force and moment
CN103808322A (en) * 2014-03-04 2014-05-21 韩涛 Compound type six degrees of freedom input method and system in inertia measurement LED (light-emitting diode) detection
CN104570732A (en) * 2014-12-15 2015-04-29 浙江大学 Online generation method for walking mode of humanoid robot
CN105956617A (en) * 2016-04-27 2016-09-21 西北工业大学 Joint estimation method for people and bike attitude in bike-like sport activity
CN106492439A (en) * 2016-10-08 2017-03-15 武汉新英赛健康科技有限公司 A kind of work-break exercises methods of marking
CN106650591A (en) * 2016-10-08 2017-05-10 武汉新英赛健康科技有限公司 Work-break exercise scoring system
CN107456173A (en) * 2016-06-06 2017-12-12 北京小米移动软件有限公司 Barrier crossing method and device
CN107728635A (en) * 2017-11-13 2018-02-23 北京赛曙科技有限公司 A kind of automatic balancing arrangement and method of motor cycle type robot
CN107883951A (en) * 2017-10-19 2018-04-06 福建海图智能科技有限公司 A kind of computational methods and terminal of underwater robot 3 d pose
CN108345211A (en) * 2017-01-23 2018-07-31 深圳市祈飞科技有限公司 Biped anthropomorphic robot and its non-linear gait planning method and control method
CN108646733A (en) * 2018-04-27 2018-10-12 杭州艾豆智能科技有限公司 A kind of mobile robot and its antidote of automatic straightening
CN108897220A (en) * 2018-07-12 2018-11-27 上海硅族智能科技有限公司 A kind of adaptive stabilizing balance control method and system and biped anthropomorphic robot
CN109645960A (en) * 2019-01-15 2019-04-19 浙江大学 The physiological parameter generating device and its method of anthropomorphic robot
CN109991989A (en) * 2017-12-29 2019-07-09 深圳市优必选科技有限公司 Dynamic balance method and device of robot in idle state and storage medium
CN111956449A (en) * 2020-08-10 2020-11-20 河海大学常州校区 Exoskeleton rehabilitation treatment device for shoulder-elbow joint injury and control system thereof
CN112631277A (en) * 2020-12-08 2021-04-09 中山大学 Balance control method and system for four-legged robot standing posture conversion
CN112617811A (en) * 2021-01-04 2021-04-09 杭州风行医疗器械有限公司 Joint biomechanics measurement and evaluation device and method
WO2022199058A1 (en) * 2021-03-22 2022-09-29 深圳市优必选科技股份有限公司 Processing method and apparatus for robot stride length, robot control device, and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09273938A (en) * 1996-04-04 1997-10-21 Takao Yamaguchi Hybrid inertial navigation apparatus
WO1999020980A1 (en) * 1997-10-20 1999-04-29 Churyo Engineering Co., Ltd. Position measurement method for mobile station
CN101100059A (en) * 2007-07-27 2008-01-09 北京工业大学 Flexible double-wheel self-balancing robot attitude detecting method
CN101403619A (en) * 2008-10-21 2009-04-08 北京理工大学 Foot gesture real-time detection system and method for human-imitated robot
CN101672650A (en) * 2009-09-29 2010-03-17 北京航空航天大学 Orienting and locating navigation system in circumstance of electromagnetic interference

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09273938A (en) * 1996-04-04 1997-10-21 Takao Yamaguchi Hybrid inertial navigation apparatus
WO1999020980A1 (en) * 1997-10-20 1999-04-29 Churyo Engineering Co., Ltd. Position measurement method for mobile station
CN101100059A (en) * 2007-07-27 2008-01-09 北京工业大学 Flexible double-wheel self-balancing robot attitude detecting method
CN101403619A (en) * 2008-10-21 2009-04-08 北京理工大学 Foot gesture real-time detection system and method for human-imitated robot
CN101672650A (en) * 2009-09-29 2010-03-17 北京航空航天大学 Orienting and locating navigation system in circumstance of electromagnetic interference

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102670217A (en) * 2012-05-04 2012-09-19 嘉兴市制衡精仪有限公司 Wearable sensor measuring device and method for lower limb joint acting force and moment
CN103808322A (en) * 2014-03-04 2014-05-21 韩涛 Compound type six degrees of freedom input method and system in inertia measurement LED (light-emitting diode) detection
CN104570732A (en) * 2014-12-15 2015-04-29 浙江大学 Online generation method for walking mode of humanoid robot
CN104570732B (en) * 2014-12-15 2017-02-08 浙江大学 Online generation method for walking mode of humanoid robot
CN105956617A (en) * 2016-04-27 2016-09-21 西北工业大学 Joint estimation method for people and bike attitude in bike-like sport activity
CN105956617B (en) * 2016-04-27 2019-05-17 西北工业大学 A kind of people's vehicle posture combined estimation method for class cycling
CN107456173B (en) * 2016-06-06 2020-12-04 北京小米移动软件有限公司 Obstacle crossing method and device
CN107456173A (en) * 2016-06-06 2017-12-12 北京小米移动软件有限公司 Barrier crossing method and device
CN106650591A (en) * 2016-10-08 2017-05-10 武汉新英赛健康科技有限公司 Work-break exercise scoring system
CN106492439A (en) * 2016-10-08 2017-03-15 武汉新英赛健康科技有限公司 A kind of work-break exercises methods of marking
CN106650591B (en) * 2016-10-08 2019-12-31 温胜洁 Workshop operation scoring system
CN108345211A (en) * 2017-01-23 2018-07-31 深圳市祈飞科技有限公司 Biped anthropomorphic robot and its non-linear gait planning method and control method
CN107883951A (en) * 2017-10-19 2018-04-06 福建海图智能科技有限公司 A kind of computational methods and terminal of underwater robot 3 d pose
CN107728635B (en) * 2017-11-13 2020-10-09 北京赛曙科技有限公司 Automatic balancing device and method for motorcycle type robot
CN107728635A (en) * 2017-11-13 2018-02-23 北京赛曙科技有限公司 A kind of automatic balancing arrangement and method of motor cycle type robot
CN109991989A (en) * 2017-12-29 2019-07-09 深圳市优必选科技有限公司 Dynamic balance method and device of robot in idle state and storage medium
CN108646733A (en) * 2018-04-27 2018-10-12 杭州艾豆智能科技有限公司 A kind of mobile robot and its antidote of automatic straightening
CN108897220A (en) * 2018-07-12 2018-11-27 上海硅族智能科技有限公司 A kind of adaptive stabilizing balance control method and system and biped anthropomorphic robot
CN108897220B (en) * 2018-07-12 2020-09-29 上海硅族智能科技有限公司 Self-adaptive stable balance control method and system and biped humanoid robot
CN109645960A (en) * 2019-01-15 2019-04-19 浙江大学 The physiological parameter generating device and its method of anthropomorphic robot
CN109645960B (en) * 2019-01-15 2023-05-23 浙江大学 Physiological parameter generating device and method for humanoid robot
CN111956449A (en) * 2020-08-10 2020-11-20 河海大学常州校区 Exoskeleton rehabilitation treatment device for shoulder-elbow joint injury and control system thereof
CN112631277A (en) * 2020-12-08 2021-04-09 中山大学 Balance control method and system for four-legged robot standing posture conversion
CN112617811A (en) * 2021-01-04 2021-04-09 杭州风行医疗器械有限公司 Joint biomechanics measurement and evaluation device and method
WO2022199058A1 (en) * 2021-03-22 2022-09-29 深圳市优必选科技股份有限公司 Processing method and apparatus for robot stride length, robot control device, and storage medium

Also Published As

Publication number Publication date
CN102121828B (en) 2012-12-19

Similar Documents

Publication Publication Date Title
CN102121828A (en) Method for estimating body posture angle of humanoid robot in real time
CN106705968B (en) Indoor inertial navigation algorithm based on attitude identification and step size model
CN110986939B (en) Visual inertia odometer method based on IMU (inertial measurement Unit) pre-integration
CN111024064B (en) SINS/DVL combined navigation method for improving Sage-Husa adaptive filtering
CN104764451A (en) Target posture tracking method based on inertia and geomagnetic sensor
Jin et al. The adaptive Kalman filter based on fuzzy logic for inertial motion capture system
CN107941217A (en) A kind of robot localization method, electronic equipment, storage medium, device
CN106643715A (en) Indoor inertial navigation method based on bp neural network improvement
Marwan et al. Calibration method for articulated industrial robots
CN111086001A (en) State estimation method and system for multi-modal perception of foot robot
US11685049B2 (en) Robot localization using variance sampling
CN112326990A (en) Method and system for measuring speed of working vehicle
CN111618861A (en) Double-follow-up intelligent arm control method based on four-axis structure
CA3086559C (en) Method for predicting a motion of an object, method for calibrating a motion model, method for deriving a predefined quantity and method for generating a virtual reality view
CN108507572A (en) A kind of attitude orientation error correcting method based on MEMS Inertial Measurement Units
CN110737194B (en) Multi-MEMS sensor combination attitude measurement method
CN112650217B (en) Robot trajectory tracking strategy dynamic optimization method based on evaluation function
CN111679681B (en) Underwater robot motion path tracking method based on sliding mode robust control
CN117091592A (en) Gesture resolving method, gesture resolving device, and computer storage medium
CN116338719A (en) Laser radar-inertia-vehicle fusion positioning method based on B spline function
CN108107882B (en) Automatic calibration and detection system of service robot based on optical motion tracking
CN115077517A (en) Low-speed unmanned vehicle positioning method and system based on fusion of stereo camera and IMU
CN113190025B (en) Motion control method suitable for variable-structure underwater robot
CN114789439A (en) Slope positioning correction method and device, robot and readable storage medium
CN113358121A (en) Electrically-driven insect configuration foot type robot foot-to-ground slip estimation method

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20110713

Assignee: ZHEJIANG GUOZI ROBOT TECHNOLOGY Co.,Ltd.

Assignor: Zhejiang University

Contract record no.: 2013330000081

Denomination of invention: Method for estimating body posture angle of humanoid robot in real time

Granted publication date: 20121219

License type: Exclusive License

Record date: 20130422

LICC Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20121219

Termination date: 20211221