JP5913623B2 - センシング装置及びセンシング方法 - Google Patents

Description

この発明は、アナライトの濃度を連続的又は間欠的に定量するセンシング装置及びセンシング方法に関する。

例えば、被検体の体内にセンサ部を埋め込み、アナライトとしての血中グルコースの濃度を連続的又は間欠的に定量する持続血糖モニタシステム｛ＣＧＭ（Continuous Glucose Monitoring）システムともいう｝が近時開発されている。センサ部を介して計測信号を取得する際、通常、電気ノイズ又は光量ノイズを含む各種ノイズがこの計測信号に混入する。そこで、グルコースの濃度の定量精度を向上させるため、ノイズ成分を効果的に除去するフィルタリングに関する技術が種々提案されている。

特開２００５−１３１３７０号公報では、時間領域でのフィルタリング・アルゴリズム（特に、カルマンフィルタ）を用いて、ノイズ成分を除去する方法が提案されている。より詳細には、誤差共分散行列を信号差パラメータ（例えば、標準偏差）の関数で定義することで、フィルタ係数が動的に最適化される旨が記載されている。

ところで、装置の小型化や消費電力の削減を図るため、電気回路による処理演算量を極力少なくしたいという設計上の要望がある。しかし、特開２００５−１３１３７０号公報に記載の方法では処理演算量が大きくなる傾向があり、上記した要望に応えられないという問題があった。

本発明は上記した問題を解決するためになされたもので、比較的に簡便な構成である周波数領域でのフィルタを用いて、計測信号からノイズ成分を効果的に除去しつつも、アナライトの濃度の時間変化に対する追従性を維持可能なセンシング装置及びセンシング方法を提供することを目的とする。

本発明に係るセンシング装置は、アナライトの濃度を連続的又は間欠的に定量する装置であって、前記アナライトの濃度に相関する計測信号を逐次取得するセンサ部と、複数種類のフィルタを備えており、該複数種類のフィルタのうち１種類のフィルタを介すことで、前記センサ部により取得された前記計測信号の時系列に対して周波数領域でのフィルタ処理を施すフィルタ処理部と、前記フィルタ処理に使用される前記１種類のフィルタを前記計測信号の時間変化量に応じて切り替えるフィルタ切替部とを備える。

このように、周波数領域でのフィルタ処理に使用される１種類のフィルタを計測信号の時間変化量に応じて切り替えるフィルタ切替部を設けたので、前記計測信号の時間変化、及び前記フィルタ処理に起因する位相遅延特性を併せて考慮したフィルタを適時に選択可能である。これにより、比較的に簡便な構成である周波数領域でのフィルタを用いて、計測信号からノイズ成分を効果的に除去しつつも、アナライトの濃度の時間変化に対する追従性を維持できる。

また、前記フィルタ処理部は、前記計測信号の時系列に対して恒等変換を施す恒等変換フィルタを少なくとも備えており、前記フィルタ切替部は、前記時間変化量が閾値よりも大きい場合に前記恒等変換フィルタに切り替えることが好ましい。時間変化量が閾値よりも大きい場合にフィルタ処理による位相遅延を発生させないようにしたので、アナライトの濃度の時間変化に対する追従性を確保できる。

更に、前記フィルタ処理部は、遮断周波数以下の帯域での位相遅延量の平均値が異なる２種類のフィルタを少なくとも備えており、前記フィルタ切替部は、前記時間変化量が大きい場合に前記位相遅延量の平均値が小さいフィルタに切り替えると共に、前記時間変化量が小さい場合に前記位相遅延量の平均値が大きいフィルタに切り替えることが好ましい。時間変化量が大きい場合にフィルタ処理による位相遅延を発生させないようにしたので、アナライトの濃度の時間変化に対する追従性を確保できる。また、時間変化量が小さい場合には、上記した追従性がそれほど要求されないので、計測信号からノイズ成分を一層効果的に除去できる。

本発明に係るセンシング方法は、アナライトの濃度を連続的又は間欠的に定量する方法であって、センサを用いて、前記アナライトの濃度に相関する計測信号を逐次取得する取得ステップと、複数種類のフィルタのうち１種類のフィルタを介すことで、前記センサにより取得された前記計測信号の時系列に対して周波数領域でのフィルタ処理を施す処理ステップと、前記フィルタ処理に使用される前記１種類のフィルタを前記計測信号の時間変化量に応じて切り替える切替ステップとを備える。

また、前記複数種類のフィルタには、前記計測信号の時系列に対して恒等変換を施す恒等変換フィルタが少なくとも含まれており、前記切替ステップでは、前記時間変化量が閾値よりも大きい場合に前記恒等変換フィルタに切り替えることが好ましい。

