[go: up one dir, main page]

WO2010055871A1 - 姿勢検出装置の補正パラメーター作成方法、姿勢検出装置の補正パラメーター作成用装置及び姿勢検出装置 - Google Patents

姿勢検出装置の補正パラメーター作成方法、姿勢検出装置の補正パラメーター作成用装置及び姿勢検出装置 Download PDF

Info

Publication number
WO2010055871A1
WO2010055871A1 PCT/JP2009/069249 JP2009069249W WO2010055871A1 WO 2010055871 A1 WO2010055871 A1 WO 2010055871A1 JP 2009069249 W JP2009069249 W JP 2009069249W WO 2010055871 A1 WO2010055871 A1 WO 2010055871A1
Authority
WO
WIPO (PCT)
Prior art keywords
axis
sensor
detection
correction
angular velocity
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.)
Ceased
Application number
PCT/JP2009/069249
Other languages
English (en)
French (fr)
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.)
Miyazaki Epson Corp
Original Assignee
Epson Toyocom 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 Epson Toyocom Corp filed Critical Epson Toyocom Corp
Priority to US13/126,446 priority Critical patent/US20110202300A1/en
Priority to CN2009801453275A priority patent/CN102216790B/zh
Publication of WO2010055871A1 publication Critical patent/WO2010055871A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/02Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
    • G01B21/04Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness by measuring coordinates of points
    • G01B21/045Correction of measurements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C19/00Gyroscopes; Turn-sensitive devices using vibrating masses; Turn-sensitive devices without moving masses; Measuring angular rate using gyroscopic effects
    • G01C19/56Turn-sensitive devices using vibrating masses, e.g. vibratory angular rate sensors based on Coriolis forces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
    • G01C25/005Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P21/00Testing or calibrating of apparatus or devices covered by the preceding groups
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F3/00Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
    • G06F3/01Input arrangements or combined input and output arrangements for interaction between user and computer
    • G06F3/011Arrangements for interaction with the human body, e.g. for user immersion in virtual reality
    • G06F3/012Head tracking input arrangements

