WO2006028013A1 - 加速度推定方法、加速度推定装置及び移動端末装置 - Google Patents

加速度推定方法、加速度推定装置及び移動端末装置 Download PDF

Info

Publication number
WO2006028013A1
WO2006028013A1 PCT/JP2005/016112 JP2005016112W WO2006028013A1 WO 2006028013 A1 WO2006028013 A1 WO 2006028013A1 JP 2005016112 W JP2005016112 W JP 2005016112W WO 2006028013 A1 WO2006028013 A1 WO 2006028013A1
Authority
WO
WIPO (PCT)
Prior art keywords
acceleration
component
frequency
rotation
low
Prior art date
Application number
PCT/JP2005/016112
Other languages
English (en)
French (fr)
Inventor
Naomi Nishikata
Hirohisa Kusuda
Original Assignee
Vodafone K.K.
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 Vodafone K.K. filed Critical Vodafone K.K.
Priority to JP2006527204A priority Critical patent/JPWO2006028013A1/ja
Priority to EP05781378A priority patent/EP1788397A4/en
Publication of WO2006028013A1 publication Critical patent/WO2006028013A1/ja
Priority to US11/706,239 priority patent/US8521467B2/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR 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/017Gesture based interaction, e.g. based on a set of recognized hand gestures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F1/00Details not covered by groups G06F3/00 - G06F13/00 and G06F21/00
    • G06F1/16Constructional details or arrangements
    • G06F1/1613Constructional details or arrangements for portable computers
    • G06F1/1626Constructional details or arrangements for portable computers with a single-body enclosure integrating a flat display, e.g. Personal Digital Assistants [PDAs]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F1/00Details not covered by groups G06F3/00 - G06F13/00 and G06F21/00
    • G06F1/16Constructional details or arrangements
    • G06F1/1613Constructional details or arrangements for portable computers
    • G06F1/1633Constructional details or arrangements of portable computers not specific to the type of enclosures covered by groups G06F1/1615 - G06F1/1626
    • G06F1/1684Constructional details or arrangements related to integrated I/O peripherals not covered by groups G06F1/1635 - G06F1/1675
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR 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/03Arrangements for converting the position or the displacement of a member into a coded form
    • G06F3/033Pointing devices displaced or positioned by the user, e.g. mice, trackballs, pens or joysticks; Accessories therefor
    • G06F3/038Control and interface arrangements therefor, e.g. drivers or device-embedded control circuitry
    • G06F3/0383Signal control means within the pointing device
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V30/00Character recognition; Recognising digital ink; Document-oriented image-based pattern recognition
    • G06V30/10Character recognition
    • G06V30/22Character recognition characterised by the type of writing
    • G06V30/228Character recognition characterised by the type of writing of three-dimensional handwriting, e.g. writing in the air
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2200/00Indexing scheme relating to G06F1/04 - G06F1/32
    • G06F2200/16Indexing scheme relating to G06F1/16 - G06F1/18
    • G06F2200/161Indexing scheme relating to constructional details of the monitor
    • G06F2200/1614Image rotation following screen orientation, e.g. switching from landscape to portrait mode
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2200/00Indexing scheme relating to G06F1/04 - G06F1/32
    • G06F2200/16Indexing scheme relating to G06F1/16 - G06F1/18
    • G06F2200/163Indexing scheme relating to constructional details of the computer
    • G06F2200/1637Sensing arrangement for detection of housing movement or orientation, e.g. for controlling scrolling or cursor movement on the display of an handheld computer
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04MTELEPHONIC COMMUNICATION
    • H04M2250/00Details of telephonic subscriber devices
    • H04M2250/12Details of telephonic subscriber devices including a sensor for measuring a physical value, e.g. temperature or motion