更に、前記複数種類のフィルタには、遮断周波数以下の帯域での位相遅延量の平均値が異なる２種類のフィルタが少なくとも含まれており、前記切替ステップでは、前記時間変化量が大きい場合に前記位相遅延量の平均値が小さいフィルタに切り替えると共に、前記時間変化量が小さい場合に前記位相遅延量の平均値が大きいフィルタに切り替えることが好ましい。

本発明に係るセンシング装置及びセンシング方法によれば、周波数領域でのフィルタ処理に使用される１種類のフィルタを計測信号の時間変化量に応じて切り替えるようにしたので、前記計測信号の時間変化、及び前記フィルタ処理に起因する位相遅延特性を併せて考慮したフィルタを適時に選択可能である。これにより、比較的に簡便な構成である周波数領域でのフィルタを用いて、計測信号からノイズ成分を効果的に除去しつつも、アナライトの濃度の時間変化に対する追従性を維持できる。

第１及び第２実施形態に共通するセンシング装置の概略ブロック図である。 第１実施形態に係るセンサ制御回路のブロック図である。 図２に示す第１フィルタの回路構成図である。 図４Ａは、図２に示す第１フィルタのフィルタ係数を示す図である。図４Ｂは、図４Ａのフィルタ係数に応じたフィルタ特性を示すグラフである。 図５Ａ及び図５Ｂは、血中グルコースの濃度変化及びその定量結果を表すグラフである。 第１実施形態に係るセンシング装置の動作説明に供されるフローチャートである。 複数の標本点から時間変化量を算出する一例を表す概略説明図である。 ３通りのフィルタ処理を施して定量した場合における、グルコースの定量誤差率の累積ヒストグラムである。 第２実施形態に係るセンサ制御回路のブロック図である。 図１０Ａは、図９に示す第２フィルタのフィルタ係数を示す図である。図１０Ｂは、図１０Ａのフィルタ係数に応じたフィルタ特性を示すグラフである。 第２実施形態に係るセンシング装置の動作説明に供されるフローチャートである。 ４通りのフィルタ処理を施して定量した場合における、グルコースの定量誤差率の累積ヒストグラムである。

以下、本発明に係るセンシング方法について、センシング装置との関係において好適な実施形態を挙げ、添付の図面を参照しながら説明する。

［第１及び第２実施形態に共通するセンシング装置１０の構成］
先ず、第１及び第２実施形態に共通するセンシング装置１０の構成について、図１の概略ブロック図を参照しながら説明する。

図１に示すように、センシング装置１０は、センサ部（センサ）１２と、センサ制御回路１４（センサ制御回路８０）と、演算部１６と、電源回路１８と、ＲＯＭ（Read Only Memory）２０と、ＲＡＭ（Random Access Memory）２２と、クロック発生器２４と、入力部２６と、表示器２８とを基本的に備える。

センサ部１２は、アナライトの濃度に相関する信号（以下、計測信号Ｓ）を取得する。センサ部１２として、サンプリング間隔Ｔｓを容易に変更可能な光学センサ（例えば、蛍光センサ）を適用することが好ましい。なお、センサ部１２の形態はこれに限定されるものではなく、例えば、血糖値をグルコースオキシダーゼ（ＧＯＤ）等の酵素を用いた酵素電極法等による電気的（電気化学方式）に測定するセンサを適用してもよい。

センサ制御回路１４は、センサ部１２を駆動制御することで、所望のタイミングで計測信号Ｓを取得可能である。センサ制御回路１４は、計測信号Ｓとしての電流値（アナログ信号）を電圧値に変換すると共に、この電圧値を量子化してデジタル信号に変換する。センサ制御回路１４は、このアナログ信号又はデジタル信号に対して所定のフィルタ処理を施すことで、計測信号Ｓに混入するノイズ成分を除去する。

演算部１６は、ＣＰＵ（Central Processing Unit）、ＭＰＵ（Micro-Processing Unit）等で構成されており、ＲＯＭ２０に予め記憶されたプログラムを読み出し、後述する各種信号処理を実行する。演算部１６は、センサ制御回路１４から取得した信号値Ｓｆ（ｋ）に基づいてアナライトの濃度を定量する濃度定量部２９として機能する。

電源回路１８は、演算部１６を含むセンシング装置１０内の各構成要素に電力を供給する。ＲＡＭ２２は、センサ部１２を介して入力された計測信号Ｓの他、本発明に係るセンシング方法を実施するために必要な各種データを読み出し又は書込み可能である。クロック発生器２４は、所定周期でクロック信号を発生し、演算部１６側に供給する。これにより、演算部１６は、信号値Ｓｆ（ｋ）の取得タイミングを制御可能である。

