JP6690343B2 - Measuring device, measuring method, and measuring program - Google Patents

Measuring device, measuring method, and measuring program Download PDF

Info

Publication number
JP6690343B2
JP6690343B2 JP2016059469A JP2016059469A JP6690343B2 JP 6690343 B2 JP6690343 B2 JP 6690343B2 JP 2016059469 A JP2016059469 A JP 2016059469A JP 2016059469 A JP2016059469 A JP 2016059469A JP 6690343 B2 JP6690343 B2 JP 6690343B2
Authority
JP
Japan
Prior art keywords
timing
signal
acceleration
acceleration signal
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2016059469A
Other languages
Japanese (ja)
Other versions
JP2017169837A (en
Inventor
雅行 広浜
雅行 広浜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Casio Computer Co Ltd
Original Assignee
Casio Computer Co Ltd
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 Casio Computer Co Ltd filed Critical Casio Computer Co Ltd
Priority to JP2016059469A priority Critical patent/JP6690343B2/en
Publication of JP2017169837A publication Critical patent/JP2017169837A/en
Application granted granted Critical
Publication of JP6690343B2 publication Critical patent/JP6690343B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、測定装置、測定方法、及び測定プログラムに関する。   The present invention relates to a measuring device, a measuring method, and a measuring program.

従来、人体に加速度センサを取り付けた状態で歩行したときの、加速度センサの進行方向に対する出力結果と、鉛直方向に対する出力結果とから、例えば歩行時における着地タイミングや離地タイミング、接地時間等に関するデータを測定することのできる測定方法が開示されている(例えば特許文献1参照)。   Conventionally, when walking with an acceleration sensor attached to a human body, output data for the traveling direction of the acceleration sensor and output results for the vertical direction, for example, data regarding landing timing, take-off timing, ground contact time, etc. during walking. There is disclosed a measuring method capable of measuring (see, for example, Patent Document 1).

特開2012−179114号公報JP, 2012-179114, A

しかしながら、上記特許文献1に開示されている測定方法では、走行時において、特に走り始めなどの速度が遅い状態や、スピードに乗り速度が速くなった状態では、加速度センサからの前後方向の出力結果の波形が安定せず、離地に関するデータを誤検出してしまう場合があった。
そこで、本発明の課題は、走行時においても離地に関するデータを精度良く測定することである。
However, in the measuring method disclosed in Patent Document 1, the output result of the acceleration sensor in the front-rear direction is low when the vehicle is running, especially when the speed is low, such as when the vehicle starts to run, or when the speed is fast. There was a case where the waveform of was not stable and erroneously detected the data regarding the departure.
Therefore, an object of the present invention is to accurately measure data regarding takeoff even during traveling.

前記課題を解決するために、本発明の一の態様によれば、
生体に装着され、前記生体の動きに対応する加速度を検出して前記加速度に対応する加速度信号を出力する加速度センサと、
前記生体が進行方向に移動しているときに前記加速度信号に基づいて、前記生体の足の動きに関する少なくとも一つの動作タイミングを取得する制御部と、
を備え、
前記制御部は、
前記加速度信号に基づいて、前記進行方向の加速度信号を第1加速度信号として算出するとともに、上下方向の加速度信号を第2加速度信号として算出し、
前記第1加速度信号の第1期間の移動平均を用いて前記第1加速度信号を平滑化した第1信号と、前記第1加速度信号の前記第1期間より短い第2期間の移動平均を用いて前記第1加速度信号を平滑化した第2信号と、前記第2加速度信号を二度積分した第3信号と、を算出し、
前記第1信号、前記第2信号及び前記第3信号に基づいて、前記生体の足が離地したタイミングと推定される離地タイミングを前記動作タイミングとして取得することを特徴とする測定装置が提供される。
In order to solve the above problems, according to one aspect of the present invention,
An acceleration sensor that is attached to a living body, detects an acceleration corresponding to the movement of the living body, and outputs an acceleration signal corresponding to the acceleration;
Based on the acceleration signal when the living body is moving in the traveling direction, a control unit that acquires at least one operation timing regarding the movement of the foot of the living body,
Equipped with
The control unit is
Based on the acceleration signal, the acceleration signal in the traveling direction is calculated as a first acceleration signal, and the acceleration signal in the vertical direction is calculated as a second acceleration signal,
A first signal obtained by smoothing the first acceleration signal using a moving average of the first acceleration signal in a first period, and a moving average of the first acceleration signal in a second period shorter than the first period are used. Calculating a second signal obtained by smoothing the first acceleration signal and a third signal obtained by integrating the second acceleration signal twice,
A measurement apparatus, characterized in that, based on the first signal, the second signal, and the third signal, a take-off timing estimated to be a timing at which the foot of the living body is taken off is acquired as the operation timing. To be done.

また、本発明の他の態様によれば
定装置における測定方法であって、
前記測定装置は、生体に装着され、前記生体の動きに対応する加速度を検出して前記加速度に対応する加速度信号を出力する加速度センサと、前記生体が進行方向に移動しているときに前記加速度信号に基づいて、前記生体の足の動きに関する少なくとも一つの動作タイミングを取得する制御部と、を有し、
前記加速度信号に基づいて、前記進行方向の加速度信号を第1加速度信号として算出するとともに、上下方向の加速度信号を第2加速度信号として算出し、
前記第1加速度信号の第1期間の移動平均を用いて前記第1加速度信号を平滑化した第1信号と、前記第1加速度信号の前記第1期間より短い第2期間の移動平均を用いて前記第1加速度信号を平滑化した第2信号と、前記第2加速度信号を二度積分した第3信号と、を算出し、
前記第1信号、前記第2信号及び前記第3信号に基づいて、前記生体の足が離地したタイミングと推定される離地タイミングを前記動作タイミングとして取得することを特徴とする測定方法が提供される。
According to another aspect of the present invention ,
A measuring method in measurement device,
The measuring device is mounted on a living body, detects an acceleration corresponding to the movement of the living body, and outputs an acceleration signal corresponding to the acceleration; and an acceleration sensor when the living body is moving in a traveling direction. A control unit that acquires at least one operation timing related to the movement of the legs of the living body based on a signal ,
Based on the acceleration signal, the acceleration signal in the traveling direction is calculated as a first acceleration signal, and the acceleration signal in the vertical direction is calculated as a second acceleration signal,
A first signal obtained by smoothing the first acceleration signal using a moving average of the first acceleration signal in a first period, and a moving average of the first acceleration signal in a second period shorter than the first period are used. Calculating a second signal obtained by smoothing the first acceleration signal and a third signal obtained by integrating the second acceleration signal twice,
A measurement method, characterized in that, based on the first signal, the second signal, and the third signal, a takeoff timing estimated to be a timing at which the foot of the living body is taken off is obtained as the operation timing. To be done.

また、本発明の他の態様によれば
定装置における測定プログラムであって、
前記測定装置は、生体に装着され、前記生体の動きに対応する加速度を検出して前記加速度に対応する加速度信号を出力する加速度センサと、前記生体が進行方向に移動しているときに前記加速度信号に基づいて、前記生体の足の動きに関する少なくとも一つの動作タイミングを取得する制御部と、を有し、
前記加速度信号に基づいて、前記進行方向の加速度信号を第1加速度信号として算出させるとともに、上下方向の加速度信号を第2加速度信号として算出させ、
前記第1加速度信号の第1期間の移動平均を用いて前記第1加速度信号を平滑化した第1信号と、前記第1加速度信号の前記第1期間より短い第2期間の移動平均を用いて前記第1加速度信号を平滑化した第2信号と、前記第2加速度信号を二度積分した第3信号と、を算出させ、
前記第1信号、前記第2信号及び前記第3信号に基づいて、前記生体の足が離地したタイミングと推定される離地タイミングを前記動作タイミングとして取得させることを特徴とする測定プログラムが提供される。
According to another aspect of the present invention ,
A measurement program in the measurement device,
The measuring device is mounted on a living body, detects an acceleration corresponding to the movement of the living body, and outputs an acceleration signal corresponding to the acceleration; and an acceleration sensor when the living body is moving in a traveling direction. A control unit that acquires at least one operation timing related to the movement of the legs of the living body based on a signal ,
Based on the acceleration signal, the acceleration signal in the traveling direction is calculated as a first acceleration signal, and the acceleration signal in the vertical direction is calculated as a second acceleration signal,
A first signal obtained by smoothing the first acceleration signal using a moving average of the first acceleration signal in a first period, and a moving average of the first acceleration signal in a second period shorter than the first period are used. Calculating a second signal obtained by smoothing the first acceleration signal and a third signal obtained by integrating the second acceleration signal twice,
A measurement program, characterized in that, based on the first signal, the second signal, and the third signal, a take-off timing estimated to be a timing at which the foot of the living body takes off is acquired as the operation timing. To be done.

本発明によれば、走行時においても離地に関するデータを精度良く測定することができる。   According to the present invention, it is possible to accurately measure data regarding takeoff even during traveling.

本実施形態に係る測定装置をユーザが装着した状態を示す説明図である。It is explanatory drawing which shows the state which the user mounted the measuring device which concerns on this embodiment. 本実施形態に係る本体部の主制御構成を示すブロック図である。It is a block diagram which shows the main control structure of the main-body part which concerns on this embodiment. 図1の測定装置で実行される測定処理の流れを示すフローチャートである。3 is a flowchart showing a flow of measurement processing executed by the measuring device of FIG. 1. 図1の測定装置で実行される着地タイミング設定処理の流れを示すフローチャートである。6 is a flowchart showing a flow of landing timing setting processing executed by the measuring apparatus of FIG. 1. 前後方向、左右方向、上下方向それぞれの加速度信号の波形の一部分を例示的に抜き出して示すグラフである。It is a graph which extracts and shows a part of waveform of the acceleration signal of each of the front-back direction, the left-right direction, and the up-down direction. 図1の測定装置で実行される離地タイミング設定処理の流れを示すフローチャートである。3 is a flowchart showing a flow of a takeoff timing setting process executed by the measuring apparatus of FIG. 1. 前後方向、上下方向それぞれの加速度信号の波形の一部分を例示的に抜き出して示すグラフである。It is a graph which extracts and shows a part of waveform of the acceleration signal of each of the front-back direction and the up-down direction. 図7に示すグラフのうち、28.9秒から29.1秒の間の一部分を拡大したグラフである。It is the graph which expanded a part between 28.9 seconds and 29.1 seconds among the graphs shown in FIG. 前後方向、上下方向それぞれの加速度信号の波形の一部分を例示的に抜き出して示すグラフである。It is a graph which extracts and shows a part of waveform of the acceleration signal of each of the front-back direction and the up-down direction. 前後方向、上下方向それぞれの加速度信号の波形の一部分を例示的に抜き出して示すグラフである。It is a graph which extracts and shows a part of waveform of the acceleration signal of each of the front-back direction and the up-down direction. 前後方向、上下方向それぞれの加速度信号の波形の一部分を例示的に抜き出して示すグラフである。It is a graph which extracts and shows a part of waveform of the acceleration signal of each of the front-back direction and the up-down direction. フォースプレートにより測定された接地時間と図1の測定装置により算出された接地時間の散布図である。FIG. 2 is a scatter diagram of the contact time measured by the force plate and the contact time calculated by the measuring device of FIG. 1. 光学測定装置により測定された接地時間と図1の測定装置により算出された接地時間の散布図である。FIG. 2 is a scatter diagram of the contact time measured by the optical measuring device and the contact time calculated by the measuring device of FIG. 1.

以下、本発明に係る測定装置1について説明する。
なお、以下に述べる実施形態には、本発明を実施するために技術的に好ましい種々の限定が付されているが、本発明の範囲を以下の実施形態及び図示例に限定するものではない。
Hereinafter, the measuring device 1 according to the present invention will be described.
It should be noted that various technically preferable limitations for carrying out the present invention are attached to the embodiments described below, but the scope of the present invention is not limited to the following embodiments and illustrated examples.

図1は、本実施形態に係る測定装置1をユーザが装着した状態を示す説明図である。
図1に示すように、測定装置1は、本体部2と、ベルト部3とを有しており、ベルト部3によって、ユーザの腰の位置で本体部2が固定されている。ここで、左右方向をX軸とし、前後方向をY軸とし、上下方向をZ軸とする。X軸においては左手方向を正、右手方向を負とする。Y軸においては進行方向逆向きを正とし、進行方向を負とする。Z軸においては上方向を正、下方向を負とする。
FIG. 1 is an explanatory diagram showing a state in which a user wears the measuring apparatus 1 according to the present embodiment.
As shown in FIG. 1, the measuring device 1 has a main body 2 and a belt 3, and the main body 2 is fixed by the belt 3 at the waist of the user. Here, the left-right direction is the X axis, the front-back direction is the Y axis, and the up-down direction is the Z axis. On the X-axis, the left-hand direction is positive and the right-hand direction is negative. On the Y-axis, the direction opposite to the traveling direction is positive and the traveling direction is negative. On the Z-axis, the upward direction is positive and the downward direction is negative.

図2は、本実施形態に係る本体部2の主制御構成を示すブロック図である。
図2に示すように本体部2は、加速度センサ4、通信部5、表示部6、操作部7、及びこれらを制御する制御部8を備えて構成されている。
加速度センサ4は、互いに直交する3軸方向の加速度を測定し、測定した加速度に対応する3軸方向の加速度信号を制御部8に出力するようになっている。
通信部5は、取得したデータを制御部8による制御に基づいて外部の情報端末に出力するものであり、例えばUSB端子などの有線式の通信部や、Bluetooth(登録商標)などの近距離無線通信規格を採用した通信部である。
表示部6は、取得したデータを制御部8による制御に基づいて表示するものであり、例えば液晶パネルなどである。
操作部7は、電源のON/OFFを切り替える電源ボタン(図示省略)、データ取得の開始/停止を指示する開始/停止ボタン(図示省略)、表示内容を切り替える表示切替ボタン(図示省略)等を備えており、この操作部7からの指示に基づいて制御部8は各部を制御するようになっている。
FIG. 2 is a block diagram showing the main control configuration of the main body unit 2 according to this embodiment.
As shown in FIG. 2, the main body 2 is configured to include an acceleration sensor 4, a communication unit 5, a display unit 6, an operation unit 7, and a control unit 8 that controls these.
The acceleration sensor 4 measures the accelerations in the directions of the three axes orthogonal to each other, and outputs the acceleration signals in the directions of the three axes corresponding to the measured accelerations to the control unit 8.
The communication unit 5 outputs the acquired data to an external information terminal under the control of the control unit 8. For example, a wired communication unit such as a USB terminal or a short-range wireless communication such as Bluetooth (registered trademark). It is a communication unit that adopts a communication standard.
The display unit 6 displays the acquired data under the control of the control unit 8, and is, for example, a liquid crystal panel or the like.
The operation unit 7 includes a power button (not shown) for switching the power on / off, a start / stop button (not shown) for instructing start / stop of data acquisition, a display switch button (not shown) for switching display contents, and the like. The control section 8 controls each section based on an instruction from the operation section 7.

制御部8は、例えば、CPU、ROM、RAM(いずれも図示せず)からなり、ROMに記録された処理プログラムをRAMに展開してCPUによりこの処理プログラムを実行するものである。
具体的には、制御部8は、加速度センサ4から出力される3軸方向の加速度信号を所定のサンプリング周期(例えば、200Hz)でサンプリングして、RAM中に記憶し、サンプリングした加速度信号を周知の手法により前後方向、左右方向、上下方向それぞれの加速度信号に変換し、変換後の前後方向、左右方向、上下方向のそれぞれの時間軸に対する波形を作成する。そして制御部8は、前後方向、左右方向、上下方向それぞれの加速度信号の波形に基づいて、ランニング動作におけるサイクルの各局面のデータを算出する。
The control unit 8 is composed of, for example, a CPU, a ROM, and a RAM (none of which is shown). The control program is expanded in the RAM and the CPU executes the process program.
Specifically, the control unit 8 samples the acceleration signals in the three-axis directions output from the acceleration sensor 4 at a predetermined sampling cycle (for example, 200 Hz), stores them in the RAM, and publicizes the sampled acceleration signals. Is converted into acceleration signals in the front-rear direction, the left-right direction, and the up-down direction, and the converted waveforms for the front-rear direction, the left-right direction, and the up-down direction are created. Then, the control unit 8 calculates data of each phase of the cycle in the running motion based on the waveforms of the acceleration signals in the front-rear direction, the left-right direction, and the up-down direction.

次に、本実施形態に係る測定装置1で実行される上記各局面のデータの測定処理について説明する。この測定処理によって本発明に係る測定方法が実行される。
図3は、測定処理の流れを示すフローチャートである。
なお、この処理ではユーザが所定距離だけ走った分の前後方向、左右方向、上下方向それぞれの加速度信号が、処理の実行前に取得されている場合を例示して説明する。例えば、開始/停止ボタンが操作されることによってデータ取得が開始されてから、ユーザが所定距離だけ走り終わり、開始/停止ボタンが再度操作されデータ取得が停止されると、制御部8は上記測定処理に関するプログラムを読み出して実行する。
図5、及び図7〜11は、前後方向、左右方向、及び上下方向の加速度信号の波形の一部分を例示的に抜き出して示すグラフである。以下の説明では、上記測定処理に加速度信号の波形がどう用いられているかを、図5、及び図7〜11を参照して説明する。また、本実施形態では、図5、及び図7〜11に示す部分にのみ上記測定処理が施される場合を例示して説明するが、当該測定処理は、前後方向、左右方向、及び上下方向それぞれの加速度信号の波形の全体を対象として実行される。
Next, the measurement processing of the data of each aspect described above, which is executed by the measuring apparatus 1 according to the present embodiment, will be described. The measuring method according to the present invention is executed by this measuring process.
FIG. 3 is a flowchart showing the flow of measurement processing.
Note that, in this process, an explanation will be given by exemplifying a case where the acceleration signals in the front-rear direction, the left-right direction, and the up-down direction corresponding to the amount of travel of the user by a predetermined distance are acquired before the process is executed. For example, when data acquisition is started by operating the start / stop button, the user finishes running for a predetermined distance, and the start / stop button is operated again to stop data acquisition, the control unit 8 causes the measurement to be performed. Read and execute the program related to the process.
5 and FIGS. 7 to 11 are graphs showing, as an example, a part of the waveform of the acceleration signal in the front-rear direction, the left-right direction, and the up-down direction. In the following description, how the waveform of the acceleration signal is used in the measurement process will be described with reference to FIGS. 5 and 7 to 11. In the present embodiment, the case where the measurement process is performed only on the portions shown in FIGS. 5 and 7 to 11 will be described as an example. It is executed for the entire waveform of each acceleration signal.

図3に示すように、上記測定処理が実行されると、制御部8は、先ず、着地タイミング設定処理を行う(ステップS1)。ここで、着地タイミングとは、ランニング動作におけるサイクルの一局面であり、走っているユーザの足(一方の足)が地面に触れたタイミングを指す。なお、着地タイミング設定処理の詳細については後述する。   As shown in FIG. 3, when the measurement process is executed, the control unit 8 first performs a landing timing setting process (step S1). Here, the landing timing is one phase of the cycle in the running motion, and refers to the timing when the running user's foot (one foot) touches the ground. The details of the landing timing setting process will be described later.

次いで、制御部8は、最下点タイミングの設定を行う(ステップS2)。具体的には、制御部8は、後述する高さ位置波形Tの極小値Tmin(図5参照)を示すタイミングを最下点タイミングとして設定する。ここで、最下点タイミングとは、ランニング動作におけるサイクルの一局面であり、走っているユーザの足(一方の足)が地面に触れている間において、腰が最も低い位置を通過したタイミングを指す。 Next, the control unit 8 sets the lowest point timing (step S2). Specifically, the control unit 8 sets the timing indicating the minimum value T min (see FIG. 5) of the height position waveform T described later as the lowest point timing. Here, the lowest point timing is one phase of the cycle in the running motion, and is the timing at which the waist passes the lowest position while the running user's foot (one foot) touches the ground. Point to.

次いで、制御部8は、離地タイミング設定処理を行う(ステップS3)。ここで、離地タイミングとは、ランニング動作におけるサイクルの一局面であり、走っているユーザの足(一方の足)が地面から離れたタイミングを指す。なお、離地タイミング設定処理の詳細については後述する。   Next, the control unit 8 performs a takeoff timing setting process (step S3). Here, the take-off timing is one phase of the cycle in the running motion, and refers to the timing at which the running user's foot (one foot) leaves the ground. The details of the departure timing setting process will be described later.

次いで、制御部8は、最高点タイミングの設定を行う(ステップS4)。具体的には、制御部8は、後述する高さ位置波形Tの極大値Tmax(図5参照)を示すタイミングを最高点タイミングとして設定する。ここで、最高点タイミングとは、ランニング動作におけるサイクルの一局面であり、走っているユーザの足(一方の足)が地面から離れた後、腰が最も高い位置を通過したタイミングを指す。 Next, the control unit 8 sets the highest point timing (step S4). Specifically, the control unit 8 sets the timing indicating the maximum value T max (see FIG. 5) of the height position waveform T described later as the highest point timing. Here, the highest point timing is one phase of a cycle in a running motion, and refers to the timing when the running user's foot (one foot) leaves the ground and then the waist passes the highest position.

次いで、制御部8は、着地タイミング設定処理によって設定された着地タイミングと、離地タイミング設定処理によって設定された離地タイミングと、の差分の時間を接地時間として算出する(ステップS5)。
これにより、ランニング動作におけるサイクルの各局面のデータ(着地タイミング、最下点タイミング、離地タイミング、最高点タイミング、及び接地時間)が算出されて、制御部8は測定処理を終了する。
Next, the control unit 8 calculates the time difference between the landing timing set by the landing timing setting process and the takeoff timing set by the takeoff timing setting process as the ground contact time (step S5).
As a result, data of each phase of the cycle in the running motion (landing timing, lowest point timing, takeoff timing, highest point timing, and ground contact time) is calculated, and the control unit 8 ends the measurement process.

次に、上記の測定処理における着地タイミング設定処理(ステップS1)の詳細について説明する。
図4は、着地タイミング設定処理の流れを示すフローチャートである。
図5は、前後方向、左右方向、及び上下方向の加速度信号の波形の一部分を例示的に抜き出して示すグラフである。以下の説明では、着地タイミング設定処理に加速度信号の波形がどう用いられているかを、図5を参照して説明する。また、本実施形態では、図5に示す部分にのみ着地タイミング設定処理が施される場合を例示して説明するが、着地タイミング設定処理は、前後方向、左右方向、及び上下方向それぞれの加速度信号の波形の全体を対象として実行される。
Next, details of the landing timing setting process (step S1) in the above measurement process will be described.
FIG. 4 is a flowchart showing the flow of landing timing setting processing.
FIG. 5 is a graph exemplarily showing a part of the waveforms of the acceleration signals in the front-rear direction, the left-right direction, and the up-down direction. In the following description, how the waveform of the acceleration signal is used in the landing timing setting process will be described with reference to FIG. In the present embodiment, the case where the landing timing setting process is performed only on the portion shown in FIG. 5 will be described as an example. However, the landing timing setting process is performed in the front-rear direction, the left-right direction, and the vertical direction acceleration signal. Is executed for the entire waveform of.

図4に示すように、着地タイミング設定処理が実行されると、制御部8は、先ず、上下方向の加速度信号(第2加速度信号)AccZ1に対して例えば移動平均等の周知の平滑化処理を行う(ステップS11)。
次いで、制御部8は、平滑化後の上下方向の加速度信号AccZ2の波形の極大値Zmaxを求め、この極大値Zmaxを示すタイミングを基準タイミング(第四基準タイミング)として時間軸を分割する(ステップS12)。図5では、極大値Zmaxに基づいて線P1,P2,P3で時間軸を分割している。この線P1と,線P2の間の分割領域、線P2と線P3との間の分割領域を第一分割領域R1と称す。
As shown in FIG. 4, when the landing timing setting process is executed, the control unit 8 first performs a known smoothing process such as a moving average on the vertical acceleration signal (second acceleration signal) AccZ1. Perform (step S11).
Next, the control unit 8 obtains a maximum value Z max of the waveform of the acceleration signal AccZ2 in the vertical direction after smoothing, and divides the time axis with the timing indicating the maximum value Z max as the reference timing (fourth reference timing). (Step S12). In FIG. 5, the time axis is divided by lines P1, P2, and P3 based on the maximum value Z max . The divided area between the line P1 and the line P2 and the divided area between the line P2 and the line P3 are referred to as a first divided area R1.

次いで、制御部8は、平滑化前の上下方向の加速度信号AccZ1を第一分割領域R1毎に二度積分することで、加速度センサ4の高さ位置を表す高さ位置波形Tを求める(ステップS13)。
次いで、制御部8は、各第一分割領域R1内にある高さ位置波形Tの極大値Tmaxを示すタイミングを基準タイミング(第一基準タイミング)として時間軸を分割する。図5では、極大値Tmaxに基づいて線P4,P5,P6で時間軸を分割している。この線P4,P5,P6を歩の切れ目(例えば奇数歩目と偶数歩目の境界)に設定する(ステップS14)。この線P4,P5,P6で分割された領域を第二分割領域と称す。なお、説明の便宜上、以下においては、連続する二つの第二分割領域のうち、先(時間的に前)の方を「先の第二分割領域R21」とし、後(時間的に後)の方を「後の第二分割領域R22」とする。
Next, the control unit 8 obtains a height position waveform T representing the height position of the acceleration sensor 4 by integrating the vertical acceleration signal AccZ1 before smoothing twice for each first divided region R1 (step S13).
Next, the control unit 8 divides the time axis with the timing indicating the maximum value T max of the height position waveform T in each first divided region R1 as the reference timing (first reference timing). In FIG. 5, the time axis is divided by lines P4, P5, and P6 based on the maximum value T max . The lines P4, P5 and P6 are set as step breaks (for example, boundaries between odd steps and even steps) (step S14). The area divided by the lines P4, P5 and P6 is referred to as a second divided area. Note that, for convenience of explanation, in the following, of the two consecutive second divided areas, the earlier (temporally earlier) one is referred to as the “first second divided area R21” and the latter (temporally later) One is referred to as a “second divided region R22 later”.

次いで、制御部8は、各第二分割領域R21,R22のうち前半部分(例えば歩の切れ目(線P4,P5)から、平滑化後の上下方向の加速度信号AccZ2の波形の極大値位置(線P2,P3))内において、前後方向の加速度信号(第1加速度信号)AccYの波形の正の極大値を探索する(ステップS15)。
図5では、先の第二分割領域R21に正の極大値Ymaxが一つあり、後の第二分割領域R22に正の極大値Ymaxが二つある場合を例示している。
Next, the control unit 8 controls the maximum value position (line of the waveform of the smoothed vertical acceleration signal AccZ2 from the first half of each of the second divided regions R21 and R22 (for example, a step break (lines P4 and P5)). P2, P3)), a positive maximum value of the waveform of the acceleration signal (first acceleration signal) AccY in the front-rear direction is searched (step S15).
In Figure 5, there is one positive maximum value Y max on the second divisional region R21 previous, positive maximum value Y max on the second divisional region R22 of the subsequent illustrates the case where there are two.