Definitions

  • Acceleration estimation method acceleration estimation device, and mobile terminal device
  • the present invention relates to an acceleration estimation method, an acceleration estimation device, and a mobile terminal device. More specifically, the present invention contributes to translation in a moving object, and contributes to translation contribution acceleration and rotation.
  • the present invention relates to an acceleration estimation method for estimating the acceleration, an acceleration estimation device using the acceleration estimation method, and a mobile terminal device including the acceleration estimation device.
  • Patent Document 1 Japanese Unexamined Patent Application Publication No. 2002-169645
  • the above-described conventional technique is excellent as a technique for using a measurement result by a sensor in a mobile phone device or the like.
  • this conventional example when drawing a character trajectory, it is necessary for the user to perform a translation operation on the two-dimensional plane with almost no rotation on the two-dimensional plane. This is difficult for the user to operate.
  • an acceleration sensor mounted on a mobile phone or the like used near the ground surface is usually a static sensor. Gravity acceleration is detected in the stationary state. For this reason, with the mobile phone, etc. in a stationary state, the state in which only gravitational acceleration is detected is the initial state, and by constantly monitoring the change in acceleration, the gravitational direction (vertically downward) is identified and the gravitational acceleration is detected. However, it is difficult to specify the direction of gravity with high accuracy after a long time in the initial state force. For this reason, only the detection result by the acceleration sensor installed in a mobile terminal device with a variable posture, etc., can be used to easily subtract the detection result force gravitational acceleration and accurately derive the acceleration that contributes to the motion. .
  • the present invention has been made in view of the above circumstances, and estimates the translation contribution acceleration contributing to the translation and the rotation contribution acceleration contributing to the rotation of the moving object with improved accuracy. It is an object of the present invention to provide an acceleration estimation method and an acceleration estimation apparatus that can perform the above-described process.
  • the present invention provides a mobile terminal device capable of estimating the translation contribution acceleration contributing to its own translation and the rotation contribution acceleration contributing to the rotation with improved accuracy. With the goal.
  • the first axis uniquely defined in the moving object, the second axis orthogonal to the first axis, and the orthogonal to the first axis and the second axis Within the third axis.
  • This is an acceleration estimation method for estimating the translation contribution acceleration that contributes to the translation of the moving object and the rotation contribution acceleration that contributes to the rotation in the measurement results of acceleration along at least two axial directions.
  • a measurement result collecting step for collecting a measurement result of acceleration along the at least two axes and a measurement result of a posture angle that is a rotation angle of a reference posture force of the moving object; and along the at least two axes;
  • An acceleration frequency separation step for separating a low-frequency acceleration component and a high-frequency acceleration component in a measured acceleration result;
  • a posture angle change rate deriving step for deriving a time change rate of the posture angle;
  • a translational rotation acceleration estimation step for estimating the translation contribution acceleration and the rotation contribution acceleration based on a low-frequency acceleration component and a time change rate of the posture angle. It is a method.
  • the first axis that is uniquely defined in the moving object, the second axis that is orthogonal to the first axis, the third axis that is orthogonal to the first axis and the second axis Measurement results of acceleration along at least two of the axis directions and rotation angle from the reference posture, that is, the measurement result of the posture angle of the moving object are collected.
  • the low frequency acceleration component and the high frequency acceleration component in each acceleration measurement result are separated through the acceleration frequency separation step.
  • Such frequency separation can be achieved, for example, by extracting the low-frequency acceleration component from each of the periodically collected acceleration measurement results by the moving average method and subtracting the low-frequency acceleration component extracted from each acceleration measurement result. This is done by extracting the acceleration component.
  • the posture angle change rate deriving step the time change rate of the posture angle is derived.
  • the high-frequency acceleration component, the low-frequency acceleration component, the posture angle, and the posture are calculated in the translational rotation acceleration estimation process.
  • the translation contribution acceleration and the rotation contribution acceleration are estimated based on the time change rate of the angle.
  • the gravitational acceleration is subtracted from the low-frequency acceleration component using the attitude angle, the magnitude of the rotation contribution acceleration is evaluated by the time change rate of the attitude angle, and the translation contribution acceleration and the general frequency of the rotation contribution acceleration are calculated. Considering characteristics and the like, the translation contribution acceleration and the rotation contribution acceleration are estimated.
  • the acceleration estimation method of the present invention contributes to translation in a moving object.
  • the translation contribution acceleration and the rotation contribution acceleration contributing to the rotation can be estimated with improved accuracy.
  • the translational rotation acceleration estimation step estimates that the high-frequency acceleration component does not include the rotation-contributing acceleration component but includes only the translation-contributing acceleration component.
  • Component estimation step if the absolute value of the time change rate of the posture angle is less than a predetermined value, the subtraction of gravitational acceleration from the low frequency acceleration component does not include the component of the rotation contribution acceleration, When it is estimated that only the translation contribution acceleration component is included, and the absolute value of the time change rate of the posture angle is not less than the predetermined value, the gravitational acceleration is subtracted from the low frequency acceleration component, And a low frequency component estimation step of estimating that only the rotation contribution acceleration component is included but not the translation contribution acceleration component.
  • the high frequency acceleration component does not include the rotation contribution acceleration component, but only the translation contribution acceleration component. It is estimated to be. This is because, in a mobile terminal device such as a mobile phone device, the detection result of the acceleration sensor mounted on the mobile phone device or the like rarely changes suddenly in the case of manual rotation. . In addition, in the case of manual translational motion, the translational contribution acceleration usually has a high frequency component.
  • the low frequency component estimation step next, in the low frequency component estimation step, it is determined whether or not the force is such that the absolute value of the time change rate of the posture angle is less than a predetermined value. If this determination is affirmative, it can be considered that the low-frequency acceleration component force gravitational acceleration minus many translation-contributing accelerations. Therefore, it can be assumed that the low-frequency acceleration component minus the gravitational acceleration does not include the rotation-contributing acceleration component, but only the translation-contributing acceleration component.
  • the low-frequency acceleration component minus the gravitational acceleration includes a lot of rotational contribution acceleration. Therefore, the subtraction of the low-frequency acceleration component force and gravitational acceleration does not include the translation contribution acceleration component, but only the rotation contribution acceleration component.
  • the low-frequency acceleration component contributes to rotation according to the magnitude of the time change rate of the posture angle.
  • the choice between acceleration and translation-contributing acceleration is done with reasonableness. Therefore, the translation contribution acceleration and the rotation contribution acceleration in the moving object can be estimated with improved accuracy.
  • a low frequency component of the time change rate of the posture angle is derived, and the translational acceleration estimation step includes the high frequency acceleration.
  • the component does not include the rotation contribution acceleration component, but only the translation contribution acceleration component, and a high frequency component estimation step; based on the low frequency component of the time change rate of the posture angle, A rotation-contributing low-frequency component estimating step for estimating a rotation-contributing low-frequency acceleration component that is a component of the rotation-contributing acceleration in the low-frequency acceleration component; and the rotation-contributing low-frequency acceleration component and the heavy acceleration from the low-frequency acceleration component; A translation contribution low frequency component estimation step for estimating a translation contribution low frequency acceleration component which is a low frequency component of the translation contribution acceleration. That.
  • a low frequency component of the time change rate of the posture angle is derived.
  • the derivation of the low-frequency component in the time change rate of the posture angle is performed by, for example, using the moving average method for each change amount of the posture angle measurement result collected periodically.
  • the high frequency acceleration component does not include the rotation contribution acceleration component, but includes only the translation contribution acceleration component. It is estimated to be.
  • the rotational contribution low frequency estimation step in the rotational contribution low frequency estimation step, based on the low frequency component of the time change rate of the posture angle, the rotational contribution low in the rotational contribution acceleration component in the low frequency acceleration component.
  • a frequency acceleration component is estimated.
  • the force obtained by subtracting the rotation contribution low frequency acceleration component and the gravitational acceleration from the low frequency acceleration component is the translation contribution low frequency acceleration component, which is the low frequency component of the translation contribution acceleration. It is estimated that.
  • the low-frequency acceleration component is estimated as only the translation contribution acceleration or the rotation contribution acceleration based on the magnitude of the time change rate of the posture angle, It is possible to estimate the translation contribution acceleration contributing to the translation in the moving object and the rotation contribution acceleration contributing to the rotation with improved accuracy.
  • the translational calorie velocity estimating step includes the posture angle change step.
  • a rotation-contributing high-frequency component estimation step for estimating a rotation-contributing high-frequency acceleration component that is a high-frequency component of the rotation-contributing acceleration in the high-frequency acceleration component based on a high-frequency component of the time change rate of the high-frequency acceleration component;
  • a translation-contributing high-frequency component in which the rotation-contributing high-frequency acceleration component is subtracted from the translation-contributing high-frequency acceleration component, which is a high-frequency component of the translation-contributing acceleration, and a low-frequency component of the time change rate of the posture angle On the basis of the rotation offset for estimating the rotation-contributing low-frequency acceleration component that is a component of the rotation-contributing acceleration in the low-frequency acceleration component.
  • a low-frequency component estimation step a value obtained by subtracting the rotation-contributing low-frequency acceleration component and gravity acceleration from the low-frequency acceleration component is estimated to be a translation-contributing low-frequency acceleration component that is a low-frequency component of the translation-contributing acceleration
  • the translation-contributing low-frequency component estimation step can be omitted.
  • the low frequency component and the high frequency component in the time change rate of the posture angle are separated.
  • This frequency separation of the rate of change of posture angle over time for example, derives low-frequency components by using the moving average method for each change amount of posture angle measurement results collected periodically and measures the posture angle. This is done by subtracting the derived low-frequency component from the amount of change in the result to calculate the high-frequency component.
  • the rotation contribution high frequency component estimation step Based on the high frequency component of the time change rate of the posture angle, the rotation that is the high frequency component of the rotation contribution acceleration in the high frequency acceleration component A contributing high frequency acceleration component is estimated. Then, in the translation contribution high-frequency component estimation step, the force obtained by subtracting the rotation contribution high-frequency acceleration component from the high-frequency acceleration component is estimated to be the translation contribution low-frequency high-acceleration component that is the high-frequency component of the translation-contribution acceleration.
  • the rotation contribution low frequency component estimation step Based on the low-frequency component of the time change rate of the posture angle, a rotation-contributing low-frequency acceleration component that is a component of the rotation-contributing acceleration in the low-frequency acceleration component is estimated. Then, in the translation contribution low frequency component estimation process, the rotation contribution low frequency acceleration component and the gravitational acceleration are subtracted from the low frequency acceleration component, which is the low frequency component of the translation contribution acceleration. Presumed to be a component.
  • the translation contribution acceleration and the rotation contribution acceleration are estimated only by the low-frequency acceleration component, the translation contribution contributes to the translation contribution acceleration and rotation.
  • the rotational contribution acceleration can be estimated with improved accuracy.
  • the translational rotational acceleration estimation step includes the low-frequency acceleration estimation step.
  • the absolute value of the time change rate of the posture angle is equal to or greater than the predetermined value, it is estimated that the rotation contribution acceleration component is not included and only the translation contribution acceleration component is included.
  • a high frequency component estimating step of estimating that only the rotation contribution acceleration component is included without including the translation contribution acceleration component.
  • a component that can be evaluated as a direct current component is extracted as the low frequency acceleration component.
  • the low frequency component estimation step It is estimated that the low-frequency acceleration component is derived only from the gravitational acceleration. This is because, in the case of manual movement of a mobile phone device or the like, both rotational acceleration and translational acceleration usually change momentarily, and the acceleration component that can be evaluated as a DC component is constant. This is because it is considered to be derived from gravity acting on the mobile phone device. As a result, the high-frequency acceleration component is estimated to be the result of subtracting the contribution of the acceleration measurement result force gravity acceleration.
  • the posture angle is determined in the high frequency component estimation step. It is determined whether or not the absolute value of the time change rate is less than a predetermined value. If this determination is affirmative, it can be considered that the high-frequency acceleration component contains a lot of translation-contributing acceleration. Therefore, it can be assumed that the high-frequency acceleration component does not include the rotation contribution acceleration component, but includes only the translation contribution acceleration component.
  • the high-frequency acceleration component includes a lot of rotational contribution acceleration. Therefore, it is estimated that the high-frequency acceleration component does not include the translation contribution acceleration component, but only the rotation contribution acceleration component.
  • the high-frequency acceleration component obtained by subtracting the gravitational acceleration contribution from the acceleration measurement result according to the magnitude of the time change rate of the posture angle should be the rotational contribution acceleration or the translation contribution.
  • the choice of acceleration should be made with reasonableness. Therefore, the translation contribution acceleration and the rotation contribution acceleration in the moving object can be estimated with improved accuracy.
  • the translational rotation acceleration estimation step includes the low-frequency acceleration estimation step.
  • a low-frequency component estimation step that estimates that the frequency acceleration component is derived only from gravitational acceleration; and based on the time rate of change of the posture angle, and estimates the rotation-contributing acceleration component that contributes to rotation in the high-frequency acceleration component
  • a translation contribution component estimation step of estimating that the rotation contribution acceleration component obtained by subtracting the rotation contribution acceleration component from the high frequency acceleration component is the translation contribution component estimation step.
  • a component that can be evaluated as a direct current component is extracted as the low frequency acceleration component.
  • the low-frequency acceleration component is derived only from the gravitational acceleration in the low-frequency component estimation step.
  • the high-frequency acceleration component is estimated to be the result of subtracting the contribution of the acceleration measurement result force gravity acceleration.
  • the rotation contribution component estimation step the rotation contribution component estimation step! The rotational contribution acceleration component is estimated based on the time change rate of. Then, in the translation contribution component estimation step, it is estimated that the translation contribution acceleration component is obtained by subtracting the rotation contribution acceleration component from the high frequency acceleration component.
  • the high-frequency acceleration component is estimated to be either the translation contribution acceleration or the rotation contribution acceleration based on the magnitude of the time change rate of the posture angle, it contributes to the translation in the moving object.
  • the translation contribution acceleration and the rotation contribution acceleration contributing to the rotation can be estimated with improved accuracy.
  • the first axis uniquely defined in the moving object, the second axis orthogonal to the first axis, and the first axis and the second axis are orthogonal.
  • An acceleration estimation device comprising: a measurement result collecting means for collecting a measurement result of acceleration along the at least two axes and a measurement result of a posture angle that is a rotation angle of a reference posture force of the moving object; Acceleration frequency separation means for separating low-frequency acceleration components and high-frequency acceleration components in the measurement results of acceleration along at least two axes; attitude angle change rate deriving means for deriving the time change rate of the attitude angle; Said high frequency acceleration component, front Based on the low-frequency acceleration components and the time rate of change of the attitude angle, the translation and rotation acceleration estimation means for estimating and said rotational contribution acceleration and the translational contribution acceleration; a pressurizing speed estimation device comprising a.
  • the measurement result collecting means includes a first axis that is uniquely defined in the moving object, a second axis that is orthogonal to the first axis, a third axis that is orthogonal to the first axis and the second axis. Collect the measurement results of acceleration along at least two of the axes, and the rotation angle from the reference posture, that is, the posture angle of the moving object. Subsequently, the acceleration frequency separation means separates the low frequency acceleration component and the high frequency acceleration component in each acceleration measurement result.
  • the translation / rotation acceleration estimation means performs the time of the high-frequency acceleration component, the low-frequency acceleration component, the posture angle, and the posture angle. Based on the rate of change, the translation contribution acceleration and the rotation contribution acceleration are estimated.
  • low-frequency using the attitude angle Translational acceleration is calculated by subtracting gravitational acceleration from the acceleration component, evaluating the magnitude of rotational contribution acceleration based on the time change rate of the posture angle, and taking into account the general frequency characteristics of translational contribution acceleration and rotational contribution acceleration. And the rotation contribution acceleration are estimated.
  • the acceleration estimation device of the present invention using the acceleration estimation method of the present invention described above, the translation contribution acceleration contributing to the translation of the moving object and the rotation contribution acceleration contributing to the rotation are calculated. Can be estimated. Therefore, according to the acceleration estimation apparatus of the present invention, it is possible to estimate the translation contribution acceleration and the rotation contribution acceleration contributing to the rotation of the moving object with improved accuracy.
  • the translational rotational acceleration estimation means estimates that the high-frequency acceleration component does not include the rotation-contributing acceleration component but includes only the translation-contributing acceleration component. If the absolute value of the time change rate of the posture angle is less than a predetermined value, the component-estimating means does not include the rotation-contributing acceleration component when the gravitational acceleration is subtracted from the low-frequency acceleration component. When it is estimated that only the translation contribution acceleration component is included, and the absolute value of the time change rate of the posture angle is not less than the predetermined value, the gravitational acceleration is subtracted from the low frequency acceleration component, Low-frequency component estimation means for estimating that only the rotation contribution acceleration component is included and not including the translation contribution acceleration component.
  • the high-frequency component estimation means in the translational rotational acceleration estimation means estimates that the high-frequency acceleration component does not include the rotation contribution acceleration component but includes only the translation contribution acceleration component. Further, the low frequency component estimating means in the translational rotation acceleration estimating means determines whether or not the absolute value of the time change rate of the posture angle is less than a predetermined value. If this determination is affirmative, the low-frequency component estimation means subtracts the gravitational acceleration from the low-frequency acceleration component and does not include the rotational contribution acceleration component, but only the translation contribution acceleration component. Estimated to include. On the other hand, if the above judgment is negative, the low-frequency component estimation means does not include the translation contribution acceleration component, but the rotation contribution acceleration component, in which the low-frequency acceleration component force gravity acceleration is subtracted. Is assumed to contain only
  • the low-frequency acceleration component contributes to rotation according to the magnitude of the time change rate of the posture angle.
  • the choice between acceleration and translation-contributing acceleration is done with reasonableness. Therefore, the translation contribution acceleration and the rotation contribution acceleration in the moving object can be estimated with improved accuracy.
  • the posture angle change rate deriving unit derives a low frequency component of the time change rate of the posture angle
  • the translational rotation acceleration estimating unit includes the high frequency acceleration component.
  • High-frequency component estimating means for estimating that it includes only the component of translation-contributed acceleration and does not include the component of rotational contribution acceleration; and the low-frequency component based on the low-frequency component of the time change rate of the posture angle
  • Rotation contribution low frequency component estimating means for estimating the rotation contribution low frequency acceleration component which is a component of acceleration in the frequency acceleration component; and subtracting the rotation contribution low frequency acceleration component and gravitational acceleration from the low frequency acceleration component
  • a translation-contributing low-frequency component estimator that estimates a translation-contributing low-frequency acceleration component that is a low-frequency component of the translation-contributing acceleration.
  • the posture angle change rate deriving means derives a low frequency component of the time change rate of the posture angle.
  • the high frequency component estimation means in the translational rotation acceleration estimation means estimates that the high frequency acceleration component does not include the rotation contribution acceleration component but includes only the translation contribution acceleration component as in the case described above.
  • the rotation contribution low frequency component estimation means in the translational rotation acceleration estimation means is based on the low frequency component of the time change rate of the posture angle, and the rotation contribution low frequency which is a component of the rotation contribution acceleration in the low frequency acceleration component. Estimate the acceleration component.
  • the translation-contributing low-frequency component estimation means is a translation-contributing low-frequency acceleration component that is a low-frequency component of the force-contributing acceleration obtained by subtracting the rotation-contributing low-frequency acceleration component and the heavy acceleration from the low-frequency acceleration component. presume.
  • the accuracy is improved. It is possible to estimate the translation contribution acceleration that contributes to translation in a moving object and the rotation contribution acceleration that contributes to rotation.
  • the posture angle change rate deriving means includes the posture
  • the translational rotation acceleration estimating means separates the low frequency component and the high frequency component in the time change rate of the angle, and the translation rotational acceleration estimation means determines the rotation contribution acceleration of the high frequency acceleration component based on the high frequency component of the time change rate of the posture angle.
  • Rotation contribution high frequency component estimation means for estimating a rotation contribution high frequency acceleration component which is a high frequency component; and a translation contribution high frequency which is a high frequency component of the translation contribution acceleration obtained by subtracting the rotation contribution high frequency acceleration component from the high frequency acceleration component
  • a translation-contributing high-frequency component estimator that estimates an acceleration component; and a rotation-contributing low-frequency acceleration component that is a component of the rotation-contributing acceleration in the low-frequency acceleration component based on the low-frequency component of the time change rate of the posture angle.
  • a rotation-contributing low-frequency component estimator for estimating a component; and the rotation-contributing low-frequency acceleration component from the low-frequency acceleration component.
  • the translational contribution acceleration of translational contribution low frequency component estimation means estimates that the translational contribute low frequency acceleration component which is a low-frequency component; can be configured to include a.
  • the posture angle change rate deriving means separates a low frequency component and a high frequency component in the time change rate of the posture angle.
  • the rotational contribution high-frequency component estimation means in the translational rotational acceleration estimation step estimates a rotational contribution high-frequency acceleration component that is a high-frequency component of the rotational contribution acceleration in the high-frequency acceleration component based on the high-frequency component of the time change rate of the posture angle.
  • the translation-contributing high-frequency component estimating means estimates a value obtained by subtracting the rotation-contributing high-frequency acceleration component from the high-frequency acceleration component as a translation-contributing low-circumference high-acceleration component that is a high-frequency component of the translation-contributing acceleration.
  • the rotation contribution low frequency component estimation means in the translational rotation acceleration estimation means is based on the low frequency component of the time change rate of the posture angle, and the rotation contribution low frequency which is the rotation contribution acceleration component in the low frequency acceleration component. Estimate the acceleration component.
  • the translation-contributing low-frequency component estimation means subtracts the rotation-contributing low-frequency acceleration component and the heavy acceleration from the low-frequency acceleration component to obtain a translation-contributing low-frequency acceleration component that is a low-frequency component of the translation-contributing acceleration. Estimated.
  • the acceleration frequency separation means extracts a component that can be evaluated as a direct current component as the low-frequency acceleration component
  • the translational rotation acceleration estimation means includes the low-frequency acceleration estimation means.
  • Low-frequency component estimation means that estimates that the acceleration component is derived only from gravitational acceleration; and the high-frequency acceleration component contributes to rotation when the absolute value of the time change rate of the posture angle is less than a predetermined value. If the absolute value of the time change rate of the posture angle is greater than or equal to the predetermined value, it is estimated that the translation contribution acceleration component is not included but the rotation contribution acceleration component is included.
  • High-frequency component estimating means for estimating that only the rotational contribution acceleration component is included without including the contribution acceleration component.
  • the acceleration frequency separation means extracts a component that can be evaluated as a direct current component as the low frequency acceleration component.
  • the low frequency component estimating means estimates that the low frequency acceleration component is derived only from the heavy acceleration.
  • the high-frequency acceleration component is estimated to be the result of subtracting the contribution of the acceleration measurement result force gravity acceleration.
  • the high frequency component estimating means determines whether or not the absolute value of the time change rate of the posture angle is less than a predetermined value. If this determination is affirmative, the high-frequency component estimating means estimates that the high-frequency acceleration component does not include the component of rotation contribution acceleration but only the component of translation contribution acceleration. On the other hand, if the above judgment is negative, the high frequency component estimating means estimates that the high frequency acceleration component does not include the translation contribution acceleration component but only the rotation contribution acceleration component.
  • the acceleration frequency separation means extracts a component that can be evaluated as a direct current component as the low-frequency acceleration component, and the translational rotation acceleration estimation means includes the low-frequency acceleration estimation means.
  • Low-frequency component estimating means that estimates that the acceleration component is derived only from gravitational acceleration; and a rotation-contributing component that estimates a rotation-contributing acceleration component that contributes to rotation in the high-frequency acceleration component based on the time change rate of the posture angle An estimation means; a force obtained by subtracting the rotation contribution acceleration component from the high-frequency acceleration component; and a translation contribution component estimation means for estimating the translation contribution acceleration component.
  • the acceleration frequency separation means extracts a component that can be evaluated as a direct current component as the low frequency acceleration component.
  • the low frequency component estimation means estimates that the low frequency acceleration component is derived only from the heavy acceleration.
  • the high-frequency acceleration component is estimated to be the result of subtracting the contribution of the acceleration measurement result force gravity acceleration.
  • the rotational contribution component estimation means estimates the rotational contribution acceleration component based on the time change rate of the posture angle. Then, the translation contribution component estimation means estimates the subtraction of the rotation contribution acceleration component from the high frequency acceleration component as the translation contribution acceleration component.
  • the high-frequency acceleration component is estimated to be either the translation contribution acceleration or the rotation contribution acceleration based on the magnitude of the time change rate of the posture angle, it contributes to the translation in the moving object.
  • the translation contribution acceleration and the rotation contribution acceleration contributing to the rotation can be estimated with improved accuracy.
  • a mobile terminal device is a mobile terminal device that is capable of operating while moving, and includes the acceleration estimation device according to the present invention.
  • the acceleration estimation device of the present invention is used to estimate the translation contribution acceleration contributing to the translation of itself and the rotation contribution acceleration contributing to the rotation. Therefore, according to the mobile terminal apparatus of the present invention, the translation contribution acceleration contributing to its own translation and the rotation contribution acceleration contributing to the rotation can be estimated with improved accuracy. The invention's effect
  • the translation contribution acceleration contributing to translation and the rotation contribution acceleration contributing to rotation in the moving object are accurately determined.
  • the effect that it can improve and can estimate is produced.
  • FIG. 1A is a front view schematically showing an external configuration of a mobile phone according to a first embodiment.
  • FIG. 1B is a rear view schematically showing the external configuration of the mobile phone according to the first embodiment.
  • FIG. 2 is a functional block diagram for explaining the configuration of the mobile phone device of FIGS. 1A and 1B.
  • FIG. 3 is a diagram for explaining the configuration of software executed by the control unit in FIG. 2.
  • FIG. 4 is a functional block diagram for explaining the configuration of the measurement data processing program of FIG. 3.
  • FIG. 5 is a sequence diagram for explaining measurement data processing in the cellular phone device of FIGS. 1A and 1B.
  • FIG. 6 is a flowchart for explaining the processing in the measurement data processing step in FIG.
  • FIG. 7 is a flowchart for explaining processing in an estimation step of translation contribution acceleration and rotation contribution acceleration in FIG.
  • FIG. 8 is a functional block diagram for explaining a configuration of a measurement data processing program in the second embodiment.
  • FIG. 9 is a flowchart for explaining processing in a measurement data processing step in the second embodiment.
  • FIG. 10 is a flowchart for explaining processing in an estimation step of translation contribution acceleration and rotation contribution acceleration in FIG. 9.
  • FIG. 11 is a functional block diagram for explaining a configuration of a measurement data processing program in the third embodiment.
  • FIG. 12 is a flowchart for explaining processing in a measurement data processing step in the third embodiment.
  • FIG. 13 is a flowchart for explaining processing in the translation contribution acceleration and rotation contribution acceleration estimation step in FIG.
  • FIG. 14 is a functional block diagram for explaining a configuration of a measurement data processing program in the fourth embodiment.
  • FIG. 15 is a flowchart for explaining processing in a measurement data processing step in the fourth embodiment.
  • FIG. 16 is a flowchart for explaining the processing in the translation contribution acceleration and rotation contribution acceleration estimation step in FIG.
  • FIG. 17 is a functional block diagram for explaining a configuration of a measurement data processing program in the fifth embodiment.
  • FIG. 18 is a flowchart for explaining processing in a measurement data processing step in the fifth embodiment.
  • FIG. 19 is a flowchart for explaining processing in the translation contribution acceleration and rotation contribution acceleration estimation steps in FIG. 18.
  • FIGS. 1A and 1B schematically show an external configuration of a mobile phone device 10 which is a mobile terminal device.
  • FIG. 1A shows a front view of the appearance of the mobile phone device 10
  • FIG. 1B shows a rear view of the appearance of the mobile phone device 10.
  • the mobile phone device 10 includes (a) a mobile phone body 11, (b) a numeric keypad for inputting a telephone number, and an operation mode.
  • An operation unit 12 having function keys for inputting various instructions such as switching to the control unit 21 (see Fig. 2) to be described later, and (c) a liquid crystal display for displaying operation guidance, operation status, received message, etc.
  • a display unit 13 having a device.
  • the mobile phone device 10 includes (d) a call speaker 14 that reproduces an audio signal that has also been sent by the other party during a call, and (e) a sound input during sound collection.
  • a microphone 15 for inputting voice during a call
  • a guidance speaker 16 for generating ringtones and guidance sounds
  • the antenna 17 is provided.
  • a control unit 21 that performs overall control of the operation of the mobile phone device 10 and a base station via a GO antenna 17 are provided. And (iii) read-only memory (ROM) element and random access memory (RAM) element for storing programs executed by the control unit 21 and various data. And a storage unit 23.
  • the mobile phone body 11 includes a sensor unit 25 for measuring (iv) the attitude angle of the mobile phone device 10 and the acceleration acting on the mobile phone device 10.
  • the storage unit 23 includes therein a temporary storage area 24 for temporarily storing the collected measurement data.
  • the sensor unit 25 detects the attitude angle of the mobile phone device 10 and the acceleration acting on the mobile phone device 10, and outputs the voltage value of the voltage signal from the sensor unit 26 as an analog voltage signal. And an analog-to-digital (A / D) converter 27 for converting values. Then, the digital data output from the AZD conversion 27 is notified to the control unit 21 as measurement data.
  • the sensor unit 25 starts to operate in response to a measurement start command from the control unit 21 and stops operating in response to a measurement stop command from the control unit 21.
  • the row direction in the matrix-like key arrangement of the operation unit 12 is the X axis and the column direction is the Y axis direction.
  • the XY plane is parallel to the horizontal plane.
  • the rotation angle of the reference posture force that is, the pitch angle ⁇ , roll angle ⁇ , and angle of the reference posture relative to the reference posture, with the + Y direction being the true south direction
  • the electromagnetic quantity corresponding to Y is detected.
  • the sensor unit 25 detects the pitch angle ⁇ , roll angle ⁇ , single angle ⁇ , and X direction detected at each measurement point.
  • the detection results corresponding to the velocity ⁇ and the acceleration in the heel direction are combined into a set of raw measurement data.
  • the control unit 21 is notified.
  • the control unit 21 includes a central processing unit (CPU), a digital signal processing unit (DSP), and the like, performs various data processing to realize a general mobile phone function, and performs the above-described processing. It is now possible to control the operation of other components!
  • the configuration of software such as a program executed by the control unit 21 is as shown in FIG.
  • the software in the control unit 21 (i) realizes a call function, a mail function, a character input function, etc., which are basic functions as a mobile phone, and controls the various hardware resources described above.
  • the processing unit 31 and (ii) an application 33 for providing users with various contents such as games are also configured.
  • the application 33 is an application that uses the measurement result of the sensor unit 25.
  • the basic processing unit 31 includes a measurement data processing program 35A.
  • This measurement data processing program 35A collects raw measurement data from the sensor unit 25 in response to a measurement data processing command from the application 33, and converts it into translation-contributing acceleration and rotational motion that contribute to the translational motion of the mobile phone device 10.
  • a process such as estimating the contributing rotation acceleration is performed and stored in the temporary storage area 24 in the storage unit 23. Then, in response to the measurement data request from the application 33, the processed measurement data stored in the temporary storage area 24 in the storage unit 23 is sent to the application 33.
  • the measurement data processing program 35A includes (i) a measurement data processing control unit 49 that performs overall control of the entire measurement data processing program 35A and the operation of the sensor unit 25, and (ii) ) A data collection unit 41 is provided that receives raw measurement data from the sensor unit 25, converts it into physical quantity data such as acceleration and attitude angle, and stores it in the temporary storage area 24 as measurement physical quantity data.
  • the measurement data processing program 35A includes (iii) an acceleration frequency separation unit 42 that separates a low frequency acceleration component and a high frequency acceleration component in each acceleration data in the measurement physical quantity data, and (iv)
  • a posture angle change rate deriving unit 43 for deriving the time change rate of each posture angle data is provided.
  • the acceleration frequency separation unit 42 obtains a low-frequency acceleration component by reading a predetermined number (for example, 10 pieces) of recent data and calculating a moving average value. And The high frequency acceleration component is obtained by subtracting the low frequency acceleration component from the acceleration data at the time point when it is estimated to have the calculated low frequency acceleration component. Further, the attitude angle change rate deriving unit 43 calculates the difference between the posture angle immediately before and the posture angle immediately before it is estimated to have the calculated low-frequency acceleration component as described above. The time change rate of the posture angle at the time when it is estimated to have the low-frequency acceleration component is derived.
  • the measurement data processing program 35A (V) translates the translation contribution acceleration and the rotation contribution acceleration based on the high-frequency acceleration component, the low-frequency acceleration component, the posture angle and the time change rate of the posture angle.
  • Rotation acceleration estimation unit 45A and (vi) a data notification unit 46 that notifies the application 33 of the latest translation contribution acceleration, rotation contribution acceleration, and attitude angle in response to a measurement data request from the application 33.
  • the translational rotation acceleration estimation unit 45A includes a high frequency component estimation unit 51 and a low frequency component estimation unit 52.
  • the high frequency component estimator 51 estimates that the high frequency acceleration component includes only the translation contribution caloric velocity. This is because in the mobile phone device 10, in the case of manual rotational movement, the detection result of the acceleration sensor mounted on the mobile phone device or the like rarely changes abruptly. Moreover, in the case of manual translational movement, the translational contribution acceleration usually has a high frequency component.
  • the low frequency component estimation unit 52 estimates that the low frequency acceleration component force gravity acceleration includes only the translation contribution acceleration. . This is because if the time change rate of the posture angle is less than the predetermined value, the low-frequency acceleration component force minus the gravitational acceleration includes a lot of translation contribution acceleration, and the rotation contribution acceleration includes a slight force. It depends on what is considered.
  • the low frequency component estimation unit 52 estimates that the low frequency acceleration component force gravity acceleration includes only the rotation contribution acceleration. . This is because if the time change rate of the posture angle is greater than or equal to a predetermined value, the low-frequency acceleration component force minus the gravitational acceleration includes a lot of rotational contribution acceleration, and the translation contribution acceleration includes a slight force. It depends on what is considered.
  • the predetermined value is determined by a prior experiment or evaluation.
  • processing of measurement data in the mobile phone device 10 configured as described above will be described mainly with reference to FIGS. 5 to 7 and with reference to other drawings as appropriate.
  • the application 33 sends a measurement data processing start command to the measurement data processing program 35A.
  • the measurement data processing control unit 49 receives this measurement data processing start command.
  • the measurement data processing control unit 49 analyzes the command from the abrasion 33 and recognizes that it is a measurement data processing start command.
  • the measurement data processing control unit 49 determines whether or not the sensor unit 25 is operating. If this determination is negative, the measurement data processing control unit 49 sends a sensor operation start command to the sensor unit 25.
  • the sensor unit 25 starts the sensor operation, and after the detection result by the sensor unit 26 is digitally input by the AZD conversion 27, the sensor unit 25 periodically performs measurement as raw data.
  • the data is output to the control unit 21 (more specifically, to the data collection unit 4 1) (for example, 1 msec cycle).
  • step S11A the measurement data processing program 35A starts the measurement data processing in step S11A.
  • step S11A as shown in FIG. 6, first, in step S21, new measurement raw data from the data collection unit 41 1S sensor unit 25 that has received a measurement data collection start command from the measurement data processing control unit 49 is obtained. Judge whether the force received. If this determination is negative, the process of step S21 is repeated.
  • step S22 the data collection unit 41 calculates physical quantity data that directly represents the magnitude of the acceleration and the posture angle based on the newly received measurement raw data. Then, the data collection unit 41 stores the calculated physical quantity data in the temporary storage area 24 as new measured physical quantity data. When the storage of new measurement physical quantity data in the temporary storage area 24 is thus completed, the data collection unit 41 notifies the measurement data processing control unit 49 to that effect.
  • the storage area for measurement physical quantity data in the temporary storage area 24 has a size capable of storing a set of measurement physical quantity data of a predetermined number or more, which will be described later. It is a ring buffer.
  • step S23 the measurement data processing control unit 49 receives the force that the acceleration frequency separation unit 42 has received a predetermined number or more of measurement raw data for calculating the moving average, ie, the predetermined number or more. It is determined whether the measured physical quantity data is a force stored in the temporary storage area 24 or not. If this determination is negative, the process proceeds to step S21. Then, the processes in steps S21 to S23 are repeated until a positive determination is made in step S23. If a positive determination is made in step S23, the process proceeds to step S24. In step S11A, once a positive determination is made in step S23, a state in which a predetermined number or more of measured physical data is stored in the temporary storage area 24 continues thereafter. Therefore, if step S23 is executed later, a positive determination will be made.
  • step S24 the acceleration frequency separation unit 42 that has received a command from the measurement data processing control unit 49 separates the high-frequency acceleration component and the low-frequency acceleration component based on the acceleration data in the measurement physical quantity data. .
  • the acceleration frequency separation unit 42 first obtains the latest power from the temporary storage area 24 for each of the X-axis direction acceleration data and the Y-axis direction acceleration data. Read a predetermined number in order.
  • the acceleration frequency separation unit 42 obtains a low-frequency acceleration component by calculating a moving average for each of the X-axis direction acceleration data and the Y-axis direction acceleration data.
  • the acceleration frequency separation unit 42 obtains a high-frequency acceleration component at the time point by subtracting the low-frequency acceleration from acceleration data at the time point when the calculated moving average value is regarded as the low-frequency acceleration component.
  • the acceleration frequency separation unit 42 sends the separated high frequency acceleration component and low frequency acceleration component to the translational rotation acceleration estimation unit 45A.
  • step S25A the attitude angle change rate deriving unit 43 that has received a command from the measurement data processing control unit 49 performs the above-described low-frequency acceleration component based on the attitude angle data in the measured physical quantity data.
  • the posture angle change rate deriving unit 43 reads the posture angle data immediately before and after the time point from the temporary storage area 24. Subsequently, the posture angle change rate deriving unit 43 subtracts the corresponding posture angle data value immediately before the current time point from each posture angle data value immediately after the current time point.
  • the attitude angle change rate deriving unit 43 divides the subtraction result by the time difference between the time immediately after the time and the time immediately before.
  • the posture angle change rate deriving unit 43 derives the time change rate of each posture angle at the time, that is, the angular velocity.
  • the posture angle change rate deriving unit 43 also reads out the posture angle data when there is posture angle data at the time point. On the other hand, if there is no posture angle data at the time, the average value of the posture angle data immediately before and after the time is calculated as the posture angle at the time.
  • the posture angle change rate deriving unit 43 sends the posture angle derived as described above and the time change rate of each posture angle to the translational rotation acceleration estimating unit 45A. Note that step S24 and step S25A may be executed in parallel.
  • step S26A the translational rotation acceleration estimation unit 45A performs the high-frequency acceleration component and low-frequency acceleration component from the acceleration frequency separation unit 42, and the posture angle change rate deriving unit 43 from the posture angle.
  • the translation contribution acceleration and the rotation contribution acceleration are estimated based on the time change rate of the posture angle.
  • the high-frequency component estimation unit 51 of the translational rotation acceleration estimation unit 45A estimates that all the high-frequency acceleration components are translation-contributing accelerations, as shown in FIG.
  • step S32 the low-frequency component estimation unit 52 of the translational rotation acceleration estimation unit 45A subtracts the gravitational acceleration component from the low-frequency acceleration component force based on the posture angle received from the posture angle change rate deriving unit 43. After that, it is determined whether at least one of the time change rates of the posture angle received from the posture angle change rate deriving unit 43 is not less than a predetermined value. If this determination result is affirmative, the process proceeds to step S33. In this step S33, the low frequency component estimation unit 52 estimates that all of the subtractions of the low frequency acceleration component force and gravitational acceleration are rotational contribution accelerations. On the other hand, if the determination result in step S32 is negative, the process proceeds to step S34.
  • the low-frequency component estimation unit 52 estimates that all the subtractions of the low-frequency acceleration component force and gravitational acceleration are translation-contributing accelerations. Then, the translational rotation acceleration estimation unit 45A adds the high-frequency acceleration component and the low-frequency acceleration component force minus the gravitational acceleration component, and estimates that the addition result is a translation-contributing acceleration. Is assumed to be 0.
  • the translation rotation acceleration estimation unit 45A sets one set of the estimated translation contribution acceleration, the rotation contribution calorie velocity, and the posture angle in the temporary storage area 24. Stored as notification data.
  • step S26A the process proceeds to step S21 as shown in FIG. Thereafter, the measurement data processing program 35A repeats steps S21 to S26A until a stop request for measurement data processing is received from the application 33.
  • the application 33 sends a data request with the application identifier as a parameter to the measurement data processing program 35A.
  • the measurement data processing control unit 49 receives this data request, analyzes a command from the application 33, and recognizes that the data request is from the application 33. Then, the measurement data processing control unit 49 instructs the data notification unit 46 to notify the application 33 of notification data. Upon receiving this command, the data notification unit 46 receives the data in step S12.
  • step S12 first, the data notification unit 46 sends notification data from the temporary storage area 24.
  • step S12 the data notification unit 46 notifies the application 33 of the read notification data.
  • the application 33 issues the latest translation contribution acceleration and rotation contribution acceleration estimates, and estimates of these by issuing a data request at the desired time. It is possible to obtain the posture angle at the time corresponding to the value.
  • step SI 1 A is the same as the data notification processing in step S 12. It is executed in parallel during the line and continues even after the execution of the data notification process in step S12
  • the application 33 can issue the latest translation-contributing acceleration and rotation by issuing a data request as many times as desired (in Fig. 5, the example of the second case is shown). It is possible to obtain the estimated value of the contribution acceleration and the posture angle at the time point corresponding to these estimated values.
  • the measurement data processing in step S 11 A ends when the application 33 issues a data collection stop request when desired.
  • the measurement data processing program 35A (more specifically, the measurement data processing control unit 49) receiving the request sends a sensor operation stop command to the sensor unit 25. . As a result, the operation of the sensor unit 25 is stopped.
  • the data collection unit 41 has the X axis, the Y axis, and the Z axis that are orthogonally defined and are uniquely defined in the mobile phone device 10. Collect the measurement results of the acceleration along the direction of each axis and the measurement results of the posture angle. Subsequently, the acceleration frequency separation unit 42 separates the low-frequency acceleration component and the high-frequency acceleration component in each acceleration measurement result. The posture angle change rate deriving unit 43 derives the time change rate of the posture angle. Then, the translational rotation acceleration estimation unit 45A estimates the translation contribution acceleration and the rotation contribution acceleration based on the high frequency acceleration component, the low frequency acceleration component, the posture angle and the time change rate of the posture angle.
  • the high frequency component estimation unit 51 in the translational rotation acceleration estimation unit 45A estimates that the high frequency acceleration component does not include the rotation contribution acceleration component but includes only the translation contribution acceleration component.
  • the low frequency component estimation unit 52 in the translational rotation acceleration estimation unit 45A determines whether or not the absolute value of the time change rate of the posture angle is greater than or equal to a predetermined value. If this determination is negative, the low-frequency component estimation unit 52 subtracting the gravitational acceleration from the low-frequency acceleration component does not include the rotational contribution acceleration component, but only the translation contribution acceleration component. It is estimated that it contains. On the other hand, if the determination is affirmative, the low-frequency component estimation unit 52 does not include the translation contribution acceleration component in the subtraction of the low-frequency acceleration component force gravitational acceleration, and the rotation contribution acceleration. It is estimated that only the component of is included.
  • the cellular phone device 10 of the present embodiment is different from the first embodiment in that a measurement data processing program 35B is provided instead of the measurement data processing program 35A.
  • a measurement data processing program 35B is provided instead of the measurement data processing program 35A.
  • the measurement data processing program 35B is replaced with the posture angle change rate deriving unit 43 as compared with the measurement data processing program 35A of the first embodiment.
  • a posture angle change rate deriving unit 43B for calculating a low-frequency component of the rate is provided, and that a translational rotation acceleration estimation unit 45B is provided instead of the translational rotation acceleration estimation unit 45A.
  • this translational rotation acceleration estimation unit 45B estimates the rotation-contributing low frequency acceleration component based on the low-frequency posture angle time change rate component instead of the low-frequency component estimation unit 52.
  • step S11B the measurement data of step S11B (see Fig. 9) is used instead of the measurement data processing of step S11A in the case of the first embodiment. The difference is in the processing.
  • step S11B as shown in FIG. 9, first, in the same manner as in the first embodiment, in steps S21 to S24, the data collection unit 41 collects measurement raw data and calculates measurement physical quantity data.
  • the acceleration frequency separation unit 42 determines whether the high frequency acceleration component and the low frequency acceleration component are based on the acceleration data in the measured physical quantity data. Isolate.
  • step S25B based on the attitude angle data in the attitude angle change rate deriving unit 43B force measurement physical quantity data received the command from the measurement data processing control unit 49,
  • the low frequency component of the time change rate of each posture angle at the time when the above low frequency acceleration component is calculated is derived.
  • the posture angle change rate deriving unit 43B reads the posture angle data from the temporary storage area 24.
  • the posture angle change rate deriving unit 43B obtains the change rate of the posture angle at each time point when the acceleration data used to calculate the low-frequency acceleration component described above is obtained, and the change rate of the posture angle in the first embodiment described above.
  • the calculation is performed in the same manner as the calculation of.
  • the posture angle change rate deriving unit 43B calculates the average of the time change rates of the calculated posture angle. As a result, the low frequency component of the time change rate of each posture angle at the time when the low frequency acceleration component is calculated is derived.
  • the posture angle change rate deriving unit 43B when reading the posture angle data, reads the posture angle data at the time point in the same manner as the posture angle change rate deriving unit 43 in the first embodiment. If it exists, the attitude angle data is also read. On the other hand, when there is no posture angle data at the time, the average value of the posture angle data immediately before and after the time is calculated as the posture angle at the time.
  • the posture angle change rate deriving unit 43B sends the posture angle derived as described above and the low-frequency component of the time change rate of each posture angle to the translational rotation acceleration estimating unit 45B. Note that step S24 and step S25B may be executed in parallel.
  • step S26B the translational rotational acceleration estimation unit 45B performs the high and low frequency acceleration components from the acceleration frequency separation unit 42, and the posture angle and posture angle time from the posture angle change rate deriving unit 43B. Based on the low frequency component of the rate of change, the translation contribution acceleration and the rotation contribution acceleration are estimated. As shown in FIG. 10, in the estimation, the high-frequency component estimation unit 51 of the translational rotation acceleration estimation unit 45B, as shown in FIG. Translational contribution Estimated acceleration.
  • step S42 the rotation contribution low frequency acceleration estimation unit 53 of the translational rotation acceleration estimation unit 45B receives the low frequency component of the posture angle change rate received from the posture angle change rate deriving unit 43B, and Based on the known moment of inertia of the mobile phone device 10 and the position of the center of gravity, the rotation-contributing low-frequency acceleration, which is a low-frequency component of the acceleration that contributes to the rotation. Estimate components. Subsequently, after the translation contribution low-frequency acceleration estimation unit 54 of the translational rotation acceleration estimation unit 45B subtracts the low-frequency acceleration component force gravitational acceleration component based on the posture angle that also received the posture angle change rate deriving unit 43B force in step S43.
  • the translation-contributing low-frequency acceleration estimation unit 54 estimates the thus calculated translation-contributing low-frequency acceleration component, which is a low-frequency component of the acceleration contributing to the translation.
  • the translational rotational acceleration estimation unit 45B adds the high-frequency acceleration component and the translation-contributing low-frequency acceleration component, and estimates the addition result as the translation-contributing acceleration. Further, the translational rotation acceleration estimation unit 45B estimates the rotation-contributing low frequency acceleration component as the rotation-contributing acceleration.
  • the translation rotation acceleration estimation unit 45B stores the estimated translation contribution acceleration and the rotation contribution acceleration, and the posture angle in the temporary storage area 24. Store as notification data.
  • step S26B proceeds to step S21. Thereafter, until the measurement data processing stop request is received from the application 33, the measurement data processing program 35B repeats steps S21 to S26B.
  • the application 33 issues the latest translation contribution acceleration and rotation contribution acceleration estimates, and estimates of these by issuing a data request at the desired time. It is possible to obtain the posture angle at the time corresponding to the value.
  • step S11B The measurement data processing ends when the application 33 issues a data collection stop request when desired.
  • the measurement data processing program 35B (more specifically, the measurement data processing control unit 49) that has received the request issues a sensor operation stop command to the sensor unit 25. Send. As a result, the operation of the sensor unit 25 is stopped.
  • the data collection unit 41 has the X axis, the Y axis, and the Z axis that are orthogonally defined and are defined uniquely in the mobile phone device 10. Collect the measurement results of the acceleration along the direction of each axis and the measurement results of the posture angle. Subsequently, the acceleration frequency separation unit 42 separates the low-frequency acceleration component and the high-frequency acceleration component in each acceleration measurement result. The posture angle change rate deriving unit 43B derives a low frequency component of the time change rate of the posture angle. And translation rotation acceleration estimation part
  • 45B estimates the translation contribution acceleration and the rotation contribution acceleration based on the high frequency acceleration component, the low frequency acceleration component, and the low frequency component of the posture angle and the time change rate of the posture angle.
  • the high frequency component estimation unit 51 in the translational rotation acceleration estimation unit 45B estimates that the high frequency acceleration component does not include the rotation contribution acceleration component but only the translation contribution acceleration component.
  • the rotational contribution low-frequency component estimation unit 53 in the translational rotational acceleration estimation unit 45B estimates the rotational contribution low-frequency component
  • the translational contribution low-frequency component estimation unit 54 in the translational rotational acceleration estimation unit 45B Estimate the low-frequency acceleration component.
  • the translation contribution acceleration contributing to the translation of the mobile phone device 10 and the rotation contribution acceleration contributing to the rotation can be estimated with improved accuracy.
  • the cellular phone device 10 of the present embodiment is different from the second embodiment in that a measurement data processing program 35C is provided instead of the measurement data processing program 35B. Less than
  • the measurement data processing program 35C is a time of posture angle in place of the posture angle change rate deriving unit 43B as compared to the measurement data processing program 35B of the second embodiment.
  • the only difference is that a posture angle change rate deriving unit 43C for deriving a high frequency component and a low frequency component of the rate of change is provided, and that a translational rotational acceleration estimating unit 45C is provided instead of the translational rotational acceleration estimating unit 45B.
  • this translational rotational acceleration estimation unit 45C replaces the high frequency component estimation unit 51 with a rotational contribution that estimates the rotational contribution high frequency acceleration component based on the high frequency attitude angle time change rate component.
  • the only difference is that it comprises a high-frequency component estimating unit 55 and a translation-contributing high-frequency component estimating unit 56 that estimates a translation-contributing high-frequency acceleration component based on the high-frequency acceleration component and the rotation-contributing high-frequency acceleration component.
  • step S 11C the measurement data processing of step S 11C (see FIG. 12) is used instead of the measurement data processing of step S 11B in the case of the second embodiment. Is different.
  • step S11C as shown in FIG. 12, first, in the same manner as in the second embodiment, in steps S21 to S24, the data collection unit 41 collects measurement raw data and calculates measurement physical quantity data. When a predetermined number or more sets of measured physical quantity data are obtained by storing in the temporary storage area 24, the acceleration frequency separation unit 42 performs high frequency acceleration component and low frequency acceleration component based on the acceleration data in the measured physical quantity data. And are separated.
  • step S25C the attitude angle change rate deriving unit 43C that has received a command from the measurement data processing control unit 49 determines that the low-frequency acceleration component described above is based on the attitude angle data in the measured physical quantity data.
  • the high-frequency component and low-frequency component of the time change rate for each posture angle at the calculated time are separated.
  • the posture angle change rate deriving unit 43C reads the posture angle data from the temporary storage area 24.
  • the posture angle change rate deriving unit 43C calculates the change rate of the posture angle at each time point when the acceleration data used for the calculation of the low-frequency acceleration component described above is obtained in the same manner as in the second embodiment. To do.
  • posture angle change rate deriving unit 43C calculates the average of the time change rates of the calculated posture angle. As a result, each posture angle at the time when the low-frequency acceleration component is calculated is calculated. The low frequency component of the time change rate is derived. Then, the posture angle change rate deriving unit 43C subtracts the low frequency component of the time change rate of the posture angle at the time when the low frequency acceleration component is calculated. Thereby, the high frequency component of the time change rate of each posture angle at the time when the low frequency acceleration component is calculated is derived.
  • the posture angle change rate deriving unit 43C like the posture angle change rate deriving unit 43 in the case of the first embodiment, reads the posture angle data when the posture angle data is read out. If it exists, the attitude angle data is also read. On the other hand, when there is no posture angle data at the time, the average value of the posture angle data immediately before and after the time is calculated as the posture angle at the time.
  • the attitude angle change rate deriving unit 43C sends the attitude angle derived as described above and the high frequency component and the low frequency component of the time change rate of each attitude angle to the translational rotation acceleration estimating unit 45C. Note that step S24 and step S25C may be executed in parallel.
  • step S26C the translational rotation acceleration estimation unit 45C force acceleration frequency separation unit 42 high frequency acceleration component and low frequency acceleration component, and posture angle change rate deriving unit 43C posture angle
  • the translation contribution acceleration and the rotation contribution acceleration are estimated based on the high frequency component and low frequency component of the time change rate of the posture angle.
  • the rotation contribution high frequency component estimation unit 55 of the translational rotation acceleration estimation unit 45C changes the posture angle received from the posture angle change rate deriving unit 43C.
  • the rotational contribution high-frequency acceleration component that is the high-frequency component of the acceleration contributing to the rotation is estimated.
  • the translation contribution high frequency acceleration estimation unit 56 of the translation rotation acceleration estimation unit 45C receives from the rotation contribution high frequency acceleration estimation unit 55 from the high frequency acceleration component which also received the posture angle change rate deriving unit 43C force in step S52.
  • the high-frequency acceleration component that contributed to rotation is subtracted.
  • the translation-contributing high-frequency component estimator 56 estimates what is calculated in this way as a translation-contributing high-frequency acceleration component that is a high-frequency component of acceleration that contributes to translation.
  • step S53 the rotation contribution low frequency component estimation unit 53 of the translational rotation acceleration estimation unit 45C reduces the acceleration contributing to the rotation in the same manner as in the second embodiment.
  • the rotation contribution low frequency acceleration component which is a frequency component is estimated.
  • step S54 the translation contribution low frequency component estimation unit 54 of the translational rotation acceleration estimation unit 45C contributes to the translation as a low frequency component of acceleration as in the case of the second embodiment. Translation contribution Estimates the low-frequency acceleration component.
  • the translational rotation acceleration estimation unit 45C adds the rotation-contributing high-frequency acceleration component and the rotation-contributing low-frequency acceleration component, and estimates the addition result as the rotation-contributing acceleration.
  • the translational rotation acceleration estimation unit 45C adds the translation contribution high frequency acceleration component and the translation contribution low frequency acceleration component, and estimates the addition result as the translation contribution acceleration.
  • the translation rotation acceleration estimation unit 45C stores the estimated translation contribution acceleration and rotation contribution acceleration and posture angle in a temporary storage area 24 as a set of notification data. Store as.
  • step S26C ends, the process proceeds to step S21 as shown in FIG. Thereafter, the measurement data processing program 35C repeats steps S21 to S26C until a measurement data processing stop request is received from the application 33.
  • the measurement data processing is performed as in the first and second embodiments.
  • the application 33 issues the latest translation contribution acceleration and rotation contribution acceleration estimates, and estimates of these by issuing a data request at the desired time. It is possible to obtain the posture angle at the time corresponding to the value.
  • the measurement data processing in step S11C requires the application 33 to stop collecting data when desired. It ends by issuing a request.
  • the measurement data processing program 35C (more specifically, the measurement data processing control unit 49) that has received the request sends a sensor operation stop command to the sensor unit 25. As a result, the operation of the sensor unit 25 stops.
  • the data collection unit 41 includes the X axis, the Y axis, and the Z axis that are defined uniquely in the mobile phone device 10, and that are the X axis and the Y axis. Collect the measurement results of the acceleration along the direction of each axis and the measurement results of the posture angle. Subsequently, the acceleration frequency separation unit 42 separates the low-frequency acceleration component and the high-frequency acceleration component in each acceleration measurement result. Also, the posture angle change rate deriving unit 43C separates the low frequency component and the high frequency component of the time change rate of the posture angle. Then, the translational rotation acceleration estimation unit 45C calculates the translation contribution acceleration and the rotation contribution acceleration based on the high frequency acceleration component, the low frequency acceleration component, the posture angle, and the high frequency component and low frequency component of the time change rate of the posture angle. presume.
  • the rotation contribution high frequency component estimation unit 55 in the translational rotation acceleration estimation unit 45C estimates the rotation contribution high frequency component
  • the translation contribution high frequency component estimation unit 56 in the translation rotation acceleration estimation unit 45C Estimate translation contribution high and low frequency acceleration components.
  • the rotation contribution low frequency component estimation unit 53 in the translational rotation acceleration estimation unit 45C estimates the rotation contribution low frequency component
  • the translation contribution low frequency component estimation unit 54 in the translation rotation acceleration estimation unit 45C reduces the translation contribution. Estimate the frequency acceleration component.
  • the translation contribution acceleration contributing to the translation of the mobile phone device 10 and the rotation contribution acceleration contributing to the rotation can be estimated with improved accuracy.
  • the cellular phone device 10 of the present embodiment is different from the first embodiment in that a measurement data processing program 35D is provided instead of the measurement data processing program 35A.
  • the measurement data processing program 35D includes an acceleration frequency separation unit 42D instead of the acceleration frequency separation unit 42, as compared with the measurement data processing program 35A of the first embodiment.
  • a translational rotation acceleration estimation unit 45D is provided instead of the translational rotation acceleration estimation unit 45A.
  • the acceleration frequency separation unit 42D reads a predetermined number (for example, 20 pieces) of recent data, calculates a moving average value, and obtains a low-frequency acceleration component that can be evaluated as a DC component. Then, the acceleration data force at the time when it is estimated to have the calculated low-frequency acceleration component is obtained by subtracting the low-frequency acceleration component.
  • a predetermined number for example, 20 pieces
  • the translational rotation acceleration estimation unit 45D includes a high frequency component estimation unit 51D and a low frequency component estimation unit 52D.
  • the low-frequency component estimation unit 52D estimates that the low-frequency acceleration component includes only the gravitational acceleration component. This is because, in the case of manual operation of the mobile phone device 10, both the rotational acceleration and the translational acceleration usually change little by little, and the acceleration component that can be evaluated as a direct current component is This is due to the fact that it is derived from the gravity that is constantly acting on mobile phone devices.
  • the high frequency component estimation unit 51D estimates that the high frequency acceleration component includes only the translation contribution acceleration. This is because, when the rate of change of the attitude angle with time is less than the predetermined value, the high-frequency acceleration component contains a large amount of translation contribution acceleration, and it is considered that the rotation contribution acceleration does not include a force.
  • the high frequency component estimation unit 51D estimates that the high frequency acceleration component includes only the rotation contribution acceleration. This is because when the time change rate of the posture angle is greater than or equal to a predetermined value, the high-frequency acceleration component includes a lot of rotational contribution acceleration and the translation contribution acceleration does not include a slight force. .
  • step S11D In the mobile phone device 10 of the present embodiment configured as described above, instead of the measurement data processing in step S11A in the case of the first embodiment, the measurement data processing in step S11D (see Fig. 15). Is different.
  • step S11D as shown in FIG. 15, first, in the same manner as in the first embodiment, data collection is performed in steps S21 to S23D.
  • the collecting unit 41 collects the measured raw data, calculates the measured physical quantity data, and stores the measured physical quantity data in the temporary storage area 24, whereby it is determined whether or not the force includes a predetermined number or more sets of measured physical quantity data.
  • the predetermined number is larger than that in the first embodiment, and a value capable of extracting a low frequency component of acceleration change that can be evaluated as a direct current component described later is employed. This predetermined number is determined based on experience, results of experiments and evaluations.
  • step 24D the acceleration frequency separation unit 42 that has received a command from the measurement data processing control unit 49 separates the high-frequency acceleration component and the low-frequency acceleration component based on the acceleration data in the measurement physical quantity data. .
  • the acceleration frequency separation unit 42D calculates a moving average to obtain a low-frequency acceleration component that can be evaluated as a DC component in the same manner as in the first embodiment. Ask.
  • the acceleration frequency separation unit 42D obtains the high frequency acceleration component at the time point by subtracting the low frequency acceleration from the acceleration data at the time point when the calculated moving average value is regarded as the low frequency acceleration component.
  • the acceleration frequency separation unit 42D sends the separated high frequency acceleration component and low frequency acceleration component to the translational rotation acceleration estimation unit 45D.
  • step S25D as in step 25A in the first embodiment, the posture angle change rate deriving unit 43 performs the above-described low-frequency acceleration based on the posture angle data in the measured physical quantity data.
  • the time change rate of each posture angle at the time when the component was calculated is derived.
  • step S26D the translational rotation acceleration estimation unit 45D force acceleration frequency separation unit 42 high frequency acceleration component and low frequency acceleration component, and posture angle change rate deriving unit 43D posture angle and posture angle time change Based on the rate, the translation contribution acceleration and the rotation contribution acceleration are estimated.
  • the low-frequency component estimator 52D of the translational rotation acceleration estimator 45A determines that all of the low-frequency acceleration components are gravitational acceleration components, as shown in FIG. It is estimated that the translation contribution acceleration or the rotation contribution acceleration is not included.
  • step S62 the high-frequency component estimation unit 51D of the translational rotation acceleration estimation unit 45A determines whether at least one of the time change rates of the posture angle received from the posture angle change rate deriving unit 43 is greater than or equal to a predetermined value. judge. If this determination result is affirmative, the process proceeds to step S63. In step S63, the high frequency component estimation unit 51D estimates that all of the high frequency acceleration components are rotational contribution accelerations.
  • step S64 the high frequency component estimation unit 51D estimates that all of the high frequency acceleration components are translational contribution accelerations.
  • the translation rotation acceleration estimation unit 45D sets the estimated translation contribution acceleration, the rotation contribution calorie velocity, and the posture angle as one set in the temporary storage area 24. Stored as notification data.
  • step S26D the process proceeds to step S21 as shown in FIG. Thereafter, the measurement data processing program 35D repeats steps S21 to S26D until a measurement data processing stop request is received from the application 33.
  • the application 33 issues the latest translation contribution acceleration and rotation contribution acceleration estimates, and estimates of these by issuing a data request at the desired time. It is possible to obtain the posture angle at the time corresponding to the value.
  • the measurement data processing in step S11D is performed when the application 33 issues a data collection stop request when desired. End by line.
  • the measurement data processing program 35D (more specifically, the measurement data processing control unit 49) that receives the request issues a sensor operation stop command to the sensor unit 25. Send. As a result, the operation of the sensor unit 25 is stopped.
  • the data collection unit 41 has the X axis, the Y axis, and the Z axis that are orthogonally defined and are uniquely defined in the mobile phone device 10. Collect the measurement results of the acceleration along the direction of each axis and the measurement results of the posture angle. Subsequently, the acceleration frequency separation unit 42D separates the low frequency acceleration component that can be evaluated as the direct current component in each acceleration measurement result and the other high frequency acceleration components. The posture angle change rate deriving unit 43 derives a low frequency component of the time change rate of the posture angle. Then, the translational rotation acceleration estimation unit 45D estimates the translation contribution acceleration and the rotation contribution acceleration based on the high frequency acceleration component, the low frequency acceleration component, the posture angle and the low frequency component of the time change rate of the posture angle.
  • the low-frequency component estimation unit 52D in the translational rotation acceleration estimation unit 45D estimates that the low-frequency acceleration component includes only the gravitational acceleration component. Further, the high frequency component estimator 55D in the translational rotation acceleration estimator 45D determines whether or not the absolute value of the time change rate of the posture angle is greater than or equal to a predetermined value. When this determination is negative, the high frequency component estimation unit 51D estimates that the high frequency acceleration component does not include the rotation contribution acceleration component but includes only the translation contribution acceleration component. On the other hand, when the determination is affirmative, the high frequency component estimation unit 51D estimates that the high frequency acceleration component does not include the translation contribution acceleration component but includes only the rotation contribution acceleration component.
  • the translation contribution acceleration contributing to the translation of the mobile phone device 10 and the rotation contribution acceleration contributing to the rotation can be estimated with improved accuracy.
  • the mobile phone device 10 of the present embodiment has a measurement data processing program.
  • a difference is that a measurement data processing program 35E is provided instead of the ram 35D.
  • this difference will be mainly described.
  • the measurement data processing program 35E is different from the measurement data processing program 35D of the fourth embodiment in place of the translational rotation acceleration estimation unit 45D.
  • This translational rotation acceleration estimation unit 45E is different from the translational rotation acceleration estimation unit 45D in that it replaces the high frequency component estimation unit 51D and estimates the rotation contribution component estimation unit 55 based on the posture angle time change rate.
  • E is provided with a translation contribution component estimation unit 56E that estimates the translation contribution acceleration component based on the high frequency acceleration component and the rotation contribution acceleration component.
  • step S11E the measurement data processing of step SI 1E (see FIG. 18) Is different.
  • the data collection unit 41 collects measurement raw data and calculates measurement physical quantity data.
  • the acceleration frequency separation unit 42D force Based on the acceleration data in the measured physical quantity data, the high-frequency acceleration component and the low-frequency acceleration component ( The DC component and the component that can be evaluated are separated.
  • step 25D the attitude angle change rate deriving unit 43, which has received a command from the measurement data processing control unit 49, determines the attitude in the measured physical quantity data based on the attitude angle data in the measured physical quantity data. Based on the angle data, the time change rate of each posture angle at the time when the low frequency acceleration component is calculated is derived.
  • step S26E the translational rotation acceleration estimation unit 45E force acceleration frequency separation unit 42D high frequency acceleration component and low frequency acceleration component, and posture angle change rate deriving unit 43 posture angle and Based on the low-frequency component of the time change rate of the posture angle, the translation contribution acceleration and the rotation contribution acceleration are estimated.
  • the low frequency component estimation unit 52D of the translational rotation acceleration estimation unit 45E It is estimated that all are gravitational acceleration components and do not include translation-contribution acceleration or rotation-contribution acceleration.
  • step S72 the rotation contribution acceleration estimation unit 55E of the translational rotation acceleration estimation unit 45E receives the change rate of the posture angle received from the posture angle change rate deriving unit 43, and the mobile phone that is known in advance. Based on the moment of inertia of the telephone device 10 and the position of the center of gravity, a rotation contribution acceleration component that is an acceleration component contributing to rotation is estimated. Subsequently, the translation contribution acceleration estimation unit 56E of the translation rotation acceleration estimation unit 45E subtracts the rotation contribution low frequency acceleration component received from the rotation contribution acceleration estimation unit 55E in step S73. The translation contribution acceleration estimation unit 56E thus calculated is estimated to be a translation contribution low frequency acceleration component that is an acceleration component contributing to translation.
  • the translation rotation acceleration estimation unit 45E sets the estimated translation contribution acceleration, the rotation contribution calorie velocity, and the posture angle as one set in the temporary storage area 24. Stored as notification data.
  • step S26E ends, the process proceeds to step S21 as shown in FIG. Thereafter, the measurement data processing program 35E repeats steps S21 to S26E until a measurement data processing stop request is received from the application 33.
  • the application 33 issues the latest translation contribution acceleration and rotation contribution acceleration estimates, and estimates of these by issuing a data request at the desired time. It is possible to obtain the posture angle at the time corresponding to the value.
  • the measurement data processing in step S11E is performed when the application 33 issues a data collection stop request when desired. End by line.
  • the measurement data processing program 35E (more specifically, the measurement data processing control unit 49) that has received the request issues a sensor operation stop command to the sensor unit 25. Send. As a result, the operation of the sensor unit 25 is stopped.
  • the data collection unit 41 has the X axis, the Y axis, and the Z axis that are orthogonally defined and are defined uniquely in the mobile phone device 10. Collect the measurement results of the acceleration along the direction of each axis and the measurement results of the posture angle. Subsequently, the acceleration frequency separation unit 42D separates a low frequency acceleration component that can be evaluated as a direct current component and a high frequency acceleration component in each acceleration measurement result. Also, the posture angle change rate deriving unit 43 derives the time change rate of the posture angle.
  • the translation / rotation acceleration estimation unit 45E calculates the translation-contribution acceleration and the rotation-contribution acceleration. Is estimated.
  • the low-frequency component estimation unit 52D in the translational rotation acceleration estimation unit 45E estimates that the low-frequency acceleration component includes only the gravitational acceleration component. Further, the rotation contribution component estimation unit 55E in the translational rotation acceleration estimation unit 45E estimates the rotation contribution acceleration component, and the translation contribution component estimation unit 56E force in the translation rotation acceleration estimation unit 45E estimates the translation contribution acceleration component.
  • the translation contribution acceleration contributing to the translation of the mobile phone device 10 and the rotation contribution acceleration contributing to the rotation can be estimated with improved accuracy.
  • the sensor unit 25 is mounted inside the mobile phone device 10.
  • the sensor unit 25 is arranged outside the mobile phone device 10, and the sensor unit 25 and the mobile phone device 10 are connected via an interface port (not shown) for connecting the external device of the mobile phone device 10. Even so,
  • the sensor unit 26 of the sensor unit 25 has the pitch angle ⁇ , the roll angle ⁇ , the angle ⁇ , the X-direction acceleration ⁇ , and the heel-direction acceleration ⁇ .
  • the wear Furthermore, it can also be set as the sensor which detects physical quantities other than a posture angle and acceleration.
  • the row direction in the matrix-like key arrangement of the operation unit 12 is the X axis
  • the column direction is the Y axis direction
  • the XY plane is a horizontal plane.
  • the reference posture is the posture where the + Y direction is the true south direction, but other postures can be used as the reference posture. In this case, it is necessary to measure in advance the relationship between the reference posture and the vertical direction, or to obtain the relationship between the posture and the vertical direction at the time when each piece of measurement data is obtained.
  • the sensor unit 25 may be disposed on the operation unit side. It may be arranged on the display unit side.
  • the axial direction serving as a measurement reference in the sensor unit 25 can be determined in accordance with the position where the sensor unit 25 is arranged and the model of the mobile phone.
  • the present invention is applied to a mobile phone device.
  • a portable game machine a car navigation device, a PDA (Personal Digital Assistance), etc.
  • the present invention can be applied to other types of mobile terminal devices.
  • the acceleration estimation method and the acceleration estimation apparatus of the present invention estimate the translation contribution acceleration contributing to the translation of the moving object and the rotation contribution acceleration contributing to the rotation. Can be applied.
  • the mobile terminal device of the present invention is applied to a mobile terminal device that estimates translation-contributing acceleration contributing to its own translation and rotation-contributing acceleration contributing to its own rotation, and uses the estimation result. be able to.

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Human Computer Interaction (AREA)
  • Computer Hardware Design (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Multimedia (AREA)
  • Navigation (AREA)
  • Regulating Braking Force (AREA)

Abstract

 移動物体に搭載されたセンサにより、移動物体の加速度と姿勢角とを計測し、計測結果を収集する(ステップS21,S22)。引き続き、加速度の計測結果それぞれにおける低周波加速度成分と高周波加速度成分とを分離する(ステップS24)とともに、姿勢角の時間変化率を算出導出する(ステップS25A)。こうして、高周波加速度成分及び低周波加速度成分の抽出、並びに姿勢角の時間変化率の導出が終了すると、高周波加速度成分、低周波加速度成分、姿勢角及び姿勢角の時間変化率に基づいて、並進寄与加速度と回転寄与加速度とを推定する(ステップS26A)。この結果、移動物体における並進に寄与している並進寄与加速度と回転に寄与している回転寄与加速度とが、精度良く推定される。

Description

明 細 書
加速度推定方法、加速度推定装置及び移動端末装置
技術分野
[0001] 本発明は、加速度推定方法、加速度推定装置及び移動端末装置に係り、より詳し くは、移動物体における並進に寄与して 、る並進寄与加速度と回転に寄与して 、る 回転寄与加速度とを推定する加速度推定方法と、当該加速度推定方法を使用する 加速度推定装置と、当該加速度推定装置を備える移動端末装置に関するものであ る。
背景技術
[0002] 従来から、携帯電話装置を代表とする移動しつつ動作可能な移動端末装置が、広 く普及している。力かる移動端末装置、特に携帯電話装置に関する技術の進歩は目 覚しいものがあり、必須機能である移動通信網を介した通信機能に加えて、様々な 種類のゲーム等のアプリケーションも実装されるようになって 、る。
[0003] また、移動端末装置においては、自身に作用している加速度を計測するためのセ ンサを実装し、当該センサによる計測データをアプリケーションで利用することが提案 されている。例えば、加速度を計測するセンサを備え、当該センサによって計測され た加速度のデータを用いて、その移動軌跡を求め、その移動軌跡を入力文字として 認識することができる携帯電話装置が提案されて 、る (特許文献 1参照:以下、「従来 例」と呼ぶ)。
特許文献 1 :特開 2002— 169645号公報
発明の開示
発明が解決しょうとする課題
[0004] 上述した従来例の技術は、携帯電話装置等におけるセンサによる計測結果の利用 技術としては優れたものである。しカゝしながら、この従来例では、文字軌跡を描く場合 には、利用者が、携帯電話装置等をほぼ 2次元平面上でほぼ回転を伴わずに並進 移動操作を行うことが必要であり、利用者にとって操作しずらさを伴うものであった。
[0005] また、地表近辺で使用される携帯電話等に装備される加速度センサは、通常、静 止状態において重力加速度を検知する。このため、携帯電話等を静止状態にして、 重力加速度のみが検知された状態を初期状態として、加速度変化を常に監視するこ とにより、重力方向 (鉛直下方)を特定し、重力加速度を検知結果から差し引くことを 可能とすることも考えられるが、初期状態力も長い時間を経過した後には、重力方向 を精度良く特定することは困難であった。このため、姿勢が可変な携帯端末装置等 に装備される加速度センサによる検知結果のみ力 では、検知結果力 重力加速度 を簡易に差し引き、運動に寄与する加速度を精度良く導出することができな力つた。
[0006] そこで、地磁気等を利用して携帯電話装置等の姿勢角を検知するセンサを併設す ることにより、重力方向を精度良く特定することが考えられる。こうして重力方向を特 定することにより、加速度センサによる検知結果力 重力加速度を差し引くことはでき る力 並進運動に寄与する加速度と回転に寄与する加速度とを分離することは相変 わらず困難である。
[0007] こうした並進運動に寄与する加速度と回転に寄与する加速度との分離は、従来例 のような移動軌跡の検出を利用した文字入力において大きな効果を奏するのみなら ず、ゲーム等における利用者の操作環境の充実という観点からも大きな可能を有す るものである。そこで、現在、携帯電話装置等の移動端末装置における並進運動に 寄与する加速度と回転に寄与する加速度とを分離するための技術が強く要請されて いる。
[0008] 本発明は、上記の事情を鑑みてなされたものであり、移動物体における並進に寄 与している並進寄与加速度と回転に寄与している回転寄与加速度とを精度を向上し て推定することができる加速度推定方法及び加速度推定装置を提供することを目的 とする。
[0009] また、本発明は、自身の並進に寄与している並進寄与加速度と回転に寄与してい る回転寄与加速度とを精度を向上して推定することができる移動端末装置を提供す ることを目的とする。
課題を解決するための手段
[0010] 本発明は、第 1の観点力 すると、移動物体において固有に定義される第 1軸、前 記第 1軸と直交する第 2軸、並びに前記第 1軸及び前記第 2軸と直交する第 3軸の内 の少なくとも 2つの軸方向に沿った加速度の計測結果における前記移動物体の並進 に寄与して 、る並進寄与加速度と、回転に寄与して 、る回転寄与加速度とを推定す る加速度推定方法であって、前記少なくとも 2つの軸に沿った加速度の計測結果と、 前記移動物体の基準姿勢力 の回転角である姿勢角の計測結果とを収集する計測 結果収集工程と;前記少なくとも 2つの軸に沿った加速度の計測結果における低周 波加速度成分と高周波加速度成分とを分離する加速度周波数分離工程と;前記姿 勢角の時間変化率を導出する姿勢角変化率導出工程と;前記高周波加速度成分、 前記低周波加速度成分及び前記姿勢角の時間変化率に基づいて、前記並進寄与 加速度と前記回転寄与加速度とを推定する並進回転加速度推定工程と;を備える加 速度推定方法である。
[0011] この加速度推定方法では、計測結果収集工程において、移動物体において固有 に定義される第 1軸、第 1軸と直交する第 2軸、第 1軸及び第 2軸と直交する第 3軸の 内の少なくとも 2つの軸方向に沿った加速度の計測結果、並びに基準姿勢からの回 転角、すなわち移動物体の姿勢角の計測結果が収集される。引き続き、加速度周波 数分離工程にぉ 、て、加速度の計測結果それぞれにおける低周波加速度成分と高 周波加速度成分とが分離される。こうした周波数分離は、例えば、周期的に収集され る加速度の計測結果それぞれについて移動平均法により低周波加速度成分を抽出 するとともに、加速度の計測結果それぞれから抽出された低周波加速度成分を差し 引いて高周波加速度成分を抽出することにより行われる。次いで、姿勢角変化率導 出工程において、姿勢角の時間変化率が導出される。
[0012] こうして、高周波加速度成分及び低周波加速度成分の抽出、並びに姿勢角の時間 変化率の導出が終了すると、並進回転加速度推定工程において、高周波加速度成 分、低周波加速度成分、姿勢角及び姿勢角の時間変化率に基づいて、並進寄与加 速度と回転寄与加速度とが推定される。この結果、姿勢角を利用して低周波加速度 成分から重力加速度を差し引くとともに、姿勢角の時間変化率による回転寄与加速 度の大きさの評価、並びに並進寄与加速度及び回転寄与加速度の一般的な周波数 特性の考慮等がされて、並進寄与加速度と回転寄与加速度とが推定される。
[0013] したがって、本発明の加速度推定方法によれば、移動物体における並進に寄与し ている並進寄与加速度と回転に寄与している回転寄与加速度とを精度を向上して推 定することができる。
[0014] 本発明の加速度推定方法では、前記並進回転加速度推定工程が、前記高周波加 速度成分は、前記回転寄与加速度の成分を含まず、前記並進寄与加速度の成分の みを含むと推定する高周波成分推定工程と;前記姿勢角の時間変化率の絶対値が 所定の値未満の場合には、前記低周波加速度成分から重力加速度を差し引いたも のは、前記回転寄与加速度の成分を含まず、前記並進寄与加速度の成分のみを含 むと推定し、前記姿勢角の時間変化率の絶対値が前記所定の値以上の場合には、 前記低周波加速度成分から重力加速度を差し引 、たものは、前記並進寄与加速度 の成分を含まず、前記回転寄与加速度の成分のみを含むと推定する低周波成分推 定工程と;を備免ることとすることができる。
[0015] この場合には、並進回転加速度推定工程では、まず、高周波成分推定工程にお!、 て、高周波加速度成分は、回転寄与加速度の成分を含まず、並進寄与加速度の成 分のみを含むと推定される。これは、携帯電話装置等の携帯端末装置では、手動に よる回転運動の場合には、携帯電話装置等に搭載された加速度センサの検出結果 が急激に変化することは非常に稀であることによる。また、手動による並進運動の場 合には、並進寄与加速度が高周波成分を有することが通常であることによる。
[0016] 並進回転加速度推定工程では、次に、低周波成分推定工程において、姿勢角の 時間変化率の絶対値が所定の値未満である力否かが判定される。この判定が肯定 的であった場合には、低周波加速度成分力 重力加速度を差し引いたものには、並 進寄与加速度が多く含まれると考えることができる。そこで、低周波加速度成分から 重力加速度を差し引いたものには、回転寄与加速度の成分は含まれず、並進寄与 加速度の成分のみを含むと推定することとして ヽる。
[0017] 一方、上記の判定が否定的であった場合には、低周波加速度成分から重力加速 度を差し引いたものには、回転寄与加速度が多く含まれると考えることができる。そこ で、低周波加速度成分力 重力加速度を差し引いたものには、並進寄与加速度の 成分は含まれず、回転寄与加速度の成分のみを含むと推定される。
[0018] この結果、姿勢角の時間変化率の大きさに応じて、低周波加速度成分を回転寄与 加速度とすべきか、又は、並進寄与加速度とすべきかの二者択一が合理性をもって 行われる。したがって、移動物体における並進寄与加速度と回転寄与加速度とを精 度を向上して推定することができる。
[0019] また、本発明の加速度推定方法では、前記姿勢角変化率導出工程において、前 記姿勢角の時間変化率の低周波成分を導出するとともに、前記並進加速度推定ェ 程が、前記高周波加速度成分は、前記回転寄与加速度の成分を含まず、前記並進 寄与加速度の成分のみを含むと推定する高周波成分推定工程と;前記姿勢角の時 間変化率の低周波成分に基づ 、て、前記低周波加速度成分における前記回転寄 与加速度の成分である回転寄与低周波加速度成分を推定する回転寄与低周波成 分推定工程と;前記低周波加速度成分から前記回転寄与低周波加速度成分及び重 力加速度を差し引いたもの力 前記並進寄与加速度の低周波成分である並進寄与 低周波加速度成分であると推定する並進寄与低周波成分推定工程と;を備えること とすることができる。
[0020] この場合には、姿勢角変化率導出工程において、姿勢角の時間変化率の低周波 成分が導出される。力かる姿勢角の時間変化率における低周波成分の導出は、例え ば、周期的に収集される姿勢角の計測結果の変化量それぞれについて移動平均法 を用いることにより行われる。
[0021] また、並進回転加速度推定工程では、まず、高周波成分推定工程において、上記 の場合と同様に、高周波加速度成分は、回転寄与加速度の成分を含まず、並進寄 与加速度の成分のみを含むと推定される。
[0022] 並進回転加速度推定工程では、次に、回転寄与低周波推定工程において、姿勢 角の時間変化率の低周波成分に基づいて、低周波加速度成分における回転寄与 加速度の成分である回転寄与低周波加速度成分が推定される。そして、並進寄与 低周波成分推定工程にお!ヽて、低周波加速度成分から回転寄与低周波加速度成 分及び重力加速度を差し引いたもの力 並進寄与加速度の低周波成分である並進 寄与低周波加速度成分であると推定される。
[0023] したがって、上記の姿勢角の時間変化率の大きさに基づいて、低周波加速度成分 を並進寄与加速度又は回転寄与加速度のいずれかのみと推定する場合と比べて、 移動物体における並進に寄与している並進寄与加速度と回転に寄与している回転 寄与加速度とを、精度を向上して推定することができる。
[0024] 本発明の加速度推定方法では、前記姿勢角変化率導出工程において、前記姿勢 角の時間変化率における低周波成分と高周波成分を分離するとともに、前記並進カロ 速度推定工程が、前記姿勢角の時間変化率の高周波成分に基づいて、前記高周 波加速度成分における前記回転寄与加速度の高周波成分である回転寄与高周波 加速度成分を推定する回転寄与高周波成分推定工程と;前記高周波加速度成分か ら前記回転寄与高周波加速度成分を差し引いたものが、前記並進寄与加速度の高 周波成分である並進寄与高周波加速度成分であると推定する並進寄与高周波成分 推定工程と;前記姿勢角の時間変化率の低周波成分に基づいて、前記低周波加速 度成分における前記回転寄与加速度の成分である回転寄与低周波加速度成分を 推定する回転寄与低周波成分推定工程と;前記低周波加速度成分から前記回転寄 与低周波加速度成分及び重力加速度を差し引いたものが、前記並進寄与加速度の 低周波成分である並進寄与低周波加速度成分であると推定する並進寄与低周波成 分推定工程と;を備免ることとすることができる。
[0025] この場合には、姿勢角変化率導出工程において、姿勢角の時間変化率における 低周波成分と高周波成分が分離される。かかる姿勢角の時間変化率の周波数分離 は、例えば、周期的に収集される姿勢角の計測結果の変化量それぞれについて移 動平均法を用いることにより低周波成分を導出するともに、姿勢角の計測結果の変 化量から、導出された低周波成分を差し引いて高周波成分を算出することにより行わ れる。
[0026] また、並進回転加速度推定工程では、回転寄与高周波成分推定工程にお!、て、 姿勢角の時間変化率の高周波成分に基づいて、高周波加速度成分における回転 寄与加速度の高周波成分である回転寄与高周波加速度成分が推定される。そして 、並進寄与高周波成分推定工程において、高周波加速度成分から回転寄与高周波 加速度成分を差し引いたもの力 並進寄与加速度の高周波成分である並進寄与低 周高加速度成分であると推定される。
[0027] また、並進回転加速度推定工程では、回転寄与低周波成分推定工程にぉ 、て、 姿勢角の時間変化率の低周波成分に基づいて、低周波加速度成分における回転 寄与加速度の成分である回転寄与低周波加速度成分が推定される。そして、並進 寄与低周波成分推定工程にお!ヽて、低周波加速度成分から回転寄与低周波加速 度成分及び重力加速度を差し引いたものが、並進寄与加速度の低周波成分である 並進寄与低周波加速度成分であると推定される。
[0028] したがって、低周波加速度成分のみにおいて並進寄与加速度及び回転寄与加速 度を推定する上記の場合と比べて、移動物体における並進に寄与して 、る並進寄与 加速度と回転に寄与して 、る回転寄与加速度とを、精度を向上して推定することが できる。
[0029] また、本発明の加速度推定方法では、前記加速度周波数分離工程において、前 記低周波加速度成分として、直流成分と評価できる成分を抽出するとともに、前記並 進回転加速度推定工程が、前記低周波加速度成分が重力加速度にのみ由来すると 推定する低周波成分推定工程と;前記姿勢角の時間変化率の絶対値が所定の値未 満の場合には、前記高周波加速度成分は、回転に寄与する回転寄与加速度の成分 を含まず、前記並進寄与加速度の成分のみを含むと推定し、前記姿勢角の時間変 化率の絶対値が前記所定の値以上の場合には、前記高周波成分は、前記並進寄 与加速度の成分を含まず、前記回転寄与加速度の成分のみを含むと推定する高周 波成分推定工程と;を備えるようにすることができる。
[0030] この場合には、加速度周波数分離工程では、低周波加速度成分として、直流成分 と評価できる成分を抽出する。
[0031] 並進回転加速度推定工程では、低周波成分推定工程にお!ヽて、低周波加速度成 分が重力加速度にのみ由来すると推定する。これは、携帯電話装置等の手動による 運動の場合には、回転加速度及び並進加速度の双方ともに、僅かではあっても刻々 変化することが通常であり、直流成分と評価できる加速度成分は、定常的に携帯電 話装置等に作用している重力に由来すると考えられるからである。この結果、高周波 加速度成分は、加速度の計測結果力 重力加速度の寄与分が差し引かれたものと なっていると推定されることになる。
[0032] また、並進回転加速度推定工程では、高周波成分推定工程において、姿勢角の 時間変化率の絶対値が所定の値未満である力否かが判定される。この判定が肯定 的であった場合には、高周波加速度成分には並進寄与加速度が多く含まれると考え ることができる。そこで、高周波加速度成分には、回転寄与加速度の成分は含まれ ず、並進寄与加速度の成分のみが含まれると推定することとして ヽる。
[0033] 一方、上記の判定が否定的であった場合には、高周波加速度成分には、回転寄 与加速度が多く含まれると考えることができる。そこで、高周波加速度成分には、並 進寄与加速度の成分は含まれず、回転寄与加速度の成分のみが含まれると推定さ れる。
[0034] この結果、姿勢角の時間変化率の大きさに応じて、加速度の計測結果から重力加 速度の寄与分が差し引かれた高周波加速度成分を回転寄与加速度とすべきか、又 は、並進寄与加速度とすべきかの二者択一が合理性をもって行われる。したがって、 移動物体における並進寄与加速度と回転寄与加速度とを精度を向上して推定する ことができる。
[0035] また、本発明の加速度推定方法では、前記加速度周波数分離工程において、前 記低周波加速度成分として、直流成分と評価できる成分を抽出するとともに、前記並 進回転加速度推定工程が、前記低周波加速度成分が重力加速度にのみ由来すると 推定する低周波成分推定工程と;前記姿勢角の時間変化率に基づ!、て、前記高周 波加速度成分における回転に寄与する回転寄与加速度成分を推定する回転寄与 成分推定工程と;前記高周波加速度成分から前記回転寄与加速度成分を差し引い たものが、前記並進寄与加速度成分であると推定する並進寄与成分推定工程と;を 備えるようにすることができる。
[0036] この場合には、加速度周波数分離工程にぉ 、て、低周波加速度成分として、直流 成分と評価できる成分を抽出する。
[0037] 並進回転加速度推定工程では、低周波成分推定工程にお!ヽて、低周波加速度成 分が重力加速度にのみ由来すると推定する。この結果、高周波加速度成分は、加速 度の計測結果力 重力加速度の寄与分が差し引かれたものとなっていると推定され ることになる。
[0038] また、並進回転加速度推定工程では、回転寄与成分推定工程にお!、て、姿勢角 の時間変化率に基づいて、回転寄与加速度成分が推定される。そして、並進寄与成 分推定工程において、高周波加速度成分から回転寄与加速度成分を差し引いたも のが、並進寄与加速度成分であると推定される。
[0039] したがって、上記の姿勢角の時間変化率の大きさに基づいて、高周波加速度成分 を並進寄与加速度又は回転寄与加速度のいずれかのみと推定する場合と比べて、 移動物体における並進に寄与している並進寄与加速度と回転に寄与している回転 寄与加速度とを、精度を向上して推定することができる。
[0040] 本発明は、第 2の観点力 すると、移動物体において固有に定義される第 1軸、前 記第 1軸と直交する第 2軸、並びに前記第 1軸及び前記第 2軸と直交する第 3軸の内 の少なくとも 2つの軸方向に沿った加速度の計測結果における前記移動物体の並進 に寄与して 、る並進寄与加速度と、回転に寄与して 、る回転寄与加速度とを推定す る加速度推定装置であって、前記少なくとも 2つの軸に沿った加速度の計測結果と、 前記移動物体の基準姿勢力 の回転角である姿勢角の計測結果とを収集する計測 結果収集手段と;前記少なくとも 2つの軸に沿った加速度の計測結果における低周 波加速度成分と高周波加速度成分とを分離する加速度周波数分離手段と;前記姿 勢角の時間変化率を導出する姿勢角変化率導出手段と;前記高周波加速度成分、 前記低周波加速度成分及び前記姿勢角の時間変化率に基づいて、前記並進寄与 加速度と前記回転寄与加速度とを推定する並進回転加速度推定手段と;を備える加 速度推定装置である。
[0041] この加速度推定装置では、計測結果収集手段が、移動物体において固有に定義 される第 1軸、第 1軸と直交する第 2軸、第 1軸及び第 2軸と直交する第 3軸の内の少 なくとも 2つの軸方向に沿った加速度の計測結果、並びに基準姿勢からの回転角、 すなわち移動物体の姿勢角の計測結果を収集する。引き続き、加速度周波数分離 手段が、加速度の計測結果それぞれにおける低周波加速度成分と高周波加速度成 分とを分離する。こうして、高周波加速度成分及び低周波加速度成分の抽出、並び に姿勢角の時間変化率の導出が終了すると、並進回転加速度推定手段が、高周波 加速度成分、低周波加速度成分、姿勢角及び姿勢角の時間変化率に基づいて、並 進寄与加速度と回転寄与加速度とを推定する。この結果、姿勢角を利用して低周波 加速度成分から重力加速度を差し引くとともに、姿勢角の時間変化率による回転寄 与加速度の大きさの評価、並びに並進寄与加速度及び回転寄与加速度の一般的な 周波数特性の考慮等がされて、並進寄与加速度と回転寄与加速度とが推定される。
[0042] すなわち、本発明の加速度推定装置では、上述した本発明の加速度推定方法を 使用して、移動物体の並進に寄与している並進寄与加速度と回転に寄与している回 転寄与加速度とを推定することができる。したがって、本発明の加速度推定装置によ れば、移動物体における並進に寄与して 、る並進寄与加速度と回転に寄与して 、る 回転寄与加速度とを精度を向上して推定することができる。
[0043] 本発明の加速度推定装置では、前記並進回転加速度推定手段は、前記高周波加 速度成分が、前記回転寄与加速度の成分を含まず、前記並進寄与加速度の成分の みを含むと推定する高周波成分推定手段と;前記姿勢角の時間変化率の絶対値が 所定の値未満の場合には、前記低周波加速度成分から重力加速度を差し引いたも のは、前記回転寄与加速度の成分を含まず、前記並進寄与加速度の成分のみを含 むと推定し、前記姿勢角の時間変化率の絶対値が前記所定の値以上の場合には、 前記低周波加速度成分から重力加速度を差し引 、たものは、前記並進寄与加速度 の成分を含まず、前記回転寄与加速度の成分のみを含むと推定する低周波成分推 定手段と;を備える構成とすることができる。
[0044] この場合には、並進回転加速度推定手段における高周波成分推定手段は、高周 波加速度成分が、回転寄与加速度の成分を含まず、並進寄与加速度の成分のみを 含むと推定する。また、並進回転加速度推定手段における低周波成分推定手段が、 姿勢角の時間変化率の絶対値が所定の値未満である力否かを判定する。この判定 が肯定的であった場合には、低周波成分推定手段が、低周波加速度成分から重力 加速度を差し引いたものには、回転寄与加速度の成分は含まれず、並進寄与加速 度の成分のみを含むと推定する。一方、上記の判定が否定的であった場合には、低 周波成分推定手段が、低周波加速度成分力 重力加速度を差し引いたものには、 並進寄与加速度の成分は含まれず、回転寄与加速度の成分のみを含むと推定する
[0045] この結果、姿勢角の時間変化率の大きさに応じて、低周波加速度成分を回転寄与 加速度とすべきか、又は、並進寄与加速度とすべきかの二者択一が合理性をもって 行われる。したがって、移動物体における並進寄与加速度と回転寄与加速度とを精 度を向上して推定することができる。
[0046] また、本発明の加速度推定装置では、前記姿勢角変化率導出手段が、前記姿勢 角の時間変化率の低周波成分を導出するとともに、前記並進回転加速度推定手段 は、前記高周波加速度成分が、前記回転寄与加速度の成分を含まず、前記並進寄 与加速度の成分のみを含むと推定する高周波成分推定手段と;前記姿勢角の時間 変化率の低周波成分に基づ 、て、前記低周波加速度成分における前記回転寄与 加速度の成分である回転寄与低周波加速度成分を推定する回転寄与低周波成分 推定手段と;前記低周波加速度成分から前記回転寄与低周波加速度成分及び重力 加速度を差し引いたもの力 前記並進寄与加速度の低周波成分である並進寄与低 周波加速度成分であると推定する並進寄与低周波成分推定手段と;を備える構成と することができる。
[0047] この場合には、姿勢角変化率導出手段が、姿勢角の時間変化率の低周波成分を 導出する。また、並進回転加速度推定手段における高周波成分推定手段が、上記 の場合と同様に、高周波加速度成分は、回転寄与加速度の成分を含まず、並進寄 与加速度の成分のみを含むと推定する。
[0048] また、並進回転加速度推定手段における回転寄与低周波成分推定手段が、姿勢 角の時間変化率の低周波成分に基づいて、低周波加速度成分における回転寄与 加速度の成分である回転寄与低周波加速度成分を推定する。そして、並進寄与低 周波成分推定手段が、低周波加速度成分から回転寄与低周波加速度成分及び重 力加速度を差し引いたもの力 並進寄与加速度の低周波成分である並進寄与低周 波加速度成分であると推定する。
[0049] したがって、上記の姿勢角の時間変化率の大きさに基づいて、低周波加速度成分 を並進寄与加速度又は回転寄与加速度のいずれかのみと推定する場合と比べて、 精度を向上して、移動物体における並進に寄与して 、る並進寄与加速度と回転に寄 与して 、る回転寄与加速度とを推定することができる。
[0050] また、本発明の加速度推定装置では、前記姿勢角変化率導出手段が、前記姿勢 角の時間変化率における低周波成分と高周波成分を分離するとともに、前記並進回 転加速度推定手段が、前記姿勢角の時間変化率の高周波成分に基づいて、前記 高周波加速度成分における前記回転寄与加速度の高周波成分である回転寄与高 周波加速度成分を推定する回転寄与高周波成分推定手段と;前記高周波加速度成 分から前記回転寄与高周波加速度成分を差し引いたものが、前記並進寄与加速度 の高周波成分である並進寄与高周波加速度成分であると推定する並進寄与高周波 成分推定手段と;前記姿勢角の時間変化率の低周波成分に基づいて、前記低周波 加速度成分における前記回転寄与加速度の成分である回転寄与低周波加速度成 分を推定する回転寄与低周波成分推定手段と;前記低周波加速度成分から前記回 転寄与低周波加速度成分及び重力加速度を差し引いたものが、前記並進寄与加速 度の低周波成分である並進寄与低周波加速度成分であると推定する並進寄与低周 波成分推定手段と;を備える構成とすることができる。
[0051] この場合には、姿勢角変化率導出手段が、姿勢角の時間変化率における低周波 成分と高周波成分とを分離する。また、並進回転加速度推定工程における回転寄与 高周波成分推定手段が、姿勢角の時間変化率の高周波成分に基づいて、高周波 加速度成分における回転寄与加速度の高周波成分である回転寄与高周波加速度 成分を推定する。そして、並進寄与高周波成分推定手段が、高周波加速度成分から 回転寄与高周波加速度成分を差し引いたものを、並進寄与加速度の高周波成分で ある並進寄与低周高加速度成分であると推定する。
[0052] また、並進回転加速度推定手段における回転寄与低周波成分推定手段が、姿勢 角の時間変化率の低周波成分に基づいて、低周波加速度成分における回転寄与 加速度の成分である回転寄与低周波加速度成分を推定する。そして、並進寄与低 周波成分推定手段が、低周波加速度成分から回転寄与低周波加速度成分及び重 力加速度を差し引いたものを、並進寄与加速度の低周波成分である並進寄与低周 波加速度成分であると推定する。
[0053] したがって、低周波加速度成分のみにおいて並進寄与加速度及び回転寄与加速 度を推定する上記の場合と比べて、更に精度を向上して、移動物体における並進に 寄与している並進寄与加速度と回転に寄与している回転寄与加速度とを推定するこ とがでさる。
[0054] また、本発明の加速度推定装置では、前記加速度周波数分離手段が、前記低周 波加速度成分として、直流成分と評価できる成分を抽出するとともに、前記並進回転 加速度推定手段が、前記低周波加速度成分が重力加速度にのみ由来すると推定す る低周波成分推定手段と;前記姿勢角の時間変化率の絶対値が所定の値未満の場 合には、前記高周波加速度成分は、回転に寄与する回転寄与加速度の成分を含ま ず、前記並進寄与加速度の成分のみを含むと推定し、前記姿勢角の時間変化率の 絶対値が前記所定の値以上の場合には、前記高周波成分は、前記並進寄与加速 度の成分を含まず、前記回転寄与加速度の成分のみを含むと推定する高周波成分 推定手段と;を備えるようにすることができる。
[0055] この場合には、加速度周波数分離手段が、低周波加速度成分として、直流成分と 評価できる成分を抽出する。
[0056] 並進回転加速度推定手段では、低周波成分推定手段が、低周波加速度成分が重 力加速度にのみ由来すると推定する。この結果、高周波加速度成分は、加速度の計 測結果力 重力加速度の寄与分が差し引かれたものとなっていると推定されることに なる。
[0057] また、並進回転加速度推定手段では、高周波成分推定手段が、姿勢角の時間変 化率の絶対値が所定の値未満であるか否かを判定する。この判定が肯定的であった 場合には、高周波成分推定手段は、高周波加速度成分には、回転寄与加速度の成 分は含まれず、並進寄与加速度の成分のみが含まれると推定する。一方、上記の判 定が否定的であった場合には、高周波成分推定手段は、高周波加速度成分には、 並進寄与加速度の成分は含まれず、回転寄与加速度の成分のみを含まれると推定 する。
[0058] この結果、姿勢角の時間変化率の大きさに応じて、加速度の計測結果から重力加 速度の寄与分が差し引かれた高周波加速度成分を回転寄与加速度とすべきか、又 は、並進寄与加速度とすべきかの二者択一が合理性をもって行われる。したがって、 移動物体における並進寄与加速度と回転寄与加速度とを精度を向上して推定する ことができる。 [0059] また、本発明の加速度推定装置では、前記加速度周波数分離手段が、前記低周 波加速度成分として、直流成分と評価できる成分を抽出するとともに、前記並進回転 加速度推定手段が、前記低周波加速度成分が重力加速度にのみ由来すると推定す る低周波成分推定手段と;前記姿勢角の時間変化率に基づいて、前記高周波加速 度成分における回転に寄与する回転寄与加速度成分を推定する回転寄与成分推 定手段と;前記高周波加速度成分から前記回転寄与加速度成分を差し引いたもの 力 前記並進寄与加速度成分であると推定する並進寄与成分推定手段と;を備える ようにすることができる。
[0060] この場合には、加速度周波数分離手段が、低周波加速度成分として、直流成分と 評価できる成分を抽出する。
[0061] 並進回転加速度推定手段では、低周波成分推定手段が、低周波加速度成分が重 力加速度にのみ由来すると推定する。この結果、高周波加速度成分は、加速度の計 測結果力 重力加速度の寄与分が差し引かれたものとなっていると推定されることに なる。
[0062] また、並進回転加速度推定手段では、回転寄与成分推定手段が、姿勢角の時間 変化率に基づいて、回転寄与加速度成分を推定する。そして、並進寄与成分推定 手段が、高周波加速度成分から回転寄与加速度成分を差し引いたものを、並進寄 与加速度成分であると推定する。
[0063] したがって、上記の姿勢角の時間変化率の大きさに基づいて、高周波加速度成分 を並進寄与加速度又は回転寄与加速度のいずれかのみと推定する場合と比べて、 移動物体における並進に寄与している並進寄与加速度と回転に寄与している回転 寄与加速度とを、精度を向上して推定することができる。
[0064] 本発明の移動端末装置は、移動しつつ動作が可能な移動端末装置において、本 発明の加速度推定装置を備えることを特徴とする移動端末装置である。この移動端 末装置では、本発明の加速度推定装置を使用して、自身の並進に寄与している並 進寄与加速度と回転に寄与している回転寄与加速度とを推定する。したがって、本 発明の移動端末装置によれば、自身の並進に寄与している並進寄与加速度と回転 に寄与している回転寄与加速度とを精度を向上して推定することができる。 発明の効果
[0065] 以上説明したように、本発明の加速度推定方法及び加速度推定装置によれば、移 動物体における並進に寄与している並進寄与加速度と回転に寄与している回転寄 与加速度とを精度を向上して推定することができるという効果を奏する。
[0066] また、本発明の移動端末装置によれば、自身の並進に寄与している並進寄与加速 度と回転に寄与している回転寄与加速度とを精度を向上して推定することができる。 図面の簡単な説明
[0067] [図 1A]第 1実施形態に係る携帯電話の外観構成を概略的に示す正面図である。
[図 1B]第 1実施形態に係る携帯電話の外観構成を概略的に示す背面図である。
[図 2]図 1A及び図 1Bの携帯電話装置の構成を説明するための機能ブロック図であ る。
[図 3]図 2の制御部で実行されるソフトウェアの構成を説明するための図である。
[図 4]図 3の計測データ処理プログラムの構成を説明するための機能ブロック図であ る。
[図 5]図 1A及び図 1Bの携帯電話装置における計測データ処理を説明するためのシ 一ケンス図である。
[図 6]図 5における計測データ処理ステップにおける処理を説明するためのフローチ ヤートである。
[図 7]図 6における並進寄与加速度及び回転寄与加速度の推定ステップにおける処 理を説明するためのフローチャートである。
[図 8]第 2実施形態における計測データ処理プログラムの構成を説明するための機能 ブロック図である。
[図 9]第 2実施形態における計測データ処理ステップにおける処理を説明するための フローチャートである。
[図 10]図 9における並進寄与加速度及び回転寄与加速度の推定ステップにおける処 理を説明するためのフローチャートである。
[図 11]第 3実施形態における計測データ処理プログラムの構成を説明するための機 能ブロック図である。 [図 12]第 3実施形態における計測データ処理ステップにおける処理を説明するため のフローチャートである。
[図 13]図 12における並進寄与加速度及び回転寄与加速度の推定ステップにおける 処理を説明するためのフローチャートである。
[図 14]第 4実施形態における計測データ処理プログラムの構成を説明するための機 能ブロック図である。
[図 15]第 4実施形態における計測データ処理ステップにおける処理を説明するため のフローチャートである。
[図 16]図 15における並進寄与加速度及び回転寄与加速度の推定ステップにおける 処理を説明するためのフローチャートである。
[図 17]第 5実施形態における計測データ処理プログラムの構成を説明するための機 能ブロック図である。
[図 18]第 5実施形態における計測データ処理ステップにおける処理を説明するため のフローチャートである。
[図 19]図 18における並進寄与加速度及び回転寄与加速度の推定ステップにおける 処理を説明するためのフローチャートである。
発明を実施するための最良の形態
[0068] 《第 1実施形態》
以下、本発明の第 1実施形態を、図 1A〜図 7を参照して説明する。
[0069] 図 1A及び図 1Bには、移動端末装置である携帯電話装置 10の外観構成が概略的 に示されている。ここで、図 1Aには、携帯電話装置 10の外観の正面図が示され、図 1Bには、携帯電話装置 10の外観の背面図が示されている。
[0070] 図 1A、図 IBに総合的に示されるように、携帯電話装置 10は、(a)携帯電話本体 1 1と、(b)電話番号を入力するためのテンキー、及び、動作モードの切替等の各種指 令を、後述する制御部 21 (図 2参照)に入力するためのファンクションキーを有する操 作部 12と、(c)操作案内、動作状況、受信メッセージ等を表示する液晶表示装置を 有する表示部 13とを備えている。また、携帯電話装置 10は、(d)通話時に通信相手 力も送られてきた音声信号を再生する通話用スピーカ 14と、(e)集音時に音を入力 したり、通話時に音声を入力したりするためのマイクロフォン 15と、(f)着信音や案内 音を発生するための案内用スピーカ 16と、(g)基地局との間で無線信号を授受する ためのアンテナ 17とを備えている。
[0071] 携帯電話本体 11の内部には、図 2に示されるように、(i)携帯電話装置 10全体の 動作を統括制御する制御部 21と、 GOアンテナ 17を介して、基地局との間で無線信 号の送受信を行う送受信部 22と、 (iii)制御部 21で実行されるプログラムや各種デー タを格納する、読出専用メモリ (ROM)素子やランダムアクセスメモリ (RAM)素子を 有する記憶部 23とを備えている。また、携帯電話本体 11の内部には、(iv)携帯電話 装置 10の姿勢角及び携帯電話装置 10に作用する加速度を計測するためのセンサ ユニット 25を備えている。
[0072] 記憶部 23は、その内部に、収集した計測データを一時的に記憶する一時記憶領 域 24を含んでいる。
[0073] センサユニット 25は、携帯電話装置 10の姿勢角及び携帯電話装置 10に作用する 加速度を検知し、アナログ電圧信号として出力するセンサ部 26と、センサ部 26から 電圧信号の電圧値をデジタル値に変換するアナログデジタル (A/D)変換器 27とを 備えている。そして、 AZD変翻27から出力されるデジタルデータが、計測データ として制御部 21へ通知される。また、センサユニット 25は、制御部 21からの計測開始 指令により動作を開始し、制御部 21からの計測停止指令により動作を停止するように なっている。
[0074] センサ部 26は、操作部 12のマトリクス状のキー配列における行方向を X軸とし、列 方向を Y軸方向とするとともに、本実施形態においては、 XY面が水平面と平行であ り、 +Y方向が真南方向である姿勢を基準姿勢として、当該基準姿勢力 の回転角、 すなわち当該基準姿勢に対する姿勢角であるピッチ角 Θ 、ロール角 Θ 及びョ一角
X Y
Θ
Zに対応する電磁気量を検知する。また、センサ部 26は、 X方向の加速度 )
X及 ひ Ύ方向の加速度 )
Y に対応する電磁気量を検知する。そして、センサユニット 25 は、各計測時点において検知されたピッチ角 Θ 、ロール角 Θ 、ョ一角 Θ 、 X方向加
X Y Z
速度 α 及び Υ方向加速度 α に対応する検知結果を 1組にし、計測生データとして
X Υ
制御部 21へ通知するようになって 、る。 [0075] 制御部 21は、中央処理装置 (CPU)、デジタル信号処理装置 (DSP)等を備えて おり、一般的な携帯電話機能を実現するために、様々なデータ処理を行うとともに、 上述した他の構成要素の動作制御を行うようになって!/、る。この制御部 21にお ヽて 実行されるプログラム等のソフトウェアの構成は、図 3に示されるようになつている。
[0076] すなわち、制御部 21におけるソフトウェアは、(i)携帯電話としての基本機能である 通話機能、メール機能、文字入力機能等を実現するとともに、上述した各種のハード ウェア資源の制御を行う基本処理部 31と、 (ii)ゲーム等の様々なコンテンツを利用者 に提供するためのアプリケーション 33と力も構成されている。ここで、アプリケーション 33は、センサユニット 25による計測結果を利用するアプリケーションであるものとする
[0077] 基本処理部 31は、計測データ処理プログラム 35Aを備えている。この計測データ 処理プログラム 35Aは、アプリケーション 33からの計測データ処理指令に応じて、セ ンサユニット 25からの計測生データを収集し、携帯電話装置 10の並進運動に寄与 する並進寄与加速度と回転運動に寄与する回転寄与加速度を推定する等の処理を して、記憶部 23内の一時記憶領域 24に格納する。そして、アプリケーション 33から の計測データ要求に応じて、記憶部 23内の一時記憶領域 24に記憶された処理済 計測データをアプリケーション 33へ送る。
[0078] 計測データ処理プログラム 35Aは、図 4に示されるように、(i)計測データ処理プロ グラム 35Aの全体及びセンサユニット 25の動作を統括制御する計測データ処理制 御部 49と、(ii)センサユニット 25からの計測生データを受け、加速度や姿勢角という 物理量データに変換し、計測物理量データとして、一時記憶領域 24に格納するデ ータ収集部 41とを備えている。また、計測データ処理プログラム 35Aは、(iii)計測物 理量データ内の加速度データそれぞれにおける低周波加速度成分と高周波加速度 成分とを分離させる加速度周波数分離部 42と、 (iv)計測物理量データ内の姿勢角 データそれぞれにおける時間変化率を導出する姿勢角変化率導出部 43とを備えて いる。
[0079] ここで、加速度周波数分離部 42は、所定個(例えば、 10個)の最近のデータを読 み出して、移動平均値を算出することにより、低周波加速度成分を求める。そして、 算出された低周波加速度成分を有すると推定される時点における加速度データから 当該低周波加速度成分を差し引くことにより、高周波加速度成分を求める。また、姿 勢角変化率導出部 43は、上記の算出された低周波加速度成分を有すると推定され る時点の直前の姿勢角と直後の姿勢角との差を算出することにより、上記の算出され た低周波加速度成分を有すると推定される時点における姿勢角の時間変化率を導 出する。
[0080] また、計測データ処理プログラム 35Aは、(V)高周波加速度成分、低周波加速度 成分、姿勢角及び姿勢角の時間変化率に基づいて、並進寄与加速度と回転寄与加 速度とを推定する並進回転加速度推定部 45Aと、 (vi)アプリケーション 33からの計 測データ要求に応じて、最新の並進寄与加速度、回転寄与加速度及び姿勢角をァ プリケーシヨン 33へ通知するデータ通知部 46とを備える。
[0081] 並進回転加速度推定部 45Aは、高周波成分推定部 51と、低周波成分推定部 52と を備えている。ここで、高周波成分推定部 51は、高周波加速度成分が並進寄与カロ 速度のみを含むと推定する。これは、携帯電話装置 10では、手動による回転運動の 場合には、携帯電話装置等に搭載された加速度センサの検出結果が急激に変化す ることは非常に稀であることによる。また、手動による並進運動の場合には、並進寄与 加速度が高周波成分を有することが通常であることによる。
[0082] また、低周波成分推定部 52は、姿勢角の時間変化率が所定値未満である場合に は、低周波加速度成分力 重力加速度を差し引いたものは並進寄与加速度のみを 含むと推定する。これは、姿勢角の時間変化率が所定値未満である場合には、低周 波加速度成分力 重力加速度を差し引いたものには、並進寄与加速度が多く含まれ 、回転寄与加速度は少しし力含まれな 、と考えられることによる。
[0083] 一方、低周波成分推定部 52は、姿勢角の時間変化率が所定値以上である場合に は、低周波加速度成分力 重力加速度を差し引いたものは回転寄与加速度のみを 含むと推定する。これは、姿勢角の時間変化率が所定値以上である場合には、低周 波加速度成分力 重力加速度を差し引いたものには、回転寄与加速度が多く含まれ 、並進寄与加速度は少しし力含まれな 、と考えられることによる。
[0084] なお、当該所定値は、事前の実験や評価により定められる。 [0085] 次に、上記のように構成された携帯電話装置 10における計測データの処理につい て、主に図 5〜図 7を参照し、適宜他の図面を参照して説明する。
[0086] 携帯電話装置 10における計測データの処理では、図 5に示されるように、まず、ァ プリケーシヨン 33が、計測データ処理開始指令を、計測データ処理プログラム 35A へ向けて送る。計測データ処理プログラム 35Aでは、計測データ処理制御部 49が、 この計測データ処理開始指令を受ける。そして、計測データ処理制御部 49は、アブ リケーシヨン 33からの指令を解析し、計測データ処理開始指令であることを認識する 。引き続き、計測データ処理制御部 49は、センサユニット 25が動作中である力否か を判定する。この判定が否定的であった場合には、計測データ処理制御部 49は、セ ンサユニット 25へ向けて、センサ動作開始指令を送る。
[0087] なお、センサ動作開始指令を受けると、センサユニット 25は、センサ動作を開始し、 センサ部 26による検知結果を AZD変翻27によりデジタルィ匕した後に、計測生デ ータとして定期的(例えば、 1msec周期)に制御部 21 (より詳しくは、データ収集部 4 1)へ向けて出力する。
[0088] この後、計測データ処理プログラム 35Aは、ステップ S 11Aの計測データ処理を開 始する。このステップ S11Aでは、図 6に示されるように、まず、ステップ S21において 、計測データ処理制御部 49から計測データ収集開始指令を受けたデータ収集部 41 1S センサユニット 25からの新たな計測生データを受けた力否かを判定する。この判 定が否定的であつた場合には、ステップ S 21の処理が繰り返される。
[0089] そして、データ収集部 41がセンサユニット 25からの新たな計測生データを受け、ス テツプ S21において肯定的な判定がなされると、処理はステップ S 22へ移行する。こ のステップ S22では、データ収集部 41が、新たに受けた計測生データに基づいて、 加速度や姿勢角の大きさを直接的に表す物理量データを算出する。そして、データ 収集部 41が、算出された物理量データを新たな計測物理量データとして一時記憶 領域 24に格納する。こうして新たな計測物理量データの一時記憶領域 24への格納 が完了すると、データ収集部 41は、その旨を計測データ処理制御部 49へ通知する 。なお、一時記憶領域 24における計測物理量データ用記憶領域は、後述する所定 数以上の数の組の計測物理量データを記憶できる大きさを有するとともに、いわゆる リングバッファとされている。
[0090] 次に、ステップ S23において、計測データ処理制御部 49が、加速度周波数分離部 42が移動平均を算出するための所定数以上の計測生データを受けた力、すなわち 、当該所定数以上の計測物理量データが一時記憶領域 24に格納された力否かを判 定する。この判定が否定的であった場合には、処理はステップ S21へ移行する。そし て、ステップ S23において肯定な判定がなされるまで、ステップ S21〜S23の処理を 繰り返す。そして、ステップ S23において肯定的な判定がなされると、処理はステップ S24へ移行する。なお、ステップ S 11 Aの処理においては、ステップ S23において一 度肯定的な判定がなされると、その後は一時記憶領域 24には所定数以上の計測物 理データが格納されている状態が継続するので、後にステップ S23が実行された場 合には、肯定的な判定がなされることになる。
[0091] ステップ S24では、計測データ処理制御部 49からの指令を受けた加速度周波数分 離部 42が、計測物理量データにおける加速度データに基づいて、高周波加速度成 分と低周波加速度成分とを分離する。かかる高周波加速度成分と低周波加速度成 分との分離に際して、加速度周波数分離部 42は、まず、 X軸方向加速度データ及び Y軸方向加速度データそれぞれにつ 、て、一時記憶領域 24から最新のもの力も順 に所定数を読み出す。引き続き、加速度周波数分離部 42は、 X軸方向加速度デー タ及び Y軸方向加速度データそれぞれについて、移動平均を算出することにより、低 周波加速度成分を求める。次に、加速度周波数分離部 42は、算出された移動平均 値が低周波加速度成分とされる時点の加速度データから当該低周波加速度を差し 引くことにより、当該時点における高周波加速度成分を求める。
[0092] こうして、当該時点の X軸方向加速度データ及び Y軸方向加速度データそれぞれ における高周波加速度成分と低周波加速度成分とが分離される。加速度周波数分 離部 42は、分離された高周波加速度成分と低周波加速度成分とを並進回転加速度 推定部 45Aへ送る。
[0093] 次に、ステップ S25Aにおいて、計測データ処理制御部 49からの指令を受けた姿 勢角変化率導出部 43が、計測物理量データにおける姿勢角データに基づいて、上 記の低周波加速度成分が算出された時点における姿勢角それぞれの時間変化率を 導出する。力かる姿勢角それぞれの時間変化率の導出に際しては、まず、姿勢角変 化率導出部 43が、当該時点の直前及び直後の姿勢角データを一時記憶領域 24か ら読み出す。引き続き、姿勢角変化率導出部 43は、当該時点の直後における姿勢 角データ値それぞれから、当該時点の直前における対応する姿勢角データ値を減 算する。そして、姿勢角変化率導出部 43は、減算結果を、当該時点の直後時点と直 前時点との時間差で除算する。こうして、姿勢角変化率導出部 43は、当該時点にお ける姿勢角それぞれの時間変化率、すなわち角速度を導出する。
[0094] なお、姿勢角変化率導出部 43は、上記の姿勢角データ読み出しの際に、当該時 点における姿勢角データが存在する場合には、その姿勢角データも合わせて読み 出す。一方、当該時点における姿勢角データが存在しない場合には、上記の当該時 点の直前と直後の姿勢角データの平均値を算出し、当該時点における姿勢角とする
[0095] 姿勢角変化率導出部 43は、以上のようにして導出された姿勢角、並びに姿勢角そ れぞれの時間変化率を、並進回転加速度推定部 45Aへ送る。なお、ステップ S24と ステップ S25Aとは、並行して実行してもよい。
[0096] 次 、で、ステップ S26Aにお 、て、並進回転加速度推定部 45Aが、加速度周波数 分離部 42からの高周波加速度成分及び低周波加速度成分、並びに姿勢角変化率 導出部 43からの姿勢角及び姿勢角の時間変化率に基づいて、並進寄与加速度と 回転寄与加速度とを推定する。力かる推定に際しては、図 7に示されるように、ステツ プ S31において、並進回転加速度推定部 45Aの高周波成分推定部 51が、高周波 加速度成分は、全てが並進寄与加速度であると推定する。
[0097] 引き続き、ステップ S32において、並進回転加速度推定部 45Aの低周波成分推定 部 52が、姿勢角変化率導出部 43から受けた姿勢角に基づいて低周波加速度成分 力も重力加速度分を差し引 、た後に、姿勢角変化率導出部 43から受けた姿勢角の 時間変化率の少なくとも 1つが所定値以上力否かを判定する。この判定結果が肯定 的であった場合には、処理はステップ S33へ移行する。このステップ S33では、低周 波成分推定部 52が、低周波加速度成分力 重力加速度分を差し引いたものは、全 てが回転寄与加速度であると推定する。 [0098] 一方、ステップ S32における判定結果が否定的であった場合には、処理はステップ S34へ移行する。このステップ S34では、低周波成分推定部 52が、低周波加速度成 分力 重力加速度分を差し引いたものは、全てが並進寄与加速度であると推定する 。そして、並進回転加速度推定部 45Aは、高周波加速度成分と、低周波加速度成 分力 重力加速度分を差し引いたものとを加算し、加算結果が並進寄与加速度であ ると推定するとともに、回転寄与加速度は 0であると推定する。
[0099] こうして、並進寄与加速度及び回転寄与加速度が推定されると、並進回転加速度 推定部 45Aは、一時記憶領域 24に、推定された並進寄与加速度及び回転寄与カロ 速度、並びに姿勢角を 1組の通知用データとして格納する。こうして、ステップ S26A が終了すると、図 6に示されるように、処理はステップ S21へ移行する。以後、計測デ ータ処理の停止要求をアプリケーション 33から受けるまで、計測データ処理プロダラ ム 35Aは、ステップ S21〜S26Aを繰り返す。
[0100] 図 5に戻り、上記のステップ S11Aの計測データ処理における所望の時点で、ァプ リケーシヨン 33が、アプリケーション識別子をパラメータとするデータ要求を、計測デ ータ処理プログラム 35Aへ向けて送る。計測データ処理プログラム 35Aでは、計測 データ処理制御部 49が、このデータ要求を受け、アプリケーション 33からの指令を 解析し、アプリケーション 33からのデータ要求であることを認識する。そして、計測デ ータ処理制御部 49は、データ通知部 46に対して、アプリケーション 33への通知用デ ータの通知指令を行う。この指令を受けたデータ通知部 46は、ステップ S12のデー
1 タ通知処理を行う。
[0101] ステップ S12では、まず、データ通知部 46が、一時記憶領域 24から通知用データ
1
を読み出す。引き続き、データ通知部 46は、読み出された通知用データをアプリケ ーシヨン 33へ通知する。こうして、ステップ S12の処理が終了する。
1
[0102] 以上のようにして、計測データ処理が行われることにより、アプリケーション 33は、所 望の時点におけるデータ要求の発行により、最新の並進寄与加速度及び回転寄与 加速度の推定値、並びにこれらの推定値に対応する時点における姿勢角を取得す ることがでさる。
[0103] なお、ステップ SI 1 Aの計測データ処理は、ステップ S 12のデータ通知処理の実 行中も並行して実行されるし、ステップ S12のデータ通知処理の実行終了後も継続
1
して実行される。このため、アプリケーション 33は、所望の時点で何度でも(図 5にお いては、 2度の場合の例が示されている)データ要求を発行することにより、最新の並 進寄与加速度及び回転寄与加速度の推定値、並びにこれらの推定値に対応する時 点における姿勢角を取得することができる。
[0104] また、ステップ S 11 Aの計測データ処理は、アプリケーション 33が、所望の時に、デ ータ収集停止要求を発行することにより終了する。また、アプリケーション 33がセンサ 動作停止要求を発行すると、これを受けた計測データ処理プログラム 35A (より詳しく は、計測データ処理制御部 49)は、センサユニット 25へ向けて、センサ動作停止指 令を送る。この結果、センサユニット 25の動作が停止する。
[0105] 以上説明したように、本第 1実施形態では、データ収集部 41が、携帯電話装置 10 において固有に定義される互いに直交する X軸、 Y軸及び Z軸の内の X軸及び Y軸 それぞれに沿った方向の加速度の計測結果、並びに姿勢角の計測結果を収集する 。引き続き、加速度周波数分離部 42が、加速度の計測結果それぞれにおける低周 波加速度成分と高周波加速度成分とを分離する。また、姿勢角変化率導出部 43が 姿勢角の時間変化率を導出する。そして、並進回転加速度推定部 45Aが、高周波 加速度成分、低周波加速度成分、姿勢角及び姿勢角の時間変化率に基づいて、並 進寄与加速度と回転寄与加速度とを推定する。
[0106] ここで、並進回転加速度推定部 45Aにおける高周波成分推定部 51が、高周波加 速度成分は、回転寄与加速度の成分を含まず、並進寄与加速度の成分のみを含む と推定する。また、並進回転加速度推定部 45Aにおける低周波成分推定部 52が、 姿勢角の時間変化率の絶対値が所定の値以上である力否かを判定する。この判定 が否定的であった場合には、低周波成分推定部 52が、低周波加速度成分から重力 加速度を差し引いたものには、回転寄与加速度の成分は含まれず、並進寄与加速 度の成分のみを含むと推定する。一方、当該判定が肯定的であった場合には、低周 波成分推定部 52が、低周波加速度成分力 重力加速度を差し引いたものには、並 進寄与加速度の成分は含まれず、回転寄与加速度の成分のみを含むと推定する。
[0107] したがって、移動物体である携帯電話装置 10の並進に寄与している並進寄与加速 度と回転に寄与している回転寄与加速度とを精度を向上して推定することができる。
[0108] 《第 2実施形態》
以下、本発明の第 2実施形態を説明する。なお、本実施形態の説明においては、 上述した第 1実施形態と同一又は同等の要素については同一の符号を付し、重複す る説明を省略する。
[0109] 本実施形態の携帯電話装置 10は、第 1実施形態と比べて、計測データ処理プログ ラム 35Aに代えて、計測データ処理プログラム 35Bを備える点が相違している。以下 、この相違点に主に着目して説明する。
[0110] 計測データ処理プログラム 35Bは、図 8に示されるように、第 1実施形態の計測デ ータ処理プログラム 35Aと比べて、姿勢角変化率導出部 43に代えて、姿勢角の時間 変化率の低周波成分を算出する姿勢角変化率導出部 43Bを備えるとともに、並進回 転加速度推定部 45Aに代えて、並進回転加速度推定部 45Bを備える点のみが相違 する。この並進回転加速度推定部 45Bは、並進回転加速度推定部 45Aと比べて、 低周波成分推定部 52に代えて、低周波姿勢角時間変化率成分に基づいて、回転 寄与低周波加速度成分を推定する回転寄与低周波成分推定部 53と、低周波加速 度成分、姿勢角及び回転寄与低周波加速度成分に基づいて、並進寄与低周波加 速度成分を推定する並進寄与低周波成分推定部 54とを備える点のみが相違してい る。
[0111] 上記のように構成された本実施形態の携帯電話装置 10では、第 1実施形態の場合 におけるステップ S 11 Aの計測データ処理に代えて、ステップ S11B (図 9参照)の計 測データ処理を行う点が相違している。このステップ S11Bでは、図 9に示されるよう に、まず、第 1実施形態と同様にして、ステップ S21〜S24において、データ収集部 4 1が、計測生データを収集し、計測物理量データを算出して一時記憶領域 24に格納 することにより、所定数以上の組の計測物理量データが揃うと、加速度周波数分離部 42が、計測物理量データにおける加速度データに基づいて、高周波加速度成分と 低周波加速度成分とを分離する。
[0112] 次に、ステップ S25Bにおいて、計測データ処理制御部 49からの指令を受けた姿 勢角変化率導出部 43B力 計測物理量データにおける姿勢角データに基づいて、 上記の低周波加速度成分が算出された時点における姿勢角それぞれの時間変化 率の低周波成分を導出する。力かる姿勢角それぞれの時間変化率の導出に際して は、まず、姿勢角変化率導出部 43Bが、姿勢角データを一時記憶領域 24から読み 出す。引き続き、姿勢角変化率導出部 43Bは、上述した低周波加速度成分の算出 に利用した加速度データが得られた各時点における姿勢角の変化率を、上述した第 1実施形態における姿勢角の変化率の算出の場合と同様にして、算出する。
[0113] そして、姿勢角変化率導出部 43Bは、算出された姿勢角の時間変化率の平均を算 出する。これにより、低周波加速度成分が算出された時点における姿勢角それぞれ の時間変化率の低周波成分が導出される。
[0114] なお、姿勢角変化率導出部 43Bは、第 1実施形態の場合の姿勢角変化率導出部 4 3と同様に、上記の姿勢角データ読み出しの際に、当該時点における姿勢角データ が存在する場合には、その姿勢角データも合わせて読み出す。一方、当該時点にお ける姿勢角データが存在しない場合には、上記の当該時点の直前と直後の姿勢角 データの平均値を算出し、当該時点における姿勢角とする。
[0115] 姿勢角変化率導出部 43Bは、以上のようにして導出された姿勢角、並びに姿勢角 それぞれの時間変化率の低周波成分を、並進回転加速度推定部 45Bへ送る。なお 、ステップ S24とステップ S25Bとは、並行して実行してもよい。
[0116] 次いで、ステップ S26Bにおいて、並進回転加速度推定部 45Bが、加速度周波数 分離部 42からの高周波加速度成分及び低周波加速度成分、並びに姿勢角変化率 導出部 43Bからの姿勢角及び姿勢角の時間変化率の低周波成分に基づいて、並 進寄与加速度と回転寄与加速度とを推定する。力かる推定に際しては、図 10に示さ れるように、ステップ S41において、並進回転加速度推定部 45Bの高周波成分推定 部 51が、第 1実施形態の場合と同様にして、高周波加速度成分は、全てが並進寄与 加速度であると推定する。
[0117] 次に、ステップ S42において、並進回転加速度推定部 45Bの回転寄与低周波加 速度推定部 53が、姿勢角変化率導出部 43Bから受けた姿勢角の変化率の低周波 成分、並びに、予め知られている携帯電話装置 10の慣性モーメント及び重心位置に 基づ 、て、回転に寄与して 、る加速度の低周波成分である回転寄与低周波加速度 成分を推定する。引き続き、並進回転加速度推定部 45Bの並進寄与低周波加速度 推定部 54が、ステップ S43において、姿勢角変化率導出部 43B力も受けた姿勢角 に基づいて低周波加速度成分力 重力加速度分を差し引いた後に、更に回転寄与 低周波加速度推定部 53から受けた回転寄与低周波加速度成分を差し引く。こうして 算出されたものを、並進寄与低周波加速度推定部 54は、並進に寄与している加速 度の低周波成分である並進寄与低周波加速度成分であると推定する。
[0118] 次 ヽで、並進回転加速度推定部 45Bは、高周波加速度成分と並進寄与低周波加 速度成分とを加算し、加算結果を並進寄与加速度と推定する。また、並進回転加速 度推定部 45Bは、回転寄与低周波加速度成分を回転寄与加速度と推定する。こうし て、並進寄与加速度及び回転寄与加速度が推定されると、並進回転加速度推定部 45Bは、一時記憶領域 24に、推定された並進寄与加速度及び回転寄与加速度、並 びに姿勢角を 1組の通知用データとして格納する。こうして、ステップ S26Bが終了す ると、図 9に示されるように、処理はステップ S21へ移行する。以後、計測データ処理 の停止要求をアプリケーション 33から受けるまで、計測データ処理プログラム 35Bは 、ステップ S21〜S26Bを繰り返す。
[0119] 上記のステップ S11Bの計測データ処理における所望の時点で、アプリケーション 3 3が、アプリケーション識別子をパラメータとするデータ要求を発行すると、第 1実施形 態の場合と同様に、計測データ処理プログラム 35Bでは、計測データ処理制御部 49 力 このデータ要求を受け、アプリケーション 33からの指令を解析し、アプリケーショ ン 33からのデータ要求であることを認識する。そして、計測データ処理制御部 49は、 データ通知部 46に対して、アプリケーション 33への通知用データの通知指令を行う 。この指令を受けたデータ通知部 46は、ステップ S12 (j = l, 2, · · ·)のデータ通知処 理を行う(図 5参照)。
[0120] 以上のようにして、計測データ処理が行われることにより、アプリケーション 33は、所 望の時点におけるデータ要求の発行により、最新の並進寄与加速度及び回転寄与 加速度の推定値、並びにこれらの推定値に対応する時点における姿勢角を取得す ることがでさる。
[0121] また、本実施形態においても、第 1実施形態の場合と同様にして、ステップ S11Bの 計測データ処理は、アプリケーション 33が、所望の時に、データ収集停止要求を発 行することにより終了する。また、アプリケーション 33がセンサ動作停止要求を発行す ると、これを受けた計測データ処理プログラム 35B (より詳しくは、計測データ処理制 御部 49)は、センサユニット 25へ向けて、センサ動作停止指令を送る。この結果、セ ンサユニット 25の動作が停止する。
[0122] 以上説明したように、本第 2実施形態では、データ収集部 41が、携帯電話装置 10 において固有に定義される互いに直交する X軸、 Y軸及び Z軸の内の X軸及び Y軸 それぞれに沿った方向の加速度の計測結果、並びに姿勢角の計測結果を収集する 。引き続き、加速度周波数分離部 42が、加速度の計測結果それぞれにおける低周 波加速度成分と高周波加速度成分とを分離する。また、姿勢角変化率導出部 43B が姿勢角の時間変化率の低周波成分を導出する。そして、並進回転加速度推定部
45Bが、高周波加速度成分、低周波加速度成分、姿勢角及び姿勢角の時間変化率 の低周波成分に基づ!ヽて、並進寄与加速度と回転寄与加速度とを推定する。
[0123] ここで、並進回転加速度推定部 45Bにおける高周波成分推定部 51が、高周波加 速度成分は、回転寄与加速度の成分を含まず、並進寄与加速度の成分のみを含む と推定する。また、並進回転加速度推定部 45Bにおける回転寄与低周波成分推定 部 53が、回転寄与低周波加速度成分を推定するとともに、並進回転加速度推定部 4 5Bにおける並進寄与低周波成分推定部 54が、並進寄与低周波加速度成分を推定 する。
[0124] したがって、携帯電話装置 10の並進に寄与している並進寄与加速度と回転に寄与 している回転寄与加速度とを精度を向上して推定することができる。
[0125] 《第 3実施形態》
以下、本発明の第 3実施形態を説明する。なお、本実施形態の説明においては、 上述した第 1又は第 2実施形態と同一又は同等の要素については同一の符号を付し
、重複する説明を省略する。
[0126] 本実施形態の携帯電話装置 10は、第 2実施形態と比べて、計測データ処理プログ ラム 35Bに代えて、計測データ処理プログラム 35Cを備える点が相違している。以下
、この相違点に主に着目して説明する。 [0127] 計測データ処理プログラム 35Cは、図 11に示されるように、第 2実施形態の計測デ ータ処理プログラム 35Bと比べて、姿勢角変化率導出部 43Bに代えて、姿勢角の時 間変化率の高周波成分と低周波成分と導出する姿勢角変化率導出部 43Cを備える とともに、並進回転加速度推定部 45Bに代えて、並進回転加速度推定部 45Cを備え る点のみが相違する。この並進回転加速度推定部 45Cは、並進回転加速度推定部 45Bと比べて、高周波成分推定部 51に代えて、高周波姿勢角時間変化率成分に基 づいて、回転寄与高周波加速度成分を推定する回転寄与高周波成分推定部 55と、 高周波加速度成分及び回転寄与高周波加速度成分に基づ!、て、並進寄与高周波 加速度成分を推定する並進寄与高周波成分推定部 56とを備える点のみが相違して いる。
[0128] 上記のように構成された本実施形態の携帯電話装置 10では、第 2実施形態の場合 におけるステップ S 11Bの計測データ処理に代えて、ステップ S 11C (図 12参照)の 計測データ処理を行う点が相違している。このステップ S11Cでは、図 12に示される ように、まず、第 2実施形態と同様にして、ステップ S21〜S24において、データ収集 部 41が、計測生データを収集し、計測物理量データを算出して一時記憶領域 24に 格納することにより、所定数以上の組の計測物理量データが揃うと、加速度周波数分 離部 42が、計測物理量データにおける加速度データに基づいて、高周波加速度成 分と低周波加速度成分とを分離する。
[0129] 次に、ステップ S25Cにおいて、計測データ処理制御部 49からの指令を受けた姿勢 角変化率導出部 43Cが、計測物理量データにおける姿勢角データに基づいて、上 記の低周波加速度成分が算出された時点における姿勢角それぞれの時間変化率の 高周波成分と低周波成分とを分離する。かかる姿勢角それぞれの時間変化率の高 周波成分と低周波成分との分離に際しては、まず、姿勢角変化率導出部 43Cが、姿 勢角データを一時記憶領域 24から読み出す。引き続き、姿勢角変化率導出部 43C は、第 2実施形態の場合と同様にして、上述した低周波加速度成分の算出に利用し た加速度データが得られた各時点における姿勢角の変化率を算出する。
[0130] 次いで、姿勢角変化率導出部 43Cは、算出された姿勢角の時間変化率の平均を 算出する。これにより、低周波加速度成分が算出された時点における姿勢角それぞ れの時間変化率の低周波成分が導出される。そして、姿勢角変化率導出部 43Cは、 低周波加速度成分が算出された時点における姿勢角の変化率力 姿勢角の時間変 化率の低周波成分を差し引く。これにより、低周波加速度成分が算出された時点に おける姿勢角それぞれの時間変化率の高周波成分が導出される。
[0131] なお、姿勢角変化率導出部 43Cは、第 1実施形態の場合の姿勢角変化率導出部 43と同様に、上記の姿勢角データ読み出しの際に、当該時点における姿勢角デー タが存在する場合には、その姿勢角データも合わせて読み出す。一方、当該時点に おける姿勢角データが存在しない場合には、上記の当該時点の直前と直後の姿勢 角データの平均値を算出し、当該時点における姿勢角とする。
[0132] 姿勢角変化率導出部 43Cは、以上のようにして導出された姿勢角、並びに姿勢角 それぞれの時間変化率の高周波成分及び低周波成分を、並進回転加速度推定部 4 5Cへ送る。なお、ステップ S24とステップ S25Cとは、並行して実行してもよい。
[0133] 次 、で、ステップ S26Cにお 、て、並進回転加速度推定部 45C力 加速度周波数 分離部 42からの高周波加速度成分及び低周波加速度成分、並びに姿勢角変化率 導出部 43Cからの姿勢角、姿勢角の時間変化率の高周波成分及び低周波成分に 基づいて、並進寄与加速度と回転寄与加速度とを推定する。力かる推定に際しては 、図 13に示されるように、ステップ S51において、並進回転加速度推定部 45Cの回 転寄与高周波成分推定部 55が、姿勢角変化率導出部 43Cから受けた姿勢角の変 化率の高周波成分、並びに、予め知られている携帯電話装置 10の慣性モーメント及 び重心位置に基づいて、回転に寄与している加速度の高周波成分である回転寄与 高周波加速度成分を推定する。
[0134] 引き続き、並進回転加速度推定部 45Cの並進寄与高周波加速度推定部 56が、ス テツプ S52において、姿勢角変化率導出部 43C力も受けた高周波加速度成分から、 回転寄与高周波加速度推定部 55から受けた回転寄与高周波加速度成分を差し引 く。こうして算出されたものを、並進寄与高周波成分推定部 56は、並進に寄与してい る加速度の高周波成分である並進寄与高周波加速度成分であると推定する。
[0135] 次に、ステップ S53において、並進回転加速度推定部 45Cの回転寄与低周波成 分推定部 53が、第 2実施形態の場合と同様にして、回転に寄与している加速度の低 周波成分である回転寄与低周波加速度成分を推定する。引き続き、ステップ S54に おいて、並進回転加速度推定部 45Cの並進寄与低周波成分推定部 54が、第 2実施 形態の場合と同様にして、並進に寄与して 、る加速度の低周波成分である並進寄与 低周波加速度成分を推定する。
[0136] 次いで、並進回転加速度推定部 45Cは、回転寄与高周波加速度成分と回転寄与 低周波加速度成分とを加算し、この加算結果を回転寄与加速度と推定する。また、 並進回転加速度推定部 45Cは、並進寄与高周波加速度成分と並進寄与低周波加 速度成分を加算し、この加算結果を並進寄与加速度と推定する。こうして、並進寄与 加速度及び回転寄与加速度が推定されると、並進回転加速度推定部 45Cは、一時 記憶領域 24に、推定された並進寄与加速度及び回転寄与加速度、並びに姿勢角 を 1組の通知用データとして格納する。こうして、ステップ S26Cが終了すると、図 12 に示されるように、処理はステップ S21へ移行する。以後、計測データ処理の停止要 求をアプリケーション 33から受けるまで、計測データ処理プログラム 35Cは、ステップ S21〜S26Cを繰り返す。
[0137] 上記のステップ S11Cの計測データ処理における所望の時点で、アプリケーション 3 3が、アプリケーション識別子をパラメータとするデータ要求を発行すると、第 1及び第 2実施形態の場合と同様に、計測データ処理プログラム 35Cでは、計測データ処理 制御部 49が、このデータ要求を受け、アプリケーション 33からの指令を解析し、アブ リケーシヨン 33からのデータ要求であることを認識する。そして、計測データ処理制 御部 49は、データ通知部 46に対して、アプリケーション 33への通知用データの通知 指令を行う。この指令を受けたデータ通知部 46は、ステップ S12 (j = l, 2, · · ·)のデ ータ通知処理を行う(図 5参照)。
[0138] 以上のようにして、計測データ処理が行われることにより、アプリケーション 33は、所 望の時点におけるデータ要求の発行により、最新の並進寄与加速度及び回転寄与 加速度の推定値、並びにこれらの推定値に対応する時点における姿勢角を取得す ることがでさる。
[0139] また、本実施形態においても、第 1及び第 2実施形態の場合と同様にして、ステップ S11Cの計測データ処理は、アプリケーション 33が、所望の時に、データ収集停止要 求を発行することにより終了する。また、アプリケーション 33がセンサ動作停止要求を 発行すると、これを受けた計測データ処理プログラム 35C (より詳しくは、計測データ 処理制御部 49)は、センサユニット 25へ向けて、センサ動作停止指令を送る。この結 果、センサユニット 25の動作が停止する。
[0140] 以上説明したように、本第 3実施形態では、データ収集部 41が、携帯電話装置 10 において固有に定義される互いに直交する X軸、 Y軸及び Z軸の内の X軸及び Y軸 それぞれに沿った方向の加速度の計測結果、並びに姿勢角の計測結果を収集する 。引き続き、加速度周波数分離部 42が、加速度の計測結果それぞれにおける低周 波加速度成分と高周波加速度成分とを分離する。また、姿勢角変化率導出部 43C が姿勢角の時間変化率の低周波成分と高周波成分とを分離する。そして、並進回転 加速度推定部 45Cが、高周波加速度成分、低周波加速度成分、姿勢角、並びに姿 勢角の時間変化率の高周波成分及び低周波成分に基づいて、並進寄与加速度と 回転寄与加速度とを推定する。
[0141] ここで、並進回転加速度推定部 45Cにおける回転寄与高周波成分推定部 55が、 回転寄与高周波加速度成分を推定するとともに、並進回転加速度推定部 45Cにお ける並進寄与高周波成分推定部 56が、並進寄与高低周波加速度成分を推定する。 また、並進回転加速度推定部 45Cにおける回転寄与低周波成分推定部 53が、回転 寄与低周波加速度成分を推定するとともに、並進回転加速度推定部 45Cにおける 並進寄与低周波成分推定部 54が、並進寄与低周波加速度成分を推定する。
[0142] したがって、携帯電話装置 10の並進に寄与している並進寄与加速度と回転に寄与 している回転寄与加速度とを精度を向上して推定することができる。
[0143] 《第 4実施形態》
以下、本発明の第 4実施形態を説明する。なお、本実施形態の説明においては、 上述した第 1実施形態と同一又は同等の要素については同一の符号を付し、重複す る説明を省略する。
[0144] 本実施形態の携帯電話装置 10は、第 1実施形態と比べて、計測データ処理プログ ラム 35Aに代えて、計測データ処理プログラム 35Dを備える点が相違している。以下 、この相違点に主に着目して説明する。 [0145] 計測データ処理プログラム 35Dは、図 14に示されるように、第 1実施形態の計測デ ータ処理プログラム 35Aと比べて、加速度周波数分離部 42に代えて、加速度周波 数分離部 42Dを備えるとともに、並進回転加速度推定部 45Aに代えて、並進回転加 速度推定部 45Dを備える点のみが相違する。
[0146] 加速度周波数分離部 42Dは、所定個(例えば、 20個)の最近のデータを読み出し て、移動平均値を算出することにより、直流成分と評価できる低周波加速度成分を求 める。そして、算出された低周波加速度成分を有すると推定される時点における加速 度データ力 当該低周波加速度成分を差し引くことにより、高周波加速度成分を求 める。
[0147] 並進回転加速度推定部 45Dは、高周波成分推定部 51Dと、低周波成分推定部 5 2Dとを備えている。ここで、低周波成分推定部 52Dは、低周波加速度成分には重力 加速度の成分のみが含まれると推定する。これは、携帯電話装置 10の手動による運 動の場合には、回転加速度及び並進加速度の双方ともに、僅かではあっても刻々変 化することが通常であり、直流成分と評価できる加速度成分は、定常的に携帯電話 装置等に作用している重力に由来すると考えられることによる。
[0148] また、高周波成分推定部 51Dは、姿勢角の時間変化率が所定値未満である場合 には、高周波加速度成分が、並進寄与加速度のみを含むと推定する。これは、姿勢 角の時間変化率が所定値未満である場合には、高周波加速度成分には、並進寄与 加速度が多く含まれ、回転寄与加速度が少しし力含まれないと考えられることによる
[0149] 一方、高周波成分推定部 51Dは、姿勢角の時間変化率が所定値以上である場合 には、高周波加速度成分は回転寄与加速度のみを含むと推定する。これは、姿勢角 の時間変化率が所定値以上である場合には、高周波加速度成分には、回転寄与加 速度が多く含まれ、並進寄与加速度は少しし力含まれな 、と考えられることによる。
[0150] 上記のように構成された本実施形態の携帯電話装置 10では、第 1実施形態の場合 におけるステップ S 11 Aの計測データ処理に代えて、ステップ S11D (図 15参照)の 計測データ処理を行う点が相違している。このステップ S11Dでは、図 15に示される ように、まず、第 1実施形態と同様にして、ステップ S21〜S23Dにおいて、データ収 集部 41が、計測生データを収集し、計測物理量データを算出して一時記憶領域 24 に格納することにより、所定数以上の組の計測物理量データが揃った力否かが判定 される。なお、本実施形態においては、当該所定数が第 1実施形態の場合よりも多く 、後述する直流成分と評価できる加速度変化の低周波成分を抽出が可能な値が採 用される。この所定数は、経験や実験や評価の結果に基づいて定められる。
[0151] 引き続き、ステップ 24Dにおいて、計測データ処理制御部 49からの指令を受けた 加速度周波数分離部 42が、計測物理量データにおける加速度データに基づいて、 高周波加速度成分と低周波加速度成分とを分離する。かかる高周波加速度成分と 低周波加速度成分との分離に際して、加速度周波数分離部 42Dは、第 1実施形態 の場合と同様にして、移動平均を算出することにより、直流成分と評価できる低周波 加速度成分を求める。次に、加速度周波数分離部 42Dは、算出された移動平均値 が低周波加速度成分とされる時点の加速度データから当該低周波加速度を差し引く ことにより、当該時点における高周波加速度成分を求める。
[0152] こうして、当該時点の X軸方向加速度データ及び Y軸方向加速度データそれぞれ における高周波加速度成分と低周波加速度成分とが分離される。加速度周波数分 離部 42Dは、分離された高周波加速度成分と低周波加速度成分とを並進回転加速 度推定部 45Dへ送る。
[0153] 次に、ステップ S25Dにおいて、第 1実施形態におけるステップ 25Aの場合と同様 にして、姿勢角変化率導出部 43が、計測物理量データにおける姿勢角データに基 づいて、上記の低周波加速度成分が算出された時点における姿勢角それぞれの時 間変化率を導出する。
[0154] 次いで、ステップ S26Dにおいて、並進回転加速度推定部 45D力 加速度周波数 分離部 42からの高周波加速度成分及び低周波加速度成分、並びに姿勢角変化率 導出部 43Dからの姿勢角及び姿勢角の時間変化率に基づいて、並進寄与加速度と 回転寄与加速度とを推定する。力かる推定に際しては、図 16に示されるように、ステ ップ S61において、並進回転加速度推定部 45Aの低周波成分推定部 52Dが、低周 波加速度成分は、全てが重力加速度成分であり、並進寄与加速度又は回転寄与加 速度を含まないと推定する。 [0155] 引き続き、ステップ S62において、並進回転加速度推定部 45Aの高周波成分推定 部 51Dが、姿勢角変化率導出部 43から受けた姿勢角の時間変化率の少なくとも 1 つが所定値以上力否かを判定する。この判定結果が肯定的であった場合には、処 理はステップ S63へ移行する。このステップ S63では、高周波成分推定部 51Dが、 高周波加速度成分の全てが回転寄与加速度であると推定する。
[0156] 一方、ステップ S62における判定結果は否定的であった場合には、処理はステップ S64へ移行する。このステップ S64では、高周波成分推定部 51Dが、高周波加速度 成分の全てが並進寄与加速度であると推定する。
[0157] こうして、並進寄与加速度及び回転寄与加速度が推定されると、並進回転加速度 推定部 45Dは、一時記憶領域 24に、推定された並進寄与加速度及び回転寄与カロ 速度、並びに姿勢角を 1組の通知用データとして格納する。こうして、ステップ S26D が終了すると、図 15に示されるように、処理はステップ S21へ移行する。以後、計測 データ処理の停止要求をアプリケーション 33から受けるまで、計測データ処理プログ ラム 35Dは、ステップ S21〜S26Dを繰り返す。
[0158] 上記のステップ S11Dの計測データ処理における所望の時点で、アプリケーション 33が、アプリケーション識別子をパラメータとするデータ要求を発行すると、第 1実施 形態の場合と同様に、計測データ処理プログラム 35Dでは、計測データ処理制御部 49力 このデータ要求を受け、アプリケーション 33からの指令を解析し、アプリケーシ ヨン 33からのデータ要求であることを認識する。そして、計測データ処理制御部 49は 、データ通知部 46に対して、アプリケーション 33への通知用データの通知指令を行 う。この指令を受けたデータ通知部 46は、ステップ S12 (j = l, 2, · · ·)のデータ通知 処理を行う(図 5参照)。
[0159] 以上のようにして、計測データ処理が行われることにより、アプリケーション 33は、所 望の時点におけるデータ要求の発行により、最新の並進寄与加速度及び回転寄与 加速度の推定値、並びにこれらの推定値に対応する時点における姿勢角を取得す ることがでさる。
[0160] また、本実施形態においても、第 1実施形態の場合と同様にして、ステップ S11Dの 計測データ処理は、アプリケーション 33が、所望の時に、データ収集停止要求を発 行することにより終了する。また、アプリケーション 33がセンサ動作停止要求を発行す ると、これを受けた計測データ処理プログラム 35D (より詳しくは、計測データ処理制 御部 49)は、センサユニット 25へ向けて、センサ動作停止指令を送る。この結果、セ ンサユニット 25の動作が停止する。
[0161] 以上説明したように、本第 4実施形態では、データ収集部 41が、携帯電話装置 10 において固有に定義される互いに直交する X軸、 Y軸及び Z軸の内の X軸及び Y軸 それぞれに沿った方向の加速度の計測結果、並びに姿勢角の計測結果を収集する 。引き続き、加速度周波数分離部 42Dが、加速度の計測結果それぞれにおける直 流成分と評価できる低周波加速度成分と、それ以外の高周波加速度成分とを分離 する。また、姿勢角変化率導出部 43が姿勢角の時間変化率の低周波成分を導出す る。そして、並進回転加速度推定部 45Dが、高周波加速度成分、低周波加速度成 分、姿勢角及び姿勢角の時間変化率の低周波成分に基づいて、並進寄与加速度と 回転寄与加速度とを推定する。
[0162] ここで、並進回転加速度推定部 45Dにおける低周波成分推定部 52Dが、低周波 加速度成分は、重力加速度成分のみを含むと推定する。また、並進回転加速度推 定部 45Dにおける高周波成分推定部 55Dが、姿勢角の時間変化率の絶対値が所 定の値以上である力否かを判定する。この判定が否定的であった場合には、高周波 成分推定部 51Dが、高周波加速度成分には、回転寄与加速度の成分は含まれず、 並進寄与加速度の成分のみを含むと推定する。一方、当該判定が肯定的であった 場合には、高周波成分推定部 51Dが、高周波加速度成分には、並進寄与加速度の 成分は含まれず、回転寄与加速度の成分のみを含むと推定する。
[0163] したがって、携帯電話装置 10の並進に寄与している並進寄与加速度と回転に寄与 している回転寄与加速度とを精度を向上して推定することができる。
[0164] 《第 5実施形態》
以下、本発明の第 5実施形態を説明する。なお、本実施形態の説明においては、 上述した第 4実施形態と同一又は同等の要素については同一の符号を付し、重複す る説明を省略する。
[0165] 本実施形態の携帯電話装置 10は、第 4実施形態と比べて、計測データ処理プログ ラム 35Dに代えて、計測データ処理プログラム 35Eを備える点が相違している。以下 、この相違点に主に着目して説明する。
[0166] 計測データ処理プログラム 35Eは、図 17に示されるように、第 4実施形態の計測デ ータ処理プログラム 35Dと比べて、並進回転加速度推定部 45Dに代えて、並進回転 加速度推定部 45Eを備える点のみが相違する。この並進回転加速度推定部 45Eは 、並進回転加速度推定部 45Dと比べて、高周波成分推定部 51Dに代えて、姿勢角 時間変化率に基づいて、回転寄与加速度成分を推定する回転寄与成分推定部 55 Eと、高周波加速度成分及び回転寄与加速度成分に基づいて、並進寄与加速度成 分を推定する並進寄与成分推定部 56Eとを備える点のみが相違している。
[0167] 上記のように構成された本実施形態の携帯電話装置 10では、第 4実施形態の場合 におけるステップ S 11Dの計測データ処理に代えて、ステップ SI 1E (図 18参照)の 計測データ処理を行う点が相違している。このステップ S11Eでは、図 18に示される ように、まず、第 4実施形態と同様にして、ステップ S21〜S24Dにおいて、データ収 集部 41が、計測生データを収集し、計測物理量データを算出して一時記憶領域 24 に格納することにより、所定数以上の組の計測物理量データが揃うと、加速度周波数 分離部 42D力 計測物理量データにおける加速度データに基づいて、高周波加速 度成分と低周波加速度成分 (直流成分と評価できる成分)とを分離する。そして、ステ ップ 25Dにお 、て、計測データ処理制御部 49からの指令を受けた姿勢角変化率導 出部 43が、計測物理量データにおける姿勢角データに基づいて、計測物理量デー タにおける姿勢角データに基づいて、上記の低周波加速度成分が算出された時点 における姿勢角それぞれの時間変化率を導出する。
[0168] 次 、で、ステップ S26Eにお 、て、並進回転加速度推定部 45E力 加速度周波数 分離部 42Dからの高周波加速度成分及び低周波加速度成分、並びに姿勢角変化 率導出部 43からの姿勢角及び姿勢角の時間変化率の低周波成分に基づいて、並 進寄与加速度と回転寄与加速度とを推定する。力かる推定に際しては、図 19に示さ れるように、ステップ S71において、並進回転加速度推定部 45Eの低周波成分推定 部 52Dが、第 4実施形態の場合と同様にして、低周波加速度成分は、全てが重力加 速度成分であり、並進寄与加速度又は回転寄与加速度を含まないと推定する。 [0169] 次に、ステップ S72において、並進回転加速度推定部 45Eの回転寄与加速度推 定部 55Eが、姿勢角変化率導出部 43から受けた姿勢角の変化率、並びに、予め知 られている携帯電話装置 10の慣性モーメント及び重心位置に基づいて、回転に寄 与している加速度成分である回転寄与加速度成分を推定する。引き続き、並進回転 加速度推定部 45Eの並進寄与加速度推定部 56Eが、ステップ S73において、回転 寄与加速度推定部 55Eから受けた回転寄与低周波加速度成分を差し引く。こうして 算出されたものを、並進寄与加速度推定部 56Eは、並進に寄与している加速度成分 である並進寄与低周波加速度成分であると推定する。
[0170] こうして、並進寄与加速度及び回転寄与加速度が推定されると、並進回転加速度 推定部 45Eは、一時記憶領域 24に、推定された並進寄与加速度及び回転寄与カロ 速度、並びに姿勢角を 1組の通知用データとして格納する。こうして、ステップ S26E が終了すると、図 18に示されるように、処理はステップ S21へ移行する。以後、計測 データ処理の停止要求をアプリケーション 33から受けるまで、計測データ処理プログ ラム 35Eは、ステップ S21〜S26Eを繰り返す。
[0171] 上記のステップ S11Eの計測データ処理における所望の時点で、アプリケーション 3 3が、アプリケーション識別子をパラメータとするデータ要求を発行すると、第 1実施形 態の場合と同様に、計測データ処理プログラム 35Eでは、計測データ処理制御部 49 力 このデータ要求を受け、アプリケーション 33からの指令を解析し、アプリケーショ ン 33からのデータ要求であることを認識する。そして、計測データ処理制御部 49は、 データ通知部 46に対して、アプリケーション 33への通知用データの通知指令を行う 。この指令を受けたデータ通知部 46は、ステップ S12 (j = l, 2, · · ·)のデータ通知処 理を行う(図 5参照)。
[0172] 以上のようにして、計測データ処理が行われることにより、アプリケーション 33は、所 望の時点におけるデータ要求の発行により、最新の並進寄与加速度及び回転寄与 加速度の推定値、並びにこれらの推定値に対応する時点における姿勢角を取得す ることがでさる。
[0173] また、本実施形態においても、第 1実施形態の場合と同様にして、ステップ S11Eの 計測データ処理は、アプリケーション 33が、所望の時に、データ収集停止要求を発 行することにより終了する。また、アプリケーション 33がセンサ動作停止要求を発行す ると、これを受けた計測データ処理プログラム 35E (より詳しくは、計測データ処理制 御部 49)は、センサユニット 25へ向けて、センサ動作停止指令を送る。この結果、セ ンサユニット 25の動作が停止する。
[0174] 以上説明したように、本第 5実施形態では、データ収集部 41が、携帯電話装置 10 において固有に定義される互いに直交する X軸、 Y軸及び Z軸の内の X軸及び Y軸 それぞれに沿った方向の加速度の計測結果、並びに姿勢角の計測結果を収集する 。引き続き、加速度周波数分離部 42Dが、加速度の計測結果それぞれにおける直 流成分と評価できる低周波加速度成分と、高周波加速度成分とを分離する。また、 姿勢角変化率導出部 43が姿勢角の時間変化率を導出する。そして、並進回転加速 度推定部 45Eが、高周波加速度成分、低周波加速度成分、姿勢角及び姿勢角の時 間変化率の低周波成分に基づ!/、て、並進寄与加速度と回転寄与加速度とを推定す る。
[0175] ここで、並進回転加速度推定部 45Eにおける低周波成分推定部 52Dが、低周波 加速度成分は、重力加速度成分のみを含むと推定する。また、並進回転加速度推 定部 45Eにおける回転寄与成分推定部 55E力 回転寄与加速度成分を推定すると ともに、並進回転加速度推定部 45Eにおける並進寄与成分推定部 56E力 並進寄 与加速度成分を推定する。
[0176] したがって、携帯電話装置 10の並進に寄与している並進寄与加速度と回転に寄与 している回転寄与加速度とを精度を向上して推定することができる。
[0177] なお、上記の第 1〜第 5の各実施形態では、センサユニット 25を携帯電話装置 10 の内部に実装した。これに対して、センサユニット 25を携帯電話装置 10の外部に配 置し、携帯電話装置 10の不図示の外部機器接続用インタフェースポートを介して、 センサユニット 25と携帯電話装置 10とを接続するようにしてもょ 、。
[0178] また、上記の第 1〜第 5の各実施形態では、センサユニット 25のセンサ部 26を、ピ ツチ角 Θ 、ロール角 Θ 、ョ一角 θ 、X方向加速度 α 及び Υ方向加速度 α を検知
X Υ Ζ X Υ
する 5軸センサとした力 更に Ζ方向加速度 α の検知を行う 6軸センサとすることもで
Ζ
きる。さらに、姿勢角や加速度以外の物理量を検出するセンサとすることもできる。 [0179] また、上記の第 1〜第 5の各実施形態では、操作部 12のマトリクス状のキー配列に おける行方向を X軸とし、列方向を Y軸方向とするとともに、 XY面が水平面と平行で あり、 +Y方向が真南方向である姿勢を基準姿勢としたが、他の姿勢を基準姿勢とす ることもできる。この場合には、基準姿勢と鉛直方向との関係を予め計測したり、計測 データのそれぞれが得られた時点における姿勢と鉛直方向との関係を求めたりする ことが必要となる。
[0180] また、上記の第 1〜第 5の各実施形態では、キーが配列された操作部 12と表示部 1 3との位置関係が固定的ないわゆるストレート型の携帯電話装置とした。これに対して 、いわゆるクラムシェル型ゃレボルバ型のように、操作部と表示部との位置関係が可 変な携帯電話装置の場合には、センサユニット 25を操作部側に配設してもよいし、 表示部側に配設してもよい。また、センサユニット 25における計測の基準となる軸方 向は、センサユニット 25の配設位置や、携帯電話の機種に対応して決めることができ る。
[0181] また、上記の第 1〜第 5の各実施形態では、携帯電話装置に対して本発明を適用 したが、携帯型ゲーム機、カーナビゲーシヨン装置、 PDA (Personal Digital Assistanc e)等の他の種類の移動端末装置に対しても本発明が適用できるのは、勿論である。 産業上の利用可能性
[0182] 以上説明したように、本発明の加速度推定方法及び加速度推定装置は、移動物体 の並進に寄与している並進寄与加速度と回転に寄与している回転寄与加速度とを推 定する場合に適用することができる。また、本発明の移動端末装置は、自身の並進 に寄与している並進寄与加速度と、自身の回転に寄与している回転寄与加速度とを 推定し、推定結果を利用する移動端末装置に適用することができる。

