WO2007105648A1 - 体動検出装置、体動検出方法および体動検出プログラム - Google Patents

体動検出装置、体動検出方法および体動検出プログラム Download PDF

Info

Publication number
WO2007105648A1
WO2007105648A1 PCT/JP2007/054710 JP2007054710W WO2007105648A1 WO 2007105648 A1 WO2007105648 A1 WO 2007105648A1 JP 2007054710 W JP2007054710 W JP 2007054710W WO 2007105648 A1 WO2007105648 A1 WO 2007105648A1
Authority
WO
WIPO (PCT)
Prior art keywords
peak position
user
body motion
motion detection
peak
Prior art date
Application number
PCT/JP2007/054710
Other languages
English (en)
French (fr)
Inventor
Kenichi Makino
Akane Sano
Motoyuki Takai
Original Assignee
Sony Corporation
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 Sony Corporation filed Critical Sony Corporation
Priority to CN2007800070178A priority Critical patent/CN101394787B/zh
Priority to US12/279,557 priority patent/US7983872B2/en
Priority to EP07738195A priority patent/EP1994883A4/en
Priority to BRPI0708620-2A priority patent/BRPI0708620A2/pt
Publication of WO2007105648A1 publication Critical patent/WO2007105648A1/ja

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/1123Discriminating type of movement, e.g. walking or running
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/22Ergometry; Measuring muscular strength or the force of a muscular blow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A61B5/7207Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
    • A61B5/7214Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts using signal cancellation, e.g. based on input of two identical physiological sensors spaced apart, or based on two signals derived from the same sensor, for different optical wavelengths
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C22/00Measuring distance traversed on the ground by vehicles, persons, animals or other moving solid bodies, e.g. using odometers, using pedometers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C22/00Measuring distance traversed on the ground by vehicles, persons, animals or other moving solid bodies, e.g. using odometers, using pedometers
    • G01C22/006Pedometers
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06MCOUNTING MECHANISMS; COUNTING OF OBJECTS NOT OTHERWISE PROVIDED FOR
    • G06M3/00Counters with additional facilities
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02Details of sensors specially adapted for in-vivo measurements
    • A61B2562/0219Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches

Definitions

  • Body motion detection device body motion detection method, and body motion detection program
  • the present invention relates to an apparatus, a method, and a program for detecting and utilizing the body movement of a user (user) such as a pedometer.
  • acceleration is detected by a plurality of body motion sensors having different detection directions, and an action axis is determined / selected by analyzing a signal pattern of each sensor.
  • a method for detecting gait by axis signal analysis is disclosed.
  • Japanese Patent Laid-Open No. 2005-140533 has a biaxial or triaxial acceleration sensor orthogonal to each other, and estimates the motion direction of the combined vector force body motion of the sensor signal of each axis, and the estimated motion direction
  • a method for detecting body movement by analyzing the signal component of the above is disclosed.
  • body movements other than walking means movements other than the vertical direction that occur in the user's body in response to walking, and a pedometer for measuring the number of steps of the user is used by the user.
  • Various movements such as the movement when the hand is swung horizontally and the pendulum movement that occurs when the pedometer for measuring the number of steps of the user is hung on the user's neck by the hanging strap, etc. Movement can be considered.
  • the user's walking pitch walking speed or the number of steps per unit time (number of times)
  • the user's movement state can be accurately grasped, and the number of steps can be further measured. It is useful to be able to do it accurately.
  • the present invention makes it possible to accurately detect a user's vertical body motion (motion) without being affected by various noises, and to make a user's use as necessary.
  • the purpose is to make it possible to accurately detect and use one body motion (motion) pitch.
  • the body motion detection device according to the invention of claim 1 provides:
  • An acceleration sensor adapted to be worn on the user's body
  • Vertical component extraction means for extracting a vertical component of acceleration from the detection output of the acceleration sensor force, and the like;
  • Decomposition means for decomposing the vertical component extracted by the vertical component extraction means into a high frequency component and a low frequency component
  • Detecting means for detecting a peak position candidate based on the low frequency component of the vertical component decomposed by the decomposing means; If the ratio of the energy of the low-frequency component and the energy of the high-frequency component in a predetermined range including the peak position candidate detected by the detecting means is smaller than a predetermined value, the peak position candidate is determined as a peak position.
  • the vertical component of the user's body motion is extracted from the detected output of the acceleration sensor force by the vertical component extraction means.
  • the extracted vertical component is decomposed into a high-frequency component and a low-frequency component by a decomposing means.
  • the low-frequency component obtained by the decomposition includes a component corresponding to the user's body movement, and is used by the detection means to detect a candidate body movement peak position.
  • the ratio of the energy of the high frequency component to the energy of the low frequency component is obtained.
  • High frequency components are likely to be mixed with noise.
  • the force for identifying the peak position candidate in the predetermined range as the peak position conversely,
  • the peak position candidates in the predetermined range are excluded from the peak position candidates. Based on the peak position candidates identified in this way, the body motion of the user is detected by the body motion detection means.
  • the body motion detection device according to claim 2 of the present application is the body motion detection device according to claim 1,
  • Determining means for determining a comparison target range for a predetermined range including each peak position candidate specified by the specifying means, and determining the peak position candidate as a peak position if they match; Prepared,
  • the body movement detecting means detects a user's body movement based on the peak position determined by the determining means.
  • a paired section is determined and the waveforms are compared, If the waveforms in both sections match, the candidate for the peak position of the comparison source is determined as the peak position. If they do not match, they are excluded from the peak position.
  • the body motion detection device according to the invention of claim 3 is the body motion detection device according to claim 1 or claim 2,
  • the step interval is estimated by analyzing the time series pattern that is the candidate power of the plurality of peak positions specified by the specifying means or the time series pattern that is the plurality of peak position forces determined by the determining means. Interval estimation means; and
  • a discriminating means for discriminating at least three states of “stationary”, “walking Z running”, and “undefined” based on the step interval estimated by the interval estimating means as a user's operation state;
  • the step estimation interval of the peak position candidates specified by the specifying means (the interval between the specified peak position candidates) by the interval estimation means. )
  • the step interval of the peak position determined by the determining means (interval between the determined peak positions) is estimated, and based on the estimated step interval, the user's operating state is set to “static It is determined by the determination means whether it is in the state of “walking Z running” or “indefinite”.
  • the body motion detection device of the invention of claim 11 is the body motion detection device of claim 1
  • the acceleration sensor is multi-axis
  • the vertical component extraction means calculates a gravitational acceleration vector from an acceleration vector that is a detection output from the multi-axis acceleration sensor, and calculates the acceleration vector from the multi-axis acceleration sensor and the calculated gravitational acceleration vector. It is characterized in that the vertical component of acceleration is extracted by calculation using.
  • the vertical component extraction means calculates the gravitational acceleration vector from the detection output (acceleration vector) from the multi-axis acceleration sensor, and multi-axis The vertical component of the acceleration is extracted by calculation using the detection output from the acceleration sensor and the calculated gravitational acceleration vector.
  • the vertical component of the detection output (acceleration vector) can be logically and accurately extracted regardless of the manner in which the multi-axis acceleration sensor is mounted on the user's body. To be able to.
  • FIG. 1 is a block diagram for explaining a basic configuration of a body motion detection device according to an embodiment.
  • FIG. 2 is a diagram for explaining an acceleration vector an, a gravitational acceleration vector g, and a vertical component vn of the acceleration vector an when a three-axis acceleration sensor is used.
  • FIG. 3 is a diagram for explaining a case where the vertical component vn of the acceleration vector an is obtained in consideration of the deviation of the gravitational acceleration vector g when a triaxial acceleration sensor is used.
  • FIG. 4 is a diagram for explaining a case where a horizontal component hn of an acceleration vector an is obtained when a triaxial acceleration sensor is used.
  • FIG. 5 is a diagram for explaining an acceleration vector an, a gravitational acceleration vector g, and a vertical component vn of the acceleration vector an when a biaxial acceleration sensor is used.
  • FIG. 6 is a diagram for explaining an example of a graph of acceleration data (A), acceleration vector length (B), vertical component (C), and water weight (D).
  • FIG. 7 is a diagram showing a graph of an acceleration signal detected when the body motion detection device of the embodiment is worn on a user's waist and used.
  • FIG. 8 is a diagram showing a graph of acceleration signals detected when the body movement detection device of the embodiment is used in a pocket of a user's pants.
  • FIG. 9 is a diagram for explaining a calculation method for obtaining the energy ratio d.
  • Fig. 10 shows the vertical component X of the acceleration vector detected by the multi-axis acceleration sensor 1.
  • FIG. 11 is a diagram for explaining the transition of the operation state of the user.
  • FIG. 12 is a diagram showing waveforms when step detection (determination of the peak position) is performed from the detection output from the acceleration sensor 1 in the body motion detection device shown in FIG.
  • FIG. 13 is a diagram for explaining an example of a calculation formula for a reference pitch Ps.
  • FIG. 14 is a diagram for explaining an example when a peak position detection failure occurs in the “walking / running” state.
  • FIG. 15 is a block diagram for explaining a pedometer to which one embodiment of the present invention is applied.
  • FIG. 16 is a flowchart for explaining vertical component extraction processing.
  • FIG. 17 is a flowchart for explaining peak position candidate detection and identification processing.
  • FIG. 18 is a flowchart for explaining processing for determining the peak position of the intermediate force specified as the candidate for the peak position and counting the number of steps of the user based on this.
  • FIG. 19 is a flowchart for explaining a process for mainly preventing a count omission near the head.
  • FIG. 20 is a flowchart for explaining body movement type estimation processing and step count counting processing.
  • FIG. 21 is a flowchart following FIG.
  • FIG. 22 is a block diagram for explaining a sound reproducing device 200 to which one embodiment of the present invention is applied.
  • FIG. 23 is a flowchart for explaining processing when an appropriate playlist is selected in accordance with the operation state of the user.
  • FIG. 1 is a block diagram for explaining the basic configuration of the body motion detection device of this embodiment.
  • the body motion detection device of this embodiment includes an acceleration sensor 1, a vertical component extraction unit 2, a high frequency Z low frequency decomposition unit 3, a peak detection Z determination processing unit 4, It consists of a step position analysis unit 5.
  • the body motion detection device of this embodiment is used by being worn on the user's body.
  • the acceleration sensor 1 detects an acceleration corresponding to the user's body movement at every fixed timing, and supplies this detection output to the vertical component extraction unit 2.
  • the vertical component extraction unit 2 extracts, from the detection output from the acceleration sensor 1, the vertical component of the acceleration vector including the component corresponding to the body movement in the vertical direction when the user performs a motion such as walking or running. This is supplied to the high-frequency Z low-frequency decomposition unit 3.
  • the analog Z-digital (AZD) conversion of the acceleration vector detected by the acceleration sensor 1 may be performed on the acceleration sensor 1 side, the vertical component extraction unit 2 side, or the acceleration sensor. It may be performed between 1 and the vertical component extraction unit 2.
  • the high-frequency Z low-frequency decomposition unit 3 determines the vertical component of the acceleration vector from the vertical component extraction unit 2 according to the high-frequency component in which noise is likely to be mixed and the vertical movement of the user. It is decomposed into low-frequency components including the components, and each of them is supplied to the peak detection Z judgment processing unit 4.
  • the peak detection Z judgment processing unit 4 receives the high-frequency component and low-frequency component of the vertical component of the acceleration vector of the high-frequency Z low-frequency decomposition unit force, and detects the body motion based on the detected peak and the detected peak. It is a part to perform.
  • the peak detection Z determination processing unit 4 detects a peak position candidate based on the low frequency component of the vertical components of the acceleration vector supplied from the high frequency Z low frequency decomposition unit 3,
  • the peak position is specified as the candidate for the peak position of the vertical component.
  • the peak detection Z determination processing unit 4 of the body motion detection device of this embodiment also performs waveform matching processing. I am doing so.
  • the peak detection Z determination processing unit 4 sets a predetermined range including a peak position candidate for each identified peak position candidate, and sets the waveform for each set predetermined range as a target. Are matched with a predetermined range of waveforms, and if a match is found, a peak position candidate included in the waveform to be compared is determined as a peak position.
  • the matching is performed in order to reliably detect periodic movements such as walking and running of the user. If the waveforms of predetermined sections including candidates for different peak positions are similar, peaking is performed. This is because it is possible to determine that the candidate for the position of the mark occurs periodically.
  • the peak position determined in this way is surely generated in accordance with the user's vertical body movement, the user's vertical direction is determined according to the determined peak position.
  • the body movement can be detected with high accuracy.
  • a step position analysis unit 5 is provided.
  • the step position analysis unit 5 receives information indicating the peak position determined by the peak detection Z determination processing unit 4 and analyzes the information, and analyzes the body movement pitch (walking) And running tempo). Then, as will be described in detail later, it is possible to perform different control depending on the detected body movement pitch of the user! /
  • pitch means the speed and number of times when the same thing is repeated or when things are performed at regular intervals. Therefore, the term “body motion pitch” means the speed and number of times of body motion, and when the body motion is walking, it means the speed of walking (walking speed) or the number of steps per unit time.
  • tempo is used as a term synonymous with “pitch”. May be used.
  • the term “tempo” originally means the speed (music progression speed) specified on the score when the music is played. Therefore, the “playback tempo” of a song is the speed at which music data is played, and means the number of beats per minute (BPM).
  • the "body motion tempo (motion tempo)" of the user is the speed of the body motion (motion).
  • the user's body motion (motion) is walking (walking) or running.
  • (Running) is the number of steps per minute, and if the movement is a jump, it is the number of jumps per minute. Means the minimum unit of movement (one movement (body movement)).
  • the body motion detection device of the present invention is based on the detection output from the acceleration sensor 1 that is adapted to be worn on the user's body !, and accurately detects the body motion in the vertical direction according to the user's walking and running. It will be able to be detected.
  • the speed sensor 1 a single-axis (single-axis) sensor or a multi-axis sensor such as two or three axes can be used.
  • a uniaxial acceleration sensor is used as the acceleration sensor 1, in order to be able to detect a user's vertical body movement, there are some restrictions on the mounting position and mounting direction.
  • a single-axis acceleration sensor when attached to an arm or leg, it may be affected by the swing of the arm or leg, and there is a possibility that it may be affected other than in the vertical direction depending on the wearing position. .
  • the acceleration sensor 1 when a single-axis acceleration sensor is used as the acceleration sensor 1, for example, it is used in order to detect the body movement in the vertical direction according to the user's walking and running as accurately as possible. On the waist of the person so that the direction of acceleration detection is vertical. There is a restriction that an acceleration sensor must be mounted. However, if this restriction is observed, the detection output from the uniaxial acceleration sensor can be used as the vertical component of acceleration. In this case, the uniaxial acceleration sensor itself has a function as the vertical component extraction unit 2.
  • the mounting position and mounting direction of the acceleration sensor 1 can be made flexible. It becomes possible.
  • a multi-axis acceleration sensor is used as the acceleration sensor 1, it is necessary to extract a vertical component from the multi-axis detection output.
  • a triaxial acceleration sensor is used as the acceleration sensor 1. For this reason, a vertical component extraction unit 2 is provided after the acceleration sensor 1.
  • the vertical component extraction unit 2 (1) uses the detection output from the multi-axis acceleration sensor 1 to estimate the gravitational acceleration vector in the gravitational field, and (2) based on the gravitational acceleration vector estimation result. Then, a process of extracting a signal component in the vertical direction from the detection output of the same acceleration sensor 1 is performed.
  • the acceleration sensor 1 for the user who does not need to estimate the action axis by extracting the vertical component by comprehensively using the detection outputs for all the axes from the multi-axis acceleration sensor 1 It is possible to accurately detect the user's vertical movement without being affected by the mounting position and mounting direction.
  • Acceleration sensor 1 has three axes, X-axis, Y-axis, and Z-axis, and acceleration vector an at a certain point n obtained from acceleration sensor 1 is expressed as axn (X-axis, as shown in equation (1 1) in Fig. 2. Component), ayn (Y-axis component), azn (Z-axis component). Acceleration vector (acceleration vector data series) an force shown in Fig. 2 (1-1) An force Gravitational acceleration vector g is estimated and body motion is also detected.
  • the gravitational acceleration vector g is estimated using the acceleration vector an, as shown in equation (1-2) in Fig. 2, gx (X-axis component), gy (Y-axis component), gz (Z Axis component).
  • the vertical component vn of the acceleration vector an can be obtained by the calculation shown in Equation (1 3) in FIG. That is, the vertical component vn of the acceleration vector an is the absolute value (magnitude) of the gravitational acceleration vector g, as shown in Equation (1 3) in Fig. 2, as the product of the inner product of the gravitational acceleration vector g and the acceleration vector an. ).
  • the vertical component vn can be accurately obtained by calculation from the acceleration vector an detected by the triaxial acceleration sensor 1 and the gravitational acceleration vector g obtained from the acceleration vector.
  • the body motion in the vertical direction of the user can be detected accurately by using the detection output from the 3-axis acceleration sensor 1 in a comprehensive manner and separating only the vertical component by numerical calculation. This invention has been made.
  • the acceleration vector an obtained from the detection output from the three-axis acceleration sensor 1 can be decomposed into a lead direct component and a horizontal component
  • the horizontal direction can be analyzed by analyzing not only the vertical component but also the horizontal component. It is also possible to detect body movements.
  • the length hn of the horizontal vector can be obtained from (3-1) in Fig. 4 and (3-2) in Fig. 4.
  • the force described using the case of using the triaxial acceleration sensor 1 as an example is not limited to this.
  • the basic concept of the present invention can be applied to the case of using a biaxial acceleration sensor in the same manner as the case of using a triaxial acceleration sensor.
  • the acceleration vector an and the gravitational acceleration vector g detected by the two-axis acceleration sensor are expressed as (4 1) (acceleration vector) in Fig. 5 and (4 2) in Fig. 5 (gravity acceleration). If it is expressed as a vector), the vertical component can be calculated according to the equation (1-3) in Fig. 2, as in the case of using a 3-axis acceleration sensor.
  • the acceleration vector an is expressed in accordance with equation (4 4) in FIG.
  • the vertical component ⁇ xn can be decomposed into the horizontal component ⁇ yn orthogonal to it, and ⁇ xn and the vertical component vn match
  • FIG. 6 shows a case in which acceleration data is acquired for 4 seconds at a sampling frequency of 50 Hz when the user is exercising, such as walking, using a three-axis acceleration sensor worn by the user.
  • the obtained acceleration data and the data obtained when the components are decomposed according to the basic concept of the present invention described above are shown in a graph.
  • the horizontal axis is time (milliseconds)
  • the vertical axis is gravitational acceleration (G).
  • G gravitational acceleration
  • FIGS. 7, 8, 10, 12, and 14 described later the horizontal axis represents time (milliseconds) and the vertical axis represents gravity acceleration (G).
  • FIG. 6A is a graph of acceleration data from the triaxial acceleration sensor
  • FIG. 6B is a graph of the length (magnitude) of the calculated acceleration vector
  • 6C is a vertical component graph obtained by calculation from the triaxial acceleration data by the method described with reference to FIGS.
  • FIG. 6D is a horizontal component graph obtained by calculating the triaxial acceleration data force by the method described with reference to FIGS.
  • the graph shown in FIG. 6 shows that the user can
  • the figure shows the case where the horizontal movement occurs mainly at the 80th, 100th, and 170th samples of the force that was mainly moving in the vertical direction, and it exists as a noise component.
  • acceleration data is divided into a vertical component (Fig. 6C) and a horizontal component (Fig. 6D) according to the basic concept of the present invention. If the components are decomposed, the noise of the horizontal component can be removed from the vertical component, so that the user's vertical body movement can be detected accurately. Of course, since the noise of the vertical component can be removed from the horizontal component, the user's horizontal body movement can be accurately detected. In other words, noise reduction can be achieved by decomposing components.
  • the force is not limited to the force that extracts the vertical component of the acceleration vector by calculation.
  • the length (magnitude) of the acceleration vector of each axis is obtained by calculation, and the longest acceleration vector is used as the vertical component of the acceleration vector that most reflects the user's body movement.
  • the mounting position of the multi-axis acceleration sensor 1 with respect to the user In some cases, the mounting direction is limited to some extent. However, when the vertical component of the acceleration vector is extracted by computation, there is no restriction on the mounting position and mounting direction for the user of the multi-axis acceleration sensor 1 as described above. In addition, the degree of freedom of wearing of the body motion detection device for the user can be improved.
  • the vertical component vn extracted in the lead direct component extraction unit 2 is represented by a function x (n).
  • the vertical component x (n) of the acceleration vector from the lead direct component extraction unit 2 is supplied to the high-frequency Z low-frequency decomposition unit 3.
  • the high-frequency Z low-frequency decomposition unit 3 is configured, for example, as an LPF (Low Pass Filter), and divides the vertical component x (n) of the acceleration vector into bands, so that the high-frequency component xh (n) and the low-frequency component xl ( n) and minutes Release.
  • LPF Low Pass Filter
  • the vertical component x (n) is separated into the high-frequency component xh (n) and the low-frequency component xl (n) as described above.
  • the high-frequency component xh (n) and the low-frequency component x Kn) decomposed in this way in the high-frequency Z low-frequency decomposition unit 3 are supplied to the peak detection Z determination processing unit 4.
  • the peak detection Z determination processing unit 4 detects a peak position candidate based on the low frequency component xl (n) of the vertical component x (n) of the acceleration vector, and also detects the peak position candidate. Based on the ratio of the energy of the low-frequency component and the high-frequency component for each predetermined range using as a reference, the candidate for the peak position is specified.
  • FIG. 7 and FIG. 8 are graphs showing acceleration signals in each band detected when the body motion detection device of this embodiment is worn on different parts of the user.
  • FIG. 7 is a diagram showing a graph of acceleration signals detected when the body motion detection device of this embodiment is used while being fixed to the user's waist
  • FIG. It is a figure which shows the graph of the acceleration signal detected when the body movement detection apparatus of this embodiment is put in a user's pants pocket and used.
  • FIG. 7 and FIG. 8 are measured by sampling the acceleration signal detected by the three-axis acceleration sensor 1 during walking of the user at 50 Hz.
  • the graphs in the top row (Figs. 7A and 8A) for both Fig. 7 and Fig. 8 are the vertical components extracted from the three-axis acceleration vector force, that is, x (n), as shown in the middle row (Figs. 7B and 8B).
  • the graph is the low-frequency component xl (n) in the vertical component x (n), and the graph in the bottom row (Figs. 7C and 8C) is the high-frequency component xh (n) of the vertical component x (n).
  • the body motion detection device of this embodiment when the body motion detection device of this embodiment is worn on the user's waist and used, only the component corresponding to the user's vertical body motion is efficiently measured. Therefore, body movements such as walking and running can be detected by threshold judgment or peak detection of the waveform of the low frequency component xl (n) of the vertical component X (n). Process, a pedometer that can accurately count the number of steps can be realized.
  • the detection output of the triaxial acceleration sensor 1 when the body motion detection device of this embodiment is used in a user's pocket is used.
  • the vibration component (noise component) other than the periodic vertical motion according to walking is the low-frequency component xl (n) ⁇ ko shown in the middle, and the high-frequency component xh (n) shown in the bottom
  • the component corresponding to the vertical movement of the user's body by walking which is a periodic motion, appears prominently in the low-frequency component xl (n), and the noise component is the high-frequency component xh (n) Appears prominently. Therefore, for each peak position candidate detected as a peak position candidate, a predetermined region having a predetermined time width is determined before and after the peak position candidate.
  • the front M sample interval and the subsequent M sample interval (M is 1 or more)
  • a predetermined range may be determined such as (integer).
  • the minimum value is detected as a peak value because it is relatively insensitive to noise, and the corresponding peak position is specified as a peak position candidate! / RU
  • the energy ratio d between the energy of the low-frequency component xl (n) and the energy of the high-frequency component xh (n) is determined for each predetermined region including the peak position candidates thus determined.
  • FIG. 9 is a diagram for explaining a calculation method for obtaining the energy ratio d.
  • the energy eh of the high-frequency component xh (n) of the vertical component X (n) of the acceleration beta can be obtained from (5-1) in Fig. 9.
  • the energy el of the low-frequency component xl (n) of the vertical component x (n) of the acceleration vector can be obtained by (5-2) in FIG.
  • nl represents the start position of the predetermined region including the peak position candidate
  • n2 represents the predetermined region including the peak position candidate. Indicate the end position.
  • the energy ratio d between the energy of the low-frequency component xl (n) and the energy of the high-frequency component xh (n) can be obtained by the equation (5-3) in FIG.
  • a position detected as a peak position candidate included in the predetermined region is specified as a peak position candidate.
  • the position specified as the peak position candidate was determined to have reached a peak according to the user's movement such as walking or running, and walking and running were performed by counting the peak position candidates. The number of steps of the user can be accurately counted.
  • the influence of vibration components other than walking and running can be eliminated.
  • the low frequency component xl (n) of the vertical component x (n) of the acceleration vector is detected.
  • the peak position candidates are detected based on the energy ratio d between the energy el of the low frequency component xl (n) and the energy eh of the high frequency component xh (n) based on Even if this is done, it may still be misjudged.
  • FIG. 10 shows a low-frequency component xl (vertical component x (n) of the acceleration vector detected by the multi-axis acceleration sensor 1 due to the influence other than the periodic movement of the user such as walking and running.
  • n It is a figure which shows the graph (irregular waveform) when an irregular signal mixes. For example, if an irregular motion other than walking or running occurs, such as falling, dropping the body motion detection device, or jumping by the user, the low-frequency component xl (n) As shown in the graph, irregular body movement components other than walking and running may be mixed.
  • irregular body motion components other than walking and running are the vertical component X of the acceleration vector.
  • the vertical component x (n) (Fig. 10A), low-frequency component xl (n) (Fig. 10B), high-frequency component xh ( n)
  • positions A, B, and C can be detected as peak positions.
  • pairs are made such as position A and position B, position B and position C, position A and position C, and matching is performed for the waveforms around each peak position between the pairs.
  • the waveform in the vicinity of position A, the waveform in the vicinity of position B, and the waveform in the vicinity of position C are not similar to each other. It is possible to judge that the movement is different from the movement of walking or running, which is a normal movement, and not the movement of walking or running.
  • the vertical component x (n), the low frequency component xl (n), and the high frequency component xh (n) of the acceleration vector are detected. If a predetermined section including the position specified as a peak position candidate is set and matching processing is performed with the waveform of one or more predetermined sections before or after it, A position identified as a candidate peak position can be determined as a formal peak position.
  • the minimum value is taken as the peak value in the graph of the low frequency component X 1 (n) so that the influence of noise is small.
  • the similarity between waveforms including adjacent peak positions is very high, so that a position specified as a candidate for each peak position can be determined as a peak position.
  • the waveform including every other peak position candidate of the peak position is taken. A position identified as a candidate for each peak position where the degree of similarity is very high can be determined as the peak position.
  • the peak position candidates are specified, and other predetermined sections to be compared (matched) are determined based on the waveforms of the predetermined sections including the candidates, and the waveforms of the two predetermined sections are compared.
  • the peak position candidate included in the predetermined section to be compared can be determined as the peak position.
  • predetermined sections to be compared can be determined as appropriate, for example, between adjacent predetermined sections, every other predetermined section, or every other predetermined section. . Or, compare with the waveform of one or more predetermined sections in the previous (past direction), and if the similarity between two or more predetermined sections is high, the peak position candidate of the comparison source section is The position and number of the predetermined sections to be compared can be arbitrarily set, such as determining the peak position.
  • the peak detection wrinkle determination processing unit 4 simply has the information (waveform) of the low frequency component xl (n) of the vertical component X ( ⁇ ) of the acceleration vector.
  • the peak position candidates are detected from the above, and the energy ratio d between the low band component energy el and the high band component energy eh is taken into account to identify the peak position candidates.
  • waveform comparison waveform matching
  • the functions of the acceleration sensor 1, the vertical component extraction unit 2, the high-frequency Z low-frequency decomposition unit 3, and the peak detection Z determination processing unit 4 described above enable the acceleration sensor 1 to be worn on the user's body. Based on this detection output, it is possible to accurately detect vertical body movements that occur in the user's body as the user walks or runs. By counting the body movements thus detected, it is possible to accurately count the number of steps of the user.
  • the user's operating state may be a stationary state that is not just a walking or running state, or an indeterminate state that is neither stationary nor operating. Can be. For this reason, if the user's movement state can be accurately grasped, the body movement is counted even when the body movement is counted, such as when the number of steps is counted, or when the movement is stationary or indefinite. It is possible to further prevent malfunctions.
  • the step position analysis unit 5 is supplied with information indicating the peak position determined by the peak detection Z determination processing unit 4, and receives it. Based on this, it is possible to accurately grasp the movement state of the user and accurately detect the body movement pitch when the user performs periodic exercise such as walking and running. So that you can!
  • FIG. 11 is a diagram for explaining the transition of the operation state of the user. As shown in Fig. 11, it is assumed that there are three user movement states: “undefined”, “still”, and “walking Z-run”. Here, “stationary” means that there is no movement of the user's body, “walking / running” means that the user is walking or running, and “undefined” means a state other than the above two states. The The initial state starts from “undefined”. In addition, when transitioning from the “Standing” state to the “Walking Z Running” state, or when transitioning from the “Walking Z Running” state to the “Standing” state, the “Undefined” state shall be used.
  • the operation state of the user is determined.
  • the user's operating state is determined by determining the transition states of transition A to transition G. [0089] That is, transition A is performed when the “undefined” state is maintained, transition B is performed when the “undefined” state is changed to the “walking Z traveling” state, and transition is performed when the “walking Z traveling” state is maintained.
  • Each transition state is determined based on the peak position determined by the peak detection Z determination processing unit 4 as described above.
  • “determined peak position” is also referred to as “detected step”.
  • the peak detection force determined by the peak detection Z determination processing unit 4 is regarded as a step corresponding to the user's vertical body movement.
  • the time interval between adjacent steps is calculated for several new steps among the detected steps (determined peak positions). If the time interval is constant, it is the state of transition B that transitions from “undefined” to “walking / running”, and it is determined that the user's operating state is “walking Z running”. .
  • the reference step interval (reference pitch Ps) is calculated in advance. Note that the reference step interval is an average value of intervals between steps used to determine that it is transition B, as will be described later.
  • an error range in which the interval between the newly detected step and the immediately preceding step is an integer multiple of the reference step interval (reference pitch Ps). If it is within the range, it is determined that the state is transition C, and the “walking / running” state is maintained. In addition, after determining that it is transition B or transition C, if the condition for determining that it is transition C in (3) above is not met, it is the state of transition D and the user's action The state is determined to be an “undefined” state.
  • whether or not the step interval is constant may be determined based on whether there is a variance or standard deviation of the step interval, Judgment is made based on whether the difference between the maximum and minimum values is below the threshold.
  • the reference step interval may be the average of the time intervals of the steps used for the determination as described above, or may be used for the determination, or the intermediate value of the time intervals of the steps may be used.
  • FIG. 12 shows the functions of the acceleration sensor 1, the vertical component extraction unit 2, the high-frequency Z low-frequency decomposition unit 3, and the peak detection Z determination processing unit 4 of the body motion detection apparatus shown in FIG.
  • FIG. 5 is a diagram showing waveforms when step detection (determination of peak position) is performed from the detection output from the acceleration sensor 1; Also in FIG. 12, the peak is taken on the minimum value side.
  • the waveform shown in Fig. 12A shows the waveform when four steps (determined peak positions) are detected in the "indeterminate" state force.
  • Fig. 12B shows the next step detected from the state shown in Fig. 12A. The waveform when this is done is shown.
  • An example of state transition will be specifically described with reference to FIG.
  • the step position analysis unit 5 calculates each step interval of adjacent step sections Tl, ⁇ 2, ⁇ 3 and compares them. .
  • the intervals of each of the step sections Tl, ⁇ 2, and ⁇ 3 are almost constant.
  • the state is a transition ⁇ transition from ⁇ indefinite '' to ⁇ walking ⁇ running ''. be able to.
  • the step position analyzing unit 5 calculates the reference pitch (reference step interval) Ps by the equation (6-1) shown in FIG.
  • the equation (6-1) shown in FIG. 13 is an equation for calculating the reference pitch Ps in the state of FIG. 12A. The section to be added and the number of sections in the denominator change.
  • step interval during walking or running in the “walking Z running” state may use the reference pitch as a representative value, or the average value of step intervals within a certain time interval, You can use the average value of a certain number of step intervals!
  • the “walking” state and the “running” state are not distinguished from each other, and the “walking Z running” state is understood as one state. However, it is not limited to this. Try to understand the “walking” and “running” states as different states.
  • the "walking" state and the “running” state may be defined as completely different motion states, but the state transition is based on FIG.
  • A above-mentioned conditions
  • B standards
  • the step interval and the reference pitch (reference step interval) Ps are not directly compared, and an integer multiple of the reference pitch Ps is used as a reference. Judgment was made. By doing so, it is useful to prevent unintended state transitions when there is a step detection failure.
  • the output Z determination processing unit 4 detects the peak position candidates, identifies the peak position candidates in consideration of the energy in the high and low bands, and further, at every predetermined interval including the peak position candidates.
  • the peak position is determined by comparing the waveforms.
  • Fig. 14 shows that in the "walking / running" state, the two steps (peak positions) indicated by the X mark (the cross mark) have been missed (count missing) and are detected most recently.
  • the interval between the step (peak position) and the step (peak position) detected immediately before is 3 steps.
  • the distance between the newly detected step and the immediately preceding step is determined based on the condition. If it is within a certain error range with respect to an integer multiple of the reference pitch Ps), it is determined that the state is transition C, and the “walking / running” state is maintained.
  • the body motion detection device of this embodiment includes the acceleration sensor 1, the vertical component extraction unit 2, the high-frequency Z low-frequency decomposition unit 3, the peak detection Z determination processing unit 4, and the step position analysis unit 5 As these parts function organically, it is possible to accurately detect and grasp body movements such as walking and running of the user.
  • the present invention is divided into a case where it is applied to a pedometer and a case where it is applied to a sound reproducing device.
  • FIG. 15 is a block diagram for explaining the pedometer 100 of this embodiment.
  • the pedometer 100 includes a three-axis acceleration sensor 101 connected to a control unit 110 via an AZD conversion unit 102, a display unit 103, and an operation unit 104. Are connected and formed.
  • the three-axis acceleration sensor 101 is formed by arranging a single-axis acceleration sensor on three axes of the X axis, the Y axis, and the Z axis orthogonal to each other, or an orthogonal three axis acceleration sensor. Although it is formed by encapsulating in one package, it is possible to use a gap.
  • the detection output (analog output) for each of the X, Y, and Z axes from the 3-axis acceleration sensor 101 is supplied to the AZD conversion unit 102, where it can be processed by the control unit 110. Is converted to digital data and supplied to the control unit 110.
  • the control unit 110 controls each unit of the pedometer 100 of this embodiment. As shown in FIG. 15, a CPU (Central Processing Unit) 111, a ROM (Read Only Memory) 11 2, and a RAM (Random Access Memory) 113 and nonvolatile memory 114 are connected via a CPU bus 115 to form a microcomputer.
  • a CPU Central Processing Unit
  • ROM Read Only Memory
  • RAM Random Access Memory
  • the CPU 111 is a main body of processing and control in the control unit 110, such as executing various programs, forming control signals to be supplied to each unit, and performing various calculations. is there.
  • the ROM 112 stores and holds various programs executed by the CPU 111 and data necessary for processing.
  • the RAM 113 is used mainly as a work area, for example, temporarily storing the results of the process, and the numerical acceleration data from the acceleration sensor 101 supplied through the AZD conversion unit 102 is used. It is also used as a buffer memory for storing and holding.
  • the non-volatile memory 114 is a memory that does not lose its stored data even when the power of EEPROM (Electrically Erasable and Programmable ROM) or flash memory is turned off. Data that should be retained, for example, set parameters Memorize and hold programs that have been added.
  • the display unit 103 connected to the control unit 110 includes a display control circuit, such as an LCD (Liquid Crystal Display), an organic EL (Electro Luminescence) display, and a CRT (Cathode-Ray Tube).
  • the display element is displayed, and the count value of the number of steps and various guidance information are displayed according to the control of the control unit 110.
  • the display unit 103 receives the display data from the control unit 110, forms a video signal to be supplied to the display data force display element, and supplies the display signal to the display element, thereby supplying the control unit 110 with the display signal. Display information corresponding to the display data from is displayed on the display screen of the display element.
  • the operation unit 104 includes a reset key, various function keys, and the like, and can receive an operation input from the user and supply an electric signal corresponding to the input to the control unit 110. is there. In response to this, the control unit 110 controls each unit so as to perform processing according to the user's instruction.
  • control unit 110 performs vertical component extraction unit 2, high-frequency Z low-frequency decomposition unit 3, peak detection Z determination processing in the body motion detection device shown in FIG.
  • control unit 110 performs vertical component extraction unit 2, high-frequency Z low-frequency decomposition unit 3, peak detection Z determination processing in the body motion detection device shown in FIG.
  • the function as the step measuring unit that performs step measurement is also realized.
  • control unit 110 of the pedometer 100 extracts the vertical component based on the detection output from the acceleration sensor 101, as described with reference to FIG. Perform low-frequency decomposition, detect peak position candidates, identify peak position candidates, determine peak positions by waveform matching, and measure the number of steps based on the determined peak positions.
  • the user's movement state is accurately grasped, and the body movement pitch such as walking and running is accurately grasped. If this is the case, then the user's step count can be measured more accurately based on the body movement pitch.
  • FIG. 16 is a flowchart for explaining the vertical component extraction processing executed mainly in the control unit 110 of the pedometer 100 of this embodiment shown in FIG.
  • the CPU 111 of the control unit 110 When power is turned on to the pedometer 100 in this embodiment and an instruction input for measuring the number of steps is received through the operation unit 104, the CPU 111 of the control unit 110 first executes the processing shown in FIG. To do.
  • the control unit 110 starts the process of acquiring acceleration data (acceleration vector) supplied through the AZD conversion unit 102 and buffering it in the RAMI 13 (step S101), and uses the acquired acceleration vector.
  • acceleration data acceleration vector
  • step S101 the initial value of the gravitational acceleration vector g is calculated (step S102).
  • step S103 the sample counter c is reset (zero cleared) (step S103).
  • the processing from step S101 to step S103 corresponds to so-called initial processing after the power is turned on.
  • control unit 110 determines whether or not the value of sample counter c is greater than a predetermined value N (step S 104).
  • the recalculation of the gravitational acceleration vector g is performed every N samples (N is an integer of 1 or more).
  • step S104 When it is determined in step S104 that the sample counter c is larger than the specified value N, the control unit 110 performs recalculation processing of the gravitational acceleration vector g (step S105). Thereafter, the value 0 is set to the sample counter c (step S106). That is, the process of step S106 is a reset process of the sample counter c. Then, using the gravitational acceleration vector g obtained by the latest recalculation and the acceleration vector an, the vertical component vn is obtained by calculation as described with reference to FIG. 2 or FIG. 3 (step S107).
  • step S104 If it is determined in step S104 that the value of the sample counter c is not greater than the value N, the control unit 110 determines the initial value of the gravitational acceleration beta g calculated in step S102 and the latest value. Using the acceleration vector an, as described above, the vertical component vn is obtained by calculation as described with reference to FIG. 2 or FIG. 3 (step S107).
  • control unit 110 obtains the latest acceleration data (step S108), adds 1 to sample counter c (step S109), Repeat these processes. In this way, the vertical component of the acceleration vector obtained in step S107 is used for the high-frequency Z low-frequency separation process.
  • the recalculation process of the gravitational acceleration vector g in step S105 and the calculation process of the initial value of the gravitational acceleration vector in step S102 are basically the same process, and the acceleration vector The gravitational acceleration vector g is calculated as the estimated value of the gravitational acceleration by taking the moving average of the data of each axis.
  • step S102 and step S105 The processing performed in step S102 and step S105 will be specifically described.
  • the X-axis gravitational acceleration vector gx is obtained by adding the X-axis acceleration data axn for each sample from the M + 1 sample before nl before the current sample position to obtain the total value.
  • the gravitational acceleration vector gx on the X axis is obtained.
  • step S107 the latest (1-3) equation shown in FIG. 2 or the (2-3) equation shown in FIG.
  • the vertical component vn can be obtained (extracted) based on the acceleration vector an and the gravitational acceleration vector g.
  • the control unit 110 executes the processing shown in FIG. 16 so as to realize the function as the vertical component extraction unit 2.
  • the gravity acceleration is estimated every N samples in order to reduce the amount of calculation.
  • the gravitational acceleration may be estimated for each sample.
  • the calculation processing for estimating the gravitational acceleration is not limited to taking a moving average of the data of each axis. For example, a least square method may be used.
  • the vertical component of the acceleration vector extracted by the processing shown in FIG. 16 is separated into a high frequency component and a low frequency component, and as described above, for example, the acceleration caused by walking or running is separated. If the component of the band of 2Hz to 4Hz containing the main component is extracted as the low frequency component, and the component of the band higher than 4Hz is extracted as the high frequency component, the vertical component is the high frequency component of the predetermined band. And a low-frequency component in a predetermined band.
  • FIG. 17 is a flowchart for explaining peak position candidate detection and identification processing executed mainly in the control unit 110 of the pedometer 100 of this embodiment shown in FIG.
  • the processing shown in FIG. 17 is performed by separating the low frequency component xl (n) and the high frequency component xh (which are obtained by separating the vertical component X (n) extracted by the vertical component extraction processing described with reference to FIG. n), for example, every sampling.
  • the control unit 110 performs a process of searching for (detecting) a peak position candidate based on the low frequency component xl (n) separated from the vertical component x (n) (step S1).
  • step S1 the latest low-frequency component data is compared with the immediately preceding low-frequency component data for each sampling, and the immediately preceding low-frequency component data is larger. If the immediately preceding low frequency component data is greater than or equal to a predetermined value, the immediately preceding low frequency component data is detected as a peak value candidate, and the peak value position is detected as a peak position candidate.
  • control unit 110 determines whether or not a peak position candidate has been detected (step S 202), and when determining that it has not been detected, ends the processing shown in FIG.
  • the processing shown in FIG. 17 is executed again at the sampling timing.
  • step S202 If it is determined in step S202 that the peak position candidate has been detected, the control unit 110 uses the low-frequency component xl (n) in the vicinity of the peak position candidate to determine the peak.
  • the energy el of the low frequency component xl (n) in the predetermined range including the position candidate is calculated according to the equation (5-2) in FIG. 9 (step S203).
  • the control unit 110 uses the high frequency component xh (n) in the vicinity of the peak position candidate to plot the energy eh of the high frequency component xh (n) in a predetermined range including the peak position candidate. Calculate according to equation (5-1) (Step S204).
  • the control unit 110 calculates the energy el of the low frequency component xl (n) calculated in step S203 and the high frequency calculated in step S204!
  • An energy ratio d of the component xh (n) to the energy eh is calculated (step S205), and it is determined whether or not the calculated energy ratio d is smaller than a predetermined threshold D (step S206).
  • step S206 if it is determined that the energy ratio d is not smaller than the predetermined threshold value D, that is, the high frequency component xh (n) is noisy, the peak position candidate Since there is a high possibility of false detection, the processing shown in FIG. 17 is terminated without specifying the detected peak position candidate as the peak position candidate, and at the next sampling timing, The process shown in FIG. 17 is executed again.
  • step S206 if it is determined that the energy ratio d is smaller than the predetermined threshold D, that is, the high frequency component xh (n) has little noise, it is determined as a peak position candidate. Since the detected one has high credibility as a peak position, the detected peak position candidate is specified as the peak position candidate, and the peak position candidate is defined as, for example, RAM. Record in the peak position recording buffer such as 113 (step S207), and then terminate the processing shown in FIG. 17 and execute the processing shown in FIG. 17 again at the next sampling timing. Like that.
  • the control unit 110 of the pedometer 100 performs the low-frequency component xl (obtained by dividing the vertical component X (n) extracted from the acceleration vector an n) detects the peak position candidates, and then the energy el of the low frequency component xl (n) and the energy of the high frequency component xh (n) in a predetermined range including the position detected as the peak position candidate. On the basis of eh, only a highly reliable candidate for the peak position is specified as a candidate for the peak position.
  • processing shown in FIG. 17 realizes the peak detection function (peak position candidate detection and identification function) among the functions as the peak detection Z determination processing unit realized by the control unit 110. Is.
  • FIG. 18 is a process executed mainly by the control unit 110 of the pedometer 100 of this embodiment shown in FIG. 15, and determines the intermediate peak position identified as the peak position candidate, It is a flowchart for demonstrating the process which counts a user's step count based on this.
  • the processing shown in FIG. 18 performs the true peak position by performing waveform matching between predetermined sections including the peak position candidates based on the peak position candidates specified by the processing shown in FIG.
  • the user's steps are counted accurately by determining the peak position and counting the determined peak position.
  • the peak position candidates identified by the process shown in FIG. 17 are stored in a peak position recording buffer such as RAMI 13 by the process of step S207 of the process shown in FIG. Therefore, the control unit 110 searches for (searches for) a peak position candidate whose peak position recording buffer power is not checked (step S301). As will be described later, whether or not the unticked force is determined is a force that can be determined by turning on or off the checked flag that is set for each identified peak position candidate. Candidates are candidates for unchecked peak positions. [0150] Then, in step S301, control unit 110 determines whether or not the unchecked peak position candidate pl is powerful (step S302).
  • step S302 If it is determined in step S302 that the unchecked peak position candidate pi is powerful, the peak position candidate p 1 is within the past N seconds from the peak position candidate p 1. Matching processing is performed on a waveform in a predetermined range including the waveform in a predetermined range including other peak position candidates (step S303).
  • step S303 matching with one or more waveforms in a predetermined range of force depending on the sampling frequency is performed. Further, the waveform used for matching may be performed based on the waveform of the low-frequency component xl (n). Of course, the waveform of the vertical component x (n) may be used as the waveform to be compared.
  • control unit 110 determines whether or not past peak position candidate p2 that matches unchecked peak position candidate pi has a strong power (step S304). ). That is, in the determination process in step S304, whether or not the waveform in the predetermined range including the candidate for the unchecked peak position pi and the waveform in the predetermined range including the past peak position candidate having a high degree of similarity is strong or not. Judging.
  • step S304 If it is determined in step S304 that the past peak position candidate p2 that matches the unchecked peak position candidate pi is strong, the unchecked peak position candidate pi is used as the peak position. And the value “1” is added to the step count, and the counted flag of the unchecked peak position candidate pi is turned on (step S305).
  • step S306 it is determined whether or not the counted flag of the past peak position candidate p2 found as a match this time is on. If it is determined in step S30 6 that the counted flag of the past peak position candidate p2 is turned on, the unchecked peak position candidate pi The checked flag is turned on (step S307), the processing shown in FIG. 18 is terminated, and the next execution timing is awaited.
  • step S306 If it is determined in the determination process in step S306 that the counted flag of the past peak position candidate p2 is not turned on, the past peak position candidate p2 is set. About the previous peak position candidate p2, it was determined as a peak position even though there was no previous matching peak, but it was matched with the new peak position candidate pl. In addition to adding “1”, turn on the counted flag of the past peak position candidate p2 (step S308), proceed to step S307, and check the flag for the unchecked peak position candidate pl. Is turned on (step S307), the processing shown in FIG. 18 is terminated, and the next execution timing is awaited.
  • step S302 it is determined that the unchecked peak position candidate pl is not found, and in the determination process of step S304, it matches the unchecked peak position candidate pl.
  • the process shown in FIG. 18 is terminated and the next execution timing is awaited.
  • the peak position is determined by the process shown in FIG. 18, and the determined peak position is counted.
  • the number of steps according to the user's walking and running can be counted accurately.
  • the number of steps according to the user's walking and running is accurately counted by counting the peak position determined by the processing of Fig. 18, but the target accuracy is If it can be ensured, the number of steps according to the user's walking or running may be counted by counting the candidates for the specified peak position by the processing of FIG.
  • the processing shown in FIG. 18 realizes a peak position determination function among the functions as the peak detection Z determination processing unit 4 realized by the control unit 110 and measures the number of steps as a pedometer.
  • the function is realized.
  • the control unit 110 searches for (searches for) a peak position candidate pi that is N seconds before the peak position recording buffer power and is not yet counted (step S401).
  • step S401 is a process of finding a past uncounted peak position candidate. In short, it is a process of searching for a head peak position candidate. Then, the control unit 110 determines whether or not the uncounted peak position candidate pi is powerful (step S402).
  • step S402 If it is determined in step S402 that an uncounted peak position candidate pi has been found, an uncounted peak position candidate for the range from the uncounted peak position candidate pi to the present is targeted.
  • a waveform in a predetermined range including pi and a waveform in a predetermined range including a peak position candidate existing in the range from the uncounted peak position candidate p1 to the present are matched (step S403).
  • step S303 shown in FIG. 18 is a waveform matching process directed from the present to the past direction, whereas the process of step S403 is a waveform matching process from the past to the present.
  • control unit 110 determines whether or not the peak position candidate p2 in the current direction that matches the unchecked peak position candidate pi has the power (step). S404). In other words, in the determination process of step S404, it is uncounted! /, A predetermined range of waveforms including the peak position candidate pi, and the current direction (time passage direction) as seen from the peak position candidate pi. It is determined whether or not a waveform having a predetermined range including a candidate of a peak position having a high similarity is found.
  • step S404 when it is determined that the peak position candidate p2 in the current direction that matches the unchecked peak position candidate pi is powerful, the value "1" is added to the step count. At the same time, the counted flag of the unchecked peak position candidate p 1 is turned on (step S405).
  • control unit 110 determines whether or not the counted flag for peak position candidate p2 in the current direction is on (step S406). In the determination process of step S406, if it is determined that the counted flag for the peak position candidate p2 in the current direction is not on, the peak position candidate p2 in the current direction is also determined as a peak position and is counted as a step count. In addition to adding “1”, V is turned on for the peak position candidate p2 in the current direction (step S407).
  • step S406 if it is determined that the counted flag for the peak position candidate p2 in the current direction is on, or after the process of step S407, the control unit 110 finishes the processing shown in FIG. 19 and waits for the next execution timing.
  • the waveform near the first peak is later in time (future).
  • the waveform near the peak it is possible to judge whether or not the peak force is truly, and when it is a peak, it can be counted.
  • the peak position candidates are identified and counted, or the peak true position of the peak position candidate is determined and counted. It is possible to measure the number of steps during exercise such as walking and running. However, peaks that cannot be matched may occur due to noise.
  • the function as the step position analysis unit 5 realized by the control unit 110 is used to determine the user's movement state (body movement type according to the movement state).
  • the user's movement state body movement type according to the movement state.
  • FIG. 20 and FIG. 21 are flowcharts for explaining body motion type estimation processing and step count counting processing executed mainly in the control unit 110 of the pedometer 100 of this embodiment.
  • the peak position information used in the processing shown in FIG. 20 and FIG. 21 includes the peak position candidate data specified by the processing in FIG. 17 and recorded in the peak position recording buffer. Data is used to perform peak position determination processing by waveform matching, body motion type estimation processing, and step count processing.
  • FIGS. 20 and 21 are executed.
  • the processing of FIGS. 20 and 21 is realized by the function of the peak detection Z determination processing unit 4 and the function of the step position analysis unit 5.
  • control unit 110 searches for (searches for) a peak position candidate that has not been checked for peak position recording buffer power such as RAMI 13 (step S501).
  • a peak position candidate that has not been checked for peak position recording buffer power such as RAMI 13 (step S501).
  • whether or not the unticked force is determined is a force that can be determined by turning on and off the checked flag set for each identified peak position candidate.
  • the most recently recorded peak position candidate is Candidate for unchecked peak position.
  • control unit 110 determines whether or not the unchecked peak position candidate pl is powerful (step S502). If it is determined in step S502 that the unchecked peak position candidate pi is strong, the value “0 (zero)” is set in the stationary timer C (step S503). Then, the control unit 110, within a range within the past N seconds from the peak position candidate p 1, the waveform in a predetermined range including the peak position candidate p 1 and the predetermined range including other peak position candidates. Matching with the waveform is performed (step S504).
  • step S504 the force depending on the sampling frequency is matched with one or more waveforms in a predetermined range in the past. Further, the waveform used for matching may be performed based on the waveform of the low-frequency component xl (n). Of course, the waveform of the vertical component x (n) may be used as the waveform to be compared.
  • control unit 110 determines whether or not past peak position candidate p2 that matches unchecked peak position candidate pi has a strong power (step S505). ). That is, in the determination process of step S505, whether or not the waveform in the predetermined range including the candidate for the peak position with high similarity and the waveform in the predetermined range including the past peak position candidate having a high degree of similarity is the power of the determination. Judging.
  • step S505 when it is determined that the past peak position candidate p2 that matches the unchecked peak position candidate pi is strong, the unchecked peak position The position candidate pl is determined as the peak position, the value “1” is added to the step count, and the counted flag of the unchecked peak position candidate pi is turned on (step S506).
  • step S507 it is determined whether or not the counted flag of the past peak position candidate p2 found as a match this time is turned on. If it is determined in step S507 that the counted flag of the past peak position candidate p2 is on, the process proceeds to the process shown in FIG.
  • step S507 If it is determined in step S507 that the counted flag of the past peak position candidate p2 is not turned on, the previous peak position candidate p2 matches the previous peak. Although it has not existed, it has been matched with the new peak position candidate pl, so even the past peak position candidate p2 is determined as a peak position and the value ⁇ 1 '' is added to the step count.
  • the counted flag of the past peak position candidate p2 is turned on (step S508), and the process proceeds to FIG.
  • step S502 If it is determined in step S502 shown in FIG. 20 that an unchecked peak position candidate P1 is not found, the value "1" is added to the stationary timer C (step S509). ), It is determined whether or not the stationary timer C is greater than the reference value (step S510).
  • step S510 when it is determined that the stationary timer C has become larger than the reference value, the user's operation state is set to the "stationary" state (step S511). The processing in FIG. 21 is terminated and the next execution timing is awaited. Also, if it is determined in step S510 that the stationary timer C is greater than the reference value, the processing shown in FIGS. 20 and 21 is terminated without doing anything. It will wait for the next execution timing.
  • step S507 when it is determined that the counted flag of the past peak position candidate p2 is on, and after the process of step S508, In FIG. 21, the processing shown in FIG. 21 is performed.
  • the control unit 110 first turns on the checked flag for the unchecked peak position candidate pl (step S 512).
  • step S 512 control unit 110 identifies the past peak position p 1 ′ closest to the peak position p 1 with the latest checked flag turned on by the process of step S512 (step S 513 )
  • the interval TO between the peak position pi and the nearest peak position ⁇ is calculated (step S514).
  • the three peak intervals TO, Tl, and ⁇ 2 sandwiched between the four most recent peak positions can be maintained. I have to. Of course, more peak intervals may be retained and used.
  • control unit 110 determines whether or not the latest determination result of the user's motion state is the “walking / running” state (step S515). If it is determined in step S515 that the user's motion state force S is not in the [walking ⁇ running] state, as described above, the three intervals ⁇ ⁇ ⁇ , Tl, ⁇ 2 determined according to the four most recent peak positions are as described above.
  • the average value Ta is calculated (step S516), and the average value Ta is calculated for each of the peak intervals TO, Tl, and ⁇ 2, and this value is the reference value for all peak intervals. It is determined whether it is smaller than (step S517).
  • step S517 it is determined whether or not the peak intervals TO, Tl, and ⁇ ⁇ ⁇ ⁇ 2 have substantially the same pitch, and the force is such that walking and running, which are periodic body movements, are performed. It is processing.
  • the control unit 110 uses the average value Ta as a reference.
  • it is determined (specified) that the user's movement state is the “walking Z running” state (step S 518).
  • control unit 110 shifts peak intervals TO, Tl, and ⁇ 2 (step S519) o That is, peak interval T1 is shifted to peak interval ⁇ 2, and peak interval ⁇ 0 is shifted to peak interval T1. To do.
  • the peak interval TO may be initialized.
  • step S519 After the process of step S519 and in the determination process of step S517, in the determination process of step S517, the difference between each of the peak intervals TO, Tl, ⁇ 2 and the average value Ta is determined as the reference.
  • the control unit 110 ends the processing shown in FIGS. 20 and 21 and waits for the next execution timing.
  • step S5 the user's movement state is "walking Z running"
  • the control unit 110 calculates the value r by dividing the latest peak interval TO by the reference pitch Ps (step S520).
  • the processing in step S520 is processing that considers the case where what should be detected as a peak value is not detected as a peak value due to the influence of noise or the like. This is the process of calculating how many times the peak interval TO force is equal to the reference pitch Ps.
  • step S521 an integer R closest to the value r calculated in step S520 is calculated (step S521). For example, if the value r is “0.1” or “0.2”, the integer R is “0” and the value r is “0.9” or “1.1” The integer R is “1”, and if the value r is “1.9” or “2.1”, the integer R is “2”. Then, control unit 110 determines whether or not the absolute value of the value obtained by subtracting integer R from value r is smaller than a predetermined reference value (step S522).
  • the determination processing in step S522 is processing for determining whether or not the latest peak interval TO force is a multiple of the reference pitch Ps.
  • the latest peak interval TO is set to the reference pitch Ps. It is determined that it is not an integral multiple, and it is determined (specified) that the user's operating state is in the “undefined” state (step S523). Then, the control unit 110 ends the processing shown in FIGS. 20 and 21 and waits for the next execution timing.
  • step S522 If it is determined in step S522 that the absolute value of the value r force obtained by subtracting the integer R is smaller than a predetermined reference value, the latest peak interval is determined. It is determined that TO is an integer multiple of the reference pitch Ps, and it is determined whether or not the integer R1 is greater than the value “1” (step S524).
  • the determination process in step S524 is a process for determining whether or not the peak interval T 0 is more than twice the reference pitch Ps.
  • step S524 When it is determined in step S524 that the integer R is not larger than 1, the latest peak interval TO is not an interval more than twice the reference pitch, so the control unit 110 does nothing.
  • the processing shown in FIG. 20 and FIG. 21 ends, and the next execution timing is waited.
  • step S524 when it is determined in step S524 that the integer R is greater than 1 Since the latest peak interval TO is a section that is at least twice the reference pitch, the control unit 110 adds a value obtained by subtracting 1 from the integer R to the count number (step S525). The processing shown in FIG. 21 is terminated and the next execution timing is awaited.
  • the latest peak interval TO is an interval that is an integral multiple of the reference pitch
  • the number of steps is measured according to the strong peak position that cannot be measured. It is not limited to this.
  • the latest peak interval T 0 is an interval that is an integral multiple of the reference pitch
  • the amplitude of the low frequency component of the vertical component in the interval or the amplitude of the vertical component itself is also taken into consideration and an integral multiple of the reference pitch. Try to determine whether or not the peak position exists at the correct position.
  • a reference value is set for the amplitude at an integer multiple of the reference pitch, and if it is greater than or equal to this reference value, it is recognized as a peak value and peak position, and is smaller than the reference value. Is not recognized as a peak value and a peak position. It is possible to prescribe an appropriate value for the reference value by conducting an experiment.
  • FIG. 22 is a block diagram for explaining the sound reproducing device 200 of this embodiment.
  • the sound reproduction apparatus 200 of this embodiment is a portable type that can use a relatively large-capacity recording medium for storing music data (music data). belongs to.
  • a recording medium using various storage media such as a hard disk, a magneto-optical disk such as an MD (mini disk (registered trademark)), an optical disk such as a CD or a DVD, a memory card, or a semiconductor memory is considered. It is done. However, here, in order to simplify the explanation, it is assumed that the recording medium for storing content data such as music data is a hard disk. To do.
  • a three-axis acceleration sensor 201 is connected to a control unit 210 via an AZD conversion unit 202, and relatively As a recording medium having a large storage capacity, a music database (hereinafter referred to as music DB) 203 and a playlist storage unit 204 are connected.
  • music DB music database
  • a speaker 206 is connected to the control unit 210 via an audio reproduction processing unit 205, an operation unit 221 as a user interface, and a display unit 222 are connected, and an external interface ( Hereinafter, it is referred to as an external IZF.)
  • the input / output terminal 224 is connected via the 223.
  • the three-axis acceleration sensor 201 has a single-axis acceleration sensor arranged on the three axes X, Y, and Z that are orthogonal to each other. Any of those formed by encapsulating three orthogonal acceleration sensors in a single package may be used.
  • the X axis, Y axis, and Z axis detection outputs (analog output) from the 3-axis acceleration sensor 201 are supplied to the AZD conversion unit 202, where they are converted into digital data that can be processed by the control unit 210. It is converted and supplied to the control unit 210.
  • the control unit 210 controls each unit of the sound reproducing device 200 according to this embodiment.
  • the CPU 211, the ROM 212, the RAM 213, and the nonvolatile memory 214 power CPU bus 215 is connected to form a microcomputer.
  • the CPU 211 is a main body of processing and control in the control unit 210, such as executing various programs, forming control signals to be supplied to each unit, and performing various calculations. is there.
  • the ROM 212 stores and holds various programs executed by the CPU 211 and data necessary for processing.
  • the RAM 213 is mainly used as a work area, such as temporarily storing intermediate results of processing, and numerically input acceleration data from the acceleration sensor 201 supplied through the AZD conversion unit 202. It is also used as a buffer memory for storing and holding.
  • Non-volatile memory 214 for example, even if power to EEPROM or flash memory is turned off This is a memory that does not lose its stored data, and stores and holds data that should be retained even when the power is turned off, such as set parameters and added programs.
  • the music DB 203 and the playlist storage unit 204 are both hard disks as described above.
  • the music DB stores and holds a large number of music data that can be played back by the sound playback device 200 of this embodiment, for example, data compressed according to a predetermined data compression method.
  • the playlist storage unit 204 has a stationary playlist instructing music to be played when the user is in the "still” state, and a playback when the user is in the "walking" state.
  • a walking playlist that indicates the music to be played and a running playlist that indicates the music to be played when the user is in the “running” state are stored and held.
  • the playlist stored in the playlist storage unit 204 is different from the stationary playlist, the walking playlist, and the traveling playlist.
  • the playlist type added to each playlist, etc. can be distinguished by information.
  • there are multiple playlists at rest, playlists at walking, and playlists at running, and each playlist has a song ID (music identifier) such as a file name that identifies one or more songs to be played. are registered in advance in the order of reproduction, for example.
  • the music DB 203 and the power that separately shows the playlist storage unit 204 may be formed on the same hard disk.
  • the audio reproduction processing unit 205 receives supply of the music data read from the music DB 203 by the control unit 210, and compresses and decompresses the music data to restore the original music data before data compression. At the same time, an analog audio signal in a format to be supplied to the speaker 206 is formed from the restored music data and supplied to the speaker 206. As a result, sound corresponding to the music data to be played is emitted from the speaker 206.
  • the audio reproduction processing unit 205 of this embodiment is also connected to a headphone terminal.
  • the audio reproduction processing unit 205 The audio signal processed in 205 is supplied not to the speaker 206 but via the headphone terminal to the headphones connected thereto, and passes through the headphones. So that you can listen to the playback sound.
  • the operation unit 221 includes a play key, a stop key, a fast-forward key, a fast-rewind key, and various other function keys.
  • the operation unit 221 accepts user input and controls electrical signals accordingly.
  • the unit 210 can be supplied.
  • the control unit 210 controls each unit so that processing according to a user instruction can be performed.
  • the display unit 222 includes a display control circuit, and includes a display element such as an LCD (Liquid Crystal Display organic EL (Electro Luminescence) display, CRT (Cathode-Ray Tube), etc.
  • a display element such as an LCD (Liquid Crystal Display organic EL (Electro Luminescence) display, CRT (Cathode-Ray Tube), etc.
  • Various guidance information, etc. are displayed in accordance with the control, etc.
  • the display unit 222 receives display data from the control unit 210 and receives a video signal to be supplied from the display data to the display element. By forming and supplying this to the display element, display information corresponding to the display data from the control unit 210 is displayed on the display screen of the display element.
  • the external I / F 223 and the input / output terminal 224 can be connected to an external device such as a personal computer. Then, music data and playlists are supplied from an external device through the input / output terminal 224 and the external IZF 223, and the music data is stored in the music DB 203, and the playlist can be stored in the playlist storage unit 204. I can do it.
  • the music data stored in the music DB 203 is output to the external device through the control unit 210, the external IZF 223, and the input / output terminal 224, and backed up to the recording medium of the external device.
  • a list of music data that can be stored in the music DB 203 is displayed on the display unit 222, and the target music data is selected through the operation unit 221, and the resting playlist, walking playlist, and running are selected. By inputting another time playlist, the playlist can also be created in the sound reproduction apparatus 200 of this embodiment.
  • the music data selected through the operation unit 221 is read from the music DB 203 by the control unit 210, and this is the audio reproduction process. Supplied to the unit 205. As a result, the music instructed by the user through the operation unit 221 can be reproduced and listened to.
  • the music data is read from the music DB 203 by the control unit 210 according to the instructed playlist, and is supplied to the audio reproduction processing unit 205. Music can also be played according to the playlist.
  • the control unit 210 includes a vertical component extraction unit 2, a high-frequency Z low-frequency decomposition unit 3, and a peak detection Z determination processing unit 4 in the body motion detection device shown in FIG.
  • the user's movement pitch (movement tempo) is accurately grasped, and the resting playlist, walking playlist, and running according to the user's movement state. The power of on-the-fly playlists so that the appropriate one can be selected automatically.
  • control unit 210 of the sound reproduction device 200 of this embodiment based on the detection output from the acceleration sensor 201, extracts the vertical component, increases the vertical component as described with reference to FIG. It has functions to perform processing such as region Z low frequency decomposition, detection of peak position candidates and identification of peak position candidates, and determination of peak positions by waveform matching.
  • FIG. 23 is a flowchart for explaining processing when an appropriate playlist is selected in accordance with the user's operation state in the sound reproduction device of this embodiment.
  • step S601 The processing shown in Fig. 23 is mainly performed in control unit 210 of sound reproduction device 200 of this embodiment.
  • the control unit 210 sets the control unit 110 of the pedometer 100 shown in FIG.
  • body motion type estimation processing is executed (step S601).
  • step S601 (1) as described with reference to FIG. 16, the vertical component is extracted from the acceleration vector from the triaxial acceleration sensor 201, and (2) the extracted vertical component is converted to the high frequency component. And (3) using these, as described with reference to FIG. 17, identify peak position candidates, and (4) use the information on the identified peak position candidates.
  • the processing shown in FIGS. 20 and 21 is executed to grasp the user's operation state and operation pitch.
  • control unit 210 determines whether or not the user's operating state has changed force (step S602).
  • the previous determination result and the current determination result are held as the determination result of the operation state, and when the two are different, it can be determined that the operation state has changed.
  • step S602 If it is determined that the operation state has changed! /, After the determination process in step S602, the process shown in FIG. 23 is terminated, and the next execution timing is awaited. .
  • step S603 it is determined whether or not the state after the change, that is, whether or not the current determination result is a “static” state.
  • step S603 when it is determined that the current determination result is in the "still” state, the stationary playlist stored in the playlist storage unit 204 is used, and the stationary play is performed.
  • the music (music) is played according to the list (step S604). Thereafter, the processing shown in FIG. 23 is terminated and the next execution timing is awaited.
  • step S603 when it is determined that the current determination result is not in the “still” state, it is determined whether or not the current determination result force is “walking Z running” state. Tape S605). If it is determined in the determination process of step S605 that the state is the “walking Z running” state, the reference pitch Ps obtained in the process of step S601 is smaller than the specified value in the process of FIGS. It is determined whether or not (step S606). [0230] The specified value used in the determination process of step S606 is for determining whether the user's movement state 1S is walking or running. If the reference pitch Ps is not smaller than the specified value, The reference pitch is not yet fast enough and can be determined to be a walking state. Conversely, if the reference pitch is smaller than the specified value, it can be determined that the reference pitch is sufficiently fast.
  • step S606 when it is determined in step S606 that the reference pitch PS is not smaller than the specified value, it is determined that the user's motion state is the "walking" state, and the playlist storage unit 204 And a music playlist (music) is reproduced in accordance with the walking playlist (step S607). Thereafter, the control unit 210 ends the processing shown in FIG. 23 and waits until the next execution timing.
  • step S606 determines that the reference pitch PS is smaller than the specified value, it is determined that the user's operation state is the "running" state, and the playlist The travel playlist stored in the storage unit 204 is used, and music (music) is reproduced according to the travel playlist (step S608). Thereafter, the control unit 210 ends the processing shown in FIG. 23 and waits for the next execution timing.
  • the sound reproduction device 200 appropriately determines the user's operation state, automatically switches to the playlist according to the user's operation state, and performs the user's body.
  • the music corresponding to the movement can be played.
  • each of the stationary playlist, the walking playlist, and the traveling playlist may be prepared.
  • the order of use should be determined in advance and used according to that order, or the frequency of use may be stored for each playlist, and the less frequently used may be used, On the contrary, a frequently used one may be used.
  • each of the playlist when stationary, the playlist when walking, and the playlist when traveling! The user may specify what to use in advance.
  • the processing of Fig. 19 is performed prior to the processing of Figs. 20 and 21, so that the accurate number of steps can be measured from the beginning of the step count measurement processing. Can be.
  • the triaxial acceleration sensors 101 and 201 are used. However, it is not limited to this. A uniaxial acceleration sensor may be used, or a biaxial acceleration sensor may be used. However, when using a 1-axis or 2-axis acceleration sensor, it is necessary to place it so that the main components of body movement can be picked up as signals.
  • a triaxial acceleration sensor is used, and the triaxial acceleration vector force is a force that extracts a vertical component.
  • the signal with the strongest signal strength among the three axes may be selected as the vertical component.
  • the force that determines the pitch period of the body motion by detecting the peak position of the vertical component of the acceleration vector is not limited.
  • a switch is placed on the sole and a switch is closed by stepping on the foot while walking or running, or a strain gauge is attached to the sole. It is also possible to use a means for detecting depression by arranging the
  • the body movement pitch (walking tempo and running tempo) of the user can be accurately measured, so that the control unit 210 performs the sound reproduction processing unit.
  • the playback tempo of the music to be played back can be controlled according to the body movement tempo of the user.
  • the body motion detection device, the body motion detection method, and the body motion detection program according to the above-described embodiments are used in a fitness club or the like that can be used only with a pedometer and a portable sound reproduction device.
  • a fitness club or the like that can be used only with a pedometer and a portable sound reproduction device.
  • a training machine such as an internal running machine or suitable for a stationary sound reproduction device. Can be applied.
  • FIGS. 16 to 21 and FIG. 23 are those to which the method and the program according to the present invention are applied.
  • the method according to the present invention can be used, and the program according to the present invention can be realized by creating a program according to the flowcharts shown in FIGS.
  • the method according to the present invention can be realized by defining the processing in the control unit so as to realize the function of each unit shown in FIG. 1. Similarly, the function of each unit shown in FIG. 1 is realized. By creating a program to be executed by the control unit as described above, the program according to the present invention can be realized.
  • the present invention it is possible to remove noise and accurately detect a user's vertical body movement. That is, it is possible to realize an apparatus, method, and program for detecting body movement that is resistant to noise (not easily affected by noise).

Abstract

 3軸の加速度センサ(1)からの使用者の体動に応じた加速度ベクトルから鉛直成分抽出手段(2)により鉛直成分を抽出し、この鉛直成分を高域/低域分解部(3)により高域成分と低域成分とに分解する。これらをピーク検出/判定処理部(4)において用いて、使用者の鉛直方向の動作のピーク位置の候補を検出すると共に、高域成分と低域成分のエネルギー比に基づいてピーク位置の候補を特定し、各ピーク位置の候補を含む所定範囲についての波形マッチングを行うことにより、ピーク位置を決定して、これに基づき体動を検出すると共に、ステップ位置解析部(5)において体動ピッチを検出する。これにより、種々のノイズの影響を受けることなく、ユーザーの鉛直方向の体動を正確に検出すると共に、必要に応じてユーザーの体動ピッチをも正確に検出できるようにする。

Description

明 細 書
体動検出装置、体動検出方法および体動検出プログラム
技術分野
[0001] この発明は、例えば、歩数計などの使用者 (ユーザー)の体動を検出して利用する 装置、方法、プログラムに関する。
背景技術
[0002] 自己の健康管理などのために、歩数計を用いて自己の運動量を把握することが広 く行われている。従来の歩数計では、正確に歩数をカウントするためには、その構造 上、本体の装着位置や装着方向に制約があるものが多力つた。しかし、近年におい ては、使い勝手をより向上させるために、装着位置や装着方向を自由に設定できる 体動検出装置が提案される。その中でも、多軸の加速度センサを用いて姿勢推定を 行いつつ、同じセンサを用いて体動検出を行う方法および装置が提案されている。こ れは、姿勢の推定に角度センサ等を用いる方法と比較して、安価に実現できるという 禾 IJ点がある。
[0003] 例えば、特開 2004— 141669号公報には、互いに検出方向が異なる複数の体動 センサにより加速度を検出し、各センサの信号パターンを解析することによって作用 軸を判定/選択し、作用軸の信号解析によって歩行検出を行う方法が開示されてい る。また、特開 2005— 140533号公報には、互いに直交する 2軸もしくは 3軸の加速 度センサを搭載し、各軸のセンサ信号の合成ベクトル力 体動の運動方向を推定し、 推定した運動方向の信号成分を解析することによって体動を検出する方法が開示さ れている。
発明の開示
[0004] し力しながら、特開 2004— 141669号公報に開示された方法では、複数のセンサ 力も測定に適したセンサを作用軸として一つだけ選択するため、検出する使用者の 体動の方向(運動方向)と選定した作用軸とがー致しないケースもあり、歩行以外の 体動の影響を受けて正確な歩数計測ができなくなる可能性があると考えられる。また
、特開 2005— 140533号公報に開示された方法では、信号波形の性質は考慮され ず、信号の強度のみを問題としているために、歩行以外の体動もカウントしてしまう可 能性があると考えられる。
[0005] ここで、歩行以外の体動とは、歩行に応じて使用者の身体に発生する鉛直方向以 外の動きを意味し、使用者の歩数を計測するための歩数計が使用者の手に持たれ て水平方向に振られた場合の動きや例えば使用者の歩数を計測するための歩数計 が吊り紐によって使用者の首に掛けられて 、る場合に生じる振り子運動など、様々な 動きが考えられる。
[0006] このため、歩行以外の使用者の体動の影響を受けることなぐ歩行時の使用者の鉛 直方向の体動のみを正確に検出して、これをカウントできるようにし、使用者の歩数を 正確に計測できるようにすることが望まれる。さらに、使用者の体動を詳しく見ると、静 止している状態の場合もあれば、静止状態から歩行状態へ、あるいは、歩行状態か ら静止状態への過渡期の状態もある。また、動作状態であっても歩行 (ウォーキング) 状態と走行 (ランニング)状態なども異なる。このため、使用者の歩行ピッチ (歩行速 度、あるいは、単位時間当たりの歩数(回数))を正確に把握することができれば、使 用者の動作状態を正確に把握でき、歩数の計測をさらに正確に行うことができるなど 役立つことになる。
[0007] 以上のことにかんがみ、この発明は、種々のノイズの影響を受けることなぐユーザ 一の鉛直方向の体動(動作)を正確に検出できるようにすると共に、必要に応じてュ 一ザ一の体動(動作)ピッチをも正確に検出して利用できるようにすることを目的とす る。
[0008] 上記課題を解決するため、請求項 1に記載の発明の体動検出装置は、
使用者の身体に装着するようにされる加速度センサと、
前記加速度センサ力ゝらの検出出力から加速度の鉛直成分を抽出する鉛直成分抽 出手段と、
前記鉛直成分抽出手段によって抽出された前記鉛直成分を高域成分と低域成分 とに成分分解する分解手段と、
前記分解手段によって分解された前記鉛直成分の前記低域成分に基づいてピー ク位置の候補を検出する検出手段と、 前記検出手段によって検出されたピーク位置の候補を含む所定範囲における前記 低域成分のエネルギーと前記高域成分のエネルギーとの比が所定値よりも小さい場 合に前記ピーク位置の候補をピーク位置の候補として特定する特定手段と、 前記特定手段によって特定されるピーク位置の候補に基づ 、て、前記使用者の体 動を検出する体動検出手段と
を備えることを特徴とする。
[0009] この請求項 1に記載の発明の体動検出装置によれば、加速度センサ力 の検出出 力から鉛直成分抽出手段により使用者の体動の鉛直成分が抽出される。この抽出さ れた鉛直成分は、分解手段によって高域成分と低域成分とに分解される。分解され て得られた低域成分は使用者の体動に応じた成分を含むものであり、検出手段によ り用いられて体動のピーク位置の候補が検出される。
[0010] そして、検出された各ピーク位置の候補を含む所定の範囲毎に、低域成分のエネ ルギ一に対する高域成分のエネルギーの比を求める。高域成分にはノイズが混入し 易い。このため、低域成分のエネルギーに対する高域成分のエネルギーの比が所定 値より小さい場合には、その所定範囲のピーク位置の候補をピーク位置として特定す る力 逆に、低域成分のエネルギーに対する高域成分のエネルギーの比が所定値よ り大きい場合には、その所定範囲のピーク位置の候補は、ピーク位置の候補からは 除外する。このようにして特定されたピーク位置の候補に基づいて、体動検出手段に より使用者の体動が検出される。
[0011] これにより、ノイズなどの影響を受けることなぐ使用者の鉛直方向の体動を精度よく 検出することができるよう〖こされる。
[0012] また、この出願の請求項 2に記載の発明の体動検出装置は、請求項 1に記載の体 動検出装置であって、
前記特定手段によって特定された各ピーク位置の候補を含む所定範囲について、 比較対象の範囲を定めて波形を比較し、一致した場合にそのピーク位置の候補をピ ーク位置として決定する決定手段を備え、
前記体動検出手段は、前記決定手段によって決定されたピーク位置に基づいて、 使用者の体動を検出することを特徴とする。 [0013] この請求項 2に記載の発明の体動検出装置によれば、特定手段によって特定され た各ピーク位置の候補を含む所定区間について、ペアとなる区間を定めて波形を比 較し、両区間の波形が一致した場合には、比較元のピーク位置の候補をピーク位置 として決定し、一致しな 、場合にはピーク位置からは除外する。
[0014] これによつて、ノイズが混入することによってたまたまピークが生じた場合を除去して 、真に使用者の鉛直方向の体動に応じてピーク位置のみを抽出、特定して、使用者 の体動を正確に検出することができるようにされる。
[0015] また、請求項 3に記載の発明の体動検出装置は、請求項 1または請求項 2に記載 の体動検出装置であって、
前記特定手段によって特定された複数のピーク位置の候補力 なる時系列パター ンについて、あるいは、前記決定手段によって決定された複数のピーク位置力 なる 時系列パターンについて解析することにより、ステップ間隔を推定する間隔推定手段 と、
前記間隔推定手段において推定された前記ステップ間隔に基づいて、少なくとも「 静止」、「歩行 Z走行」、「不定」の 3状態を使用者の動作状態として判別する判別手 段と
を備えることを特徴とする。
[0016] この請求項 3に記載の発明の体動検出装置によれば、間隔推定手段によって、特 定手段によって特定されたピーク位置の候補のステップ間隔 (特定されたピーク位置 の候補間の間隔)、あるいは、決定手段によって決定されたピーク位置のステップ間 隔 (決定されたピーク位置間の間隔)が推定され、この推定されたステップ間隔に基 づいて、使用者の動作状態が、「静止」の状態にあるのか、「歩行 Z走行」の状態に あるのか、「不定」の状態にあるのかが判定手段により判定される。
[0017] これにより、使用者の動作状態を正確に把握し、使用者の動作状態に応じた機器 制御が可能になると共に、使用者力 ^歩行 Z走行」状態にある場合においては、その ステップ間隔を正確に把握することもできるようにされる。
[0018] また、請求項 11に記載の発明の体動検出装置は、請求項 1に記載の体動検出装 置であって、 前記加速度センサは多軸のものであり、
前記鉛直成分抽出手段は、前記多軸の加速度センサからの検出出力である加速 度ベクトルから重力加速度ベクトルを算出し、前記多軸の加速度センサからの前記 加速度ベクトルと算出された前記重力加速度ベクトルとを用いて演算により加速度の 鉛直成分を抽出するものであることを特徴とする。
[0019] この請求項 11に記載の発明の体動検出装置によれば、鉛直成分抽出手段は、多 軸の加速度センサからの検出出力(加速度ベクトル)から重力加速度ベクトルを算出 し、多軸の加速度センサからの検出出力と算出した重力加速度ベクトルとを用いた演 算により、加速度の鉛直成分を抽出する。
[0020] これにより、多軸の加速度センサが使用者の身体にどのような態様で装着されてい ても、その検出出力(加速度ベクトル)の鉛直成分を論理的に、かつ、正確に抽出す ることができるようにされる。
図面の簡単な説明
[0021] [図 1]図 1は、実施の形態の体動検出装置の基本的な構成について説明するための ブロック図である。
[図 2]図 2は、 3軸の加速度センサを用いた場合における、加速度ベクトル an、重力加 速度ベクトル g、加速度ベクトル anの鉛直成分 vnにつ 、て説明するための図である。
[図 3]図 3は、 3軸の加速度センサを用いた場合における、重力加速度ベクトル gの偏 角を考慮して加速度ベクトル anの鉛直成分 vnを求める場合を説明するための図で ある。
[図 4]図 4は、 3軸の加速度センサを用いた場合における、加速度ベクトル anの水平 成分 hnを求める場合を説明するための図である。
[図 5]図 5は、 2軸の加速度センサを用いた場合における、加速度ベクトル an、重力加 速度ベクトル g、加速度ベクトル anの鉛直成分 vnにつ 、て説明するための図である。
[図 6]図 6は、加速度データ (A)、加速度ベクトルの長さ(B)、鉛直成分 (C)、水平成 分 (D)のグラフの一例を説明するための図である。
[図 7]図 7は、実施の形態の体動検出装置を使用者の腰部に装着して用いた場合に 検出される加速度信号のグラフを示す図である。 [図 8]図 8は、実施の形態の体動検出装置を使用者のズボンのポケットに入れて用い た場合に検出される加速度信号のグラフを示す図である。
[図 9]図 9は、エネルギー比 dを求めるための計算方式を説明するための図である。
[図 10]図 10は、多軸の加速度センサ 1により検出される加速度ベクトルの鉛直成分 X
(n)の低域成分 xl (n)にも不規則な信号が混入した場合のグラフを示す図である。
[図 11]図 11は、使用者の動作状態の遷移を説明するための図である。
[図 12]図 12は、図 1に示した体動検出装置において加速度センサ 1からの検出出力 からステップ検出(ピーク位置の決定)をした場合の波形を示す図である。
[図 13]図 13は、基準ピッチ Psの算出式の一例を説明するための図である。
[図 14]図 14は、「歩行/走行」状態において、ピーク位置の検出漏れが生じた場合の 例を説明するための図である。
[図 15]図 15は、この発明の一実施の形態が適用された歩数計を説明するためのプロ ック図である。
[図 16]図 16は、鉛直成分抽出処理を説明するためのフローチャートである。
[図 17]図 17は、ピーク位置の候補の検出及び特定処理を説明するためのフローチヤ ートである。
[図 18]図 18は、ピーク位置の候補として特定されたものの中力もピーク位置を決定し 、これに基づいて使用者の歩数を計数する処理を説明するためのフローチャートで ある。
[図 19]図 19は、主に先頭付近のカウント漏れを防止する処理を説明するためのフロ 一チャートである。
[図 20]図 20は、体動種別の推定処理と歩数カウント処理を説明するためのフローチ ヤートである。
[図 21]図 21は、図 20に続くフローチャートである。
[図 22]図 22は、この発明の一実施の形態が適用された音響再生装置 200を説明す るためのブロック図である。
[図 23]図 23は、使用者の動作状態に応じて適切なプレイリストを選択する場合の処 理につ 、て説明するためのフローチャートである。 発明を実施するための最良の形態
[0022] 以下、図を参照しながら、この発明による装置、方法、プログラムの一実施の形態に ついて説明する。
[0023] [体動検出装置の基本構成について]
図 1は、この実施の形態の体動検出装置の基本的な構成について説明するための ブロック図である。図 1に示すように、この実施の形態の体動検出装置は、加速度セ ンサ 1と、鉛直成分抽出部 2と、高域 Z低域分解部 3と、ピーク検出 Z判定処理部 4と 、ステップ位置解析部 5とからなるものである。
[0024] この実施の形態の体動検出装置は、使用者の身体に装着されて使用されるもので ある。加速度センサ 1は、一定のタイミング毎に使用者の体動に応じた加速度を検出 し、この検出出力を鉛直成分抽出部 2に供給する。鉛直成分抽出部 2は、加速度セ ンサ 1からの検出出力から、使用者が歩行や走行といった運動を行った場合の鉛直 方向の体動に応じた成分を含む加速度べクトルの鉛直成分を抽出し、これを高域 Z 低域分解部 3に供給する。なお、加速度センサ 1において検出された加速度ベクトル のアナログ Zデジタル (AZD)変換は、加速度センサ 1側で行ってもよいし、鉛直成 分抽出部 2側で行ってもよいし、あるいは、加速度センサ 1と鉛直成分抽出部 2との間 で行うようにしてもよい。
[0025] 高域 Z低域分解部 3は、鉛直成分抽出部 2からの加速度ベクトルの鉛直成分を、ノ ィズが混入し易い高域成分と、使用者の鉛直方向の体動に応じた成分を含む低域 成分とに分解し、そのそれぞれをピーク検出 Z判定処理部 4に供給する。ピーク検出 Z判定処理部 4は、高域 Z低域分解部力 の加速度ベクトルの鉛直成分の高域成 分と低域成分の供給を受けて、ピーク検出と検出したピークに基づく体動検出とを行 う部分である。
[0026] すなわち、ピーク検出 Z判定処理部 4は、高域 Z低域分解部 3から供給された加速 度ベクトルの鉛直成分のうちの低域成分に基づいて、ピーク位置の候補を検出し、 当該ピーク位置の候補を含む所定範囲における低域成分のエネルギーに対する高 域成分のエネルギーの比が予め決められた値よりも小さい場合に、当該ピーク位置 を鉛直成分のピーク位置の候補として特定する。 [0027] このように、低域成分のエネルギーに対する高域成分のエネルギーの比を用いる のは、後述もするように、高域成分にはノイズがのりやすぐノイズの混入によって発 生したピークを除外するためである。このようにして特定されたピーク位置の候補に基 づいて使用者の鉛直方向の体動を比較的に精度よく検出することが可能である。し かし、より精度よく使用者の鉛直方向の体動を検出できるようにするために、この実施 の形態の体動検出装置のピーク検出 Z判定処理部 4は、波形マッチング処理をも行 うようにしている。
[0028] すなわち、ピーク検出 Z判定処理部 4は、特定した各ピーク位置の候補毎に、ピー ク位置の候補を含む所定範囲を設定し、この設定した所定範囲毎の波形を対象とし て他の所定範囲の波形とのマッチングを行 、、マッチングが取れた場合に比較対象 の波形に含まれるピーク位置の候補をピーク位置として決定する。このようにマツチン グを行うのは、使用者の歩行や走行といった周期的な運動を確実に検出するためで あり、異なるピーク位置の候補を含む所定区間同士の波形が類似していれば、ピー ク位置の候補が周期的に発生していると判別することが可能になるためである。
[0029] このようにして決定されたピーク位置は、使用者の鉛直方向の体動に応じて生じた ものであることは間違いないので、決定されたピーク位置に従って、使用者の鉛直方 向の体動を精度よく検出することができるようにされる。さらに、この実施の形態の体 動検出装置においては、ステップ位置解析部 5が設けられている。
[0030] このステップ位置解析部 5は、ピーク検出 Z判定処理部 4において決定されたピー ク位置を示す情報の提供を受けてこれを解析し、使用者の歩行や走行の体動ピッチ (歩行や走行のテンポ)を検出することができるようにしている。そして、詳しくは後述 もするように、検出した使用者の体動ピッチに応じて異なる制御を行うようにするなど のことができるようにして!/、る。
[0031] なお、「ピッチ」という文言は、同じことを繰り返したり、また一定の間隔でことを行つ たりする場合において、その速度や回数を意味する。したがって、「体動ピッチ」という 文言は、体動の速さや回数を意味し、体動が歩行であるときには、歩行の速さ(歩行 速度)あるいは単位時間当たりの歩数を意味する。
[0032] そして、この明細書においては、「ピッチ」と同義の文言として、「テンポ」という文言 を用いる場合もある。「テンポ」という文言は、本来、楽曲の演奏にあたって譜面に指 定された速度 (楽曲進行の速度)を意味するものである。したがって、楽曲の「再生テ ンポ」といったときには、音楽データの再生時の速さであり、 1分間あたりの柏の数 (B PM: Beat Per Minutes)を意味する。
[0033] また、使用者の「体動テンポ (動作テンポ)」と 、つたときには、体動(動作)の速度で あって、例えば、使用者の体動 (動作)が歩行 (ウォーキング)や走行 (ランニング)な どの場合には、 1分間あたりの歩数であり、また、動作がジャンプの場合には、 1分間 当たりのジャンプの回数であるというように、 1分間当たりの使用者の計数可能な最小 の動作の単位(1つの動作 (体動))の数を意味する。このように、使用者の体動(動作
)につ 、て用いる「ピッチ」と 、う文言と「テンポ」 t 、う文言とは、この明細書にお!、て はほぼ同義の文言として用いる。
[0034] [体動検出装置の各部の機能、動作について]
以下に、この実施の形態の体動検出装置を構成する各部毎の機能や動作につい て、より詳細に説明する。
[0035] [加速度センサ 1について]
まず、加速度センサ 1について説明する。この発明の体動検出装置は、使用者の 身体に装着するようにされる加速度センサ 1からの検出出力に基づ!、て、使用者の 歩行や走行に応じた鉛直方向の体動を正確に検出できるようにするものである。カロ 速度センサ 1としては、 1軸(単軸)のものや、 2軸や 3軸などの多軸のものを用いるこ とが可能である。
[0036] 加速度センサ 1として、 1軸の加速度センサを用いた場合には、使用者の鉛直方向 の体動を検出できるようにするために、装着位置や装着方向につ 、てある程度制約 を受ける。例えば、腕や足などに 1軸の加速度センサを装着した場合には、腕や足の 振りの影響を受けることが考えられるなど、装着する位置によって鉛直方向以外の影 響を受ける可能性がある。
[0037] このため、加速度センサ 1として 1軸の加速度センサを用いる場合には、使用者の 歩行、走行に応じた鉛直方向の体動をできるだけ正確に検出できるようにするために 、例えば、使用者の腰部に、加速度の検出方向が鉛直方向となるように当該 1軸の 加速度センサを装着しなければならないといった制約が生じる。しかし、このような制 約を遵守した場合には、 1軸の加速度センサからの検出出力は加速度の鉛直成分と して用いることができる。この場合、 1軸の加速度センサ自体が鉛直成分抽出部 2とし ての機能を有することになる。
[0038] これに対して、加速度センサ 1として互いに直交する軸により構成される 2軸や 3軸 の加速度センサを用いることにより、加速度センサ 1の装着位置や装着方向に柔軟 性を持たせることが可能になる。しかし、加速度センサ 1として多軸の加速度センサを 用いる場合には、多軸の検出出力から鉛直成分を抽出する必要が生じる。この実施 の形態の体動検出装置において、加速度センサ 1は、例えば 3軸の加速度センサが 用いられている。このため、加速度センサ 1の後段に鉛直成分抽出部 2を設けている
[0039] [鉛直成分抽出部 2について]
鉛直成分抽出部 2は、(1)当該多軸の加速度センサ 1からの検出出力を用いて、そ の重力場における重力加速度ベクトルを推定し、 (2)この重力加速度ベクトルの推定 結果に基づ 、て、同じ加速度センサ 1の検出出力から鉛直方向の信号成分を抽出 する処理を行う。
[0040] このように、多軸の加速度センサ 1からの全ての軸についての検出出力を総合的に 用いて、鉛直成分を抽出することによって、作用軸を推定する必要がなぐユーザー に対する加速度センサ 1の装着位置や装着方向に左右されること無ぐユーザーの 鉛直方向の運動を正確に検出することができるようにしている。
[0041] ここで、加速度センサ 1として、 3軸の加速度センサを用いる場合について具体的に 説明する。加速度センサ 1を X軸、 Y軸、 Z軸の 3軸とし、当該加速度センサ 1から得ら れるある時点 nの加速度ベクトル anを図 2の(1 1)式に示すように、 axn (X軸成分) 、 ayn (Y軸成分)、 azn(Z軸成分)で表すものとする。この図 2の(1— 1)式に示した 加速度ベクトル (加速度ベクトルのデータ系列) an力 重力加速度ベクトル gの推定を 行うと共に、体動検出をも行うようにする。
[0042] 具体的には、重力加速度ベクトル gの推定は、より簡便には加速度ベクトル anの各 軸の移動平均値を算出し、その平均ベクトルを重力加速度ベクトル gとする方法があ る。この場合、体動による信号成分の影響を少なくするために、十分長い区間で移動 平均の計算を行うことが望ましい。また、加速度ベクトル anの各軸の値を、最小二乗 法などを用いて解析し、重力加速度ベクトル gの算出を行うようにする方法を用いるよ うにしてもよい。
[0043] 加速度ベクトル anを用いて重力加速度ベクトル gを推定した結果を、図 2の(1— 2) 式に示すように、 gx(X軸成分)、 gy(Y軸成分)、 gz (Z軸成分)で表すものとする。こ の場合、加速度ベクトル anの鉛直成分 vnは、図 2の(1 3)式に示す演算により求め ることができる。すなわち、加速度ベクトル anの鉛直成分 vnは、図 2の(1 3)式に示 したように、重力加速度ベクトル gの内積と加速度ベクトル anの積を重力加速度べタト ル gの絶対値 (大きさ)で割り算することにより求めることができる。
[0044] このように、 3軸の加速度センサ 1により検出される加速度ベクトル anと当該加速度 ベクトルから求められる重力加速度ベクトル gとから、鉛直成分 vnを演算により正確に 求めることができる。すなわち、 3軸の加速度センサ 1からの検出出力を総合的に用 い、その中から鉛直成分のみを数値計算により分離することにより、ユーザーの鉛直 方向の体動を正確に検出できるという着想に基づいてこの発明はなされている。
[0045] また、重力加速度ベクトル gの 3次元空間における偏角を求めたうえで、加速度べク トル anを回転させることによつても同様の計算ができる。すなわち、重力加速度べタト ル gについての偏角 0、 φを、図 3の(2— 1)、 (2— 2)に示すときに、図 3の(2— 3) 式により計算される xnは、加速度ベクトル anの鉛直成分であり、 a'xnは、鉛直成 分 vnと一致する。また、ベクトル ynとベクトル a'znの内積は、重力加速度ベクトル g を法線ベクトルとする平面への加速度ベクトル anの正射影である。
[0046] つまり、 3軸の加速度センサ 1からの検出出力により得られる加速度ベクトル anを鉛 直成分と水平成分に成分分解できるため、鉛直成分のみならず、水平成分を解析す ることによって水平方向の体動をも検出することができる。具体的には、水平ベクトル の長さ hnは、図 4の(3— 1)式や図 4の(3— 2)式によって求めることができる。
[0047] このように、重力加速度ベクトルの偏角を考慮した演算式を用いる場合には、ユー ザ一の鉛直方向の体動と水平方向の体動とを、比較的に簡単に、し力も正確に求め ることがでさる。 [0048] なお、ここでは 3軸の加速度センサ 1を用いる場合を例にして説明した力 これに限 るものではない。この発明の基本的なコンセプトは、 2軸の加速度センサを用いる場 合にも 3軸の加速度センサを用いる場合と同様にして適用可能である。
[0049] すなわち、 2軸の加速度センサにより検出する加速度ベクトル anと重力加速度べク トル gとを、図 5の(4 1)式 (加速度ベクトル)、図 5の(4 2)式(重力加速度ベクトル )のように表すものとすると、 3軸の加速度センサを用いる場合と同様に、図 2の(1— 3 )式にしたがって鉛直成分を計算することができる。
[0050] また、重力加速度ベクトル gの偏角 0を、図 5の(4 3)式に示すように表すものとす ると、図 5の(4 4)式にしたがって、加速度ベクトル anについて、その鉛直成分^ x nとこれに直交する水平成分^ ynに成分分解でき、^ xnと鉛直成分 vnとは一致する
[0051] このように、図 5に示した各式、および、図 2に示した(1 3)式により、 2軸の加速度 センサを用いた場合であっても、ユーザーの鉛直方向の体動を正確に検出すること 力 Sできるし、また、重力加速度ベクトル gの偏角 Θを考慮するようにした場合には、ュ 一ザ一の水平方向の体動をも正確に検出することが可能となる。
[0052] 図 6は、ユーザーに装着するようにした 3軸の加速度センサにより、当該ユーザーが 歩行など運動を行っている場合に、サンプリング周波数 50Hzで 4秒間、加速度デー タを取得した場合における、この取得した加速度データと、これを上述したこの発明 の基本的なコンセプトにしたがって成分分解した場合に得られるデータとをグラフ化 して示したものである。図 6において、横軸は時間(ミリ秒)、縦軸は重力加速度 (G) である。また、後述する図 7、図 8、図 10、図 12、図 14においても、横軸は時間(ミリ 秒)、縦軸は重力加速度 (G)である。
[0053] すなわち、図 6Aは、 3軸の加速度センサからの加速度データのグラフであり、図 6B は、 3軸の加速度データ力 計算した加速度ベクトルの長さ(大きさ)のグラフであり、 図 6Cは、図 2〜図 4を用いて説明した方法によって、 3軸の加速度データから計算し て得られた鉛直成分のグラフである。また、図 6Dは、図 2〜図 4を用いて説明した方 法によって、 3軸の加速度データ力も計算して得られた水平成分のグラフである。
[0054] そして、この図 6に示したグラフは、加速度データの検出時において、ユーザーは 主に鉛直方向の運動を行っていた力 約 80サンプル目、約 100サンプル目および 1 70サンプル目に水平方向の運動が発生し、これがノイズ成分として存在している場 合を示している。
[0055] しかし、ノイズ成分は水平成分に存在するものであるために、この発明の基本的な コンセプトにしたがって、加速度データ (加速度ベクトル)を鉛直成分(図 6C)と水平 成分(図 6D)とに成分分解すれば、鉛直成分からは水平成分のノイズを除去すること ができるので、ユーザーの鉛直方向の体動を正確に検出できることが分かる。もちろ ん、水平成分からは鉛直成分のノイズを除去することができるので、ユーザーの水平 方向の体動を正確に検出できる。すなわち、成分分解することにより、ノイズの低減 効果がある。
[0056] なお、ここでは、演算により加速度ベクトルの鉛直成分を抽出するようにした力 こ れに限るものではない。例えば、各軸の加速度ベクトルの長さ(大きさ)を演算により 求めて、一番長さのある加速度ベクトルを使用者の体動を一番反映した加速度べタト ルの鉛直成分として用いるようにしたり、あるいは、使用者の体動を一番反映している と推測される軸の検出出力(加速度ベクトル)を鉛直成分として用いるようにしたりして ちょい。
[0057] しかし、加速度ベクトルの長さに応じて鉛直成分を特定したり、あるいは、鉛直成分 に応じた軸を推定したりする場合には、多軸の加速度センサ 1の使用者に対する装 着位置や装着方向がある程度制限される場合もある。しかし、演算によって加速度べ タトルの鉛直成分を抽出するようにした場合には、上述もしたように、多軸の加速度セ ンサ 1の使用者に対する装着位置や装着方向の制約を受けることは無いので、体動 検出装置の使用者に対する装着の自由度を向上させることができる。
[0058] [高域 Z低域分解部 3について]
次に、高域 Z低域分解部 3の機能及び動作について説明する。上述したように、鉛 直成分抽出部 2において抽出された鉛直成分 vnを関数 x(n)で表すものとする。鉛 直成分抽出部 2からの加速度ベクトルの鉛直成分 x(n)は、高域 Z低域分解部 3に供 給される。高域 Z低域分解部 3は、例えば LPF (Low Pass Filter)の構成とされ、加速 度ベクトルの鉛直成分 x(n)を帯域分割し、高域成分 xh(n)と低域成分 xl(n)とに分 離する。
[0059] このとき、高域 Z低域分解部 3の LPFとしての特性としては、歩行や走行による加 速度の主要成分が含まれる 2Hz〜4Hzが通過域となって 、ることが望まし 、。また、 鉛直成分 X (n)と高域成分 xh (n)および低域成分 xl (n)の位相は揃って 、ることが望 ましい。
[0060] このように、鉛直成分 x(n)について、高域成分 xh (n)と低域成分 xl (n)とに分離す るのは、上述もしたように、低域成分 xl(n)には使用者の鉛直方向の体動に応じて変 化する成分が多く含まれ、高域成分 xh (n)にはノイズ成分が混入し易いためである。 このようにして高域 Z低域分解部 3にお 、て分解された高域成分 xh (n)と低域成分 x Kn)とは、ピーク検出 Z判定処理部 4に供給される。
[0061] [ピーク検出 Z判定処理部 4の機能と動作にっ 、て]
ピーク検出 Z判定処理部 4は、上述もしたように、加速度ベクトルの鉛直成分 x (n) の低域成分 xl(n)に基づいて、ピーク位置の候補を検出すると共に、このピーク位置 の候補を基準に用いた所定範囲毎に低域成分のエネルギーと高域成分のエネルギ 一の成分比に基づ!、てピーク位置の候補を特定する。
[0062] このように、低域成分のエネルギーと高域成分のエネルギーの成分比を用いる理 由について説明する。図 7、図 8は、この実施の形態の体動検出装置を使用者の異 なる部位に装着した場合に検出される各帯域の加速度信号のグラフを示す図である 。具体的には、図 7は、この実施の形態の体動検出装置を使用者の腰部に固定して 装着して用いた場合に検出される加速度信号のグラフを示す図であり、図 8は、この 実施の形態の体動検出装置を使用者のズボンのポケットに入れて用いた場合に検 出される加速度信号のグラフを示す図である。
[0063] 図 7、図 8とも、使用者の歩行時における 3軸の加速度センサ 1において検出された 加速度信号を 50Hzでサンプリングすることにより計測したものである。また、図 7、図 8とも最上段(図 7A、図 8A)のグラフが 3軸の加速度ベクトル力 抽出した鉛直成分、 すなわち x (n)であり、中段(図 7B、図 8B)に示したグラフが鉛直成分 x(n)の内の低 域成分 xl (n)であり、最下段(図 7C、図 8C)のグラフが鉛直成分 x (n)の高域成分 xh (n)である。 [0064] そして、図 7に示したように、この実施の形態の体動検出装置を使用者の腰部に装 着した場合における 3軸の加速度センサ 1の検出出力においては、中段に示した低 域成分 xl (n)に周期的な変化が集中し、最下段に示した高域成分 xh (n)には周期 的な変化はほとんど無いことからほぼ歩行による上下動に対応する成分のみが加速 度信号として計測されて!ヽることが分かる。
[0065] したがって、この実施の形態の体動検出装置を使用者の腰部に装着して使用する ようにした場合には、使用者の鉛直方向の体動に応じた成分のみを効率よく計測す ることができるので、鉛直成分 X (n)の低域成分 xl (n)の波形を閾値判定あるいはピ ーク検出することによって歩行や走行などの体動を検出することができ、これをカウン ト処理すれば、正確に歩数のカウント (計数)が可能な歩数計が実現出来る。
[0066] これとは対称的に、図 8に示したように、この実施の形態の体動検出装置を使用者 のポケットに入れて用いるようにした場合における 3軸の加速度センサ 1の検出出力 においては、歩行に応じた周期的な上下動以外の振動成分 (ノイズ成分)が、中段に 示した低域成分 xl (n)〖こも、また、最下段に示した高域成分 xh (n)にも現れており、 特に、最下段に示した高域成分 xh (n)に多く現れている。したがって、低域成分につ いて閾値判定やピーク検出を行うようにしても、歩行以外のノイズ成分を歩行に応じ た体動として誤検出してしまう可能性が非常に高くなる。
[0067] ところで、図 8において、区間 A、区間 B、区間 Cに示すように、中段に示した低域成 分のグラフにおいて、振幅が 1Gを下回っている区間 (領域)を区切った場合、周期的 な運動である歩行による上下動に対応する波形は、区間 Aの波形と区間 Bの波形で あり、区間 Cの波形は、周期的なものではなぐ歩行以外の振動、すなわちノイズ成 分である。そして、図 8に示したように、区間 A及び区間 Bでは高域成分が弱いのに 対して、区間 Cでは高域成分が強い。
[0068] このように、周期的な運動である歩行による使用者の体の上下動に応じた成分は、 低域成分 xl (n)に顕著に現れ、ノイズ成分は高域成分 xh (n)に顕著に現れる。そこ で、ピーク位置の候補として検出したそれぞれのピーク位置の候補毎に、ピーク位置 の候補の前後に予め決められた時間幅を有する所定領域を定める。
[0069] 例えば、ピーク位置の候補の前 Mサンプル区間と後 Mサンプル区間(Mは 1以上の 整数)というように所定範囲を定めればよい。なお、図 7、図 8においては、比較的にノ ィズの影響を受け難 、極小値をピーク値として検出し、これに対応するピーク位置を ピーク位置の候補として特定するようにして!/、る。
[0070] このようにして定めたピーク位置の候補を含む各所定領域にぉ 、て、低域成分 xl ( n)のエネルギーと高域成分 xh(n)のエネルギーとのエネルギー比 dを求める。図 9は 、エネルギー比 dを求めるための計算方式を説明するための図である。加速度べタト ルの鉛直成分 X (n)の高域成分 xh (n)のエネルギー ehは、図 9の(5— 1)式によって 求めることができる。また、加速度ベクトルの鉛直成分 x(n)の低域成分 xl (n)のエネ ルギー elは、図 9の(5— 2)式によって求めることができる。
[0071] なお、図 9の(5— 1)式と(5— 2)式において、 nlはピーク位置の候補を含む所定 領域の開始位置を示し、 n2はピーク位置の候補を含む所定領域の終了位置を示し て 、る。そして、低域成分 xl (n)のエネルギーと高域成分 xh (n)のエネルギーとのェ ネルギー比 dは、図 9の(5— 3)式によって求めることができる。
[0072] このようにして求めたエネルギー比 dが、予め決められた閾値 Dよりも小さい場合に 、その所定領域に含まれるピーク位置の候補として検出された位置をピーク位置の 候補として特定する。すなわち、ピーク位置の候補として特定した位置は、使用者の 歩行あるいは走行といった運動に応じてピークになったものと判断し、このピーク位 置の候補をカウントすることにより、歩行や走行を行った使用者の歩数を正確にカウ ントすることができる。すなわち、ある閾値 Dを定めたうえで、エネルギー比 d<閾値 D のときに限って歩行あるいは走行と判定することにより、歩行や走行以外の振動成分 の影響を排除できる。
[0073] し力しながら、上述のように、歩行や走行といった使用者の体動に応じた変化を検 出するために、加速度ベクトルの鉛直成分 x (n)の低域成分 xl(n)に基づいて、ピー ク位置の候補を検出し、低域成分 xl(n)のエネルギー elと高域成分 xh(n)のェネル ギー ehとのエネルギー比 dを考慮してピーク位置の候補を特定するようにしても、ま だ誤判定をする場合があると考えられる。
[0074] 図 10は、歩行や走行といった使用者の周期的な体動以外の影響により、多軸の加 速度センサ 1により検出される加速度ベクトルの鉛直成分 x (n)の低域成分 xl (n)にも 不規則な信号が混入した場合のグラフ (不規則波形)を示す図である。例えば、転ん だり、体動検出装置を落としたり、あるいは、使用者がジャンプするなど、歩行や走行 以外の不規則な動作が発生した場合には、図 10の中段の低域成分 xl(n)のグラフ に表されるように、歩行や走行以外の不規則な体動成分が混入する場合がある。
[0075] このように、歩行や走行以外の不規則な体動成分が、加速度ベクトルの鉛直成分 X
(n)の低域成分 xl (n)に混入していた場合には、これを周期的な運動である歩行や 走行に応じた体動と誤判定をしてしまう場合があると考えられる。そこで、ある区間に お!、て信号波形が似て 、るかどうかを比較して、波形が似て 、れば歩行や走行であ ると判定をするようにする。
[0076] 例えば、歩行や走行以外の不規則な動作が発生した場合の加速度ベクトルの鉛直 成分 x (n) (図 10A)、低域成分 xl (n) (図 10B)、高域成分 xh (n) (図 IOC)の例のグ ラフである図 10において、ピーク位置として、位置 A、位置 B、位置 Cを検出すること ができる。これらに関して、位置 Aと位置 B、位置 Bと位置 C、位置 Aと位置 Cというよう にペアを作り、そのペア間において、それぞれのピーク位置の周辺の波形について マッチングを取るようにする。
[0077] この場合、図 10からも明らかなように、位置 Aの近傍の波形、位置 Bの近傍の波形 、位置 Cの近傍の波形については、相互に類似するものがないために、周期的な体 動である歩行や走行の体動とは異なるものと判断し、歩行や走行時の体動ではな 、 ものとして判断することが可能である。
[0078] これに対して、図 7や図 8に示したように、加速度ベクトルの鉛直成分 x (n)、低域成 分 xl (n)、高域成分 xh (n)が検出された場合において、ピーク位置の候補として特定 された位置を含む所定区間を設定し、その前、あるいは、後ろの 1つ以上の所定区 間の波形とマッチング処理した場合には、類似する場合があるためにピーク位置の 候補として特定された位置を、正式なピーク位置として決定することができる。
[0079] 具体的には、図 7に示した例の場合には、ノイズの影響が少な 、ように、低域成分 X 1 (n)のグラフにお 、て極小値をピーク値として取るようにすれば、隣り合うピーク位置 を含む波形の類似度は非常に高 、ので、各ピーク位置の候補として特定された位置 をピーク位置として決定することができる。また、図 8に示した例の場合においても、ノ ィズの影響が少ないように、低域成分 xl (n)のグラフにおいて極小値をピーク値とし て取るようにしているので、ピーク位置の候補の 1つおきの当該ピーク位置を含む波 形の類似度が非常に高ぐ各ピーク位置の候補として特定された位置をピーク位置と して決定することができる。
[0080] このように、ピーク位置の候補を特定し、その候補を含む所定区間の波形にっ 、て 、比較する(マッチングを取る)他の所定区間を定め、両所定区間同士の波形を比較 して類似度が高い場合に、比較対象の所定区間に含まれるピーク位置の候補をピー ク位置として決定することができる。
[0081] なお、比較する他の所定区間は、例えば、隣り合う所定区間同士、あるいは、 1つ おきの所定区間同士、あるいは、 2つおきの所定区間同士というように、適宜定めるこ とができる。あるいは、前 (過去方向)の 1つ以上の所定区間の波形と比較し、 2っ以 上の所定区間との間で類似度が高い場合に、比較元の所定区間のピーク位置の候 補をピーク位置と決定するなど、比較する所定区間の位置や数は任意に設定するこ とが可能である。
[0082] このように、この実施の形態の体動検出装置において、ピーク検出 Ζ判定処理部 4 は、単に加速度ベクトルの鉛直成分 X (η)の低域成分 xl (n)の情報 (波形)からピーク 位置の候補を検出し、低域成分のエネルギー elと高域成分のエネルギー ehとのエネ ルギー比 dを考慮してピーク位置の候補を特定するに止まらず、ピーク位置の候補を 含む所定区間同士の波形比較 (波形マッチング)を行うことによって、歩行や走行と 言った使用者の周期的な体動を正確に検出することができるようにしている。
[0083] そして、使用者の歩行や走行などの周期的な体動を正確に検出することができるの で、使用者の体動をカウントすることにより、正確に使用者の歩行時や走行時におけ る歩数のカウントが可能な歩数計を実現することができる。なお、精度を若干落として よければ、上述もしたように、波形マッチングを行うことなぐエネルギー比 dを用いて 特定するようにしたピーク位置の候補に従って歩数をカウントしてもよ 、。
[0084] [ステップ位置解析部 5の機能と動作にっ 、て]
上述した加速度センサ 1、鉛直成分抽出部 2、高域 Z低域分解部 3、ピーク検出 Z 判定処理部 4の各部の機能により、使用者の身体に装着される加速度センサ 1から の検出出力に基づいて、使用者の歩行や走行に応じて使用者の身体に生じる鉛直 方向の体動を正確に検出することができる。このようにして検出した体動をカウントす れば、使用者の歩数を正確にカウントすることも可能となる。
[0085] し力しながら、使用者の動作状態としては、歩行や走行を行っている動作状態ばか りではなぐ静止している状態もあれば、静止状態でもなく動作態でもない不定な状 態である場合もある。このため、使用者の動作状態を正確に把握できれば、歩数を力 ゥントするなどのように、体動をカウントする場合においても、静止状態や不定状態に あるときにおいても体動のカウントを行うといった誤動作をさらに防止することが可能 である。
[0086] そして、歩行や走行など運動を行っている動作状態にあることが分ければ、その期 間においては使用者の体動を正確に検出し、使用者の体動ピッチ (体動テンポ)を 正確に把握し、使用者の他動ピッチに合わせた機器制御を行うことも可能となる。そ こで、この実施の形態の体動検出装置において、ステップ位置解析部 5は、ピーク検 出 Z判定処理部 4にお 、て決定されたピーク位置を示す情報の供給を受けて、これ に基づき、使用者の動作状態を正確に把握することができると共に、使用者が歩行 や走行などの周期的な運動を行って 、る場合にぉ 、ては、その体動ピッチをも正確 に検出することができるようにして!/、る。
[0087] 図 11は、使用者の動作状態の遷移を説明するための図である。図 11に示すように 、使用者の動作状態として、「不定」、「静止」、「歩行 Z走行」の 3つの状態があるもの とする。ここで、「静止」は使用者の体の動きが全く無い状態を、「歩行/走行」は使用 者が歩行または走行している状態を、「不定」は上記 2つの状態以外の状態を意味す る。また、初期状態は「不定」から開始するものとする。また、「静止」状態から「歩行 Z 走行」状態に遷移する場合、あるいは、「歩行 Z走行」状態から「静止」状態に遷移す る場合には、「不定」状態を経由するものとする。
[0088] そして、ステップ位置解析部 5においては、ピーク検出 Z判定処理部 4からの決定 されたピーク位置情報に基づいて、まず、使用者の動作状態を判定する。ここでは、 図 11において、遷移 A〜遷移 Gのそれぞれの遷移状態を判別することにより、使用 者の動作状態を判別する。 [0089] すなわち、「不定」状態が維持される場合を遷移 A、「不定」状態から「歩行 Z走行」 状態に変化する場合を遷移 B、「歩行 Z走行」状態が維持される場合を遷移 C、「歩 行 Z走行」状態から「不定」状態に遷移する場合を遷移 D、「不定」状態から「静止」 状態に遷移する場合を遷移 E、「静止」状態から「不定」状態に遷移する場合を遷移 F、「静止」状態が維持される場合を遷移 Gというように、各状態が維持される場合に っ 、ても 1つの遷移状態として動作状態の遷移を細かく把握するようにして!/、る。
[0090] 次に、各遷移状態の判別のための条件について説明する。各遷移状態の判別は、 上述もしたように、ピーク検出 Z判定処理部 4にお 、て決定されたピーク位置に基づ いて行う。以下においては、「決定されたピーク位置」を「検出されたステップ」とも言う 。すなわち、ピーク検出 Z判定処理部 4において決定されたピーク位置力 使用者の 鉛直方向の体動に応じたステップとみなすようにしている。そして、(1)以下に説明す る遷移 Bである場合の条件も、遷移 Cである場合の条件も満たさないときには、遷移 Aの状態であると判定し、「不定」状態を維持する。
[0091] また、(2)遷移 Aの状態であると判定した後に、検出されたステップ (決定されたピ ーク位置)のうち新しい方力も数個について、隣接するステップの時間間隔を計算し 、時間間隔が一定である場合には、「不定」から「歩行/走行」へ遷移する遷移 Bの状 態であり、使用者の動作状態としては、「歩行 Z走行」状態になったと判定する。この 場合、基準ステップ間隔 (基準ピッチ Ps)を算出しておく。なお、基準ステップ間隔は 、後述もするが、遷移 Bであることの判定に用いたステップ間の間隔の平均値である。
[0092] また、(3)遷移 Bまたは遷移 Cと判定した後に、新しく検出されたステップとその直前 のステップとの間隔が、基準ステップ間隔 (基準ピッチ Ps)の整数倍に対してある誤差 範囲内である場合には、遷移 Cの状態であると判定し、「歩行/走行」状態を維持する 。また、遷移 Bまたは遷移 Cであると判定した後において、上述した(3)の遷移 Cであ ると判定するための条件を満たさない場合には、遷移 Dの状態であり、使用者の動作 状態としては、「不定」状態であると判定する。
[0093] また、(5)遷移 Aまたは遷移 Dであると判定した後において、一定時間以上、ステツ プが検出されない場合 (ピーク位置が決定されない場合)には、「不定」から「静止」へ 遷移する遷移 Eの状態であり、使用者の動作状態としては「静止」状態であると判定 する。また、(6)遷移 Eまたは遷移 Gであると判定した後において、ステップの検出が された場合 (ピーク位置が決定された場合)には、「静止」から「不定」へ遷移する遷移 Fの状態であり、使用者の動作状態としては、「不定」状態になったと判定する。
[0094] また、(7)遷移 Eまたは遷移 Gと判定された後において、ステップが検出されない場 合 (ピーク位置が決定されない場合)には、遷移 Gであると判定し、使用者の動作状 態としては。「静止」状態を維持する。
[0095] なお、上述した(2)の遷移 Bか否かの判定において、ステップ間隔が一定かどうか の判定基準は、ステップ間隔の分散や標準偏差があるかどうかで判定してもよいし、 最大値と最小値の差が閾値以下であるかどうかで判定するのでもよ 、。また基準ステ ップ間隔は、上述もしたように判定に用いたステップの時間間隔の平均を用いたり、 ある ヽは、判定に用 、たステップの時間間隔の中間値を用 、たりしてもよ 、。
[0096] 図 12は、図 1に示した体動検出装置の加速度センサ 1、鉛直成分抽出部 2、高域 Z低域分解部 3、ピーク検出 Z判定処理部 4の各部の機能を用いて、加速度センサ 1からの検出出力からステップ検出(ピーク位置の決定)をした場合の波形を示す図 である。この図 12においても、ピークは極小値側で取るようにしている。図 12Aに示 す波形は、「不定」状態力 4つのステップ (決定されたピーク位置)が検出された場 合の波形を示す、図 12Bは、図 12Aの状態から、その直後のステップが検出された 場合の波形を示している。このような図 12を用いて、状態遷移の例について具体的 に説明する。
[0097] まず、図 12Aに示した波形において、最初は「不定」の状態であったとする。この後 、 4つのステップが検出(2つのピーク位置が決定)されているため、ステップ位置解 析部 5は、隣接するステップ区間 Tl、 Τ2、 Τ3の各ステップ間隔を算出し、それらを 比較する。この例の場合、ステップ区間 Tl、 Τ2、 Τ3の各区間の間隔はほぼ一定で あるので、この場合には、「不定」から「歩行 Ζ走行」へ遷移する遷移 Βの状態である と判定することができる。
[0098] このとき、ステップ位置解析部 5は、図 13に示す (6— 1)式により、基準ピッチ (基準 ステップ間隔) Psを算出する。なお、図 13に示した(6— 1)式は、図 12Aの状態にあ るときの基準ピッチ Psを算出するための式であり、ステップ区間が増えれば、分子の 加算される区間と、分母の区間数が変わることになる。
[0099] そして、図 12Aに示した状態から、図 12Bに示した状態に変化した場合には、最新 のステップ区間 T4のステップ間隔と、図 13の(6— 1 )式によって算出した基準ピッチ Psとを比較する。この場合、ステップ区間 T4のステップ間隔と基準ピッチ Psとはほぼ 同じであり、差分は誤差範囲内であると判定されることにより、遷移 Cの状態であると 判定され、「歩行 Z走行」の状態が維持されることになる。
[oioo] なお、「歩行 Z走行」状態における歩行時や走行時におけるピッチの周期 (ステップ 間隔)は、基準ピッチを代表値として用いてもよいし、ある時間区間内のステップ間隔 の平均値や、一定数のステップ間隔の平均値を用いてもよ!、。
[oioi] また、上述の説明にお 、ては、「歩行」状態と「走行」状態とを区別せずに、「歩行 Z 走行」状態というように 1つの状態として把握するようにした。しかし、これに限るもので はな ヽ。「歩行」状態と「走行」状態とをそれぞれ別の状態として把握できるようにして ちょい。
[0102] 歩行と走行を分ける基準として、上述した(2)の遷移 Bの判定の条件にカ卩えて、(A )ステップ間隔がある一定未満の場合は「走行」状態と判定し、それ以外の場合には「 歩行」と判定するようにしたり、あるいは、(B)加速度波形のピーク値が一定以上の場 合は「走行」状態であると判定し、それ以外の場合には「歩行」状態であると判定した りすることが可能である。もちろん、上述のように、上述 (A)または(B)の何れかの条 件を用いてもょ 、し、上述の (A)と(B)との両方の条件用いるようにしてもよ!、。
[0103] また、上述したように、「歩行」状態と「走行」状態とを全く別の動作状態として定義し てもよいが、状態遷移は図 11にしたがうものとして、「歩行 Z走行」状態の属性とし、 上述した (A)、(B)の一方または両方の条件 (基準)を用いて、「歩行」と「走行」とを 分けて把握するようにしてもょ 、。
[0104] また、遷移 Cであるか否かの判定にぉ 、ては、ステップ間隔と基準ピッチ(基準ステ ップ間隔) Psとをそのまま比較するのではなぐ基準ピッチ Psの整数倍を基準に判定 するようにした。このようにすることによって、ステップの検出漏れがあった場合の意図 しない状態遷移を防止することに役立つ。
[0105] 例えば、この実施の形態の体動検出装置においては、上述もしたように、ピーク検 出 Z判定処理部 4において、ピーク位置の候補を検出して、高域と低域のエネルギ 一を考慮してピーク位置の候補を特定し、さらに、このピーク位置の候補を含む所定 間隔毎に波形比較を行うことによって、ピーク位置を決定するようにしている。
[0106] しかし、実際には使用者の歩行が続行しているのにも関わらず、一時的に加速度 波形が乱れると周囲の波形とのマッチングが取れず、検出漏れが起こる可能性があ る。このため、上述もしたように、ステップ位置解析部 5において、遷移 Cの判定を行う 場合に、「基準ピッチの整数倍を基準に判定」することにより、前段のピーク検出 Z判 定処理部 4において、歩行ステップ検出漏れ (ピーク位置の決定漏れ)が起こっても 安定した状態推定を行うことが可能となる。
[0107] 例えば、図 14は、「歩行/走行」状態において、 X印(ばつ印)で示した 2つのステツ プ (ピーク位置)の検出漏れ (カウント漏れ)が生じ、一番最近に検出されたステップ( ピーク位置)とその直前に検出されたステップ (ピーク位置)との間隔が 3ステップ分に なる場合の例を示している。
[0108] この図 14に示した例の場合には、遷移 Cである力否かを判断するために条件にお いて、新しく検出されたステップとその直前のステップとの間隔力 基準ステップ間隔 (基準ピッチ Ps)の整数倍に対してある誤差範囲内である場合には、遷移 Cの状態で あると判断され、「歩行/走行」の状態は維持されることになる。
[0109] 従って、単純に検出されたステップを 1つづつ積算するのではなぐステップ (ピーク 位置)が検出されない場合であっても、「歩行 Z走行」状態が維持されている場合に おいては、基準ピッチに対する倍数を積算の対象とすることで、より正確な歩数計測 も可能となる。したがって、この発明を歩数計に適用することで、歩数計の精度向上 にも役立つ。
[0110] このように、この実施の形態の体動検出装置は、加速度センサ 1、鉛直成分抽出部 2、高域 Z低域分解部 3、ピーク検出 Z判定処理部 4、ステップ位置解析部 5の各部 が有機的に機能することによって、使用者の歩行や走行といった体動を正確に検出 、これを把握することができる。
[0111] [具体的な装置への適用 ]
次に、この発明による装置、方法、プログラムを具体的な装置に適用した場合につ いて説明する。以下においては、この発明を歩数計に適用した場合と、音響再生装 置に適用した場合とに分けて説明する。
[0112] [歩数計への適用について]
まず、この発明の装置、方法、プログラムの一実施の形態を歩数計に適用した場合 について説明する。図 15は、この実施の形態の歩数計 100を説明するためのブロッ ク図である。図 15に示すように、この実施の形態の歩数計 100は、制御部 110に対し て、 3軸の加速度センサ 101が AZD変換部 102を介して接続されると共に、表示部 103、操作部 104が接続されて形成されたものである。
[0113] 3軸の加速度センサ 101は、 1軸の加速度センサを互いに直交する X軸、 Y軸、 Z軸 の 3軸に配置して形成されたもの、あるいは、直交する 3軸の加速度センサがひとつ のパッケージに封入されて形成されたものの 、ずれを用いてもょ ヽ。 3軸の加速度セ ンサ 101からの X軸、 Y軸、 Z軸のそれぞれの検出出力(アナログ出力)は、 AZD変 換部 102に供給され、ここで制御部 110にお 、て処理可能な形式のデジタルデータ に変換されて制御部 110に供給される。
[0114] 制御部 110は、この実施の形態の歩数計 100の各部を制御するものであり、図 15 に示すように、 CPU (Central Processing Unit) 111、 ROM (Read Only Memory) 11 2、 RAM (Random Access Memory) 113、不揮発性メモリ 114が、 CPUバス 115を通 じて接続されてマイクロコンピュータの構成とされたものである。
[0115] ここで、 CPU111は、各種のプログラムを実行し、各部に供給する制御信号を形成 したり、各種の演算を行ったりするなど、制御部 110における処理や制御の主体とな るものである。 ROM112は、 CPU111が実行する各種のプログラムや処理に必要に なるデータを記憶保持する。
[0116] また、 RAM113は、処理の途中結果を一時記憶するなど、主に作業領域として用 いられるものであり、 AZD変換部 102を通じて供給される加速度センサ 101からの 数値ィ匕された加速度データを記憶保持するバッファメモリなどとしても用いられる。不 揮発性メモリ 114は、例えば、 EEPROM (Electrically Erasable and Programmable R OM)やフラッシュメモリなどの電源が落とされても記憶データが消滅することのな!/、メ モリであり、電源が落とされても保持しておくべきデータ、例えば設定されたパラメ一 タゃ追加されたプログラムなどを記憶保持する。
[0117] また、制御部 110に接続された表示部 103は、表示制御回路を含み、例えば、 LC D (Liquid Crystal Display)、有機 EL (Electro Luminescence)ディスプレイ、 CRT (Ca thode-Ray Tube)等の表示素子を備え、制御部 110の制御に応じて歩数のカウント 値や種々のガイダンス情報等を表示する。具体的には、表示部 103は、制御部 110 力 表示データの供給を受けて、この表示データ力 表示素子に供給する映像信号 を形成し、これを表示素子に供給することによって、制御部 110からの表示データに 応じた表示情報を表示素子の表示画面に表示する。
[0118] また、操作部 104は、リセットキーや各種のファンクションキーなどを備え、ユーザー 力ゝらの操作入力を受け付けて、これに応じた電気信号を制御部 110に供給すること ができるものである。これに応じて、制御部 110は各部を制御し、ユーザーの指示に 応じた処理を行うことができるようにして 、る。
[0119] そして、この実施の形態の歩数計 100において、制御部 110が、図 1に示す体動検 出装置における鉛直成分抽出部 2、高域 Z低域分解部 3、ピーク検出 Z判定処理部 4、ステップ位置解析部 5としての機能を実現すると共に、ステップ計測を行うステップ 計測部としての機能をも実現するようにしている。
[0120] すなわち、この実施の形態の歩数計 100の制御部 110は、加速度センサ 101から の検出出力に基づいて、図 1を用いて説明したように、鉛直成分抽出、鉛直成分の 高域 Z低域の分解、ピーク位置の候補の検出とピーク位置の候補の特定、波形マツ チングによるピーク位置の決定を行って、決定したピーク位置に基づ 、て歩数を計 測する。
[0121] さらには、決定したピーク位置に基づいて、使用者の動作状態を正確に把握すると 共に、歩行や走行などの体動ピッチを正確に把握し、使用者の動作状態が「歩行 Z 走行」である場合にぉ 、て、体動ピッチに基づ 、て使用者の歩数をより正確に計測 することちでさるよう〖こして ヽる。
[0122] 以下、図 16〜図 21のフローチャートを参照しながら、図 15に示したこの実施の形 態の歩数計の主に制御部 110において行われる種々の処理について詳細に説明す る。 [0123] [鉛直成分抽出処理について]
図 16は、図 15に示したこの実施の形態の歩数計 100の主に制御部 110において 実行される鉛直成分抽出処理を説明するためのフローチャートである。この実施の形 態の歩数計 100に電源が投入され、操作部 104を通じて歩数の計測を実行するよう にする指示入力を受け付けると、制御部 110の CPU111は、まず、図 16に示す処理 を実行する。
[0124] 制御部 110は、 AZD変換部 102を通じて供給される加速度データ (加速度べタト ル)を取得して RAMI 13にバッファリングする処理を開始し (ステップ S101)、取得し た加速度ベクトルを用いて重力加速度ベクトル gの初期値を算出する (ステップ S102 )。そして、サンプルカウンタ cをリセット(ゼロクリア)する(ステップ S 103)。このステツ プ S101〜ステップ S103までの処理が電源投入後のいわゆる初期処理に相当する
[0125] そして、制御部 110は、サンプルカウンタ cの値は、予め決められた値 Nより大きい か否かを判断する (ステップ S 104)。この実施の形態においては、演算量を削減する ために、重力加速度ベクトル gの再計算は、 N (Nは、 1以上の整数)サンプルおきに 行うようにしている。
[0126] ステップ S 104の判断処理において、サンプルカウンタ cが規定値 Nよりも大きいと 判断したときには、制御部 110は、重力加速度ベクトル gの再計算処理を行って (ステ ップ S105)、この後、サンプルカウンタ cに値 0をセットする(ステップ S 106)。すなわ ち、ステップ S 106の処理は、サンプルカウンタ cのリセット処理である。そして、最新 に再計算して得た重力加速度ベクトル gと、加速度ベクトル anとを用いて、図 2または 図 3を用いて説明したように演算により鉛直成分 vnを求める (ステップ S107)。
[0127] また、制御部 110は、ステップ S104において、サンプルカウンタ cの値は、値 Nより 大きくないと判断した場合には、ステップ S102において算出した重力加速度べタト ル gの初期値と最新の加速度ベクトル anとを用いて、上述もしたように、図 2または図 3を用いて説明したように演算により鉛直成分 vnを求める (ステップ S107)。
[0128] そして、ステップ S107の処理の後、制御部 110は、最新の加速度データを取得し( ステップ S108)、サンプルカウンタ cに 1加算して(ステップ S109)、ステップ S104か らの処理を繰り返す。このようにして、ステップ S 107で求められた加速度ベクトルの 鉛直成分が、高域 Z低域分離処理に用いられる。
[0129] なお、この実施の形態において、ステップ S 105の重力加速度ベクトル gの再計算 処理と、ステップ S 102の重力加速度ベクトルの初期値の算出処理は、基本的に同じ 処理であり、加速度ベクトルの各軸のデータの移動平均を取ることによって、重力加 速度の推定値として、重力加速度ベクトル gを算出する。
[0130] ステップ S 102、ステップ S 105で行われる処理を具体的に説明する。ここでは、現 在のサンプル位置を nlとして、過去 Mサンプル分の加速度データから重力加速度を 求める場合を考える。この場合において、 X軸の重力加速度ベクトル gxを求めには、 現在のサンプル位置よりも M+ 1サンプル前から nlまでの各サンプルにおける X軸 の加速度データ axnを加算して合計値を求め、この合計値を値 Mで割り算することに より、 X軸の重力加速度ベクトル gxが求められる。
[0131] 同様に、 Y軸の重力加速度ベクトル gyを求める場合には、現在のサンプル位置より も M+ 1サンプル前から nlまでの各サンプルにおける Y軸の加速度データ aynをカロ 算して合計値を求め、この合計値を値 Mで割り算することにより、 Y軸の重力加速度 ベクトル gyが求められる。また、 Z軸の重力加速度ベクトル gzを求める場合には、現 在のサンプル位置よりも M+ 1サンプル前から nlまでの各サンプルにおける Z軸の加 速度データ aznを加算して合計値を求め、この合計値を値 Mで割り算することにより、 Z軸の重力加速度ベクトル gzが求められる。
[0132] このようにして重力加速度ベクトルを算出する場合には、動きによる加速度成分が 平均化されてキャンセルされるように十分長い区間で平均を取るようにすることが望ま しい。しかし、あまり長すぎると装置の傾きに追従しなくなりため(装置の傾きを正しく 反映させることができなくなるために)、例えば、数秒間程度に設定するのが適当であ る。
[0133] そして、ステップ S107の鉛直成分の抽出処理においては、上述もしたように、図 2 に示した (1 3)式により、あるいは、図 3に示した(2— 3)式により、最新の加速度べ タトル anと重力加速度ベクトル gとに基づいて、鉛直成分 vnを求める(抽出する)こと ができる。 [0134] このように、制御部 110は、図 16に示した処理を実行することによって、鉛直成分 抽出部 2としての機能を実現するようにして 、る。
[0135] なお、図 16に示したように、この実施の形態の歩数計 100においては、重力加速 度の推定は、演算量削減のために Nサンプルおきに行うようにした力 これに限るも のではない。重力加速度の推定についても毎サンプルごとに行ってもよい。重力加 速度の推定のための演算処理は、各軸のデータの移動平均をとるものに限るもので はない。例えば、最小二乗法などを用いてもよい。
[0136] [鉛直成分の高域 Z低域分離処理にっ 、て]
そして、図 16に示した処理により抽出された加速度ベクトルの鉛直成分を高域成分 と低域成分とに分離するのであるが、これを上述もしたように、例えば、歩行や走行に よる加速度の主要成分が含まれる 2Hz〜4Hzの帯域の成分を低域成分として抽出し 、 4Hzより高い帯域の成分を高域成分として抽出するようにすれば、鉛直成分を予め 決められた帯域の高域成分と予め決められた帯域の低域成分とに分離することがで きる。
[0137] 図 16に示した処理により抽出される鉛直成分 x(n)には、使用者の歩行運動に伴う 上下動に対応したピークが現れ、特に、鉛直成分 x (n)の低域成分 xl(n)には、使用 者の歩行運動に伴う上下動に対応したピークが顕著に現れる。このため、鉛直成分 X (n)を分離することにより得られた低域成分 xl(n)に基づいて、ピーク位置の候補を 検出し、低域成分 xl (n)のエネルギーと高域成分 xh (n)のエネルギーとの比に基づ いて、ピーク位置の候補を特定する。
[0138] [ピーク位置の候補の検出及び特定処理について]
図 17は、図 15に示したこの実施の形態の歩数計 100の主に制御部 110において 実行されるピーク位置の候補の検出及び特定処理を説明するためのフローチャート である。この図 17に示す処理は、図 16を用いて説明した鉛直成分抽出処理により抽 出された鉛直成分 X (n)を分離することにより得られる低域成分 xl (n)と高域成分 xh ( n)とに基づいて、例えば、サンプリング毎に実行される。
[0139] まず、制御部 110は、上述したように、鉛直成分 x (n)から分離された低域成分 xl(n )に基づいてピーク位置の候補を探索 (検出)する処理を行う(ステップ S201)。具体 的には、低域成分 xl (x)について、例えば、サンプリング毎に、最新に取得した低域 成分データとその直前の低域成分データとを比較し、直前の低域成分データの方が 大きぐかつ、当該直前の低域成分データが所定値以上の場合には、当該直前の低 域成分データをピーク値の候補として検出すると共に、そのピーク値の位置をピーク 位置の候補として検出するようにする。
[0140] そして、制御部 110は、ピーク位置の候補が検出されたか否かを判断し (ステップ S 202)、検出されていないと判断したときには、この図 17に示す処理を終了して、次 のサンプリングのタイミングで、再度、この図 17に示す処理を実行するようにする。
[0141] また、ステップ S 202の判断処理において、ピーク位置の候補が検出されたと判断 したときには、制御部 110は、ピーク位置の候補近傍の低域成分 xl (n)を用いて、ピ ーク位置の候補を含む所定範囲における低域成分 xl(n)のエネルギー elを図 9の(5 — 2)式に従って算出する (ステップ S203)。また、同様にして、制御部 110は、ピー ク位置の候補近傍の高域成分 xh (n)を用いて、ピーク位置の候補を含む所定範囲 における高域成分 xh (n)のエネルギー ehを図 9の(5— 1)式に従って算出する (ステ ップ S 204)。
[0142] そして、制御部 110は、図 9の(5— 3)式に従って、ステップ S203において算出し た低域成分 xl (n)のエネルギー elとステップ S204にお!/、て算出した高域成分 xh (n) のエネルギー ehとのエネルギー比 dを算出し (ステップ S205)、算出したエネルギー 比 dが、予め決められた閾値 Dより小さいか否かを判断する (ステップ S206)。
[0143] ステップ S206の判断処理において、エネルギー比 dが、予め決められた閾値 Dより 小さくない、すなわち、高域成分 xh (n)にノイズが多いと判断した場合には、ピーク位 置の候補として検出したものは誤検出の可能性が高いので、検出したピーク位置の 候補をピーク位置の候補としては特定せずに、この図 17に示す処理を終了して、次 のサンプリングのタイミングで、再度、この図 17に示す処理を実行するようにする。
[0144] ステップ S206の判断処理において、エネルギー比 dが、予め決められた閾値 Dより 小さい、すなわち、高域成分 xh(n)にノイズが少ないと判断した場合には、ピーク位 置の候補として検出したものはピーク位置としての信憑性が高 、ので、検出したピー ク位置の候補をピーク位置の候補として特定して、ピーク位置の候補を例えば RAM 113などのピーク位置記録バッファに記録するようにし (ステップ S207)、この後、こ の図 17に示す処理を終了して、次のサンプリングのタイミングで、再度、この図 17に 示す処理を実行するようにする。
[0145] このようにして、この実施の形態の歩数計 100の制御部 110は、加速度ベクトル an カゝら抽出した鉛直成分 X (n)を帯域分割することにより得られた低域成分 xl (n)に基 づ 、てピーク位置の候補を検出し、さらにピーク位置の候補として検出した位置を含 む所定範囲における低域成分 xl (n)のエネルギー elと高域成分 xh(n)のエネルギー ehとに基づ 、て、ピーク位置の候補として信憑性の高 、ものだけをピーク位置の候 補として特定するようにして 、る。
[0146] なお、この図 17に示した処理は、制御部 110が実現するピーク検出 Z判定処理部 としての機能の内のピーク検出機能 (ピーク位置の候補の検出と特定機能)を実現す るものである。
[0147] [波形マッチングと歩数カウント処理にっ ヽて]
図 18は、図 15に示したこの実施の形態の歩数計 100の主に制御部 110において 実行される処理であって、ピーク位置の候補として特定されたものの中力 ピーク位 置を決定し、これに基づいて使用者の歩数を計数する処理を説明するためのフロー チャートである。
[0148] この図 18に示す処理は、図 17に示した処理により特定されたピーク位置の候補に 基づいて、ピーク位置の候補を含む所定区間同士の波形マッチングを行うことにより 、真のピーク位置を決定し、決定したピーク位置をカウントすることによって、使用者 の歩数を正確にカウントするものである。
[0149] 図 17に示した処理により特定されたピーク位置の候補は、図 17に示した処理のス テツプ S207の処理により、例えば、 RAMI 13などのピーク位置記録バッファに格納 される。このため、制御部 110は、ピーク位置記録バッファ力 未チェックのピーク位 置の候補を探索 (検索)する (ステップ S301)。未チ ック力否かは、後述もするように 、特定されたピーク位置の候補毎に設定されるチェック済みフラグのオン Zオフによ つて判別できる力 通常、最新に記録されたピーク位置の候補が未チェックのピーク 位置の候補となる。 [0150] そして、制御部 110は、ステップ S301の処理において、未チェックのピーク位置の 候補 plが見つ力つた力否かを判断する(ステップ S302)。ステップ S302の判断処理 において、未チェックのピーク位置の候補 piが見つ力つたと判断したときには、ピー ク位置の候補 p 1から過去 N秒以内の範囲にぉ 、て、ピーク位置の候補 p 1を含む所 定範囲の波形と、他のピーク位置の候補を含む所定範囲の波形とのマッチング処理 を行う(ステップ S 303)。
[0151] なお、ステップ S303の処理においては、サンプリング周波数にもよる力 過去の 1 つ以上の所定範囲の波形とのマッチングが行われることになる。また、マッチングに 用いられる波形は、低域成分 xl(n)の波形に基づいて行うようにすればよい。もちろ ん、比較する波形としては、鉛直成分 x (n)の波形を用いるようにしてもよい。
[0152] そして、制御部 110は、ステップ S303の処理に基づいて、未チェックのピーク位置 の候補 piとマッチする過去のピーク位置の候補 p2が見つ力 た力否かを判断する( ステップ S304)。すなわち、ステップ S304の判断処理においては、未チェックのピー ク位置の候補 piを含む所定範囲の波形と類似度の高い過去のピーク位置の候補を 含む所定範囲の波形が見つ力つた力否かを判断している。
[0153] ステップ S304の判断処理において、未チェックのピーク位置の候補 piとマッチす る過去のピーク位置の候補 p2が見つ力つたと判断したときには、未チェックのピーク 位置の候補 piをピーク位置として決定して歩数カウントに値「1」を加算すると共に、 未チェックのピーク位置の候補 piのカウント済みフラグをオンにする(ステップ S305)
[0154] そして、今回、マッチするものとして見つかった過去のピーク位置の候補 p2のカウ ント済みフラグはオンになっているか否かを判断する(ステップ S306)。ステップ S30 6の判断処理にお!、て、過去のピーク位置の候補 p2のカウント済みフラグがオンにな つて 、ると判断した場合には、未チェックのピーク位置の候補 piにつ 、てのチェック 済みフラグをオンにして (ステップ S307)、この図 18に示す処理を終了し、次の実行 タイミングを待つ。
[0155] また、ステップ S306の判断処理において、過去のピーク位置の候補 p2のカウント 済みフラグはオンになっていないと判断したときには、過去のピーク位置の候補 p2に ついて、それ以前にマッチするピークは存在しな力つたものの、新たにピーク位置の 候補 plとマッチしたので、過去のピーク位置の候補 p2につ 、てもピーク位置として 決定して歩数カウントに値「1」を加算すると共に、過去のピーク位置の候補 p2のカウ ント済みフラグをオンにして (ステップ S308)、ステップ S307の処理に進み、未チエツ クのピーク位置の候補 plについてのチェック済みフラグをオンにして (ステップ S307 )、この図 18に示す処理を終了し、次の実行タイミングを待つ。
[0156] また、ステップ S302の判断処理において、未チェックのピーク位置の候補 plが見 つからないと判断した場合と、ステップ S304の判断処理において、未チェックのピー ク位置の候補 plとマッチする過去のピーク位置の候補 P2が見つ力 ないと判断した 場合には、この図 18に示す処理を終了し、次の実行タイミングを待つことになる。
[0157] このように、図 17を用いて説明した処理により、ピーク位置の候補を特定した後に おいて、この図 18に示した処理により、ピーク位置を決定し、決定したピーク位置を カウントすることにより、使用者の歩行や走行に応じた歩数を正確にカウントすること ができる。
[0158] なお、ここでは、図 18の処理により、決定したピーク位置をカウントすることにより、 使用者の歩行や走行に応じた歩数を正確にカウントするようにしたが、 目的とする精 度が確保できれば、図 17の処理により、特定したピーク位置の候補をカウントするこ とにより、使用者の歩行や走行に応じた歩数をカウントするようにしてもよい。
[0159] なお、この図 18に示した処理は、制御部 110が実現するピーク検出 Z判定処理部 4としての機能の内のピーク位置の決定機能を実現するとともに、歩数計としての歩 数計測の機能を実現するものである。
[0160] [先頭のカウント漏れ防止処理について]
ところで、図 18に示した処理の場合、静止状態力も歩き始めた時など過去 N秒以 内にピークが無い場合についてはカウント漏れが生じてしまう。そこで、図 19のように 過去のピーク力 現在のピークに対してマッチングを取る処理を追加することによつ てカウント漏れを防ぐことができる。
[0161] すなわち、使用者の動作 (運動)開始直後の歩数のカウント漏れを防止する場合に は、図 19に示す処理を、図 18に示した処理前段または後段に追加する。この場合、 制御部 110は、ピーク位置記録バッファ力 現在よりも N秒前のピーク位置の候補で あって、未カウントのピーク位置の候補 piを探索 (検索)する(ステップ S401)。
[0162] このステップ S401の処理は、過去の未カウントのピーク位置の候補を見つける処 理であり、端的に言えば、先頭のピーク位置の候補を検索する処理である。そして、 制御部 110は、未カウントのピーク位置の候補 piが見つ力つた力否かを判断する (ス テツプ S402)。
[0163] ステップ S402の判断処理において、未カウントのピーク位置の候補 piが見つかつ たと判断したときには、未カウントのピーク位置の候補 piから現在にかけての範囲を 対象として、未カウントのピーク位置の候補 piを含む所定範囲の波形と、未カウント のピーク位置の候補 p 1から現在にかけての範囲に存在するピーク位置の候補を含 む所定範囲の波形とのマッチングを行う(ステップ S403)。
[0164] 図 18に示したステップ S303の処理が現在から過去方向に向力う波形マッチング 処理であるのに対して、ステップ S403の処理は、過去から現在に向かう波形マッチ ング処理である。
[0165] そして、制御部 110は、ステップ S403の処理に基づいて、未チェックのピーク位置 の候補 piとマッチする現在方向のピーク位置の候補 p2が見つ力 た力否かを判断 する(ステップ S404)。すなわち、ステップ S404の判断処理においては、未カウント となって!/、るピーク位置の候補 piを含む所定範囲の波形と、ピーク位置の候補 piか ら見て現在方向(時間の経過方向)に類似度の高いピーク位置の候補を含む所定範 囲の波形が見つかった力否かを判断して 、る。
[0166] ステップ S404の判断処理において、未チェックのピーク位置の候補 piとマッチす る現在方向のピーク位置の候補 p2が見つ力つたと判断したときには、歩数カウントに 値「 1」を加算すると共に、未チェックのピーク位置の候補 p 1のカウント済みフラグをォ ンにする(ステップ S405)。
[0167] これにより、過去にピーク位置の候補が存在しないために、ピーク位置として決定で きな力つた過去のピーク位置の候補を、それよりも後の波形と比較することにより類似 度の高 、波形が存在した場合に、過去のピーク位置の候補をピーク位置として決定 して、これをカウントすることができる。 [0168] さらに、制御部 110は、現在方向のピーク位置の候補 p2についてのカウント済みフ ラグがオンか否かを判断する(ステップ S406)。ステップ S406の判断処理において、 現在方向のピーク位置の候補 p2についてのカウント済みフラグがオンでないと判断 したときには、その現在方向のピーク位置の候補 p2についてもピーク位置として決定 して、歩数カウントに値「1」を加算すると共に、現在方向のピーク位置の候補 p2につ V、てのカウント済みフラグをオンにする(ステップ S407)。
[0169] そして、ステップ S406の判断処理において、現在方向のピーク位置の候補 p2につ いてのカウント済みフラグがオンであると判断した場合、または、ステップ S407の処 理の後においては、制御部 110はこの図 19に示す処理を終了し、次の実行タイミン グを待つことになる。
[0170] このような図 19に示した処理を、図 18に示した処理にカ卩えて実行することにより、 例えば、一番初めのピーク近傍の波形については、時間的に後(未来)のピーク近傍 の波形とマッチングを取ることにより、真にピーク力否かを判断し、ピークであるときに はこれをカウントすることができる。
[0171] [体動種別の推定処理と歩数カウント処理につ!、て]
上述したように、ピーク位置の候補を特定して、これをカウントしたり、あるいは、ピー ク位置の候補の中力 真のピーク位置を決定して、これをカウントしたりすることにより 、使用者の歩行や走行などの運動時における歩数を計測することが可能である。し かし、ノイズなどの影響で、マッチングが取れないピークが生ずることがある。
[0172] このため、この実施の形態の歩数計 100においては、制御部 110が実現するステツ プ位置解析部 5としての機能によって、使用者の動作状態 (動作状態に応じた体動 種別)をも正確に把握し、歩行や走行を行っている場合には、マッチングが取れない ピークが存在する期間においても、使用者の歩数を正確にカウントすることができるよ うにしている。
[0173] 図 20、図 21は、この実施の形態の歩数計 100の主に制御部 110において実行さ れる体動種別の推定処理と歩数カウント処理を説明するためのフローチャートである 。この図 20、図 21に示す処理において用いられるピーク位置の情報としては、図 17 の処理により特定されピーク位置記録バッファに記録されたピーク位置の候補のデ ータを用い、波形マッチングによるピーク位置の決定処理と、体動種別の推定処理と 、歩数カウント処理とを行うようにしているものである。
[0174] すなわち、図 16、図 17の処理に続いて、この図 20、図 21に示す処理が実行される 。この図 20、図 21の処理は、ピーク検出 Z判定処理部 4の機能と、ステップ位置解 析部 5の機能とにより実現されるものである。
[0175] そして、制御部 110は、 RAMI 13などのピーク位置記録バッファ力 未チェックの ピーク位置の候補を探索 (検索)する (ステップ S501)。未チ ック力否かは、後述も するように、特定されたピーク位置の候補毎に設定されるチヱック済みフラグのオン Z オフによって判別できる力 通常、最新に記録されたピーク位置の候補が未チェック のピーク位置の候補となる。
[0176] そして、制御部 110は、ステップ S501の処理において、未チェックのピーク位置の 候補 plが見つ力つた力否かを判断する(ステップ S502)。ステップ S502の判断処理 において、未チェックのピーク位置の候補 piが見つ力つたと判断したときには、静止 タイマ Cに値「0 (ゼロ)」をセットする (ステップ S503)。そして、制御部 110は、ピーク 位置の候補 p 1から過去 N秒以内の範囲にぉ 、て、ピーク位置の候補 p 1を含む所定 範囲の波形と、他のピーク位置の候補を含む所定範囲の波形とのマッチング処理を 行う(ステップ S 504)。
[0177] なお、ステップ S504の処理においては、サンプリング周波数にもよる力 過去の 1 つ以上の所定範囲の波形とのマッチングが行われることになる。また、マッチングに 用いられる波形は、低域成分 xl(n)の波形に基づいて行うようにすればよい。もちろ ん、比較する波形としては、鉛直成分 x (n)の波形を用いるようにしてもよい。
[0178] そして、制御部 110は、ステップ S504の処理に基づいて、未チェックのピーク位置 の候補 piとマッチする過去のピーク位置の候補 p2が見つ力 た力否かを判断する( ステップ S505)。すなわち、ステップ S505の判断処理においては、未チェックのピー ク位置の候補 piを含む所定範囲の波形と類似度の高い過去のピーク位置の候補を 含む所定範囲の波形が見つ力つた力否かを判断している。
[0179] ステップ S505の判断処理において、未チェックのピーク位置の候補 piとマッチす る過去のピーク位置の候補 p2が見つ力つたと判断したときには、未チェックのピーク 位置の候補 plをピーク位置として決定して歩数カウントに値「1」を加算すると共に、 未チェックのピーク位置の候補 piのカウント済みフラグをオンにする(ステップ S506)
[0180] そして、今回、マッチするものとして見つかった過去のピーク位置の候補 p2のカウ ント済みフラグはオンになっているか否かを判断する(ステップ S507)。ステップ S50 7の判断処理にお!、て、過去のピーク位置の候補 p2のカウント済みフラグがオンにな つていると判断した場合には、図 21に示す処理に進む。
[0181] また、ステップ S507の判断処理において、過去のピーク位置の候補 p2のカウント 済みフラグはオンになっていないと判断したときには、過去のピーク位置の候補 p2に ついて、それ以前にマッチするピークは存在しな力つたものの、新たにピーク位置の 候補 plとマッチしたので、過去のピーク位置の候補 p2につ 、てもピーク位置として 決定して歩数カウントに値「1」を加算すると共に、過去のピーク位置の候補 p2のカウ ント済みフラグをオンにして (ステップ S508)、図 21の処理に進む。
[0182] また、図 20に示したステップ S502の判断処理において、未チェックのピーク位置 の候補 P 1が見つからな 、と判断したときには、静止タイマ Cに値「 1」を加算し (ステツ プ S509)、静止タイマ Cが基準値よりも大きくなつた力否かを判断する (ステップ S51 0)。
[0183] ステップ S510の判断処理において、静止タイマ Cが基準値よりも大きくなつたと判 断したときには、使用者の動作状態は「静止」状態であると設定し (ステップ S511)、 この図 20、図 21の処理を終了し、次の実行タイミングを待つことになる。また、ステツ プ S510の判断処理にぉ 、て、静止タイマ Cが基準値よりも大きくなつて 、な 、と判断 したときには、何もせずに、この図 20、図 21に示す処理を終了し、次の実行タイミン グを待つことになる。
[0184] そして、上述したように、ステップ S507の判断処理にお!、て、過去のピーク位置の 候補 p2のカウント済みフラグがオンになっていると判断した場合と、ステップ S508の 処理の後においては、図 21に示す処理を行うことになる。図 21に示す処理において は、制御部 110は、まず、未チェックのピーク位置の候補 plについてのチェック済み フラグをオンにする (ステップ S 512)。 [0185] 次に、制御部 110は、ステップ S512の処理により、最新にチェック済みフラグがォ ンになつたピーク位置 p 1に最も近 、過去のピーク位置 p 1 'を特定し (ステップ S 513) 、ピーク位置 piとその直近のピーク位置 ρΓとの間隔 TOを算出する (ステップ S514) 。この実施の形態の歩数計 100においては、図 12を用いて説明したように、直近の 4 つのピーク位置のそれぞれによって挟まれた 3つのピーク間隔 TO、 Tl、 Τ2を保持す ることができるようにしている。もちろん、より多くのピーク間隔を保持して用いるように してちよい。
[0186] そして、制御部 110は、最新の使用者の動作状態の判定結果が [歩行 Ζ走行]状 態か否かを判断する (ステップ S515)。ステップ S515の判断処理において、使用者 の動作状態力 S [歩行 Ζ走行]状態ではないと判断したときには、上述したように、直近 の 4つのピーク位置に応じて決まる 3つの間隔 ΤΟ、 Tl、 Τ2の平均値 Taを算出し (ス テツプ S516)、全てのピーク間隔 TO、 Tl、 Τ2のそれぞれについて、平均値 Taを減 算した値を求め、この求めた値が全てのピーク間隔について、基準値よりも小さいか 否かを判断する (ステップ S517)。
[0187] このステップ S517の判断処理は、ピーク間隔 TO、 Tl、 Τ2のそれぞれがほぼ同じ ピッチであり、周期的な体動である歩行や走行が行われるようになった力否かを判断 する処理である。ステップ S517の判断処理において、ピーク間隔 ΤΟ、 Tl、 Τ2のそ れぞれと平均値 Taとの差分のそれぞれが、基準値以下であると判断したときには、 制御部 110は、平均値 Taを基準ピッチ Psとするとともに、使用者の動作状態を「歩行 Z走行」状態にあると判定 (特定)する (ステップ S 518)。
[0188] そして、制御部 110は、ピーク間隔 TO、 Tl、 Τ2のシフト処理を行う(ステップ S519 ) οすなわち、ピーク間隔 Τ2にピーク間隔 T1を移行し、ピーク間隔 T1にピーク間隔 Τ 0を移行する。なお、ピーク間隔 TOについては、初期化を行うようにしてもよい。
[0189] このステップ S519の処理の後、および、ステップ S517の判断処理において、ステ ップ S517の判断処理において、ピーク間隔 TO、 Tl、 Τ2のそれぞれと平均値 Taと の差分のそれぞれが、基準値以下ではないと判断したときには、制御部 110は、この 図 20、図 21に示す処理を終了し、次の実行タイミングを待つことになる。
[0190] また、ステップ S515の判断処理において、使用者の動作状態が [歩行 Z走行]状 態であると判断したときには、制御部 110は、最新のピーク間隔 TOを基準ピッチ Psで 割り算することにより値 rを算出する(ステップ S520)。このステップ S520の処理は、 図 14を用いて説明したように、ノイズ等の影響により、本来、ピーク値として検出され るべきものがピーク値として検出されない場合を考慮する処理であり、最新に検出し たピーク間隔 TO力 基準ピッチ Psの何倍になるかを計算する処理である。
[0191] そして、ステップ S520において算出した値 rに最も近い整数 Rを算出する (ステップ S521)。例えば、値 rが「0. 1」や「0. 2」などである場合には、整数 Rは「0」であり、値 rが「0. 9」や「1. 1」などである場合には、整数 Rは「1」であり、また、値 rが「1. 9」や「 2. 1」であれば、整数 Rは「2」になる。そして、制御部 110は、値 rから整数 Rを減算 することにより得られる値の絶対値力 予め決められる基準値よりも小さいか否かを判 断する (ステップ S522)。
[0192] このステップ S522の判断処理は、最新のピーク間隔 TO力 基準ピッチ Psの整数 倍か否かを判断する処理である。ステップ S522の判断処理において、値 rから整数 Rを減算することにより得られる値の絶対値が、予め決められる基準値よりも小さくな いと判断したときには、最新のピーク間隔 TOが、基準ピッチ Psの整数倍ではないと 判断し、使用者の動作状態を「不定」状態にあると判定 (特定)する (ステップ S523) 。そして制御部 110は、この図 20、図 21に示す処理を終了し、次の実行タイミングを 待つことになる。
[0193] また、ステップ S522の判断処理において、値 r力も整数 Rを減算することにより得ら れる値の絶対値が、予め決められる基準値よりも小さいと判断したときには、最新のピ ーク間隔 TOが、基準ピッチ Psの整数倍であると判断し、整数 R1は、値「1」より大きい か否かを判断する(ステップ S 524)。このステップ S524の判断処理は、ピーク間隔 T 0が基準ピッチ Psの 2倍以上力否かを判断する処理である。
[0194] ステップ S524の判断処理において、整数 Rは 1より大きくないと判断したときには、 最新のピーク間隔 TOは、基準ピッチの 2倍以上の区間ではないので、制御部 110は 、何もすることなぐこの図 20、図 21に示す処理を終了し、次の実行タイミングを待つ ことになる。
[0195] 一方、ステップ S524の判断処理において、整数 Rは 1より大きいと判断したときに は、最新のピーク間隔 TOは、基準ピッチの 2倍以上の区間であるので、制御部 110 は、カウント数に、整数 Rから 1減算した値を加算して (ステップ S525)、この図 20、図 21に示す処理を終了し、次の実行タイミングを待つことになる。
[0196] このように、図 20、図 21に示した処理により、使用者の動作状態をも正確に判別し つつ、歩行や走行を行っている場合には、マッチングが取れないピークが存在する 期間においても、使用者の歩数を正確にカウントすることができ、しかも、「静止」状態 、「不定」状態、「歩行」状態、「走行」状態の 4つの動作状態 (動作状態に応じた体動 種別)をリアルタイムに判別することもできるようにして ヽる。
[0197] なお、図 20、図 21を用いて説明した処理において、最新のピーク間隔 TOが、基準 ピッチの整数倍の区間であるときには、計測できな力つたピーク位置に応じて、歩数 の計測を行うようにした力 これに限るものではない。すなわち、最新のピーク間隔 T 0が、基準ピッチの整数倍の区間であるときには、その区間における鉛直成分の低域 成分の振幅、あるいは、鉛直成分自体の振幅をも考慮し、基準ピッチの整数倍の位 置に確実にピーク位置が存在するカゝ否かを判別するようにしてもょ 、。
[0198] 具体的には、基準ピッチの整数倍の位置の振幅に対する基準値を設けておき、こ の基準値以上であれば、ピーク値及びピーク位置として認識し、基準値より小さい場 合には、ピーク値及びピーク位置として認識しないようにすればよい。なお、基準値 は、実験を行うことにより、適切な値を説定することが可能である。
[0199] [音響再生装置への適用について]
次に、この発明の装置、方法、プログラムの一実施の形態を音響再生装置に適用 した場合について説明する。図 22は、この実施の形態の音響再生装置 200を説明 するためのブロック図である。この実施の形態の音響再生装置 200は、後述もするよ うに、楽曲データ (音楽データ)を記憶するための比較的に大容量の記録媒体を用 V、ることができるようにされた携帯型のものである。
[0200] 記録媒体としては、ハードディスク、 MD (ミニディスク (登録商標) )などの光磁気デ イスク、 CDや DVDなどの光ディスク、メモリカード、半導体メモリなど、種々の記憶媒 体を用いるものが考えられる。しかし、ここでは、説明を簡単にするため、楽曲データ などのコンテンツデータを記憶する記録媒体はハードディスクであるものとして説明 する。
[0201] 図 22に示すように、この実施の形態の音響再生装置 200は、制御部 210に対して 、 3軸の加速度センサ 201が AZD変換部 202を介して接続されると共に、比較的に 記憶容量の大きな記録媒体として、楽曲データベース(以下、楽曲 DBという。) 203 と、プレイリスト格納部 204とが接続される。
[0202] また、制御部 210に対しては、音声再生処理部 205を介してスピーカ 206が接続さ れると共に、ユーザインターフェースとしての操作部 221と、表示部 222とが接続され 、さらに外部インターフェース(以下、外部 IZFという。) 223を介して入出力端子 224 が接続された構成とされて 、る。
[0203] 3軸の加速度センサ 201は、上述した歩数計 100の加速度センサ 101の場合と同 様に、 1軸の加速度センサを互いに直交する X軸、 Y軸、 Z軸の 3軸に配置して形成さ れたもの、あるいは、直交する 3軸の加速度センサがひとつのパッケージに封入され て形成されたもののいずれを用いてもよい。 3軸の加速度センサ 201からの X軸、 Y 軸、 Z軸のそれぞれの検出出力(アナログ出力)は、 AZD変換部 202に供給され、こ こで制御部 210において処理可能な形式のデジタルデータに変換されて制御部 21 0に供給される。
[0204] 制御部 210は、この実施の形態の音響再生装置 200の各部を制御するものであり 、上述した歩数計 100の帘 IJ御咅 10と同様に、 CPU211、 ROM212, RAM213, 不揮発性メモリ 214力 CPUバス 215を通じて接続されてマイクロコンピュータの構 成とされたものである。
[0205] ここで、 CPU211は、各種のプログラムを実行し、各部に供給する制御信号を形成 したり、各種の演算を行ったりするなど、制御部 210における処理や制御の主体とな るものである。 ROM212は、 CPU211が実行する各種のプログラムや処理に必要に なるデータを記憶保持する。
[0206] また、 RAM213は、処理の途中結果を一時記憶するなど、主に作業領域として用 いられるものであり、 AZD変換部 202を通じて供給される加速度センサ 201からの 数値ィ匕された加速度データを記憶保持するバッファメモリなどとしても用いられる。不 揮発性メモリ 214は、例えば、 EEPROMやフラッシュメモリなどの電源が落とされても 記憶データが消滅することのないメモリであり、電源が落とされても保持しておくべき データ、例えば設定されたパラメータや追加されたプログラムなどを記憶保持する。
[0207] そして、楽曲 DB203、プレイリスト格納部 204は、上述もしたように、いずれもハード ディスクである。楽曲 DBには、この実施の形態の音響再生装置 200において再生可 能な多数の楽曲データ力 例えば所定のデータ圧縮方式に従ってデータ圧縮され て記憶保持されている。
[0208] また、プレイリスト格納部 204には、使用者が「静止」状態にあるときに再生すべき楽 曲を指示する静止時プレイリストと、使用者が「歩行」状態にあるときに再生すべき楽 曲を指示する歩行時プレイリストと、使用者が「走行」状態にあるときに再生すべき楽 曲を指示する走行時プレイリストとが記憶保持されている。
[0209] すなわち、プレイリスト格納部 204に格納されているプレイリストは、静止時プレイリ スト、歩行時プレイリスト、走行時プレイリストの別力 例えば、各プレイリストに付加さ れるプレイリスト種別等の情報によって区別することができるようにされている。また、 静止時プレイリスト、歩行時プレイリスト、走行時プレイリストのそれぞれは複数存在し 、各プレイリストには、再生すべき 1曲以上の楽曲を特定するファイル名などの楽曲 I D (楽曲識別子)が例えば再生順に予め登録されて形成されたものである。
[0210] なお、図 22においては、構成を明確にするために、楽曲 DB203と、プレイリスト格 納部 204とを別個に示した力 これらは同じハードディスク上に形成される場合ももち ろんある。
[0211] 音声再生処理部 205は、制御部 210により楽曲 DB203から読み出された楽曲デ ータの供給を受けて、これを圧縮解凍処理してデータ圧縮前の元の楽曲データを復 元すると共に、復元した楽曲データからスピーカ 206に供給する形式のアナログ音声 信号を形成し、これをスピーカ 206に供給する。これにより、再生するようにされた楽 曲データに応じた音声がスピーカ 206から放音される。
[0212] なお、図示しないが、この実施の形態の音声再生処理部 205には、ヘッドホン端子 も接続するようにされており、当該ヘッドホン端子にヘッドホンが接続された場合には 、音声再生処理部 205において処理された音声信号は、スピーカ 206にではなくへ ッドホン端子を介して、これに接続されたヘッドホンに供給され、当該ヘッドホンを通 じて再生音声を聴取することができるようにして 、る。
[0213] また、操作部 221は、再生キー、停止キー、早送りキー、早戻しキー、その他、各種 のファンクションキーなどを備え、ユーザー力 の操作入力を受け付けて、これに応じ た電気信号を制御部 210に供給することができるものである。これに応じて、制御部 2 10は各部を制御し、ユーザーの指示に応じた処理を行うことができるようにしている。
[0214] また、表示部 222は、表示制御回路を含み、例えば、 LCD (Liquid Crystal Display 有機 EL (Electro Luminescence)ディスプレイ、 CRT (Cathode- Ray Tube)等の表 示素子を備え、制御部 210の制御に応じて種々のガイダンス情報等を表示する。具 体的には、表示部 222は、制御部 210から表示データの供給を受けて、この表示デ ータから表示素子に供給する映像信号を形成し、これを表示素子に供給することに よって、制御部 210からの表示データに応じた表示情報を表示素子の表示画面に表 示する。
[0215] また、外部 I/F223、入出力端子 224を通じて、パーソナルコンピュータなどの外 部機器に接続することができるようにされている。そして、入出力端子 224、外部 IZF 223を通じて、外部機器から楽曲データやプレイリストの供給を受け、楽曲データに ついては楽曲 DB203に格納し、プレイリストについては、プレイリスト格納部 204に 格納することができるようにして 、る。
[0216] 逆に、制御部 210の制御により、楽曲 DB203に格納されている楽曲データを制御 部 210、外部 IZF223、入出力端子 224を通じて外部機器に出力し、外部機器の記 録媒体にバックアップを形成したり、同様にして、プレイリスト格納部 204に記憶保持 されているプレイリストのノ ックアップを外部機器の記録媒体に作成したりすることもで きるようにしている。
[0217] また、表示部 222に楽曲 DB203に格納されていえる楽曲データの一覧リストを表 示させ、操作部 221を通じて目的する楽曲データを選択すると共に、静止時プレイリ スト、歩行時プレイリスト、走行時プレイリストの別を入力することにより、この実施の形 態の音響再生装置 200においてもプレイリストを作成することができる。
[0218] そして、この実施の形態の音響再生装置 200は、操作部 221を通じて選択された 楽曲データが制御部 210により楽曲 DB203から読み出され、これが音声再生処理 部 205に供給される。これにより、使用者から操作部 221を通じて指示された楽曲を 再生し、その楽曲を聴取することができるようにされる。
[0219] また、操作部を通じて用いるプレイリストを指示した場合は、その指示したプレイリス に従って、制御部 210により楽曲 DB203から楽曲データが読み出され、音声再生処 理部 205に供給するようにされ、プレイリストに従って楽曲を再生することもできるよう にしている。
[0220] さらに、プレイリストを自動選択するモードが選ばれている場合には、使用者の動作 状態に応じて適切なプレイリストを選択することもできるようにしている。この実施の形 態の音響再生装置 200において、制御部 210が、図 1に示す体動検出装置におけ る鉛直成分抽出部 2、高域 Z低域分解部 3、ピーク検出 Z判定処理部 4、ステップ位 置解析部 5としての機能を実現すると共に、使用者の動作ピッチ (動作テンポ)を正確 に把握し、使用者の動作状態に応じて、静止時プレイリスト、歩行時プレイリスト、走 行時プレイリストの中力 適切なものを自動選択することができるようにして 、る。
[0221] すなわち、この実施の形態の音響再生装置 200の制御部 210は、加速度センサ 2 01からの検出出力に基づいて、図 1を用いて説明したように、鉛直成分抽出、鉛直 成分の高域 Z低域の分解、ピーク位置の候補の検出とピーク位置の候補の特定、波 形マッチングによるピーク位置の決定などの処理を行う機能を有している。
[0222] そして、決定したピーク位置に基づいて、使用者の動作状態を正確に把握すると共 に、歩行や走行などの体動ピッチを正確に把握し、体動ピッチをも考慮して、使用者 の動作状態に応じて適切なプレイリストを自動的に選択し、この選択したプレイリスト を用いて楽曲の再生を行うことができるようにして 、る。
[0223] 以下、図 23のフローチャートを参照しながら、この実施の形態の音響再生装置にお いて、プレイリストを自動選択するモードが選ばれている場合において、使用者の動 作状態に応じて適切なプレイリストを選択する場合の処理について説明する。図 23 は、この実施の形態の音響再生装置において、使用者の動作状態に応じて適切な プレイリストを選択する場合の処理について説明するためのフローチャートである。
[0224] 図 23に示す処理は、主に、この実施の形態の音響再生装置 200の制御部 210に おいて行われる。まず、制御部 210は、図 15に示した歩数計 100の制御部 110の場 合と同様に、体動種別の推定処理を実行する (ステップ S601)。このステップ S601 の処理は、(1)図 16を用いて説明したように、 3軸の加速度センサ 201からの加速度 ベクトルから鉛直成分を抽出すると共に、 (2)この抽出した鉛直成分を高域成分と低 域成分とに分離し、(3)これらを用いて、図 17を用いて説明したように、ピーク位置の 候補を特定して、(4)この特定したピーク位置の候補の情報を用いて、図 20、図 21 に示した処理を実行することにより、使用者の動作状態と動作ピッチを把握する処理 である。
[0225] なお、図 20、図 21に示した処理においては、歩数をカウントする処理を行う必要は ない。ステップ S506〜ステップ S508の処理、および、ステップ S525の処理を行う必 要はない。
[0226] そして、制御部 210は、使用者の動作状態が変化した力否かを判断する (ステップ S602)。この判断処理は、例えば、動作状態の判定結果として、前回の判定結果と 今回の判定結果とを保持しておき、両者が異なる場合に動作状態が変化したと判断 することができる。
[0227] ステップ S602の判断処理にぉ 、て、動作状態は変化して!/、な 、と判断したときに は、この図 23に示す処理を終了し、次の実行タイミングまで待つことになる。ステップ S602の判断処理において、動作状態が変化したと判断したときには、変化後の状 態、すなわち、今回の判定結果が、「静止」状態である力否かを判断する (ステップ S 603)。
[0228] ステップ S603の判断処理において、今回の判定結果が「静止」状態であると判断 したときには、プレイリスト格納部 204に格納されている静止時プレイリストを用いるよ うにし、当該静止時プレイリストに従って楽曲(音楽)を再生する (ステップ S604)。こ の後、この図 23に示す処理を終了し、次の実行タイミングまで待つことになる。
[0229] また、ステップ S603の判断処理にぉ 、て、今回の判定結果が「静止」状態ではな いと判断したときには、今回の判定結果力 ^歩行 Z走行」状態か否かを判断する (ス テツプ S605)。ステップ S605の判断処理において、「歩行 Z走行」状態であると判 断したときには、ステップ S601の処理であって、図 20、図 21の処理において求めら れた基準ピッチ Psが規定値よりも小さ 、か否かを判断する (ステップ S606)。 [0230] このステップ S606の判断処理において用いられる規定値は、使用者の動作状態 1S 歩行か走行かを判断するためのものであり、基準ピッチ Psが規定値よりも小さくな い場合には、基準ピッチはまだ十分に早くない状態であり、歩行状態であると判断す ることができる。逆に、基準ピッチが規定値よりも小さい場合には、基準ピッチは十分 に早く走行状態であると判断することができる。
[0231] このため、ステップ S606の判断処理において、基準ピッチ PSが規定値よりも小さく ないと判断したときには、使用者の動作状態は、「歩行」状態であると判断し、プレイリ スト格納部 204に格納されて 、る歩行時プレイリストを用いるようにし、当該歩行時プ レイリストに従って楽曲(音楽)を再生する (ステップ S607)。この後、制御部 210は、 図 23に示す処理を終了し、次の実行タイミングまで待つことになる。
[0232] また、ステップ S606の判断処理にぉ 、て、基準ピッチ PSが規定値よりも小さ 、と判 断したときには、使用者の動作状態は、「走行」状態であると判断し、プレイリスト格納 部 204に格納されている走行時プレイリストを用いるようにし、当該走行時プレイリスト に従って楽曲(音楽)を再生する (ステップ S608)。この後、制御部 210は、図 23に 示す処理を終了し、次の実行タイミングまで待つことになる。
[0233] このように、この実施の形態の音響再生装置 200は、使用者の動作状態を適切に 判定し、使用者の動作状態に応じたプレイリストに自動的に切り換えて、使用者の体 動に応じた楽曲を再生することができるようにされる。
[0234] なお、上述もしたように、静止時プレイリスト、歩行時プレイリスト、走行時プレイリスト のそれぞれが複数用意されている場合もある。このような場合には、使用する順番を 予め定めておき、その順番に従って使用するようにしたり、各プレイリスト毎に使用頻 度を記憶しておき、使用頻度の低いものを用いるようにしたり、逆に、使用頻度の高 いものを用いるようにしたりしてもよい。また、静止時プレイリスト、歩行時プレイリスト、 走行時プレイリストのそれぞれにつ!/、て用いるものを予め使用者が指定するようにし てもよい。
[0235] また、上述した歩数計 100においても、図 20、図 21の処理に先立って、図 19の処 理を行うことにより、歩数の計測処理の初めから、正確の歩数の計測を行うようにする ことができる。 [0236] また、上述した歩数計 100、音響再生装置 200においては、 3軸の加速度センサ 1 01、 201を用いるものとして説明した。しかし、これに限るものではない。 1軸の加速 度センサを用いてもよいし、 2軸の加速度センサを用いてもよい。しかし、 1軸または 2 軸の加速度センサを用いる場合は、体動の主要成分が信号として拾えるような向きに 配置される必要がある。
[0237] しかし、 3軸の加速度センサを用い、上述もしたように、鉛直成分の抽出、高域成分 と低域成分のエネルギー比の考慮といった処理を行うことによって、加速度センサの 装着位置や装着方向に制限を生じさせることが無いので、 3軸の加速度センサを用 、ることが好まし!/、。
[0238] また、上述した実施の形態においては、 3軸の加速度センサを用い、 3軸の加速度 ベクトル力も鉛直成分を抽出するようにした力 これに限るものではない。例えば、 3 軸の内で最も信号強度の強!ヽ軸の信号を鉛直成分とみなして選択するようにしても よい。
[0239] また、上述した歩数計 100、音響再生装置 200にお 、ては、加速度ベクトルの鉛直 成分のピーク位置を検出していくことにより、体動のピッチ周期を求めるようにした力 これに限るものではな 、。ピーク位置 (ステップ位置)を特定するための手段としては 、靴底にスィッチを配置し歩行または走行時足を踏み込むことでスィッチの接点が閉 じるような機構を用いたり、靴底にひずみゲージを配置したりして踏み込みを検出す る手段を用いてもよい。
[0240] また、上述した実施の形態の音響再生装置 200においては、使用者の体動ピッチ ( 歩行テンポや走行テンポ)を正確に計測することができるので、制御部 210が音声再 生処理部 205を制御することにより、使用者の体動テンポに応じて、再生する楽曲の 再生テンポを制御するようにすることもできる。
[0241] このように、使用者の体動テンポに応じて、再生する楽曲の再生テンポを制御する ことにより、使用者は再生される楽曲に対して一体感を感じ、歩行や走行などの運動 を効率よく行ったり、無理なく継続したりすることができるようにされる。
[0242] また、上述した実施の形態の体動検出装置、体動検出方法、体動検出プログラム は、歩数計、携帯型の音響再生装置だけでなぐフィットネスクラブ等で用いられる室 内用のランニングマシンなどのトレーニングマシンに適用したり、据え置き型の音響 再生装置に適当したりするなど、使用者の鉛直方向の体動を検出して利用する必要 性のある種々の電子機器に適用することができる。
[0243] また、図 16〜図 21、図 23は、この発明による、方法、プログラムが適用されたもの であり、図 16〜図 21、図 23に示したフローチャートに従って処理を行うことにより、こ の発明による方法を用いることができるし、また、図 16〜図 21、図 23に示したフロー チャートに従ったプログラムを作成することにより、この発明によるプログラムを実現す ることがでさる。
[0244] また、図 1に示した各部の機能を実現するように制御部での処理を規定することに よりこの発明による方法を実現でき、同様に、図 1に示した各部の機能を実現するよう に制御部で実行するプログラムを作成することにより、この発明によるプログラムを実 現することができる。
[0245] この発明によれば、ノイズを除去し使用者の鉛直方向の体動を正確に検出すること ができる。すなわち、ノイズに強い(ノイズの影響を受け難い)体動を検出するための 装置、方法、プログラムを実現できる。

Claims

請求の範囲
[1] 使用者の身体に装着するようにされる加速度センサと、
前記加速度センサ力ゝらの検出出力から加速度の鉛直成分を抽出する鉛直成分抽 出手段と、
前記鉛直成分抽出手段によって抽出された前記鉛直成分を高域成分と低域成分 とに成分分解する分解手段と、
前記分解手段によって分解された前記鉛直成分の前記低域成分に基づいてピー ク位置の候補を検出する検出手段と、
前記検出手段によって検出されたピーク位置の候補を含む所定範囲における前記 低域成分のエネルギーと前記高域成分のエネルギーとの比が所定値よりも小さい場 合に前記ピーク位置の候補をピーク位置の候補として特定する特定手段と、 前記特定手段によって特定されるピーク位置の候補に基づ 、て、前記使用者の体 動を検出する体動検出手段と
を備えることを特徴とする体動検出装置。
[2] 請求項 1に記載の体動検出装置であって、
前記特定手段によって特定された各ピーク位置の候補毎に、ピーク位置の候補を 含む所定範囲を設定し、設定された所定範囲毎の波形を対象として他の所定範囲 の波形とのマッチングを行い、類似度が高い場合に、対象となった波形に含まれるピ ーク位置の候補をピーク位置として決定する決定手段を備え、
前記体動検出手段は、前記決定手段によって決定されたピーク位置に基づいて、 使用者の体動を検出することを特徴とする体動検出装置。
[3] 請求項 1または請求項 2に記載の体動検出装置であって、
前記特定手段によって特定された複数のピーク位置の候補力 なる時系列パター ンについて、あるいは、前記決定手段によって決定された複数のピーク位置力 なる 時系列パターンについて解析することにより、ステップ間隔を推定する間隔推定手段 と、
前記間隔推定手段において推定された前記ステップ間隔に基づいて、少なくとも「 静止」、「歩行 Z走行」、「不定」の 3状態を使用者の動作状態として判別する判別手 段と
を備えることを特徴とする体動検出装置。
[4] 請求項 3に記載の体動検出装置であって、
判定手段は、ピーク位置が全く検出されない「静止」状態でもなぐステップ間隔が 等間隔でもない「不定」の状態にあるときに、前記間隔推定手段によって、数ステップ の間、ステップ間隔が等間隔であることが推定された場合に、使用者の動作状態を「 不定」から「歩行 Z走行」に変更することができるものであり、
前記判定手段によって、使用者の動作状態が「不定」から「歩行 Z走行」に変更さ れた場合に、前記間隔推定手段によって推定されたステップ間隔に基づいて、基準 ピッチを算出する基準ピッチ算出手段と
を備えることを特徴とする体動検出装置。
[5] 請求項 4に記載の体動検出装置であって、
前記判定手段は、前記間隔推定手段により推定された最新のステップ間隔と、前 記基準ピッチとを比較し、その差が、予め設定された誤差範囲に収まっているか、基 準ピッチの整数倍に対して所定の誤差範囲に収まっている場合には、「歩行 Z走行」 の状態を維持することを特徴とする体動検出装置。
[6] 請求項 4に記載の体動検出装置であって、
前記間隔推定手段によって推定されたステップ間隔に基づ!/、て、使用者のステツ プをカウントするカウント手段を備え、
前記判別手段によって使用者の動作状態が「歩行 Z走行」であると判別されて 、る 場合に、前記カウント手段は、前記ステップ間隔の前記基準ピッチに対する倍数に基 づいて、カウントの補間を行うことを特徴とする体動検出装置。
[7] 請求項 4に記載の体動検出装置であって、
前記判別手段により、使用者の動作状態が「歩行 Z走行」の状態にあると判別され ている場合において、前記間隔推定手段によって推定された前記ステップ間隔に基 づ ヽて、使用者の動作状態が「歩行」か「走行」かを判定する判定手段を備えたことを 特徴とする体動検出装置。
[8] 請求項 1に記載の体動検出装置であって、 前記判別手段により、使用者の動作状態が「歩行 Z走行」の状態にあると判別され て 、る場合にぉ 、て、前記取得手段によって取得された前記低域成分のピーク値と 、前記間隔推定手段によって推定された前記ステップ間隔とに基づいて、使用者の 動作状態が「歩行」か「走行」かを判定する判定手段を備えたことを特徴とする体動検 出装置。
[9] 請求項 3に記載の体動検出装置であって、
使用者の動作状態に対応して形成された 1以上の楽曲再生リストを記憶保持する 楽曲再生リスト保持手段と、
前記判別手段によって判別された使用者の動作状態に応じて、前記楽曲再生リス ト保持手段力 楽曲の再生に用いる楽曲再生リストを選択する選択手段と
を備えることを特徴とする体動検出装置。
[10] 請求項 1に記載の体動検出装置であって、
楽曲データの再生手段と、
前記体動検出手段によって検出される使用者の体動の周期に応じて、前記再生手 段における楽曲の再生速度を制御する再生制御手段と
を備え、楽曲再生装置としての機能を備えることを特徴とする体動検出装置。
[11] 請求項 1に記載の体動検出装置であって、
前記加速度センサは多軸のものであり、
前記鉛直成分抽出手段は、前記多軸の加速度センサからの検出出力である加速 度ベクトルから重力加速度ベクトルを算出し、前記多軸の加速度センサからの前記 加速度ベクトルと算出された前記重力加速度ベクトルとを用いて演算により加速度の 鉛直成分を抽出するものであることを特徴とする体動検出装置。
[12] 使用者の身体に装着するようにされる前記加速度センサ力 の検出出力から加速 度の鉛直成分を抽出する鉛直成分抽出工程と、
前記鉛直成分抽出工程において抽出した前記鉛直成分を高域成分と低域成分と に成分分解する分解工程と、
前記分解工程にお!ヽて分解した前記鉛直成分の前記低域成分に基づ!ヽてピーク 位置の候補を検出する検出工程と、 前記検出工程において検出したピーク位置の候補を含む所定範囲における前記 低域成分のエネルギーと前記高域成分のエネルギーとの比が所定値よりも小さい場 合に前記ピーク位置の候補をピーク位置の候補として特定する特定工程と、 前記特定工程にぉ 、て特定したピーク位置の候補に基づ!、て、前記使用者の体 動を検出する体動検出工程と
を有することを特徴とする体動検出方法。
[13] 請求項 12に記載の体動検出方法であって、
前記特定工程にぉ 、て特定した各ピーク位置の候補毎に、ピーク位置の候補を含 む所定範囲を設定し、設定した所定範囲毎の波形を対象として他の所定範囲の波 形とのマッチングを行い、類似度が高い場合に、対象となった波形に含まれるピーク 位置の候補をピーク位置として決定する決定工程を有し、
前記体動検出工程においては、前記決定工程において決定したピーク位置に基 づいて、使用者の体動を検出することを特徴とする体動検出方法。
[14] 請求項 12または請求項 13に記載の体動検出方法であって、
前記特定工程において特定した複数のピーク位置の候補力 なる時系列パターン について、あるいは、前記決定工程において決定した複数のピーク位置からなる時 系列パターンについて解析することにより、ステップ間隔を推定する間隔推定工程と 前記間隔推定工程において推定した前記ステップ間隔に基づいて、少なくとも「静 止」、「歩行 Z走行」、「不定」の 3状態を使用者の動作状態として判別する判別工程 と
を有することを特徴とする体動検出方法。
[15] 請求項 12に記載の体動検出方法であって、
使用者の身体に装着するようにされる前記加速度センサは多軸のものであり、 前記鉛直成分抽出工程においては、前記多軸の加速度センサからの検出出力で ある加速度ベクトル力 重力加速度ベクトルを算出し、前記多軸の加速度センサから の前記加速度ベクトルと算出された前記重力加速度ベクトルとを用いて演算により加 速度の鉛直成分を抽出することを特徴とする体動検出方法。
[16] 使用者の身体に装着するようにされる加速度センサを備え、当該加速度センサから の検出出力を用いて使用者の体動を検出する体動検出装置に搭載されたコンビュ ータに、
前記加速度センサ力ゝらの検出出力から加速度の鉛直成分を抽出する鉛直成分抽 出ステップと、
前記鉛直成分抽出ステップにおいて抽出した前記鉛直成分を高域成分と低域成 分とに成分分解する分解ステップと、
前記分解ステップにお ヽて分解した前記鉛直成分の前記低域成分に基づ!ヽてピ ーク位置の候補を検出する検出ステップと、
前記検出ステップにおいて検出したピーク位置の候補を含む所定範囲における前 記低域成分のエネルギーと前記高域成分のエネルギーとの比が所定値よりも小さい 場合に前記ピーク位置の候補をピーク位置の候補として特定する特定ステップと、 前記特定ステップにお 、て特定したピーク位置の候補に基づ 、て、前記使用者の 体動を検出する体動検出ステップと
を実行させることを特徴とする体動検出プログラム。
[17] 請求項 16に記載の体動検出プログラムであって、
前記特定ステップにお 、て特定した各ピーク位置の候補毎に、ピーク位置の候補 を含む所定範囲を設定し、設定した所定範囲毎の波形を対象として他の所定範囲の 波形とのマッチングを行い、類似度が高い場合に、対象となった波形に含まれるピー ク位置の候補をピーク位置として決定する決定ステップを備え、
前記体動検出ステップにお 、ては、前記決定ステップにお 、て決定したピーク位 置に基づいて、使用者の体動を検出することを特徴とする体動検出プログラム。
[18] 請求項 16または請求項 17に記載の体動検出プログラムであって、
前記特定ステップにおいて特定した複数のピーク位置の候補力 なる時系列バタ ーンについて、あるいは、前記決定ステップにおいて決定した複数のピーク位置から なる時系列パターンについて解析することにより、ステップ間隔を推定する間隔推定 ステップと、
前記間隔推定ステップにおいて推定した前記ステップ間隔に基づいて、少なくとも「 静止」、「歩行 Z走行」、「不定」の 3状態を使用者の動作状態として判別する判別ス テツプと
を有することを特徴とする体動検出プログラム。
請求項 16に記載の体動検出プログラムであって、
使用者の身体に装着するようにされる前記加速度センサは多軸のものであり、 前記鉛直成分抽出ステップにおいては、前記多軸の加速度センサからの検出出力 である加速度ベクトルから重力加速度ベクトルを算出し、前記多軸の加速度センサか らの前記加速度ベクトルと算出された前記重力加速度ベクトルとを用いて演算により 加速度の鉛直成分を抽出することを特徴とする体動検出プログラム。
PCT/JP2007/054710 2006-03-14 2007-03-09 体動検出装置、体動検出方法および体動検出プログラム WO2007105648A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN2007800070178A CN101394787B (zh) 2006-03-14 2007-03-09 人体运动检测设备和人体运动检测方法
US12/279,557 US7983872B2 (en) 2006-03-14 2007-03-09 Body movement detector, body movement detection method and body movement detection program
EP07738195A EP1994883A4 (en) 2006-03-14 2007-03-09 BODY MOTION DETECTOR, METHOD AND PROGRAM FOR DETECTING BODY MOVEMENT
BRPI0708620-2A BRPI0708620A2 (pt) 2006-03-14 2007-03-09 aparelho de detecção de movimento do corpo, método de detecção de movimento de corpo, e, programa de detecção de movimento de corpo

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2006-069343 2006-03-14
JP2006069343A JP4904861B2 (ja) 2006-03-14 2006-03-14 体動検出装置、体動検出方法および体動検出プログラム

Publications (1)

Publication Number Publication Date
WO2007105648A1 true WO2007105648A1 (ja) 2007-09-20

Family

ID=38509471

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2007/054710 WO2007105648A1 (ja) 2006-03-14 2007-03-09 体動検出装置、体動検出方法および体動検出プログラム

Country Status (10)

Country Link
US (1) US7983872B2 (ja)
EP (1) EP1994883A4 (ja)
JP (1) JP4904861B2 (ja)
KR (1) KR20080104123A (ja)
CN (1) CN101394787B (ja)
BR (1) BRPI0708620A2 (ja)
MY (1) MY150291A (ja)
RU (1) RU2442534C2 (ja)
TW (1) TW200821551A (ja)
WO (1) WO2007105648A1 (ja)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008073077A (ja) * 2006-09-19 2008-04-03 Mitsubishi Chemicals Corp データ処理装置,データ処理方法及びデータ処理プログラム
CN102056541A (zh) * 2008-06-12 2011-05-11 全球动力股份有限公司 运动机能减退和/或运动机能亢进状态的检测
EP2341821A1 (en) * 2008-09-23 2011-07-13 Koninklijke Philips Electronics N.V. Power measurement method and apparatus
CN103679512A (zh) * 2013-12-19 2014-03-26 深圳市宇恒互动科技开发有限公司 根据用户需求推送广告的可分离控制系统与装置
JP2017504440A (ja) * 2014-01-30 2017-02-09 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. ユーザの動きの測定値における歩行の検出の改善
JPWO2017175720A1 (ja) * 2016-04-08 2018-12-20 シャープ株式会社 行動判定装置および行動判定方法
US10292635B2 (en) 2013-03-01 2019-05-21 Global Kinetics Pty Ltd System and method for assessing impulse control disorder
US10736577B2 (en) 2014-03-03 2020-08-11 Global Kinetics Pty Ltd Method and system for assessing motion symptoms
RU2736005C1 (ru) * 2020-01-22 2020-11-11 Федеральное Государственное Бюджетное Образовательное Учреждение Высшего Образования "Тюменский Государственный Медицинский Университет" Министерства Здравоохранения Российской Федерации Способ центильного распределения и парциальной оценки уровня двигательной активности у детей старшего дошкольного и младшего школьного возраста

Families Citing this family (120)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7867141B2 (en) * 2004-07-21 2011-01-11 Panasonic Electric Works Co., Ltd. Physical activity measuring system
US8055469B2 (en) * 2006-03-03 2011-11-08 Garmin Switzerland Gmbh Method and apparatus for determining the attachment position of a motion sensing apparatus
US8206325B1 (en) 2007-10-12 2012-06-26 Biosensics, L.L.C. Ambulatory system for measuring and monitoring physical activity and risk of falling and for automatic fall detection
JP5416896B2 (ja) * 2007-11-15 2014-02-12 シチズン・システムズ株式会社 体動検出装置
JP2009125506A (ja) * 2007-11-27 2009-06-11 Panasonic Electric Works Co Ltd 歩容改善支援システム
US7877226B2 (en) * 2008-03-03 2011-01-25 Idt Technology Limited Apparatus and method for counting exercise repetitions
JP5176622B2 (ja) * 2008-03-18 2013-04-03 オムロンヘルスケア株式会社 歩数計
JP5608958B2 (ja) * 2008-03-31 2014-10-22 富士通株式会社 移動方向算出装置、移動方向算出プログラムおよび移動方向算出方法
US9482755B2 (en) 2008-11-17 2016-11-01 Faro Technologies, Inc. Measurement system having air temperature compensation between a target and a laser tracker
JP5155117B2 (ja) * 2008-11-21 2013-02-27 ラピスセミコンダクタ株式会社 重力軸判定装置及びこれを用いたモバイル端末装置
WO2010073684A1 (ja) * 2008-12-26 2010-07-01 オムロンヘルスケア株式会社 歩数検出システム及び歩数検出方法及び活動量計
JP4739439B2 (ja) * 2009-04-02 2011-08-03 株式会社タニタ 体動検出装置及び体動検出方法
JP5417970B2 (ja) * 2009-04-28 2014-02-19 富士通株式会社 歩数計及び歩数計数方法
US8744121B2 (en) 2009-05-29 2014-06-03 Microsoft Corporation Device for identifying and tracking multiple humans over time
CN101594208B (zh) * 2009-06-24 2014-04-09 中兴通讯股份有限公司 一种配置预编码矩阵的方法
EP2281612B1 (en) * 2009-08-04 2016-10-05 IDT Technology Limited Apparatus and method for counting exercise repetitions
CN102548474B (zh) * 2009-09-30 2014-09-03 三菱化学株式会社 针对体动信号的信息处理方法、信息处理系统及信息处理装置
KR101046604B1 (ko) * 2009-11-19 2011-07-06 한국산업기술시험원 주파수 디텍팅 기능을 구비한 진동 측정 장치 및 그 동작 방법
CN101879066A (zh) * 2010-03-08 2010-11-10 北京英福生科技有限公司 一种运动监测仪和运动健康数据的监测与传输方法
US8979665B1 (en) 2010-03-22 2015-03-17 Bijan Najafi Providing motion feedback based on user center of mass
JP5617299B2 (ja) * 2010-03-25 2014-11-05 オムロンヘルスケア株式会社 活動量計、制御プログラムおよび活動種類特定方法
JP5757742B2 (ja) * 2010-03-25 2015-07-29 セイコーインスツル株式会社 電子機器、歩数計、およびプログラム
CN101894252B (zh) * 2010-03-29 2012-12-05 天津大学 基于3轴加速度传感器信号的步行运动分类方法
US9772394B2 (en) 2010-04-21 2017-09-26 Faro Technologies, Inc. Method and apparatus for following an operator and locking onto a retroreflector with a laser tracker
US9377885B2 (en) 2010-04-21 2016-06-28 Faro Technologies, Inc. Method and apparatus for locking onto a retroreflector with a laser tracker
US9400170B2 (en) 2010-04-21 2016-07-26 Faro Technologies, Inc. Automatic measurement of dimensional data within an acceptance region by a laser tracker
US8724119B2 (en) 2010-04-21 2014-05-13 Faro Technologies, Inc. Method for using a handheld appliance to select, lock onto, and track a retroreflector with a laser tracker
US8422034B2 (en) 2010-04-21 2013-04-16 Faro Technologies, Inc. Method and apparatus for using gestures to control a laser tracker
US8537371B2 (en) 2010-04-21 2013-09-17 Faro Technologies, Inc. Method and apparatus for using gestures to control a laser tracker
US8619265B2 (en) 2011-03-14 2013-12-31 Faro Technologies, Inc. Automatic measurement of dimensional data with a laser tracker
JP5459179B2 (ja) * 2010-05-14 2014-04-02 カシオ計算機株式会社 移動状態判別装置、移動状態判別方法およびプログラム
TWI494078B (zh) * 2010-06-10 2015-08-01 Univ Nat Cheng Kung 健康照護監測系統
TWI490011B (zh) * 2010-07-01 2015-07-01 Ind Tech Res Inst 運動分析系統與方法
JP2012019927A (ja) * 2010-07-14 2012-02-02 Nippon Telegr & Teleph Corp <Ntt> 周期推定装置、周期推定方法、及びプログラム
JP5216825B2 (ja) * 2010-09-03 2013-06-19 京セラ株式会社 携帯機器、歩数カウント方法、重力方向検出方法、歩数カウントプログラムおよび重力方向検出プログラム
JP2012063208A (ja) * 2010-09-15 2012-03-29 Fujitsu Ltd 状態判定装置および状態判定プログラム
JPWO2012036135A1 (ja) * 2010-09-17 2014-02-03 三菱化学株式会社 情報処理方法、情報処理装置、出力装置、情報処理システム、情報処理用プログラムおよび同プログラムを記録したコンピュータ読み取り可能な記録媒体
TWI486880B (zh) 2010-09-20 2015-06-01 Ind Tech Res Inst 運動模式之自動辨識方法
JP2012107992A (ja) * 2010-11-17 2012-06-07 Sony Corp 歩行状況検出装置、歩行状況検出方法及び歩行状況検出プログラム
US8753275B2 (en) 2011-01-13 2014-06-17 BioSensics LLC Intelligent device to monitor and remind patients with footwear, walking aids, braces, or orthotics
BR112013020417A2 (pt) * 2011-02-10 2016-10-18 Dorsavi Pty Ltd aparelho e método para classificar a orientação de um corpo de um mamífero
US8666693B2 (en) * 2011-02-24 2014-03-04 Qualcomm Incorporated Low average velocity pedestrial motion identification
JP5364744B2 (ja) * 2011-03-01 2013-12-11 日本電信電話株式会社 加速度検出装置、方法及びプログラム
GB2518544A (en) 2011-03-03 2015-03-25 Faro Tech Inc Target apparatus and method
JP5779943B2 (ja) * 2011-03-31 2015-09-16 ソニー株式会社 情報処理装置、情報処理方法、プログラム、及び記録媒体
US9164173B2 (en) 2011-04-15 2015-10-20 Faro Technologies, Inc. Laser tracker that uses a fiber-optic coupler and an achromatic launch to align and collimate two wavelengths of light
JP2014516409A (ja) 2011-04-15 2014-07-10 ファロ テクノロジーズ インコーポレーテッド レーザトラッカの改良位置検出器
US9686532B2 (en) 2011-04-15 2017-06-20 Faro Technologies, Inc. System and method of acquiring three-dimensional coordinates using multiple coordinate measurement devices
US9482529B2 (en) 2011-04-15 2016-11-01 Faro Technologies, Inc. Three-dimensional coordinate scanner and method of operation
CN102263811A (zh) * 2011-04-28 2011-11-30 华南理工大学 智能健康呵护机
CN102290045B (zh) * 2011-05-13 2013-05-01 北京瑞信在线系统技术有限公司 一种控制音乐节奏的方法、装置及移动终端
DE102011079187A1 (de) * 2011-07-14 2013-01-17 Bayerische Motoren Werke Aktiengesellschaft Verfahren und Vorrichtung zur Schritterkennung eines Fußgängers für ein tragbares Endgerät
US8914037B2 (en) 2011-08-11 2014-12-16 Qualcomm Incorporated Numerically stable computation of heading without a reference axis
US20130046505A1 (en) * 2011-08-15 2013-02-21 Qualcomm Incorporated Methods and apparatuses for use in classifying a motion state of a mobile device
JP2013050308A (ja) * 2011-08-30 2013-03-14 Sony Corp 情報処理装置、情報処理方法、プログラム、及び記録媒体
JP2014223098A (ja) * 2011-09-27 2014-12-04 テルモ株式会社 活動量計及びその処理方法
JP5740285B2 (ja) * 2011-10-31 2015-06-24 株式会社東芝 歩行分析装置及び歩行分析プログラム
US9569439B2 (en) 2011-10-31 2017-02-14 Elwha Llc Context-sensitive query enrichment
CN102566796B (zh) * 2011-12-07 2015-03-11 东蓝数码有限公司 智能终端上基于线加速度传感器的人机交互检测及处理方法
JP5358656B2 (ja) * 2011-12-08 2013-12-04 京セラ株式会社 携帯電子機器
CN103181766B (zh) * 2011-12-30 2015-05-06 华为终端有限公司 一种人体监护方法及装置
US9638507B2 (en) 2012-01-27 2017-05-02 Faro Technologies, Inc. Measurement machine utilizing a barcode to identify an inspection plan for an object
DE102012201498A1 (de) * 2012-02-02 2013-08-08 Robert Bosch Gmbh Bedienvorrichtung und Verfahren zum Betrieb einer Bedienvorrichtung
CN104144643A (zh) * 2012-03-02 2014-11-12 日本电气株式会社 动作判定装置、动作判定系统和动作判定方法
CN102818913B (zh) * 2012-07-31 2013-11-27 宋子健 一种人体运动信息检测装置及其检测方法
CN102836543A (zh) * 2012-08-28 2012-12-26 洪德伟 一种健身运动识别方法及相关装置
US20140074431A1 (en) * 2012-09-10 2014-03-13 Apple Inc. Wrist Pedometer Step Detection
US10447844B2 (en) 2012-12-14 2019-10-15 Apple Inc. Method and apparatus for automatically setting alarms and notifications
US10564178B2 (en) 2012-12-21 2020-02-18 Qualcomm Incorporated Swing compensation in step detection
US9210566B2 (en) * 2013-01-18 2015-12-08 Apple Inc. Method and apparatus for automatically adjusting the operation of notifications based on changes in physical activity level
WO2014129169A1 (ja) 2013-02-22 2014-08-28 旭化成株式会社 保持状態変化検知装置、保持状態変化検知方法、およびプログラム
US9041914B2 (en) 2013-03-15 2015-05-26 Faro Technologies, Inc. Three-dimensional coordinate scanner and method of operation
US9311789B1 (en) 2013-04-09 2016-04-12 BioSensics LLC Systems and methods for sensorimotor rehabilitation
US10352724B1 (en) 2013-05-03 2019-07-16 Apple Inc. Calibration factors for step frequency bands
JP6111837B2 (ja) * 2013-05-10 2017-04-12 オムロンヘルスケア株式会社 歩行姿勢計およびプログラム
CN103301618A (zh) * 2013-05-12 2013-09-18 安徽工程大学 长跑辅助训练装置
US9086281B1 (en) * 2013-05-21 2015-07-21 Amazon Technologies Using accelerometer data to determine movement direction
CN103310192B (zh) * 2013-06-06 2016-05-11 南京邮电大学 一种基于单轴加速度传感器的运动行为识别方法
CN104297519B (zh) * 2013-07-19 2017-11-28 广州三星通信技术研究有限公司 人体运动姿态识别方法和移动终端
CN103488285A (zh) * 2013-08-22 2014-01-01 英华达(南京)科技有限公司 互动装置及方法
CN104580660B (zh) * 2013-10-12 2017-07-21 深圳市汇顶科技股份有限公司 一种移动智能终端及其计步方法、系统
JP6277674B2 (ja) * 2013-10-30 2018-02-14 富士通株式会社 情報処理装置、情報処理装置の作動方法及びプログラム
JP6424424B2 (ja) * 2013-12-02 2018-11-21 Tdk株式会社 行動推定装置及び活動量計
CN103727959B (zh) 2013-12-31 2016-09-14 歌尔声学股份有限公司 计步方法及装置
CN103767710B (zh) * 2013-12-31 2015-12-30 歌尔声学股份有限公司 人体运动状态监视方法和装置
TWI502167B (zh) * 2014-02-25 2015-10-01 Acer Inc 計步方法及其電子裝置
CN104880198A (zh) * 2014-02-27 2015-09-02 宏碁股份有限公司 计步方法及其电子装置
US11100234B2 (en) * 2014-06-13 2021-08-24 Hitachi Systems, Ltd. Work recording apparatus, system, program, and method preventing confidential information leaks
US9395174B2 (en) 2014-06-27 2016-07-19 Faro Technologies, Inc. Determining retroreflector orientation by optimizing spatial fit
JP6365031B2 (ja) * 2014-07-07 2018-08-01 オムロンヘルスケア株式会社 活動量測定装置、活動量の測定方法、活動量の測定プログラム
US10652696B2 (en) * 2014-07-30 2020-05-12 Trusted Positioning, Inc. Method and apparatus for categorizing device use case for on foot motion using motion sensor data
CN104121925A (zh) * 2014-08-08 2014-10-29 沈迪 一种计步方法和装置
SG11201704170TA (en) * 2014-11-27 2017-06-29 Razer (Asia-Pacific) Pte Ltd Step counter devices and step counting methods
US10197416B2 (en) * 2015-01-21 2019-02-05 Quicklogic Corporation Multiple axis wrist worn pedometer
JP6557995B2 (ja) * 2015-03-02 2019-08-14 富士通株式会社 計測プログラム、計測装置及び計測方法
RU2593983C1 (ru) * 2015-03-13 2016-08-10 Общество С Ограниченной Ответственностью "Хилби" Способ определения вида двигательной активности человека и устройство для его осуществления
JP6044670B2 (ja) * 2015-04-23 2016-12-14 Tdk株式会社 歩行分析装置
WO2016176301A2 (en) * 2015-04-30 2016-11-03 Honeywell International Inc. System for integrating multiple sensor data to predict a fall risk
JP6382775B2 (ja) * 2015-08-17 2018-08-29 日本電信電話株式会社 歩数計測装置、歩数計測方法、およびプログラム
TWI588743B (zh) * 2015-10-19 2017-06-21 Portable electronic device
JP6505614B2 (ja) * 2016-01-15 2019-04-24 株式会社日立製作所 トレーニング分類システム、トレーニング分類方法およびトレーニング分類サーバ
CN105953814B (zh) * 2016-06-21 2019-01-11 Oppo广东移动通信有限公司 一种计步方法、计步装置及移动终端
BR112019003561B1 (pt) 2016-08-31 2022-11-22 Apple Inc Sistema e método para melhorar uma exatidão de um dispositivo para ser usado junto ao corpo e determinar um movimento do braço de um usuário
US11896368B2 (en) 2016-08-31 2024-02-13 Apple Inc. Systems and methods for determining swimming metrics
US10041800B2 (en) 2016-09-23 2018-08-07 Qualcomm Incorporated Pedestrian sensor assistance in a mobile device during typical device motions
CN106491138B (zh) * 2016-10-26 2019-04-09 歌尔科技有限公司 一种运动状态检测方法及装置
CN106598537A (zh) * 2016-11-16 2017-04-26 上海斐讯数据通信技术有限公司 移动终端音乐播放的控制方法、系统及移动终端
US10120455B2 (en) 2016-12-28 2018-11-06 Industrial Technology Research Institute Control device and control method
JP6911418B2 (ja) * 2017-03-17 2021-07-28 カシオ計算機株式会社 表示装置、電子時計、表示方法及びプログラム
US11051720B2 (en) * 2017-06-01 2021-07-06 Apple Inc. Fitness tracking for constrained-arm usage
CN107273857B (zh) * 2017-06-19 2021-03-02 深圳市酷浪云计算有限公司 运动动作的识别方法及装置、电子设备
US11363970B2 (en) 2017-10-10 2022-06-21 Hunter Cronin Hand-held dexterity testing apparatus
US10854066B2 (en) 2018-04-12 2020-12-01 Apple Inc. Methods and systems for disabling sleep alarm based on automated wake detection
JP7080722B2 (ja) * 2018-05-17 2022-06-06 パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカ 検知方法、検知装置及び検知システム
CN109870172B (zh) * 2019-02-25 2020-12-18 广州市香港科大霍英东研究院 计步检测方法、装置、设备及存储介质
US11937904B2 (en) 2019-09-09 2024-03-26 Apple Inc. Detecting the end of cardio machine activities on a wearable device
CN112945257A (zh) * 2019-12-11 2021-06-11 瑞昱半导体股份有限公司 计步装置及方法
CN111166322A (zh) * 2020-01-17 2020-05-19 哈尔滨工业大学 心电信号特征波提取方法
CN111504433A (zh) * 2020-04-23 2020-08-07 永康龙飘传感科技有限公司 一种用于评估移动时动物的体重的方法及装置
JP7173102B2 (ja) * 2020-07-01 2022-11-16 カシオ計算機株式会社 情報処理装置、情報処理方法及びプログラム

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1142220A (ja) * 1997-04-17 1999-02-16 Kaajiopeeshingu Res Lab:Kk 体動検出方法及びその装置
JP2004141669A (ja) 2000-10-16 2004-05-20 Omron Healthcare Co Ltd 体動検出装置
JP2005140533A (ja) 2003-11-04 2005-06-02 Tanita Corp 体動測定装置
JP2006175206A (ja) * 2004-11-29 2006-07-06 National Institute Of Advanced Industrial & Technology 身体状態検出装置、その検出方法及び検出プログラム
JP2006293861A (ja) * 2005-04-13 2006-10-26 Matsushita Electric Works Ltd 歩数計

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0678902A (ja) * 1992-09-04 1994-03-22 Res:Kk 移動手段判別法
JPH11142220A (ja) 1997-11-13 1999-05-28 Canon Inc 流動体量測定装置及び現像装置
CN1137661C (zh) * 2000-03-06 2004-02-11 郑红 分布式电磁人体运动检测装置
JP2002200059A (ja) * 2000-12-28 2002-07-16 Terumo Corp 体動検出装置及びそれを用いた医療装置及び体動検出方法
KR100786703B1 (ko) * 2004-07-24 2007-12-21 삼성전자주식회사 가속도 센서를 이용한 운동량 측정장치 및 방법
JP4364157B2 (ja) * 2005-04-22 2009-11-11 トレックス・セミコンダクター株式会社 落下検出装置
JP4967368B2 (ja) * 2006-02-22 2012-07-04 ソニー株式会社 体動検出装置、体動検出方法および体動検出プログラム

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1142220A (ja) * 1997-04-17 1999-02-16 Kaajiopeeshingu Res Lab:Kk 体動検出方法及びその装置
JP2004141669A (ja) 2000-10-16 2004-05-20 Omron Healthcare Co Ltd 体動検出装置
JP2005140533A (ja) 2003-11-04 2005-06-02 Tanita Corp 体動測定装置
JP2006175206A (ja) * 2004-11-29 2006-07-06 National Institute Of Advanced Industrial & Technology 身体状態検出装置、その検出方法及び検出プログラム
JP2006293861A (ja) * 2005-04-13 2006-10-26 Matsushita Electric Works Ltd 歩数計

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008073077A (ja) * 2006-09-19 2008-04-03 Mitsubishi Chemicals Corp データ処理装置,データ処理方法及びデータ処理プログラム
CN102056541A (zh) * 2008-06-12 2011-05-11 全球动力股份有限公司 运动机能减退和/或运动机能亢进状态的检测
US11596327B2 (en) 2008-06-12 2023-03-07 Global Kinetics Pty Ltd Detection of hypokinetic and hyperkinetic states
US9826921B2 (en) 2008-06-12 2017-11-28 Global Kinetics Corporation Limited Detection of hypokinetic and hyperkinetic states
US9730617B2 (en) 2008-09-23 2017-08-15 Koninklijke Philips N.V. Power measurement method and apparatus
CN102164532A (zh) * 2008-09-23 2011-08-24 皇家飞利浦电子股份有限公司 功率测量方法和装置
EP2341821A1 (en) * 2008-09-23 2011-07-13 Koninklijke Philips Electronics N.V. Power measurement method and apparatus
US10292635B2 (en) 2013-03-01 2019-05-21 Global Kinetics Pty Ltd System and method for assessing impulse control disorder
CN103679512A (zh) * 2013-12-19 2014-03-26 深圳市宇恒互动科技开发有限公司 根据用户需求推送广告的可分离控制系统与装置
JP2017504440A (ja) * 2014-01-30 2017-02-09 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. ユーザの動きの測定値における歩行の検出の改善
US11039760B2 (en) 2014-01-30 2021-06-22 Koninklijke Philips N.V. Detection of walking in measurements of the movement of a user
US10736577B2 (en) 2014-03-03 2020-08-11 Global Kinetics Pty Ltd Method and system for assessing motion symptoms
JPWO2017175720A1 (ja) * 2016-04-08 2018-12-20 シャープ株式会社 行動判定装置および行動判定方法
RU2736005C1 (ru) * 2020-01-22 2020-11-11 Федеральное Государственное Бюджетное Образовательное Учреждение Высшего Образования "Тюменский Государственный Медицинский Университет" Министерства Здравоохранения Российской Федерации Способ центильного распределения и парциальной оценки уровня двигательной активности у детей старшего дошкольного и младшего школьного возраста

Also Published As

Publication number Publication date
CN101394787A (zh) 2009-03-25
TW200821551A (en) 2008-05-16
MY150291A (en) 2013-12-31
KR20080104123A (ko) 2008-12-01
EP1994883A1 (en) 2008-11-26
JP4904861B2 (ja) 2012-03-28
TWI325491B (ja) 2010-06-01
EP1994883A4 (en) 2012-02-22
US7983872B2 (en) 2011-07-19
JP2007244495A (ja) 2007-09-27
CN101394787B (zh) 2010-06-09
US20090240461A1 (en) 2009-09-24
BRPI0708620A2 (pt) 2011-06-07
RU2008136904A (ru) 2010-03-20
RU2442534C2 (ru) 2012-02-20

Similar Documents

Publication Publication Date Title
WO2007105648A1 (ja) 体動検出装置、体動検出方法および体動検出プログラム
JP4967368B2 (ja) 体動検出装置、体動検出方法および体動検出プログラム
JP4811046B2 (ja) コンテンツの再生装置、オーディオ再生機およびコンテンツの再生方法
JP4389821B2 (ja) 体動検出装置、コンテンツ再生装置、体動検出方法およびコンテンツ再生方法
US20110054833A1 (en) Processing motion sensor data using accessible templates
EP1583073A1 (en) Audio generating method and apparatus based on motion
US20210068713A1 (en) Detecting swimming activities on a wearable device
JP6788350B2 (ja) 多軸慣性センサを用いた手首装着型歩数計
US20120303319A1 (en) Pedometer device and step detection method
US20150120022A1 (en) Training supporting apparatus and system for supporting training of walking and/or running
JP2018523200A (ja) ウェアラブルコンピューティングデバイスの、微細運動に基づく入力ジェスチャ制御のための技術
US20070157798A1 (en) Apparatus and method for reproducing audio signal
JP2013252366A (ja) 動作判定方法、携帯端末およびプログラム
JPWO2017179423A1 (ja) 運動計測装置、情報処理装置及び運動計測方法
JP6940813B2 (ja) 歩数計測装置、歩数計測方法及びプログラム
EP3979908A1 (en) Method and apparatus for calibrating a user activity model used by a mobile device
Skoglund et al. Activity tracking using ear-level accelerometers
CN116114265A (zh) 耳机和用于识别耳机是否插入到用户的耳朵中的方法

Legal Events

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

Ref document number: 07738195

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 6343/DELNP/2008

Country of ref document: IN

WWE Wipo information: entry into national phase

Ref document number: 1020087019852

Country of ref document: KR

WWE Wipo information: entry into national phase

Ref document number: 12279557

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 200780007017.8

Country of ref document: CN

WWE Wipo information: entry into national phase

Ref document number: 2007738195

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2008136904

Country of ref document: RU

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: PI0708620

Country of ref document: BR

Kind code of ref document: A2

Effective date: 20080905