次いで、制御部8は、各第二分割領域R21,R22の前半部分にある正の極大値Ymaxが一つであるか否かを判断し、一つである場合にはその正の極大値Ymaxを特定しステップS18に移行し、二つ以上である場合にはステップS17に移行する(ステップS16)。
ここで、着地時においてはその着地動作による衝撃によって減速するために、前後方向の加速度信号AccYの波形には正の値にピーク(極大値)が生じることになる。この正の極大値Ymaxが、進行方向に対して減速を示す極値である。
しかし、走り方や速度によっては同様のピークが複数生じることもある。換言すると、着地時においては前後方向の加速度信号AccYの波形に少なくとも一つは正の極大値が生じることになる。このため、ステップS16では、正の極大値が一つであるか否かを判断している。なお、本実施形態では、前後方向の加速度信号AccYは、進行方向逆向きを正とし、進行方向を負としているため、前後方向の加速度信号AccYの波形における減速を示す極値が正の極大値Ymaxとなっているが、正負が逆の場合は減速を示す極値が負の極小値となる。
Next, the control unit 8 determines whether or not there is one positive maximum value Y max in the first half of each of the second divided regions R21 and R22, and when it is one, the positive maximum value Y max. Ymax is specified, the process proceeds to step S18, and if it is two or more, the process proceeds to step S17 (step S16).
Here, at the time of landing, deceleration is caused by the impact due to the landing operation, so that a positive peak (maximum value) occurs in the waveform of the acceleration signal AccY in the front-rear direction. This positive maximum value Y max is an extreme value that indicates deceleration in the traveling direction.
However, multiple similar peaks may occur depending on the running style and speed. In other words, at the time of landing, at least one of the waveforms of the acceleration signal AccY in the front-rear direction has a positive maximum value. Therefore, in step S16, it is determined whether or not there is one positive maximum value. In this embodiment, the acceleration signal AccY in the front-rear direction is positive in the reverse direction of the traveling direction and negative in the traveling direction. Therefore, the extreme value indicating the deceleration in the waveform of the acceleration signal AccY in the front-rear direction is a positive maximum value. Although it is Y max , the extreme value indicating deceleration becomes a negative minimum value when the positive and negative signs are opposite.