Claims

請求の範囲
[1] 移動物体において固有に定義される第 1軸、前記第 1軸と直交する第 2軸、並びに 前記第 1軸及び前記第 2軸と直交する第 3軸の内の少なくとも 2つの軸方向に沿った 加速度の計測結果における前記移動物体の並進に寄与している並進寄与加速度と 、回転に寄与して 、る回転寄与加速度とを推定する加速度推定方法であって、 前記少なくとも 2つの軸に沿った加速度の計測結果と、前記移動物体の基準姿勢 力 の回転角である姿勢角の計測結果とを収集する計測結果収集工程と; 前記少なくとも 2つの軸に沿った加速度の計測結果における低周波加速度成分と 高周波加速度成分とを分離する加速度周波数分離工程と;
前記姿勢角の時間変化率を導出する姿勢角変化率導出工程と;
前記高周波加速度成分、前記低周波加速度成分、前記姿勢角及び前記姿勢角の 時間変化率に基づ!、て、前記並進寄与加速度及び前記回転寄与加速度を推定す る並進回転加速度推定工程と;を備える加速度推定方法。
[2] 前記並進回転加速度推定工程は、
前記高周波加速度成分は、前記回転寄与加速度の成分を含まず、前記並進寄 与加速度の成分のみを含むと推定する高周波成分推定工程と;
前記姿勢角の時間変化率の絶対値が所定の値未満の場合には、前記低周波加 速度成分力 重力加速度を差し引 ヽたものは、前記回転寄与加速度の成分を含ま ず、前記並進寄与加速度の成分のみを含むと推定し、前記姿勢角の時間変化率の 絶対値が前記所定の値以上の場合には、前記低周波加速度成分から重力加速度 を差し引いたものは、前記並進寄与加速度の成分を含まず、前記回転寄与加速度 の成分のみを含むと推定する低周波成分推定工程と;を備えることを特徴とする請求 項 1に記載の加速度推定方法。
[3] 前記姿勢角変化率導出工程では、前記姿勢角の時間変化率の低周波成分を導 出し、
前記並進加速度推定工程は、
前記高周波加速度成分は、前記回転寄与加速度の成分を含まず、前記並進寄 与加速度の成分のみを含むと推定する高周波成分推定工程と; 前記姿勢角の時間変化率の低周波成分に基づいて、前記低周波加速度成分に おける前記回転寄与加速度の成分である回転寄与低周波加速度成分を推定する回 転寄与低周波成分推定工程と;
前記低周波加速度成分から前記回転寄与低周波加速度成分及び重力加速度を 差し引いたものが、前記並進寄与加速度の低周波成分である並進寄与低周波加速 度成分であると推定する並進寄与低周波成分推定工程と;を備えることを特徴とする 請求項 1に記載の加速度推定方法。
[4] 前記姿勢角変化率導出工程では、前記姿勢角の時間変化率における低周波成分 と高周波成分を分離し、
前記並進加速度推定工程は、
前記姿勢角の時間変化率の高周波成分に基づいて、前記高周波加速度成分に おける前記回転寄与加速度の高周波成分である回転寄与高周波加速度成分を推 定する回転寄与高周波成分推定工程と;
前記高周波加速度成分から前記回転寄与高周波加速度成分を差し引いたもの 力 前記並進寄与加速度の高周波成分である並進寄与高周波加速度成分であると 推定する並進寄与高周波成分推定工程と;
前記姿勢角の時間変化率の低周波成分に基づいて、前記低周波加速度成分に おける前記回転寄与加速度の低周波成分である回転寄与低周波加速度成分を推 定する回転寄与低周波成分推定工程と;
前記低周波加速度成分から前記回転寄与低周波加速度成分及び重力加速度を 差し引いたものが、前記並進寄与加速度の低周波成分である並進寄与低周波加速 度成分であると推定する並進寄与低周波成分推定工程と;を備えることを特徴とする 請求項 1に記載の加速度推定方法。
[5] 前記加速度周波数分離工程では、前記低周波加速度成分として、直流成分と評 価できる成分を抽出し、
前記並進回転加速度推定工程は、
前記低周波加速度成分が重力加速度にのみ由来すると推定する低周波成分推 定工程と; 前記姿勢角の時間変化率の絶対値が所定の値未満の場合には、前記高周波加 速度成分は、回転に寄与する回転寄与加速度の成分を含まず、前記並進寄与加速 度の成分のみを含むと推定し、前記姿勢角の時間変化率の絶対値が前記所定の値 以上の場合には、前記高周波成分は、前記並進寄与加速度の成分を含まず、前記 回転寄与加速度の成分のみを含むと推定する高周波成分推定工程と;を備えること を特徴とする請求項 1に記載の加速度推定方法。
[6] 前記加速度周波数分離工程では、前記低周波加速度成分として、直流成分と評 価できる成分を抽出し、
前記並進回転加速度推定工程は、
前記低周波加速度成分が重力加速度にのみ由来すると推定する低周波成分推 定工程と;
前記姿勢角の時間変化率に基づいて、前記高周波加速度成分における回転に 寄与する回転寄与加速度成分を推定する回転寄与成分推定工程と;
前記高周波加速度成分から前記回転寄与加速度成分を差し引いたものが、前記 並進寄与加速度成分であると推定する並進寄与成分推定工程と;を備えることを特 徴とする請求項 1に記載の加速度推定方法。
[7] 移動物体において固有に定義される第 1軸、前記第 1軸と直交する第 2軸、並びに 前記第 1軸及び前記第 2軸と直交する第 3軸の内の少なくとも 2つの軸方向に沿った 加速度の計測結果における前記移動物体の並進に寄与している並進寄与加速度と 、回転に寄与して 、る回転寄与加速度とを推定する加速度推定装置であって、 前記少なくとも 2つの軸に沿った加速度の計測結果と、前記移動物体の基準姿勢 からの回転角である姿勢角の計測結果とを収集する計測結果収集手段と; 前記少なくとも 2つの軸に沿った加速度の計測結果における低周波加速度成分と 高周波加速度成分とを分離する加速度周波数分離手段と;
前記姿勢角の時間変化率を導出する姿勢角変化率導出手段と;
前記高周波加速度成分、前記低周波加速度成分、前記姿勢角及び前記姿勢角の 時間変化率に基づ!、て、前記並進寄与加速度及び前記回転寄与加速度を推定す る並進回転加速度推定手段と;を備える加速度推定装置。
[8] 前記並進回転加速度推定手段は、
前記高周波加速度成分が、回転に寄与する回転寄与加速度の成分を含まず、 前記並進寄与加速度の成分のみを含むと推定する高周波成分推定手段と;
前記姿勢角の時間変化率の絶対値が所定の値未満の場合には、前記低周波加 速度成分力 重力加速度を差し引 ヽたものは、前記回転寄与加速度の成分を含ま ず、前記並進寄与加速度の成分のみを含むと推定し、前記姿勢角の時間変化率の 絶対値が前記所定の値以上の場合には、前記低周波加速度成分から重力加速度 を差し引いたものは、前記並進寄与加速度の成分を含まず、前記回転寄与加速度 の成分のみを含むと推定する低周波成分推定手段と;を備えることを特徴とする請求 項 7に記載の加速度推定装置。
[9] 前記姿勢角変化率導出手段は、前記姿勢角の時間変化率の低周波成分を導出し 前記並進回転加速度推定手段は、
前記高周波加速度成分が、回転に寄与する回転寄与加速度の成分を含まず、 前記並進寄与加速度の成分のみを含むと推定する高周波成分推定手段と;
前記姿勢角の時間変化率の低周波成分に基づいて、前記低周波加速度成分に おける前記回転寄与加速度の成分である回転寄与低周波加速度成分を推定する回 転寄与低周波成分推定手段と;
前記低周波加速度成分から前記回転寄与低周波加速度成分及び重力加速度を 差し引いたものが、前記並進寄与加速度の低周波成分である並進寄与低周波加速 度成分であると推定する並進寄与低周波成分推定手段と;を備えることを特徴とする 請求項 7に記載の加速度推定装置。
[10] 前記姿勢角変化率導出手段は、前記姿勢角の時間変化率における低周波成分と 高周波成分を分離し、
前記並進回転加速度推定手段は、
前記姿勢角の時間変化率の高周波成分に基づいて、前記高周波加速度成分に おける前記回転寄与加速度の高周波成分である回転寄与高周波加速度成分を推 定する回転寄与高周波成分推定手段と; 前記高周波加速度成分から前記回転寄与高周波加速度成分を差し引いたもの 力 前記並進寄与加速度の高周波成分である並進寄与高周波加速度成分であると 推定する並進寄与高周波成分推定手段と;
前記姿勢角の時間変化率の低周波成分に基づいて、前記低周波加速度成分に おける前記回転寄与加速度の低周波成分である回転寄与低周波加速度成分を推 定する回転寄与低周波成分推定手段と;
前記低周波加速度成分から前記回転寄与低周波加速度成分及び重力加速度を 差し引いたものが、前記並進寄与加速度の低周波成分である並進寄与低周波加速 度成分であると推定する並進寄与低周波成分推定手段と;を備えることを特徴とする 請求項 7に記載の加速度推定装置。
[11] 前記加速度周波数分離手段は、前記低周波加速度成分として、直流成分と評価 できる成分を抽出し、
前記並進回転加速度推定手段は、
前記低周波加速度成分が重力加速度にのみ由来すると推定する低周波成分推 定手段と;
前記姿勢角の時間変化率の絶対値が所定の値未満の場合には、前記高周波加 速度成分は、回転に寄与する回転寄与加速度の成分を含まず、前記並進寄与加速 度の成分のみを含むと推定し、前記姿勢角の時間変化率の絶対値が前記所定の値 以上の場合には、前記高周波成分は、前記並進寄与加速度の成分を含まず、前記 回転寄与加速度の成分のみを含むと推定する高周波成分推定手段と;を備えること を特徴とする請求項 7に記載の加速度推定装置。
[12] 前記加速度周波数分離手段は、前記低周波加速度成分として、直流成分と評価 できる成分を抽出し、
前記並進回転加速度推定手段は、
前記低周波加速度成分が重力加速度にのみ由来すると推定する低周波成分推 定手段と;
前記姿勢角の時間変化率に基づいて、前記高周波加速度成分における回転に 寄与する回転寄与加速度成分を推定する回転寄与成分推定手段と; 前記高周波加速度成分から前記回転寄与加速度成分を差し引いたものが、前記 並進寄与加速度成分であると推定する並進寄与成分推定手段と;を備えることを特 徴とする請求項 7に記載の加速度推定装置。
移動しつつ動作が可能な移動端末装置において、
請求項 7に記載の加速度推定装置を備える、ことを特徴とする移動端末装置。
PCT/JP2005/016112 2004-09-07 2005-09-02 加速度推定方法、加速度推定装置及び移動端末装置 WO2006028013A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2006527204A JPWO2006028013A1 (ja) 2004-09-07 2005-09-02 加速度推定方法、加速度推定装置及び移動端末装置
EP05781378A EP1788397A4 (en) 2004-09-07 2005-09-02 WATER POUCH FOR ULTRASONIC PROBE USED IN ULTRASONIC DIAGNOSTIC DEVICE
US11/706,239 US8521467B2 (en) 2004-09-07 2007-02-15 Method for acceleration estimation, acceleration estimating apparatus, and mobile terminal

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2004259613 2004-09-07
JP2004-259613 2004-09-07

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US11/706,239 Continuation US8521467B2 (en) 2004-09-07 2007-02-15 Method for acceleration estimation, acceleration estimating apparatus, and mobile terminal