入力部２６は、演算部１６での演算に供される各種情報（例えば、定量間隔Ｔｄ）を入力可能に設けられている。例えば、押圧式ボタンであってもよいし、表示器２８に組み込まれたタッチパネルであってもよい。表示器２８は、演算部１６により定量されたアナライトの濃度（以下、定量濃度ともいう）に関する各種情報を可視化して表示する。表示器２８は、モノクロ又はカラー表示可能な表示モジュールであり、液晶パネル、有機ＥＬ（Electro-Luminescence）、無機ＥＬパネル等で構成されてもよい。

なお、センサ部１２は、酵素センサ、グルコースセンサ、ｐＨセンサ、免疫センサ、又は微生物センサ等、多様な用途に適用可能である。また、センサ部１２の構成は、本構成に限られることなく種々の構成を採り得る。例えば、物理的に分離されたセンサ制御回路１４（８０）及び演算部１６の間を無線で通信可能に設けることで、センサ部１２を被検体の体内に完全に埋め込んだ状態で間欠的又は連続的に定量可能である。なお、無線通信の際には、近距離通信用の規格（例えば、「IEEE 802.15.6」で規定するボディエリアネットワーク等）を適用してもよい。

＜第１実施形態＞
続いて、第１実施形態に係るセンサ制御回路１４の構成及び動作について、図２〜図７を参照しながら説明する。なお、本明細書では、アナライトとしてグルコースを用いた場合の定量動作を中心に説明する。

［センサ制御回路１４のブロック図］
図２は、第１実施形態に係るセンサ制御回路１４（図１参照）のブロック図である。

センサ制御回路１４は、センサ部１２からの計測信号Ｓを入力する信号入力部３０と、アナログ信号である計測信号Ｓをデジタル信号である原信号値Ｓ（ｋ）に変換するアナログ・デジタル変換部（以下、ＡＤＣ３２という）と、複数種類のフィルタの中から１種類のフィルタを択一的に切り替えるフィルタ切替部３４と、原信号値Ｓ（ｋ）に対して周波数領域でのフィルタ処理を施すフィルタ処理部３６と、直近の原信号値Ｓ（ｋ）を一時的に格納するバッファメモリ４０と、フィルタ切替部３４が備える複数種類のフィルタを切り替える変数（以下、切替変数Ｖｓ）を算出する切替変数算出部４２と、を備える。

フィルタ切替部３４のスイッチ４４は、第１端子４６ａ又は第２端子４６ｂのいずれか一方に接続された状態（オン状態）、或いは、第１端子４６ａ及び第２端子４６ｂのいずれにも接続されていない状態（オフ状態）を採り得る。

フィルタ処理部３６は、周波数領域でのフィルタリング処理を実行するための第１フィルタ４８を備える。第１フィルタ４８は、フィルタ切替部３４の第１端子４６ａ側に接続されている。一方、フィルタ切替部３４の第２端子４６ｂ側には、フィルタが接続されていない。換言すれば、フィルタ処理部３６は、原信号値Ｓ（ｋ）に対して恒等変換を施す恒等変換フィルタ４９を更に備えている。

図３は、第１フィルタ４８の回路構成図である。第１フィルタ４８は、５つの乗算器５１、５２、５３、５４、５５、４つの加算器５６、５７、５８、５９、並びに４つの遅延器６０、６１、６２、６３で構成される。すなわち、第１フィルタ４８は、タップ数が５であるＦＩＲ（Finite Impulse Response）フィルタに相当する。以下、乗算器５１〜５５に設定された乗数（以下、フィルタ係数）をそれぞれ順番に、ｈ０、ｈ１、ｈ２、ｈ３、ｈ４と表記する。なお、ＦＩＲフィルタは周知の電気回路であるため、各演算器の機能及び接続関係についての説明を省略する。

図４Ａは、図２に示す第１フィルタ４８のフィルタ係数を示す図である。具体的には、乗算器５１（以下、図３参照）のフィルタ係数をｈ０＝０．１５９、乗算器５２のフィルタ係数をｈ１＝０．２２０、乗算器５３のフィルタ係数をｈ２＝０．２４３、乗算器５４のフィルタ係数をｈ３＝０．２２０、乗算器５５のフィルタ係数をｈ４＝０．１５９、にそれぞれ設定する。

図４Ｂは、図４Ａのフィルタ係数に応じたフィルタ特性を示すグラフである。