ステップS17では、制御部8は、左右方向の加速度信号AccXの波形に基づいて複数の正の極大値Ymaxから一つ特定する。例えば、走行による着地時においては、片足だけが着地することになるため、上体は左右方向に傾くことになる。このとき、無意識のうちに体勢を整えるために、左右方向に上体が揺れ動くことになる。つまり、左右方向の加速度信号AccXの波形においては、着地時に左右のブレを示す波形が現れることになる。この波形を基にすれば、着地時を特定することが可能である。
具体的には、制御部8は、後の第二分割領域R22の前半部分内において、所定間隔内に左右方向の加速度信号AccXの波形の極値Xm1,Xm2,Xm3が三つ存在し、なおかつ当該三つの極値Xm1,Xm2,Xm3のうち隣接する極値Xm1,Xm2,Xm3の差分が所定値以上であると、左右のブレを示す波形として認定する。そして、制御部8は三つの極値Xm1,Xm2,Xm3のうち最初に発生した極値Xm1に近い正の極大値Ymaxを特定する。
なお、「所定間隔」や、「所定値」については、実験やシミュレーションなどにより得られた値が用いられる。具体的には所定間隔としては40〜100msの範囲に収まる値が好ましく、70msがより好ましい。また、所定値としては5〜15m/sの範囲に収まる値が好ましく、10m/sがより好ましい。
In step S17, the control unit 8 identifies one from a plurality of positive maximum values Y max based on the waveform of the acceleration signal AccX in the left-right direction. For example, when landing due to running, only one foot will land, so the upper body leans in the left-right direction. At this time, in order to adjust the posture unconsciously, the upper body sways in the left-right direction. That is, in the waveform of the acceleration signal AccX in the left-right direction, a waveform showing a left-right blur appears when landing. Based on this waveform, it is possible to specify the landing time.
Specifically, the control unit 8 has three extreme values X m1 , X m2 , and X m3 of the waveform of the acceleration signal AccX in the left-right direction within a predetermined interval in the first half portion of the subsequent second divided region R22. However, if the difference between the adjacent extreme values X m1 , X m2 , and X m3 among the three extreme values X m1 , X m2 , and X m3 is equal to or more than a predetermined value, it is recognized as a waveform indicating left and right blur. Then, the control unit 8 specifies the positive maximum value Y max close to the first generated extreme value X m1 among the three extreme values X m1 , X m2 , and X m3 .
As the "predetermined interval" and the "predetermined value", values obtained by experiments and simulations are used. Specifically, the predetermined interval is preferably a value within the range of 40 to 100 ms, more preferably 70 ms. Further, preferably a value within the range of 5 to 15 m / s 2 as the predetermined value, 10 m / s 2 is more preferable.

次いで、制御部8は、特定した正の極大値Ymaxの位置(時間)を着地タイミングとして設定し(ステップS18)、着地タイミング設定処理を終了する。 Next, the control unit 8 sets the position (time) of the identified positive maximum value Y max as the landing timing (step S18), and ends the landing timing setting process.

次に、上記の測定処理における離地タイミング設定処理(ステップS3)の詳細について説明する。
図6は、離地タイミング設定処理の流れを示すフローチャートである。
図7、図9、図10は、前後方向、及び上下方向の加速度信号の波形の一部分を例示的に抜き出して示すグラフである。
図8は、図7に示すグラフのうち、28.9秒から29.1秒の間の一部分を拡大したグラフである。
以下の説明では、離地タイミング設定処理に加速度信号の波形がどう用いられているかを、図7〜図10を参照して説明する。また、本実施形態では、図7〜図10に示す部分にのみ離地タイミング設定処理が施される場合を例示して説明するが、離地タイミング設定処理は、前後方向、及び上下方向それぞれの加速度信号の波形の全体を対象として実行される。
Next, details of the departure timing setting process (step S3) in the above measurement process will be described.
FIG. 6 is a flowchart showing the flow of the departure timing setting process.
FIG. 7, FIG. 9, and FIG. 10 are graphs showing, as an example, a part of the waveform of the acceleration signal in the front-rear direction and the waveform in the up-down direction.
FIG. 8 is a graph obtained by enlarging a part of the graph shown in FIG. 7 between 28.9 seconds and 29.1 seconds.
In the following description, how the waveform of the acceleration signal is used in the departure timing setting process will be described with reference to FIGS. In the present embodiment, the case where the take-off timing setting process is performed only on the portions shown in FIGS. 7 to 10 will be described as an example, but the take-off timing setting process is performed in the front-rear direction and the vertical direction. It is executed for the entire waveform of the acceleration signal.

