CN101473193A - 姿势角检测装置与姿势角检测方法 - Google Patents

姿势角检测装置与姿势角检测方法 Download PDF

Info

Publication number
CN101473193A
CN101473193A CNA200780023332XA CN200780023332A CN101473193A CN 101473193 A CN101473193 A CN 101473193A CN A200780023332X A CNA200780023332X A CN A200780023332XA CN 200780023332 A CN200780023332 A CN 200780023332A CN 101473193 A CN101473193 A CN 101473193A
Authority
CN
China
Prior art keywords
acceleration
component
angle
moving body
moved
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
CNA200780023332XA
Other languages
English (en)
Other versions
CN101473193B (zh
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.)
Toyota Motor Corp
Original Assignee
Toyota Motor Corp
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 Toyota Motor Corp filed Critical Toyota Motor Corp
Publication of CN101473193A publication Critical patent/CN101473193A/zh
Application granted granted Critical
Publication of CN101473193B publication Critical patent/CN101473193B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C9/00Measuring inclination, e.g. by clinometers, by levels
    • G01C9/02Details
    • G01C9/08Means for compensating acceleration forces due to movement of instrument
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/12Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to parameters of the vehicle itself, e.g. tyre models
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P15/00Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration
    • G01P15/18Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration in two or more dimensions
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • B60W40/105Speed
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • B60W40/107Longitudinal acceleration
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • B60W40/109Lateral acceleration
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • B60W40/11Pitch movement
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • B60W40/112Roll movement
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • B60W40/114Yaw movement

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Mathematical Physics (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Navigation (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Gyroscopes (AREA)

Abstract

提供通过对来自加速度传感器的输出值进行补正而高精度地检测移动体的移动中的姿势角的姿势角检测装置及其方法。其特征在于,包括:加速度传感器(2),其测定施加在移动体(7)上的加速度;横摆率传感器(3),其测定移动体(7)的横摆率;速度传感器(4),其测定移动体(7)的速度;移动成分加速度计算单元(5),其从速度计算实际加速度,从横摆率和速度计算离心力,并计算作为实际加速度与离心力的合力的移动成分加速度;和姿势角计算单元(6),其从用移动成分加速度对加速度进行补正而得到的重力加速度计算姿势角。

Description

姿势角检测装置与姿势角检测方法
技术领域
本发明涉及移动体的姿势角检测装置及其方法。
背景技术
作为得知车辆等移动体本身的位置、速度的装置(导航装置)的一种,有不需要地面上的设备的辅助的惯性导航装置(INS,Inertial NavigationSystem)。使用于该INS的惯性传感器包括检测速度的变化的加速度传感器与检测姿势的变化的陀螺传感器。INS使用惯性传感器并用计算机将从出发点开始的速度、姿势的变化累积,由此能够得知自己的位置、速度。
通常,为了确定移动体的姿势状态,使用将移动体的前后(轴)方向设为x(侧倾,roll)轴、将左右(轴)方向设为y(纵倾,pitch)轴、将上下(轴)方向设为z(横摆,yaw)轴的正交3轴坐标系。在移动体被放置在水平的地面上时,移动体的横摆轴大致朝向地球的铅直轴(Z轴),侧倾轴以及纵倾轴与地球的局部水平面(X-Y平面)大致平行。而且,将侧倾轴从局部水平面X轴倾斜的角度设为纵倾角,将纵倾轴从局部水平面Y轴倾斜的角度设为侧倾角、将侧倾轴或者纵倾轴在X-Y平面上偏离的角设为横摆角,通过这3个角度表示姿势角。为了测定这3个角度,使用3轴式陀螺传感器。
例如,专利文献1的发明对3轴式陀螺传感器与加速度传感器的输出值分别进行加权,计算姿势角。陀螺传感器的输出结果为了姿势角检测而被积分。但是,陀螺传感器具有随着时间的经过而向相同方向输出误差的现象(漂移,drift),如果对包含误差的值进行积分,则误差会被累积起来。而且,陀螺传感器的价格较高,所以为了以低成本满足功能,具有通过加速度传感器检测重力(gravity)来计算姿势角的方法。
专利文献2的发明基于来自2个加速度传感器的输出结果检测移动体的姿势角。但是,加速度传感器的输出包含伴随着移动体的移动的加速度,所以在专利文献2所示的发明中,难以检测出正确的姿势角。
专利文献1:特开平8-178687号公报
专利文献2:特开2003-139536号公报
发明内容
本发明是鉴于上述课题而进行的,其目的在于提供通过对来自加速度传感器的输出值进行补正而高精度地检测移动体的移动中的姿势角的姿势角检测装置及其方法。
由于在仅从施加在移动体上的加速度计算姿势角时,移动体的移动中的姿势角难以被正确地检测,所以本发明者着眼于从移动体的速度或角速度求得与移动体的移动有关的加速度的分力(离心力等),完成了本发明。
本发明的姿势角检测装置,其特征在于,包括:加速度传感器,该加速度传感器测定施加于移动体的加速度;横摆率传感器,该横摆率传感器测定该移动体的横摆率;速度传感器,该速度传感器测定该移动体的速度;移动成分加速度计算单元,该移动成分加速度计算单元从该速度计算实际加速度,从该横摆率和该速度计算离心力,并计算作为该实际加速度与该离心力的合力的移动成分加速度;和姿势角计算单元,该姿势角计算单元从用该移动成分加速度对该加速度进行补正而得到的重力加速度计算姿势角。
在本发明的姿势角检测装置中,加速度传感器对于侧倾轴方向、纵倾轴方向、横摆轴方向这3个方向测定施加在移动体上的加速度,横摆率传感器测定移动体的对于横摆轴方向的旋转的横摆率(角速度)。而且,通过从移动体的速度计算求得与移动(运动)有关的加速度(运动量),设为实际加速度。同样,使用移动体的速度以及横摆率计算与移动体有关的离心力。实际加速度与离心力是由于移动体的移动而产生的加速度(移动成分加速度),包含于由加速度传感器得到的加速度。因此,通过从加速度中除去移动成分加速度,能够将与移动体的移动有关的运动量除去而得到重力加速度。从而能够使用重力加速度求得姿势角。
在本发明的姿势角检测装置中,通过检测移动体的正交的3个轴的加速度与移动体的速度,只要使用用于检测横摆角的1个陀螺传感器,便能够求得移动体的姿势角。因此,能够减少高价的陀螺传感器的个数。但是,陀螺传感器的个数并不限定于1个,通过在正交的3个轴上使用并且从加速度与速度检测移动体的姿势角,能够得到更正确的检测结果。另外,通过不将速度传感器以及加速度传感器的输出结果直接使用于姿势角的检测,而是计算并从加速度传感器的加速度中去除伴随着移动体的移动而产生的加速度(运动量),能够计算正确的姿势角。
优选的是,在本发明的姿势角检测装置中所使用的姿势角计算单元包括:重力加速度计算单元,该重力加速度计算单元计算从加速度中除去移动成分加速度之后的重力加速度;侧倾角计算单元,该侧倾角计算单元从重力加速度计算移动体的侧倾角;和纵倾角计算单元,该纵倾角计算单元从重力加速度计算移动体的纵倾角。移动体的姿势角通过侧倾角与纵倾角表示。通过重力加速度计算单元可求得各方向的重力加速度,所以能够从纵倾轴方向的重力加速度计算侧倾角,从侧倾轴方向的重力加速度计算纵倾角。由于能够使用该重力分力对侧倾角与纵倾角求取移动体的姿势角,所以能够检测更正确的姿势角。
优选的是,在本发明的姿势角检测装置中所使用的移动成分加速度计算单元包括:前后轴加速度计算单元,该前后轴加速度计算单元对速度进行时间微分而计算实际加速度,作为移动成分加速度的前后轴方向的分力;左右轴加速度计算单元,该左右轴加速度计算单元从横摆率和速度计算离心力,作为移动成分加速度的左右轴方向的分力;和上下轴加速度计算单元,该上下轴加速度计算单元计算离心力的在移动体的上下轴方向上的分力,作为移动成分加速度的上下轴方向的分力。移动成分加速度是与移动体的移动(运动)有关的运动量,表示与移动体的正交的3个轴相对应的分力。因此,由于通过对移动体的正交的3个轴中的每个求得并从关于移动体的加速度中去除,来检测移动体的姿势角,所以能够求得更正确的重力加速度。
优选的是,在本发明的姿势角检测装置中所使用的上下轴加速度计算单元进行的离心力的分力的计算基于由姿势角计算单元计算的侧倾角或者纵倾角进行;姿势角计算单元基于计算出的离心力的分力计算侧倾角以及纵倾角中的至少一方。在上下轴加速度计算单元,基于侧倾角或者纵倾角进行离心力的分力的计算,通过使离心力的分力反映到纵倾角以及侧倾角的计算,能够更正确地求得移动体的姿势角。
本发明的姿势角检测方法,其特征在于,具有:测定步骤,该测定步骤测定施加于移动体的加速度、该移动体的横摆率和该移动体的速度;移动成分加速度计算步骤,该移动成分加速度计算步骤从该速度计算实际加速度,从该横摆率和该速度计算离心力,并计算作为该实际加速度与该离心力的合力的移动成分加速度;和姿势角计算步骤,该姿势角计算步骤从用该移动成分加速度对该加速度进行补正而得到的重力加速度计算姿势角。
本发明的姿势角检测方法中,对于侧倾轴方向、纵倾轴方向、横摆轴方向这3个方向由加速度传感器测定施加在移动体上的加速度,横摆率传感器测定移动体的对于横摆轴方向的旋转的横摆率。而且,通过从移动体的速度计算求得与移动有关的加速度,设为实际加速度。同样,使用移动体的速度计算与移动体有关的离心力。实际加速度与离心力是由于移动体的移动而产生的加速度(移动成分加速度),包含于由加速度传感器得到的加速度。因此,通过从加速度中除去移动成分加速度,能够将与移动体的移动有关的运动量除去而得到重力加速度。从而能够使用重力加速度求得姿势角。
在本发明的姿势角检测方法中,能够通过检测移动体的速度以及正交的3个轴的加速度,求得移动体的姿势角,所以只需要1个检测横摆角的陀螺传感器,能够减少高价的陀螺传感器的个数。另外,通过不将速度传感器以及加速度传感器的输出结果直接使用于姿势角的检测,而计算并从加速度传感器的加速度中去除伴随着移动体的移动而产生的加速度,能够计算正确的姿势角。
优选的是,在本发明的姿势角检测方法中所使用的姿势角计算步骤包括:重力加速度计算步骤,该重力加速度计算步骤计算从加速度中除去移动成分加速度之后的重力加速度;侧倾角计算步骤,该侧倾角计算步骤从重力加速度计算移动体的侧倾角;和纵倾角计算步骤,该纵倾角计算步骤从重力加速度计算移动体的纵倾角。移动体的姿势角通过侧倾角与纵倾角表示。通过重力加速度计算单元,求得各方向的重力加速度,所以能够从纵倾轴方向的重力加速度计算侧倾角,从侧倾角轴方向的重力加速度计算纵倾角。由于能够使用该重力分力对侧倾角与纵倾角计算移动体的姿势角,所以能够检测更正确的姿势角。
优选的是,在本发明的姿势角检测方法中所使用的移动成分加速度计算步骤包括:前后轴加速度计算步骤,该前后轴加速度计算步骤对速度进行时间微分而计算实际加速度,作为移动成分加速度的前后轴方向的分力;左右轴加速度计算步骤,该左右轴加速度计算步骤从横摆率和速度计算离心力,作为移动成分加速度的左右轴方向的分力;和上下轴加速度计算步骤,该上下轴加速度计算步骤计算离心力的在移动体的上下方向上的分力,作为移动成分加速度的上下方向的分力。移动成分加速度是与移动体的移动(运动)有关的运动量,表示与移动体的正交的3个轴相对应的分力。因此,由于通过对移动体的正交的3个轴中的每个求得并从涉及移动体的加速度中去除,来检测移动体的姿势角,所以能够求得更正确的重力加速度。
优选的是,在本发明的姿势角检测方法中所使用的重力加速度计算步骤是这样的步骤,将重力加速度分为前后轴方向、左右轴方向以及上下轴方向的分力进行计算,基于在侧倾角计算步骤以及纵倾角计算步骤中的一方中计算出的侧倾角或者纵倾角计算各分力中的一部分。而且,优选的是,侧倾角计算步骤以及纵倾角计算步骤中的一方基于各分力中的残余部分计算侧倾角或者纵倾角,侧倾角计算步骤以及纵倾角计算步骤中的另一方基于各分力中的一部分乃至全部计算侧倾角或者纵倾角。在本发明的姿势角检测方法中由于计算侧倾角以及纵倾角,从而基于侧倾角或者纵倾角计算通过重力加速度计算步骤计算的重力加速度的各轴方向的分力。由此,通过使重力加速度的分力反映到纵倾角以及侧倾角的计算,能够更正确地求得移动体的姿势角。
根据本发明的姿势角检测装置及其方法,不必在正交的3个轴中的所有轴上使用高价的陀螺传感器,使用1个便能起作用,所以能够以较低的成本检测移动体的姿势角。另外,陀螺传感器的输出使用角速度(横摆率),在计算中不使用包含由陀螺传感器的漂移引起的误差的值,所以能够检测出正确的姿势角。
而且,不是将加速度传感器的输出值直接用于计算,而是求得去除了与移动体的运动有关的加速度的重力加速度,所以即使移动体移动也能够正确地检测姿势角。
附图说明
图1是本实施例的姿势角检测装置的结构图;
图2是表示本实施例的姿势角检测装置以及方法中所使用的绝对轴与相对轴的图;
图3是本实施例的姿势角检测方法的流程图;
图4是表示本实施例的移动体的顺侧倾(normal bank)状态的图;
图5是表示本实施例的移动体的逆侧倾(reverse bank)状态的图;
图6是表示本实施例的姿势角检测方法中所使用的判定侧倾(bank)与转弯的流程图。
符号说明
1:姿势角检测装置
2:加速度传感器
3:横摆率(yaw rate)传感器
4:速度传感器
5:移动加速度成分计算单元
6:姿势角计算单元
7:移动体
51:前后轴加速度计算单元
52:左右轴加速度计算单元
53:上下轴加速度计算单元
61:重力加速度计算单元
62:侧倾角计算单元
63:纵倾角计算单元
71:纵倾角(pitch angle)
72:侧倾角(roll angle)
73:横摆角(yaw angle)
81:移动体7的上下轴(横摆轴)
82:移动体7的前后轴(侧倾轴)
83:移动体7的左右轴(纵倾轴)
91:地球的铅直轴(Z轴)
92:局部水平面X轴
93:局部水平面Y轴
94:地球的水平面(X-Y平面)
具体实施方式
下面,使用实施例对本发明进行具体说明。
(实施例)
图1是本实施例的姿势角检测装置的结构图。
姿势角检测装置1包括:加速度传感器2,横摆率传感器3,速度传感器4,移动加速度成分计算单元5和姿势角计算单元6。移动加速度成分计算单元5和姿势角计算单元6作为计算机上的逻辑而被实现。
安装姿势角检测装置1的移动体7如图2所示,在将移动体7放置在水平的地面上时,移动体7的上下轴(横摆轴)81大致朝向地球的铅直轴(Z轴)91方向,移动体7的前后轴(侧倾轴)82以及左右轴(纵倾轴)83与地球的水平面(X-Y平面)94大致平行。将移动体7的侧倾轴82从局部水平面X轴92倾斜的角度设为纵倾角71,将纵倾轴83从局部水平面Y轴93倾斜的角度设为侧倾角72,将侧倾轴82或者纵倾轴83在X-Y平面上偏离的角设为横摆角73,通过这3个角度表示移动体7的姿势角。
加速度传感器2包括分别测定侧倾轴82、纵倾轴83以及横摆轴81即正交的这3个轴的加速度的Ax21、Ay22、Ax23。加速度传感器2测定移动体的加速度,向移动加速度成分计算单元5以及姿势角计算单元6输出。加速度传感器2被设置成:正交的3个轴中的侧倾轴82与移动体7的前后方向一致、纵倾轴83与移动体7的左右方向一致、横摆轴81与移动体7的上下方向一致,与移动体7的旋转同时地进行旋转。
横摆率传感器3是测定作为移动体7的上下轴的横摆轴81的角速度(横摆率)并向移动加速度成分计算单元5输出的陀螺传感器。速度传感器4是测定移动体7的速度并向移动加速度成分计算单元5输出的单元。在移动体为车辆的情况下,测定车轮速度。
移动成分加速度计算单元5是基于从横摆率传感器3以及速度传感器4输入的横摆率以及速度对正交的3个轴中的每个计算与移动体7的移动有关的加速度的单元,包括:前后轴加速度计算单元51,左右轴加速度计算单元52和上下轴加速度计算单元53。移动成分加速度计算单元5将计算结果向姿势角计算单元6输出。
前后轴加速度计算单元51是计算与移动体7的移动有关的前后方向的加速度的单元。左右轴加速度计算单元52是计算与移动体7的移动有关的左右方向的加速度的单元。上下轴加速度计算单元53是计算与移动体7的移动有关的上下方向的加速度的单元。
姿势角计算单元6包括:重力加速度计算单元61,侧倾角计算单元62和纵倾角计算单元63。重力加速度计算单元61是从由加速度传感器2输入的加速度计算使用由移动成分加速度计算单元5输入的移动成分加速度而进行了补正的重力加速度的单元。侧倾角计算单元62是基于从重力加速度计算单元61输入的重力加速度计算侧倾角72的单元。纵倾角计算单元63是基于从重力加速度计算单元61输入的重力加速度计算纵倾角73的单元。
图3是通过流程图表示在本实施例的姿势角检测装置1中所使用的姿势角检测方法(下面,称作姿势角检测例程)的图。本例程在姿势角检测装置1起动后,在预定期间中每一预定时间执行一次。
在本例程的处理开始后,加速度传感器2测定施加在移动体7上的加速度,横摆率传感器3测定横摆率(ω),速度传感器4测定速度(Vw),向移动成分加速度计算单元5以及姿势角计算单元6输出(测定步骤S1)。加速度传感器2通过Ax21测定侧倾轴方向的加速度(Ax_sensor)、通过Ay22测定纵倾轴方向的加速度(Ay_sensor)、通过Az23测定横摆轴方向的加速度(Az_sensor)。
接下来,移动成分加速度计算单元5根据横摆率以及速度进行与移动体7的移动(运动)有关的加速度(移动成分加速度)的计算(移动成分加速度计算步骤)。
然后,姿势角计算单元6用重力加速度进行侧倾角(Roll)72、纵倾角(Pitch)71的计算(姿势角计算步骤),其中该重力加速度是用从移动成分加速度计算单元5输入的移动成分加速度对由加速度传感器2测定的加速度进行补正而求得的。
移动成分加速度计算单元5具有前后轴加速度计算步骤S21,在该步骤中,前后轴加速度计算单元51对速度(Vw)进行时间微分而计算实际加速度,设为移动成分加速度的前后方向的分力(Ax_cor)。而且,具有左右轴加速度计算步骤22,在该步骤中,左右轴加速度计算单元52根据横摆率(ω)和速度(Vw)计算离心力,设为移动成分加速度的左右方向的分力(Ay_cor);并具有上下轴加速度计算步骤23,在该步骤中,上下轴加速度计算单元53计算离心力在移动体7的上下方向上的分力而设为移动成分加速度的上下方向的分力(Az_cor)。
离心力(Ay_cor)可以通过式1求得。
式1:Ay_cor=速度(Vw)×横摆率(ω)
移动成分加速度的上下方向的分力(Az_cor)基于通过姿势角计算单元6的姿势角计算步骤计算出的Roll而计算。
姿势角计算单元6具有:重力加速度计算单元61计算用移动成分加速度对加速度进行补正后的重力加速度的重力加速度计算步骤S31,根据重力加速度计算移动体7的侧倾角(Roll)的侧倾角计算步骤S32,和根据重力加速度计算移动体7的纵倾角(Pitch)的纵倾角计算步骤S33。
在步骤S32中,根据移动体7倾斜的方向,Ay_sensor在图4所示的顺侧倾时与图5所示的逆侧倾时通过不同的公式进行计算。
在顺侧倾时:
式2:Ay_sensor=Ay_cor-g×sin(Roll)
g=该地域的重力加速度
在逆侧倾时:
式3:Ay_sensor=Ay_cor+g×sin(Roll)
g=该地域的重力加速度
式2以及式3的g·sin(Roll)是移动成分加速度的左右方向的分力。在为了求得左右方向的分力而对式2以及式3进行移项时,根据移动体7向哪个方向倾斜,对于在所测定的Ay_sensor上施加离心力还是从其减去离心力而言是不同的。因此,如图6所示,将Ay_sensor的大小与离心力的大小相比,在离心力较小时为顺侧倾,在离心力较大时为逆侧倾。在顺侧倾时将1设为侧倾标志(B_flag)的值,在逆侧倾时将-1设为侧倾标志(B_flag)的值。另外,对于移动体7的转弯(turn),根据左转弯和右转弯的不同,Ay_sensor与离心力的力的作用方向也不同,所以如图6所示,通过ω进行判断,在右转弯时将1设为方向标志(RL_flag)的值,在左转弯时将-1设为方向标志(RL_flag)的值。重力加速度g为常数,所以只要能够求得左右方向的分力,就能够计算Roll(侧倾角)。
式4
Roll=sin-1((|Ay_sensor|-|Vw×ω|)/g)×RL_flag×B_flag
接下来,在S33中,根据重力加速度来计算Pitch(纵倾角),其中该重力加速度是通过移动成分加速度的正交的3个轴的各分力(Ax_cor、Ay_cor、Az_cor)对加速度(Ax_sensor、Ay_sensor、Az_sensor)进行补正而得到的。
式5
Pitch=tan-1((Ax_sensor-Ax_cor)/((Ay_sensor-Ay_cor)2+
            (Az_sensor-Az_cor)2))1/2)
移动成分加速度的正交的3个轴的各分力(Ax_cor、Ay_cor、Az_cor)通过移动成分加速度计算单元5的步骤S21、S22、S23而求得。
根据本实施例的姿势角检测装置1及其方法,通过使用加速度传感器与速度传感器,不必在正交的3个轴中的所有轴上使用高价的陀螺传感器,能够以较低的成本检测移动体的姿势角。另外,也不会将包含伴随着陀螺传感器的漂移而产生的误差的值用于计算,所以能够检测出正确的姿势角。
而且,不是将加速度传感器的输出值直接用于计算,而是通过去除与移动体的运动有关的加速度而求得重力加速度,所以即使移动体移动也能够正确地检测姿势角。
上面,对本发明的优选的实施例进行了说明,但本发明并不限定于上述实施例。例如,如果移动体是具有车轮的汽车等,则作为移动体的速度测定车轮速度,在要求取纵倾角(Pitch)时,通过对移动成分加速度的上下方向的分力的求取方法追加根据纵倾角和速度计算的移动体的加速度,能够检测更高精度的姿势角。另外,如果移动体是不具有车轮的飞机等,则可以考虑测定螺旋桨的旋转速度作为移动体的速度。
本发明的姿势角检测装置与姿势角检测方法能够使用于工业、例如汽车工业、飞机工业等领域。

Claims (8)

1.一种姿势角检测装置,其特征在于,包括:
加速度传感器,该加速度传感器测定施加于移动体的加速度;
横摆率传感器,该横摆率传感器测定该移动体的横摆率;
速度传感器,该速度传感器测定该移动体的速度;
移动成分加速度计算单元,该移动成分加速度计算单元从该速度计算实际加速度,从该横摆率和该速度计算离心力,并计算作为该实际加速度与该离心力的合力的移动成分加速度;和
姿势角计算单元,该姿势角计算单元从用该移动成分加速度对该加速度进行补正而得到的重力加速度计算姿势角。
2.如权利要求1所述的姿势角检测装置,其中,所述姿势角计算单元包括:重力加速度计算单元,该重力加速度计算单元计算从所述加速度中除去所述移动成分加速度之后的所述重力加速度;侧倾角计算单元,该侧倾角计算单元从该重力加速度计算该移动体的侧倾角;和纵倾角计算单元,该纵倾角计算单元从该重力加速度计算该移动体的纵倾角。
3.如权利要求1或2所述的姿势角检测装置,其中,所述移动成分加速度计算单元包括:前后轴加速度计算单元,该前后轴加速度计算单元对所述速度进行时间微分而计算所述实际加速度,作为所述移动成分加速度的前后轴方向的分力;左右轴加速度计算单元,该左右轴加速度计算单元从所述横摆率和该速度计算离心力,作为所述移动成分加速度的左右轴方向的分力;和上下轴加速度计算单元,该上下轴加速度计算单元计算该离心力的在该移动体的上下轴方向上的分力,作为所述移动成分加速度的上下轴方向的分力。
4.如权利要求3所述的姿势角检测装置,其中:
所述上下轴加速度计算单元进行的所述离心力的分力的计算基于由所述姿势角计算单元计算的所述侧倾角或者所述纵倾角进行;
该姿势角计算单元基于计算出的该离心力的分力计算该侧倾角以及该纵倾角中的至少一方。
5.一种姿势角检测方法,其特征在于,具有:
测定步骤,该测定步骤测定施加于移动体的加速度、该移动体的横摆率和该移动体的速度;
移动成分加速度计算步骤,该移动成分加速度计算步骤从该速度计算实际加速度,从该横摆率和该速度计算离心力,并计算作为该实际加速度与该离心力的合力的移动成分加速度;和
姿势角计算步骤,该姿势角计算步骤从用该移动成分加速度对该加速度进行补正而得到的重力加速度计算姿势角。
6.如权利要求5所述的姿势角检测方法,其中,所述姿势角计算步骤包括:重力加速度计算步骤,该重力加速度计算步骤计算从所述加速度中除去所述移动成分加速度之后的所述重力加速度;侧倾角计算步骤,该侧倾角计算步骤从该重力加速度计算该移动体的侧倾角;和纵倾角计算步骤,该纵倾角计算步骤从该重力加速度计算该移动体的纵倾角。
7.如权利要求5或6所述的姿势角检测方法,其中,所述移动成分加速度计算步骤包括:前后轴加速度计算步骤,该前后轴加速度计算步骤对所述速度进行时间微分而计算所述实际加速度,作为所述移动成分加速度的前后轴方向的分力;左右轴加速度计算步骤,该左右轴加速度计算步骤从所述横摆率和该速度计算离心力,作为所述移动成分加速度的左右轴方向的分力;和上下轴加速度计算步骤,该上下轴加速度计算步骤计算该离心力的在该移动体的上下轴方向上的分力,作为所述移动成分加速度的上下轴方向的分力。
8.如权利要求7所述的姿势角检测方法,其中:
所述重力加速度计算步骤是这样的步骤:将所述重力加速度分为前后轴方向、左右轴方向以及上下轴方向的分力进行计算,基于在所述侧倾角计算步骤以及所述纵倾角计算步骤中的一方中计算出的所述侧倾角或者所述纵倾角计算各分力中的一部分;
该侧倾角计算步骤以及该纵倾角计算步骤中的一方基于各分力中的残余部分计算该侧倾角或者该纵倾角,该侧倾角计算步骤以及该纵倾角计算步骤中的另一方基于各分力中的一部分乃至全部计算该侧倾角或者该纵倾角。
CN200780023332XA 2006-06-23 2007-06-25 姿势角检测装置与姿势角检测方法 Expired - Fee Related CN101473193B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP173964/2006 2006-06-23
JP2006173964A JP4321554B2 (ja) 2006-06-23 2006-06-23 姿勢角検出装置と姿勢角検出方法
PCT/JP2007/062707 WO2007148818A1 (ja) 2006-06-23 2007-06-25 姿勢角検出装置と姿勢角検出方法

Publications (2)

Publication Number Publication Date
CN101473193A true CN101473193A (zh) 2009-07-01
CN101473193B CN101473193B (zh) 2012-05-09

Family

ID=38833551

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200780023332XA Expired - Fee Related CN101473193B (zh) 2006-06-23 2007-06-25 姿势角检测装置与姿势角检测方法

Country Status (6)

Country Link
US (1) US8200452B2 (zh)
EP (1) EP2034270B1 (zh)
JP (1) JP4321554B2 (zh)
KR (1) KR101049362B1 (zh)
CN (1) CN101473193B (zh)
WO (1) WO2007148818A1 (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102506822A (zh) * 2011-11-30 2012-06-20 北京天地玛珂电液控制系统有限公司 一种矿用倾角传感器
CN102648115A (zh) * 2009-11-25 2012-08-22 康蒂特米克微电子有限公司 用于在行驶的汽车中估算侧倾角的方法
CN102730000A (zh) * 2011-03-31 2012-10-17 比亚迪股份有限公司 车辆动态质心的计算方法、横摆力矩的计算方法和系统
CN103171561A (zh) * 2013-03-25 2013-06-26 广州市雄兵汽车电器有限公司 汽车姿态检测方法
CN105035091A (zh) * 2015-08-10 2015-11-11 北京经纬恒润科技有限公司 一种车身姿势变化的检测方法及装置
CN105799708A (zh) * 2014-09-22 2016-07-27 现代自动车株式会社 通过使用重力加速度传感器估计道路坡度的设备和方法
CN108627155A (zh) * 2018-03-30 2018-10-09 华南农业大学 一种农业机械非直线作业向心加速度估算及倾角融合方法
CN110192001A (zh) * 2016-11-09 2019-08-30 萨兰达有限公司 用于可旋转元件的传感器
CN110832342A (zh) * 2017-06-30 2020-02-21 古野电气株式会社 航迹预测装置、航迹预测方法及航迹预测程序
CN111886481A (zh) * 2018-05-30 2020-11-03 松下知识产权经营株式会社 位置调整装置以及具备其的投影系统

Families Citing this family (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5039349B2 (ja) * 2006-09-27 2012-10-03 ヤマハ発動機株式会社 鞍乗型車両の姿勢検知装置、エンジン制御装置、及び鞍乗型車両
JP4946562B2 (ja) * 2007-03-23 2012-06-06 日本電気株式会社 姿勢検出装置および方法
US9157737B2 (en) * 2008-06-11 2015-10-13 Trimble Navigation Limited Altimeter with calibration
US7856336B2 (en) 2008-06-11 2010-12-21 Trimble Navigation Limited Forward-looking altitude detector
US8150651B2 (en) * 2008-06-11 2012-04-03 Trimble Navigation Limited Acceleration compensated inclinometer
US8566034B1 (en) * 2008-06-11 2013-10-22 Trimble Navigation Limited Inclinometer with calibration
EP2140915B1 (en) 2008-06-30 2019-03-06 Nintendo Co., Ltd. Orientation calculation apparatus, storage medium having orientation calculation program stored therein, game apparatus, and storage medium having game program stored therein
EP2140916B1 (en) 2008-06-30 2018-10-31 Nintendo Co., Ltd. Coordinate calculation apparatus and storage medium having coordinate calculation program stored therein
JP5455191B2 (ja) * 2009-04-08 2014-03-26 任天堂株式会社 姿勢算出装置、姿勢算出プログラム、姿勢算出システム、および姿勢算出方法
US8209065B2 (en) * 2009-01-21 2012-06-26 The Boeing Company Steering logic for spacecraft slew maneuvers
EP2228109B1 (en) 2009-03-09 2021-03-24 Nintendo Co., Ltd. Information processing apparatus, storage medium having information processing program stored therein, information processing system, and display range control method
JP5438382B2 (ja) * 2009-05-27 2014-03-12 任天堂株式会社 情報処理プログラムおよび情報処理装置
JP5352922B2 (ja) * 2009-10-07 2013-11-27 多摩川精機株式会社 角速度算出装置、ナビゲーション装置
JP5393432B2 (ja) 2009-12-21 2014-01-22 ヤマハ発動機株式会社 ロール角推定装置および輸送機器
KR101240578B1 (ko) * 2010-03-31 2013-03-08 (주)블루포인트 차량의 경사각 측정방법 및 그 측정장치
JP2012141219A (ja) * 2010-12-28 2012-07-26 Pioneer Electronic Corp 傾斜角検出装置、方法、プログラムおよび記録媒体
JP6094026B2 (ja) * 2011-03-02 2017-03-15 セイコーエプソン株式会社 姿勢判定方法、位置算出方法及び姿勢判定装置
JP5931198B2 (ja) 2012-08-01 2016-06-08 三菱電機株式会社 傾斜角検出装置および車載機器
JP5996572B2 (ja) * 2014-03-27 2016-09-21 本田技研工業株式会社 車体のロール角推定装置
DE102015002601A1 (de) * 2015-02-28 2016-09-01 Audi Ag Verfahren zum Erfassen einer Neigung einer Fahrbahn
JP6454857B2 (ja) * 2015-05-27 2019-01-23 多摩川精機株式会社 姿勢検出装置及び姿勢検出方法
JP6380468B2 (ja) * 2016-06-21 2018-08-29 マツダ株式会社 四輪駆動車の制御装置
KR101910515B1 (ko) * 2017-10-26 2018-10-22 엘아이지넥스원 주식회사 수중 운동체 항법 제어 장치 및 방법
FR3077792B1 (fr) * 2018-02-14 2020-10-02 Aml Systems Procede et dispositif autonomes de determination d’une assiette d’un vehicule automobile.
CN108592878A (zh) * 2018-05-17 2018-09-28 中国人民解放军海军工程大学 海洋和航空重力仪平台静态水平倾角的测量方法
CN109204317B (zh) * 2018-07-24 2020-05-22 吉林大学 轮毂驱动电动汽车纵、横和垂向力集成控制优化方法
KR20230118642A (ko) * 2021-01-18 2023-08-11 아이치 세이코우 가부시키가이샤 제어 방법 및 제어 시스템

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0783659A (ja) * 1993-09-10 1995-03-28 Mitsubishi Electric Corp 車両用傾斜角度検出装置
JPH08178687A (ja) 1994-12-26 1996-07-12 Pioneer Electron Corp 姿勢検出方法及び装置
JP2885125B2 (ja) * 1995-03-30 1999-04-19 トヨタ自動車株式会社 車両の旋回に伴って変化する運動状態量の推定方法
JPH095104A (ja) 1995-06-23 1997-01-10 Nippon Telegr & Teleph Corp <Ntt> 移動物体の三次元姿勢角測定法および三次元姿勢角計測装置
JP2904118B2 (ja) 1996-05-27 1999-06-14 日本電気株式会社 姿勢角検出装置及び方法
JP3576789B2 (ja) 1998-02-09 2004-10-13 株式会社日本自動車部品総合研究所 道路形状計測装置
US6263261B1 (en) * 1999-12-21 2001-07-17 Ford Global Technologies, Inc. Roll over stability control for an automotive vehicle
US6332104B1 (en) 1999-12-21 2001-12-18 Ford Global Technologies, Inc. Roll over detection for an automotive vehicle
JP2001356014A (ja) * 2000-06-13 2001-12-26 Mitsubishi Motors Corp 路面勾配検出装置
JP2003139536A (ja) 2001-11-05 2003-05-14 Japan Aviation Electronics Industry Ltd 方位計および方位測定方法
WO2006028013A1 (ja) * 2004-09-07 2006-03-16 Vodafone K.K. 加速度推定方法、加速度推定装置及び移動端末装置
CN100400331C (zh) * 2004-11-17 2008-07-09 丰田自动车株式会社 车辆以及车辆的控制方法
JP2007040763A (ja) * 2005-08-01 2007-02-15 Toyota Motor Corp 加速度センサの補正装置

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102648115B (zh) * 2009-11-25 2015-11-25 康蒂特米克微电子有限公司 用于在行驶的汽车中估算侧倾角的方法
CN102648115A (zh) * 2009-11-25 2012-08-22 康蒂特米克微电子有限公司 用于在行驶的汽车中估算侧倾角的方法
CN102730000A (zh) * 2011-03-31 2012-10-17 比亚迪股份有限公司 车辆动态质心的计算方法、横摆力矩的计算方法和系统
CN102730000B (zh) * 2011-03-31 2016-06-01 比亚迪股份有限公司 车辆动态质心的计算方法、横摆力矩的计算方法和系统
CN102506822A (zh) * 2011-11-30 2012-06-20 北京天地玛珂电液控制系统有限公司 一种矿用倾角传感器
CN102506822B (zh) * 2011-11-30 2013-10-23 北京天地玛珂电液控制系统有限公司 一种矿用倾角传感器
CN103171561A (zh) * 2013-03-25 2013-06-26 广州市雄兵汽车电器有限公司 汽车姿态检测方法
CN105799708A (zh) * 2014-09-22 2016-07-27 现代自动车株式会社 通过使用重力加速度传感器估计道路坡度的设备和方法
CN105799708B (zh) * 2014-09-22 2019-05-21 现代自动车株式会社 通过使用重力加速度传感器估计道路坡度的设备和方法
CN105035091A (zh) * 2015-08-10 2015-11-11 北京经纬恒润科技有限公司 一种车身姿势变化的检测方法及装置
CN110192001A (zh) * 2016-11-09 2019-08-30 萨兰达有限公司 用于可旋转元件的传感器
US11402205B2 (en) 2016-11-09 2022-08-02 Salunda Limited Sensor for a rotatable element
CN110832342A (zh) * 2017-06-30 2020-02-21 古野电气株式会社 航迹预测装置、航迹预测方法及航迹预测程序
CN108627155A (zh) * 2018-03-30 2018-10-09 华南农业大学 一种农业机械非直线作业向心加速度估算及倾角融合方法
CN111886481A (zh) * 2018-05-30 2020-11-03 松下知识产权经营株式会社 位置调整装置以及具备其的投影系统

Also Published As

Publication number Publication date
WO2007148818A1 (ja) 2007-12-27
EP2034270A4 (en) 2012-11-28
KR20090018659A (ko) 2009-02-20
EP2034270A1 (en) 2009-03-11
JP4321554B2 (ja) 2009-08-26
US20090326858A1 (en) 2009-12-31
KR101049362B1 (ko) 2011-07-13
JP2008002992A (ja) 2008-01-10
US8200452B2 (en) 2012-06-12
EP2034270B1 (en) 2014-03-19
CN101473193B (zh) 2012-05-09

Similar Documents

Publication Publication Date Title
CN101473193B (zh) 姿势角检测装置与姿势角检测方法
Ahmed et al. Accurate attitude estimation of a moving land vehicle using low-cost MEMS IMU sensors
CN107655493B (zh) 一种光纤陀螺sins六位置系统级标定方法
CN106814753B (zh) 一种目标位置矫正方法、装置及系统
CN102066875A (zh) 测斜仪测量系统和提供用于运动诱导加速度误差的校正的方法
CN108132053B (zh) 一种行人轨迹构建方法、系统及惯性测量装置
JP2011185899A (ja) 位置標定装置、位置標定装置の位置標定方法および位置標定プログラム
CN114061619B (zh) 一种基于在线标定的惯导系统姿态补偿方法
CN106153069A (zh) 自主导航系统中的姿态修正装置和方法
JP2016017796A (ja) 車輌位置計測装置及び方法
CN110426011A (zh) 车辆转向角测量系统及方法
Promrit et al. In-vehicle MEMS IMU calibration using accelerometer
CN109186634B (zh) 一种mems惯组导航性能测量方法及装置
CN111902693A (zh) 用于校准装备于车辆的陀螺仪的方法
CN109827572A (zh) 一种检测车位置预测的方法及装置
CN114964224B (zh) 一种微惯导系统误差自主抑制方法
CN112577512A (zh) 一种基于轮速融合的状态量误差确定方法及车载终端
JP5591185B2 (ja) 車両の位置姿勢測定装置
JP2009053039A (ja) 車両姿勢推定装置、車両姿勢推定方法
CN106931992B (zh) 用于检测物体翻滚的方法和装置
CN105222764B (zh) 一种对惯性角速度传感器进行地球自转补偿的方法
CN113985466A (zh) 一种基于模式识别的组合导航方法及系统
CN113959464A (zh) 一种陀螺仪辅助的加速度计现场校准方法和系统
JP6409625B2 (ja) 車両位置算出装置
JP6454857B2 (ja) 姿勢検出装置及び姿勢検出方法

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120509

CF01 Termination of patent right due to non-payment of annual fee