実線で示すグラフの横軸は周波数（単位：ｍＨｚ）であり、縦軸は振幅比（単位：なし）である。ここで、振幅比とは、入力される周期信号（正弦波形信号）の振幅に対する出力信号の振幅の比である。理想的には、信号成分を通過させる周波数帯域で１であり、信号成分を遮断させる周波数帯域で０である。本図例のフィルタ特性は、０〜ｆｃ［ｍＨｚ］の帯域では５０％以上の信号成分を通過させると共に、ｆｃ［ｍＨｚ］以上の帯域では５０％以上の信号成分を遮断させる、いわゆるローパスフィルタ型の特徴を示している。以下、ｆｃ＝０．４４［ｍＨｚ］のことを遮断周波数（カットオフ周波数）と称する。

破線で示すグラフの横軸は周波数（単位：ｍＨｚ）であり、縦軸は位相遅延量（単位：ｓｅｃ）である。ここで、位相遅延量とは、入力される周期信号（正弦波形信号）の位相に対する出力信号の位相の差であり、理想的には０である。本図例のフィルタ特性は、０〜０．９［ｍＨｚ］の帯域では相対的に大きな位相の遅延が生じると共に、０．９［ｍＨｚ］以上の帯域では位相遅延量が相対的に小さくなる特徴を示している。なお、遮断周波数ｆｃ以下の帯域Ｒ１での位相遅延量の平均値は、約６００［ｓｅｃ］である。

なお、第１フィルタ４８（又は後述する第２フィルタ８４）を決定する際、ＦＩＲフィルタ又はＩＩＲ（Infinite Impulse Response）フィルタに関する公知の設計手法を種々適用してもよい。例えば、信号の通過域に関して、低域通過フィルタ、高域通過フィルタ、帯域通過フィルタ、帯域除去フィルタ、及び全域通過フィルタのうちのいずれかを適用してもよい。また、振幅特性の形状に関して、バターワース特性、チェビシェフ特性、逆チェビシェフ特性、及び連立チェビシェフ特性（楕円特性）のうちのいずれかを適用してもよい。

［フィルタ処理と定量結果との関係］
図５Ａ及び図５Ｂは、血中グルコースの濃度変化及びその定量結果を表すグラフである。各グラフの横軸は時間（単位：ｍｉｎ）であり、縦軸はグルコースの濃度、いわゆる血糖値（単位：ｍｇ／ｄｌ）である。実線のグラフでそれぞれ示す血糖値の変化は、被検体の食事前後における体内の血糖値の時間変化を模擬している。

図５Ａにおける破線のグラフは、計測信号Ｓに対してフィルタ処理を施すことなく得られた定量結果を示す。本グラフから理解されるように、計測信号Ｓに高周波数のノイズ成分が混入することで、実際値と定量値との間に不規則な誤差が生じている。特に、信号のレベルが相対的に低いＢ領域において、ノイズ成分に起因する定量精度への影響が顕著に現われている。

図５Ｂにおける破線のグラフは、計測信号Ｓに対して第１フィルタ４８によるフィルタ処理を施して得られた定量結果を示す。本グラフから理解されるように、ノイズ成分に起因する不規則な誤差が少ない定量値が得られる。しかし、第１フィルタ４８による位相遅延が起こることで、実際値と定量値との間に乖離が生じている。特に、信号の時間変化量が大きいＡ領域において、追従性の低下に起因する定量精度への影響が顕著に現われている。

［センサ制御回路１４を含むセンシング装置１０の動作］
上記した定量誤差の発生を抑制するため、第１実施形態に係るセンシング方法ではフィルタの有無を適時に切り替える。以下、センサ制御回路１４（図２）を含むセンシング装置１０の動作について、図６のフローチャートを参照しながら詳細に説明する。初期状態では、フィルタ切替部３４のスイッチ４４はオフ状態であったとする。

ステップＳ１において、信号入力部３０は、所定のサンプリング間隔Ｔｓに従って、センサ部１２から計測信号Ｓを入力する。その後、ＡＤＣ３２は、信号入力部３０から取得したアナログ信号をデジタル信号｛以下、原信号値Ｓ（ｋ）という｝に変換する。

ステップＳ２において、ステップＳ１で入力・取得された原信号値Ｓ（ｋ）をバッファメモリ４０に一時的に格納する。

ステップＳ３において、センサ制御回路１４は、アナライトの濃度を定量する指示があったか否かを判別する。具体的には、センサ制御回路１４は、濃度の定量を指示する信号（以下、定量指示信号）を演算部１６から受信したか否かを判別する。

演算部１６は、ステップＳ１及びＳ２の実行と並列して、クロック発生器２４から入力されたクロック信号のパルス数をカウントする。そして、カウント上限値（定量間隔Ｔｄに相当する値）に到達した場合、演算部１６は、センサ制御回路１４側に定量指示信号を送出し、次のステップ（Ｓ４）に進む。