図6に示すように、離地タイミング設定処理が実行されると、制御部8は、先ず、前後方向の加速度信号(第1加速度信号)AccYに対して移動平均を用いた平滑化処理を行う(ステップS21)。
具体的には、制御部8は、前後方向の加速度信号AccYに対して、タップ数が11〜15(ここで、タップ数はサンプリング周波数が200Hzであって、サンプリング周期が5msである場合のサンプル数であって、移動平均を算出する期間が対象とする信号の前後25〜35msの期間(第1期間)である場合に対応する)の三角移動平均(第一の移動平均)を算出し、当該三角移動平均による平滑化処理後の第1信号(以下、FAccY_2と称す)を求めるとともに、タップ数が3〜5(タップ数はサンプリング周波数が200Hzである場合のサンプル数であって、移動平均を算出する期間が対象とする信号の前後5〜10msの期間(第2期間)である場合に対応する)の三角移動平均(第二の移動平均)を算出し、当該三角移動平均による平滑化処理後の第2信号(以下、FAccY_1と称する)を求める。ここで、三角移動平均は対象とする信号のノイズを低減する効果を有するものであり、周知の手法である。この三角移動平均による平滑化処理においては、タップ数(サンプル数)を多くする程、信号が平滑化される効果が大きくなる。
As illustrated in FIG. 6, when the take-off timing setting process is executed, the control unit 8 first performs a smoothing process using a moving average on the longitudinal acceleration signal (first acceleration signal) AccY. (Step S21).
Specifically, the control unit 8 sets the number of taps to 11 to 15 with respect to the acceleration signal AccY in the front-rear direction (here, the number of taps is a sampling frequency of 200 Hz and a sampling period of 5 ms. A number, which corresponds to the case where the period for calculating the moving average is a period (first period) of 25 to 35 ms before and after the signal of interest), the triangular moving average (first moving average) is calculated, The first signal (hereinafter referred to as FAccY_2) after the smoothing processing by the triangular moving average is obtained, and the number of taps is 3 to 5 (the number of taps is the number of samples when the sampling frequency is 200 Hz, (Corresponding to the case where the period for calculating is a period of 5 to 10 ms before and after the target signal (second period)) is calculated, and the triangular moving average (second moving average) is calculated. The second signal (hereinafter referred to as FAccY_1) after the smoothing processing by the moving average is obtained. Here, the triangular moving average has an effect of reducing noise of a target signal, and is a known method. In the smoothing process using the triangular moving average, the larger the number of taps (the number of samples), the greater the effect of smoothing the signal.

次いで、制御部8は、変数nを、サンプリング周期(例えば、5ms)毎のサンプリングタイミングに、時間の経過に応じて順次付した番号として、第二分割領域(上述の図5における線P4、P5間、又は、線P5、P6間の第二分割領域と同様に設定される領域であり、図7の線L1,L2間の領域)のうちの後半部分において、FAccY_2が負から正、すなわち加速から減速へと変化するタイミング(プラスゼロクロスタイミング;n=ZC)を探索する(ステップS22)。そして、制御部8は、探索して取得した上記タイミング(ZC)を第二基準タイミングとして設定する。図7〜図10では、第二基準タイミングを線LZCで示している。
ここで、プラスゼロクロスタイミング(ZC)を探索するに際して、図8の破線の丸で囲った領域に示すように、離地前にも、AccYの値がノイズや体の動き方によって負から正に急峻に変化するタイミングが発生してしまう場合がある。そして、タップ数が比較的小さいFAccY_1ではこの急峻な変化が残っていて、同様に負から正に急峻に変化する。これに対して、図8に示すように、タップ数が比較的大きいFAccY_2では、AccYやFAccY_1では負から正に急峻に変化する場合であっても、この急峻な変化が緩和されて、負から正に変化するタイミングは現れないようになる。そこで、本実施形態では、タップ数を比較的大きくしたFAccY_2を用いることによって、プラスゼロクロスタイミング(ZC)を適正に探索することができるようにしている。なお、上記のタップ数はサンプリング周波数が200Hzでサンプリング周期が5msの場合であり、タップ数は加速度信号のサンプリング周波数によって変化する。また、サンプリング数、サンプリング期間は上記の値に限定されるものではなく、上記の値の近傍であればよい。
Next, the control unit 8 sets the variable n as a number sequentially added to the sampling timing of each sampling cycle (for example, 5 ms) according to the passage of time, and the second divided region (lines P4 and P5 in FIG. 5 described above). Or a second divided area between the lines P5 and P6, which is set in the second half of the area between the lines L1 and L2 in FIG. 7, FAccY_2 is from negative to positive, that is, acceleration. A timing (plus zero crossing timing; n = ZC) at which the speed changes from to deceleration is searched (step S22). Then, the control unit 8 sets the above-mentioned timing (ZC) obtained by searching as the second reference timing. In FIGS. 7 to 10, the second reference timing is indicated by the line L ZC .
Here, when searching for plus-zero cross timing (ZC), as shown in the area surrounded by the broken line circle in FIG. 8, the value of AccY is changed from negative to positive depending on noise and the way the body moves even before takeoff. There may be a case where a sudden change timing occurs. Then, in FAccY_1 having a relatively small number of taps, this steep change remains, and similarly, the steep change from negative to positive occurs. On the other hand, as shown in FIG. 8, in the case of FAccY_2 having a relatively large number of taps, even in the case of a sharp change from negative to positive in AccY and FAccY_1, this abrupt change is relaxed, and The timing of positive change will not appear. Therefore, in the present embodiment, FAccY_2 having a relatively large number of taps is used so that the plus-zero cross timing (ZC) can be appropriately searched. The number of taps is for a sampling frequency of 200 Hz and a sampling period of 5 ms, and the number of taps changes depending on the sampling frequency of the acceleration signal. Further, the number of samplings and the sampling period are not limited to the above values, and may be in the vicinity of the above values.

次いで、制御部8は、第二分割領域(図7の線L1,L2で分割された領域)のうちの後半部分において、FAccY_2が極大値となるタイミング(LMX)を探索する(ステップS23)。そして、制御部8は、上記タイミング(LMX)を第三基準タイミングとして設定する。図7〜図10では、第三基準タイミングを線LLMXで示している。 Next, the control unit 8 searches for the timing (LMX) at which FAccY_2 has the maximum value in the latter half of the second divided area (area divided by the lines L1 and L2 in FIG. 7) (step S23). Then, the control unit 8 sets the timing (LMX) as the third reference timing. In FIGS. 7 to 10, the third reference timing is indicated by the line L LMX .

次いで、制御部8は、第二基準タイミングと第三基準タイミングとの間、すなわち線LZC,LLMXで分割された探索期間において、先ず、第二基準タイミングから1サンプリング周期後のサンプリングタイミング(n=ZC+1、ステップS24)から第三基準タイミング(LMX)の間のサンプリング周期毎のFAccY_1の値の、FAccY_1の1サンプリング周期前の値に対する差分(FAccY_1(n)−FAccY_1(n−1))を算出し、当該差分の値が予め定められた閾値TH未満であるか否かを判定する(ステップS25)。ここで、ステップS25の閾値処理は、例えば、ノイズや体の固有の動き等の影響によって、第二基準タイミングの後、極大値が一気に現れる場合や、途中で変曲点が見られる波形の後に極大値が現れる場合に、離地タイミングを正確に特定できるようにすることを目的としてなされている。また、当該閾値処理では、AccYのデータに含まれているノイズの影響を抑制するため、FAccY_1を用いている。ここで、閾値THは、0m/sに近い値(例えば、0.5m/s)とする。なお、閾値THは、FAccY_1の極大値に応じて段階的に設定しても良い。例えば、FAccY_1の極大値が5m/s未満の場合は0.5m/sとし、当該極大値が5m/s以上の場合は1.0m/sとする。 Then, the control unit 8 firstly, between the second reference timing and the third reference timing, that is, in the search period divided by the lines L ZC and L LMX , first, the sampling timing (sampling timing 1 sampling period after the second reference timing) ( n = ZC + 1, the difference between the value of FAccY_1 for each sampling cycle from the step S24) to the third reference timing (LMX) with respect to the value of FAccY_1 one sampling cycle before (FAccY_1 (n) -FAccY_1 (n-1)). Is calculated and it is determined whether or not the value of the difference is less than a predetermined threshold TH (step S25). Here, the threshold value processing of step S25 is performed, for example, after the second reference timing, when the maximum value appears at once due to the influence of noise or peculiar movement of the body, or after the waveform where an inflection point is seen in the middle. The purpose is to enable the departure timing to be accurately specified when the maximum value appears. Further, in the threshold processing, FAccY_1 is used in order to suppress the influence of noise included in the AccY data. Here, the threshold value TH is a value close to 0 m / s 2 (for example, 0.5 m / s 2 ). The threshold TH may be set stepwise according to the maximum value of FAccY_1. For example, when the maximum value of FAccY_1 is less than 5 m / s 2 , it is set to 0.5 m / s 2, and when the maximum value is 5 m / s 2 or more, it is set to 1.0 m / s 2 .

ステップS25において、上記差分が閾値TH未満でないと判定された場合(ステップS25;NO)、制御部8は、上記変数nを1インクリメントし(ステップS26)、このときの変数nが示すサンプリングタイミングが第三基準タイミング(LMX)であるか否かを判定する(ステップS27)。
ステップS27において、このときの変数nが示すサンプリングタイミングが第三基準タイミング(LMX)でないと判定された場合(ステップS27;NO)、ステップS25へ戻り、制御部8は、ステップS25以降の処理を繰り返し行う。
一方、ステップS27において、このときの変数nが示すサンプリングタイミングが第三基準タイミング(LMX)であると判定された場合(ステップS27;YES)、すなわち、第二基準タイミング(ZC)と第三基準タイミング(LMX)との間において、上記差分が閾値TH未満であるという条件が満たされなかった場合、制御部8は、図9に示すように、当該第三基準タイミング(LMX)を離地タイミングとして設定し(ステップS28)、離地タイミング設定処理を終了する。
When it is determined in step S25 that the difference is not less than the threshold value TH (step S25; NO), the control unit 8 increments the variable n by 1 (step S26), and the sampling timing indicated by the variable n at this time is It is determined whether it is the third reference timing (LMX) (step S27).
When it is determined in step S27 that the sampling timing indicated by the variable n at this time is not the third reference timing (LMX) (step S27; NO), the process returns to step S25, and the control unit 8 performs the processing of step S25 and thereafter. Repeat.
On the other hand, when it is determined in step S27 that the sampling timing indicated by the variable n at this time is the third reference timing (LMX) (step S27; YES), that is, the second reference timing (ZC) and the third reference timing. When the condition that the difference is less than the threshold value TH is not satisfied with respect to the timing (LMX), the control unit 8 sets the third reference timing (LMX) to the departure timing as shown in FIG. 9. (Step S28), and the take-off timing setting process ends.

また、ステップS25において、上記差分が閾値TH未満であると判定された場合(ステップS25;YES)、制御部8は、このときの変数nが示すサンプリングタイミングにおける上下方向の加速度信号AccZの値が所定値C1未満であり、且つ、FAccY_1の値が特定値C2以上であるか否かを判定する(ステップS29)。
ステップS29の判定処理でAccZの値が所定値C1未満であるか否かを判定しているのは、足が地面に設置している間は、AccZの値が正側すなわち上向きになっており、AccZの値が所定値C1以上である場合(例えば、図9の破線の丸で囲った領域参照)は、足はまだ地面に接していると想定されるためである。ここで、所定値C1は、重力加速度9.80665m/sよりも若干大きい値(例えば、12m/s)とする。
また、ステップS29の判定処理でFAccY_1の値が特定値C2以上であるか否かを判定しているのは、FAccY_1の極大値に対して、FAccY_1の値が非常に小さい場合(図10の破線の丸で囲った領域参照)は、そのタイミングにおいてAccYの上昇率が0になっていたとしても、まだ足が地面から離れていないと想定されるためであり、離地タイミングでは、比較的大きな速度の変化(減少)が生じると想定されるためである。ここで、特定値C2は、第二基準タイミング(ZC)から第二分割領域の後端(例えば、図7の線L2参照)までの範囲におけるFAccY_1の値の最大値に所定の係数(例えば、40〜60%)を掛けた値とする。
Further, when it is determined in step S25 that the difference is less than the threshold value TH (step S25; YES), the control unit 8 determines that the value of the vertical acceleration signal AccZ at the sampling timing indicated by the variable n at this time. It is determined whether the value is less than the predetermined value C1 and the value of FAccY_1 is the specific value C2 or more (step S29).
In the determination process of step S29, it is determined whether or not the value of AccZ is less than the predetermined value C1 because the value of AccZ is positive, that is, upward while the foot is installed on the ground. , AccZ is equal to or larger than the predetermined value C1 (for example, refer to a region surrounded by a broken line circle in FIG. 9), it is assumed that the foot is still in contact with the ground. The predetermined value C1 is slightly larger than the gravitational acceleration 9.80665 m / s 2 (e.g., 12m / s 2) to.
Further, in the determination processing of step S29, it is determined whether the value of FAccY_1 is the specific value C2 or more when the value of FAccY_1 is extremely small with respect to the maximum value of FAccY_1 (broken line in FIG. 10). This is because even if the rate of increase of AccY is 0 at that timing, it is assumed that the foot is not yet off the ground, and at the time of takeoff, it is relatively large. This is because it is assumed that the speed will change (decrease). Here, the specific value C2 is a predetermined coefficient (for example, the maximum value of FAccY_1 in the range from the second reference timing (ZC) to the rear end of the second divided region (for example, see the line L2 in FIG. 7). 40 to 60%).

ステップS29において、上記判定条件を満たしていないと判定された場合(ステップS29;NO)、ステップS26へ移行し、制御部8は、それ以降の処理を行う。
一方、ステップS29において、上記判定条件を満たしていると判定された場合(ステップS29;YES)、制御部8は、このときの変数nが示すサンプリングタイミング(図8及び図10参照)を離地タイミングとして設定し(ステップS30)、離地タイミング設定処理を終了する。
When it is determined in step S29 that the above determination conditions are not satisfied (step S29; NO), the process proceeds to step S26, and the control unit 8 performs the subsequent processing.
On the other hand, when it is determined in step S29 that the above determination condition is satisfied (step S29; YES), the control unit 8 departs from the sampling timing (see FIGS. 8 and 10) indicated by the variable n at this time. The timing is set (step S30), and the take-off timing setting process ends.

上述したように、制御部8は、着地タイミング設定処理(ステップS1)、最下点タイミングの設定処理(ステップS2)、離地タイミング設定処理(ステップS3)、及び最高点タイミングの設定処理(ステップS4)を実行することにより、例えば、図11に示すように、一方の足が地面に触れた接地時(着地タイミング)、この足の接地中に腰が最も低い位置を通過した最下点時(最下点タイミング)、この足が地面から離れた離地時(離地タイミング)、及び、この足が地面から離れた後、腰が最も高い位置を通過した最高点時(最高点タイミング)、並びに、逆の足が地面に触れた逆足接地時(着地タイミング)、この足の接地中に腰が最も低い位置を通過した逆足最下点時(最下点タイミング)、この足が地面から離れた逆足離地時(離地タイミング)、及び、この足が地面から離れた後、腰が最も高い位置を通過した逆足最高点時(最高点タイミング)を検出することができる。   As described above, the control unit 8 sets the landing timing setting process (step S1), the lowest point timing setting process (step S2), the departure timing setting process (step S3), and the highest point timing setting process (step S3). By executing S4), for example, as shown in FIG. 11, when one foot touches the ground and touches down (landing timing), the lowest point when the waist passes the lowest position while touching the foot (Bottom point timing), when this foot leaves the ground (takeoff timing), and when this foot has left the ground and is at the highest point when the waist has passed the highest position (highest point timing) , And, when the reverse foot touches the ground, the reverse foot touches the ground (landing timing), and when the hip passes the lowest position while touching the foot, the reverse foot at the lowest point (lowest point timing) Reverse foot off the ground When (away land timing), and, after the foot leaves the ground, it is possible to detect the time Gyakuashi highest point has passed the highest position waist (highest point timing).

