(57)【特許請求の範囲】
【請求項１】
身体の各部位をそれぞれマッサージする複数のエアバッグと、前記身体の各部位がマッサージされる毎に心拍の脈波波形を検出する心拍検出手段とを備え、前記エアバッグの動作条件である施療要素に基づいてそのエアバッグを動作させるマッサージ装置であって、
前記心拍検出手段によって検出された脈波波形に基づいて前記身体の各部位がマッサージされる毎に心拍周期を算出する心拍周期検出手段と、
前記心拍周期検出手段により算出された心拍周期に基づいて、単位時間内の心拍周期の平均である平均心拍周期と、単位時間内の心拍周期の変化の割合である心拍変化率と、単位時間内の心拍周期の長さの変化の大きさである心拍変動度とをそれぞれ求める心拍特徴検出手段と、
前記平均心拍周期と前記心拍変化率と前記心拍変動度とから前記身体の各部位毎における快適度を算出する施療要素変更手段と、
該施療要素変更手段によって算出された快適度に応じて施療要素を変更する動作制御部と、を備え、
その各部位毎に変更された施療要素に基づいて前記各エアバッグを動作させることを特徴とするマッサージ装置。
【請求項２】
前記心拍周期検出手段は、前記脈波波形毎の振幅の一定割合を閾値とし、この閾値と前記脈波波形とのクロスポイントを心拍周期の基準として心拍周期を算出することを特徴とする請求項１に記載のマッサージ装置。
【請求項３】
前記施療要素の変更は、前記各部位毎に行うマッサージについて施療頻度の増減と施療時間の増減と施療の強さの増減のうち少なくともいずれか１つを行うこと又は施療を実施から非実施に変えることを特徴とする請求項１に記載のマッサージ装置。
【請求項４】
前記各エアバッグを予め設定した順序で動作させてマッサージ部位を順じ変えてマッサージを行っていく自動マッサージモードにおいて、
前記心拍特徴検出手段は、その最初のサイクルで前記身体の各部位毎に前記平均心拍周期と前記心拍変化率と前記心拍変動度とを検出し、
前記施療要素変更手段は、その各部位毎の前記平均心拍周期と前記心拍変化率と前記心拍変動度とに基づいて快適度を算出していくことを特徴とする請求項１に記載のマッサージ装置。
【請求項５】
前記平均心拍周期と前記心拍変化率と前記心拍変動度とに基づいて各部位毎に求められた快適度を表示することを特徴とする請求項１に記載のマッサージ装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この発明は、心地よい部位を自動的に探してその部位をマッサージするマッサージ装置に関する。
【０００２】
【従来の技術】
従来から、特開２００１−１９０６１６号公報や特開平９−１９２１８６号公報に示すマッサージ機が知られている。
【０００３】
前者のマッサージ機は、リモコン装置に設けたセンサ部から皮膚温、心拍周期、血圧などの生体情報を検出し、この検出した生態情報に基づいて使用者の体調に合わせたマッサージを行うようにしたものである。
【０００４】
また、後者のマッサージ機は、心拍数,呼吸数や脳波等の生体情報を検出する検出装置と、動作速度が可変となっているマッサージ機構と、検出装置の検出出力に応じてマッサージ機構の動作速度を変更する制御部とからなり、生体情報を参照しつつマッサージ機構の動作速度を変更していくものである。
【０００５】
【発明が解決しようとする課題】
しかしながら、このような従来のマッサージ機では、人体の各部位毎に心地よさを検出してマッサージを行う構成になっていない。このため、心地よい部位だけをマッサージする場合、使用者は操作ボタンを操作して心地よい部位をいちいち選択していかなければならないという問題があった。
【０００６】
この発明の目的は、自動的に心地よい部位を探してその部位をマッサージすることのできるマッサージ装置を提供することにある。
【０００７】
【課題を解決するための手段】
この発明は、上記目的を達成するため、身体の各部位をそれぞれマッサージする複数のエアバッグと、前記身体の各部位がマッサージされる毎に心拍の脈波波形を検出する心拍検出手段とを備え、前記エアバッグの動作条件である施療要素に基づいてそのエアバッグを動作させるマッサージ装置であって、
前記心拍検出手段によって検出された脈波波形に基づいて前記身体の各部位がマッサージされる毎に心拍周期を算出する心拍周期検出手段と、
前記心拍周期検出手段により算出された心拍周期に基づいて、単位時間内の心拍周期の平均である平均心拍周期と、単位時間内の心拍周期の変化の割合である心拍変化率と、単位時間内の心拍周期の長さの変化の大きさである心拍変動度とをそれぞれ求める心拍特徴検出手段と、
前記平均心拍周期と前記心拍変化率と前記心拍変動度とから前記身体の各部位毎における快適度を算出する施療要素変更手段と、
該施療要素変更手段によって算出された快適度に応じて施療要素を変更する動作制御部と、を備え、
その各部位毎に変更された施療要素に基づいて前記各エアバッグを動作させることを特徴とする。
【０００８】
【実施の形態】
以下、この発明に係るマッサージ装置の実施の形態を図面に基づいて説明する。
【０００９】
図１に示すマッサージ装置１０は、座部１１と、この座部１１の後ろに設けた背もたれ部１２と、座部１１の前に設けたオットマン１３と、座部１１の両側に設けたアームレスト部１４,１５と、座部１１の下側に設けられた給排気装置１００とを有している。背もたれ部１２はリクライニングが可能となっており、モータＭ１（図３参照）を備えた図示しない電動機構により背もたれ部１２が傾動するようになっている。
【００１０】
背もたれ部１２の前面には、首肩用のエアバック（施療子）２０と、背筋用のエアバック（施療子）２１と、背中用の一対のエアバック（施療子）２２,２３と、腰用のエアバック（施療子）２４とが設けられている。座部１１の上面には、臀部用のエアバック（施療子）２５と、もも用のエアバック（施療子）２６とが設けられている。また、オットマン１３には左右に一対のエアバック（施療子）２７,２８がそれぞれ設けられている。首肩用のエアバック２０は、図示しない電動機構により上下にスライド移動できるようになっており、この電動機構はモータＭ２（図３参照）を備えている。
【００１１】
各エアバック２０〜２８は、座部１１,背もたれ部１２,オットマン１３に取り付けられたカバーＫによって覆われている。
【００１２】
アームレスト部１４の上面の前側には表示部３０が設けられており、この表示部３０には、図２に示すように、オン・オフスイッチ３１と、全身コースと上半身コースと下半身コース等を選択する選択スイッチ３２と、エアバック２０を上にスライド移動させる上スイッチ３３Ｓ１と、エアバック２０を下にスライド移動させる下スイッチ３３Ｓ２と、背もたれ部１２を倒していく倒スイッチ３４Ｓ１と、背もたれ部１２を起こしていく起スイッチ３４Ｓ２と、各エアバック２０〜２８の膨・縮の速さを調整するパルススイッチ３５と、マッサージの強さを調整する強さスイッチ３６と、快適お探しモードと自動コースモードなどを選択する選択スイッチ３７と、納得スイッチ３８と、足同時スイッチ３９等とが設けられている。各スイッチ３１〜３９はタッチパネルで構成され、そのタッチパネルに触れることにより各スイッチ３１〜３９が操作されるようになっている。
【００１３】
また、表示部３０には、首／肩,背中,背筋,腰,しり,もも,脚の心地よさを示す４つのランプ画像Ｌ１ａ〜Ｌ１ｄ,…Ｌ７ａ〜Ｌ７ｄがそれぞれ表示され、その心地よさに応じてランプ画像Ｌ１ａ〜Ｌ１ｄ,…Ｌ７ａ〜Ｌ７ｄが点灯表示するようになっている。
【００１４】
また、表示部３０には、その心地よさを示す棒グラフ画像４０が表示される他に人体を示す人体画像４１と、各エアバック２０〜２８に対応したエアバック画像Ｇ２０〜Ｇ２８とが表示され、各エアバック画像Ｇ２０〜Ｇ２８にはタッチパネルが設けられており、これらタッチパネルに触れることにより各エアバック画像Ｇ２０〜Ｇ２８に対応した各エアバック２０〜２８が選択され、その選択されたエアバックだけを膨縮させることができるようになっている。
【００１５】
給排気装置１００は、図３に示すように、ポンプ１０１と複数の電磁弁１０２〜１０９と制御装置１１０とを有している。
【００１６】
電磁弁１０２は、開成されるとポンプ１０１とエアーバック２０とを連通してポンプ１０１からエアーバック２０にエアーを給気するようになっている。電磁弁１０２が閉成されるとエアーバック２０が大気に開放されてエアーバック２０のエアは排気されていくようになっている。
【００１７】
同様に、各電磁弁１０３〜１０９は、開成されるとポンプ１０１とエアーバック２１〜２８とを連通してポンプ１０１からエアーバック２１〜２８にエアーを給気するようになっており、電磁弁１０３〜１０９が閉成されるとエアーバック２１〜２８が大気に開放されてエアーバック２１〜２８のエアは排気されていくようになっている。
【００１８】
制御装置１１０は、ＣＰＵ等を有しており、表示部３０のスイッチ操作と心拍検出部２００の検出信号に基づいてポンプ１０１,電磁弁１０２〜１０９,モータＭ１,Ｍ２等を制御する。
【００１９】
また、制御装置１１０は、図４に示すように演算部１２０と動作制御部１３０とを有している。演算部１２０は、心拍検出部２００の検出信号から心拍周期を算出する心拍周期算出手段（心拍周期検出手段）１２１と、この心拍周期算出手段１２１が算出した心拍周期から平均心拍周期を算出する平均心拍周期算出手段（心拍特徴検出手段）１２２と、その心拍周期から心拍変化率を算出する心拍変化率算出手段（心拍特徴検出手段）１２３と、その心拍周期から心拍変動度を算出する心拍変動度算出手段（心拍特徴検出手段）１２４と、平均心拍周期算出手段１２２が算出した平均心拍周期と心拍変化率算出手段１２３が算出した心拍変化率と心拍変動度算出手段１２４が算出した心拍変動度とのうちの少なくとも１つから快適度を算出する快適度算出手段（施療要素変更手段：施療要素算出手段）１２５と、快適度算出手段１２５が算出する快適度とこの快適度に対応した各部位とを記憶するメモリ（記憶手段）１２７とを備えている。１２６は心拍周期を算出するための閾値を設定する閾値設定手段である。この閾値設定手段１２６は必ずしも必要とするものではない。
【００２０】
動作制御部１３０は、メモリ１２７に記憶された各部位の快適度に応じて、各エアバック２０〜２８の膨縮回数（動作条件：施療要素）、膨縮時間（動作条件：施療要素）、膨張度によるマッサージの強さ（動作条件：施療要素、）などを求めてポンプ１０１および電磁弁１０２〜１０９を制御していく。また、動作制御部１３０は、表示部３０のスイッチ操作に基づいてモータＭ１,Ｍ２を制御したり、表示部３０の画像表示を制御するようになっている。
【００２１】
心拍検出部２００は、耳朶に装着して耳朶の血流に応じた血流信号を出力する脈波センサ（心拍センサ：心拍検出手段）２０１と、血流信号を増幅する増幅回路２０２と、増幅した血流信号をＡ／Ｄ変換するＡ／Ｄ変換回路２０３とを備えている。
【００２２】
脈波センサ２０１は、耳朶に向けて光りを発光する発光素子と、その耳朶を透過する光を受光する受光センサ等とから構成されている。ここでは、脈波センサ２０１は耳朶の血流を検出するが、指などの血流を検出するものでもよい。また、受光素子から十分な大きさの血流信号が得られる場合には増幅回路２０２は不要である。また、Ａ／Ｄ変換回路２０３も制御装置１１０の演算部１２０にＡ／Ｄ変換回路があれば不要となる。
【００２３】
そして、心拍検出部２００と演算部１２０とでリラックス状態推定装置が構成されている。すなわち、脈波センサ２０１と心拍周期算出手段１２１と平均心拍周期算出手段１２２と心拍変化率算出手段１２３と心拍変動度算出手段１２４と快適度算出段１２５とでリラックス状態推定装置が構成されている。このリラックス状態推定装置は、平均心拍周期算出手段１２２と心拍変化率算出手段１２３と心拍変動度算出手段１２４のうち少なくともいずれか１つの算出手段を備えていればよい。
【００２４】
ここで、脈波波形の一例を図５に示す。図５において、横軸が時間、縦軸が  透過光量（すなわち血流量）を示している。血流量は心臓の鼓動で変化するものであるから、脈波波形の１つの山が心臓一拍分とみなすことができる。
【００２５】
従って、例えば脈波波形のピークを検知し、次のピークまでの時間を計測すれば心拍の周期が求められることになる。尚、図５ではピーク検出を脈波波形の  山のピークで行っているが、谷（バレー）を検知しても同様に心拍周期を求める事ができる。また、ピーク検知を行う以外にもゼロクロス法等で心拍周期を求めてもよい。これら心拍周期は心拍周期算出手段１２１が求めていく。
【００２６】
このようにして、得られた瞬時心拍周期から、評価対象単位時間内の心拍周期を１つのグループとみなし、そのグループ内での変化の仕方を捉える事で、マッサージ装置の使用者の心地よさの程度、すなわちリラックス状態（マッサージ効果）を数値として推定することができる。
【００２７】
ここで  評価対象単位時間とは、同じマッサージを続けている時間であり、発明者らの実験によれば５秒から１０秒程度がもっともマッサージ効果の推定が  しやすいことが判った。これは、これよりも短い時間では情報量（心拍数）が  少なく、正確に推定できず、逆に長い時間では使用者のマッサージに対する感覚が変わるため、その感覚によって心拍周期が変化したのかマッサージ効果によって心拍周期が変わったのかが判断できないからである。
【００２８】
次に、このマッサージ効果の数値化について説明する。
【００２９】
マッサージ効果の推定は、平均心拍周期、心拍変化率、心拍変動度の３つの  値のうち少なくとも１つを算出する事で実現できる。以下、その具体的計算方法を示す。
【００３０】
まず、対象単位時間内（対象刺激単位内）の心拍周期の平均値である平均心拍周期は下記の（１）式で求めることができる。
【００３１】
平均心拍周期＝（Σｙi）／ｎ        …（１）
ここで、（１）式のｙiは対象単位時間内のi番目の心拍周期で、ｎは対象時間内のデータ数（心拍周期の数）である。
【００３２】
演算部１２０の平均心拍周期算出手段１２２は、心拍周期算出手段１２１が算出した心拍周期ｙiに基づいて（１）式から平均心拍周期を算出していく。
【００３３】
この平均心拍周期の値が大きいほどリラックスの度合（心地よさの度合）が大きいと判断する。
【００３４】
心拍変化率は、評価対象単位時間内で心拍周期が長くなってゆくのか、逆に短くなっていくのかを表すものであり、ある評価対象単位時間（評価対象刺激単位）内で横軸に時間を縦軸に心拍周期をとった図６のグラフＦ１を直線近似した場合のその直線Ｆ２の傾きで定義されるものである。
【００３５】
この傾き（心拍変化率）は、（２）式によって求められる。
【００３６】
心拍変化率＝｛ｎΣｘiｙi−(Σｘi)(Σｙi)｝／｛ｎΣｘi2−(Σｘi)2｝…(２)
ここで、ｎはデータ数、ｙiはi番目の心拍周期、ｘiはi番目までの経過時間、すなわちこの経過時間はｙ１からｙiまでの心拍周期の累積時間を示す。
【００３７】
この心拍変化率が大きいほどリラックスしている度合いが大きいと判断する。
【００３８】
この心拍変化率は、演算部１２０の心拍変化率算出手段１２３が心拍周期算出手段１２１の心拍周期ｙiに基づいて（２）式から算出するものである。
【００３９】
次に、心拍変動度の算出方法について説明する。
【００４０】
図７はあるマッサージ刺激を１分間継続したときの心拍周期を、横軸：時間、縦軸：心拍周期としてプロットしたグラフを示すものであり、＋印は心拍周期そのものを（生データ）示す。●印は５つの生データの平均値を算出したものであり、この５つの生データの取り方は、１番目から５番目までの平均と、２番目から６番目までの平均と、…平均とを取っていくものであり、５つの生データを１つづつずらしながら順番に取って、その各々の平均を求めた値（第２移動平均値）である。
【００４１】
また、図７において、○印は９つの生データの平均値を算出したものであり、９つの生データを１つづつずらしながら順番に取って、その各々の平均を求めた値（第１移動平均値）である。
【００４２】
この＋印の生データは、マッサージ刺激終了後にマッサージ使用者に対するインタビューを行い、気持ちよかったという回答が有ったもののデータを示したものであり、気持ちよかった場合には、図７に示すように心拍周期が周期的に変動することが発明者らの実験で明らかになった。
【００４３】
これは、心拍周期の移動平均値、心拍周期の心拍変化率や周期的変動度の大きさからリラックス状態の度合いを検出することができることを意味するものである。
【００４４】
そこで、この周期的変動の大きさを心拍変動度として算出し、この心拍変動度からマッサージ効果、すなわちリラックス状態を推定すればより的確なリラックス状態の推定が可能になる。
【００４５】
ところで、移動平均は、ローパスフィルターと同等の効果があり、平均をとる数（ｎ）を増やすと、波形がだんだん滑らかになる。そこで、ｎ＝９の場合とｎ＝５の場合とをそれぞれ算出し、それぞれ全体のばらつき、すなわち標準偏差を算出し、その差を算出すれば周期的な変動のある場合のみ値が大きくなり、周期変動の大きさを数値化する事が可能である。
【００４６】
つまり、周期変動の大きさ（心拍変動度）は次式で求める事が可能である。
【００４７】
0004017915.tif 000002
ただし、ｎ１＞ｎ２
尚、周期的変動の大きさを求めるには、フーリエ変換等の周波数解析手法を使うことも可能であるが、上の方法の様に単位データ数の異なる２種類の移動平均値を算出して、それぞれのデータ群の標準偏差を算出する方法で心拍変動度の計算を行うことで、フーリエ変換などの複雑な計算を行うことなく、周期的変動の大きさを求めることが出来る。この心拍変動度が大きいほどリラックスしている度合いが大きいと判断する。
【００４８】
図７から心拍周期の周期的変動の周波数は０.１Ｈｚ程度（周期が約１０秒）であることが分かり、評価対象単位時間（同じマッサージを行う時間）が短い場合、例えば１０秒程度では、山１つが収まるだけである。そこで、心拍変動を数値化する方法として次のような方法も適用する事が可能である。
【００４９】
図８は、ある評価対象単位時間内で横軸に時間を縦軸に心拍周期をとったグラフを示した例である。
【００５０】
マッサージの効果に応じて、周期的な変動を起こす例は、短い時間（図８では１０秒）のデータを見ると、いったん心拍周期が短くなった後に長くなる。すなわち、谷型のグラフとなる（山型の場合もある）。そこで、この谷型の分布を放物線近似して、この近似放物線の２乗係数をマッサージ効果（リラックス状態）の推定に用いる事が可能である。
【００５１】
この様に、評価対象単位時間内の心拍周期の時系列データを放物線近似した場合の２次係数として心拍変動度を求める事で、短い時間で的確にマッサージ効果、すなわちリラックス状態を推定することができる。
【００５２】
近似放物線の２乗係数は次の（４）式で求められる。
【００５３】
0004017915.tif 000003
ただし、ｎはデータの個数
ｙiはi番目の心拍周期
ｚiは  ｘi−ｘ
ｘiは  i番目までの累積時間
ｘは  平均心拍周期
心拍変動度は、演算部１２０の心拍変動度算出手段１２４が心拍周期算出手段１２１の心拍周期ｙiに基づいて（３）式あるいは（４）式から算出するものである。
【００５４】
図１２は実験結果の一例を示したものである。この実験はエアーバッグ式マッサージ装置で、もも、脚、首の三箇所を各一分間マッサージした際の心拍周期を測定したものであり、３つの刺激の前後（各一分間）に比較対象用に安静（マッサージなし）の心拍も測定している。
【００５５】
横軸が時間、縦軸が心拍周期を示している。実験において被験者（マッサージを受ける人）はグリップスイッチを保持しており、自身が気持ちよいと感じている場合にボタンを押すように指示してあり、実験システムはボタンを押している間だけプロットを黒丸から白丸に変えるように作成されている。つまり図１２の白丸プロット部分は被験者自身が気持ちよいと感じている時を指している。
【００５６】
図１２をみると脚のマッサージ中（ほぼ全域で白丸プロット）に心拍周期が１０秒前後の周期で変動を繰り返している事がわかる。また、同一の領域で緩やかな右上がりをしている事がわかる。さらに首、もも、のマッサージ中も含めて、白丸プロット部分は全体的に黒丸プロットに比べて、心拍周期が長い傾向が読み取れる。これは、心拍周期の長さからリラックス状態の度合いを検出することができることを示すものである。
【００５７】
発明者は、同様な実験を繰り返すことで、心拍平均周期、心拍変化率、心拍変動度が刺激（マッサージ）を受ける人のリラックス状態（心地よいという感覚）と相関があることを見い出した。すなわち、心拍平均周期、心拍変化率、心拍変動度が大きいほどリラックス状態の度合いが大きいことが分かった。
【００５８】
また、図１２の例では同一の人物で心拍平均周期と心拍変化率、心拍変動度の３つに変化が現れたが、人によってはこのうち１つあるいは２つだけがリラックスした場合に変化する事もあった。そこで、どのようなタイプの人が使っても同じようにリラックス状態（リラックス度）の推定が行えるように下記の（５）式を定義した。
【００５９】
この（５）式により、リラックスした場合に例えば心拍周期が変化する（遅くなる）だけのタイプの人でも、他の刺激部分との差をとることによって当該刺激のリラックス状態が相対的に定量化できるようになった。尚、（５）式の重み定数Ｃ１、Ｃ２、Ｃ３は多数の実験を繰り返し、その結果を用い快適度が各被験者の印象に最もよく合うように最適化した値を用いた。また、年代や性別毎にカテゴライズしたデータを用いて各カテゴリー毎に最適化を行う事も可能である。この場合は、刺激を受ける人がどの分類に属するかの情報を与える事でより厳密にリラックス状態の推定を行う事が可能である。
【００６０】
リラックス状態（心地よさ）＝Ｃ１×平均心拍周期＋Ｃ２×心拍変化率＋Ｃ３×（心拍変動度）-2…（５）
リラックス状態である快適度は、演算部１２０の快適度算出手段１２５によって（５）式に基づいて算出する。
【００６１】
この様に算出したリラックス状態（マッサージ効果の推定結果：心地よさ）をもとに、効果の高かったところを優先的にマッサージする事で使用者が満足できるマッサージを実現する事が可能となる。
【００６２】
ところで、心拍周期を検出する際に、ピーク検出の精度不足やノイズ等の影響で正しく検出できない場合がある。また、図９の例に見られるように心地よさを感じていても、心拍周期は一般にばらつきが大きく、直線近似や放物線近似がうまくできない場合がある（図９の●印は心拍周期の生データであり、○は５つの生データの平均値を算出した移動平均値である）。
【００６３】
そこで、生データから移動平均を算出し、そのデータを用いてその後の解析を行うと、ノイズ等の影響が小さくなり、より的確にマッサージ効果の推定が可能になる。
【００６４】
移動平均を算出する単位データ数があまり大きくなると、周期的な変動が埋もれてしまう。また、逆に小さすぎるとノイズに対するフィルター効果が得られない。発明者らの実験結果では、単位データ数は５前後がもっとも適していた（図９の○印）。ただし、心拍周期は脈波波形の注目点（先に述べた例では山のピーク）を順次カウントして求めているため、仮に注目点検知に誤差があって、心拍周期が現実より短く検出されると、次の心拍はその分長く検出されることになる。したがって、注目点のずれが大きい場合にはそれが相殺されるように、単位データ数は偶数の値をとることが望ましい。それ以外の場合は、所定位置の心拍周期の前後同数が取れるように単位データを奇数とすればよい。
【００６５】
先に述べた例では、脈波波形から心拍周期を求める方法として、脈波の山もしくは谷のピークを検出し、そのピークと次のピークまでの時間を測定する方法について説明したが、それ以外に２つの別な方法でも可能である。以下、その方法について説明する。
【００６６】
図１０に示すように、脈波波形上に一定の閾値を設定し、その閾値と脈波波形の例えば立ち上がり波形がクロスする点を検知し、クロス点から次のクロス点までを心拍周期とする事ができる。
【００６７】
一般に、山または谷のピーク付近はなだらかな変化をしているため、時間軸方向でどのポイントが真のピークかが判定しにくい場合がある。これに比べて急峻な立ち上がり部分でのクロスポイントは正確に時間を定める事が可能である。
【００６８】
図１１は、もう一つの方法を示すものであり、一般に、脈波波形は時間軸方向のみでなく、振幅方向にも変動する。そこで、谷側のピークと山側のピークから当該波形の一波分の振幅を求め、この振幅に対し一定の割合（図では５０％）を閾値とし、それと脈波波形のクロスするポイントを求める事でも心拍周期を定義することができる。
【００６９】
この場合、一定の閾値を用いる場合に比べて、振幅方向の変動の影響が小さくなるというメリットを有している。
【００７０】
ここで、脈波波形の振幅方向の変動は、心臓からの血流量の変化以外に、血管の緊張度合いなどの影響を受けているものと考えられる。この血管の緊張の度合いは自律神経系の活動状態によっても変化するものであり、例えばマッサージ効果によって使用者がリラックスできたとすると、脈波の振幅に変化が現れる事も起こりえる。
【００７１】
また、この振幅方向の変動はマッサージ椅子の使用者の年齢などによっても  変わると考えられる。従って、これらの心拍周期の検出方法のうち、その使用者に適したものを用いればよい。
【００７２】
発明者らの実験によれば、脈波振幅の変動に応じて閾値レベルを変化させる変動型閾値設定方法を用いてマッサージ効果を推定した結果がマッサージ椅子の使用者自身の感覚に合うと答える人が多かった。
【００７３】
先に述べたように、マッサージ装置にはもみ玉を機械的に動作させてツボを刺激するメカ式のものとエアーバッグを収縮させて刺激を得るエアーバッグ式がある。また、肩や腰などの上半身はメカ式で、ふくらはぎなどの下半身はエアーバッグで挟み込むように刺激するハイブリット型のものもある。
【００７４】
ところで、図９の測定例はエアーバッグ式のマッサージ装置を用いて心拍周期を測定したものであるが、一旦短くなってその後長くなる谷型の変動となっている１０秒から２０秒までの期間では、マッサージ後の使用者に対するインタビューで、使用者が気持ちよかったと答えた部分である。この一旦短くなって長くなるタイミングがエアーバッグの収縮のタイミングとほぼ一致することが実験的に確認されている。
【００７５】
すなわち、心拍周期が一旦短くなる区間ではエアーバッグは膨張し、心拍周期が増加している区間ではエアーバッグが収縮している。これは、エアーバッグの刺激により末梢血管が収縮した事の影響と考えられるが、同じ刺激でも気持ちよいと答えなかった場合には谷型にならない傾向が強いことも分かった。すなわち、実際に使用者が気持ちよいと感じた場合にのみエアーバッグの収縮と同期した
変動が見られる傾向が強い。
【００７６】
一方、メカ式のマッサージ機の場合には、一般にもみ玉が数ヘルツから十数ヘルツで振動する刺激であるから、このような谷型の変動は見られなかった。
【００７７】
従って、本マッサージ効果の推定方法は、心拍変動度による推定が有効になるという意味で、エアーバッグ式のマッサージ機において特に有効である。いうまでも無く、ハイブリットタイプのエアーバッグの収縮による下半身マッサージにも同様に効果が高い。
［動  作］
次に、上記のように構成されるマッサージ装置１０の動作を図１３に示すフロー図を参照しながら説明する。
【００７８】
先ず、電源プラグＰをコンセント(図示せず)に接続する。この接続により表示部３０に図２に示す画面が表示される。そして、使用者は例えば選択スイッチ３７を指で触れると、快適自動モードと自動コースとお好みコースとが順次選択されていく。この選択スイッチ３７により快適自動モードを選択し（ステップ１）、脈波センサ２０１を耳朶に装着し、オン・オフスイッチ３１をオンすると、ステップ２で脚用のエアバック２７,２８が膨縮されていく。ステップ３ではエアバック２６が膨縮され、ステップ４ではエアバック２５が膨縮されていく。同様にして、ステップ５ないしステップ９ではエアバック２４〜２０,２７,２８が順次膨縮されていく。
【００７９】
それぞれのエアバック２０〜２８の膨縮の回数は１回であり、膨縮の時間は５秒から１０秒の間である。
【００８０】
各エアバック２０〜２８の膨縮時の血流信号（心拍信号）を脈波センサ２０１が出力し、心拍周期算出手段１２１がその心拍信号から心拍周期を逐一算出していく。そして、平均心拍周期算出手段１２２がその逐一算出される心拍周期から平均心拍周期を算出し、心拍変化率算出手段１２３がその逐一算出される心拍周期から心拍変化率を算出し、心拍変動度算出手段１２４がその逐一算出される心拍周期から心拍変動度を算出する。
【００８１】
快適度算出手段１２５は、算出された平均心拍周期と心拍変化率と心拍変動度とから各部位であるもも,しり,腰,背筋,背中,首／肩,脚の快適度（リラックス状態）を算出し、この算出した快適度から快適度のランク（施療評価）を判定していく。すなわち、快適度算出手段１２５は各部位毎に快適度のランクを求めていく。この算出された各部位毎の快適度と判定したランクはメモリ１２７に記憶されていく。
【００８２】
そして、ステップ１０では、メモリ１２７に記憶された各部位毎の快適度のランク、すなわち施療評価に応じてランプ画像Ｌ１ａ〜Ｌ１ｄ,…Ｌ７ａ〜Ｌ７ｄがそれぞれ点灯表示される。また、その快適度に応じて棒グラフ画像４０が表示される。
【００８３】
使用者は、ランプ画像Ｌ１ａ〜Ｌ１ｄ,…Ｌ７ａ〜Ｌ７ｄの点灯表示を見て納得できれば納得スイッチ３８をタッチする。ステップ１１では納得スイッチ３８がタッチされたか否かが判断され、タッチされればイエスと判断されてステップ１２へ進む。ステップ１２では、判定結果、すなわちメモリ１２７に記憶された各部位毎の快適度のランクに応じて各エアバックが膨縮されていく。
【００８４】
各エアバック２０〜２８の膨縮回数は快適度のランクに応じて多くなる。例えば、ランプ画像Ｌ１ａ〜Ｌ７ａが点灯すれば膨縮回数は１回、ランプ画像Ｌ１ｂ〜Ｌ７ｂが点灯すれば膨縮回数は２回、ランプ画像Ｌ１ｃ〜Ｌ７ｃが点灯すれば膨縮回数は３回、ランプ画像Ｌ１ｄ〜Ｌ７ｄが点灯すれば膨縮回数は４回となる。
【００８５】
例えば、ランプ画像Ｌ１ａとランプ画像Ｌ２ａ,Ｌ２ｂとランプ画像Ｌ３ａ〜Ｌ３ｃとランプ画像Ｌ４ａ〜Ｌ４ｄとランプ画像Ｌ６ａとランプ画像Ｌ７ａとが点灯した場合、例えば最初にエアバック,２７,２８が１回膨縮され、この後エアバック２６が１回膨縮される。そして、エアバック２５は膨縮されず（非実施）、エアバック２４が４回膨縮され、この後エアバック２１が３回膨縮される。そして、エアバック２２,２３が２回膨縮され、エアバック２０が１回膨縮されていく。そして、これら動作が繰り返し行われていく。
【００８６】
このように、自動的に心地よい部位を探し、その探した各部位の快適度に応じて各エアバック２０〜２８の膨張回数が多く行われることにより、気持ちよいと感じる部分だけが重点的にマッサージされることになる。このため、マッサージ効果を上げることができ、使用者は気持ちよさを感じて十分にリラックスすることができることになる。また、心地よい部位を自動的に探していくので、使用者は操作ボタンを操作して心地よい部位をいちいち選択していく必要がない。
【００８７】
ステップ１３では、マッサージを行っている時間が予め設定されている時間に達したか否かが判断され、ノーであればステップ１２へ戻り、イエスであれば終了する。
【００８８】
ところで、ランプ画像Ｌ１ａ〜Ｌ１ｄ,…Ｌ７ａ〜Ｌ７ｄの点灯表示の数に納得の行かない場合には、背中や腰等の表示部（各部にタッチパネルが設けられている。）をタッチして、各部分の点灯表示の増減を行なって使用者の納得状態にする（ステップ１４）。各部分の点灯表示の増減により、メモリ１２７に記憶された各部位毎の快適度のランクが書き換えられる。
【００８９】
そして、納得スイッチ３８をタッチすればステップ１１からステップ１２へ進み、ステップ１２では変更した判定結果、すなわちメモリ１２７に書き換えられた各部位毎の快適度のランクに応じて上記と同様に各エアバック２０〜２８が膨縮されていく。
【００９０】
選択スイッチ３７により自動コースを選択した場合、選択スイッチ３２により全身コースを選択し、オン・オフスイッチ３１をオンすると、予め設定されている順序で各エアバック２０〜２８が膨縮されていく。同様に、上半身,下半身コースを選択すれば上半身または下半身のエアバックが順番に膨縮されていく。
【００９１】
選択スイッチ３７によりお好みコースを選択した場合には、マッサージしたい部分の各エアバック画像Ｇ２０〜Ｇ２８をタッチすれば、このタッチした各エアバック画像Ｇ２０〜Ｇ２８に対応した各エアバック２０〜２８が選択される。そして、オン・オフスイッチ３１をオンすれば、その選択されたエアバックだけが膨縮していく。
【００９２】
上記実施形態では、各部位の快適度に応じて各エアバック２０〜２８の膨張回数を多くしているが、エアバック２０〜２８の膨縮の時間を長くしたり、マッサージの強さを強くしたりしてもよい。また、各部位の快適度に応じて各エアバック２０〜２８で行うマッサージの頻度を多くしてもよく、膨張回数,膨縮時間,強さ,頻度のうち少なくとも２つ以上のものを行うようにしてもよい。
【００９３】
また、上記実施形態では椅子式のマッサージ装置１０について説明したが、マット式のマッサージ装置であってもよく、また、エアバックの替わりに例えばローラ等を使用する機械式のマッサージ装置であってもよい。
【００９４】
上記実施形態では、（５）式からリラックス状態を求めているが、平均心拍周期と心拍変化率と心拍変動度のうち少なくともいずれか１つを使用してリラックス状態を求めてもよい。また、移動平均値からリラックス状態を求めてもよい。
【００９５】
また、リラックスモードと施療モードの２つのモードを設けておき、リラックスモードが設定されたら上述の動作を行い、施療モードが設定されたら上記と逆の動作を行うようにしてもよい。すなわち、快適度が低いほどエアバックの膨張回数や膨縮時間,強さ,頻度等を多くするものである。
【００９６】
これは、例えば肩が凝りすぎている場合、その部分をマッサージすると心地よさを感じずに、逆に痛みなどを感じてしまうからである。すなわち、その患部を施療する場合、快適度が低いほどエアバックの膨張回数等を多くすれば、患部の部位の施療を行うことができることになる。
【００９７】
【発明の効果】
以上説明したように、この発明によれば、自動的に心地よい部位を探してその部位をマッサージすることができる。
【図面の簡単な説明】
【図１】この発明に係る椅子式のマッサージ装置を示した斜視図である。
【図２】図１に示したマッサージ装置の表示部を示した説明図である。
【図３】図１に示したマッサージ装置の給排気装置の構成を示したブロック図である。
【図４】図３に示す制御装置の構成を示したブロック図である。
【図５】脈波波形の一例を示した説明図である。
【図６】心拍周期と直線近似した直線とを示したグラフである。
【図７】マッサージ刺激を１分間継続したときの心拍周期を示したグラフである。
【図８】マッサージ刺激を行ったときの心拍周期を示した他の例のグラフである。
【図９】心拍周期と移動平均値とを示したグラフである。
【図１０】心拍周期の他の求め方を示した説明図である。
【図１１】心拍周期の別な他の求め方を示した説明図である。
【図１２】マッサージ箇所と心拍周期とマッサージを行わなかった場合の心拍周期を示したグラフである。
【図１３】マッサージ装置の動作を示したフロー図である。
【符号の説明】
２０〜２８    エアバック（施療子）
１２２        平均心拍周期算出手段（心拍特徴検出手段）
１２３        心拍変化率算出手段（心拍特徴検出手段）
１２４        心拍変動度算出手段（心拍特徴検出手段）
１２５        快適度算出手段（施療要素変更手段）
１３０        動作制御部
２０１        脈波センサ（心拍検出手段）
【図１】
0004017915.tif 000004
【図２】
0004017915.tif 000005
【図３】
0004017915.tif 000006
【図４】
0004017915.tif 000007
【図５】
0004017915.tif 000008
【図６】
0004017915.tif 000009
【図７】
0004017915.tif 000010
【図８】
0004017915.tif 000011
【図９】
0004017915.tif 000012
【図１０】
0004017915.tif 000013
【図１１】
0004017915.tif 000014
【図１２】
0004017915.tif 000015
【図１３】
0004017915.tif 000016