Publications (1)

Publication Number Publication Date
WO2006028013A1 true WO2006028013A1 (ja) 2006-03-16

Family

ID=36036298

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2005/016112 WO2006028013A1 (ja) 2004-09-07 2005-09-02 加速度推定方法、加速度推定装置及び移動端末装置

Country Status (4)

Country Link
US (1) US8521467B2 (ja)
EP (1) EP1788397A4 (ja)
JP (1) JPWO2006028013A1 (ja)
WO (1) WO2006028013A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011047879A (ja) * 2009-08-28 2011-03-10 Oki Semiconductor Co Ltd 運動検出装置、電子機器、運動検出方法及びプログラム

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006028013A1 (ja) 2004-09-07 2006-03-16 Vodafone K.K. 加速度推定方法、加速度推定装置及び移動端末装置
JP4321554B2 (ja) * 2006-06-23 2009-08-26 トヨタ自動車株式会社 姿勢角検出装置と姿勢角検出方法
CN101711364A (zh) * 2007-04-13 2010-05-19 领航网络有限公司 确定移动物体的旋转半径的力传感设备和方法
US9047850B1 (en) 2007-10-17 2015-06-02 David Wiley Beaty Electric instrument music control device with magnetic displacement sensors
US7923623B1 (en) * 2007-10-17 2011-04-12 David Beaty Electric instrument music control device with multi-axis position sensors
US8370106B2 (en) * 2008-07-07 2013-02-05 Keynetik, Inc. Spatially aware inference logic
JP5736106B2 (ja) * 2009-05-19 2015-06-17 古野電気株式会社 移動状態検出装置
US8284847B2 (en) * 2010-05-03 2012-10-09 Microsoft Corporation Detecting motion for a multifunction sensor device
EP2392991A1 (en) * 2010-06-02 2011-12-07 Fraunhofer-Gesellschaft zur Förderung der Angewandten Forschung e.V. Hand-held pointing device, software cursor control system and method for controlling a movement of a software cursor
KR102336879B1 (ko) * 2015-05-20 2021-12-08 삼성전자주식회사 화면을 표시하는 전자 장치, 그 제어 방법

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS625185A (ja) * 1985-07-02 1987-01-12 Nissan Motor Co Ltd 横加速度検出装置
JPH0728591A (ja) * 1993-05-13 1995-01-31 Toshiba Corp 空間操作マウスシステム及び空間動作パターン入力方法
JPH0943269A (ja) * 1995-07-28 1997-02-14 Omron Corp 加速度トランスデューサ
JPH1021000A (ja) * 1996-06-28 1998-01-23 Sumitomo Metal Ind Ltd 信号入力装置
JP2002169645A (ja) * 2000-12-01 2002-06-14 Toshiba Corp 移動通信端末
JP2004102738A (ja) * 2002-09-10 2004-04-02 Sony Corp 電子機器

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6270766A (ja) * 1985-09-25 1987-04-01 Nissan Motor Co Ltd 車両揺動運動検出装置
US5452698A (en) * 1990-05-07 1995-09-26 Robert Bosch Gmbh Device for suppressing discontinuous motion of a moving motor vehicle
JP2878498B2 (ja) * 1991-09-13 1999-04-05 株式会社日立製作所 角加速度検出装置
US5225637A (en) * 1991-10-01 1993-07-06 Kurta Corporation Position resolving system
US5598187A (en) 1993-05-13 1997-01-28 Kabushiki Kaisha Toshiba Spatial motion pattern input system and input method
US5627465A (en) * 1995-10-25 1997-05-06 Honeywell Inc. Rotational position sensor with mechanical adjustment of offset and gain signals
US5902968A (en) * 1996-02-20 1999-05-11 Ricoh Company, Ltd. Pen-shaped handwriting input apparatus using accelerometers and gyroscopes and an associated operational device for determining pen movement
JPH10224509A (ja) 1997-02-03 1998-08-21 Matsushita Graphic Commun Syst Inc 通信端末装置
US6360602B1 (en) * 1999-07-29 2002-03-26 Litton Systems, Inc. Method and apparatus reducing output noise in a digitally rebalanced accelerometer
GB2358108A (en) * 1999-11-29 2001-07-11 Nokia Mobile Phones Ltd Controlling a hand-held communication device
TW577975B (en) * 2000-07-25 2004-03-01 American Gnc Corp Core inertial measurement unit
US6725173B2 (en) * 2000-09-02 2004-04-20 American Gnc Corporation Digital signal processing method and system thereof for precision orientation measurements
JP4198875B2 (ja) * 2000-11-30 2008-12-17 株式会社東芝 移動通信端末
JP3775258B2 (ja) * 2001-07-31 2006-05-17 株式会社デンソー フィルタリング方法及びフィルタ機能を有するa/d変換装置
JP4032681B2 (ja) * 2001-08-27 2008-01-16 株式会社デンソー 同期検波方法及び装置並びにセンサ信号検出装置
US20050212760A1 (en) * 2004-03-23 2005-09-29 Marvit David L Gesture based user interface supporting preexisting symbols
WO2006028013A1 (ja) 2004-09-07 2006-03-16 Vodafone K.K. 加速度推定方法、加速度推定装置及び移動端末装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS625185A (ja) * 1985-07-02 1987-01-12 Nissan Motor Co Ltd 横加速度検出装置
JPH0728591A (ja) * 1993-05-13 1995-01-31 Toshiba Corp 空間操作マウスシステム及び空間動作パターン入力方法
JPH0943269A (ja) * 1995-07-28 1997-02-14 Omron Corp 加速度トランスデューサ
JPH1021000A (ja) * 1996-06-28 1998-01-23 Sumitomo Metal Ind Ltd 信号入力装置
JP2002169645A (ja) * 2000-12-01 2002-06-14 Toshiba Corp 移動通信端末
JP2004102738A (ja) * 2002-09-10 2004-04-02 Sony Corp 電子機器

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011047879A (ja) * 2009-08-28 2011-03-10 Oki Semiconductor Co Ltd 運動検出装置、電子機器、運動検出方法及びプログラム