次に、本実施形態の測定装置1による接地時間の算出の精度を検証するため、本実施形態による接地時間の算出値と、フォースプレートにより測定された接地時間、及びオプトジャンプ(登録商標)光学測定装置により測定された接地時間とを比較した結果を、図12及び図13を用いて説明する。ここで、フォースプレートとは、歩行や走行など動作の力学的計測を行う装置である。また、オプトジャンプ(登録商法)光学測定装置(以下、光学測定装置と称する)とは、走行路に沿って一定間隔で設けられた、走行路に直交する方向に光を照射する複数の送信部と、送信部から照射された光を受ける複数の受信部と、を有して、当該光の遮断(接地)と接続(離地)を感知する装置である。   Next, in order to verify the accuracy of the calculation of the ground contact time by the measuring device 1 of the present embodiment, the calculated value of the ground contact time according to the present embodiment, the ground contact time measured by the force plate, and Optojump (registered trademark) optics. The result of comparison with the contact time measured by the measuring device will be described with reference to FIGS. 12 and 13. Here, the force plate is a device that mechanically measures movements such as walking and running. Further, the opt-jump (registered trademark) optical measuring device (hereinafter, referred to as an optical measuring device) is a plurality of transmitting units which are provided at regular intervals along the traveling road and emit light in a direction orthogonal to the traveling road. And a plurality of receiving units for receiving the light emitted from the transmitting unit, and detecting the interruption (grounding) and the connection (departure) of the light.

図12は、フォースプレートにより測定された接地時間と本実施形態の測定装置1により算出された接地時間の散布図を示している。なお、図中の[Cor]は相関係数、[MeanErr]は平均誤差、[MAE]は平均絶対誤差、[ErrStd]は誤差の標準偏差を表している。
図12に示すように、接地時間の相関係数[Cor]は0.896であり、本実施形態の測定装置1により算出された接地時間において高い相関が認められた。
また、フォースプレートにより測定された接地時間と、本実施形態の測定装置1とは離地タイミングの設定方法が異なる従来の測定装置により算出された接地時間と、の相関係数は0.856であり、本実施形態の測定装置1により算出された接地時間は、従来に比べて、より高い相関が得られて測定精度が向上しているといえる。
ここで、従来の測定装置は、着地タイミングの設定方法に関しては、本実施形態の測定装置1と同じ方法を採用している。一方、離地タイミングに関しては、第二分割領域(図7の線L1,L2で分割された領域)のうちの後半部分で、前後方向の加速度信号AccYの値が負から正、つまり加速から減速へと変化した後の前後方向の加速度信号AccYの波形における最初の正の極大値を探索し、当該最初の正の極大値が存在した場合には、その位置(時間)を離地タイミングとする一方、最初の正の極大値が存在しなかった場合には、第二分割領域の92〜97%の範囲に収まる所定地点を離地タイミングとしている。
つまり、本実施形態の測定装置1と従来の測定装置とでは、接地時間を算出する際、離地タイミングの値だけが異なっていることから、接地時間の相関が従来に比べて高い本実施形態の測定装置1では離地タイミングの検出の精度が向上しているといえる。
FIG. 12 shows a scatter diagram of the contact time measured by the force plate and the contact time calculated by the measuring apparatus 1 of the present embodiment. In the figure, [Cor] represents a correlation coefficient, [MeanErr] represents an average error, [MAE] represents an average absolute error, and [ErrStd] represents a standard deviation of the error.
As shown in FIG. 12, the correlation coefficient [Cor] of the contact time was 0.896, and a high correlation was recognized in the contact time calculated by the measuring apparatus 1 of the present embodiment.
Further, the correlation coefficient between the ground contact time measured by the force plate and the ground contact time calculated by the conventional measuring device having a method of setting the departure timing different from that of the measuring device 1 of the present embodiment is 0.856. Therefore, it can be said that the ground contact time calculated by the measuring apparatus 1 of the present embodiment has a higher correlation than that in the conventional case and the measurement accuracy is improved.
Here, the conventional measuring apparatus employs the same method as the measuring apparatus 1 of the present embodiment with respect to the landing timing setting method. On the other hand, regarding the take-off timing, the value of the acceleration signal AccY in the front-rear direction is negative to positive, that is, from acceleration to deceleration, in the latter half of the second divided area (area divided by lines L1 and L2 in FIG. 7). The first positive maximum value in the waveform of the acceleration signal AccY in the front-rear direction after the change to is searched for, and if the first positive maximum value exists, the position (time) is set as the departure timing. On the other hand, when the first positive maximum value does not exist, the predetermined point within the range of 92 to 97% of the second divided area is set as the departure timing.
That is, since the measuring apparatus 1 of the present embodiment and the conventional measuring apparatus differ only in the value of the takeoff timing when calculating the ground contact time, the correlation of the ground contact time is higher than that in the conventional embodiment. It can be said that the accuracy of detection of the departure timing is improved in the measuring apparatus 1 of 1.

図13は、光学測定装置により測定された接地時間と本実施形態の測定装置1により算出された接地時間の散布図を示している。
図13に示すように、接地時間の相関係数[Cor]は0.892であり、本実施形態の測定装置1により算出された接地時間において高い相関が認められた。
また、光学測定装置により算出された接地時間と、上述の従来の測定装置により算出された接地時間と、の相関係数は0.756であり、この場合においても、本実施形態の測定装置1により算出された接地時間は、従来に比べて、より高い相関が得られて測定精度が向上しているといえる。
つまり、図13の散布図から導出された相関係数と、上述の従来の測定装置における相関係数と、の関係からしても接地時間の相関が従来に比べて高い本実施形態の測定装置1では離地タイミングの検出の精度が向上しているといえる。
FIG. 13 is a scatter diagram of the contact time measured by the optical measuring device and the contact time calculated by the measuring device 1 of the present embodiment.
As shown in FIG. 13, the correlation coefficient [Cor] of the contact time is 0.892, and a high correlation was recognized in the contact time calculated by the measuring apparatus 1 of the present embodiment.
Further, the correlation coefficient between the ground contact time calculated by the optical measuring device and the ground contact time calculated by the above-mentioned conventional measuring device is 0.756, and even in this case, the measuring device 1 of the present embodiment is also used. It can be said that the contact time calculated by the above has a higher correlation than that in the conventional case and the measurement accuracy is improved.
That is, even from the relationship between the correlation coefficient derived from the scatter diagram of FIG. 13 and the correlation coefficient in the above-described conventional measurement apparatus, the correlation of the ground contact time is higher than that of the conventional measurement apparatus. In No. 1, it can be said that the accuracy of detection of the departure timing is improved.

以上のように、本実施形態によれば、加速度センサ4により得られた前後方向、上下方向それぞれの加速度信号の波形に基づいて離地タイミングが求められているので、走行時においても離地に関するデータの一つである離地タイミングを精度良く測定することができる。   As described above, according to the present embodiment, the take-off timing is determined based on the waveforms of the acceleration signals in the front-rear direction and the vertical direction obtained by the acceleration sensor 4, so that the take-off timing can be related to the take-off even during traveling. It is possible to accurately measure the take-off timing, which is one of the data.

また、離地タイミング設定処理におけるステップS25及びステップS29の判定条件を満たすサンプルが存在しなかった場合には、第三基準タイミング(LMX)を離地タイミングとするので、かかる場合でも離地タイミングを推測して設定することができる。   Further, when there is no sample that satisfies the determination conditions of step S25 and step S29 in the takeoff timing setting process, the third reference timing (LMX) is set as the takeoff timing. It can be guessed and set.

また、前後方向の加速度信号AccYに対して、タップ数を15とする三角移動平均による平滑化処理を行い、当該処理後のFAccY_2のデータを用いることによって、プラスゼロクロスタイミング(ZC)を探索するに際して、AccYやFAccY_1のデータから検出される離地前の瞬間的なプラスゼロクロスタイミングを除外することができ、適正なプラスゼロクロスタイミング(ZC)を探索することができるようになる。   In addition, when the acceleration signal AccY in the front-rear direction is subjected to smoothing processing by a triangular moving average with 15 taps and the data of FAccY_2 after the processing is used, a plus-zero cross timing (ZC) is searched for. , AccY and FAccY_1, the instantaneous plus-zero cross timing before takeoff detected from the data can be excluded, and an appropriate plus-zero cross timing (ZC) can be searched for.

また、前後方向の加速度信号AccYに対して、タップ数を3とする三角移動平均による平滑化処理を行い、当該処理後のFAccY_1のデータを用いることによって、離地タイミング設定処理のステップS25の閾値処理において、例えば、ノイズや体の固有の動き等の影響によって、第二基準タイミングの後、極大値が一気に現れるタイミングや、途中で変曲点が見られる波形の後に極大値が現れるタイミングを離地タイミングの候補から除外することができる。   Further, the acceleration signal AccY in the front-rear direction is subjected to smoothing processing by a triangular moving average with the number of taps set to 3, and the data of FAccY_1 after the processing is used to determine the threshold value in step S25 of the takeoff timing setting processing. In the processing, for example, due to the influence of noise or peculiar motion of the body, the timing at which the maximum value appears suddenly after the second reference timing, or the timing at which the maximum value appears after the waveform with an inflection point in the middle, is separated. Can be excluded from ground timing candidates.