一方、上記したカウント上限値に到達していない場合、演算部１６は、定量指示信号を送出することなく、パルス数のカウントを継続する。すなわち、ステップＳ１に戻り、以下ステップＳ１及びＳ２を順次繰り返す。

ところで、サンプリング間隔Ｔｓは、センサ制御回路１４側に計測信号Ｓを入力する時間の間隔であると共に、定量間隔Ｔｄは、演算部１６側でアナライトの濃度を定量する時間間隔である。すなわち、定量間隔Ｔｄは、サンプリング間隔Ｔｓと異なるパラメータであることから、サンプリング間隔Ｔｓと同一の又は異なる値を取ってもよい。また、定量間隔Ｔｄがサンプリング間隔Ｔｓに等しい場合、センサ制御回路１４は、定量指示信号を受け付けた後に原信号値Ｓ（ｋ）を取得・格納してもよい。

ステップＳ４において、切替変数算出部４２は、ステップＳ２で順次格納された原信号値Ｓ（ｋ）の時系列に基づいて、計測信号Ｓの時間変化量を表すパラメータである切替変数Ｖｓを算出する。ここで、時間変化量とは、直近の複数の標本点から推定される計測信号Ｓの変動傾向（トレンド）を意味する。

図７に示すように、現在の標本点７５から時間的に近い順に、標本点７４、７３、７２、７１が既に得られていたとする。この場合、過去の標本点７１〜７４のみならず標本点７５も併せ用いて、破線で図示する回帰直線７６を求めることができる。そして、切替変数算出部４２は、この回帰直線７６の勾配（時間の１次微分）の絶対値を切替変数Ｖｓとして算出する。

なお、上記勾配を算出する手法は、加重平均、最小二乗法を含む種々の最適化手法を用いることができる。また、トレンドの推定に供される標本点の個数は５つに限られず、演算量、処理時間等を総合的に考慮して適宜決定してもよい。更に、切替変数算出部４２は、回帰直線７６の勾配に限られず、例えば、隣接する標本点を結ぶ直線勾配の統計値（例えば、平均値）、近似曲線における曲率（時間の２次微分）等を用いて、切替変数Ｖｓを算出してもよい。

ステップＳ５において、フィルタ切替部３４は、ステップＳ４で算出された切替変数Ｖｓに応じて、フィルタ切替部３４のスイッチ４４を切り替える。具体的には、フィルタ切替部３４は、予め設定された閾値Ｖｓ^＊との大小関係によってスイッチ４４の切替状態を決定する。

Ｖｓ≦Ｖｓ^＊を満たす場合、フィルタ切替部３４は、スイッチ４４を一定時間だけ第１端子４６ａ側に切り替える（ステップＳ６）。そうすると、原信号値Ｓ（ｋ）は、スイッチ４４、第１端子４６ａ及び第１フィルタ４８を介して、センサ制御回路１４の外部に出力される。

Ｖｓ＞Ｖｓ^＊を満たす場合、フィルタ切替部３４は、スイッチ４４を一定時間だけ第２端子４６ｂ側に切り替える（ステップＳ７）。そうすると、原信号値Ｓ（ｋ）は、スイッチ４４及び第２端子４６ｂを介して、センサ制御回路１４の外部に出力される。以下、用語の区別を明確にするため、フィルタ処理部３６を通過し、センサ制御回路１４から出力された原信号値Ｓ（ｋ）を「信号値Ｓｆ（ｋ）」と称する。

なお、フィルタ切替部３４は、閾値Ｖｓ^＊による判別処理の際にデッドバンド（不感帯）を設けてもよい。これにより、時系列的な判別結果のゆらぎを抑制可能であり、フィルタの切替制御を安定的に実行できる。

ステップＳ８において、濃度定量部２９は、ＲＡＭ２２から読み出した定量係数等を用いて、センサ制御回路１４から出力された信号値Ｓｆ（ｋ）に基づく濃度の定量を行う。ここで、濃度の定量方法は、センサ部１２における検出方式、材質、感度特性、個体差等に適した種々の手法を採ることができる。

ステップＳ９において、表示器２８は、ステップＳ８における定量結果を表示する。表示処理に先立ち、演算部１６は、得られた定量結果のうち、表示器２８に表示させる可視情報（以下、定量可視情報という。）を決定した後、その定量可視情報に応じた制御信号を表示器２８側に供給する。なお、定量可視情報として、定量値のみならず、例えば、トレンド、定量の成否、定量時刻、診断結果等が挙げられる。

ステップＳ１０において、演算部１６は、この一連の定量動作の終了指示があったか否かを判別する。終了指示がなかったと判別された場合、ステップＳ１に戻り、以下ステップＳ１〜Ｓ９の動作を同様に繰り返す。一方、終了指示があった場合、センシング装置１０は、アナライトの定量動作を終了する。このようにして、演算部１６は、所定の定量間隔Ｔｄに従って、各定量時点における濃度の時系列データを得る。

［第１実施形態に係るセンシング方法により得られる作用効果］
以下、第１実施形態に係るセンシング方法により得られる作用効果について、図８を参照しながら説明する。より詳細には、図５Ａ及び図５Ｂに示す血中グルコースの濃度変化（実線のグラフ）の下、異なる形態のフィルタ処理を用いてそれぞれ定量した結果を比較する。いずれの場合も、サンプリング間隔Ｔｓ＝５［ｍｉｎ］、定量間隔Ｔｄ＝５［ｍｉｎ］として計測・定量した。また、回帰直線７６（図７）の勾配の絶対値に関する閾値を、Ｖｓ^＊＝０．３［１／ｍｉｎ］に設定した。

図８は、３通りのフィルタ処理を施して定量した場合における、グルコースの定量誤差率の累積ヒストグラムである。ヒストグラムの横軸は定量値の誤差率（単位：％）であり、縦軸は累積頻度（単位：％）である。

なお、本図中における「フィルタなし」は、スイッチ４４が第２端子４６ｂに常時接続された場合での定量結果（図５Ａに示す破線のグラフ）に対応する。「フィルタ固定」は、スイッチ４４が第１端子４６ａに常時接続された場合での定量結果（図５Ｂに示す破線のグラフ）に対応する。「フィルタ切替」は、図６のフローチャートに従ってスイッチ４４を適時切り替えた場合での定量結果に対応する。

本図から理解されるように、誤差率が１０％以下の頻度は、「フィルタ切替」＞「フィルタなし」＞「フィルタ固定」の順番で高くなっている。そして、誤差率が２０％以下の頻度は、「フィルタ切替」＞「フィルタ固定」＞「フィルタなし」の順番で高くなっている。このように、「フィルタ切替」は、「フィルタなし」及び「フィルタ固定」と比べて、濃度の定量誤差が有意に少ないと結論付けられる。

＜第２実施形態＞
続いて、第２実施形態に係るセンサ制御回路８０の構成及び動作について、図９〜図１２を参照しながら説明する。なお、第１実施形態と同様の構成については、同一の参照符号を付すると共にその説明を省略する。

［センサ制御回路８０のブロック図］
図９は、第２実施形態に係るセンサ制御回路８０（図１参照）のブロック図である。センサ制御回路８０は、センサ制御回路１４（図２）と略同様の構成を採るが、フィルタ処理部３６に代わって別の構成のフィルタ処理部８２を備える。

フィルタ処理部８２は、第１フィルタ４８と、該第１フィルタ４８と同じ回路構成（図３参照）である第２フィルタ８４とを備える。第２フィルタ８４は、フィルタ切替部３４の第２端子４６ｂ側に接続されている。

図１０Ａは、図９に示す第２フィルタ８４のフィルタ係数を示す図である。具体的には、乗算器５１（以下、図３参照）のフィルタ係数をｈ０＝０．３０１、乗算器５２のフィルタ係数をｈ１＝０．３９８、乗算器５３のフィルタ係数をｈ２＝０．３０１、乗算器５４のフィルタ係数をｈ３＝０．０００、乗算器５５のフィルタ係数をｈ４＝０．０００、にそれぞれ設定する。このように、ｈ３＝ｈ４＝０．０００であることから、第２フィルタ８４は、実質的にタップ数が３であるＦＩＲフィルタとして機能する。

図１０Ｂは、図１０Ａのフィルタ係数に応じたフィルタ特性を示すグラフである。実線で示すグラフの横軸は周波数（単位：ｍＨｚ）であり、縦軸は振幅比（単位：なし）である。本図例のフィルタ特性は、図４Ｂと同様に、ローパスフィルタ型の特徴を示しており、遮断周波数ｆｃは、ｆｃ＝０．７４［ｍＨｚ］である。つまり、第２フィルタ８４の遮断周波数ｆｃは、第１フィルタ４８の遮断周波数ｆｃ（＝０．４４）よりも高い。

破線で示すグラフの横軸は周波数（単位：ｍＨｚ）であり、縦軸は位相遅延量（単位：ｓｅｃ）である。本図例では、０〜１．７［ｍＨｚ］の帯域で位相遅延は略一定（＝３００［ｓｅｃ］）である。この場合、遮断周波数ｆｃ以下の帯域Ｒ２での位相遅延量の平均値は約３００［ｓｅｃ］であり、第１フィルタ４８での位相遅延量の平均値（図４Ｂ参照；約６００［ｓｅｃ］）よりも小さい。

