(57)【要約】【課題】フィルタ処理技術では除去できないノイズ成分を除去できるようにし、これにより解析精度のさらなる向上を図る。【解決手段】左右の各側頭筋の表面筋電図（ｓＥＭＧ）の計測データを取得し、これらの計測データに対し２０−４５０Hzのみを通過させ他の周波数成分を除去するバンドパスフィルタ処理を行った後、特徴量を示すＲＭＳを算出する。そして、上記左右の側頭筋に対応するＲＭＳ間の相関係数（Ｒ）を算出し、この相関係数（Ｒ）がしきい値以下の場合に、上記ＲＭＳには側頭筋の運動に寄らないノイズ成分が多く残留していると見做して、上記ＲＭＳの値を事前に設定しておいたベースライン値に置換する。【選択図】図２
2019030520.tif 000002
【特許請求の範囲】
【請求項１】
連動する第１および第２の筋肉部位の表面に設置された各電極から、それぞれ前記第１および第２の筋肉部位による表面筋電図を示す第１および第２の計測信号を取得する取得部と、
前記第１および第２の計測信号に基づいてそれぞれ前記表面筋電図の特徴量を算出する特徴量算出部と、
前記算出された各特徴量間の相関の度合いを示す相関値を算出する相関値算出部と、
前記算出された相関値に基づいて前記各特徴量を補正する補正部と
を具備する筋電計測装置。
【請求項２】
前記特徴量算出部は、
前記第１および第２の計測信号からそれぞれ予め設定した周波数帯域の信号を抽出するフィルタ部と、
前記フィルタ部により抽出された各信号からそれぞれ前記特徴量を算出する算出部と
を備える、請求項１に記載の筋電計測装置。
【請求項３】
前記補正部は、
前記各特徴量を補正するための補正情報を記憶する記憶媒体と、
前記算出された相関値が予め設定したしきい値以下の場合に、前記算出された各特徴量を前記記憶された補正情報に基づいて補正する補正処理部と、
を備える、請求項１または２に記載の筋電計測装置。
【請求項４】
前記補正部は、
前記相関値が最大値をとるときの前記各特徴量を求め、当該各特徴量を前記補正情報として前記記憶媒体に記憶させる補正情報生成部を、
さらに備える、請求項３に記載の筋電計測装置。
【請求項５】
筋電計測装置が実行する筋電計測方法であって、
連動する第１および第２の筋肉部位の表面に設置された各電極から、それぞれ前記第１および第２の筋肉部位による表面筋電図を示す第１および第２の計測信号を取得する過程と、
前記第１および第２の計測信号に基づいてそれぞれ前記表面筋電図の特徴量を算出する過程と、
前記算出された各特徴量間の相関の度合いを示す相関値を算出する過程と、
前記算出された相関値に基づいて前記各特徴量を補正する過程と、
を具備する筋電計測方法。
【請求項６】
前記補正する過程は、前記算出された相関値が予め設定したしきい値以下の場合に、前記算出された各特徴量を、予め記憶媒体に記憶された補正情報に基づいて補正する、
請求項５記載の筋電計測方法。
【請求項７】
前記補正する過程は、前記相関値が最大値をとるときの前記各特徴量を求め、当該各特徴量を前記補正情報として前記記憶媒体に記憶させる過程を、
さらに備える、請求項６記載の筋電計測方法。
【請求項８】
請求項１乃至４のいずれか１項に記載の筋電計測装置が具備する前記各部としてプロセッサを機能させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、表面筋電図を計測するための筋電計測装置、方法及びプログラムに関する。
【背景技術】
【０００２】
近年、医学的研究分野およびリハビリテーションの臨床のほか、スポーツ科学や産業医学などの人間光学的分野では、筋活動を観察するために表面筋電図（surface Electromyography：ｓＥＭＧ）が用いられている。表面筋電図は、使い捨て電極（ディスポーザブル電極）や乾式電極を、観察対象の筋に近い皮膚に設置するだけで、比較的簡易に計測することが可能である。
【０００３】
筋電図の解析においては、計測信号を二乗平均平方根（Root Mean Square：ＲＭＳ）などの特徴量に変換することが一般的である。しかし、表面筋電図の計測において、体動などにより電極が浮いたりずれたりして皮膚と電極の接触状況が変化すると、その変化がノイズ成分となって計測値に含まれ、ＲＭＳ値が変化する。そのため、観察対象の筋活動を解析する際には、筋活動による電位変化とノイズ成分による電位変化とを分離して解析する必要がある。
【０００４】
そこで従来では、信号処理技術として一般的に用いられているフィルタ処理技術により、表面筋電図の計測値からノイズ成分を除去する手法が提案されている。例えば、非特許文献１では、通過周波数帯域が２０〜４５０Hzのバンドパスフィルタを用いて、体動ノイズ成分を除去する手法が提案されている。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】Carlo J. De Luca et al. “Filtering the surface EMG signal Movement artifact and baseline noise contamination”, Journal of Biomechanics, Volume 43, Issue 8, Pages 1573−1579, (2010)
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ところが、バンドパスフィルタを用いても、２０〜４５０Hz帯に含まれるノイズ成分については除去することができないため、解析精度のさらなる向上は期待できない。
【０００７】
この発明は上記事情に着目してなされたもので、その目的とするところは、フィルタ処理技術では除去できないノイズ成分を除去できるようにし、これにより解析精度のさらなる向上を図った筋電計測装置、方法およびプログラムを提供することにある。
【課題を解決するための手段】
【０００８】
上記目的を達成するために、この発明に係る筋電計測装置または方法の第１の態様は、連動する第１および第２の筋肉部位の表面に設置された各電極から、それぞれ前記第１および第２の筋肉部位による表面筋電図を示す第１および第２の計測信号を取得し、前記第１および第２の計測信号に基づいてそれぞれ前記表面筋電図の特徴量を算出し、前記算出された各特徴量間の相関の度合いを示す相関値を算出し、前記算出された相関値に基づいて前記各特徴量を補正するようにしたものである。
【発明の効果】
【０００９】
この発明の第１の態様によれば、以下のような作用効果が奏せられる。すなわち、連動する第１および第２の筋肉部位において計測された表面筋電図の各特徴量は、筋活動による成分のみを含む場合には略一致するため、その相関値は最大値をとる。これに対し、例えば筋活動によらない皮膚の動きや電極ずれ等が生じると、上記各特徴量にノイズ成分が混入し、相関値は低下する。そこで、上記各特徴量間の相関値に基づいて特徴量を補正することで、フィルタ処理技術のみでは除去できないノイズ成分を除去することが可能となり、これにより解析精度のさらなる向上を図ることができる。
【図面の簡単な説明】
【００１０】
【図１】図１は、この発明の一実施形態に係る筋電計測装置を備えたシステムの全体構成図である。
【図２】図２は、この発明の一実施形態に係る筋電計測装置の機能構成を示すブロック図である。
【図３】図３は、図２に示した筋電計測装置による筋電計測制御の手順と処理内容を示すフローチャートである。
【図４】咀嚼運動をしたときの側頭筋の筋活動を示すＲＭＳと相関係数Ｒの一例を示す図である。
【図５】歯を食いしばったときの側頭筋の筋活動を示すＲＭＳと相関係数Ｒの一例を示す図である。
【図６】瞬きをしたときの側頭筋の筋活動を示すＲＭＳと相関係数Ｒの一例を示す図である。
【図７】帽子の鍔を左右に動かしたときの側頭筋の筋活動を示すＲＭＳと相関係数Ｒの一例を示す図である。
【発明を実施するための形態】
【００１１】
以下、図面を参照してこの発明に係わる実施形態を説明する。  
［一実施形態］
（構成）
図１は、この発明の一実施形態に係る筋電計測装置を備えたシステムの全体構成図である。  
図１に示すように、ユーザの左側頭筋上および右側頭筋上の各皮膚面にはそれぞれ電極２Ｌ，２Ｒが貼付される。左側頭筋および右側頭筋は、例えば咀嚼運動を行うための筋肉の１つであり、こめかみ付近に位置している。
【００１２】
電極２Ｌ，２Ｒは、例えば電極サイズが３０mm×３０mmの布電極からなり、それぞれ左側頭筋および右側頭筋の筋活動を表す表面筋電図（ｓＥＭＧ）を、計測信号ＩＳL ,ＩＳRとして出力する。これらの計測信号ＩＳL ,ＩＳR は、例えば信号ケーブルを介して筋電計測装置１に送られる。
【００１３】
なお、電極２Ｌ，２Ｒに、例えばBluetooth（登録商標）等の小電力無線データ通信規格を採用した無線ユニットを付設し、これにより計測信号ＩＳL ,ＩＳR を筋電計測装置１に対し無線送信するように構成してもよい。また、電極２Ｌ，２Ｒの貼付位置は、左右の側頭筋上であれば任意に選択することができ、さらに電極２Ｌ，２Ｒのサイズや形状、材質、貼付手段等についても、表面筋電図を計測可能なものであれば如何なるものであってもよい。
【００１４】
筋電計測装置１は、例えばパーソナルコンピュータを用いたもので、以下のように構成される。図２はそのハードウェアおよびソフトウェアの機能構成を示すブロック図である。筋電計測装置１は、制御ユニット１０と、記憶ユニット２０と、入出力インタフェースユニット３０とを備える。
【００１５】
入出力インタフェースユニット３０は、上記電極２Ｌ，２Ｒから出力された表面筋電図（ｓＥＭＧ）を表す計測信号ＩＳL ,ＩＳR を、制御ユニット１０で取り扱うことが可能なディジタルデータに変換する機能を備える。また入出力インタフェースユニット３０は、制御ユニット１０から出力された補正後の計測データを外部デバイスへ出力する機能も備える。外部デバイスとしては、例えば表示器を備えたユーザ端末や、計測データをもとにユーザの咀嚼運動等を判定するパーソナルコンピュータやサーバ等の情報処理装置が想定される。
【００１６】
記憶ユニット２０は、例えば、記憶媒体としてＨＤＤ（Hard Disk Drive）またはＳＳＤ（Solid State Drive）等の随時書込および読出しが可能な不揮発性メモリとＲＡＭ等の揮発メモリを使用したもので、その記憶領域には、プログラム記憶部に加え、計測データ記憶部２１およびベースライン値記憶部２２が設けられている。
【００１７】
計測データ記憶部２１は、電極２Ｌ，２Ｒにより検出された左右の各側頭筋の表面筋電図を表す計測データを記憶するために使用される。ベースライン値記憶部２２は、上記表面筋電図（ｓＥＭＧ）を表す計測データにノイズ成分が含まれている場合に、当該計測データを補正するために使用するベースライン値を記憶するために使用される。
【００１８】
制御ユニット１０は、少なくとも１個のプロセッサと作業用メモリを有し、本実施形態を実施するために必要な処理機能として、ｓＥＭＧ計測データ取得制御部１１と、ＢＰＦ処理部１２と、ＲＭＳ算出部１３と、ベースライン値算出部１４と、相関係数算出部１５と、補正処理部１６とを備えている。これらの処理機能は、いずれも上記記憶ユニット２０に格納されたプログラムを上記プロセッサに実行させることにより実現される。
【００１９】
ｓＥＭＧ計測データ取得制御部１１は、上記入出力インタフェースユニット３０から、上記左右の各側頭筋の動きを表す表面筋電図（ｓＥＭＧ）に対応する計測データを所定時間分取り込み、これらの計測データを時系列に従い計測データ記憶部２１に格納するための制御を行う。
【００２０】
ＢＰＦ処理部１２は、体動等によるノイズ成分を除去するために周波数通過帯域が予め設定されたバンドパスフィルタ処理機能を有する。そして、上記計測データ記憶部２１から読み出された計測データに対し上記バンドパスフィルタ処理を行うことで、上記計測データから体動等によるノイズ成分を除去する。
【００２１】
ＲＭＳ算出部１３は、上記ＢＰＦ処理部１２によりフィルタリング処理された計測データに対し所定の時間幅を有する窓を設定し、この窓ごとに計測データの特徴量を算出する処理を行う。特徴量としては、例えば二乗平均平方根（ＲＭＳ）が用いられる。
【００２２】
ベースライン値算出部１４は、キャリブレーションフェーズにおいて、上記ＲＭＳ算出部１３により窓ごとに算出された計測データの特徴量（ＲＭＳ）の所定時間分の平均値を算出する。そして、この算出したＲＭＳの平均値をベースライン値としてベースライン値記憶部２２に記憶させる処理を行う。
【００２３】
相関係数算出部１５は、上記ベースライン値設定後の通常計測モードにおいて、上記ＲＭＳ算出部１３により窓ごとに算出された、左右の各側頭筋の表面筋電図（ｓＥＭＧ）のＲＭＳ間の相関係数（Ｒ）を算出する処理を行う。
【００２４】
補正処理部１６は、上記相関係数算出部１５により算出された相関係数（Ｒ）を予め設定されたしきい値と比較し、この比較の結果相関係数（Ｒ）がしきい値以下の場合に、このときのＲＭＳにはノイズ成分がまだ多く残留していると判断する。そして、当該ＲＭＳの値を上記ベースライン値記憶部２２に記憶されているベースライン値に基づいて補正する処理を行う。
【００２５】
（動作）
次に、以上のように構成された筋電計測装置１による表面筋電図の計測動作を説明する。  
（１）ベースライン値の設定
実際の表面筋電図の計測に先立ち、筋電計測装置１は、ＲＭＳの補正に使用するベースライン値を定めるキャリブレーション処理を以下のように実行する。
【００２６】
先ず図１に例示したように、電極２Ｌ，２Ｒをユーザの左右の側頭筋上の皮膚に貼付する。例えば、電極サイズが３０mm×３０mmの布電極２Ｌ，２Ｒを、ユーザの左右の側頭筋上の皮膚に、電極間中心距離が４０mmとなるように位置決めして貼付する。
【００２７】
そして、安静状態において、制御ユニット１０がｓＥＭＧ計測データ取得制御部１１の制御の下、上記各電極２Ｌ，２Ｒにより計測された表面筋電図を表す計測データを、５秒間にわたり入出力インタフェースユニット３０から取り込み、計測データ記憶部２１に記憶させる。なお、このとき入出力インタフェースユニット３０では、サンプリングレート１０００Hz、ゲイン倍率１０００倍で、計測信号ＩＳL ，ＩＳR のサンプリングが行われ、ディジタル信号からなる計測データに変換される。
【００２８】
次に筋電計測装置１は、ＢＰＦ処理部１２により、上記計測データ記憶部２１から各計測データを読み出して、当該計測データに対しバンドパスフィルタ処理を行う。このときバンドパスフィルタ処理の通過帯域は２０−４５０Hzに設定される。この結果、ＢＰＦ処理部１２により、上記通過帯域は２０−４５０Hz以外の帯域に含まれるノイズ成分が除去された計測データが得られる。
【００２９】
続いて、ＲＭＳ算出部１３の制御の下、上記ＢＰＦ処理部１２によりノイズ成分が除去された各計測データに対し１００msecごとに窓を設定し、これらの窓ごとに上記各計測データの二乗平均平方根（ＲＭＳ）を算出する。そして、上記各計測データについて、上記窓ごとに算出されたＲＭＳの上記５秒間分の平均値をそれぞれ算出し、この算出されたＲＭＳの平均値を、各計測データのベースライン値としてベースライン値記憶部２２に保存する。かくして、左右の各側頭筋に対応するベースライン値が得られる。
【００３０】
（２）表面筋電図の計測
上記ベースライン値の決定が終了すると、筋電計測装置１はユーザの例えば咀嚼運動の計測処理を以下のように実行する。図３はその処理手順と処理内容を示すフローチャートである。
【００３１】
筋電計測装置１は、上記ベースライン値の設定処理と同様に、先ずｓＥＭＧ計測データ取得制御部１１の制御の下、ステップＳ１１において、電極２Ｌ，２Ｒにより検出された左右の各側頭筋の表面筋電図を表す計測データを、入出力インタフェースユニット３０から時系列に従い順次取り込み、計測データ記憶部２１に記憶させる。
【００３２】
筋電計測装置１は、次にＢＰＦ処理部１２の制御の下、ステップＳ１２において、上記計測データ記憶部２１から左右の各側頭筋に対応する計測データを読み出し、これらの計測データに対し２０−４５０Hzの帯域のみを通過させて他の帯域を減衰させるフィルタリング処理を行う。このフィルタリング処理により、咀嚼運動とは無関係のノイズ成分の大部分が除去される。
【００３３】
筋電計測装置１は、続いてＲＭＳ算出部１３の制御の下、ステップＳ１３により、上記フィルタリング処理後の左右の各側頭筋に対応する計測データに対し、それぞれ１００msecの時間幅の窓を設定し、これらの窓ごとに上記各計測データの特徴量を表す二乗平均平方根（ＲＭＳ）を算出する。
【００３４】
筋電計測装置１は、次に相関係数算出部１５の制御の下、ステップＳ１４において、上記左右の各側頭筋に対応する計測データのＲＭＳに対し、５秒間の窓を設定して、この窓の期間における上記各計測データのＲＭＳ間の相関係数（Ｒ）を算出する。
【００３５】
筋電計測装置１は、続いて補正処理部１６の制御の下、先ずステップＳ１５において、上記算出された相関係数（Ｒ）をしきい値と比較し、相関係数（Ｒ）がしきい値以下であるか否かを判定する。しきい値は、例えば０．９に設定される。
【００３６】
上記比較の結果、相関係数（Ｒ）がしきい値より大きければ、上記左右の各側頭筋の計測データから算出したＲＭＳの相似性は高く、左右のＲＭＳ値を不均衡にするノイズ成分は含まれていないと判断する。そして、ステップＳ１６ではＲＭＳ値の補正を行わずに、ステップＳ１８に移行し、上記計測されたＲＭＳ値をそのまま計測結果を表すデータＯＳとして、入出力インタフェースユニット３０から図示しないユーザ端末または管理用サーバ等へ出力させる。
【００３７】
これに対し、上記比較の結果相関係数（Ｒ）がしきい値以下だったとする。この場合補正処理部１６は、ＲＭＳにノイズ成分が多く含まれており補正が必要と判断する。そして、ステップＳ１７において、ベースライン値記憶部２２からベースライン値を読み出し、上記ＲＭＳ値を上記ベースライン値に置換する。そして、ステップＳ１８に移行して、上記補正後のＲＭＳ値を計測結果を表すデータＯＳとして、入出力インタフェースユニット３０から図示しないユーザ端末または管理用サーバ等へ出力させる。
【００３８】
（３）計測結果の具体例
（３−１）咀嚼運動
咀嚼運動を行った場合、その閉口運動では左右の側頭筋がほぼ同時に同程度の収縮力で動くため、ＲＭＳの値は左右で同程度の値となり、結果的に相関係数（Ｒ）は高い値、つまり“１”に近い値になる。
図４は、咀嚼運動を５秒間行ったときの左右のＲＭＳの値の変化と、相関係数（Ｒ）の一例を示すものである。同図に示すように、ＲＭＳの値は左右で同程度の値となり、その結果相関係数（Ｒ）は“０．９７”と高い値となる。
【００３９】
（３−２）歯を食いしばった場合
歯を食いしばった場合、左右の側頭筋は連動するため、上記咀嚼運動の場合と同様に、ＲＭＳの値は左右で同程度の値となり、相関係数（Ｒ）も高い値となる。  
図５は、歯を食いしばった状態を５秒間保持したときの左右のＲＭＳの値の変化と、相関係数（Ｒ）の一例を示すものである。同図に示すように、ＲＭＳの値は左右で類似した値となり、その結果相関係数（Ｒ）は“０．９５”と高い値となる。
【００４０】
（３−３）瞬きを行った場合
これに対し瞬きを行った場合には、側頭筋は活動していないにも拘わらず、側頭筋表面部の皮膚が変形する。このため、ＲＭＳにはバンドパスフィルタ処理では除去できないノイズ成分が混入してしまう。この場合、瞬きによる計測値へのノイズ成分の混入の度合いは左右で異なり、一定の関係にならないため、相関係数は低くなる。  
図６は、瞬きを５秒間行ったときの左右のＲＭＳの値の変化と、相関係数（Ｒ）の算出結果の一例を示すものである。同図に示すように、ＲＭＳの値は左右で関連性が少なくなり、その結果相関係数（Ｒ）は“０．８３”と、咀嚼運動等を行った場合に比べて低くなる。
【００４１】
（３−４）帽子の鍔を左右に動かした場合
ユーザが帽子の鍔を持って左右に動かすと、この場合も側頭筋は活動していないにも拘わらず、電極がずれることによりノイズ成分が混入してしまう。この場合、ノイズ成分の混入の度合いは左右で異なり一定とはならないため、相関係数は低くなる。  
図７は、帽子を左右に５秒間動かしたときの左右のＲＭＳの値の変化と、相関係数（Ｒ）の算出結果の一例を示すものである。同図に示すように、ＲＭＳの値は左右で関連性が少なくなり、その結果相関係数（Ｒ）は“０．８”と、咀嚼運動を行った場合等に比べて低くなる。
【００４２】
本実施形態では、以上の計測結果に基づいて、しきい値を“０．９”に設定し、補正の要否を判定している。従って、ユーザの咀嚼運動をノイズの発生要因となる他の動きと区別して、精度良く判定することが可能となる。
【００４３】
（効果）
以上詳述したように一実施形態では、左右の各側頭筋の表面筋電図（ｓＥＭＧ）の計測データを取得し、これらの計測データに対し２０−４５０Hzのみを通過させ他の周波数成分を除去するバンドパスフィルタ処理を行った後、特徴量を示すＲＭＳを算出する。そして、上記左右の側頭筋に対応するＲＭＳ間の相関係数（Ｒ）を算出し、この相関係数（Ｒ）がしきい値以下の場合に、上記ＲＭＳには側頭筋の運動に寄らないノイズ成分が多く残留していると見做して、上記ＲＭＳの値を事前に設定しておいたベースライン値に置換するようにしている。
【００４４】
従って、フィルタ処理技術のみでは除去できないノイズ成分を除去することが可能となり、これにより解析精度のさらなる向上を図ることができる。また、バンドパスフィルタによるノイズ除去処理を併用しているため、例えば咀嚼運動以外の体動成分の影響を事前に排除することができ、これによりＲＭＳの算出および補正の要否判定をより正確に行うことが可能となる。さらに、補正用のベースライン値を事前にベースライン値記憶部２２に記憶しているため、ベースライン値を読み出すだけで簡単な処理により補正処理を行うことができる。しかも、補正用の要否判定に使用するしきい値は、ユーザごとに事前に計測したデータをもとに設定するようにしているため、ユーザごとに適切なしきい値を設定することができる。
【００４５】
［他の実施形態］
前記一実施形態では、補正処理方法として、算出された相関係数がしきい値以下の場合に、算出されたＲＭＳをベースライン値に置換するようにした。しかし、それに限らず、相関係数としきい値との差に応じてベースライン値に対する補正量を算出し、この補正量に基づいて、算出されたＲＭＳの値を補正するようにしてもよい。
【００４６】
前記一実施形態では、ベースライン値を実際の計測結果に基づいて設定するようにしたが、キャリブレーション処理を省略し、事前に用意された標準的なベースライン値を記憶するようにしてもよい。また、しきい値についても、予め固定的な設定した値を使用する以外に、咀嚼運動等の解析精度に基づいて定期的または随時更新するようにしてもよい。
【００４７】
さらに、前記一実施形態では、咀嚼運動を解析するために左右の各側頭筋の動きを計測する場合を例にとって説明した。しかし、それに限るものではなく、側頭筋の他にも奥歯付近に位置する咬筋も咀嚼運動時に左右で連動するため、同様の処理手順でノイズ成分を除去して表面筋電図（ｓＥＭＧ）を測定することができる。
【００４８】
また、咀嚼運動以外にも、屈伸運動等のように左右対称と筋肉が連動して動作する他の運動にも、この発明は適用可能である。その他、筋電計測装置の構成、表面筋電図の計測処理の手順と処理内容、しきい値や等についても、この発明の要旨を逸脱しない範囲で種々変形して実施可能である。
【００４９】
要するにこの発明は、上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合せにより種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態に亘る構成要素を適宜組み合せてもよい。
【符号の説明】
【００５０】
１…筋電計測装置、２Ｒ，２Ｌ…電極、１０…制御ユニット、１１…ｓＥＭＧ計測データ取得制御部、１２…ＢＰＦ処理部、１３…ＲＭＳ算出部、１４…ベースライン値算出部、１５…相関係数算出部、２０…記憶ユニット、２１…計測データ記憶部、２２…ベースライン値記憶部、３０…入出力インタフェースユニット。
【図１】
2019030520.tif 000003
【図２】
2019030520.tif 000004
【図３】
2019030520.tif 000005
【図４】
2019030520.tif 000006
【図５】
2019030520.tif 000007
【図６】
2019030520.tif 000008
【図７】
2019030520.tif 000009