また、加速度センサ4により得られた前後方向、上下方向それぞれの加速度信号の波形に基づいて、足が地面に触れた着地タイミング、足の着地中に腰が最も低い位置を通過した最下点タイミング、及び、足が地面から離れた後、腰が最も高い位置を通過した最高点タイミングも求められているので、ランニング動作におけるサイクルの各曲面を好適に検出することができる。   Also, based on the waveforms of the acceleration signals in the front-back direction and the vertical direction obtained by the acceleration sensor 4, the landing timing when the foot touches the ground, the lowest point timing when the hip passes the lowest position during landing of the foot. Also, since the highest point timing at which the waist has passed the highest position after the foot has left the ground is obtained, it is possible to preferably detect each curved surface of the cycle in the running motion.

なお、以上本発明の実施形態について説明したが、本発明は、かかる実施形態に限定されず、その要旨を逸脱しない範囲で、種々変形が可能であることは言うまでもない。   Although the embodiments of the present invention have been described above, it is needless to say that the present invention is not limited to the embodiments and various modifications can be made without departing from the scope of the invention.

例えば、上記実施形態では、予め取得された前後方向、左右方向、上下方向それぞれの加速度信号の波形に対して測定処理が施される場合を例示して説明したが、加速度信号を取得しながらリアルタイムで、あるいはほぼリアルタイムで測定処理を実行するようにしてもよい。   For example, in the above-described embodiment, the case where the measurement process is performed on the waveforms of the acceleration signals in the front-rear direction, the left-right direction, and the up-down direction that have been acquired in advance has been described. Alternatively, the measurement process may be executed in real time.

また、上記実施形態においては、制御部8が加速度センサ4とともに一体化されている場合を例示して説明したが、制御部は加速度センサ4と別体であってもよい。具体的には、加速度センサ4が取得した前後方向、左右方向、上下方向のそれぞれの加速度信号を外部の制御部に出力し、当該加速度信号及びその波形に基づいて算出部が着地、離地等に関するデータを算出するようにしてもよい。外部の制御部としては、例えばパソコンや携帯電話、スマートフォン、タブレット機器、リストバンド型端末などの情報端末が挙げられる。   In the above embodiment, the case where the control unit 8 is integrated with the acceleration sensor 4 has been described as an example, but the control unit may be separate from the acceleration sensor 4. Specifically, the acceleration sensor 4 outputs the acceleration signals in the front-rear direction, the left-right direction, and the up-down direction to an external control unit, and the calculation unit lands and takes off based on the acceleration signal and its waveform. You may make it calculate the data regarding. Examples of the external control unit include an information terminal such as a personal computer, a mobile phone, a smartphone, a tablet device, and a wristband type terminal.

また、上記実施形態では、人体(生体、ユーザ)の腰の角速度を所定時間毎に測定して角速度信号を出力する角速度センサを更に備え、制御部8は、角速度信号を積分することで、腰の回転量を表す回転角度の波形を求め、回転角度の波形に基づき、上述のランニング動作におけるサイクルの各局面の腰の回転量を導出するようにしても良い。
これにより、ランニング動作におけるサイクルの各局面に照らし、腰の回転量の変化を比較評価することができるようになる。
Further, in the above-described embodiment, the control unit 8 further includes an angular velocity sensor that measures the angular velocity of the waist of a human body (living body, user) at predetermined time intervals and outputs an angular velocity signal. It is also possible to obtain a rotation angle waveform that represents the rotation amount and to derive the waist rotation amount for each phase of the cycle in the running motion described above based on the rotation angle waveform.
This makes it possible to compare and evaluate changes in the amount of rotation of the waist in the light of each phase of the cycle in the running motion.

以上、本発明の実施形態を説明したが、本発明の範囲は、上述の実施の形態に限定するものではなく、特許請求の範囲に記載された発明の範囲をその均等の範囲を含む。
以下に、この出願の願書に最初に添付した特許請求の範囲に記載した発明を付記する。付記に記載した請求項の項番は、この出願の願書に最初に添付した特許請求の範囲の通りである。
Although the embodiments of the present invention have been described above, the scope of the present invention is not limited to the above-mentioned embodiments, and includes the scope of the invention described in the claims and its equivalent scope.
Hereinafter, the inventions described in the claims attached to the application of this application will be additionally described. The claim numbers listed in the appendices are as set forth in the claims initially attached to the application for this application.

〔付記〕
<請求項1>
生体に装着され、前記生体の動きに対応する加速度を検出して前記加速度に対応する加速度信号を出力する加速度センサと、
前記生体が進行方向に移動しているときに前記加速度信号に基づいて、前記生体の足の動きに関する少なくとも一つの動作タイミングを取得する制御部と、
を備え、
前記制御部は、
前記加速度信号に基づいて、前記進行方向の加速度信号を第1加速度信号として算出するとともに、上下方向の加速度信号を第2加速度信号として算出し、
前記第1加速度信号の第1期間の移動平均を用いて前記第1加速度信号を平滑化した第1信号と、前記第1加速度信号の前記第1期間より短い第2期間の移動平均を用いて前記第1加速度信号を平滑化した第2信号と、前記第2加速度信号を二度積分した第3信号と、を算出し、
前記第1信号、前記第2信号及び前記第3信号に基づいて、前記生体の足が離地したタイミングと推定される離地タイミングを前記動作タイミングとして取得することを特徴とする測定装置。
<請求項2>
前記第1期間は25〜35ms、前記第2期間は5〜10msに設定されていることを特徴とする請求項1に記載の測定装置。
<請求項3>
前記制御部は、前記加速度信号をサンプリングタイミング毎に取り込み、
前記第1信号における前記移動平均は、前記第1期間における複数の前記サンプリングタイミングでの前記加速度信号の値による三角移動平均であり、
前記第2信号における前記移動平均は、前記第2期間における複数の前記サンプリングタイミングでの前記加速度信号の値による三角移動平均であることを特徴とする請求項1又は2に記載の測定装置。
<請求項4>
前記制御部は、
前記第3信号の時間軸に対する波形が極大値を示す複数の第一基準タイミングを抽出し、
隣接する2つの前記第一基準タイミング間の後半の期間内において、前記第1信号の値が前記進行方向の値から前記進行方向の逆方向の値となった第二基準タイミングと、前記第二基準タイミングより時間的に後で、前記第1信号の時間軸に対する波形が前記進行方向の逆方向で極大値を示す第三基準タイミングと、を抽出し、
前記第二基準タイミングと前記第三基準タイミングとの間を第1探索期間として、前記第2信号の値と、前記第2信号の値と前記第2加速度信号の値との相互関係と、に基づいて前記生体の足が離地したタイミングと推定される第1タイミングを前記第1探索期間において探索し、前記探索により前記第1タイミングが求められたとき、前記第1タイミングを前記離地タイミングとして取得することを特徴とする請求項1乃至3の何れか一項に記載の測定装置。
<請求項5>
前記制御部は、
前記第1探索期間に含まれる複数の前記サンプリングタイミングにおいて、第一条件、第二条件及び第三条件の全てを満たす特定のサンプリングタイミングを前記第1タイミングとして探索し、前記探索により求められた前記特定のサンプリングタイミングを前記離地タイミングとし、
前記第一条件は、前記特定のサンプリングタイミングでの前記第2信号の値と、当該特定のサンプリングタイミングから1サンプリングタイミングだけ前のサンプリングタイミングでの前記第2信号の値との差分が閾値より小さい値であることであり、前記第二条件は、前記特定のサンプリングタイミングでの前記第2加速度信号の値が所定値より小さい値であることであり、前記第三条件は、前記特定のサンプリングタイミングでの前記第2信号の値が、前記第1探索期間内での前記第2信号の値の最大値に所定の係数を掛けた特定値より大きい値であることであることを特徴とする請求項4に記載の測定装置。
<請求項6>
前記閾値は、0.5〜1.0m/sに設定され、前記所定値は、重力加速度より大きい値に設定され、前記所定の係数は、0.4〜0.6の値に設定されていることを特徴とする請求項5に記載の測定装置。
<請求項7>
前記制御部は、前記第1探索期間における前記探索おいて、前記第1タイミングが求められなかったとき、前記第三基準タイミングを前記離地タイミングとして取得することを特徴とする請求項4乃至6の何れか一項に記載の測定装置。
<請求項8>
前記制御部は、
前記第2加速度信号に対し平滑化処理を行った第4信号を算出し、
前記第4信号の時間軸に対する波形が極大値を示す複数の第四基準タイミングを抽出し、
前記第一基準タイミングと、前記第一基準タイミングより時間的に後で、前記第一基準タイミングに最も近い前記第四タイミングとの間を第2探索期間とし、
前記第2探索期間において、前記第1加速度信号の時間軸に対する波形における前記進行方向に対する減速動作の極値のタイミングを探索し、前記探索により求められた前記極値のタイミングに基づく第2タイミングを、前記動作タイミングにおける、前記足が着地したタイミングと推定される着地タイミングとして取得し、
前記第3信号の前記波形が極小値を示すタイミングを、前記動作タイミングにおける、前記足の着地中に前記生体の腰が最も低い位置となった最下点タイミングとして取得し、
前記第一基準タイミングを、前記動作タイミングにおける、前記足が地面から離れた後に前記腰が最も高い位置となった最高点タイミングとして取得することを特徴とする請求項1乃至7の何れか一項に記載の測定装置。
<請求項9>
測定装置における測定方法であって、
前記測定装置は、生体に装着され、前記生体の動きに対応する加速度を検出して前記加速度に対応する加速度信号を出力する加速度センサを有し、
前記加速度信号に基づいて、前記進行方向の加速度信号を第1加速度信号として算出するとともに、上下方向の加速度信号を第2加速度信号として算出し、
前記第1加速度信号の第1期間の移動平均を用いて前記第1加速度信号を平滑化した第1信号と、前記第1加速度信号の前記第1期間より短い第2期間の移動平均を用いて前記第1加速度信号を平滑化した第2信号と、前記第2加速度信号を二度積分した第3信号と、を算出し、
前記第1信号、前記第2信号及び前記第3信号に基づいて、前記生体の足が離地したタイミングと推定される離地タイミングを前記動作タイミングとして取得することを特徴とする測定方法。
<請求項10>
測定装置における測定プログラムであって、
前記測定装置は、生体に装着され、前記生体の動きに対応する加速度を検出して前記加速度に対応加速度信号を出力する加速度センサを有し、
前記加速度信号に基づいて、前記進行方向の加速度信号を第1加速度信号として算出させるとともに、上下方向の加速度信号を第2加速度信号として算出させ、
前記第1加速度信号の第1期間の移動平均を用いて前記第1加速度信号を平滑化した第1信号と、前記第1加速度信号の前記第1期間より短い第2期間の移動平均を用いて前記第1加速度信号を平滑化した第2信号と、前記第2加速度信号を二度積分した第3信号と、を算出させ、
前記第1信号、前記第2信号及び前記第3信号に基づいて、前記生体の足が離地したタイミングと推定される離地タイミングを前記動作タイミングとして取得させることを特徴とする測定プログラム。
[Appendix]
<Claim 1>
An acceleration sensor that is attached to a living body, detects an acceleration corresponding to the movement of the living body, and outputs an acceleration signal corresponding to the acceleration;
Based on the acceleration signal when the living body is moving in the traveling direction, a control unit that acquires at least one operation timing regarding the movement of the foot of the living body,
Equipped with
The control unit is
Based on the acceleration signal, the acceleration signal in the traveling direction is calculated as a first acceleration signal, and the acceleration signal in the vertical direction is calculated as a second acceleration signal,
A first signal obtained by smoothing the first acceleration signal using a moving average of the first acceleration signal in a first period, and a moving average of the first acceleration signal in a second period shorter than the first period are used. Calculating a second signal obtained by smoothing the first acceleration signal and a third signal obtained by integrating the second acceleration signal twice,
A take-off timing estimated to be a take-off timing of the foot of the living body is acquired as the operation timing based on the first signal, the second signal, and the third signal.
<Claim 2>
The measuring device according to claim 1, wherein the first period is set to 25 to 35 ms and the second period is set to 5 to 10 ms.
<Claim 3>
The control unit fetches the acceleration signal at each sampling timing,
The moving average in the first signal is a triangular moving average based on the values of the acceleration signals at the plurality of sampling timings in the first period,
The measuring apparatus according to claim 1, wherein the moving average in the second signal is a triangular moving average based on values of the acceleration signal at a plurality of the sampling timings in the second period.
<Claim 4>
The control unit is
Extracting a plurality of first reference timings at which the waveform of the third signal with respect to the time axis exhibits a maximum value,
A second reference timing at which the value of the first signal becomes a value in the opposite direction of the traveling direction from a value in the traveling direction in a latter half period between two adjacent first reference timings; And a third reference timing in which the waveform of the first signal with respect to the time axis shows a maximum value in a direction opposite to the traveling direction, temporally after the reference timing,
The first search period is between the second reference timing and the third reference timing, and the value of the second signal and the correlation between the value of the second signal and the value of the second acceleration signal are set. Based on the first timing, the first timing estimated to be the timing at which the foot of the living body has departed is searched based on the first search period, and when the first timing is obtained by the search, the first timing is set to the departing timing. The measuring device according to claim 1, wherein the measuring device is obtained as
<Claim 5>
The control unit is
Among the plurality of sampling timings included in the first search period, a specific sampling timing that satisfies all of the first condition, the second condition, and the third condition is searched as the first timing, and the specific timing obtained is obtained by the search. A specific sampling timing is the takeoff timing,
The first condition is that a difference between a value of the second signal at the specific sampling timing and a value of the second signal at a sampling timing one sampling timing before the specific sampling timing is smaller than a threshold value. Is a value, the second condition is that the value of the second acceleration signal at the specific sampling timing is smaller than a predetermined value, and the third condition is the specific sampling timing. The value of the second signal in the above is larger than a specific value obtained by multiplying the maximum value of the values of the second signal in the first search period by a predetermined coefficient. Item 4. The measuring device according to item 4.
<Claim 6>
The threshold value is set to 0.5 to 1.0 m / s 2 , the predetermined value is set to a value larger than the gravitational acceleration, and the predetermined coefficient is set to a value of 0.4 to 0.6. The measuring device according to claim 5, wherein
<Claim 7>
7. The control unit acquires the third reference timing as the departure timing when the first timing is not obtained in the search in the first search period. The measuring device according to any one of 1.
<Claim 8>
The control unit is
Calculating a fourth signal obtained by performing a smoothing process on the second acceleration signal,
Extracting a plurality of fourth reference timings in which the waveform of the fourth signal with respect to the time axis exhibits a maximum value,
A second search period between the first reference timing and the fourth timing that is closest to the first reference timing after the first reference timing.
In the second search period, the timing of the extreme value of the deceleration operation with respect to the traveling direction in the waveform of the first acceleration signal with respect to the time axis is searched for, and the second timing based on the timing of the extreme value obtained by the search is set. , As the landing timing estimated to be the landing timing of the foot in the operation timing,
The timing at which the waveform of the third signal shows a minimum value is acquired as the lowest point timing at which the waist of the living body is at the lowest position during landing of the foot at the operation timing,
8. The first reference timing is acquired as a highest point timing at which the waist is at a highest position after the foot has left the ground in the operation timing. The measuring device according to 1.
<Claim 9>
A measuring method in a measuring device,
The measuring device includes an acceleration sensor that is attached to a living body, detects an acceleration corresponding to the movement of the living body, and outputs an acceleration signal corresponding to the acceleration.
Based on the acceleration signal, the acceleration signal in the traveling direction is calculated as a first acceleration signal, and the acceleration signal in the vertical direction is calculated as a second acceleration signal,
Using a first signal obtained by smoothing the first acceleration signal using a moving average of the first acceleration signal in a first period, and a moving average of the first acceleration signal in a second period shorter than the first period. Calculating a second signal obtained by smoothing the first acceleration signal and a third signal obtained by integrating the second acceleration signal twice,
A measurement method, characterized in that, based on the first signal, the second signal, and the third signal, a takeoff timing estimated to be a timing at which the foot of the living body has taken off is acquired as the operation timing.
<Claim 10>
A measuring program in a measuring device,
The measuring device is attached to a living body, and has an acceleration sensor that detects an acceleration corresponding to the movement of the living body and outputs an acceleration signal corresponding to the acceleration,
Based on the acceleration signal, the acceleration signal in the traveling direction is calculated as a first acceleration signal, and the acceleration signal in the vertical direction is calculated as a second acceleration signal,
Using a first signal obtained by smoothing the first acceleration signal using a moving average of the first acceleration signal in a first period, and a moving average of the first acceleration signal in a second period shorter than the first period. Calculating a second signal obtained by smoothing the first acceleration signal and a third signal obtained by integrating the second acceleration signal twice,
A measurement program, characterized in that, based on the first signal, the second signal, and the third signal, a take-off timing estimated to be a timing at which the foot of the living body is taken off is acquired as the operation timing.