［センサ制御回路８０を含むセンシング装置１０の動作］
上記した欠点を相互に補完するため、第２実施形態に係るセンシング方法では複数種類のフィルタを適時に切り替える。以下、センサ制御回路８０（図９）を含むセンシング装置１０の動作について、図１１のフローチャートを参照しながら説明する。ステップＳ１〜Ｓ４、ステップＳ８〜Ｓ１０に関して、図６のフローチャート（第１実施形態）と同様であるため、その説明を割愛する。

ステップＳ５において、フィルタ切替部３４は、算出された切替変数Ｖｓに応じてスイッチ４４を切り替える。Ｖｓ≦Ｖｓ^＊を満たす場合、フィルタ切替部３４は、スイッチ４４を一定時間だけ第１端子４６ａ側に切り替える（ステップＳ６）。

一方、Ｖｓ＞Ｖｓ^＊を満たす場合、フィルタ切替部３４は、スイッチ４４を一定時間だけ第２端子４６ｂ側に切り替える（ステップＳ７Ａ）。そうすると、原信号値Ｓ（ｋ）は、スイッチ４４、第２端子４６ｂ及び第２フィルタ８４を介して、センサ制御回路８０の外部に出力される。

このようにして、演算部１６は、所定の定量間隔Ｔｄに従って、各定量時点における濃度の時系列データを得る。

［第２実施形態に係るセンシング方法での定量結果］
以下、第２実施形態に係るセンシング方法により得られる作用効果について、図１２を参照しながら説明する。より詳細には、図５Ａ及び図５Ｂに示す血中グルコースの濃度変化（実線のグラフ）の下、異なる形態のフィルタ処理を用いてそれぞれ定量した結果を比較する。いずれの場合も、サンプリング間隔Ｔｓ＝５［ｍｉｎ］、定量間隔Ｔｄ＝５［ｍｉｎ］として計測・定量した。また、閾値Ｖｓ^＊＝０．３［１／ｍｉｎ］に設定した。

図１２は、４通りのフィルタ処理を施して定量した場合における、グルコースの定量誤差率の累積ヒストグラムである。ヒストグラムの横軸は定量値の誤差率（単位：％）であり、縦軸は累積頻度（単位：％）である。

なお、本図中における「フィルタなし」及び「第１フィルタ固定」は、図８の「フィルタなし」及び「フィルタ固定」にそれぞれ対応する。「第２フィルタ固定」は、スイッチ４４が第２端子４６ｂに常時接続された場合での定量結果に対応する。「フィルタ切替」は、図１１のフローチャートに従ってスイッチ４４を適時切り替えた場合での定量結果に対応する。

本図から理解されるように、誤差率が１０％以下の頻度は、「フィルタ切替」＞「第２フィルタ固定」＞「フィルタなし」＞「第１フィルタ固定」の順番で高くなっている。そして、誤差率が２０％以下の頻度は、「フィルタ切替」＞「第２フィルタ固定」＞「第１フィルタ固定」＞「フィルタなし」の順番で高くなっている。このように、「フィルタ切替」は、「フィルタなし」、「第１フィルタ固定」及び「第２フィルタ固定」と比べて、濃度の定量誤差が相対的に少ないと結論付けられる。

［本発明の効果］
以上のように、このセンシング装置１０は、アナライトの濃度に相関する計測信号Ｓを逐次取得するセンサ部１２と、複数種類のフィルタ（４８、４９、８４）のうち１種類のフィルタを介すことで、計測信号Ｓの時系列に対して周波数領域でのフィルタ処理を施すフィルタ処理部３６、８２を備える。

そして、周波数領域でのフィルタ処理に使用される１種類のフィルタを計測信号Ｓの時間変化量（例えば、切替変数Ｖｓ）に応じて切り替えるフィルタ切替部３４を設けたので、計測信号の時間変化、及びフィルタ処理に起因する位相遅延特性を併せて考慮した第１フィルタ４８等を適時に選択可能である。これにより、比較的に簡便な構成である周波数領域でのフィルタ（第１フィルタ４８等）を用いて、計測信号Ｓからノイズ成分を効果的に除去しつつも、アナライトの濃度の時間変化に対する追従性を維持できる。

また、フィルタ処理部３６は、計測信号Ｓの時系列に対して恒等変換を施す恒等変換フィルタ４９を少なくとも備えており、フィルタ切替部３４は、切替変数Ｖｓが閾値Ｖｓ^＊よりも大きい場合に恒等変換フィルタ４９に切り替えてもよい。