Definitions

  • the present invention provides a correction parameter creation method, a correction parameter creation device, and a correction function for correcting a detection value of an attitude detection device including a sensor that detects a three-axis angular velocity or acceleration into a detection value in a predetermined orthogonal coordinate system.
  • the present invention relates to an attitude detection device.
  • Japanese Patent Application Laid-Open No. 9-106322 discloses a head mount that detects the posture of the user's head so that the image displayed on the front display changes in conjunction with the movement of the head and allows the user to experience a virtual space. A display is described. An image that matches the posture angle of the user's head is displayed on the head-mounted display. In order to detect the attitude angle, an attitude detection device including an angular velocity sensor and an acceleration sensor is attached to a predetermined position of the head mounted display.
  • the posture detection device When the posture detection device is installed, if the sensor detection axis is not mounted so as to be parallel to the three axes of the Cartesian coordinate system for representing the posture angle of the head, posture detection is caused by this mounting angle error.
  • the detected value of the apparatus includes an error. Therefore, the position and angle at which the posture detection device is attached are strictly defined.
  • Equations (1) and (2) represent an angular velocity sensor correction equation and an acceleration sensor correction equation using correction parameters, respectively.
  • the function determinant (Jacobiane) J fG is a correction parameter for the angular velocity sensor, and f G (x) and f G (p) are the current correction value and the previous correction value (ideal value) of the angular velocity sensor, respectively. It is.
  • the function determinant (Jacobiane) J fA is a correction parameter for the acceleration sensor, and f A (x) and f A (p) are the correction values (according to the current and previous values of the acceleration sensor, respectively) Ideal value).
  • x and p are the current and previous detected values of the angular velocity sensor or the angular velocity sensor, respectively
  • o is a Landau symbol.
  • FIGS. 13A to 13C and FIGS. 14A to 14C show a conventional method for creating correction parameters (J fG , J fA ).
  • the attitude detection device is set in the socket 520 attached to the table 510, and the rotary arm 530 is rotated around the X axis, the Y axis, and the Z axis in the order shown in FIGS. 13A to 13C.
  • each detection value of the attitude detection device solves the simultaneous equations obtained by substituting each detection value and each ideal value into the equation (1), and create a correction parameter for the angular velocity sensor. Further, in the order shown in FIGS. 14A to 14C, the rotary arm 530 is operated to stand still in a state in which the positive directions of the X axis, the Y axis, and the Z axis are vertically upward (gravity acceleration is applied vertically downward). Each detection value of the detection device is acquired, and a correction parameter for the acceleration sensor is created by solving simultaneous equations obtained by substituting each detection value and each ideal value into Equation (2).
  • the present invention has been made in view of the above problems, and according to some aspects of the present invention, a correction parameter for correcting an error in a detected value caused by a mounting angle error of a sensor is provided. It is possible to provide a correction parameter creation method of a posture detection device that can be created at a lower cost, a correction parameter creation device that can be realized at a lower cost, and a posture detection device with a correction function.
  • the present invention provides a first sensor for detecting angular velocity or acceleration, which is mounted so that detection axes are substantially parallel to a first axis, a second axis, and a third axis that are orthogonal to each other, And a detection value of an attitude detection device that detects an attitude of an object based on detection signals of the first sensor, the second sensor, and the third sensor.
  • the second axis, and the third axis as a coordinate axis, and a correction parameter of a correction formula for correcting to a detected value in a rectangular coordinate system, wherein the rotating plate is mounted so that the upper surface is horizontal.
  • the first axis of the rectangular parallelepiped jig having a first surface, a second surface, and a third surface orthogonal to each other, and the first axis is perpendicular to the second surface And the second axis is perpendicular to the third plane
  • the correspondence between the first, second, and third axes and the X, Y, and Z axes is not particularly limited. .
  • the first axis, the second axis, and the third axis are perpendicular to the first surface, the second surface, and the third surface of the jig, respectively. It is easy to fix the posture detection device to the first surface. And if the rotating plate is installed so that the upper surface is horizontal, it is easy to fix the second surface, the third surface, and the surface facing the first surface of the jig to the upper surface of the rotating plate.
  • the first axis, the second axis, and the third axis can be parallel to the vertical direction.
  • the detection value of the acceleration sensor or the angular velocity sensor can be easily reduced in a short time by stationary or rotating the rotating plate. Can be obtained at.
  • the rotating plate is installed so that the upper surface is horizontal only once at the beginning, the rotation direction of the rotating plate is fixed, so the detection values for the first axis, the second axis, and the third axis are acquired. Setting time can be greatly shortened. Therefore, according to the present invention, it is possible to create a correction parameter for correcting an error in a detected value caused by a sensor mounting angle error at a lower cost.
  • the correction formula may use the orthogonal coordinates for the detection values of the first sensor, the second sensor, and the third sensor as the correction parameter. Including a first correction matrix, a second correction matrix, and a third correction matrix for correcting each detected value in the system, the first correction matrix, the second correction matrix, and the third correction matrix 3 obtained by multiplying each detected value of each of the first sensor, the second sensor, and the third sensor with a digital value obtained by A / D conversion as an element, respectively. It may be given as the sum of two matrices.
  • the first correction matrix, the second correction matrix, and the third correction matrix are the first sensor, the second sensor, and
  • the detection axis of the third sensor may be an inverse matrix of a rotation matrix that converts the detection axis into the first axis, the second axis, and the third axis, respectively.
  • the correction parameter creation step includes the second sensor and the third sensor based on the detection value acquired in the first detection value acquisition step. Based on the detection value acquired in the step of calculating each mounting angle error of the sensor about the first axis and the second detection value acquisition step, the first sensor and the third sensor of the third sensor. And calculating each mounting angle error about the second axis, and the third axis of the first sensor and the second sensor based on the detection value acquired in the third detection value acquisition step. And calculating the mounting angle error of the first sensor around the second axis and the mounting angle error around the third axis of the first sensor.
  • the present invention provides a first sensor that detects angular velocity or acceleration, and is attached so that detection axes are substantially parallel to a first axis, a second axis, and a third axis that are orthogonal to each other. And a detection value of an attitude detection device that detects an attitude of an object based on detection signals of the first sensor, the second sensor, and the third sensor.
  • a correction parameter creation device used to create a correction parameter of a correction formula for correcting to a detected value in a Cartesian coordinate system with the second axis and the third axis as coordinate axes.
  • a rectangular parallelepiped-shaped jig capable of fixing the posture detection device so as to be straight, and on the upper surface, the first surface, the second surface, and the second surface of the jig, respectively. It includes a rotating plate that can fix any one of them, and a rotation control unit that rotates the rotating plate at a predetermined angular velocity.
  • the correction parameter creation device since a rotating arm is not required by using a rectangular parallelepiped jig and a rotating plate, it is possible to provide a more compact and low-cost device for creating correction parameters.
  • the correction parameter creation device as described above, the correction parameter of the detection value of each sensor attached to the posture detection device can be easily acquired in a short time.
  • the present invention provides a first sensor that detects angular velocity or acceleration, and is attached so that detection axes are substantially parallel to a first axis, a second axis, and a third axis that are orthogonal to each other.
  • the first sensor, the third sensor, and the detected values of the first sensor, the second sensor, and the third sensor are coordinate axes of the first axis, the second axis, and the third axis.
  • a storage unit storing correction parameters of correction formulas for correcting detection values in the Cartesian coordinate system, and detecting signals of the first sensor, the second sensor, and the third sensor are converted into digital signals.
  • the first correction matrix, the second correction matrix, and the second correction matrix for correcting the detection values of the first sensor, the second sensor, and the third sensor to the detection values in the orthogonal coordinate system, respectively.
  • the detection value is given as the sum of three matrices obtained by the product of each of the matrices each including an A / D converted digital value as an element.
  • the function determinant (Jacobiane) in the conventional correction formula (1) and correction formula (2) is not a correction parameter that directly reflects the sensor mounting angle error, and the correction formula (1) and correction formula (2).
  • a correction value cannot be obtained if any mapping is applied to the detection value. Therefore, there is a limit to increasing the correction accuracy in the correction equations (1) and (2).
  • the present invention it is possible to directly reflect the mounting angle error of each sensor in the three correction matrices included in the correction formula calculated by the correction calculation processing unit. Further, according to the present invention, the correction formula calculated by the correction calculation processing unit does not require the previous detection value in calculating the correction value for the current detection value. Can be calculated. Therefore, according to the present invention, it is possible to realize an attitude detection device with higher correction accuracy and faster correction calculation processing.
  • the first correction matrix, the second correction matrix, and the third correction matrix are the first sensor, the second sensor, and the third sensor.
  • Each detection axis may be an inverse matrix of a rotation matrix that converts the detection axis into the first axis, the second axis, and the third axis, respectively.
  • the posture detection apparatus performs a process of sequentially selecting any one of the detection signals of the first sensor, the second sensor, and the third sensor at a predetermined cycle.
  • a signal selection processing unit may be included, and the A / D conversion processing unit may include an A / D conversion circuit that sequentially performs A / D conversion processing on the detection values selected by the signal selection processing unit.
  • FIG. 1 is a diagram illustrating an example of a configuration of a posture detection apparatus that is a target of a correction parameter creation method according to the present embodiment.
  • FIG. 2 is a perspective view of the posture detection apparatus in the present embodiment.
  • FIG. 3 is a plan view showing an example of a vibrator included in the angular velocity sensor.
  • FIG. 4 is a diagram for explaining the operation of the vibrator included in the angular velocity sensor.
  • FIG. 5 is a diagram for explaining the operation of the vibrator included in the angular velocity sensor.
  • FIG. 6 is a diagram illustrating an example of a configuration of a drive circuit and a detection circuit included in the angular velocity sensor.
  • FIG. 7A is a diagram for explaining a mounting angle error of a sensor.
  • FIG. 7B is a diagram for explaining a sensor mounting angle error.
  • FIG. 7C is a diagram for explaining a sensor mounting angle error.
  • FIG. 8 is a diagram showing the configuration of the correction parameter creation apparatus of the present embodiment.
  • FIG. 9 is a flowchart showing an example of a correction parameter creation procedure in the present embodiment.
  • FIG. 10A is a diagram for describing a correction parameter creation procedure in the present embodiment.
  • FIG. 10B is a diagram for describing a correction parameter creation procedure in the present embodiment.
  • FIG. 10C is a diagram for describing a correction parameter creation procedure in the present embodiment.
  • FIG. 11 is a diagram illustrating a configuration of the posture detection apparatus of the present embodiment.
  • FIG. 12 is a diagram illustrating another configuration of the posture detection apparatus according to the present embodiment.
  • FIG. 13A is a diagram for explaining a conventional correction parameter creation method.
  • FIG. 13B is a diagram for explaining a conventional correction parameter creation method.
  • FIG. 13C is a diagram for explaining a conventional correction parameter creation method.
  • FIG. 14A is a diagram for explaining a conventional correction parameter creation method.
  • FIG. 14B is a diagram for explaining a conventional correction parameter creation method.
  • FIG. 14C is a diagram for explaining a conventional correction parameter creation method.
  • the first axis, the second axis, and the third axis in the present invention correspond to the X axis, the Y axis, and the Z axis, respectively, but the first axis in the present invention.
  • the correspondence relationship between the second axis, the third axis, the X axis, the Y axis, and the Z axis is not limited to this, and may be an arbitrary correspondence relationship.
  • FIG. 1 is a diagram illustrating an example of the configuration of an attitude detection device that is a target of the correction parameter creation method of the present embodiment.
  • the posture detection apparatus 1 in this embodiment detects an angular velocity sensor module 2 that detects angular velocities around the X axis, Y axis, and Z axis, and detects acceleration in the X axis, Y axis, and Z axis directions.
  • the acceleration sensor module 3 is configured to be included.
  • the angular velocity sensor module 2 includes an X-axis angular velocity sensor 10a, a Y-axis angular velocity sensor 10b, and a Z-axis angular velocity sensor 10c that detect angular velocities around the X, Y, and Z axes, respectively.
  • the X-axis angular velocity sensor 10a includes a vibrator 11a, a drive circuit 20a that vibrates the vibrator 11a, and a detection circuit 30a that generates an angular velocity detection signal 38a.
  • the detection electrodes 14a and 15a of 11a and the detection circuit 30a are connected to each other.
  • the Y-axis angular velocity sensor 10b includes a vibrator 11b, a drive circuit 20b that vibrates the vibrator 11b, and a detection circuit 30b that generates an angular velocity detection signal 38b, and the drive electrodes 12b and 13b of the vibrator 11b and the drive circuit 20b.
  • the detection electrodes 14b and 15b of the vibrator 11b and the detection circuit 30b are connected to each other.
  • the Z-axis angular velocity sensor 10c includes a vibrator 11c, a drive circuit 20c that drives the vibrator 11c to vibrate, and a detection circuit 30c that generates an angular velocity detection signal 38c, and drives the drive electrodes 12c and 13c of the vibrator 11c.
  • the detection electrodes 14c and 15c of the circuit 20c and the vibrator 11c are connected to the detection circuit 30c.
  • the acceleration sensor module 3 includes an X-axis acceleration sensor 50a, a Y-axis acceleration sensor 50b, and a Z-axis acceleration sensor 50c that detect accelerations in the X-axis, Y-axis, and Z-axis directions, respectively.
  • the X-axis acceleration sensor 50a includes a vibrator 51a, a drive circuit 60a that vibrates the vibrator 51a, and a detection circuit 70a that generates an acceleration detection signal 78a.
  • the detection electrodes 54a and 55a of 51a and the detection circuit 70a are connected to each other.
  • the Y-axis acceleration sensor 50b includes a vibrator 51b, a drive circuit 60b that vibrates the vibrator 51b, and a detection circuit 70b that generates an acceleration detection signal 78b.
  • the detection electrodes 54b and 55b of the vibrator 51b and the detection circuit 70b are connected to each other.
  • the Z-axis acceleration sensor 50c includes a vibrator 51c, a drive circuit 60c that vibrates the vibrator 51c, and a detection circuit 70c that generates an acceleration detection signal 78c.
  • the detection electrodes 54c and 55c of the vibrator 51c and the detection circuit 70c are connected to each other.
  • the angular velocity sensors 10a, 10b, and 10c function as the first sensor, the second sensor, and the third sensor in the present invention, respectively.
  • the acceleration sensors 50a, 50b, and 50c function as the first sensor, the second sensor, and the third sensor in the present invention, respectively.
  • FIG. 2 is a perspective view of the posture detection apparatus in the present embodiment.
  • the angular velocity sensor module 2 and the acceleration sensor module 3 are each formed in a cubic shape (in the broad sense, a rectangular parallelepiped, the same applies hereinafter), and the rectangular parallelepiped package 4. Is housed inside.
  • the X axis, Y axis, and Z axis are determined based on the attitude detection device 1.
  • the axes perpendicular to the three orthogonal surfaces 5a, 5b, and 5c of the package 4 may be the X axis, the Y axis, and the Z axis, respectively. it can.
  • the positive directions of the X axis, the Y axis, and the Z axis can be arbitrarily determined. In this embodiment, the direction toward the tip of the arrow shown in FIG.
  • Angular Velocity Sensor Module As shown in FIG. 2, in the angular velocity sensor module 2, the angular velocity sensors 10a, 10b, and 10c are placed on the insulating substrate 80 so that the detection axes are substantially parallel to the X axis, the Y axis, and the Z axis, respectively.
  • the vibrators 11a, 11b, and 11c are mounted in the packages 82a, 82b, and 82c, respectively.
  • the periphery of the packages 82a, 82b, and 82c is covered with a resin mold material.
  • the packages 82a, 82b, and 82c include a package body 84a and a lid body 86a, a package body 84b and a lid body 86b, and a package body 84c and a lid body 86c, respectively.
  • the package bodies 84a, 84b, 84c are formed in a rectangular parallelepiped box shape by laminating and sintering a plurality of ceramic sheets.
  • the lids 86a, 86b, and 86c are formed of a glass plate, a metal plate, a ceramic sheet, and the like, and accommodate the vibrators 11a, 11b, and 11c through a bonding material such as a metal brazing material and low-melting glass.
  • the upper surface openings of the package bodies 84a, 84b, 84c are vacuum sealed.
  • the vibrators 11a, 11b, and 11c are connected to the drive circuits 20a, 20b, and 20c and the detection circuits 30a, 30b, and 30c by wiring patterns (not shown) formed on the insulating substrate 80, respectively.
  • the driving circuit 20a and the detection circuit 30a, the driving circuit 20b and the detection circuit 30b, and the driving circuit 20c and the detection circuit 30c may be integrated into three chips and accommodated in the packages 82a, 82b, and 82c, respectively. Further, the drive circuits 20a, 20b, 20c and the detection circuits 30a, 30b, 30c may be integrated on a single chip and disposed on the insulating substrate 80.
  • the detection signals 38a, 38b, and 38c from the detection circuits 30a, 30b, and 30c are output to the outside of the attitude detection device 1 via external output terminals (not shown). It has come to be.
  • FIG. 3 is a plan view showing an example of a vibrator included in the angular velocity sensor. Since all the vibrators 11a, 11b, and 11c included in the angular velocity sensors 10a, 10b, and 10c have the same structure, only the structure of the vibrator 11a is illustrated in FIG. Note that the X axis, Y axis, and Z axis in FIG. 3 indicate crystal axes, and are independent of the X axis, Y axis, and Z axis in FIG.
  • the vibrator 11a is formed of a thin plate made of a piezoelectric material such as quartz, and a driving vibrating arm 41a (driving vibrating piece in a broad sense) extends from the driving base 44a in the Y-axis direction of the crystal.
  • Drive electrodes 12a and 13a are formed on the side and top surfaces of the drive vibrating arm 41a, respectively. As shown in FIG. 1, the drive electrodes 12a and 13a are connected to the drive circuit 20a.
  • the drive base 44a is connected to the detection base 47a via a connecting arm 45a extending in the X-axis direction of the crystal.
  • the detection vibration arm 42a (detection vibration piece in a broad sense) extends from the detection base 47a in the Y-axis direction of the crystal.
  • Detection electrodes 14a and 15a are formed on the upper surface of the detection vibrating arm 42a, and an electrode 16a is formed on the side surface of the detection vibrating arm 42a. As shown in FIG. 1, the detection electrodes 14a and 15a are connected to the drive circuit 20a.
  • the electrode 16a is grounded.
  • the drive vibration arm 41a When a drive signal comprising an alternating voltage / alternating current is applied between the drive electrode 12a and the drive electrode 13a of the drive vibration arm 41a, the drive vibration arm 41a is as shown by an arrow B by the piezoelectric effect as shown in FIG. Bends and vibrates.
  • the drive vibration arm 41a is perpendicular to both the direction of the bending vibration of the arrow B and the crystal Z axis. Get Coriolis power in the direction.
  • the connecting arm 45a vibrates as indicated by an arrow C.
  • the detection vibrating arm 42a bends and vibrates like the arrow D with the connecting arm 45a in conjunction with the vibration (arrow C) of the connecting arm 45a.
  • the detection base 47a is arranged in the center, and the detection vibrating arm 42a extends from the detection base 47a in both the + Y axis and ⁇ Y axis directions. I am letting. Further, the connecting arm 45a extends from the detection base 47a in both the + X-axis and ⁇ X-axis directions, and the drive vibrating arm 41a extends from each of the connecting arms 45a in both the + Y-axis and ⁇ Y-axis directions. .
  • the Coriolis force is increased by making the distal end of the drive vibrating arm 41a a wide wide portion 43a and attaching a weight. Also, due to the weight effect, a desired resonance frequency can be obtained with a short vibrating arm. For the same reason, the distal end of the detection vibrating arm 42a is a wide portion 46a, and a weight is attached.
  • the vibrator 11a is not limited to the above-described configuration, and may be any vibrator that outputs a detection signal including an angular velocity component based on Coriolis force.
  • it may be configured to serve as both the drive vibration arm and the detection vibration arm, or may be configured such that a piezoelectric film is formed on the drive vibration arm and the detection vibration arm.
  • FIG. 6 is a diagram showing an example of the configuration of a drive circuit and a detection circuit included in the angular velocity sensor. Since the drive circuits 20a, 20b, and 20c all have the same configuration, and the detection circuits 30a, 30b, and 30c all have the same configuration, FIG. 3 illustrates only the configuration of the drive circuit 20a and the detection circuit 30a.
  • the drive circuit 20a includes a current-voltage converter (I / V converter) 21a, an AC amplifier 22a, an automatic gain control circuit (AGC) 23a, and a comparator 24a.
  • I / V converter current-voltage converter
  • AGC automatic gain control circuit
  • the current-voltage converter (I / V converter) 21a converts the input AC current into an AC voltage signal having the same frequency as the vibration frequency of the vibrator 11a and outputs the AC voltage signal.
  • the AC voltage signal output from the current-voltage converter (I / V converter) 21a is input to the AC amplifier 22a.
  • the AC amplifier a amplifies and outputs the input AC voltage signal.
  • the AC voltage signal output from the AC amplifier 22a is input to the automatic gain control circuit (AGC) 23a.
  • the automatic gain control circuit (AGC) 23a controls the gain so as to maintain the amplitude of the input AC voltage signal at a constant value, and outputs the AC voltage signal after gain control to the drive electrode 12a of the vibrator 11a.
  • the vibrator 11a vibrates by an AC voltage signal input to the drive electrode 12a.
  • the AC voltage signal amplified by the AC amplifier 22a is input to the comparator 24a, and a square wave voltage signal that switches the output level in accordance with the comparison result between the AC voltage signal and the reference voltage signal with the amplitude center of the AC voltage signal as the reference voltage.
  • a square wave voltage signal that switches the output level in accordance with the comparison result between the AC voltage signal and the reference voltage signal with the amplitude center of the AC voltage signal as the reference voltage.
  • the detection circuit 30a includes charge amplifiers 31a and 32a, a differential amplifier 33a, an AC amplifier 34a, a synchronous detection circuit 35a, a DC amplifier 36a, and an integration circuit (LPF) 37a.
  • LPF integration circuit
  • detection signals (alternating currents) of opposite phases detected by the vibrator 11a are input via the detection electrodes 12a and 13a.
  • the charge amplifiers 31a and 32a convert the input detection signal (alternating current) into an alternating voltage signal centered on the reference voltage.
  • the differential amplifier 33a differentially amplifies the output signal of the charge amplifier 31a and the output signal of the charge amplifier 32a.
  • the output signal of the differential amplifier 33a is further amplified by the AC amplifier 34a.
  • the synchronous detection circuit 35a extracts an angular velocity component by synchronously detecting the output signal of the AC amplifier 34a based on the square wave voltage signal output from the comparator 24a. For example, when the voltage level of the square wave voltage signal is higher than the reference voltage, the synchronous detection circuit 35a outputs the output signal of the AC amplifier 34a as it is, and when the voltage level of the square wave voltage signal is lower than the reference voltage, AC is detected. It can be configured as a switch circuit that inverts the output signal of the amplifier 34a with respect to the reference voltage and outputs the inverted signal.
  • the integration circuit (LPF) 37a generates an angular velocity detection signal 38a by attenuating a high frequency component from the output signal of the DC amplifier 35a and extracting a direct current component, and outputs it to the outside.
  • the acceleration sensor module 3 includes a base 90, a weight 100, and three acceleration sensors 50a, 50b, and 50c.
  • the drive circuits 60a, 60b, and 60c and the detection circuits 70a, 70b, and 70c shown in FIG. 1 are omitted, but the drive circuits 60a, 60b, and 60c are arranged at appropriate positions in the package 4, and the detection circuits 70a,
  • the detection signals 78a, 78b, 78c from 70b, 70c are output to the outside of the attitude detection device 1 via an external output terminal (not shown).
  • the base 90 is formed by making three square wall portions orthogonal to each other so as to form a cubic shape, and has three attachment surfaces 91, 92, 93 orthogonal to each other in the X-axis, Y-axis, and Z-axis directions.
  • the weight 100 is made of a cube having a predetermined mass, and has three joining surfaces 101, 102, and 103 that are orthogonal to each other.
  • the base 90 and the weight 100 are formed using an appropriate material such as an aluminum alloy.
  • the acceleration sensors 50a, 50b, and 50c are configured to include double tuning fork type vibrators 51a, 51b, and 51c formed of a thin plate of a piezoelectric material such as quartz, respectively.
  • the transducers 51a, 51b, and 51c have element base surfaces 91a, 92b, and 56c on one base end portions 56a, 56b, and 56c, respectively, such that the detection axes are substantially parallel to the X, Y, and Z axes, respectively. , 93 and is vertically supported on each wall portion of the base 90.
  • the other base end portions 57a, 57b, and 57c of the vibrators 51a, 51b, and 51c are bonded to the element bonding surfaces 101 to 103 of the weight 100 corresponding to the element mounting surfaces 91, 92, and 93, respectively.
  • the weight 100 is supported in a floating state by the vibrators 51a, 51b, and 51c from the X-axis, Y-axis, and Z-axis directions.
  • the two drive vibrating arms 58a of the vibrator 51a are provided with drive electrodes 52a and 53a on the upper and lower main surfaces and both side surfaces (not shown), and a predetermined alternating current is provided between the drive electrodes 52a and 53a by the drive circuit 60a.
  • a voltage is applied, the two drive vibrating arms 58a bend and vibrate at a predetermined frequency in directions opposite to each other, that is, in directions close to or away from each other.
  • the vibrator 51a When an external force is applied to the acceleration sensor module 3 and the acceleration in the X-axis direction is applied to the weight 100 in a state where the vibrator 51a is vibrated at a predetermined frequency, the vibrator 51a is applied to the vibrator 51a according to the magnitude and direction thereof.
  • a force to compress or pull in the longitudinal direction acts.
  • the frequency of the vibrator 51a changes so as to decrease when a compressing force is applied and to increase when a pulling force is applied. Therefore, the detection circuit 70a detects the amount of change in frequency in the vibrator 51a, and calculates the load acting on the weight 100 from the amount of change in frequency, thereby calculating the magnitude of the acceleration acting on the weight 100 in the X-axis direction.
  • the direction can be calculated.
  • the structures of the vibrators 51b and 51c are the same as the structure of the vibrator 51a, and similarly, the magnitude and direction of acceleration in the Y-axis and Z-axis directions can be calculated.
  • the drive circuits 60a, 60b, and 60c have the same configuration as the drive circuit 20a shown in FIG. 6, and the detection circuits 70a, 70b, and 70c have the same configuration as a known circuit that detects the amount of change in frequency. The description thereof is omitted.
  • correction parameter creation method 2-1 Sensor mounting angle error
  • the angular velocity sensors 10a, 10b, and 10c are ideally mounted so that the detection axes are strictly parallel to the X, Y, and Z axes, respectively.
  • the acceleration sensors 50a, 50b, and 50c are ideally mounted so that the detection axes are strictly parallel to the X axis, the Y axis, and the Z axis, respectively.
  • the X-axis angular velocity sensor 10a is actually arranged so that the detection axis is parallel to the X ′ axis that is a small angle ⁇ 2x about the Y axis and a small angle ⁇ 3x about the Z axis. Is attached.
  • the Y-axis angular velocity sensor 10b is actually parallel to the Y ′ axis whose detection axis is a small angle ⁇ 3y about the Z axis and a small angle ⁇ 1y about the X axis. As shown in FIG.
  • the Z-axis angular velocity sensor 10c is actually configured such that the detection axis has a small angle ⁇ 1z about the X axis and a small angle ⁇ 2z about the Y axis. It is attached to be parallel to That is, the mounting angle error about the Y axis and the mounting angle error about the Z axis of the X-axis angular velocity sensor 10a are ⁇ 2x and ⁇ 3x , respectively.
  • the mounting angle errors are ⁇ 3y and ⁇ 1y , respectively, and the mounting angle error around the X axis and the mounting angle error around the Y axis of the Z-axis angular velocity sensor 10c are ⁇ 3z and ⁇ 1z , respectively.
  • An arbitrary rotation in the three-dimensional Euclidean space can be represented by a combination of products of rotation matrices T 1 , T 2 , and T 3 .
  • rotation matrices T 1 , T 2 , and T 3 For example, rotated about the Z-axis by an angle theta 3, is rotated by an angle theta 2 about the Y-axis, by rotating about the X axis by an angle theta 1, X'Y'Z the XYZ coordinate system '
  • T ⁇ to be converted to the coordinate system is given by equation (4).
  • T ⁇ is referred to as a “transformation matrix”.
  • the detection axes are actually the X ′ axis, respectively, due to mounting angle errors. It is assumed that they are mounted so as to be parallel to the Y ′ axis and the Z ′ axis.
  • the relational expression (5) between the detected values G x ′, G y ′, G z ′ of the angular velocity sensors 10 a, 10 b, 10 c and the ideal values G x , G y , G z by the transformation matrix T ⁇ Is established.
  • the ideal values G x , G y , G z can be calculated from the detected values G x ′, G y ′, G z ′ of the angular velocity sensors 10a, 10b, 10c by the following equation (6).
  • T ⁇ ⁇ 1 can be obtained by any method, the detected values of the angular velocity sensors 10a, 10b, and 10c can be corrected to ideal values using the equation (6).
  • T ⁇ ⁇ 1 is referred to as a “correction matrix”.
  • ⁇ 1 , ⁇ 2 , and ⁇ 3 are directly derived, and the rotation matrices T 1 , T 2 , and T 3 are calculated using Equation (3).
  • the correction matrix T ⁇ ⁇ 1 can be obtained using the inverse matrices T 1 ⁇ 1 , T 2 ⁇ 1 , and T 3 ⁇ 1 .
  • the angular velocity sensors 10a, 10b, and 10c are rotated around the X, Y, and Z axes.
  • Three input conditions that are reflected are selected, and the detected values G x ′, G y ′, G z ′ and ideal values G x , G y , G z of the angular velocity sensors 10a, 10b, 10c corresponding to the input conditions are selected.
  • ⁇ 1 , ⁇ 2 , and ⁇ 3 can be derived by solving three simultaneous equations obtained by substituting. However, since these simultaneous equations are very complicated, ⁇ 1 , ⁇ 2 , and ⁇ 3 cannot be easily derived.
  • the transformation matrix T ⁇ is expressed as the following equation (8).
  • the transformation matrix T ⁇ can be represented by a linear sum of three base matrices J 1 , J 2 , and J 3 .
  • the detected values G x ′, G y ′, G z ′ of the X-axis angular velocity sensor 10a, the Y-axis angular velocity sensor 10b, and the Z-axis angular velocity sensor 10c are respectively set to ideal values G x , G z. It can correct
  • the correction matrices T ⁇ x ⁇ 1 , T ⁇ y ⁇ 1 , T ⁇ z ⁇ 1 in the correction formula (13) are the mounting angle errors ⁇ 2x , ⁇ 3x , ⁇ 1y , ⁇ 3y , ⁇ 1z of the angular velocity sensors 10 a, 10 b, 10 c. , ⁇ 2z is directly reflected. Further, according to the correction formula (13), if the detection values G x ′, G y ′, G z ′ of this time are obtained without requiring the previous detection values, the correction values (ideal values) G x , G are immediately obtained. y and Gz can be calculated. Therefore, according to the correction equation (13), it is possible to improve the correction accuracy and speed up the correction calculation process.
  • the correction matrices T ⁇ x ⁇ 1 , T ⁇ y ⁇ 1 , and T ⁇ z ⁇ 1 correspond to the first correction matrix, the second correction matrix, and the third correction matrix in the present invention, respectively.
  • ⁇ yy ′ can be obtained by multiplying the detection value (angular velocity around the Y ′ axis) of the Y-axis angular velocity sensor 10b by a predetermined time, by substituting ⁇ yy ′ and ⁇ x into equation (17), ⁇ 3y can be obtained.
  • ⁇ zz ′ can be obtained by multiplying the detection value (angular velocity around the Z ′ axis) of the Z-axis angular velocity sensor 10c by a predetermined time, by substituting ⁇ zz ′ and ⁇ x into equation (18), ⁇ 2z can be obtained.
  • ⁇ xx ′ can be obtained by multiplying the detection value (angular velocity around the X ′ axis) of the X-axis angular velocity sensor 10a by a predetermined time, by substituting ⁇ xx ′ and ⁇ y into the equation (19), ⁇ 3x can be obtained.
  • ⁇ zz ′ can be obtained by multiplying the detection value (angular velocity about the Z ′ axis) of the Z-axis angular velocity sensor 10c by a predetermined time, by substituting ⁇ zz ′ and ⁇ y into the equation (20), ⁇ 1z can be obtained.
  • ⁇ xx ′ can be obtained by multiplying the detection value (angular velocity around the X ′ axis) of the X-axis angular velocity sensor 10a by a predetermined time, by substituting ⁇ xx ′ and ⁇ z into Equation (21), ⁇ 2x can be obtained.
  • ⁇ yy ′ can be obtained by multiplying the detection value (angular velocity around the Y ′ axis) of the Y-axis angular velocity sensor 10b by a predetermined time, by substituting ⁇ yy ′ and ⁇ z into the equation (22), ⁇ 1y can be obtained.
  • Equation (13) The inverse matrices T ⁇ x ⁇ 1 , T ⁇ y ⁇ 1 , and T ⁇ z ⁇ 1 calculated from ⁇ 2x , ⁇ 3x , ⁇ 1y , ⁇ 3y , ⁇ 1z , and ⁇ 2z obtained as described above are expressed by Equation (13).
  • the angular velocity sensor 10a, 10b each of the detection values G x of 10c ', G y', G z ' each ideal value G x, G y, can be corrected to G z.
  • the correction calculation of Expression (13) is performed with a digital value by a CPU or a dedicated circuit. Therefore, the X ′ axis and the Y ′ axis obtained by multiplying the A / D conversion values of the detection values G x ′, G y ′, and G z ′ of the angular velocity sensors 10a, 10b, and 10c by the A / D conversion sample period ⁇ t. , Z ′ axis micro rotation angles ⁇ x ′ , ⁇ y ′ , ⁇ z ′ are expressed by the following formula (23) as micro rotation angles ⁇ x , ⁇ y , ⁇ about the X axis, Y axis, and Z axis. It is corrected to z .
  • the X-axis acceleration sensor 50a is mounted such that a small angle ⁇ 2x around the Y axis and an X ′ axis rotated by a small angle ⁇ 3x around the Z axis are detection axes ( ⁇ 2x and ⁇ 3x are mounting angle errors).
  • the transformation matrix T ⁇ x is expressed as the following equation (24).
  • the Y-axis acceleration sensor 50b is mounted so that the Y ′ axis rotated by a minute angle ⁇ 1y around the X axis and the minute angle ⁇ 3y around the Z axis becomes the detection axis ( ⁇ 1y and ⁇ 3y are
  • the conversion matrix T ⁇ y is expressed as the following equation (25).
  • the Z-axis acceleration sensor 50c is mounted such that the Z ′ axis rotated by a minute angle ⁇ 1z around the X axis and the minute angle ⁇ 2z around the Y axis becomes the detection axis ( ⁇ 1z and ⁇ 2z are
  • the conversion matrix T ⁇ z is expressed as the following equation (26).
  • Expressions (24) to (26) correspond to Expressions (10) to (12) in the angular velocity sensors 10a, 10b, and 10c.
  • Expressions (27) to (29) correspond to Expressions (14) to (16) in the angular velocity sensors 10a, 10b, and 10c. Then, the following (30) to (35) are obtained by the same method as the derivation of the expressions (17) to (22) in the angular velocity sensors 10a, 10b, and 10c.
  • Equation (36) corresponds to Equation (13) in angular velocity sensors 10a, 10b, and 10c.
  • the correction calculation of Expression (36) is performed with a digital value by the CPU or a dedicated circuit. Therefore, the X ′ axis and the Y ′ axis obtained by multiplying the A / D conversion values of the detection values A x ′, A y ′, and A z ′ of the acceleration sensors 50a, 50b, and 50c by the A / D conversion sample period ⁇ t. , Z ′ axis direction micro speeds ⁇ v x ′ , ⁇ v y ′ , ⁇ v z ′ are converted into X axis, Y axis, Z axis direction micro speeds ⁇ v x , ⁇ v y , ⁇ v z by the following equation (37). It is corrected.
  • FIG. 8 is a diagram showing the configuration of the correction parameter creation device of the present embodiment.
  • the correction parameter creation apparatus 200 is used to create a correction parameter (correction matrix) for correcting a detection value including an error caused by a mounting angle error of each sensor included in the attitude detection apparatus 1 to an ideal value.
  • the correction parameter creating apparatus 200 includes a cube jig 210, a socket 220, a rotating plate 230, a rotating motor 240, a support base 250, a cable 260, and the like.
  • the cube jig 210 is molded into a cubic shape (or a rectangular parallelepiped shape) using a material such as metal, and is chamfered so that the three surfaces 211, 212, and 213 are orthogonal to each other.
  • the socket 220 is fixed.
  • the three surfaces 211, 212, and 213 of the cube jig 210 correspond to the first surface, the second surface, and the third surface of the jig in the present invention, respectively.
  • the socket 220 includes a socket body 222 and a lid body 224 that can be opened and closed.
  • the socket body 220 can accommodate the attitude detection device 1 in a predetermined direction without a gap.
  • the cube jig 210 can fix the posture detection device 1 by setting the posture detection device 1 in the socket 220 so that the X axis, the Y axis, and the Z axis are perpendicular to the surfaces 212, 213, and 211, respectively. It has become.
  • fixing brackets (not shown) are provided on the three surfaces 214, 215, and 216 respectively facing the surfaces 211, 212, and 213 of the cube jig 210.
  • the rotating plate 230 is so small that the unevenness of the upper surface 231 is negligible, and a fixing bracket (not shown) is provided on the upper surface 231, and any fixing bracket of the cube jig 210 is used as the fixing bracket of the rotating plate 230.
  • any fixing bracket of the cube jig 210 is used as the fixing bracket of the rotating plate 230.
  • the rotation plate 230 can be adjusted in inclination, and is strictly adjusted so that the upper surface 231 of the rotation plate 230 is horizontal in a state where the data correction device 200 is installed.
  • the rotary motor 240 is attached to the support base 250, and can rotate at an angular velocity in a predetermined range clockwise or counterclockwise about the vertical direction.
  • the cable 260 is connected to a control circuit (not shown) of the rotary motor 240.
  • a control device such as a personal computer is connected to the cable 260 so that the rotational speed of the rotary motor 250 can be adjusted by an interface such as GPIB (General / Purpose / Interface / Bus).
  • the rotary motor 240 functions as a rotation control unit in the present invention.
  • the use of the cube jig 210 and the rotating plate 230 eliminates the need for the rotating arm 530 shown in FIG. 13A and the like, thereby providing a correction parameter generating apparatus 200 that is more compact and low cost. be able to.
  • correction matrices T ⁇ x ⁇ 1 , T ⁇ y ⁇ 1 , T ⁇ z ⁇ 1 , T ⁇ x ⁇ 1 , T ⁇ y ⁇ 1 , T ⁇ z -1 can be easily obtained in a short time.
  • FIG. 9 is a flowchart showing an example of a correction parameter creation procedure in the present embodiment.
  • the rotating plate 230 is installed so that the upper surface 231 is horizontal (step S10).
  • the posture detection device 1 is set in the socket 220 attached to the surface 211 of the cube jig 210 (step S12).
  • step S14 the surface 215 facing the surface 212 of the cube jig 210 is fixed to the upper surface 231 of the rotating plate 230 (step S14).
  • the correction parameter creating apparatus 200 is set as shown in FIG. 10A, and the cube jig 210 is fixed to the rotating plate 230 so that the positive direction of the X axis is vertically upward.
  • Step S18 the minute rotation angles ⁇ y ′ , ⁇ z ′ obtained by sampling the detected values G y ′, G z ′ of the angular velocity sensors 10 b, 10 c and A / D converting the sample period ⁇ t, respectively.
  • ⁇ y ′ and ⁇ z ′ correspond to ⁇ yy ′ in equation (17) and ⁇ zz ′ in equation (18), respectively.
  • step S20 the surface 216 facing the surface 213 of the cube jig 210 is fixed to the upper surface 231 of the rotating plate 230 (step S20).
  • the correction parameter creating apparatus 200 is set as shown in FIG. 10B, and the cube jig 210 is fixed to the rotating plate 230 so that the positive direction of the Y-axis is vertically upward.
  • Step S22 the detected values of the X-axis acceleration sensor 50a and the Z-axis acceleration sensor 50c are acquired while the rotating plate 230 is stationary, and ⁇ 3x and ⁇ 1z are calculated from the equations (32) and (33) ( Step S22). Since the specific process in step S22 is the same as that in step S16, description thereof is omitted.
  • Step S24 Since the specific process in step S24 is the same as that in step S18, description thereof is omitted.
  • step S26 the surface 214 facing the surface 211 of the cube jig 210 is fixed to the upper surface 231 of the rotating plate 230 (step S26).
  • the correction parameter creating apparatus 200 is set as shown in FIG. 10C, and the cube jig 210 is fixed to the rotating plate 230 so that the positive direction of the Z-axis is vertically upward.
  • Step S28 the detected values of the X-axis acceleration sensor 50a and the Y-axis acceleration sensor 50b are acquired while the rotating plate 230 is stationary, and ⁇ 2x and ⁇ 1y are calculated from the equations (34) and (35) ( Step S28). Since the specific process in step S28 is the same as that in step S16, description thereof is omitted.
  • step S30 Since the specific process in step S30 is the same as that in step S18, description thereof is omitted.
  • correction matrices T ⁇ x ⁇ 1 , T ⁇ y ⁇ 1 , T ⁇ z ⁇ 1 , T ⁇ x ⁇ 1 , T ⁇ y ⁇ 1 , T ⁇ z ⁇ 1 are created (step S32).
  • the correction matrix T ⁇ x ⁇ 1 is created by calculating the inverse matrix of the transformation matrix T ⁇ x obtained by substituting ⁇ 3x and ⁇ 2x calculated in steps S24 and S30, respectively, into Equation (10). be able to.
  • the correction matrix T ⁇ y ⁇ 1 can be created by calculating the inverse matrix of the transformation matrix T ⁇ y obtained by substituting ⁇ 3y and ⁇ 1y calculated in steps S18 and S30, respectively, into the equation (11).
  • the correction matrix T ⁇ z ⁇ 1 can be created by calculating the inverse matrix of the transformation matrix T ⁇ z obtained by substituting ⁇ 2z and ⁇ 1z calculated in steps S18 and S24, respectively, into the equation (12). it can.
  • the correction matrix T ⁇ x ⁇ 1 can be created by calculating the inverse matrix of the transformation matrix T ⁇ x obtained by substituting ⁇ 3x and ⁇ 2x calculated in steps S22 and S28, respectively, into the equation (24). it can.
  • the correction matrix T ⁇ y ⁇ 1 can be created by calculating the inverse matrix of the transformation matrix T ⁇ y obtained by substituting ⁇ 3y and ⁇ 1y calculated in steps S16 and S28 respectively into the equation (25). it can.
  • the above processing is performed by a personal computer or the like connected to the cable 260 of the correction parameter creating apparatus 200.
  • the correction parameter created using this embodiment is used in the task for the correction calculation process mounted in the user side microcomputer connected to the back
  • the detected values of the acceleration sensors 50a, 50b, and 50c are rotated by rotating the rotating plate 230.
  • the detection values of the angular velocity sensors 10a, 10b, and 10c can be easily acquired in a short time.
  • the rotating plate 230 is installed so that the upper surface 231 is horizontal only once at the beginning, the rotation direction of the rotating plate 230 is fixed, so that the detection values for the X axis, the Y axis, and the Z axis are acquired. Setting time can be greatly reduced. Therefore, according to the present embodiment, the correction matrices T ⁇ x ⁇ 1 , T ⁇ y ⁇ 1 , T ⁇ z ⁇ 1 , T ⁇ x ⁇ 1 , T ⁇ y ⁇ 1 , and T ⁇ z ⁇ 1 can be created at a lower cost. .
  • the correction matrix T ⁇ x ⁇ 1 of the X-axis angular velocity sensor 10a is created based on the detection values of the X-axis angular velocity sensor 10a during rotation around the Y-axis and the Z-axis. Therefore, a correction parameter is created in consideration of both the mounting angle error and the other axis sensitivity error for each detection axis.
  • Correction parameters created using the present embodiment include a device for performing posture detection and posture control of a moving body and a robot, a head mounted display used for virtual reality, a tracker for detecting a head posture angle, a 3D game pad, and the like Can be used to correct the detection value of a posture detection device incorporated in various electronic devices such as game machines, digital cameras, mobile phones, portable information terminals, car navigation systems and the like.
  • FIG. 11 is a diagram showing the configuration of the posture detection device of the present embodiment.
  • the posture detection apparatus 300 with a correction function includes an angular velocity sensor module 2, an acceleration sensor module 3, anti-alias filters 310a, 310b, 310c, 350a, 350b, 350c, and A / D conversion circuits 320a, 320b, 320c, 360a, 360b, 360c.
  • a correction calculation processing unit 370 and a storage unit 380 are included.
  • the anti-alias filters 310a, 310b, 310c, 350a, 350b, and 350c are disposed in front of the A / D conversion circuits 320a, 320b, 320c, 360a, 360b, and 360c, respectively, and angular velocity detection signals 38a, 38b, and 38c and acceleration detection are performed.
  • the signals 78a, 78b, and 78c are attenuated to a level that can be ignored in advance by returning to the frequency band near DC by sampling the A / D conversion circuits 320a, 320b, 320c, 360a, 360b, and 360c, respectively.
  • the anti-alias filters 310a, 310b, 310c, 360a, 360b, 360c can be configured as, for example, switched capacitor filters (Switched Capacitor Filter (SCF)).
  • SCF Switchched Capacitor Filter
  • the A / D conversion circuits 320a, 320b, 320c, 360a, 360b, and 360c include angular velocity detection signals 38a, 38b, and 38c and acceleration detection signals 78a, 78b, and 78c, and anti-alias filters 310a, 310b, 310c, 350a, 350b, and 350c. Are converted into angular velocity detection signals 322a, 322b, and 322c and acceleration detection signals 362a, 362b, and 362c having a predetermined number of bits, respectively.
  • the A / D conversion circuits 320a, 320b, 320c, 360a, 360b, and 360c function as an A / D conversion processing unit in the present invention, and are of a flash type (parallel comparison type), a pipeline type, a successive approximation type, and a delta sigma type. It can be constituted by various known types of AD conversion circuits.
  • the storage unit 380 stores an angular velocity sensor correction parameter 382 and an acceleration sensor correction parameter 384.
  • the correction parameters 382 are correction matrices T ⁇ x ⁇ 1 , T ⁇ y ⁇ 1 , T ⁇ z ⁇ 1
  • the correction parameters 384 are correction matrices T ⁇ x ⁇ 1 , T ⁇ y ⁇ 1 , T ⁇ z ⁇ 1 . .
  • the correction calculation processing unit 370 calculates the correction equation (23) based on the angular velocity detection signals 322a, 322b, and 322c and the correction parameter 382, so that the angular velocity detection signal caused by the mounting angle error of the angular velocity sensors 10a, 10b, and 10c.
  • Angular velocity detection signals 302a, 302b, and 302c in which errors of 38a, 38b, and 38c are corrected are generated.
  • the correction calculation processing unit 370 obtains a value obtained by multiplying the digital value of the angular velocity detection signals 322a, 322b, and 322c by the A / D conversion sample period ⁇ t, and the minute rotation angle ⁇ x ′ of the correction equation (23), ⁇ y ', ⁇ z' each assignment to small rotation angle ⁇ x, ⁇ y, calculates the [Delta] [theta] z, further small rotation angle ⁇ x, ⁇ y, the angular velocity corresponding to [Delta] [theta] z digital value divided by ⁇ t Detection signals 302a, 302b, and 302c are generated.
  • the correction calculation processing unit 370 calculates the correction equation (37) based on the acceleration detection signals 362a, 362b, 362c and the correction parameter 384, thereby causing the mounting angle error of the acceleration sensors 50a, 50b, 50c.
  • Angular velocity detection signals 304a, 304b, and 304c in which errors in the acceleration detection signals 78a, 78b, and 78c are corrected are generated.
  • the correction calculation processing unit 370 multiplies the digital values of the acceleration detection signals 362a, 362b, and 362c by the A / D conversion sample period ⁇ t, and the micro speeds ⁇ v x ′ and ⁇ v of the correction formula (37).
  • Acceleration detection signal 304a corresponding to a digital value obtained by substituting y ′ and ⁇ v z ′ for subtle velocities ⁇ v x , ⁇ v y , and ⁇ v z and further dividing the micro velocities ⁇ v x , ⁇ v y , and ⁇ v z by ⁇ t. , 304b, 304c are generated.
  • the correction calculation processing unit 370 can be realized as a dedicated circuit for performing the correction calculation processing, or the function of the correction calculation processing unit 370 is executed by a CPU (CentralCPUProcessing ⁇ ⁇ Unit) executing a program stored in the storage unit 380 or the like. It can also be realized.
  • a CPU CentralCPUProcessing ⁇ ⁇ Unit
  • FIG. 12 is a diagram illustrating another configuration of the attitude detection device according to the present embodiment.
  • the configurations of the angular velocity sensor module 2, the acceleration sensor module 3, the anti-alias filters 310a, 310b, 310c, 350a, 350b, 350c, and the storage unit 380 are the same as those in FIG.
  • the multiplexer 390 time-divides the signals obtained by filtering the angular velocity detection signals 38a, 38b, and 38c and the acceleration detection signals 78a, 78b, and 78c by the anti-alias filters 310a, 310b, 310c, 350a, 350b, and 350c in a predetermined cycle. Select sequentially.
  • the correction calculation processing unit 370 samples the detection signal 322 every predetermined period, and when the detection signal 322 corresponds to the angular velocity detection signals 38a, 38b, and 38c, the correction formula (23) is based on the detection signal 322 and the correction parameter 382. ) Is generated, and the detected angular velocity detection signals are generated and output as detection signals 322 in a time-sharing manner.
  • the correction calculation processing unit 370 is corrected by calculating the correction formula (37) based on the detection signal 322 and the correction parameter 384.
  • the acceleration detection signals are generated and output as detection signals 322 in a time-sharing manner.
  • the attitude detection device 300 shown in FIG. 11 or FIG. 12 if the bypass calculation mode 370 is bypassed and an output of the A / D conversion circuit 320a and the like can be output to the outside, a bypass mode is provided.
  • the correction parameters 382 and 384 can be created by using the correction parameter creation method of the present embodiment in the state set as above. That is, the posture detection apparatus 300 can also be an application target of the correction parameter creation method of the present embodiment.
  • the correction calculation processing unit 370 calculates the correction value according to the correction formula (23) and the correction formula (37) that can improve the correction accuracy and speed up the correction calculation processing. Therefore, it is possible to realize an attitude detection device with higher correction accuracy and faster correction calculation processing.
  • the posture detection device 300 of the present embodiment digitizes and outputs the sensor detection signal, it is not necessary to connect an A / D conversion circuit between the posture detection device 300 and the user side microcomputer.
  • the posture detection device 300 uses a device for performing posture detection and posture control of a moving body or robot, a head-mounted display used for virtual reality, a tracker for detecting the posture angle of the head, a 3D game pad, and the like. It can be incorporated into various electronic devices such as game machines, digital cameras, mobile phones, portable information terminals, car navigation systems.
  • the present invention is not limited to the present embodiment, and various modifications can be made within the scope of the gist of the present invention.
  • the posture detection apparatus shown in FIG. 1 includes three angular velocity sensors 10a, 10b, and 10c and three acceleration sensors 50a, 50b, and 50c, and posture detection to which the correction parameter creation method of the present invention is applied.
  • the device is not limited to this. That is, the posture detection device to which the correction parameter creation method of the present invention is applied may be configured to be able to detect orthogonal three-axis angular velocities or accelerations.
  • the posture detection device includes only the angular velocity sensors 10a, 10b, and 10c.
  • An attitude detection apparatus including only the acceleration sensors 50a, 50b, and 50c, an attitude detection apparatus including only the angular velocity sensors 10a and 10b and the acceleration sensor 50c, an attitude detection apparatus including only the angular velocity sensor 10a and the acceleration sensors 50b and 50c, and the like. include.
  • the correction parameter creating apparatus 200 shown in FIG. 8 has only one socket 220 attached to the surface 211, a plurality of sockets 220 may be attached to the surface 211. Then, by setting the posture detection device 1 in each of the plurality of sockets 220, the detection values of the posture detection devices 1 can be acquired simultaneously.
  • the detection values of the posture detection device 1 are acquired in the order of the X axis, the Y axis, and the Z axis. Can be acquired.
  • the present invention includes substantially the same configuration (for example, a configuration having the same function, method and result, or a configuration having the same purpose and effect) as the configuration described in the embodiment.
  • the invention includes a configuration in which a non-essential part of the configuration described in the embodiment is replaced.
  • the present invention includes a configuration that exhibits the same operational effects as the configuration described in the embodiment or a configuration that can achieve the same object.
  • the invention includes a configuration in which a known technique is added to the configuration described in the embodiment.
  • 1 attitude detection device 2 angular velocity sensor module, 3 acceleration sensor module, 4 package, 5a, 5b, 5c package surface, 10a, 10b, 10c angular velocity sensor, 11a, 11b, 11c vibrator, 12a, 12b, 12c drive electrode 13a, 13b, 13c drive electrode, 14a, 14b, 14c detection electrode, 15a, 15b, 15c detection electrode, 20a, 20b, 20c drive circuit, 21a current-voltage converter (I / V converter), 22a AC amplifier, 23a automatic gain control circuit (AGC) 24a comparator, 30a, 30b, 30c detection circuit, 31a, 32a charge amplifier, 33a differential amplifier, 34a AC amplifier, 35a synchronous detection circuit, 36a DC amplifier, 37a integration circuit (LPF), 3 a, 38b, 38c detection signal, 41a drive vibration arm, 42a detection vibration arm, 43a wide part, 44a drive base, 45a connection arm, 46a wide part, 47a detection base, 50a, 50b, 50c acceleration sensor, 51

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • General Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Human Computer Interaction (AREA)
  • Gyroscopes (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Position Input By Displaying (AREA)

Abstract

 回転板230を上面231が水平になるように設置し(S10)、立方体治具210の面211に、X軸(第1の軸)が面212(第2の面)に垂直になり、Y軸(第2の軸)が面213(第3の面)に垂直になり、Z軸(第3の軸)が面211(第1の面)に垂直になるように、姿勢検出装置1を固定する(S12)。そして、立方体治具の面212、213、211と対向する各面を回転板の上面に順次固定し(S14、S20、S26)、回転板を静止又は所定の角速度で回転させて姿勢検出装置の検出値を取得し(S16、S18、S22、S24、S28、S30)、補正パラメーターを作成する(S32)。

Description

姿勢検出装置の補正パラメーター作成方法、姿勢検出装置の補正パラメーター作成用装置及び姿勢検出装置
 本発明は、3軸の角速度又は加速度を検出するセンサーを含む姿勢検出装置の検出値を所定の直交座標系における検出値に補正するための補正パラメーター作成方法、補正パラメーター作成用装置及び補正機能付き姿勢検出装置に関する。
 近年、角速度センサーや加速度センサーにより物体の姿勢を検出する姿勢検出装置が様々な用途に使用されている。例えば、特開平9-106322号公報には、ユーザーの頭の姿勢を検出することにより、目前にあるディスプレイに表示された映像が頭の動きに連動して変化し、仮想空間を体験できるヘッドマウントディスプレイが記載されている。ユーザーの頭の姿勢角にあった映像がヘッドマウントディスプレイに映し出される。この姿勢角を検出するために角速度センサーや加速度センサーを備えた姿勢検出装置がヘッドマウントディスプレイの所定の位置に取り付けられる。姿勢検出装置を取り付ける場合、センサーの検出軸が頭の姿勢角を表すための直交座標系の3軸とそれぞれ平行になるように取り付けられていなければ、この取付角の誤差に起因して姿勢検出装置の検出値が誤差を含むことになる。そのため、姿勢検出装置を取り付ける位置や角度が厳密に規定される。
 しかし、姿勢検出装置の内部にセンサーが取り付けられる際に取付角にわずかなずれがあると、姿勢検出装置を取り付ける位置や角度が厳密に規定されていたとしても高精度の検出結果を得ることができない。センサーの取付角誤差が無いようにするのはコスト面から現実的ではないため、あらかじめ取付角誤差を算出して姿勢検出装置の検出値を取付角誤差に応じた補正パラメーターで補正することが行われる。式(1)及び式(2)は、それぞれ補正パラメーターを用いた角速度センサー用の補正式及び加速度センサー用の補正式を表している。
Figure JPOXMLDOC01-appb-M000001
Figure JPOXMLDOC01-appb-M000002
 式(1)において、関数行列式(ヤコビアン)JfGは角速度センサー用の補正パラメーターであり、f(x)及びf(p)はそれぞれ角速度センサーの今回及び前回の補正値(理想値)である。同様に、式(2)において、関数行列式(ヤコビアン)JfAは加速度センサー用の補正パラメーターであり、f(x)及びf(p)はそれぞれ加速度センサーの今回及び前回の補正値(理想値)である。また、式(1)、式(2)において、x及びpはそれぞれ角速度センサー又は角速度センサーの今回及び前回の検出値であり、oはランダウの記号である。
 センサーの取付角誤差は姿勢検出装置毎に異なるため、出荷テスト時等に各姿勢検出装置に対して式(1)、式(2)の補正パラメーター(JfG、JfA)が作成される。図13A~図13C及び図14A~図14Cは、補正パラメーター(JfG、JfA)を作成する従来方法を示している。従来方法では、まず、テーブル510に取り付けられたソケット520に姿勢検出装置をセットし、図13A~図13Cに示す順に、回転腕530をX軸回り、Y軸回り、Z軸回りに所定の角速度で回転させて姿勢検出装置の各検出値を取得し、各検出値と各理想値を式(1)に代入して得られる連立方程式を解いて角速度センサー用の補正パラメーターを作成する。さらに、図14A~図14Cに示す順に、回転腕530を操作してX軸、Y軸、Z軸の正方向が鉛直上向きになる状態(鉛直下向きに重力加速度が加わる状態)で静止させて姿勢検出装置の各検出値を取得し、各検出値と各理想値を式(2)に代入して得られる連立方程式を解いて加速度センサー用の補正パラメーターを作成する。
 図13A~図13C及び図14A~図14Cに示す回転腕530の操作において、X軸、Y軸、Z軸に対してテーブル510を所定の角度に正確に固定しなければ、角速度センサーおよび加速度センサーの取付角誤差が正確に反映された検出値を取得することができない。しかし、X軸、Y軸、Z軸に対してテーブル510を所定の角度に正確に固定するためには、テーブル510及び回転腕530を含む補正パラメーター作成用装置500が大がかりな装置になる傾向がある。また、X軸、Y軸、Z軸に対してテーブル510を所定の角度に正確に固定するために相当の時間がかかる。角速度センサーと加速度センサーをともに含む姿勢検出装置に対する補正パラメーターを作成するためには、図13A~図13Cに示す回転腕530の操作と図14A~図14Cに示す回転腕530の操作を別々に行う必要があるためさらに時間がかかる。そのため、従来方法では、補正パラメーターの作成にかかるコストが大きいという問題がある。
 本発明は、以上のような問題点に鑑みてなされたものであり、本発明のいくつかの態様によれば、センサーの取付角誤差に起因する検出値の誤差を補正するための補正パラメーターをより低コストで作成することができる姿勢検出装置の補正パラメーター作成方法、より低コストで実現可能な補正パラメーター作成用装置及び補正機能付きの姿勢検出装置を提供することができる。
 (1)本発明は、検出軸が互いに直交する第1の軸、第2の軸及び第3の軸とそれぞれ略平行になるように取り付けられ角速度又は加速度を検出する第1のセンサー、第2のセンサー及び第3のセンサーを含み、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの検出信号に基づいて物体の姿勢を検出する姿勢検出装置の検出値を、前記第1の軸、前記第2の軸及び前記第3の軸を座標軸とする直交座標系における検出値に補正する補正式の補正パラメーターを作成する方法であって、上面が水平になるように回転板を設置するステップと、互いに直交する第1の面、第2の面、第3の面を有する直方体形状の治具の前記第1の面に、前記第1の軸が前記第2の面に垂直になり、前記第2の軸が前記第3の面に垂直になり、前記第3の軸が前記第1の面に垂直になるように、前記姿勢検出装置を固定するステップと、前記治具の前記第2の面と対向する面を前記回転板の前記上面に固定し、前記回転板を静止又は所定の角速度で回転させて前記姿勢検出装置の検出値を取得する第1検出値取得ステップと、前記治具の前記第3の面と対向する面を前記回転板の前記上面に固定し、前記回転板を静止又は所定の角速度で回転させて前記姿勢検出装置の検出値を取得する第2検出値取得ステップと、前記治具の前記第1の面と対向する面を前記回転板の前記上面に固定し、前記回転板を静止又は所定の角速度で回転させて前記姿勢検出装置の検出値を取得する第3検出値取得ステップと、取得した検出値に基づいて、前記補正パラメーターを作成する補正パラメーター作成ステップと、を含むことを特徴とする。
 X軸、Y軸、Z軸を座標軸とする直交座標系を考えた場合、第1の軸、第2の軸、第3の軸とX軸、Y軸、Z軸の対応関係は特に限定されない。
 本発明によれば、直方体形状の治具を用いるので、第1の軸、第2の軸、第3の軸が治具のそれぞれ第1の面、第2の面、第3の面と垂直になるように第1の面に姿勢検出装置を固定することが容易である。そして、上面が水平になるように回転板を設置すれば、治具の第2の面、第3の面、第1の面と対向する面をそれぞれ回転板の上面に固定するだけで、簡単に、第1の軸、第2の軸、第3の軸をそれぞれ鉛直方向と平行にすることができる。さらに、第1の軸、第2の軸、第3の軸をそれぞれ鉛直方向と平行にした状態において、回転板を静止又は回転させることにより、加速度センサー又は角速度センサーの検出値を簡単に短時間で取得することができる。
 すなわち、最初に一度だけ、上面が水平になるように回転板を設置すれば回転板の回転方向が固定されるので、第1の軸、第2の軸、第3の軸に関する検出値を取得するためのセッティング時間を大幅に短縮することができる。従って、本発明によれば、センサーの取付角誤差に起因する検出値の誤差を補正するための補正パラメーターをより低コストで作成することができる。
 (2)本発明の姿勢検出装置の補正パラメーター作成方法において、前記補正式は、前記補正パラメーターとして前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出値を前記直交座標系における各検出値に補正するための第1の補正行列、第2の補正行列及び第3の補正行列を含み、前記第1の補正行列、前記第2の補正行列及び前記第3の補正行列の各々と、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出値がA/D変換されたデジタル値をそれぞれ要素として含む行列の各々と、の積により得られる3つの行列の和として与えられるようにしてもよい。
 (3)本発明の姿勢検出装置の補正パラメーター作成方法において、前記第1の補正行列、前記第2の補正行列及び前記第3の補正行列は、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出軸をそれぞれ前記第1の軸、前記第2の軸及び前記第3の軸に変換する回転行列の逆行列であるようにしてもよい。
 (4)本発明の姿勢検出装置の補正パラメーター作成方法において、前記補正パラメーター作成ステップは、前記第1検出値取得ステップにおいて取得した前記検出値に基づいて、前記第2のセンサー及び前記第3のセンサーの前記第1の軸回りの各取付角誤差を算出するステップと、前記第2検出値取得ステップにおいて取得した前記検出値に基づいて、前記第1のセンサー及び前記第3のセンサーの前記第2の軸回りの各取付角誤差を算出するステップと、前記第3検出値取得ステップにおいて取得した前記検出値に基づいて、前記第1のセンサー及び前記第2のセンサーの前記第3の軸回りの各取付角誤差を算出するステップと、前記第1のセンサーの前記第2の軸回りの前記取付角誤差及び前記第3の軸回りの前記取付角誤差に基づいて、前記第1の補正行列を作成するステップと、前記第2のセンサーの前記第1の軸回りの前記取付角誤差及び前記第3の軸回りの前記取付角誤差に基づいて、前記第2の補正行列を作成するステップと、前記第3のセンサーの前記第1の軸回りの前記取付角誤差及び前記第2の軸回りの前記取付角誤差に基づいて、前記第3の補正行列を作成するステップと、を含むようにしてもよい。
 (5)本発明は、検出軸が互いに直交する第1の軸、第2の軸及び第3の軸とそれぞれ略平行になるように取り付けられ角速度又は加速度を検出する第1のセンサー、第2のセンサー及び第3のセンサーを含み、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの検出信号に基づいて物体の姿勢を検出する姿勢検出装置の検出値を、前記第1の軸、前記第2の軸及び前記第3の軸を座標軸とする直交座標系における検出値に補正する補正式の補正パラメーターを作成するために使用される補正パラメーター作成用装置であって、互いに直交する第1の面、第2の面、第3の面を有し、前記第1の面に、前記第1の軸が前記第2の面に垂直になり、前記第2の軸が前記第3の面に垂直になり、前記第3の軸が前記第1の面に垂直になるように前記姿勢検出装置を固定可能な直方体形状の治具と、上面に、前記治具の前記第1の面、前記第2の面、前記第2の面とそれぞれ対向する面のいずれかを固定可能な回転板と、前記回転板を所定の角速度で回転させる回転制御部と、を含むことを特徴とする。
 本発明によれば、直方体形状の治具と回転板を用いることにより、回転腕を必要としないので、よりコンパクトかつ低コストの補正パラメーター作成用装置を提供することができる。本発明に係る補正パラメーター作成用装置を用いることにより、上記の通り、姿勢検出装置に取り付けられた各センサーの検出値の補正パラメーターを簡単に短時間で取得することができる。
 (6)本発明は、検出軸が互いに直交する第1の軸、第2の軸及び第3の軸とそれぞれ略平行になるように取り付けられ角速度又は加速度を検出する第1のセンサー、第2のセンサー及び第3のセンサーと、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出値を前記第1の軸、前記第2の軸及び前記第3の軸を座標軸とする直交座標系における検出値に補正する補正式の補正パラメーターが記憶された記憶部と、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出信号をデジタル信号に変換する処理を行うA/D変換処理部と、前記デジタル信号の各々と前記補正パラメーターに基づいて前記補正式を計算する処理を行う補正計算処理部と、を含み、前記補正式は、前記補正パラメーターとして前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出値を前記直交座標系における各検出値に補正するための第1の補正行列、第2の補正行列及び第3の補正行列を含み、前記第1の補正行列、前記第2の補正行列及び前記第3の補正行列の各々と、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出値がA/D変換されたデジタル値をそれぞれ要素として含む行列の各々と、の積により得られる3つの行列の和として与えられることを特徴とする。
 従来の補正式(1)、補正式(2)における関数行列式(ヤコビアン)はセンサーの取付角誤差を直接的に反映する補正パラメーターではなく、また、補正式(1)、補正式(2)では前回の検出値を元に関数行列式(ヤコビアン)を用いて今回の検出値を類推するため、検出値に何らかの写像を施すと補正値が得られるようになっていない。そのため、補正式(1)、補正式(2)では補正精度を上げるのに限界がある。
 本発明によれば、補正計算処理部が計算する補正式に含まれる3つの補正行列に各センサーの取付角誤差を直接的に反映させることができる。また、本発明によれば、補正計算処理部が計算する補正式は今回の検出値に対する補正値の計算において前回の検出値を必要としないので、今回の検出値が得られれば直ちに補正値を計算することができる。従って、本発明によれば、補正精度がより高く、かつ、補正計算処理がより速い姿勢検出装置を実現することができる。
 (7)本発明の姿勢検出装置において、前記第1の補正行列、前記第2の補正行列及び前記第3の補正行列は、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出軸をそれぞれ前記第1の軸、前記第2の軸及び前記第3の軸に変換する回転行列の逆行列であるようにしてもよい。
 (8)本発明の姿勢検出装置は、所定の周期で、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各前記検出信号のいずれか1つを順次選択する処理を行う信号選択処理部を含み、前記A/D変換処理部は、前記信号選択処理部が選択した検出値を順次A/D変換処理するA/D変換回路を含むようにしてもよい。
図1は、本実施形態の補正パラメーター作成方法の対象となる姿勢検出装置の構成の一例を示す図である。 図2は、本実施形態における姿勢検出装置の斜視図である。 図3は、角速度センサーに含まれる振動子の一例を示す平面図である。 図4は、角速度センサーに含まれる振動子の動作について説明するための図である。 図5は、角速度センサーに含まれる振動子の動作について説明するための図である。 図6は、角速度センサーに含まれる駆動回路及び検出回路の構成の一例を示す図である。 図7Aは、センサーの取付角誤差について説明するための図である。 図7Bは、センサーの取付角誤差について説明するための図である。 図7Cは、センサーの取付角誤差について説明するための図である。 図8は、本実施形態の補正パラメーター作成用装置の構成を示す図である。 図9は、本実施形態における補正パラメーター作成手順の一例を示すフローチャート図である。 図10Aは、本実施形態における補正パラメーター作成手順について説明するための図である。 図10Bは、本実施形態における補正パラメーター作成手順について説明するための図である。 図10Cは、本実施形態における補正パラメーター作成手順について説明するための図である。 図11は、本実施形態の姿勢検出装置の構成を示す図である。 図12は、本実施形態の姿勢検出装置の他の構成を示す図である。 図13Aは、従来の補正パラメーター作成方法について説明するための図である。 図13Bは、従来の補正パラメーター作成方法について説明するための図である。 図13Cは、従来の補正パラメーター作成方法について説明するための図である。 図14Aは、従来の補正パラメーター作成方法について説明するための図である。 図14Bは、従来の補正パラメーター作成方法について説明するための図である。 図14Cは、従来の補正パラメーター作成方法について説明するための図である。
 以下、本発明の好適な実施形態について図面を用いて詳細に説明する。なお、以下に説明する実施の形態は、特許請求の範囲に記載された本発明の内容を不当に限定するものではない。また以下で説明される構成の全てが本発明の必須構成要件であるとは限らない。
 なお、以下の説明において、本発明における第1の軸、第2の軸、第3の軸はそれぞれX軸、Y軸、Z軸に対応するものとするが、本発明における第1の軸、第2の軸、第3の軸とX軸、Y軸、Z軸の対応関係はこれに限らず、任意の対応関係とすることができる。
 1.姿勢検出装置
 1-1.姿勢検出装置の構成
 図1は、本実施形態の補正パラメーター作成方法の対象となる姿勢検出装置の構成の一例を示す図である。
 図1に示すように、本実施形態における姿勢検出装置1は、X軸、Y軸、Z軸回りの角速度を検出する角速度センサーモジュール2と、X軸、Y軸、Z軸方向の加速度を検出する加速度センサーモジュール3を含んで構成されている。
 角速度センサーモジュール2は、それぞれX軸、Y軸、Z軸回りの角速度を検出するX軸角速度センサー10a、Y軸角速度センサー10b、Z軸角速度センサー10cを含む。
 X軸角速度センサー10aは、振動子11a、振動子11aを振動させる駆動回路20a、角速度検出信号38aを生成する検出回路30aを含み、振動子11aの駆動電極12a、13aと駆動回路20a、振動子11aの検出電極14a、15aと検出回路30aがそれぞれ接続されている。
 同様に、Y軸角速度センサー10bは、振動子11b、振動子11bを振動させる駆動回路20b、角速度検出信号38bを生成する検出回路30bを含み、振動子11bの駆動電極12b、13bと駆動回路20b、振動子11bの検出電極14b、15bと検出回路30bがそれぞれ接続されている。
 同様に、Z軸角速度センサー10cは、振動子11c、振動子11cを振動させる駆動する駆動回路20c、角速度検出信号38cを生成する検出回路30cを含み、振動子11cの駆動電極12c、13cと駆動回路20c、振動子11cの検出電極14c、15cと検出回路30cがそれぞれ接続されている。
 加速度センサーモジュール3は、それぞれX軸、Y軸、Z軸方向の加速度を検出するX軸加速度センサー50a、Y軸加速度センサー50b、Z軸加速度センサー50cを含む。
 X軸加速度センサー50aは、振動子51a、振動子51aを振動させる駆動回路60a、加速度検出信号78aを生成する検出回路70aを含み、振動子51aの駆動電極52a、53aと駆動回路60a、振動子51aの検出電極54a、55aと検出回路70aがそれぞれ接続されている。
 同様に、Y軸加速度センサー50bは、振動子51b、振動子51bを振動させる駆動回路60b、加速度検出信号78bを生成する検出回路70bを含み、振動子51bの駆動電極52b、53bと駆動回路60b、振動子51bの検出電極54b、55bと検出回路70bがそれぞれ接続されている。
 同様に、Z軸加速度センサー50cは、振動子51c、振動子51cを振動させる駆動回路60c、加速度検出信号78cを生成する検出回路70cを含み、振動子51cの駆動電極52c、53cと駆動回路60c、振動子51cの検出電極54c、55cと検出回路70cがそれぞれ接続されている。
 なお、角速度センサー10a、10b、10cは、それぞれ本発明における第1のセンサー、第2のセンサー及び第3のセンサーとして機能する。同様に、加速度センサー50a、50b、50cは、それぞれ本発明における第1のセンサー、第2のセンサー及び第3のセンサーとして機能する。
 図2は、本実施形態における姿勢検出装置の斜視図である。
 図2に示すように、姿勢検出装置1において、角速度センサーモジュール2及び加速度センサーモジュール3は、それぞれ立方体(広義には直方体。以下、同じ。)の形状に形成されており、直方体形状のパッケージ4の中に収納されている。
 X軸、Y軸、Z軸は姿勢検出装置1を基準にして決定される。例えば、姿勢検出装置1を構成するパッケージ4が直方体形状である場合は、パッケージ4の直交する3つの面5a、5b、5cと垂直な軸をそれぞれX軸、Y軸、Z軸とすることができる。また、X軸、Y軸、Z軸の正方向は任意に決定することができ、本実施形態では、図2に示す矢印の先に向かう方向を各軸の正方向とする。
 1-2.角速度センサーモジュール
 図2に示すように、角速度センサーモジュール2において、角速度センサー10a、10b、10cは、検出軸がそれぞれX軸、Y軸、Z軸とほぼ平行になるように絶縁基板80の上に取り付けられ、それぞれ振動子11a、11b、11cがパッケージ82a、82b、82c内に収容されている。パッケージ82a、82b、82cの周囲は樹脂モールド材で覆われている。
 パッケージ82a、82b、82cは、それぞれパッケージ本体84aと蓋体86a、パッケージ本体84bと蓋体86b、パッケージ本体84cと蓋体86cから構成されている。パッケージ本体84a、84b、84cは、複数のセラミックシートを積層、焼結して直方体の箱状に形成されている。蓋体86a、86b、86cは、ガラス板や金属板、セラミックシート等で形成されており、金属ロウ材、低融点ガラスなどの接合材を介して、振動子11a、11b、11cがそれぞれ収容されたパッケージ本体84a、84b、84cの上面開口部を真空封止している。振動子11a、11b、11cは、絶縁基板80に形成された配線パターン(図示せず)によってそれぞれ駆動回路20a、20b、20cや検出回路30a、30b、30cと接続されている。
 駆動回路20aと検出回路30a、駆動回路20bと検出回路30b、駆動回路20cと検出回路30cは、3つのチップにIC化してそれぞれパッケージ82a、82b、82c内にそれぞれ収容されていてもよい。また、駆動回路20a、20b、20c、検出回路30a、30b、30cを1チップにIC化して絶縁基板80上に配置するようにしてもよい。
 なお、図2では図示を省略しているが、検出回路30a、30b、30cからの各検出信号38a、38b、38cは外部出力端子(図示せず)を介して姿勢検出装置1の外部に出力されるようになっている。
 図3は、角速度センサーに含まれる振動子の一例を示す平面図である。角速度センサー10a、10b、10cにそれぞれ含まれる振動子11a、11b、11cはすべて同一の構造であるため、図3では振動子11aの構造のみが図示されている。なお、図3におけるX軸、Y軸、Z軸は水晶の軸を示し、図2におけるX軸、Y軸、Z軸とは無関係である。
 振動子11aは、水晶などの圧電材料の薄板から形成され、駆動用基部44aから駆動振動腕41a(広義には、駆動用振動片)が水晶のY軸方向に延出している。駆動振動腕41aの側面及び上面にはそれぞれ駆動電極12a及び13aが形成されている。図1に示したように、駆動電極12a、13aは駆動回路20aに接続される。
 駆動用基部44aは、水晶のX軸方向に延びる連結腕45aを介して検出用基部47aに接続されている。検出振動腕42a(広義には、検出振動片)は、検出用基部47aから水晶のY軸方向に延出されている。検出振動腕42aの上面には検出電極14a及び15aが形成されており、検出振動腕42aの側面には電極16aが形成されている。図1に示したように、検出電極14a、15aは駆動回路20aに接続される。また、電極16aは接地される。
 駆動振動腕41aの駆動電極12aと駆動電極13aとの間に交番電圧/交番電流からなる駆動信号が与えられると、図4に示すように、駆動振動腕41aは圧電効果によって矢印Bのように屈曲振動する。
 ここで、図5に示すように振動子11aが水晶のZ軸を回転軸とした回転運動をすると、駆動振動腕41aは、矢印Bの屈曲振動の方向と水晶のZ軸の両方に垂直な方向にコリオリの力を得る。その結果、連結腕45aは矢印Cで示すような振動をする。そして、検出振動腕42aは、連結腕45aの振動(矢印C)に連動して、連結腕45aと矢印Dのような屈曲振動をする。
 そして、これらの屈曲振動に基づいて発生する逆圧電効果によって、検出振動腕42aの検出電極14a、15aと電極16aとの間には、それぞれ逆方向の交番電圧/交番電流が発生する。以上のようにして、振動子11aは、水晶のZ軸を検出軸としてコリオリの力に基づく角速度成分を検出し、検出電極14a、15aを介して検出信号を出力する。
 なお、図3の構成では、振動子11aのバランスを良くするために、検出用基部47aを中央に配置し、検出用基部47aから+Y軸と-Y軸の両方向に検出振動腕42aを延出させている。さらに、検出用基部47aから+X軸と-X軸の両方向に連結腕45aを延出させ、連結腕45aのそれぞれから、+Y軸と-Y軸の両方向に駆動振動腕41aを延出させている。
 また、駆動振動腕41aの先端を幅広の幅広部43aにし、さらに、錘を付けることでコリオリの力を大きくしている。また、錘効果によって、所望の共振周波数を、短い振動腕で得ることができる。同様の理由で、検出振動腕42aの先端を幅広の幅広部46aにし、さらに、錘を付けている。
 なお、振動子11aは、上述の構成に限らず、コリオリの力に基づく角速度成分を含む検出信号を出力する振動子であれば良い。例えば、駆動振動腕と検出振動腕とを兼ねる構成であっても良く、また、駆動振動腕や検出振動腕に圧電膜を形成した構成であっても良い。
 図6は、角速度センサーに含まれる駆動回路及び検出回路の構成の一例を示す図である。駆動回路20a、20b、20cはすべて同じ構成であり、検出回路30a、30b、30cはすべて同じ構成であるので、図3では駆動回路20a及び検出回路30aの構成のみを図示している。
 図6に示すように、駆動回路20aは、電流電圧変換器(I/V変換器)21a、AC増幅器22a、自動利得制御回路(AGC)23a、コンパレータ24aを含んで構成されている。
 振動子11aが振動すると、圧電効果に基づく交流電流がフィードバック信号として駆動電極13aから出力され、電流電圧変換器(I/V変換器)21aに入力される。電流電圧変換器(I/V変換器)21aは、入力された交流電流を振動子11aの振動周波数と同一の周波数の交流電圧信号に変換して出力する。
 電流電圧変換器(I/V変換器)21aから出力された交流電圧信号は、AC増幅器22aに入力される。AC増幅器aは、入力された交流電圧信号を増幅して出力する。
 AC増幅器22aから出力された交流電圧信号は自動利得制御回路(AGC)23aに入力される。自動利得制御回路(AGC)23aは、入力された交流電圧信号の振幅を一定値に保持するように利得を制御し、利得制御後の交流電圧信号を振動子11aの駆動電極12aに出力する。この駆動電極12aに入力される交流電圧信号により振動子11aが振動する。
 AC増幅器22aが増幅した交流電圧信号はコンパレータ24aに入力され、交流電圧信号の振幅中心を基準電圧として、交流電圧信号と基準電圧信号との比較結果に応じて出力レベルを切り替える方形波電圧信号を、検出回路30aの同期検波回路35aに出力する。
 図6に示すように、検出回路30aは、チャージアンプ31a、32a、差動増幅器33a、AC増幅器34a、同期検波回路35a、DC増幅器36a及び積分回路(LPF)37aを含んで構成されている。
 チャージアンプ31a、32aには、振動子11aにより検出された互いに逆位相の検出信号(交流電流)が検出電極12a、13aを介して入力される。そして、チャージアンプ31a、32aは、入力された検出信号(交流電流)を基準電圧を中心とする交流電圧信号に変換する。
 差動増幅器33aはチャージアンプ31aの出力信号とチャージアンプ32aの出力信号を差動増幅する。差動増幅器33aの出力信号は、さらにAC増幅器34aで増幅される。
 同期検波回路35aは、コンパレータ24aが出力する方形波電圧信号を基に、AC増幅器34aの出力信号を同期検波することにより角速度成分を抽出する。同期検波回路35aは、例えば、方形波電圧信号の電圧レベルが基準電圧よりも高い時はAC増幅器34aの出力信号をそのまま出力し、方形波電圧信号の電圧レベルが基準電圧よりも低い時はAC増幅器34aの出力信号を基準電圧に対して反転して出力するスイッチ回路として構成することができる。
 同期検波回路35aで抽出された角速度成分信号は、DC増幅器36aで増幅されて積分回路(LPF)37aに入力される。
 積分回路(LPF)37aは、DC増幅器35aの出力信号から高周波成分を減衰させて直流成分を抽出することにより角速度検出信号38aを生成して外部に出力する。
 1-3.加速度センサーモジュール
 図2に示すように、加速度センサーモジュール3は、ベース90、ウエイト100、3つの加速度センサー50a、50b、50cを有する。なお、図2では、図1に示した駆動回路60a、60b、60c及び検出回路70a、70b、70cの図示を省略しているがパッケージ4の中の適当な位置に配置され、検出回路70a、70b、70cからの各検出信号78a、78b、78cは外部出力端子(図示せず)を介して姿勢検出装置1の外部に出力されるようになっている。
 ベース90は、立方形を作るように3つの正方形の壁部を互いに直交させて形成され、X軸、Y軸、Z軸方向に互いに直交する3つの取付面91、92、93を有する。ウエイト100は、所定の質量を有する立方体からなり、互いに直交する3つの接合面101、102、103を有する。ベース90及びウエイト100は、例えばアルミニウム合金などの適当な材料を用いて形成される。
 本実施形態では、加速度センサー50a、50b、50cは、それぞれ、水晶などの圧電材料の薄板から形成された双音叉型の振動子51a、51b、51cを含んで構成されている。
 振動子51a、51b、51cは、検出軸がそれぞれX軸、Y軸、Z軸とほぼ平行になるように、一方の基端部56a、56b、56cがそれぞれベース90の素子取付面91、92、93に取り付けられ、ベース90の各壁部に垂直に支持されている。振動子51a、51b、51cの他方の基端部57a、57b、57cは、それぞれ素子取付面91、92、93に対応するウエイト100の素子接合面101~103に接合されている。これにより、ウエイト100が、X軸、Y軸、Z軸方向から振動子51a、51b、51cによって浮遊した状態に支持される。
 振動子51aの2つの駆動振動腕58aには、その上下主面及び両側面に駆動電極52a、53aが設けられており(図示省略)、駆動回路60aにより駆動電極52a、53a間に所定の交流電圧が印加されると、2つの駆動振動腕58aは互いに逆向きに即ち近接または離反する向きに所定の周波数で屈曲振動する。
 振動子51aを所定の周波数で振動させた状態で、加速度センサーモジュール3に外力が作用してウエイト100にX軸方向の加速度が加わると、その大きさ及び向きに対応して、振動子51aには長手方向(すなわち、X軸方向)に圧縮または引張する力が作用する。振動子51aの周波数は、圧縮する力が作用すると減少し、引張する力が作用すると増加するように変化する。従って、検出回路70aにより振動子51aにおける周波数の変化量を検出し、周波数の変化量からX軸方向に作用する荷重を算出することにより、ウエイト100に作用したX軸方向の加速度の大きさ及び向きを計算することができる。
 振動子51b、51cの構造も振動子51aの構造と同一であり、同様にしてY軸及びZ軸方向の加速度の大きさ及び向きを計算することができる。
 なお、駆動回路60a、60b、60cは図6に示した駆動回路20aと同様の構成であり、また、検出回路70a、70b、70cは周波数の変化量を検出する既知の回路と同様の構成にすることができるので、その説明を省略する。
 2.補正パラメーター作成方法
 2-1.センサー取付角誤差
 角速度センサー10a、10b、10cは、理想的には、検出軸がそれぞれX軸、Y軸、Z軸と厳密に平行になるように取り付けられる。同様に、加速度センサー50a、50b、50cは、理想的には、検出軸がそれぞれX軸、Y軸、Z軸と厳密に平行になるように取り付けられる。しかし、角速度センサー10a、10b、10c及び加速度センサー50a、50b、50cをそのように厳密に取り付けるのはコスト面から難しい。そのため、X軸角速度センサー10aは、実際には、図7Aに示すように、検出軸がY軸回りに微小角Δθ2x、Z軸回りに微小角Δθ3x回ったX’軸と平行になるように取り付けられている。同様に、Y軸角速度センサー10bは、実際には、図7Bに示すように、検出軸がZ軸回りに微小角Δθ3y、X軸回りに微小角Δθ1y回ったY’軸と平行になるように取り付けられており、Z軸角速度センサー10cは、実際には、図7Cに示すように、検出軸がX軸回りに微小角Δθ1z、Y軸回りに微小角Δθ2z回ったZ’軸と平行になるように取り付けられている。すなわち、X軸角速度センサー10aのY軸回りの取付角誤差及びZ軸回りの取付角誤差はそれぞれΔθ2x、Δθ3xであり、Y軸角速度センサー10bのZ軸回りの取付角誤差及びX軸回りの取付角誤差はそれぞれΔθ3y、Δθ1yであり、Z軸角速度センサー10cのX軸回りの取付角誤差及びY軸回りの取付角誤差はそれぞれΔθ3z、Δθ1zである。
 加速度センサー50a、50b、50cについても同様に取付角誤差が存在する。そのため、角速度センサー10a、10b、10c、加速度センサー50a、50b、50cの各検出値は理想値とずれている。
 2-2.数学的考察
 従来の補正式(1)、補正式(2)における関数行列式(ヤコビアン)はセンサーの取付角誤差を直接的に反映する補正パラメーターではなく、また、補正式(1)、補正式(2)では前回の検出値を元に関数行列式(ヤコビアン)を用いて今回の検出値を類推するため、検出値に何らかの写像を施すと補正値が得られるようになっていない。そのため、補正式(1)、補正式(2)では補正精度を上げるのに限界がある。そこで、以下では、より精度の高い補正について数学的に考察する。
 3次元ユークリッド空間において、X軸、Y軸、Z軸の回りにそれぞれ角度θの回転を施す回転行列T、T、Tは式(3)によって与えられる。
Figure JPOXMLDOC01-appb-M000003
 そして、3次元ユークリッド空間における任意の回転は回転行列T、T、Tの積の組み合わせにより表すことができる。例えば、Z軸の回りに角度θだけ回転させ、Y軸の回りに角度θだけ回転させ、X軸の回りに角度θだけ回転させることにより、XYZ座標系をX’Y’Z’座標系に変換する行列Tδは式(4)によって与えられる。以下では、Tδを「変換行列」ということにする。
Figure JPOXMLDOC01-appb-M000004
 3つの角速度センサー10a、10b、10cを、検出軸がそれぞれX軸、Y軸、Z軸と平行になるように取り付けたとき、取付角誤差のために実際には検出軸がそれぞれX’軸、Y’軸、Z’軸と平行になるように取り付けられていると仮定する。この場合、角速度センサー10a、10b、10cの各検出値G’、G’、G’と理想値G、G、Gの間には変換行列Tδによる関係式(5)が成立する。
Figure JPOXMLDOC01-appb-M000005
 従って、次の式(6)により、角速度センサー10a、10b、10cの各検出値G’、G’、G’から理想値G、G、Gを計算することができる。
Figure JPOXMLDOC01-appb-M000006
 すなわち、何らかの方法でTδ -1を得ることができれば、式(6)を用いて角速度センサー10a、10b、10cの各検出値を理想値に補正することができる。以下では、Tδ -1を「補正行列」ということにする。
 角速度センサー10a、10b、10cの取付角を光学観測できる場合は、直接θ、θ、θを導出して式(3)を用いて回転行列T、T、Tを計算し、その逆行列T -1、T -1、T -1を用いて補正行列Tδ -1を得ることができる。
 一方、光学観測ができない場合は、例えば、X軸、Y軸、Z軸を中心に角速度センサー10a、10b、10cを回転させる等、X軸、Y軸、Z軸に関する取付角誤差が検出値に反映されるような3通りの入力条件を選択し、その入力条件に対する角速度センサー10a、10b、10cの検出値G’、G’、G’及び理想値G、G、Gをそれぞれ式(6)に代入することにより得られる3つの連立方程式を解けばθ、θ、θを導出することができる。しかし、この連立方程式は非常に複雑であるためθ、θ、θを簡単に導出することができない。
 一方、仮にθが非常に小さい値であれば、次の式(7)が成り立つ。
Figure JPOXMLDOC01-appb-M000007
 従って、取付角誤差θ、θ、θが非常に小さい値であれば、変換行列Tδは次の式(8)のように表される。
Figure JPOXMLDOC01-appb-M000008
 従って、次の式(9)に示すように、変換行列Tδは3つの基底となる行列J、J、Jの線形和で表すことができる。
Figure JPOXMLDOC01-appb-M000009
 前記のように、X軸角速度センサー10aは、Y軸回りに微小角Δθ2x、Z軸回りに微小角Δθ3x回ったX’軸が検出軸になるように取り付けられている。X軸回りの取付角誤差Δθ1x=0であるので、式(9)より、変換行列Tδxは次の式(10)のように表される。
Figure JPOXMLDOC01-appb-M000010
 同様に、Y軸角速度センサー10bは、X軸回りに微小角Δθ1y、Z軸回りに微小角Δθ3y回ったY’軸が検出軸になるように取り付けられている。Y軸回りの取付角誤差Δθ2y=0であるので、式(9)より、変換行列Tδyは次の式(11)のように表される。
Figure JPOXMLDOC01-appb-M000011
 同様に、Z軸角速度センサー10cは、X軸回りに微小角Δθ1z、Y軸回りに微小角Δθ2z回ったZ’軸が検出軸になるように取り付けられている。Z軸回りの取付角誤差Δθ3z=0であるので、式(9)より、変換行列Tδzは次の式(12)のように表される。
Figure JPOXMLDOC01-appb-M000012
 式(10)、(11)、(12)によれば、何らかの方法でΔθ2x、Δθ3x、Δθ1y、Δθ3y、Δθ1z、Δθ2zを得ることができれば、変換行列Tδx、Tδy、Tδzを計算することができる。そして、変換行列Tδx、Tδy、Tδzの逆行列を計算すれば、補正行列Tδx -1、Tδy -1、Tδz -1を作成することができる。そうすると、次の式(13)より、X軸角速度センサー10a、Y軸角速度センサー10b、Z軸角速度センサー10cの各検出値G’、G’、G’をそれぞれ理想値G、G、Gに補正することができる。
Figure JPOXMLDOC01-appb-M000013
 補正式(13)における補正行列Tδx -1、Tδy -1、Tδz -1は、角速度センサー10a、10b、10cの各取付角誤差Δθ2x、Δθ3x、Δθ1y、Δθ3y、Δθ1z、Δθ2zを直接的に反映している。また、補正式(13)によれば、前回の検出値を必要とせず、今回の検出値G’、G’、G’が得られれば直ちに補正値(理想値)G、G、Gを計算することができる。従って、補正式(13)によれば、補正精度の向上及び補正計算処理の高速化を実現することができる。
 なお、補正行列Tδx -1、Tδy -1、Tδz -1は、それぞれ本発明における第1の補正行列、第2の補正行列及び第3の補正行列に相当する。
 次に、Δθ2x、Δθ3x、Δθ1y、Δθ3y、Δθ1z、Δθ2zを得るための方法について説明する。
 X軸角速度センサー10aをX軸回りに角度Δθxx、Y軸回りに角度Δθxy、Z軸回りに角度Δθxzだけ移動させたときに、X’軸回りに角度Δθxx’、Y’軸回りに角度Δθxy’、Z’軸回りに角度Δθxz’だけ移動することを示す検出値が得られる関係にあるとすると、式(10)より次の式(14)が成立する。
Figure JPOXMLDOC01-appb-M000014
 同様に、Y軸角速度センサー10bをX軸回りに角度Δθyx、Y軸回りに角度Δθyy、Z軸回りに角度Δθyzだけ移動させたときに、X’軸回りに角度Δθyx’、Y’軸回りに角度Δθyy’、Z’軸回りに角度Δθyz’だけ移動することを示す検出値が得られる関係にあるとすると、式(11)より次の式(15)が成立する。
Figure JPOXMLDOC01-appb-M000015
 同様に、Z軸角速度センサー10cをX軸回りに角度Δθzx、Y軸回りに角度Δθzy、Z軸回りに角度Δθzzだけ移動させたときに、X’軸回りに角度Δθzx’、Y’軸回りに角度Δθzy’、Z’軸回りに角度Δθzz’だけ移動することを示す検出値が得られる関係にあるとすると、式(12)より次の式(16)が成立する。
Figure JPOXMLDOC01-appb-M000016
 まず、姿勢検出装置1をY軸及びZ軸回りには回転させずにX軸回りにのみ角度Δθだけ回転させた場合、式(15)においてΔθyx=Δθ、Δθyy=0、Δθyz=0なので、Δθyx’=Δθ、Δθyz’=0とともに以下の関係式(17)を得ることができる。
Figure JPOXMLDOC01-appb-M000017
 Y軸角速度センサー10bの検出値(Y’軸回りの角速度)に所定時間を掛けることによりΔθyy’を得ることができるので、式(17)にΔθyy’とΔθを代入することにより、Δθ3yを得ることができる。
 同様に、式(16)においてΔθzx=Δθ、Δθzy=0、Δθzz=0なので、Δθzx’=Δθ、Δθzy’=0とともに以下の関係式(18)を得ることができる。
Figure JPOXMLDOC01-appb-M000018
 Z軸角速度センサー10cの検出値(Z’軸回りの角速度)に所定時間を掛けることによりΔθzz’を得ることができるので、式(18)にΔθzz’とΔθを代入することにより、Δθ2zを得ることができる。
 次に、姿勢検出装置1をX軸及びZ軸回りには回転させずにY軸回りにのみ角度Δθだけ回転させた場合、式(14)においてΔθxx=0、Δθxy=Δθ、Δθxz=0なので、Δθxy’=Δθ、Δθxz’=0とともに以下の関係式(19)を得ることができる。
Figure JPOXMLDOC01-appb-M000019
 X軸角速度センサー10aの検出値(X’軸回りの角速度)に所定時間を掛けることによりΔθxx’を得ることができるので、式(19)にΔθxx’とΔθを代入することにより、Δθ3xを得ることができる。
 同様に、式(16)においてΔθzx=0、Δθzy=Δθ、Δθzz=0なので、Δθzx’=0、Δθzy’=Δθとともに以下の関係式(20)を得ることができる。
Figure JPOXMLDOC01-appb-M000020
 Z軸角速度センサー10cの検出値(Z’軸回りの角速度)に所定時間を掛けることによりΔθzz’を得ることができるので、式(20)にΔθzz’とΔθを代入することにより、Δθ1zを得ることができる。
 最後に、姿勢検出装置1をX軸及びY軸回りには回転させずにZ軸回りにのみ角度Δθだけ回転させた場合、式(14)においてΔθxx=0、Δθxy=0、Δθxz=Δθなので、Δθxy’=0、Δθxz’=Δθとともに以下の関係式(21)を得ることができる。
Figure JPOXMLDOC01-appb-M000021
 X軸角速度センサー10aの検出値(X’軸回りの角速度)に所定時間を掛けることによりΔθxx’を得ることができるので、式(21)にΔθxx’とΔθを代入することにより、Δθ2xを得ることができる。
 同様に、式(15)においてΔθyx=0、Δθyy=0、Δθyz=Δθなので、Δθyx’=0、Δθyz’=Δθとともに以下の関係式(22)を得ることができる。
Figure JPOXMLDOC01-appb-M000022
 Y軸角速度センサー10bの検出値(Y’軸回りの角速度)に所定時間を掛けることによりΔθyy’を得ることができるので、式(22)にΔθyy’とΔθを代入することにより、Δθ1yを得ることができる。
 以上のようにして得られたΔθ2x、Δθ3x、Δθ1y、Δθ3y、Δθ1z、Δθ2zから計算される逆行列Tδx -1、Tδy -1、Tδz -1を式(13)に代入すれば、角速度センサー10a、10b、10cの各検出値G’、G’、G’をそれぞれ理想値G、G、Gに補正することができる。
 なお、実際には、式(13)の補正計算はCPU又は専用回路によりデジタル値で行われる。そのため、角速度センサー10a、10b、10cの検出値G’、G’、G’のA/D変換値にA/D変換のサンプル周期Δtを掛けて得られるX’軸、Y’軸、Z’軸回りの微小回転角Δθx’、Δθy’、Δθz’が、次の式(23)により、X軸、Y軸、Z軸回りの微小回転角Δθ、Δθ、Δθに補正される。
Figure JPOXMLDOC01-appb-M000023
 X軸加速度センサー50a、Y軸加速度センサー50b、Z軸加速度センサー50cの各検出値の補正においても、同様の理論が成り立つ。
 X軸加速度センサー50aは、Y軸回りに微小角Δφ2x、Z軸回りに微小角Δφ3x回ったX’軸が検出軸になるように取り付けられており(Δφ2x、Δφ3xが取付角誤差になる)、変換行列Tγxは次の式(24)のように表される。
Figure JPOXMLDOC01-appb-M000024
 同様に、Y軸加速度センサー50bは、X軸回りに微小角Δφ1y、Z軸回りに微小角Δφ3y回ったY’軸が検出軸になるように取り付けられており(Δφ1y、Δφ3yが取付角誤差になる)、変換行列Tγyは次の式(25)のように表される。
Figure JPOXMLDOC01-appb-M000025
 同様に、Z軸加速度センサー50cは、X軸回りに微小角Δφ1z、Y軸回りに微小角Δφ2z回ったZ’軸が検出軸になるように取り付けられており(Δφ1z、Δφ2zが取付角誤差になる)、変換行列Tγzは次の式(26)のように表される。
Figure JPOXMLDOC01-appb-M000026
 式(24)~式(26)は角速度センサー10a、10b、10cにおける式(10)~式(12)と対応している。
 X軸加速度センサー50aをX軸方向に速度Δvxx、Y軸方向に速度Δvxy、Z軸方向に速度Δvxzで移動させたときに、X’軸方向に速度Δvxx’、Y’軸方向に速度Δvxy’、Z’軸方向に速度Δvxz’で移動することを示す検出値が得られる関係にあるとすると、次の式(27)が成立する。
Figure JPOXMLDOC01-appb-M000027
 Y軸加速度センサー50bをX軸方向に速度Δvyx、Y軸方向に速度Δvyy、Z軸方向に速度Δvyzで移動させたときに、X’軸方向に速度Δvyx’、Y’軸方向に速度Δvyy’、Z’軸方向に速度Δvyz’で移動することを示す検出値が得られる関係にあるとすると、次の式(28)が成立する。
Figure JPOXMLDOC01-appb-M000028
 Z軸加速度センサー50cをX軸方向に速度Δvzx、Y軸方向に速度Δvzy、Z軸方向に速度Δvzzで移動させたときに、X’軸方向に速度Δvzx’、Y’軸方向に速度Δvzy’、Z’軸方向に速度Δvzz’で移動することを示す検出値が得られる関係にあるとすると、次の式(29)が成立する。
Figure JPOXMLDOC01-appb-M000029
 式(27)~式(29)は角速度センサー10a、10b、10cにおける式(14)~式(16)と対応している。そして、角速度センサー10a、10b、10cにおける式(17)~式(22)の導出と同様の方法により以下の(30)~式(35)が得られる。
Figure JPOXMLDOC01-appb-M000030
Figure JPOXMLDOC01-appb-M000031
Figure JPOXMLDOC01-appb-M000032
Figure JPOXMLDOC01-appb-M000033
Figure JPOXMLDOC01-appb-M000034
Figure JPOXMLDOC01-appb-M000035
 そして、Δφ2x、Δφ3x、Δφ1y、Δφ3y、Δφ1z、Δφ2zから計算される逆行列Tγx -1、Tγy -1、Tγz -1を次の式(36)に代入すれば、加速度センサー50a、50b、50cの各検出値A’、A’、A’をそれぞれ理想値A、A、Aに補正することができる。なお、式(36)は角速度センサー10a、10b、10cにおける式(13)と対応している。
Figure JPOXMLDOC01-appb-M000036
 なお、実際には、式(36)の補正計算はCPU又は専用回路によりデジタル値で行われる。そのため、加速度センサー50a、50b、50cの検出値A’、A’、A’のA/D変換値にA/D変換のサンプル周期Δtを掛けて得られるX’軸、Y’軸、Z’軸方向の微小速度Δvx’、Δvy’、Δvz’が、次の式(37)により、X軸、Y軸、Z軸方向の微小速度Δv、Δv、Δvに補正される。
Figure JPOXMLDOC01-appb-M000037
 2-3.補正パラメーター作成用装置
 図8は、本実施形態の補正パラメーター作成用装置の構成を示す図である。
 補正パラメーター作成用装置200は、姿勢検出装置1に含まれる各センサーの取付角誤差に起因する誤差を含む検出値を理想値に補正するための補正パラメーター(補正行列)を作成するために使用される。
 補正パラメーター作成用装置200は、立方体治具210、ソケット220、回転板230、回転モーター240、支持台250、ケーブル260等により構成される。
 立方体治具210は、金属等の材料により立法体形状(直方体形状でもよい)に成型され、3つの面211、212、213が互いに直交するように厳密に面取りがされており、面211にはソケット220が固着されている。立方体治具210の3つの面211、212、213は、それぞれ本発明における治具の第1の面、第2の面及び第3の面に相当する。
 ソケット220は、ソケット本体222と開閉可能な蓋体224で構成されており、ソケット本体220は姿勢検出装置1を所定の向きに隙間無く収容できるようになっている。
 立方体治具210は、姿勢検出装置1をソケット220にセットすることにより、X軸、Y軸及びZ軸がそれぞれ面212、213及び211と垂直になるように、姿勢検出装置1を固定可能になっている。また、立方体治具210の面211、212、213とそれぞれ対向する3つの面214、215、216にはそれぞれ固定金具(図示せず)が設けられている。
 回転板230は、上面231の凹凸が無視できるほど小さく、上面231には固定金具(図示せず)が設けられており、立方体治具210のいずれかの固定金具を回転板230の固定金具に接合することにより、立法体治具210の面214、215、216のいずれかを上面231に固定可能になっている。
 また、回転板230は傾きを調整することができるようになっており、データ補正装置200が設置された状態で回転板230の上面231が水平になるように厳密に調整される。
 回転モーター240は、支持台250に取り付けられており、鉛直方向を軸として時計回り又は半時計回りに所定範囲の角速度で回転できるようになっている。
 ケーブル260は回転モーター240の制御回路(図示せず)に接続されている。ケーブル260にはパーソナルコンピュータ等の制御装置(図示せず)が接続され、GPIB(General Purpose Interface Bus)等のインターフェースにより回転モーター250の回転速度が調整できるようになっている。
 なお、回転モーター240は、本発明における回転制御部として機能する。
 本実施形態によれば、立方体治具210と回転板230を用いることにより、図13A等に示した回転腕530を必要としないので、よりコンパクトかつ低コストの補正パラメーター作成用装置200を提供することができる。本実施形態に係る補正パラメーター作成用装置200を用いることにより、後述するように、補正行列Tδx -1、Tδy -1、Tδz -1、Tγx -1、Tγy -1、Tγz -1を簡単に短時間で取得することができる。
 2-4.補正パラメーター作成手順
 次に、図8に示した補正パラメーター作成用装置200を用いて、補正パラメーター(補正行列)を作成する手順の一例について説明する。
 図9は、本実施形態における補正パラメーター作成手順の一例を示すフローチャート図である。
 まず、回転板230を上面231が水平になるように設置する(ステップS10)。
 次に、立方体冶具210の面211に取り付けられているソケット220に姿勢検出装置1をセットする(ステップS12)。
 次に、立方体治具210の面212と対向する面215を回転板230の上面231に固定する(ステップS14)。これにより、補正パラメーター作成用装置200は図10Aのようにセッティングされ、X軸の正方向が鉛直上向きになるように、立方体冶具210が回転板230に固定される。
 次に、回転板230を静止させた状態でY軸加速度センサー50b、Z軸加速度センサー50cの各検出値を取得し、式(30)、式(31)よりΔφ3y、Δφ2zを計算する(ステップS16)。具体的には、加速度センサー50b、50cの各検出値A’、A’をサンプリングしてA/D変換した値にサンプル周期Δtを掛けて得られる微小速度Δvy’、Δvz’を計算する。ここで、Δvy’及びΔvz’は、それぞれ式(30)におけるΔvyy’及び式(31)におけるΔvzz’に相当する。また、加速度センサー50b、50cにはX軸の負方向に重力加速度gが加わるため、式(30)、式(31)におけるΔv=-g×Δtである。従って、式(30)、式(31)よりΔφ3y、Δφ2zを計算することができる。
 次に、回転板230を角速度ωで回転させてY軸角速度センサー10b、Z軸角速度センサー10cの各検出値を取得し、式(17)、式(18)よりΔθ3y、Δθ2zを計算する(ステップS18)。具体的には、角速度センサー10b、10cの各検出値G’、G’をサンプリングしてA/D変換した値にサンプル周期Δtを掛けて得られる微小回転角Δθy’、Δθz’を計算する。ここで、Δθy’及びΔθz’は、それぞれ式(17)におけるΔθyy’及び式(18)におけるΔθzz’に相当する。また、角速度センサー10b、10cはX軸回りに角速度ωで回転するため、式(17)、式(28)におけるΔθ=ω×Δtである。従って、式(17)、式(18)よりΔθ3y、Δθ2zを計算することができる。
 次に、立方体治具210の面213と対向する面216を回転板230の上面231に固定する(ステップS20)。これにより、補正パラメーター作成用装置200は図10Bのようにセッティングされ、Y軸の正方向が鉛直上向きになるように、立方体冶具210が回転板230に固定される。
 次に、回転板230を静止させた状態でX軸加速度センサー50a、Z軸加速度センサー50cの各検出値を取得し、式(32)、式(33)よりΔφ3x、Δφ1zを計算する(ステップS22)。ステップS22における具体的な処理は、ステップS16と同様であるため説明を省略する。
 次に、回転板230を角速度ωで回転させてX軸角速度センサー10a、Z軸角速度センサー10cの各検出値を取得し、式(19)、式(20)よりΔθ3x、Δθ1zを計算する(ステップS24)。ステップS24における具体的な処理は、ステップS18と同様であるため説明を省略する。
 次に、立方体治具210の面211と対向する面214を回転板230の上面231に固定する(ステップS26)。これにより、補正パラメーター作成用装置200は図10Cのようにセッティングされ、Z軸の正方向が鉛直上向きになるように、立方体冶具210が回転板230に固定される。
 次に、回転板230を静止させた状態でX軸加速度センサー50a、Y軸加速度センサー50bの各検出値を取得し、式(34)、式(35)よりΔφ2x、Δφ1yを計算する(ステップS28)。ステップS28における具体的な処理は、ステップS16と同様であるため説明を省略する。
 次に、回転板230を角速度ωで回転させてX軸角速度センサー10a、Y軸角速度センサー10bの各検出値を取得し、式(21)、式(22)よりΔθ2x、Δθ1yを計算する(ステップS30)。ステップS30における具体的な処理は、ステップS18と同様であるため説明を省略する。
 最後に、補正行列Tδx -1、Tδy -1、Tδz -1、Tγx -1、Tγy -1、Tγz -1を作成する(ステップS32)。具体的には、ステップS24、S30でそれぞれ計算したΔθ3x、Δθ2xを式(10)に代入して得られる変換行列Tδxの逆行列を計算することにより補正行列Tδx -1を作成することができる。同様に、ステップS18、S30でそれぞれ計算したΔθ3y、Δθ1yを式(11)に代入して得られる変換行列Tδyの逆行列を計算することにより補正行列Tδy -1を作成することができる。同様に、ステップS18、S24でそれぞれ計算したΔθ2z、Δθ1zを式(12)に代入して得られる変換行列Tδzの逆行列を計算することにより補正行列Tδz -1を作成することができる。同様に、ステップS22、S28でそれぞれ計算したΔφ3x、Δφ2xを式(24)に代入して得られる変換行列Tγxの逆行列を計算することにより補正行列Tγx -1を作成することができる。同様に、ステップS16、S28でそれぞれ計算したΔφ3y、Δφ1yを式(25)に代入して得られる変換行列Tγyの逆行列を計算することにより補正行列Tγy -1を作成することができる。同様に、ステップS16、S22でそれぞれ計算したΔφ2z、Δφ1zを式(26)に代入して得られる変換行列Tγzの逆行列を計算することにより補正行列Tγz -1を作成することができる。
 なお、以上の処理は、補正パラメーター作成用装置200のケーブル260に接続されたパーソナルコンピュータ等により行われる。そして、本実施形態を用いて作成した補正パラメーターは、例えば、姿勢検出装置1の後段に接続されるユーザー側のマイコンに実装される補正計算処理用のタスクにおいて使用される。
 本実施形態によれば、立方体治具210を用いることにより、X軸、Y軸、Z軸が立方体治具210のそれぞれ面212、面213、面211と垂直になるように、面211に姿勢検出装置1を容易に固定することができる。そして、上面231が水平になるように回転板230を設置すれば、立方体治具210の面215、216、214をそれぞれ回転板230の上面231に固定するだけで、簡単に、X軸、Y軸、Z軸をそれぞれ鉛直方向と平行にすることができる。さらに、X軸、Y軸、Z軸をそれぞれ鉛直方向と平行にした状態において、回転板230を静止させることにより加速度センサー50a、50b、50cの各検出値を、回転板230を回転させることにより角速度センサー10a、10b、10cの各検出値を、それぞれ簡単に短時間で取得することができる。
 すなわち、最初に一度だけ、上面231が水平になるように回転板230を設置すれば回転板230の回転方向が固定されるので、X軸、Y軸、Z軸に関する検出値を取得するためのセッティング時間を大幅に短縮することができる。従って、本実施形態によれば、補正行列Tδx -1、Tδy -1、Tδz -1、Tγx -1、Tγy -1、Tγz -1をより低コストで作成することができる。
 なお、本実施形態では、例えば、X軸角速度センサー10aの補正行列Tδx -1は、Y軸及びZ軸回りの回転時のX軸角速度センサー10aの検出値に基づいて作成される。そのため、各検出軸に対する取付角誤差と他軸感度誤差の両方が考慮された補正パラメーターが作成される。
 本実施形態を用いて作成した補正パラメーターは、移動体やロボットの姿勢検出、姿勢制御を行う装置、バーチャルリアリティー等に使用されるヘッドマウントディスプレイ、頭の姿勢角度を検出するトラッカー、3Dゲームパッド等を使用するゲーム機、デジタルカメラ、携帯電話機、携帯型情報端末、カーナビゲーションシステム等の種々の電子機器に組み込まれる姿勢検出装置の検出値の補正に使用することができる。
 3.補正機能付き姿勢検出装置
 図11は、本実施形態の姿勢検出装置の構成を示す図である。
 補正機能付き姿勢検出装置300は、角速度センサーモジュール2、加速度センサーモジュール3、アンチエリアスフィルター310a、310b、310c、350a、350b、350c、A/D変換回路320a、320b、320c、360a、360b、360c、補正計算処理部370、記憶部380を含んで構成されている。
 角速度センサーモジュール2及び加速度センサーモジュール3は図1及び図2と同じ構成であるため、その説明を省略する。
 アンチエリアスフィルター310a、310b、310c、350a、350b、350cは、それぞれA/D変換回路320a、320b、320c、360a、360b、360cの前段に配置され、角速度検出信号38a、38b、38c及び加速度検出信号78a、78b、78cに対して、それぞれA/D変換回路320a、320b、320c、360a、360b、360cのサンプリングによりDC付近の周波数帯域に折り返すノイズをあらかじめ無視できる程度にまで減衰させる。
 なお、アンチエリアスフィルター310a、310b、310c、360a、360b、360cは、例えば、スイッチトキャパシタフィルター(Switched Capacitor Filter(SCF))として構成することができる。
 A/D変換回路320a、320b、320c、360a、360b、360cは、角速度検出信号38a、38b、38c及び加速度検出信号78a、78b、78cがアンチエリアスフィルター310a、310b、310c、350a、350b、350cによってそれぞれフィルター処理された信号をそれぞれ所定ビット数の角速度検出信号322a、322b、322c及び加速度検出信号362a、362b、362cに変換する。A/D変換回路320a、320b、320c、360a、360b、360cは、本発明におけるA/D変換処理部として機能し、フラッシュ型(並列比較型)、パイプライン型、逐次比較型、デルタシグマ方式等の既知の様々なタイプのAD変換回路により構成することができる。
 記憶部380には、角速度センサーの補正パラメーター382と加速度センサーの補正用パラメーター384が記憶されている。具体的には、補正パラメーター382は補正行列Tδx -1、Tδy -1、Tδz -1であり、補正パラメーター384は補正行列Tγx -1、Tγy -1、Tγz -1である。
 補正計算処理部370は、角速度検出信号322a、322b、322cと補正パラメーター382に基づいて補正式(23)を計算することにより、角速度センサー10a、10b、10cの取付角誤差に起因する角速度検出信号38a、38b、38cの誤差が補正された角速度検出信号302a、302b、302cを生成する。具体的には、補正計算処理部370は、角速度検出信号322a、322b、322cのデジタル値にA/D変換のサンプル周期Δtを掛けた値を補正式(23)の微小回転角Δθx’、Δθy’、Δθz’にそれぞれ代入して微小回転角Δθ、Δθ、Δθを計算し、さらに微小回転角Δθ、Δθ、ΔθをΔtで割ったデジタル値に対応する角速度検出信号302a、302b、302cを生成する。
 同様に、補正計算処理部370は、加速度検出信号362a、362b、362cと補正パラメーター384に基づいて補正式(37)を計算することにより、加速度センサー50a、50b、50cの取付角誤差に起因する加速度検出信号78a、78b、78cの誤差が補正された角速度検出信号304a、304b、304cを生成する。具体的には、補正計算処理部370は、加速度検出信号362a、362b、362cのデジタル値にA/D変換のサンプル周期Δtを掛けた値を補正式(37)の微小速度Δvx’、Δvy’、Δvz’にそれぞれ代入して微小速度Δv、Δv、Δvを計算し、さらに微小速度Δv、Δv、ΔvをΔtで割ったデジタル値に対応する加速度検出信号304a、304b、304cを生成する。
 補正計算処理部370は補正計算処理を行う専用回路として実現することもできるし、CPU(Central Processing Unit)が記憶部380等に記憶されたプログラムを実行することにより補正計算処理部370の機能を実現することもできる。
 図12は、本実施形態の姿勢検出装置の他の構成を示す図である。
 図12において、角速度センサーモジュール2、加速度センサーモジュール3、アンチエリアスフィルター310a、310b、310c、350a、350b、350c、記憶部380の構成は図11と同様であるため、その説明を省略する。
 マルチプレクサ390は、角速度検出信号38a、38b、38c及び加速度検出信号78a、78b、78cがアンチエリアスフィルター310a、310b、310c、350a、350b、350cによってそれぞれフィルター処理された信号を所定周期で時分割に順次選択する。
 A/D変換回路320は、マルチプレクサ390により選択された信号を所定ビット数の検出信号322に変換する。A/D変換回路320及びマルチプレクサ390は、それぞれ本発明におけるA/D変換処理部及び信号選択処理部として機能する。
 補正計算処理部370は、所定周期毎に検出信号322をサンプリングし、検出信号322が角速度検出信号38a、38b、38cに対応する場合は、検出信号322と補正パラメーター382に基づいて補正式(23)を計算することにより補正された角速度検出信号をそれぞれ生成し、検出信号322として時分割に出力する。
 同様に、補正計算処理部370は、検出信号322が加速度検出信号78a、78b、78cに対応する場合は、検出信号322と補正パラメーター384に基づいて補正式(37)を計算することにより補正された加速度検出信号をそれぞれ生成し、検出信号322として時分割に出力する。
 なお、図11又、図12に示した姿勢検出装置300において、補正計算処理部370をバイパスしてA/D変換回路320a等の出力を外部に出力可能なバイパスモードを設ければ、バイパスモードに設定した状態で本実施形態の補正パラメーター作成方法を用いて補正パラメーター382、384を作成することができる。すなわち、姿勢検出装置300も本実施形態の補正パラメーター作成方法の適用対象になり得る。
 本実施形態によれば、補正計算処理部370は、前記の通り、補正精度の向上及び補正計算処理の高速化を実現可能な補正式(23)及び補正式(37)に従い補正値を計算するので、補正精度がより高く、かつ、補正計算処理がより速い姿勢検出装置を実現することができる。
 また、本実施形態によれば、姿勢検出装置300の後段に接続されるユーザー側のマイコンに補正計算処理を実装する必要がなくなるため、タスクカプセル化の見地からユーザーに受け入れられやすい。
 また、本実施形態の姿勢検出装置300はセンサー検出信号をデジタル化して出力するので、姿勢検出装置300とユーザー側のマイコンの間にA/D変換回路を接続する必要がなくなる。
 本実施形態の姿勢検出装置300は、移動体やロボットの姿勢検出、姿勢制御を行う装置、バーチャルリアリティー等に使用されるヘッドマウントディスプレイ、頭の姿勢角度を検出するトラッカー、3Dゲームパッド等を使用するゲーム機、デジタルカメラ、携帯電話機、携帯型情報端末、カーナビゲーションシステム等の種々の電子機器に組み込むことができる。
 なお、本発明は本実施形態に限定されず、本発明の要旨の範囲内で種々の変形実施が可能である。
 例えば、図1に示した姿勢検出装置は、3つの角速度センサー10a、10b、10cと3つの加速度センサー50a、50b、50cを含んでいるが、本発明の補正パラメーター作成方法が適用される姿勢検出装置はこれに限られない。すなわち、本発明の補正パラメーター作成方法が適用される姿勢検出装置は直交する3軸の角速度又は加速度を検出可能な構成であればよく、例えば、角速度センサー10a、10b、10cのみを含む姿勢検出装置、加速度センサー50a、50b、50cのみを含む姿勢検出装置、角速度センサー10a、10bと加速度センサー50cのみを含む姿勢検出装置、角速度センサー10aと加速度センサー50b、50cのみを含む姿勢検出装置等も適用対象に含まれる。
 また、例えば、図8に示した補正パラメーター作成用装置200は、面211にソケット220が1つだけ取り付けられているが、面211に複数のソケット220を取り付けるようにしてもよい。そして、複数のソケット220の各々に姿勢検出装置1をセットすることにより、各姿勢検出装置1の検出値を同時に取得することができる。
 また、例えば、図9の手順では、X軸、Y軸、Z軸の順に姿勢検出装置1の検出値を取得するようになっているが、任意の軸の順に姿勢検出装置1の検出値を取得することができる。
 本発明は、実施の形態で説明した構成と実質的に同一の構成(例えば、機能、方法及び結果が同一の構成、あるいは目的及び効果が同一の構成)を含む。また、本発明は、実施の形態で説明した構成の本質的でない部分を置き換えた構成を含む。また、本発明は、実施の形態で説明した構成と同一の作用効果を奏する構成又は同一の目的を達成することができる構成を含む。また、本発明は、実施の形態で説明した構成に公知技術を付加した構成を含む。
1 姿勢検出装置、2 角速度センサーモジュール、3 加速度センサーモジュール、4 パッケージ、5a,5b,5c パッケージの面、10a,10b,10c 角速度センサー、11a,11b,11c 振動子、12a,12b,12c 駆動電極、13a,13b,13c 駆動電極、14a,14b,14c 検出電極、15a,15b,15c 検出電極、20a,20b,20c 駆動回路、21a 電流電圧変換器(I/V変換器)、22a AC増幅器、23a 自動利得制御回路(AGC) 24a コンパレータ、30a,30b,30c 検出回路、31a,32a チャージアンプ、33a 差動増幅器、34a AC増幅器、35a 同期検波回路、36a DC増幅器、37a 積分回路(LPF)、38a,38b,38c 検出信号、41a 駆動振動腕、42a 検出振動腕、43a 幅広部、44a 駆動用基部、45a 連結腕、46a 幅広部、47a 検出用基部、50a,50b,50c 加速度センサー、51a,51b,51c 振動子、52a,52b,52c 駆動電極、53a,53b,53c 駆動電極、54a,54b,54c 検出電極、55a,55b,55c 検出電極、56a,56b,56c 基端部、57a,57b,57c 基端部、58a,58b,58c 駆動振動腕、60a,60b,60c 駆動回路、70a,70b,70c 検出回路、78a,78b,78c 検出信号、80 絶縁基板、82a,82b,82c パッケージ、84a,84b,84c パッケージ本体、86a,86b,86c 蓋体、90 ベース、91,92,93 素子取付面、100 ウエイト、101,102,103 素子接合面、200 補正パラメーター作成用装置、210 立方体治具、211,212,213,214,215,216 立方体治具の面、220 ソケット、222 ソケット本体、224 蓋体、230 回転板、231 回転板の上面、240 回転モーター、250 支持台、260 ケーブル、300 姿勢検出装置、302,302a,302b,302c 検出信号、304a,304b,304c 検出信号、310a,310b,310c アンチエリアスフィルター、310a,310b,310c A/D変換回路、320 A/D変換回路、322a,322b,322c 検出信号、350a,350b,350c アンチエリアスフィルター、360a,360b,360c A/D変換回路、362a,362b,362c 検出信号、370 補正計算処理部、380 記憶部、382 補正パラメーター、384 補正パラメーター、390 マルチプレクサ、500 補正パラメーター作成用装置、510 テーブル、520 ソケット、530 回転腕

Claims (8)

  1.  検出軸が互いに直交する第1の軸、第2の軸及び第3の軸とそれぞれ略平行になるように取り付けられ角速度又は加速度を検出する第1のセンサー、第2のセンサー及び第3のセンサーを含み、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの検出信号に基づいて物体の姿勢を検出する姿勢検出装置の検出値を、前記第1の軸、前記第2の軸及び前記第3の軸を座標軸とする直交座標系における検出値に補正する補正式の補正パラメーターを作成する方法であって、
     上面が水平になるように回転板を設置するステップと、
     互いに直交する第1の面、第2の面、第3の面を有する直方体形状の治具の前記第1の面に、前記第1の軸が前記第2の面に垂直になり、前記第2の軸が前記第3の面に垂直になり、前記第3の軸が前記第1の面に垂直になるように、前記姿勢検出装置を固定するステップと、
     前記治具の前記第2の面と対向する面を前記回転板の前記上面に固定し、前記回転板を静止又は所定の角速度で回転させて前記姿勢検出装置の検出値を取得する第1検出値取得ステップと、
     前記治具の前記第3の面と対向する面を前記回転板の前記上面に固定し、前記回転板を静止又は所定の角速度で回転させて前記姿勢検出装置の検出値を取得する第2検出値取得ステップと、
     前記治具の前記第1の面と対向する面を前記回転板の前記上面に固定し、前記回転板を静止又は所定の角速度で回転させて前記姿勢検出装置の検出値を取得する第3検出値取得ステップと、
     取得した検出値に基づいて、前記補正パラメーターを作成する補正パラメーター作成ステップと、を含むことを特徴とする姿勢検出装置の補正パラメーター作成方法。
  2.  請求項1において、
     前記補正式は、
     前記補正パラメーターとして前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出値を前記直交座標系における各検出値に補正するための第1の補正行列、第2の補正行列及び第3の補正行列を含み、前記第1の補正行列、前記第2の補正行列及び前記第3の補正行列の各々と、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出値がA/D変換されたデジタル値をそれぞれ要素として含む行列の各々と、の積により得られる3つの行列の和として与えられることを特徴とする姿勢検出装置の補正パラメーター作成方法。
  3.  請求項2において、
     前記第1の補正行列、前記第2の補正行列及び前記第3の補正行列は、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出軸をそれぞれ前記第1の軸、前記第2の軸及び前記第3の軸に変換する回転行列の逆行列であることを特徴とする姿勢検出装置の補正パラメーター作成方法。
  4.  請求項2又は3において、
     前記補正パラメーター作成ステップは、
     前記第1検出値取得ステップにおいて取得した前記検出値に基づいて、前記第2のセンサー及び前記第3のセンサーの前記第1の軸回りの各取付角誤差を算出するステップと、
     前記第2検出値取得ステップにおいて取得した前記検出値に基づいて、前記第1のセンサー及び前記第3のセンサーの前記第2の軸回りの各取付角誤差を算出するステップと、
     前記第3検出値取得ステップにおいて取得した前記検出値に基づいて、前記第1のセンサー及び前記第2のセンサーの前記第3の軸回りの各取付角誤差を算出するステップと、
     前記第1のセンサーの前記第2の軸回りの前記取付角誤差及び前記第3の軸回りの前記取付角誤差に基づいて、前記第1の補正行列を作成するステップと、
     前記第2のセンサーの前記第1の軸回りの前記取付角誤差及び前記第3の軸回りの前記取付角誤差に基づいて、前記第2の補正行列を作成するステップと、
     前記第3のセンサーの前記第1の軸回りの前記取付角誤差及び前記第2の軸回りの前記取付角誤差に基づいて、前記第3の補正行列を作成するステップと、を含むことを特徴とする姿勢検出装置の補正パラメーター作成方法。
  5.  検出軸が互いに直交する第1の軸、第2の軸及び第3の軸とそれぞれ略平行になるように取り付けられ角速度又は加速度を検出する第1のセンサー、第2のセンサー及び第3のセンサーを含み、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの検出信号に基づいて物体の姿勢を検出する姿勢検出装置の検出値を、前記第1の軸、前記第2の軸及び前記第3の軸を座標軸とする直交座標系における検出値に補正する補正式の補正パラメーターを作成するために使用される補正パラメーター作成用装置であって、
     互いに直交する第1の面、第2の面、第3の面を有し、前記第1の面に、前記第1の軸が前記第2の面に垂直になり、前記第2の軸が前記第3の面に垂直になり、前記第3の軸が前記第1の面に垂直になるように前記姿勢検出装置を固定可能な直方体形状の治具と、
     上面に、前記治具の前記第1の面、前記第2の面、前記第2の面とそれぞれ対向する面のいずれかを固定可能な回転板と、
     前記回転板を所定の角速度で回転させる回転制御部と、を含むことを特徴とする姿勢検出装置の補正パラメーター作成用装置。
  6.  検出軸が互いに直交する第1の軸、第2の軸及び第3の軸とそれぞれ略平行になるように取り付けられ角速度又は加速度を検出する第1のセンサー、第2のセンサー及び第3のセンサーと、
     前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出値を前記第1の軸、前記第2の軸及び前記第3の軸を座標軸とする直交座標系における検出値に補正する補正式の補正パラメーターが記憶された記憶部と、
     前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出信号をデジタル信号に変換する処理を行うA/D変換処理部と、
     前記デジタル信号の各々と前記補正パラメーターに基づいて前記補正式を計算する処理を行う補正計算処理部と、を含み、
     前記補正式は、
     前記補正パラメーターとして前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出値を前記直交座標系における各検出値に補正するための第1の補正行列、第2の補正行列及び第3の補正行列を含み、前記第1の補正行列、前記第2の補正行列及び前記第3の補正行列の各々と、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出値がA/D変換されたデジタル値をそれぞれ要素として含む行列の各々と、の積により得られる3つの行列の和として与えられることを特徴とする姿勢検出装置。
  7.  請求項6において、
     前記第1の補正行列、前記第2の補正行列及び前記第3の補正行列は、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各検出軸をそれぞれ前記第1の軸、前記第2の軸及び前記第3の軸に変換する回転行列の逆行列であることを特徴とする姿勢検出装置。
  8.  請求項6又は7において、
     所定の周期で、前記第1のセンサー、前記第2のセンサー及び前記第3のセンサーの各前記検出信号のいずれか1つを順次選択する処理を行う信号選択処理部を含み、
     前記A/D変換処理部は、
     前記信号選択処理部が選択した検出値を順次A/D変換処理するA/D変換回路を含むことを特徴とする姿勢検出装置。
PCT/JP2009/069249 2008-11-13 2009-11-12 姿勢検出装置の補正パラメーター作成方法、姿勢検出装置の補正パラメーター作成用装置及び姿勢検出装置 Ceased WO2010055871A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US13/126,446 US20110202300A1 (en) 2008-11-13 2009-11-12 Method for creating correction parameter for posture detecting device, device for creating correction parameter for posture detecting device, and posture detecting device
CN2009801453275A CN102216790B (zh) 2008-11-13 2009-11-12 姿势检测装置的校正参数生成方法、用于生成姿势检测装置的校正参数的装置以及姿势检测装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2008-291213 2008-11-13
JP2008291213A JP2010117260A (ja) 2008-11-13 2008-11-13 姿勢検出装置の補正パラメーター作成方法、姿勢検出装置の補正パラメーター作成用装置及び姿勢検出装置

Publications (1)

Publication Number Publication Date
WO2010055871A1 true WO2010055871A1 (ja) 2010-05-20

Family

ID=42170003

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2009/069249 Ceased WO2010055871A1 (ja) 2008-11-13 2009-11-12 姿勢検出装置の補正パラメーター作成方法、姿勢検出装置の補正パラメーター作成用装置及び姿勢検出装置

Country Status (4)

Country Link
US (1) US20110202300A1 (ja)
JP (1) JP2010117260A (ja)
CN (3) CN103257251A (ja)
WO (1) WO2010055871A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109978991A (zh) * 2019-03-14 2019-07-05 西安交通大学 基于视觉快速实现复杂构件装夹位姿误差在线测量的方法
CN115409767A (zh) * 2021-05-10 2022-11-29 北京芙睿特无限科技发展有限公司 一种包裹体积确定方法和装置
WO2024185856A1 (ja) * 2023-03-07 2024-09-12 住友精密工業株式会社 角速度センサの補正方法、角速度センサの補正プログラム、および、角速度センサシステム

Families Citing this family (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2743502B2 (ja) * 1989-08-04 1998-04-22 神鋼電機株式会社 誘導溶解炉の炉壁冷却機構
CN102607588A (zh) * 2012-01-18 2012-07-25 中国商用飞机有限责任公司 二自由度转台的校准装置
JP6071302B2 (ja) * 2012-07-26 2017-02-01 オリンパス株式会社 キャリブレーション装置およびプログラム
JP6120435B2 (ja) * 2012-11-30 2017-04-26 株式会社レールテック 曲線整正のための修正量算出システム及び修正量算出用コンピュータプログラム
KR101891980B1 (ko) * 2012-12-11 2018-08-27 한국전자통신연구원 회전체의 회전 특성 측정 방법 및 장치
CN103994845B (zh) * 2013-02-19 2019-09-10 精工爱普生株式会社 力检测装置、机械手、以及移动体
CN103292769B (zh) * 2013-06-19 2015-11-25 桂林电子科技大学 一种基于最小区域的平面倾斜度误差检测方法
CN103292770B (zh) * 2013-06-23 2015-12-23 桂林电子科技大学 一种计算圆锥体零件作用尺寸的方法
CN103292771B (zh) * 2013-06-25 2015-12-02 桂林电子科技大学 一种基于最小区域的零件圆锥度的计算方法
CN106456051A (zh) * 2014-11-21 2017-02-22 深圳迈瑞生物医疗电子股份有限公司 一种呼吸监测设备、方法和装置
US9970781B2 (en) * 2015-03-03 2018-05-15 West Virginia University Apparatus for three-axis IMU calibration with a single-axis rate table
US9651399B2 (en) * 2015-03-25 2017-05-16 Northrop Grumman Systems Corporation Continuous calibration of an inertial system
JP6454857B2 (ja) * 2015-05-27 2019-01-23 多摩川精機株式会社 姿勢検出装置及び姿勢検出方法
JP6514089B2 (ja) * 2015-11-02 2019-05-15 株式会社ソニー・インタラクティブエンタテインメント 情報処理装置、情報処理システム、および情報処理方法
JP6693214B2 (ja) * 2016-03-25 2020-05-13 セイコーエプソン株式会社 物理量検出装置、電子機器及び移動体
CN107063237A (zh) * 2016-12-14 2017-08-18 歌尔股份有限公司 一种测量物体姿态角的方法和装置
CN108332741B (zh) * 2017-01-18 2021-09-14 宏达国际电子股份有限公司 定位装置及方法
DE102017108194A1 (de) * 2017-04-18 2018-10-18 Vorwerk & Co. Interholding Gmbh Verfahren zum Betrieb eines sich selbsttätig fortbewegenden Fahrzeugs
DE102017208365A1 (de) * 2017-05-18 2018-11-22 Robert Bosch Gmbh Verfahren zur Orientierungsschätzung eines tragbaren Gerätes
CN109129523B (zh) * 2018-08-30 2021-07-06 燕山大学 基于人机交互的移动机器人实时远程控制系统
CN109697835A (zh) * 2019-01-29 2019-04-30 程家煜 一种监控坐姿和预防近视的方法与装置
JP7251385B2 (ja) * 2019-07-30 2023-04-04 セイコーエプソン株式会社 振動デバイス、電子機器および移動体
CN111486871A (zh) * 2020-04-27 2020-08-04 新石器慧通(北京)科技有限公司 传感器检测方法、装置、检测设备及可读存储介质
JP7571395B2 (ja) * 2020-06-10 2024-10-23 セイコーエプソン株式会社 慣性センサー装置及び慣性センサー装置の製造方法
JP7635550B2 (ja) * 2020-12-28 2025-02-26 セイコーエプソン株式会社 振動整流誤差補正装置、センサーモジュール及び振動整流誤差補正方法
TWI778582B (zh) * 2021-04-15 2022-09-21 桃苗汽車股份有限公司 感知器安裝角度檢查裝置及方法
DE102024118279B3 (de) * 2024-06-27 2025-07-10 Deutsches Zentrum für Luft- und Raumfahrt e.V. Rotationstisch mit verkippbarer Sensoraufnahme

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08304448A (ja) * 1995-05-11 1996-11-22 Nikon Corp 力学量検出器および力学量検出方法
JPH102914A (ja) * 1996-06-13 1998-01-06 Fujikura Ltd 加速度発生装置及びこれを用いた加速度センサ測定装置
JPH10267651A (ja) * 1997-03-28 1998-10-09 Data Tec:Kk 傾斜測定方法及び傾斜測定装置
JP2000338128A (ja) * 1999-03-19 2000-12-08 Ngk Insulators Ltd 加速度センサ素子の感度評価方法
JP2004150900A (ja) * 2002-10-29 2004-05-27 Japan Aviation Electronics Industry Ltd 姿勢角検出装置および取付け角度補正量取得方法
WO2005095998A1 (ja) * 2004-03-31 2005-10-13 National Institute Of Advanced Industrial Science And Technology 加速度を検出するセンサの横感度を計測する方法および加速度計測方法
JP2008128674A (ja) * 2006-11-16 2008-06-05 Eastman Kodak Co 角速度較正方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2353629C (en) * 1998-12-17 2005-12-27 Tokin Corporation Orientation angle detector
JP4853937B2 (ja) * 2003-04-28 2012-01-11 独立行政法人産業技術総合研究所 慣性センサの動的感度マトリックス計測装置およびその計測方法
JP2006176084A (ja) * 2004-12-24 2006-07-06 Advics:Kk 車両挙動センサの検出値補正方法
CN101105503B (zh) * 2007-06-02 2010-10-27 中北大学 捷联式惯导测量组合中加速度计装配误差标量修正方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08304448A (ja) * 1995-05-11 1996-11-22 Nikon Corp 力学量検出器および力学量検出方法
JPH102914A (ja) * 1996-06-13 1998-01-06 Fujikura Ltd 加速度発生装置及びこれを用いた加速度センサ測定装置
JPH10267651A (ja) * 1997-03-28 1998-10-09 Data Tec:Kk 傾斜測定方法及び傾斜測定装置
JP2000338128A (ja) * 1999-03-19 2000-12-08 Ngk Insulators Ltd 加速度センサ素子の感度評価方法
JP2004150900A (ja) * 2002-10-29 2004-05-27 Japan Aviation Electronics Industry Ltd 姿勢角検出装置および取付け角度補正量取得方法
WO2005095998A1 (ja) * 2004-03-31 2005-10-13 National Institute Of Advanced Industrial Science And Technology 加速度を検出するセンサの横感度を計測する方法および加速度計測方法
JP2008128674A (ja) * 2006-11-16 2008-06-05 Eastman Kodak Co 角速度較正方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109978991A (zh) * 2019-03-14 2019-07-05 西安交通大学 基于视觉快速实现复杂构件装夹位姿误差在线测量的方法
CN115409767A (zh) * 2021-05-10 2022-11-29 北京芙睿特无限科技发展有限公司 一种包裹体积确定方法和装置
WO2024185856A1 (ja) * 2023-03-07 2024-09-12 住友精密工業株式会社 角速度センサの補正方法、角速度センサの補正プログラム、および、角速度センサシステム

Also Published As

Publication number Publication date
CN103257251A (zh) 2013-08-21
JP2010117260A (ja) 2010-05-27
US20110202300A1 (en) 2011-08-18
CN102216790A (zh) 2011-10-12
CN103292766A (zh) 2013-09-11
CN102216790B (zh) 2013-06-05

Similar Documents

Publication Publication Date Title
WO2010055871A1 (ja) 姿勢検出装置の補正パラメーター作成方法、姿勢検出装置の補正パラメーター作成用装置及び姿勢検出装置
JP4893335B2 (ja) ジャイロモジュール
JP5682267B2 (ja) 角速度センサ
JP5088540B2 (ja) 検出装置、検出方法及び電子機器
JP2010117371A (ja) 姿勢検出装置
Seeger et al. Development of high-performance high-volume consumer MEMS gyroscopes
CN102192731A (zh) 物理量检测元件、物理量检测装置以及电子设备
JP5135832B2 (ja) 振動型ジャイロセンサ、制御回路、電子機器及び振動型ジャイロセンサの製造方法
JP2016099269A (ja) ジャイロセンサー、電子機器、および移動体
JP4967663B2 (ja) 振動型ジャイロセンサ、制御回路及び電子機器
CN103363982A (zh) 陀螺传感器、电子设备以及移动体
CN102200438A (zh) 物理量检测元件、物理量检测装置及电子设备
JP2017173208A (ja) 回路装置、物理量検出装置、電子機器及び移動体
US8453503B2 (en) Vibrating reed, vibrator, physical quantity sensor, and electronic apparatus
CN102331512B (zh) 二维/三维角速度检测装置、方法、姿态感知设备
JP4893838B2 (ja) ジャイロモジュール
JP2016017792A (ja) センサーモジュールおよび電子機器
JP2016153745A (ja) 回路装置、物理量検出装置、電子機器及び移動体
JP5131491B2 (ja) ジャイロモジュール
JP2015137991A (ja) 機能素子、センサーデバイス、電子機器、および移動体
US20260002782A1 (en) Physical Quantity Detection Device
JP2022049898A (ja) 振動デバイス
WO2021220574A1 (ja) センサモジュール
JP2006292506A (ja) 加速度・角速度検出装置および加速度・角速度検出方法
JP2012021964A (ja) 角速度センサ装置

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200980145327.5

Country of ref document: CN

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 09826118

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 13126446

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 09826118

Country of ref document: EP

Kind code of ref document: A1