1 測定装置
2 本体部
3 ベルト部
4 加速度センサ
5 通信部
6 表示部
7 操作部
8 制御部
1 Measuring Device 2 Main Body 3 Belt 4 Acceleration Sensor 5 Communication 6 Display 7 Operation 8 Control

Claims (10)

生体に装着され、前記生体の動きに対応する加速度を検出して前記加速度に対応する加速度信号を出力する加速度センサと、
前記生体が進行方向に移動しているときに前記加速度信号に基づいて、前記生体の足の動きに関する少なくとも一つの動作タイミングを取得する制御部と、
を備え、
前記制御部は、
前記加速度信号に基づいて、前記進行方向の加速度信号を第1加速度信号として算出するとともに、上下方向の加速度信号を第2加速度信号として算出し、
前記第1加速度信号の第1期間の移動平均を用いて前記第1加速度信号を平滑化した第1信号と、前記第1加速度信号の前記第1期間より短い第2期間の移動平均を用いて前記第1加速度信号を平滑化した第2信号と、前記第2加速度信号を二度積分した第3信号と、を算出し、
前記第1信号、前記第2信号及び前記第3信号に基づいて、前記生体の足が離地したタイミングと推定される離地タイミングを前記動作タイミングとして取得することを特徴とする測定装置。
An acceleration sensor that is attached to a living body, detects an acceleration corresponding to the movement of the living body, and outputs an acceleration signal corresponding to the acceleration;
Based on the acceleration signal when the living body is moving in the traveling direction, a control unit that acquires at least one operation timing regarding the movement of the foot of the living body,
Equipped with
The control unit is
Based on the acceleration signal, the acceleration signal in the traveling direction is calculated as a first acceleration signal, and the acceleration signal in the vertical direction is calculated as a second acceleration signal,
A first signal obtained by smoothing the first acceleration signal using a moving average of the first acceleration signal in a first period, and a moving average of the first acceleration signal in a second period shorter than the first period are used. Calculating a second signal obtained by smoothing the first acceleration signal and a third signal obtained by integrating the second acceleration signal twice,
A take-off timing estimated to be a take-off timing of the foot of the living body is acquired as the operation timing based on the first signal, the second signal, and the third signal.
前記第1期間は25〜35ms、前記第2期間は5〜10msに設定されていることを特徴とする請求項1に記載の測定装置。   The measuring device according to claim 1, wherein the first period is set to 25 to 35 ms and the second period is set to 5 to 10 ms. 前記制御部は、前記加速度信号をサンプリングタイミング毎に取り込み、
前記第1信号における前記移動平均は、前記第1期間における複数の前記サンプリングタイミングでの前記加速度信号の値による三角移動平均であり、
前記第2信号における前記移動平均は、前記第2期間における複数の前記サンプリングタイミングでの前記加速度信号の値による三角移動平均であることを特徴とする請求項1又は2に記載の測定装置。
The control unit fetches the acceleration signal at each sampling timing,
The moving average in the first signal is a triangular moving average based on the values of the acceleration signals at the plurality of sampling timings in the first period,
The measuring apparatus according to claim 1, wherein the moving average in the second signal is a triangular moving average based on values of the acceleration signal at a plurality of the sampling timings in the second period.
前記制御部は、
前記第3信号の時間軸に対する波形が極大値を示す複数の第一基準タイミングを抽出し、
隣接する2つの前記第一基準タイミング間の後半の期間内において、前記第1信号の値が前記進行方向の値から前記進行方向の逆方向の値となった第二基準タイミングと、前記第二基準タイミングより時間的に後で、前記第1信号の時間軸に対する波形が前記進行方向の逆方向で極大値を示す第三基準タイミングと、を抽出し、
前記第二基準タイミングと前記第三基準タイミングとの間を第1探索期間として、前記第2信号の値と、前記第2信号の値と前記第2加速度信号の値との相互関係と、に基づいて前記生体の足が離地したタイミングと推定される第1タイミングを前記第1探索期間において探索し、前記探索により前記第1タイミングが求められたとき、前記第1タイミングを前記離地タイミングとして取得することを特徴とする請求項1乃至3の何れか一項に記載の測定装置。
The control unit is
Extracting a plurality of first reference timings at which the waveform of the third signal with respect to the time axis exhibits a maximum value,
A second reference timing at which the value of the first signal becomes a value in the opposite direction of the traveling direction from a value in the traveling direction in a latter half period between two adjacent first reference timings; And a third reference timing in which the waveform of the first signal with respect to the time axis shows a maximum value in a direction opposite to the traveling direction, temporally after the reference timing,
The first search period is between the second reference timing and the third reference timing, and the value of the second signal and the correlation between the value of the second signal and the value of the second acceleration signal are set. Based on the first timing, the first timing estimated to be the timing at which the foot of the living body has departed is searched based on the first search period, and when the first timing is obtained by the search, the first timing is set to the departing timing. The measuring device according to claim 1, wherein the measuring device is obtained as
前記制御部は、
前記第3信号の時間軸に対する波形が極大値を示す複数の第一基準タイミングを抽出し、
隣接する2つの前記第一基準タイミング間の後半の期間内において、前記第1信号の値が前記進行方向の値から前記進行方向の逆方向の値となった第二基準タイミングと、前記第二基準タイミングより時間的に後で、前記第1信号の時間軸に対する波形が前記進行方向の逆方向で極大値を示す第三基準タイミングと、を抽出し、
前記第二基準タイミングと前記第三基準タイミングとの間を第1探索期間として、前記第2信号の値と、前記第2信号の値と前記第2加速度信号の値との相互関係と、に基づいて前記生体の足が離地したタイミングと推定される第1タイミングを前記第1探索期間において探索し、前記探索により前記第1タイミングが求められたとき、前記第1タイミングを前記離地タイミングとして取得し、
前記第1探索期間に含まれる複数の前記サンプリングタイミングにおいて、第一条件、第二条件及び第三条件の全てを満たす特定のサンプリングタイミングを前記第1タイミングとして探索し、前記探索により求められた前記特定のサンプリングタイミングを前記離地タイミングとし、
前記第一条件は、前記特定のサンプリングタイミングでの前記第2信号の値と、当該特定のサンプリングタイミングから1サンプリングタイミングだけ前のサンプリングタイミングでの前記第2信号の値との差分が閾値より小さい値であることであり、前記第二条件は、前記特定のサンプリングタイミングでの前記第2加速度信号の値が所定値より小さい値であることであり、前記第三条件は、前記特定のサンプリングタイミングでの前記第2信号の値が、前記第1探索期間内での前記第2信号の値の最大値に所定の係数を掛けた特定値より大きい値であることであることを特徴とする請求項に記載の測定装置。
The control unit is
Extracting a plurality of first reference timings at which the waveform of the third signal with respect to the time axis exhibits a maximum value,
A second reference timing at which the value of the first signal becomes a value in the opposite direction of the traveling direction from a value in the traveling direction in a latter half period between two adjacent first reference timings; And a third reference timing in which the waveform of the first signal with respect to the time axis shows a maximum value in a direction opposite to the traveling direction, temporally after the reference timing,
The first search period is between the second reference timing and the third reference timing, and the value of the second signal and the correlation between the value of the second signal and the value of the second acceleration signal are set. Based on the first timing, the first timing estimated to be the timing at which the foot of the living body has departed is searched based on the first search period, and when the first timing is obtained by the search, the first timing is set to the departing timing. Get as
Among the plurality of sampling timings included in the first search period, a specific sampling timing that satisfies all of the first condition, the second condition, and the third condition is searched as the first timing, and the specific timing obtained is obtained by the search. A specific sampling timing is the takeoff timing,
The first condition is that a difference between a value of the second signal at the specific sampling timing and a value of the second signal at a sampling timing one sampling timing before the specific sampling timing is smaller than a threshold value. Is a value, the second condition is that the value of the second acceleration signal at the specific sampling timing is smaller than a predetermined value, and the third condition is the specific sampling timing. The value of the second signal in the above is larger than a specific value obtained by multiplying the maximum value of the values of the second signal in the first search period by a predetermined coefficient. Item 3. The measuring device according to Item 3 .
前記閾値は、0.5〜1.0m/sに設定され、前記所定値は、重力加速度より大きい値に設定され、前記所定の係数は、0.4〜0.6の値に設定されていることを特徴とする請求項5に記載の測定装置。 The threshold value is set to 0.5 to 1.0 m / s 2 , the predetermined value is set to a value larger than the gravitational acceleration, and the predetermined coefficient is set to a value of 0.4 to 0.6. The measuring device according to claim 5, wherein 前記制御部は、前記第1探索期間における前記探索おいて、前記第1タイミングが求められなかったとき、前記第三基準タイミングを前記離地タイミングとして取得することを特徴とする請求項4乃至6の何れか一項に記載の測定装置。   7. The control unit acquires the third reference timing as the departure timing when the first timing is not obtained in the search in the first search period. The measuring device according to any one of 1. 前記制御部は、
前記第2加速度信号に対し平滑化処理を行った第4信号を算出し、
前記第4信号の時間軸に対する波形が極大値を示す複数の第四基準タイミングを抽出し、
前記第一基準タイミングと、前記第一基準タイミングより時間的に後で、前記第一基準タイミングに最も近い前記第四基準タイミングとの間を第2探索期間とし、
前記第2探索期間において、前記第1加速度信号の時間軸に対する波形における前記進行方向に対する減速動作の極値のタイミングを探索し、前記探索により求められた前記極値のタイミングに基づく第2タイミングを、前記動作タイミングにおける、前記足が着地したタイミングと推定される着地タイミングとして取得し、
前記第3信号の前記波形が極小値を示すタイミングを、前記動作タイミングにおける、前記足の着地中に前記生体の腰が最も低い位置となった最下点タイミングとして取得し、
前記第一基準タイミングを、前記動作タイミングにおける、前記足が地面から離れた後に前記腰が最も高い位置となった最高点タイミングとして取得することを特徴とする請求項乃至7の何れか一項に記載の測定装置。
The control unit is
Calculating a fourth signal obtained by performing a smoothing process on the second acceleration signal,
Extracting a plurality of fourth reference timings in which the waveform of the fourth signal with respect to the time axis exhibits a maximum value,
A second search period between the first reference timing and the fourth reference timing that is closest to the first reference timing after the first reference timing.
In the second search period, the timing of the extreme value of the deceleration operation with respect to the traveling direction in the waveform of the first acceleration signal with respect to the time axis is searched for, and the second timing based on the timing of the extreme value obtained by the search is set. , As the landing timing estimated to be the landing timing of the foot in the operation timing,
The timing at which the waveform of the third signal shows a minimum value is acquired as the lowest point timing at which the waist of the living body is at the lowest position during landing of the foot at the operation timing,
The first reference timing, in the operation timing, any one of claims 4 to 7 wherein the foot and wherein said that the waist is obtained as the highest point timing became highest position after leaving the ground The measuring device according to 1.
測定装置における測定方法であって、
前記測定装置は、生体に装着され、前記生体の動きに対応する加速度を検出して前記加速度に対応する加速度信号を出力する加速度センサと、前記生体が進行方向に移動しているときに前記加速度信号に基づいて、前記生体の足の動きに関する少なくとも一つの動作タイミングを取得する制御部と、を有し、
前記加速度信号に基づいて、前記進行方向の加速度信号を第1加速度信号として算出するとともに、上下方向の加速度信号を第2加速度信号として算出し、
前記第1加速度信号の第1期間の移動平均を用いて前記第1加速度信号を平滑化した第1信号と、前記第1加速度信号の前記第1期間より短い第2期間の移動平均を用いて前記第1加速度信号を平滑化した第2信号と、前記第2加速度信号を二度積分した第3信号と、を算出し、
前記第1信号、前記第2信号及び前記第3信号に基づいて、前記生体の足が離地したタイミングと推定される離地タイミングを前記動作タイミングとして取得することを特徴とする測定方法。
A measuring method in a measuring device,
The measuring device is mounted on a living body, detects an acceleration corresponding to the movement of the living body, and outputs an acceleration signal corresponding to the acceleration; and an acceleration sensor when the living body is moving in a traveling direction. A control unit that acquires at least one operation timing related to the movement of the legs of the living body based on a signal ,
Based on the acceleration signal, the acceleration signal in the traveling direction is calculated as a first acceleration signal, and the acceleration signal in the vertical direction is calculated as a second acceleration signal,
A first signal obtained by smoothing the first acceleration signal using a moving average of the first acceleration signal in a first period, and a moving average of the first acceleration signal in a second period shorter than the first period are used. Calculating a second signal obtained by smoothing the first acceleration signal and a third signal obtained by integrating the second acceleration signal twice,
A measurement method, characterized in that, based on the first signal, the second signal, and the third signal, a takeoff timing estimated to be a timing at which the foot of the living body has taken off is acquired as the operation timing.
測定装置における測定プログラムであって、
前記測定装置は、生体に装着され、前記生体の動きに対応する加速度を検出して前記加速度に対応する加速度信号を出力する加速度センサと、前記生体が進行方向に移動しているときに前記加速度信号に基づいて、前記生体の足の動きに関する少なくとも一つの動作タイミングを取得する制御部と、を有し、
前記加速度信号に基づいて、前記進行方向の加速度信号を第1加速度信号として算出させるとともに、上下方向の加速度信号を第2加速度信号として算出させ、
前記第1加速度信号の第1期間の移動平均を用いて前記第1加速度信号を平滑化した第1信号と、前記第1加速度信号の前記第1期間より短い第2期間の移動平均を用いて前記第1加速度信号を平滑化した第2信号と、前記第2加速度信号を二度積分した第3信号と、を算出させ、
前記第1信号、前記第2信号及び前記第3信号に基づいて、前記生体の足が離地したタイミングと推定される離地タイミングを前記動作タイミングとして取得させることを特徴とする測定プログラム。
A measuring program in a measuring device,
The measuring device is mounted on a living body, detects an acceleration corresponding to the movement of the living body, and outputs an acceleration signal corresponding to the acceleration; and an acceleration sensor when the living body is moving in a traveling direction. A control unit that acquires at least one operation timing related to the movement of the legs of the living body based on a signal ,
Based on the acceleration signal, the acceleration signal in the traveling direction is calculated as a first acceleration signal, and the acceleration signal in the vertical direction is calculated as a second acceleration signal,
A first signal obtained by smoothing the first acceleration signal using a moving average of the first acceleration signal in a first period, and a moving average of the first acceleration signal in a second period shorter than the first period are used. Calculating a second signal obtained by smoothing the first acceleration signal and a third signal obtained by integrating the second acceleration signal twice,
A measurement program, characterized in that, based on the first signal, the second signal, and the third signal, a take-off timing estimated to be a timing at which the foot of the living body is taken off is acquired as the operation timing.
JP2016059469A 2016-03-24 2016-03-24 Measuring device, measuring method, and measuring program Active JP6690343B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016059469A JP6690343B2 (en) 2016-03-24 2016-03-24 Measuring device, measuring method, and measuring program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016059469A JP6690343B2 (en) 2016-03-24 2016-03-24 Measuring device, measuring method, and measuring program