更に、フィルタ処理部８２は、遮断周波数ｆｃ以下の帯域での位相遅延量の平均値が異なる２種類のフィルタ（４８、８４）を少なくとも備えており、フィルタ切替部３４は、切替変数Ｖｓが大きい場合に位相遅延量の平均値が小さい第２フィルタ８４に切り替えると共に、切替変数Ｖｓが小さい場合に位相遅延量の平均値が大きい第１フィルタ４８に切り替えてもよい。

切替変数Ｖｓが大きい場合にフィルタ処理による位相遅延を発生させないようにしたので、アナライトの濃度の時間変化に対する追従性を確保できる。また、切替変数Ｖｓが小さい場合には、上記した追従性がそれほど要求されないので、計測信号Ｓからノイズ成分を一層効果的に除去できる。

なお、この発明は、上述した実施形態に限定されるものではなく、この発明の主旨を逸脱しない範囲で自由に変更できることは勿論である。

例えば、第１及び第２実施形態では、フィルタ処理部３６、８２をデジタルフィルタ回路で構成する例を示したが、アナログフィルタ回路で構成してもよい。また、デジタルフィルタを適用する場合、ハードウェア及び／又はソフトウェアで実現してもよい。更に、このフィルタ処理をソフトウェアで実現する場合、センサ制御回路１４、８０に代わって演算部１６に実行させてもよい。

Claims (6)

1. アナライトの濃度を連続的又は間欠的に定量するセンシング装置（１０）であって、
   前記アナライトの濃度に相関する計測信号を逐次取得するセンサ部（１２）と、
   複数種類のフィルタ（４８、４９、８４）を備えており、該複数種類のフィルタ（４８、４９、８４）のうち１種類のフィルタを介すことで、前記センサ部（１２）により取得された前記計測信号の時系列に対して周波数領域でのフィルタ処理を施すフィルタ処理部（３６、８２）と、
   前記フィルタ処理に使用される前記１種類のフィルタを前記計測信号の時間変化量に応じて切り替えるフィルタ切替部（３４）と
   を備えることを特徴とするセンシング装置（１０）。
2. 請求項１記載のセンシング装置（１０）において、
   前記フィルタ処理部（３６）は、前記計測信号の時系列に対して恒等変換を施す恒等変換フィルタ（４９）を少なくとも備えており、
   前記フィルタ切替部（３４）は、前記時間変化量が閾値よりも大きい場合に前記恒等変換フィルタ（４９）に切り替える
   ことを特徴とするセンシング装置（１０）。
3. 請求項１記載のセンシング装置（１０）において、
   前記フィルタ処理部（８２）は、遮断周波数以下の帯域での位相遅延量の平均値が異なる２種類のフィルタ（４８、８４）を少なくとも備えており、
   前記フィルタ切替部（３４）は、前記時間変化量が大きい場合に前記位相遅延量の平均値が小さいフィルタ（８４）に切り替えると共に、前記時間変化量が小さい場合に前記位相遅延量の平均値が大きいフィルタ（４８）に切り替える
   ことを特徴とするセンシング装置（１０）。
4. アナライトの濃度を連続的又は間欠的に定量するセンシング方法であって、
   センサ（１２）を用いて、前記アナライトの濃度に相関する計測信号を逐次取得する取得ステップと、
   複数種類のフィルタ（４８、４９、８４）のうち１種類のフィルタを介すことで、前記センサ（１２）により取得された前記計測信号の時系列に対して周波数領域でのフィルタ処理を施す処理ステップと、
   前記フィルタ処理に使用される前記１種類のフィルタを前記計測信号の時間変化量に応じて切り替える切替ステップと
   を備えることを特徴とするセンシング方法。
5. 請求項４記載のセンシング方法において、
   前記複数種類のフィルタ（４８、４９、８４）には、前記計測信号の時系列に対して恒等変換を施す恒等変換フィルタ（４９）が少なくとも含まれており、
   前記切替ステップでは、前記時間変化量が閾値よりも大きい場合に前記恒等変換フィルタ（４９）に切り替える
   ことを特徴とするセンシング方法。
6. 請求項４記載のセンシング方法において、
   前記複数種類のフィルタ（４８、４９、８４）には、遮断周波数以下の帯域での位相遅延量の平均値が異なる２種類のフィルタ（４８、８４）が少なくとも含まれており、
   前記切替ステップでは、前記時間変化量が大きい場合に前記位相遅延量の平均値が小さいフィルタ（８４）に切り替えると共に、前記時間変化量が小さい場合に前記位相遅延量の平均値が大きいフィルタ（４８）に切り替える
   ことを特徴とするセンシング方法。

Images