JP2009268831A - 酸素摂取量測定方法およびこの方法を用いる装置 - Google Patents
酸素摂取量測定方法およびこの方法を用いる装置 Download PDFInfo
- Publication number
- JP2009268831A JP2009268831A JP2008123945A JP2008123945A JP2009268831A JP 2009268831 A JP2009268831 A JP 2009268831A JP 2008123945 A JP2008123945 A JP 2008123945A JP 2008123945 A JP2008123945 A JP 2008123945A JP 2009268831 A JP2009268831 A JP 2009268831A
- Authority
- JP
- Japan
- Prior art keywords
- oxygen intake
- test object
- axis direction
- dimensional
- acceleration sensor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
【解決手段】被験者運動中の酸素摂取量のダイナミクス(時々刻々と変化する代謝量)を推定する。そのために、三次元加速度信号の瞬時振幅(Me)とbreath-by-breathで計測した酸素摂取量(Vo2)との間の伝達関数を、FIR(Finite Impulse Response)モデルを含んだ形のARX(Auto-Regressive eXogenous)モデルで解析し、このインパルス応答と三次元加速度センサの検知結果から、被験者の酸素摂取量のダイナミクスを推定する。
【選択図】 図28
Description
ここでA(q)はV−o2の自己回帰係数行列、B(q)は出力に対する入力の時間遅れ係数行列、w(k)は白色雑音である。A(q)の次数がゼロであればFIR(Finite Impulse Response)モデルとなる。次数とそれぞれの係数は最小二乗法で求める。三次元加速度信号から酸素摂取量Vo2のダイナミクスを推定するまでのデータ処理の流れは図28のようになる。
・加速度計本体100の大きさは5.9×8.3cm2でポータブル;
・加速度センサは10mm立方体であり、±2Gまで計測可能;
・計測データはメモリカード(フラッシュメモリ)に書き込む;
・計測時、三次元加速度計本体100は被験者の外部に固定;
・三次元加速度センサ部分の身体への装着方法は、X軸を前後、Y軸を左右、Z軸を上下方向として胸部剣状突起部にマジックテープ(登録商標)を用いて貼り付け;
・計測後のデータはUSBインターフェイス経由でパソコンに取り込む。
・被験者はトレッドミル上で20分間歩行運動を行う(図2参照);
・広帯域の伝達特性を得るためトレッドミルの速度は15秒おきに2.5km/hと5.5km/hでランダムに変化させる(PRBS運動);
・歩行運動中の酸素摂取量を呼吸流速計と質量分析器を用いてbreath-by-breathで計測する;
・酸素摂取量の計測と同時に三次元加速度計から10Hzのサンプリング速度で加速度信号を取り込む。
・被験者はトレッドミル上で速度2.5km/hで3分間、5.5km/hで5分間の歩行運動
(STEP運動)を、PRBS運動に引き続いて行う(図3参照);
・STEP運動時もPRBS運動同様に三次元加速度と酸素摂取量を同時計測する。
<加速度信号と呼吸データの処理>
01)加速度信号
計測したX、Y、Z方向の加速度データの二乗平均を前掲の式(1)により求めて瞬時Activity(A(t))とした。
breath-by-breathの酸素摂取量をFIとし、酸素濃度と呼気酸素濃度をFEとし、吸気流速と呼気流速をそれぞれV・ I、V・ Eとして、式(4)に示すように呼吸酸素濃度と呼吸流速の積分によって、酸素摂取量V・o2(t)を求めた。
11)計測した加速度信号の各成分を二乗平均することで瞬時Activity(A(t))求め(式(1)参照)、A(t)の瞬時振幅(AH(t))を、A(t)とそのHilbert変換(式(2)(3)参照)によって求める。
<ARXモデル>
システム同定には線形回帰モデルであるARX(Auto-Regressive eXogeneous)モデルを用いた(図4参照)。ARXモデルを次の式(5)に示す。(式(5)は前述した式(a)に対応する内容である。)
酸素摂取量ダイナミクスの推定手順は以下のとおりである。
<酸素摂取量と加速度平均瞬時振幅の関係>
酸素摂取量と加速度の瞬時振幅間における線形性の有無を調べるため、予備実験を行った。トレッドミル上で歩行速度を2.5km/h、→4.0km/h、→5.5km/hと変えて歩行運動を行ったところ、加速度の平均瞬時振幅A− H(t)と定常状態の酸素摂取量V・o2は線形関係にあることがわかった(図5参照)。
31)加速度信号
PRBS運動(図2参照)およびSTEP運動(図3参照)で計測された加速度信号の例を図6および図7に示す。図6からは歩行速度が2.5km/hから5.5km/hへ上昇するたびに加速度の振幅が大きくなっているのがわかる。また、図7からは、2.5km/hで歩行開始後の180秒目以後、5.5km/hへの歩行速度の上昇によって振幅が大きくなっているのがわかる。
データ解析後の、PRBS運動時およびSTEP運動時における酸素摂取量V・o2(t)と加速度の平均瞬時振幅A− H(t)の例を図8および図9に示す。図8からは、加速度平均瞬時振幅A− H(t)がランダムに上下するのに伴い酸素摂取量V・o2(t)も上下しているのがわかる。図9からは、STEP運動を開始してから180秒以後、加速度平均瞬時振幅A− H(t)の増加に伴い酸素摂取量V・o2(t)も増加しており、STEP運動を開始してからおよそ300秒付近で定常状態となっていることがわかる。
図10は、被験者Aについて、ARXモデルでシミュレートしたPRBS運動時の酸素摂取量V・o2(t)を例示している(図26(a)も参照)。被験者Aの場合、式(5)のA(q)とB(q)の次数na、nbは、それぞれ“0”と“43”になった(図32参照)。図11は被験者Aのシミュレーションによって得られたインパルス応答の一例を示している(図26(b)も参照)。その他の各被験者のA(q)、B(q)の次数は図32に例示されている。いずれの被験者も次数naは小さく、FIRモデルに近い伝達特性であった(次数naが“0”の被験者A〜CについてはFIRモデルをそのまま適用できる)。
図12は、システム同定によって得られたインパルス応答を用いて推定した、STEP運動時の酸素摂取量ダイナミクスの実測値と推定値の一例を示している。酸素摂取量V・o2(t)の実測値は加速度平均瞬時振幅A− H(t)の増加に比べゆっくりと増加しているが、酸素摂取量V・o2(t)の推定値もA− H(t)と同様の形を示しており、良い一致が見られた。また、実測値と推定値との間の誤差を示す残差は全体的に小さな値を示し、特定の傾向を示すことなく推定されており、加速度信号の平均瞬時振幅情報(A− H(t))を用いた酸素摂取量(V・o2(t))のダイナミクスの推定は可能であることがわかった。
図31は、三次元加速度センサ100Hが自転車のペダルに装備される場合を説明する図である。walkingやrunningと異なり、自転車による運動量測定の場合では、センサを被験者の体に装着しても正確なデータ収集は難しい。しかし、自転車のペダル部分に三次元加速度センサ100Hおよびロード測定素子(ペダルを漕ぐときの力を測定する感圧素子)を組み込むことにより、ペダルを漕ぐ力(運動負荷の大きさ)をパラメータに含めた上で、酸素摂取量を運動量として計測することができる。
(1)<酸素摂取量測定方法>
この酸素摂取量測定方法では、検査対象(被験者)の三次元空間における動きに伴う三次元方向の加速度の瞬時ベクトルを、互いに直交するx軸方向とy軸方向とz軸方向(図25の(a)(b)(c))毎に検出し、
検出したx軸方向の加速度信号(式(1)のx2(t))とy軸方向の加速度信号(式(1)のy2(t))とz軸方向の加速度信号(式(1)のz2(t))を平均化(式(1)の√(x2+y2+z2))して平均振幅(式(1)のA(t))を求め(図25の(d)のsqrt(x2+y2+z2))、
求めた平均振幅(式(1)のA(t))から、ヒルベルト変換(式(3))(積分変換の1つ)を用いて前記平均振幅の瞬時値(式(2)のAH(t)ejφ(t) の振幅AH(t)を平均化したA− H(t))(式(a)のM−e)を求め(図26の(a)ではA− H(t))、
求めた平均振幅の瞬時値(図5のA− H(t) または図24のM−e)と前記検査対象の酸素摂取量(式(a)のV−o2)(図5のV・o2または図23のsteady state Vo2)との間の対応関係(図5におけるV・o2とA− Hの比例関係、または図23におけるsteady state Vo2とsteady stateM−eの比例関係)から、前記検査対象の酸素摂取量(V・o2)を求める。
この酸素摂取量測定装置は、検査対象(被験者)の三次元空間における動きに伴う三次元方向の加速度の瞬時ベクトルを、互いに直交するx軸方向とy軸方向とz軸方向毎に検出する三次元加速度センサ(図29の100)と、
検出したx軸方向の加速度信号とy軸方向の加速度信号とz軸方向の加速度信号を平均化して平均振幅(A(t)))を求める構成(図29の210の処理機能)と、
求めた平均振幅(A(t))から、積分変換(好ましい例としてはヒルベルト変換)を用いて前記平均振幅の瞬時値(A− H(t))を求める構成(図29の210の処理機能)と、
求めた平均振幅の瞬時値(図5のA− H(t) または図24のM−e)と前記検査対象の酸素摂取量(図5のV・o2または図23のsteady state Vo2)との間の対応関係(図5におけるV・o2とA− Hの比例関係、または図23におけるsteady state Vo2とsteady stateM−eの比例関係)から、前記検査対象の酸素摂取量(V・o2)を求める構成(図29の210の処理機能)を具備する。
この発明を実施することにより、walkingあるいはrunning中の酸素摂取量のダイナミクスを三次元加速度センサの検知結果から推定することができる。そのために、三次元加速度信号の瞬時振幅とbreath-by-breathで計測した酸素摂取量間の伝達関数をFIR(Finite Impulse Response)モデルを含んだ形のARX(Auto-Regressive eXogenous)モデルで解析し、そのインパルス応答を求める。このインパルス応答の係数をフィルタとして機器にインプリメントすることにより、正確なエネルギー代謝計測装置を作ることができる。この装置は携帯型とすることが容易であり、安価に量産できる。
Claims (11)
- 検査対象の三次元空間における動きに伴う三次元方向の加速度を、互いに直交するx軸方向とy軸方向とz軸方向毎に検出し、
検出したx軸方向の加速度信号とy軸方向の加速度信号とz軸方向の加速度信号から平均振幅を求め、
求めた平均振幅から、積分変換を用いて前記平均振幅の瞬時値を求め、
求めた平均振幅の瞬時値と前記検査対象の酸素摂取量との間の対応関係から、前記検査対象の酸素摂取量を求める酸素摂取量測定方法。 - 前記対応関係が前記検査対象毎に異なる場合において、前記前記平均振幅の瞬時値を入力とし、前記検査対象の運動時酸素摂取量を出力とする伝達関数のインパルス応答を前記検査対象毎に求め、求めたインパルス応答から前記対応関係を決める請求項1に記載の方法。
- 前記積分変換としてヒルベルト変換を採用し、前記平均振幅の瞬時値と前記検査対象の酸素摂取量との間の対応関係として線形比例関係を採用する請求項1または請求項2に記載の方法。
- 検査対象の三次元空間における動きに伴う三次元方向の加速度を互いに直交するx軸方向とy軸方向とz軸方向毎に検出する三次元加速度センサと、
検出したx軸方向の加速度信号とy軸方向の加速度信号とz軸方向の加速度信号から平均振幅を求める構成と、
求めた平均振幅から、積分変換を用いて前記平均振幅の瞬時値を求める構成と、
求めた平均振幅の瞬時値と前記検査対象の酸素摂取量との間の対応関係から、前記検査対象の酸素摂取量を求める構成と
を具備した酸素摂取量測定装置。 - 請求項4に記載の前記三次元加速度センサが組み込まれ請求項4の装置により前記前記検査対象の酸素摂取量を求めるように構成された歩行計。
- 請求項4に記載の前記三次元加速度センサが組み込まれ請求項4の装置により前記前記検査対象の酸素摂取量を求めるように構成された自転車ペダル。
- 請求項4に記載の前記三次元加速度センサが組み込まれ請求項4の装置により前記前記検査対象の酸素摂取量を求めるように構成された靴。
- 請求項4に記載の前記三次元加速度センサが組み込まれ請求項4の装置により前記前記検査対象の酸素摂取量を求めるように構成されたバンド装置。
- 請求項4に記載の前記三次元加速度センサが組み込まれ請求項4の装置により前記前記検査対象の酸素摂取量を求めるように構成され前記検査対象により振り回されるように構成されたリモコン。
- 請求項4に記載の前記三次元加速度センサが組み込まれ請求項4の装置により前記前記検査対象の酸素摂取量を求めるように構成され前記検査対象が握るように構成された棒状装置。
- 請求項4に記載の前記三次元加速度センサが組み込まれ請求項4の装置により前記前記検査対象の酸素摂取量を求めるように構成された携帯電話。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008123945A JP4719768B2 (ja) | 2008-05-09 | 2008-05-09 | 酸素摂取量測定方法およびこの方法を用いる装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008123945A JP4719768B2 (ja) | 2008-05-09 | 2008-05-09 | 酸素摂取量測定方法およびこの方法を用いる装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2009268831A true JP2009268831A (ja) | 2009-11-19 |
JP4719768B2 JP4719768B2 (ja) | 2011-07-06 |
Family
ID=41435886
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008123945A Active JP4719768B2 (ja) | 2008-05-09 | 2008-05-09 | 酸素摂取量測定方法およびこの方法を用いる装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4719768B2 (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102657517A (zh) * | 2012-03-31 | 2012-09-12 | 戴晓晖 | 一种测量人体运动消耗的能量的装置 |
JP2013244357A (ja) * | 2012-05-29 | 2013-12-09 | Toyota Tsusho Corp | 活動量測定端末および活動量測定表示システム |
JP2014104260A (ja) * | 2012-11-29 | 2014-06-09 | Toyota Tsusho Corp | 活動量測定端末を用いた健康管理システム |
JP2017042627A (ja) * | 2016-11-22 | 2017-03-02 | 豊田通商株式会社 | 活動量測定端末を用いた健康管理システム |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR102029576B1 (ko) * | 2018-03-19 | 2019-10-07 | 재단법인대구경북과학기술원 | 생체 정보 추정 슈즈 및 이를 포함하는 시스템 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006263449A (ja) * | 2005-02-28 | 2006-10-05 | Jukunen Taiiku Daigaku Research Center | 脚筋力算出装置、脚筋力算出方法及びプログラム |
-
2008
- 2008-05-09 JP JP2008123945A patent/JP4719768B2/ja active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006263449A (ja) * | 2005-02-28 | 2006-10-05 | Jukunen Taiiku Daigaku Research Center | 脚筋力算出装置、脚筋力算出方法及びプログラム |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102657517A (zh) * | 2012-03-31 | 2012-09-12 | 戴晓晖 | 一种测量人体运动消耗的能量的装置 |
JP2013244357A (ja) * | 2012-05-29 | 2013-12-09 | Toyota Tsusho Corp | 活動量測定端末および活動量測定表示システム |
JP2014104260A (ja) * | 2012-11-29 | 2014-06-09 | Toyota Tsusho Corp | 活動量測定端末を用いた健康管理システム |
JP2017042627A (ja) * | 2016-11-22 | 2017-03-02 | 豊田通商株式会社 | 活動量測定端末を用いた健康管理システム |
Also Published As
Publication number | Publication date |
---|---|
JP4719768B2 (ja) | 2011-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2017386412B2 (en) | Systems and methods for real-time data quantification, acquisition, analysis, and feedback | |
KR101860654B1 (ko) | 운동 움직임 속성으로부터의 페이스 및 에너지 지출의 계산 | |
US11679300B2 (en) | Systems and methods for real-time data quantification, acquisition, analysis, and feedback | |
JP6183906B2 (ja) | 歩容推定装置とそのプログラム、転倒危険度算出装置とそのプログラム | |
RU2435521C2 (ru) | Способ детектирования и компенсации для контроля места активности на теле | |
CN107205661B (zh) | 使用来自多个设备的数据进行能量消耗计算 | |
JP6127873B2 (ja) | 歩行特徴の解析方法 | |
JP4719768B2 (ja) | 酸素摂取量測定方法およびこの方法を用いる装置 | |
Silsupadol et al. | Smartphone-based assessment of gait during straight walking, turning, and walking speed modulation in laboratory and free-living environments | |
JP6127455B2 (ja) | 歩行年齢評価方法 | |
Twomey et al. | Comparison of accelerometer-based energy expenditure estimation algorithms | |
Aubol et al. | Tibial acceleration reliability and minimal detectable difference during overground and treadmill running | |
WO2016194907A1 (ja) | 指標導出装置、ウェアラブル機器及び携帯機器 | |
Duhaylungsod et al. | Detection of gait abnormality through leg symmetry and temporal parameters | |
JP5888582B2 (ja) | エネルギー消費量呈示装置およびエネルギー消費量呈示装置の作動方法 | |
JP2020192307A (ja) | 下肢筋力評価方法、下肢筋力評価プログラム、下肢筋力評価装置及び下肢筋力評価システム | |
JP2009285269A (ja) | 人体構造の異常状態の物理的な観測と解析方法及び、同方法を用いた測定装置 | |
KR101171593B1 (ko) | 3차원 가속도 센서를 이용한 실시간 칼로리 계산방법 | |
CN108158587B (zh) | 室内人体运动量的测量方法和测量装置 | |
JP6552875B2 (ja) | 移動運動解析装置、方法及びプログラム | |
JP6303746B2 (ja) | 歩行目標の設定方法 | |
TWI542321B (zh) | Measurement of Energy Metabolism of Fat and Sugar Based on Heart Rate , Wearable electronic devices, and computer programs | |
Andreoni et al. | Smart garments for performance and training assessment in sport | |
TWI796035B (zh) | 生物力學評量系統,其生物力學感測裝置及生物力學評量平台 | |
Callaway | Quantification of performance analysis factors in front crawl swimming using micro electronics: a data rich system for swimming. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20101125 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20101207 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110207 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20110308 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110404 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140408 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |