(57)【要約】【課題】人の片麻痺を簡易に検査することが可能な検査装置を提供する。【解決手段】人の片麻痺を検査する装置が提供される。検査装置は、人の歩行時の時系列の加速度データに基づいて、人の歩行時の位置データの周波数分布を算出する周波数解析部と、ｈ＝１，２，・・・，Ｈ（Ｈは、２以上の整数）に対し、ｈ番目のピーク周波数ｆhからＨ番目のピーク周波数ｆHまでの区間、又は区間とその近傍の区間を積分区間として周波数分布を積分し、積分値ｇ（ｈ）を算出する積分部と、ｈの対数と、ｇ（ｈ）の対数との線形関係を表す直線を導出する直線導出部と、直線の傾き及び切片、並びに、人のストライド周波数ｆ1の３つの要素の中から選択される少なくとも２つの要素の値により特定される点が、少なくとも２つの要素をそれぞれ軸とする空間内において分布する位置に応じて、人が片麻痺か否かを判断する判断部とを備える。【選択図】図１
2020048683.tif 000002
【特許請求の範囲】
【請求項１】
人の片麻痺を検査する検査装置であって、
前記人の歩行時の時系列の加速度データに基づいて、前記人の歩行時の位置データの周波数分布を算出する周波数解析部と、
ｈ＝１，２，・・・，Ｈ（Ｈは、２以上の整数）に対し、ｈ番目のピーク周波数ｆhからＨ番目のピーク周波数ｆHまでの区間、又は前記区間とその近傍の区間を積分区間として前記周波数分布を積分し、積分値ｇ（ｈ）を算出する積分部と、
ｈの対数と、ｇ（ｈ）の対数との線形関係を表す直線を導出する直線導出部と、
前記直線の傾き及び切片、並びに、前記人のストライド周波数ｆ1の３つの要素の中から選択される少なくとも２つの要素の値により特定される点が、前記少なくとも２つの要素をそれぞれ軸とする空間内において分布する位置に応じて、前記人が片麻痺か否かを判断する判断部と
を備える、検査装置。
【請求項２】
前記判断部は、前記切片と前記ストライド周波数ｆ1の値により特定される点が、前記切片と前記ストライド周波数ｆ1をそれぞれ２軸とする平面内において分布する位置に応じて、前記人が片麻痺か否かを判断する、
請求項１に記載の検査装置。
【請求項３】
前記判断部は、前記傾きと前記切片の値により特定される点が、前記傾きと前記切片をそれぞれ２軸とする平面内において分布する位置に応じて、前記人が片麻痺か否かを判断する、
請求項１に記載の検査装置。
【請求項４】
請求項１から３のいずれかに記載の検査装置と、
前記検査装置に内蔵又は外付けされ、前記加速度データを計測する加速度計と
を備える、検査システム。
【請求項５】
人の片麻痺を検査するためのプログラムであって、
前記人の歩行時の時系列の加速度データに基づいて、前記人の歩行時の位置データの周波数分布を算出するステップと、
ｈ＝１，２，・・・，Ｈ（Ｈは、２以上の整数）に対し、ｈ番目のピーク周波数ｆhからＨ番目のピーク周波数ｆHまでの区間、又は前記区間とその近傍の区間を積分区間として前記周波数分布を積分し、積分値ｇ（ｈ）を算出するステップと、
ｈの対数と、ｇ（ｈ）の対数との線形関係を表す直線を導出するステップと、
前記直線の傾き及び切片、並びに、前記人のストライド周波数の３つの要素の中から選択される少なくとも２つの要素の値が、前記少なくとも２つの要素をそれぞれ軸とする空間内において分布する位置に応じて、前記人が片麻痺か否かを判断するステップと
をコンピュータに実行させる、プログラム。
【請求項６】
人の片麻痺を検査する方法であって、
前記人の歩行時の時系列の加速度データに基づいて、前記人の歩行時の位置データの周波数分布を算出するステップと、
ｈ＝１，２，・・・，Ｈ（Ｈは、２以上の整数）に対し、ｈ番目のピーク周波数ｆhからＨ番目のピーク周波数ｆHまでの区間、又は前記区間とその近傍の区間を積分区間として前記周波数分布を積分し、積分値ｇ（ｈ）を算出するステップと、
ｈの対数と、ｇ（ｈ）の対数との線形関係を表す直線を導出するステップと、
前記直線の傾き及び切片、並びに、前記人のストライド周波数の３つの要素の中から選択される少なくとも２つの要素の値が、前記少なくとも２つの要素をそれぞれ軸とする空間内において分布する位置に応じて、前記人が片麻痺か否かを判断するステップと
を含む、方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、人の片麻痺を検査する検査装置、プログラム及び方法に関する。
【背景技術】
【０００２】
脳梗塞や脳出血を含む脳卒中の検査には、通常、ＭＲＩやＣＴが使用される。しかし、このような設備は、病院等の特定の施設にしかなく、このような設備を使用する検査には、お金と時間がかかる。また、ＣＴには、放射能侵襲があるし、ＭＲＩには、磁場侵襲がある。さらに、このような検査は、実際上、毎日受けることはできない。しかし、例えば、脳卒中の中でも近年多く見られる脳梗塞の場合、発症後、４．５時間以内に血栓溶解薬を注射すれば、回復が見込まれることが知られている。よって、脳卒中は、早期発見が重要である。
【０００３】
脳卒中は、片麻痺の症状を伴うことが多くある。片麻痺の回復過程を評価する指標に、ブルンストロームステージと呼ばれる指標がある。これによれば、片麻痺の程度がＩ〜ＶＩの６段階に分けて評価され、数字が大きい程軽微であり、健常者に近い。特に最も軽微なステージＶＩの人の場合、片麻痺が発症していたとしても、周りからは分かりづらく、本人にも自覚がないことがしばしばある。従って、片麻痺の早期発見は難しいのが現状である。
【０００４】
特許文献１は、片麻痺の患者及び健常者を含む被験者の腰に加速度計を装着して腰の軌道パターンを計測し、計測された腰の軌道パターンの幾何学的特徴とブルンストロームステージとの関係性を評価している。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】西辰徳、他２名，「腰軌道の運動学的分析に基づく片麻痺歩行評価システム」，計測自動制御学会論文集，第４７巻，第１号，８−１６頁，２０１１年１月
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明は、人の片麻痺を簡易に検査することが可能な検査装置、プログラム及び方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
第１観点に係る検査装置は、人の片麻痺を検査する装置であって、前記人の歩行時の時系列の加速度データに基づいて、前記人の歩行時の位置データの周波数分布を算出する周波数解析部と、ｈ＝１，２，・・・，Ｈ（Ｈは、２以上の整数）に対し、ｈ番目のピーク周波数ｆhからＨ番目のピーク周波数ｆHまでの区間、又は前記区間とその近傍の区間を積分区間として前記周波数分布を積分し、積分値ｇ（ｈ）を算出する積分部と、ｈの対数と、ｇ（ｈ）の対数との線形関係を表す直線を導出する直線導出部と、前記直線の傾き及び切片、並びに、前記人のストライド周波数ｆ1の３つの要素の中から選択される少なくとも２つの要素の値により特定される点が、前記少なくとも２つの要素をそれぞれ軸とする空間内において分布する位置に応じて、前記人が片麻痺か否かを判断する判断部とを備える。
【０００８】
第２観点に係る検査装置は、第１観点に係る検査装置であって、前記判断部は、前記切片と前記ストライド周波数ｆ1の値により特定される点が、前記切片と前記ストライド周波数ｆ1をそれぞれ２軸とする平面内において分布する位置に応じて、前記人が片麻痺か否かを判断する。
【０００９】
第３観点に係る検査装置は、第１観点に係る検査装置であって、前記判断部は、前記傾きと前記切片の値により特定される点が、前記傾きと前記切片をそれぞれ２軸とする平面内において分布する位置に応じて、前記人が片麻痺か否かを判断する。
【００１０】
第４観点に係る検査システムは、第１観点から第３観点のいずれかに係る検査装置と、前記検査装置に内蔵又は外付けされ、前記加速度データを計測する加速度計とを備える。
【００１１】
第５観点に係るプログラムは、人の片麻痺を検査するためのプログラムであって、以下のステップをコンピュータに実行させる。
（１）前記人の歩行時の時系列の加速度データに基づいて、前記人の歩行時の位置データの周波数分布を算出するステップ。
（２）ｈ＝１，２，・・・，Ｈ（Ｈは、２以上の整数）に対し、ｈ番目のピーク周波数ｆhからＨ番目のピーク周波数ｆHまでの区間、又は前記区間とその近傍の区間を積分区間として前記周波数分布を積分し、積分値ｇ（ｈ）を算出するステップ。
（３）ｈの対数と、ｇ（ｈ）の対数との線形関係を表す直線を導出するステップ。
（４）前記直線の傾き及び切片、並びに、前記人のストライド周波数の３つの要素の中から選択される少なくとも２つの要素の値が、前記少なくとも２つの要素をそれぞれ軸とする空間内において分布する位置に応じて、前記人が片麻痺か否かを判断するステップ。
【００１２】
第６観点に係る方法は、人の片麻痺を検査するための方法であって、以下のステップを含む。
（１）前記人の歩行時の時系列の加速度データに基づいて、前記人の歩行時の位置データの周波数分布を算出するステップ。
（２）ｈ＝１，２，・・・，Ｈ（Ｈは、２以上の整数）に対し、ｈ番目のピーク周波数ｆhからＨ番目のピーク周波数ｆHまでの区間、又は前記区間とその近傍の区間を積分区間として前記周波数分布を積分し、積分値ｇ（ｈ）を算出するステップ。
（３）ｈの対数と、ｇ（ｈ）の対数との線形関係を表す直線を導出するステップ。
（４）前記直線の傾き及び切片、並びに、前記人のストライド周波数の３つの要素の中から選択される少なくとも２つの要素の値が、前記少なくとも２つの要素をそれぞれ軸とする空間内において分布する位置に応じて、前記人が片麻痺か否かを判断するステップ。
【発明の効果】
【００１３】
本発明によれば、人の片麻痺を簡易に検査することができる。その結果、例えば、ＭＲＩやＣＴ等の設備に頼らずとも、脳卒中の早期発見を助けることができる。
【図面の簡単な説明】
【００１４】
【図１】本発明の一実施形態に係る片麻痺の検査装置の構成を示す図。
【図２】本発明の一実施形態に係る片麻痺の検査方法の流れを示すフローチャート。
【図３Ａ】健常者の加速度データを示すグラフ。
【図３Ｂ】股関節異常患者の加速度データを示すグラフ。
【図３Ｃ】片麻痺の脳卒中患者の加速度データを示すグラフ。
【図４Ａ】健常者の加速度データの周波数分布ｐ（ｆ）を示すグラフ。
【図４Ｂ】股関節異常患者の加速度データの周波数分布ｐ（ｆ）を示すグラフ。
【図４Ｃ】片麻痺の脳卒中患者の加速度データの周波数分布ｐ（ｆ）を示すグラフ。
【図５Ａ】健常者の位置データの周波数分布の様々な積分区間での積分値ｇ（ｈ）を示すグラフ。
【図５Ｂ】股関節異常患者の位置データの周波数分布の様々な積分区間での積分値ｇ（ｈ）を示すグラフ。
【図５Ｃ】片麻痺の脳卒中患者の位置データの周波数分布の様々な積分区間での積分値ｇ（ｈ）を示すグラフ。
【図６Ａ】健常者についてのｈの対数と、積分値ｇ（ｈ）の対数との関係を示すグラフ。
【図６Ｂ】股関節異常患者についてのｈの対数と、積分値ｇ（ｈ）の対数との関係を示すグラフ。
【図６Ｃ】片麻痺の脳卒中患者についてのｈの対数と、積分値ｇ（ｈ）の対数との関係を示すグラフ。
【図７】健常者、股関節異常患者、並びにブルンストロームステージＩＶ、Ｖ及びＶＩの脳卒中患者を含む４１名についての、切片ｃiとストライド周波数ｆ1のプロット図。
【図８】健常者、股関節異常患者、並びにブルンストロームステージＩＶ、Ｖ及びＶＩの脳卒中患者を含む４１名についての、傾きｓiと切片ｃiのプロット図。
【図９】健常者、股関節異常患者、並びにブルンストロームステージＩＶ、Ｖ及びＶＩの脳卒中患者を含む４１名についての、傾きｓiとストライド周波数ｆ1のプロット図。
【発明を実施するための形態】
【００１５】
以下、図面を参照しつつ、本発明の一実施形態に係る片麻痺の検査装置、プログラム及び方法について説明する。
【００１６】
＜１．片麻痺の検査装置の構成＞
図１に、本実施形態に係る片麻痺の検査装置１の構成図を示す。検査装置１は、ハードウェアとしては、汎用のコンピュータであり、これに限定されないが、本実施形態では、加速度計３が内蔵されたスマートフォンから構成される。スマートフォンは、インターネットを介して所定のプログラム２をダウンロードし、インストールすることにより、片麻痺の検査を行う検査装置１として使用可能となる。プログラム２は、加速度計３により計測される人Ｐ１の歩行時の時系列の加速度データに基づいて、人Ｐ１の片麻痺を簡易に検査するアプリケーションソフトウェアであり、検査装置１に後述する処理を実行させる。検査装置１によれば、ＭＲＩやＣＴ等の大掛かりな設備に頼らずとも、簡易に片麻痺の検査を行うことができ、脳卒中の早期発見を助けることができる。検査装置１による片麻痺の検査は、医師や理学療法士等の専門家とともに行うこともできるが、自分だけで行うこともできるため、日常的な検査が可能になる。
【００１７】
図１に示すとおり、検査装置１は、表示部１０、入力部２０、記憶部３０、制御部４０及び通信部５０を有する。これらの部１０〜５０は、互いにバス線７を介して接続されている。表示部１０は、液晶ディスプレイ等で構成され、適当な画面をユーザーに対し表示する。入力部２０は、タッチパネルや操作ボタン等で構成され、検査装置１に対するユーザーからの操作を受け付ける。通信部５０は、検査装置１をインターネットや、BLUETOOTH（登録商標）のような近距離無線通信ネットワーク等に通信接続する通信インターフェースである。記憶部３０は、フラッシュメモリ等の不揮発性の記憶装置で構成され、記憶部３０内には、プログラム２が格納されている。
【００１８】
制御部４０は、ＣＰＵ、ＲＯＭ及びＲＡＭ等から構成されており、記憶部３０内に格納されているプログラム２を読み出して実行することにより、仮想的に、周波数解析部４１、積分部４２、直線導出部４３、判断部４４及び結果出力部４５として動作する。各部４１〜４５の動作については、後述する。
【００１９】
＜２．片麻痺の検査方法の流れ＞
以下、図２を参照しつつ、本実施形態に係る片麻痺の検査方法の流れについて説明する。まず、片麻痺の検査を行うに当たり、ユーザーは、入力部２０を操作して、検査装置１上でプログラム２を起動させる。これにより、図２に示す処理が開始する。なお、ここでいうユーザーとは、被験者となる人Ｐ１自身の他、医師や理学療法士等の専門家のような、人Ｐ１の片麻痺の検査結果を必要とする者の総称である。
【００２０】
次に、被験者となる人Ｐ１の身体の所定の位置、本実施形態では、図１に示すように、腰の後ろの仙骨に近い位置に、検査装置１を図視されないベルト等を介して装着する。ただし、検査装置１の装着位置はこれに限定されず、首の後ろ等、その他の部位に装着することも可能である。そして、この状態で、人Ｐ１がテスト歩行を行うと、検査装置１に含まれる加速度計３により、人Ｐ１の歩行時の時系列の加速度データが収集される（ステップＳ１）。テスト歩行では、人Ｐ１に８歩以上歩かせることが好ましい。また、歩き始めと歩き終わりの加速度データは、定常的な歩行の状態を表すデータとはならないため、後述する解析には必ずしも適さない。従って、加速度計３により収集された加速度データのうち、歩き始めと歩き終わりのデータを削除し、残りの一定速度のデータを後述する解析の対象とすることが好ましい。８歩分の加速度データであれば、通常は４秒程度の時間で収集することができ、６歩分の一定速度のデータが含まれる。
【００２１】
加速度計３により計測された加速度データは、記憶部３０内に保存される。なお、加速度計３は、３軸加速度センサであり、ここでは、左右、前後及び上下方向の３方向の加速度データが収集される。左右、前後及び上下方向は、歩行中の人Ｐ１を基準に定義される。図３Ａ〜図３Ｃは、それぞれ、健常者、股関節に異常を有する患者、及び脳卒中による片麻痺の患者の実際の加速度データのグラフを示している。
【００２２】
続いて、周波数解析部４１が、記憶部３０内の加速度データに基づいて、人Ｐ１の歩行時の位置データ（加速度計３が装着された部位の変位データ）の周波数分布を算出する（ステップＳ２）。より具体的には、周波数解析部４１は、３方向の加速度データをそれぞれフーリエ変換し、各方向について、周波数スペクトルのコサイン成分（実部）及びサイン成分（虚部）を導出する。続いて、周波数解析部４１は、これら３方向×２つの成分の合計６つの成分を６次元ベクトルとしたときのベクトル長（スカラー値、大きさ）ｐを算出する。ｐは、周波数ｆの関数であり、以下、ｐ（ｆ）と表現する。図４Ａ〜図４Ｃは、それぞれ、図３Ａ〜図３Ｃの加速度データに基づくｐ（ｆ）のグラフである。
【００２３】
続いて、周波数解析部４１は、加速度データの周波数分布の大きさであるｐ（ｆ）を、位置データの周波数分布の大きさｑ（ｆ）へと変換する。変換式は、ｑ（ｆ）＝−ｐ（ｆ）／（２πｆ）2である。すなわち、周波数ｆのサイン関数ｓｉｎ（２πｆｔ）及びコサイン関数ｃｏｓ（２πｆｔ）は、時間ｔで二回微分すると−（２πｆ）2倍になるため、−１／（２πｆ）2倍すると、二回積分したことになる。よって、以上の変換式に基づくｑ（ｆ）は、加速度計３が装着されている腰の位置を表す位置ベクトル（３次元ベクトル）の時系列データの周波数分布（６次元ベクトル）のベクトル長（スカラー値、大きさ）を表す。
【００２４】
続いて、積分部４２は、ｈ＝１，２，・・・，Ｈ（Ｈは、２以上の整数）に対し、ｈ番目のピーク周波数ｆhの近傍からＨ番目のピーク周波数ｆHの近傍までを積分区間として、周波数分布の大きさｑ（ｆ）を積分し、積分値ｇ（ｈ）を算出する（ステップＳ３）。ｇ（ｈ）は、以下の通り、定義される。本実施形態では、積分区間は、ピーク周波数ｆhからピーク周波数ｆHまでの区間よりもやや広く、同区間と同区間の近傍の区間とを含む、（ｈ−０．５）×ｆ1から（Ｈ＋０．５）×ｆ1までの区間である。また、本実施形態では、Ｈ＝６である。
【数１】
2020048683.tif 000003
【００２５】
図４Ａ〜図４Ｃの１Ｈｚ当たりにある１番目のピークは、ストライド（右一歩と左一歩）の周波数（以下、ストライド周波数という）ｆ1のピークである。人間の歩行には、通常、１Ｈｚの近傍にストライド周波数ｆ1があり、そのパワーは、左右の足で二足歩行する人間の横揺れに対応する。そして、その２倍の周波数２ｆ1（以下、ステップ周波数という）は、一歩ずつの人間の縦揺れと前後揺れに対応する。よく知られているように、一定の弾力係数を有する物体の自由振動は、ある周波数のサイン・コサイン関数になる。人間の歩行は、自由振動のように滑らかであれば、横揺れに対応するストライド周波数ｆ1と、縦揺れ及び前後揺れに対応するステップ周波数２ｆ1とのサイン・コサイン関数だけで表されるはずである。しかし、実際には、３ｆ1以上の領域にも、パワーが存在する。歩行時には、足を床から離すトゥキックや、踵を床に接触させるヒールタッチ等の衝撃が足に加わり、それらの衝撃が３ｆ1以上の領域に現れていると考えられる。そして、これらの衝撃は、健常者の場合、膝や腰でうまく吸収して軽減されると考えられるが、片麻痺の患者の場合、そうはいかない。発明者らは、以上のとおり定義されるｇ（ｈ）が、ｈ×ｆ1以上の周波数において吸収されなかった衝撃のパワーを表すものと考えた。図５Ａ〜図５Ｃは、それぞれ、図４Ａ〜図４Ｃのｐ（ｆ）に基づくｇ（ｈ）のグラフである。
【００２６】
続いて、直線導出部４３は、ｈ＝１，２，・・・，Ｈに対し、ｈの対数ｌｏｇ（ｈ）と、ｇ（ｈ）の対数ｌｏｇ（ｇ（ｈ））を算出し、ｌｏｇ（ｈ）とｌｏｇ（ｇ（ｈ））との線形関係を表す近似直線を導出する（ステップＳ４）。近似直線は、最小二乗法等により算出される。図６Ａ〜図６Ｃは、それぞれ、図５Ａ〜図５Ｃの（ｈ，ｇ（ｈ））に基づく（ｌｏｇ（ｈ），ｌｏｇ（ｇ（ｈ）））の点をプロットしたグラフである。
【００２７】
自然の海や湖の波の周期や、風や小川のせせらぎには、１／ｆ揺らぎが現れる。１／ｆ揺らぎとは、パワーが周波数ｆに反比例する揺らぎのことであり、１／ｆ揺らぎが現れるとき、周波数ｆの対数ｌｏｇ（ｆ）と、その周波数分布ａ（ｆ）の対数ｌｏｇ（ａ（ｆ））との関係は、傾きが−１の直線となる。このことに着想を得て、本発明者らは、ｈ＝１，２，・・・，Ｈに対するＨ個の（ｌｏｇ（ｈ），ｌｏｇ（ｇ（ｈ）））の点は、１本の直線に回帰すると予測し、このことを実験により確認した。この実験は、健常者２０名と、股関節異常患者１名と、片麻痺の脳卒中患者２０名（うち、ブルンストロームステージＩＶ、Ｖ、ＶＩの人がそれぞれ７名、６名、７名）の合計４１名に対して行った。図６Ａ〜図６Ｃは、このうち、健常者１名と、股関節異常患者１名と、ブルンストロームステージＶＩの片麻痺の脳卒中患者１名分のデータを表している。なお、ステージＶＩの患者は、歩行を観察したときに健常者と余り区別がつかず、ステージＶの患者は、健常者とは区別がつくが、自立して歩けるレベルであり、ステージＩＶの患者は、かろうじて歩けるレベルである。ステージＩＩＩ以下は、自立して歩くのが困難なレベルである。この実験においては、以上の４１名の各人に対し、（ｌｏｇ（ｈ），ｌｏｇ（ｇ（ｈ）））の関係について、直線に近いほど絶対値が１に近くなる相関係数を算出した。その結果、４１名のこのような相関係数は、−１〜−０．９６の範囲に収まった。なお、このときの近似直線の傾きは−３に近くなったため、人間の歩行は、１／ｆ3揺らぎであると考察される。
【００２８】
続いて、判断部４４は、ステップＳ４で導出された近似直線の切片ｃiの値と、ステップＳ２で導出されているストライド周波数ｆ1の値により特定される（ｃi，ｆ1）の点が、ｃi−ｆ1平面内において分布する位置に応じて、人Ｐ１が片麻痺か否かを判断する（ステップＳ５）。ｃi−ｆ1平面は、切片ｉ及びストライド周波数ｆ1をそれぞれ２軸とする平面領域である。図７は、上記した実験による４１名の（ｃi，ｆ1）の点の分布図である。図中、黒丸（normal）の点は、健常者のデータを表し、「＋」（ＨｉｐＯＡ）の点は、股関節異常患者のデータを表し、三角、×、四角の点は、それぞれブルンストロームステージＩＶ、Ｖ、ＶＩの片麻痺の脳卒中患者を表している。同図に示すように、片麻痺の脳卒中患者に対応する（ｃi，ｆ1）の点が分布する領域Ａ１は、ｃi−ｆ1平面内において健常者及び股関節異常患者のそれが分布する領域Ａ２から区別される。図７には、これらの領域の境界線Ｌ１が破線で示されている。境界線Ｌ１は、直線として定義することができる。判断部４４は、予め定められ、記憶部３０内に格納されている境界線Ｌ１を特定する情報に基づいて、被験者である人Ｐ１に対応する（ｃi，ｆ1）の点が、領域Ａ１又はＡ２内のいずれに存在するかを判断する。領域Ａ１内に存在する場合には、人Ｐ１に片麻痺の症状がでていると判断し、領域Ａ２内に存在する場合には、人Ｐ１に片麻痺の症状が出ていないと判断する。
【００２９】
続いて、結果出力部４５は、ステップＳ５の判断結果を示す画面を作成し、表示部１０上に表示させる（ステップＳ６）。これを見たユーザーは、人Ｐ１に片麻痺の症状が出ているかを判断し、必要に応じて必要な治療に早期に取り掛かることができる。
【００３０】
＜３．変形例＞
以上、本発明の一実施形態について説明したが、本発明は上記実施形態に限定されるものではなく、その趣旨を逸脱しない限りにおいて、種々の変更が可能である。例えば、以下の変更が可能である。また、以下の変形例の要旨は、適宜組み合わせることができる。
【００３１】
＜３−１＞
上記実施形態において、検査装置１は、スマートフォンとして実現された。しかしながら、検査装置１のハードウェア構成は、この態様に限定されず、例えば、加速度計３を万歩計（登録商標）のような専用の装置とし、これに検査装置１としてのコンピュータを外付けで接続することにより、検査システムを構成してもよい。この構成では、加速度計３により計測した加速度データを、無線又は有線通信によりコンピュータに送信し、コンピュータにおいて上記実施形態と同様に加速度データを解析すればよい。ここでいうコンピュータは、スマートフォンであってもよいし、ラップトップコンピュータ、デスクトップコンピュータ等、任意のコンピュータとすることができる。また、加速度計３、プロセッサ及び表示部を１つの筐体内に搭載した万歩計（登録商標）のような専用の装置として、検査装置１を実現することもできる。
【００３２】
＜３−２＞
上記実施形態では、切片ｃiとストライド周波数ｆ1の値により特定される（ｃi，ｆ1）の点の、ｃi−ｆ1平面内における位置に応じて、片麻痺が判断された。しかしながら、これに限らず、ステップＳ４で導出されるｌｏｇ（ｈ）とｌｏｇ（ｇ（ｈ））との線形関係を表す近似直線の傾きｓi及び切片ｃi、並びに、ストライド周波数ｆ1の３つの要素の中から選択される少なくとも２つの要素の値により特定される点が、これらの要素をそれぞれ軸とする空間内において分布する位置に応じて、片麻痺を判断してもよい。
【００３３】
図８は、上記した実験による４１名の（ｓi，ｃ1）の点の分布図である。同図に示すように、片麻痺の脳卒中患者に対応する（ｓi，ｃ1）の点が分布する領域Ｂ１も、ｓi−ｃ1平面内において健常者及び股関節異常患者のそれが分布する領域Ｂ２から区別される。図８には、これらの領域の境界線Ｌ２が破線で示されている。境界線Ｌ２は、直線として定義することができる。よって、判断部４４は、予め定められ、記憶部３０内に格納されている境界線Ｌ２を特定する情報に基づいて、被験者である人Ｐ１に対応する（ｓi，ｃ1）の点が領域Ｂ１又はＢ２内のいずれに存在するかを判断してもよい。そして、領域Ｂ１内に存在する場合には、人Ｐ１に片麻痺の症状がでていると判断し、領域Ｂ２内に存在する場合には、人Ｐ１に片麻痺の症状が出ていないと判断してもよい。
【００３４】
同様に、図９は、上記した実験による４１名の（ｓi，ｆ1）の点の分布図である。同図に示すように、片麻痺の脳卒中患者に対応する（ｓi，ｆ1）の点が分布する領域Ｃ１も、ｓi−ｆ1平面内において健常者及び股関節異常患者のそれが分布する領域Ｃ２から概ね区別される。図９には、これらの領域の境界線Ｌ３が破線で示されている。境界線Ｌ３は、直線として定義することができる。よって、判断部４４は、予め定められ、記憶部３０内に格納されている境界線Ｌ３を特定する情報に基づいて、被験者である人Ｐ１に対応する（ｓi，ｆ1）の点が領域Ｃ１又はＣ２内のいずれに存在するかを判断してもよい。そして、領域Ｃ１内に存在する場合には、人Ｐ１に片麻痺の症状がでていると判断し、領域Ｃ２内に存在する場合には、人Ｐ１に片麻痺の症状が出ていないと判断してもよい。
【００３５】
また、以上のことから分かるように、片麻痺の脳卒中患者に対応する（ｓi，ｃ1，ｆ1）の点が分布する領域も、ｓi−ｃ1−ｆ1空間内において健常者及び股関節異常患者のそれが分布する領域から区別される。よって、判断部４４は、予め定められ、記憶部３０内に格納されているこれらの領域を区別する境界を特定する情報に従って、被験者である人Ｐ１に対応する（ｓi，ｃ1，ｆ1）の点がこれらの領域のいずれに存在するかを判断し、片麻痺の症状が出ているか否かを判断してもよい。なお、境界は、平面として定義することができる。
【符号の説明】
【００３６】
１  検査装置
２  プログラム
３  加速度計
４１  周波数解析部
４２  積分部
４３  直線導出部
４４  判断部
ｃi  切片
ｓi  傾き
ｆ1  ストライド周波数
【図１】
2020048683.tif 000004
【図２】
2020048683.tif 000005
【図３Ａ】
2020048683.tif 000006
【図３Ｂ】
2020048683.tif 000007
【図３Ｃ】
2020048683.tif 000008
【図４Ａ】
2020048683.tif 000009
【図４Ｂ】
2020048683.tif 000010
【図４Ｃ】
2020048683.tif 000011
【図５Ａ】
2020048683.tif 000012
【図５Ｂ】
2020048683.tif 000013
【図５Ｃ】
2020048683.tif 000014
【図６Ａ】
2020048683.tif 000015
【図６Ｂ】
2020048683.tif 000016
【図６Ｃ】
2020048683.tif 000017
【図７】
2020048683.tif 000018
【図８】
2020048683.tif 000019
【図９】
2020048683.tif 000020