Also Published As

Publication number Publication date
JPWO2006028013A1 (ja) 2008-05-08
EP1788397A1 (en) 2007-05-23
US20070146323A1 (en) 2007-06-28
EP1788397A4 (en) 2011-11-02
US8521467B2 (en) 2013-08-27

Similar Documents

Publication Publication Date Title
WO2006028013A1 (ja) 加速度推定方法、加速度推定装置及び移動端末装置
CN108353131B (zh) 获取图像的方法和终端设备
EP1271099B1 (en) Method and arrangement for determining movement
KR100554484B1 (ko) 동작 인식이 가능한 휴대용 단말기 및 동작 인식 방법
CN109041148B (zh) 一种移动终端操作方法及移动终端
CN107742523B (zh) 语音信号处理方法、装置以及移动终端
JP4377863B2 (ja) 運動軌跡情報処理方法及び携帯情報装置
WO2007099908A1 (ja) ウェアラブル端末、および、携帯撮像収音装置、およびそれらを実現する装置、方法、プログラム
JP2012211932A (ja) 音声認識装置及び音声認識方法
JP2001159951A (ja) 情報処理装置および情報処理方法
CN105683711A (zh) 方向估计设备、方向估计系统和估计方向的方法
CN110456923B (zh) 一种姿态传感数据处理方法及电子设备
WO2006051846A1 (ja) 計測方法及び移動情報装置
CN110398197B (zh) 一种电子设备和信息处理方法
KR20120056696A (ko) 사용자 단말의 변위 추정방법 및 장치
JP5289991B2 (ja) 携帯電子機器
CN107367252B (zh) 一种相对位置变化判断方法、装置及计算机可读存储介质
JP2010028674A (ja) 携帯端末
CN108600517B (zh) 一种切换屏幕状态的方法及终端
CN108494946B (zh) 移动终端中电子罗盘的校正方法与装置
CN108062465B (zh) 一种解锁方法及移动终端
JP4845356B2 (ja) 計測データ処理方法、計測データ処理装置及び移動端末装置
CN114279425A (zh) 指南针精度提高方法、终端及计算机可读存储介质
CN115714947A (zh) 一种音频数据的处理方法、装置、移动设备以及音频系统
US20090316882A1 (en) Communication terminal, method for operating communication terminal, and computer program

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KM KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NG NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SM SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): BW GH GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LT LU LV MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

WWE Wipo information: entry into national phase

Ref document number: 2006527204

Country of ref document: JP

121 Ep: the epo has been informed by wipo that ep was designated in this application
REEP Request for entry into the european phase

Ref document number: 2005781378

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2005781378

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 11706239

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

WWP Wipo information: published in national office

Ref document number: 2005781378

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 11706239

Country of ref document: US