Publications (2)

Publication Number Publication Date
JP2017169837A JP2017169837A (en) 2017-09-28
JP6690343B2 true JP6690343B2 (en) 2020-04-28

Family

ID=59970692

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016059469A Active JP6690343B2 (en) 2016-03-24 2016-03-24 Measuring device, measuring method, and measuring program

Country Status (1)

Country Link
JP (1) JP6690343B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2019216799A (en) * 2018-06-15 2019-12-26 カシオ計算機株式会社 Running analyzer, running analyzing method, and running analyzing program
JP7222385B2 (en) * 2020-09-28 2023-02-15 カシオ計算機株式会社 Measuring device, measuring method and program

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6018705A (en) * 1997-10-02 2000-01-25 Personal Electronic Devices, Inc. Measuring foot contact time and foot loft time of a person in locomotion
JP5017062B2 (en) * 2007-11-01 2012-09-05 シチズン・システムズ株式会社 Body motion detection device
JP5733503B2 (en) * 2011-02-28 2015-06-10 国立大学法人広島大学 Measuring device, measuring method, and measuring program
US20130085677A1 (en) * 2011-09-30 2013-04-04 Apple Inc. Techniques for improved pedometer readings
JP6390303B2 (en) * 2014-09-22 2018-09-19 カシオ計算機株式会社 Measuring apparatus, measuring method and measuring program

Also Published As

Publication number Publication date
JP2017169837A (en) 2017-09-28

Similar Documents

Publication Publication Date Title
KR102131358B1 (en) User interface device and method of operation of user interface device
FI120133B (en) A wrist unit and a method for determining motion information
US10694980B2 (en) Exercise support device, exercise support method and exercise support program
US10504381B2 (en) On-running landing position evaluation method, on-running landing position evaluation apparatus, detection method, detection apparatus, running motion evaluation method, and running motion evaluation apparatus
KR101693629B1 (en) Method for position estimation of pedestrian walking on loco-motion interface and apparatus thereof
EP3289435B1 (en) User interface control using impact gestures
US20150374307A1 (en) Data analysis device, data analysis method, and data analysis program
US10048087B2 (en) Measurement apparatus, measurement method and non-transitory computer readable recording medium with program stored thereon
WO2019033587A1 (en) Smart watch-based gps drifting filtering method, and smart watch
KR20150138831A (en) Portable device for controlling external apparatus via gesture and operating method for same
JP6690343B2 (en) Measuring device, measuring method, and measuring program
JP6496996B2 (en) Exercise quantity calculation method, exercise quantity calculation device, and portable device
US9192843B2 (en) Analysis system, analysis apparatus, electronic device, analysis method, and program
Liu et al. A zero-velocity detection method with transformation on generalized likelihood ratio statistical curve
JP6390303B2 (en) Measuring apparatus, measuring method and measuring program
WO2014181602A1 (en) Walking posture meter and program
JP2015146058A (en) Information processing apparatus, information processing method, and information processing program
JP6233123B2 (en) Sensor device, sampling method and program
KR101713496B1 (en) System and method for zero-delay real time step detection utilizing an accelerometer sensor
JP7275807B2 (en) Motion analysis device, motion analysis method and program
JP5621899B2 (en) Electronics
JP5952907B2 (en) Walking style discrimination system, walking style discrimination device, and walking style discrimination method
JP7076675B2 (en) Running method discrimination device, running method discrimination method and running method discrimination program
JP6521477B2 (en) Data analysis apparatus and data analysis method, data analysis program
JP2019101191A (en) Travel data evaluation system, information processing system, and information processing server

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190214

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20191121

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20191217

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200213

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200310

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200323

R150 Certificate of patent or registration of utility model

Ref document number: 6690343

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150