(57)【要約】【課題】脈拍検出精度が向上し、かつ、発光素子及び受光素子の数を低減可能な脈拍計測装置及びその制御方法を提供する。【解決手段】脈拍計測装置１は、所定の時間間隔で光を照射する発光素子２１と、その光を生体を介して受光する共に、体動によって形成された隙間を通って外部から入射した外光を受光可能な受光素子２２と、受光素子２２が出力する受光信号から、上記時間間隔に対応する発光周波数を中心とする第１周波数範囲の第１信号成分を抽出する帯域通過フィルタ部５２と、脈動の周波数より低い第２周波数範囲の第２信号成分を抽出する低域通過フィルタ部５３と、第１信号成分に周波数解析を施して複数の周波数スペクトルを特定すると共に第２信号成分に周波数解析を施し、第１信号成分に含まれるが、第２信号成分には含まれない周波数スペクトルｆ１の周波数に基づいて脈拍数を特定する演算処理部６０とを備える。【選択図】図３
2008220722.tif 000002
【特許請求の範囲】
【請求項１】
生体の表面に装着される生体情報検出部を備えた脈拍計測装置において、
前記生体情報検出部に設けられ、前記生体の血流部位に向けて所定の時間間隔で光を照射する単一の発光素子と、
前記生体情報検出部に設けられ、前記光を前記生体を介して受光すると共に、前記生体の体動によって前記生体と前記生体情報検出部との間に形成された所定の隙間を通って外部から入射した外光を受光可能に配置されて、前記生体の脈動及び体動に応じて受光量が変動する受光信号を生成する単一の受光素子と、
前記受光信号から、前記時間間隔に対応する前記発光素子の発光周波数を含む第１周波数範囲の第１信号成分を抽出する第１フィルタと、
前記受光信号から、前記脈動の周波数より低い第２周波数範囲の第２信号成分を抽出する第２フィルタと、
前記第１信号成分に周波数解析を施して複数の周波数スペクトルを特定すると共に、前記第２信号成分に周波数解析を施して一又は複数の周波数スペクトルを特定する周波数解析部と、
前記第１信号成分に含まれるが、前記第２信号成分には含まれない周波数スペクトルを脈動スペクトルと特定し、この脈動スペクトルの周波数に基づいて脈拍数を算出する脈拍数算出部と
を備えることを特徴とする脈拍計測装置。
【請求項２】
請求項１に記載の脈拍計測装置において、
前記第１フィルタは、前記発光周波数を略中心周波数としたバンドパスフィルタであることを特徴とする脈拍計測装置。
【請求項３】
請求項１又は２に記載の脈拍計測装置において、
前記第２フィルタは、前記生体の体動周波数範囲を通過させるローパスフィルタであることを特徴とする脈拍計測装置。
【請求項４】
請求項１乃至３のいずれか一項に記載の脈拍計測装置において、
前記生体の体動によって前記受光素子への外光入射光量が変化するように外光を導入する外光導入部を備えることを特徴とする脈拍計測装置。
【請求項５】
請求項１乃至４のいずれか一項に記載の脈拍計測装置において、
前記発光周波数を１ｋＨｚを超える周波数にしたことを特徴とする脈拍計測装置。
【請求項６】
請求項１乃至５のいずれか一項に記載の脈拍計測装置において、
前記第１フィルタが抽出した信号成分の包絡線を検波する検波部を有し、前記周波数解析部は、前記検波部の出力信号に周波数解析を施して、脈動成分及び体動成分に対応する周波数スペクトルを特定することを特徴とする脈拍計測装置。
【請求項７】
請求項６に記載の脈拍計測装置において、
前記検波部は、前記第１フィルタが抽出した信号間の極大値同士或いは極小値同士をピークホールドし、それらを結んで前記包絡線を形成することを特徴とする脈拍計測装置。
【請求項８】
請求項１乃至７のいずれか一項に記載の脈拍計測装置において、
前記周波数解析部に入力する信号を、前記第１フィルタが抽出した第１信号成分と、前記第２フィルタが抽出した第２信号成分とに選択的に切り換える入力選択部を有することを特徴とする脈拍計測装置。
【請求項９】
請求項１乃至８のいずれか一項に記載の脈拍計測装置において、
前記発光素子は、単一ピーク波長の赤外光又は単一ピーク波長の赤色光を照射することを特徴とする脈拍計測装置。
【請求項１０】
請求項１乃至９のいずれか一項に記載の脈拍計測装置において、
前記第１フィルタ及び第２フィルタが抽出した信号成分を各々増幅する増幅部と、増幅後の各信号成分の信号レベルに応じて各々の増幅率を可変制御する増幅率変更部とを備えることを特徴とする脈拍計測装置。
【請求項１１】
生体の表面に装着される生体情報検出部を備えた脈拍計測装置の制御方法において、
前記生体情報検出部に設けられた単一の発光素子により、前記生体の血流部位に向けて所定の時間間隔で光を照射し、
前記生体情報検出部に設けられた単一の受光素子により、前記光を前記生体を介して受光すると共に、前記生体の体動によって前記生体と前記生体情報検出部との間に形成された所定の隙間を通って外部から入射した外光を受光して、前記生体の脈動及び体動に応じて受光量が変動する受光信号を生成し、
第１フィルタにより、前記受光信号から、前記時間間隔に対応する発光周波数を含む第１周波数範囲の第１信号成分を抽出し、
第２フィルタにより、前記受光信号から、前記脈動の周波数より低い第２周波数範囲の第２信号成分を抽出し、
前記第１信号成分に周波数解析を施して複数の周波数スペクトルを特定すると共に、前記第２信号成分に周波数解析を施して一又は複数の周波数スペクトルを特定し、
前記第１信号成分に含まれるが、前記第２信号成分には含まれない周波数スペクトルを脈動スペクトルと特定し、この脈動スペクトルの周波数に基づいて脈拍数を算出する
ことを特徴とする脈拍計測装置の制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光学式の脈拍計測装置及びその制御方法に関する。
【背景技術】
【０００２】
従来より、光学的センシングによる脈拍計測装置が知られている。この種の光学式脈拍計測装置には、フィルタを用いて計測信号内の脈拍成分を抽出するものがあり、例えば、発光素子をｄｕｔｙ駆動して受光側で増幅した後にローパスフィルタによって脈拍波形を再現する技術（例えば、特許文献１参照）、受光後に脈拍周波数帯域を通過するバンドパスフィルタによって脈拍波形を抽出する技術（例えば、特許文献２、３参照）、及び、受光した光の変調成分だけを復調して外乱光の影響を除去する技術（例えば、特許文献４参照）がある。
【０００３】
また、この種の光学式脈拍計測装置には、計測信号内の雑音（体動成分）を検出するものがあり、例えば、波長の異なる２つの発光素子と、第１波長光、第２波長光、非発光を繰り返した際の各受光信号を各々積分する複数の積分手段とを備え、非発光時の雑音レベルが許容範囲内か否かを判断する技術（例えば、特許文献５参照）、及び、脈拍信号を通過させるフィルタ（０．５Ｈｚ〜２．３Ｈｚ）と、脈拍よりも低い周波数の信号成分を通過させる雑音成分通過フィルタとを備え、雑音発生期間検出と雑音発生検出のときに脈波信号を通過させるフィルタの出力を遮断する技術（例えば、特許文献６参照）がある。さらに、この種の脈拍計測装置には、複数のセンサを用いて脈拍と体動を検出するものが提案されている（例えば、特許文献７、８、９）。
【特許文献１】実開昭５３−７６６９７号公報
【特許文献２】特開２００１−１４５６０６号公報
【特許文献３】特開２００４−２０２１９０号公報
【特許文献４】特許第２６９３９５８号公報
【特許文献５】特許第３２９１５８１号公報
【特許文献６】特開平１−２８８２３０号公報
【特許文献７】特開平７−８８０９２号公報
【特許文献８】特開２００５−２８１５７号公報
【特許文献９】特開平１１−２７６４４８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかし、従来の構成では、体動による検出部位（血流部位）との間の密着性の変化によって外乱光が侵入したり、外乱光が生体（手首や指）を通過すること等によって脈拍波形の周波数に近い雑音成分（体動成分）が存在した場合、かかる雑音成分を除去できない問題があった。
また、上記特許文献５、７〜９記載の技術は、複数の発光素子、複数の受光素子或いは専用の体動センサを備えるため、構成部品の数が多くなってしまう問題もある。
【０００５】
本発明は、上述した事情に鑑みてなされたものであり、脈拍検出精度が向上し、かつ、発光素子及び受光素子の数を低減可能な脈拍計測装置及びその制御方法を提供することにある。
【課題を解決するための手段】
【０００６】
上述課題を解決するため、本発明は、生体の表面に装着される生体情報検出部を備えた脈拍計測装置において、前記生体情報検出部に設けられ、前記生体の血流部位に向けて所定の時間間隔で光を照射する単一の発光素子と、前記生体情報検出部に設けられ、前記光を前記生体を介して受光すると共に、前記生体の体動によって前記生体と前記生体情報検出部との間に形成された所定の隙間を通って外部から入射した外光を受光可能に配置されて、前記生体の脈動及び体動に応じて受光量が変動する受光信号を生成する単一の受光素子と、前記受光信号から、前記時間間隔に対応する前記発光素子の発光周波数を含む第１周波数範囲の第１信号成分を抽出する第１フィルタと、前記受光信号から、前記脈動の周波数より低い第２周波数範囲の第２信号成分を抽出する第２フィルタと、前記第１信号成分に周波数解析を施して複数の周波数スペクトルを特定すると共に、前記第２信号成分に周波数解析を施して一又は複数の周波数スペクトルを特定する周波数解析部と、前記第１信号成分に含まれるが、前記第２信号成分には含まれない周波数スペクトルを脈動スペクトルと特定し、この脈動スペクトルの周波数に基づいて脈拍数を算出する脈拍数算出部とを備えることを特徴とする。
【０００７】
この発明によれば、第１フィルタによって、生体を通った発光素子の光だけの第１信号成分を抽出して血流の変動に表れる脈動成分及び体動成分を取り出すことができると共に、第２フィルタによって外光の受光量変化に表れる体動成分を取り出すことができる。そして、第１信号成分に含まれるが、第２信号成分には含まれない周波数スペクトルを脈動スペクトルと特定し、この脈動スペクトルの周波数に基づいて脈拍数を算出するので、脈動スペクトルを精度良く特定でき、脈拍検出精度が向上する。また、この脈拍検出には、単一の発光素子と単一の受光素子を用いるだけなので、脈拍検出精度が高いにも拘わらず、発光素子及び受光素子の数を低減することができる。
この場合、上記脈拍数算出部が、前記第１信号成分及び前記第２信号成分に含まれる同一周波数の周波数スペクトルを体動スペクトルと特定し、体動成分を精度良く検出するようにしてもよい。
【０００８】
上記構成において、前記第１フィルタは、前記発光周波数を略中心周波数としたバンドパスフィルタであることが好ましい。この構成によれば、生体を通った発光素子の光だけの第１信号成分を精度良く取り出すことができる。
また、上記構成において、前記第２フィルタは、前記生体の体動周波数範囲を通過させるローパスフィルタであることが好ましい。この構成によれば、外光の受光量変化に表れる体動成分を精度良く取り出すことができる。
【０００９】
また、上記構成において、前記生体の体動によって前記受光素子への外光入射光量が変化するように外光を導入する外光導入部を備えることが好ましい。この構成によれば、受光素子が体動に応じた入射光量で外光をより確実に受光することができる。
また、前記発光周波数を１ｋＨｚを超える周波数にすることが好ましい。この構成によれば、赤外線方式の家電リモコンで一般に使用される周波数信号を上記第１フィルタで取り除くことができる。
【００１０】
また、上記構成において、前記第１フィルタが抽出した信号成分の包絡線を検波する検波部を有し、前記周波数解析部は、前記検波部の出力信号に周波数解析を施して、脈動成分及び体動成分に対応する周波数スペクトルを特定することが好ましい。この場合、前記検波部は、前記第１フィルタが抽出した信号間の極大値同士或いは極小値同士をピークホールドし、それらを結んで前記包絡線を形成してもよい。
【００１１】
また、上記構成において、前記周波数解析部に入力する信号を、前記第１フィルタが抽出した第１信号成分と、前記第２フィルタが抽出した信号成分とに選択的に切り換える入力選択部を有することが好ましい。この構成によれば、一つの周波数解析部を、第１信号成分及び第２信号成分に周波数解析を施すものとして共用することができ、構成部品を低減することができる。また、上記構成において、前記発光素子は、単一ピーク波長の赤外光又は単一ピーク波長の赤色光を照射することが好ましい。この構成によれば、体内で吸収されにくい光を照射することができる。
【００１２】
また、上記構成において、前記第１フィルタ及び第２フィルタが抽出した信号成分を各々増幅する増幅部と、増幅後の各信号成分の信号レベルに応じて各々の増幅率を可変制御する増幅率変更部とを備えることが好ましい。この構成によれば、周波数解析処理するのに適切な信号レベルに調整することができる。
【００１３】
また、本発明は、生体の表面に装着される生体情報検出部を備えた脈拍計測装置の制御方法において、前記生体情報検出部に設けられた単一の発光素子により、前記生体の血流部位に向けて所定の時間間隔で光を照射し、前記生体情報検出部に設けられた単一の受光素子により、前記光を前記生体を介して受光すると共に、前記生体の体動によって前記生体と前記生体情報検出部との間に形成された所定の隙間を通って外部から入射した外光を受光して、前記生体の脈動及び体動に応じて受光量が変動する受光信号を生成し、第１フィルタにより、前記受光信号から、前記時間間隔に対応する発光周波数を含む第１周波数範囲の第１信号成分を抽出し、第２フィルタにより、前記受光信号から、前記脈動の周波数より低い第２周波数範囲の第２信号成分を抽出し、前記第１信号成分に周波数解析を施して複数の周波数スペクトルを特定すると共に、前記第２信号成分に周波数解析を施して一又は複数の周波数スペクトルを特定し、前記第１信号成分に含まれるが、前記第２信号成分には含まれない周波数スペクトルを脈動スペクトルと特定し、この脈動スペクトルの周波数に基づいて脈拍数を算出することを特徴とする。
【００１４】
この発明によれば、第１フィルタによって、生体を通った発光素子の光だけの第１信号成分を抽出して血流の変動に表れる脈動成分及び体動成分を取り出すことができると共に、第２フィルタによって外光の受光量変化に表れる体動成分を取り出すことができる。そして、第１信号成分に含まれるが、第２信号成分には含まれない周波数スペクトルを脈動スペクトルと特定し、この脈動スペクトルの周波数に基づいて脈拍数を算出するので、脈動スペクトルを精度良く特定でき、脈拍検出精度が向上する。また、単一の発光素子と単一の受光素子を用いるだけなので、脈拍検出精度が高いにも拘わらず、発光素子及び受光素子の数を低減することができる。
【００１５】
また、本発明は、以上説明した脈拍計測装置及びその制御方法に適用する他、この発明を実施するための制御プログラムを電気通信回線を介して一般ユーザに配布したり、そのようなプログラムを、磁気記録媒体、光記録媒体、半導体記録媒体といった、コンピュータに読み取り可能な記録媒体に格納して一般ユーザに配布する、といった態様でも実施され得る。
【発明の効果】
【００１６】
本発明によれば、第１フィルタによって、生体を通った発光素子の光だけの第１信号成分を抽出して血流の変動に表れる脈動成分及び体動成分を取り出すと共に、第２フィルタによって外光の受光量変化に表れる体動成分を取り出し、第１信号成分に含まれるが、第２信号成分には含まれない周波数スペクトルを脈動スペクトルと特定し、この脈動スペクトルの周波数に基づいて脈拍数を算出するので、脈拍検出精度が向上し、かつ、発光素子及び受光素子の数を低減することができる。
【発明を実施するための最良の形態】
【００１７】
以下、図面を参照して本発明の実施形態を詳述する。
＜第１実施形態＞
図１は、本発明の第１実施形態に係る脈拍計測装置を示す図である。この脈拍計測装置１は、腕時計型に構成され、装置本体１０と、この装置本体１０の６時位置及び１２時位置から延びてユーザの腕（手首）に巻回されるリストバンド（帯状体）１１とを備え、ジョギングやランニング等の運動を行う際にユーザが容易に装着可能に構成されている。
【００１８】
装置本体１０は、脈拍数や時刻等の各種情報を表示する表示部１５と、ユーザが各種指示を行うための操作部として機能する複数の操作スイッチ１６Ａ、１６Ｂ、１６Ｃとを有し、その内部に、生体情報検出部２０等を構成する各種電気部品や当該脈拍計測装置１の動作電力を供給する電池等を備えている。
操作スイッチ１６Ａは、脈拍計測を開始／停止させる操作子であり、例えば、操作スイッチ１６Ａを操作する毎に、動作モードが脈拍計測モード（図１参照）と、時刻表示モードとに切り替わる。また、操作スイッチ１６Ｂは、時刻設定等の各種セットアップの開始指示や動作モードの変更指示を行う操作子であり、操作スイッチ１６Ｃは、セットアップのクリアや図示せぬライトをオンさせる操作子である。なお、図示の例では、操作スイッチ１６Ａ〜１６Ｃが押下式スイッチの場合を示しているが、静電式やメンブレン電極を用いたタッチ式スイッチでもよい。
【００１９】
生体情報検出部２０は、ユーザの腕（手首：血流部位）Ｘの表面に配置され、図２（Ａ）に示すように、一個の発光素子２１と、一個の受光素子２２とを備えて構成される光学式センサである。発光素子２１は、７００ｎｍ以上の光Ｌ１を出射するＬＥＤ、具体的には、単一ピーク波長の可視光に比べて体内で吸収されにくい赤外光、又は、単一ピーク波長の赤色光を出射するＬＥＤが適用される。この発光素子２１は、図２（Ｂ）に示すように、装置本体１０の裏側に設けられた透明カバー（図示略）を介してユーザの腕Ｘに向けて光Ｌ１を照射する。
【００２０】
受光素子２２は、発光素子２１の波長領域内の光及び外光Ｌ３を受光可能なフォトトランジスタ（例えば、受光波長領域が７００ｎｍ以上のフォトトランジスタ）が適用され、太陽光の直射によっても受光素子２２が電気的に飽和しないレベルに光を抑制する光学フィルタ２２Ａを備えている。この受光素子２２は、図２（Ａ）に示すように、装置本体１０の筐体の中心を外した位置に設けられ、図２（Ｂ）に示すように、装置本体１０の裏側に設けられた透明カバー（図示略）を介して腕Ｘの血管Ｋ等で反射された反射光Ｌ２を受光し、受光量に応じた信号レベルの受光信号Ｓ１を出力する。
【００２１】
ここで、発光素子２１から照射された光Ｌ１は、血管Ｋを流れる血流に応じて吸光度が変化するため、その反射光Ｌ２（透過光）の光量が血流に応じて変化する。この場合、血管Ｋ内の血流は、脈動だけでなく、体動（腕振り等の身体自体の動きや手首等の関節を曲げたときに変化する皮膚表面の動きがある）によっても変化するため、血管Ｋに届いた光Ｌ１の反射光Ｌ２は、脈動及び体動に応じて光量が変化することになる。
受光素子２２の受光波長領域を７００ｎｍ以上にした理由は、発光素子２１の光Ｌ１の反射光を確実に受光し、かつ、太陽光や照明光等の外光Ｌ３についても受光可能にするためである。なお、本実施形態では、反射光を受光する反射型に構成する場合について説明するが、これに限らず、この生体情報検出部２０により指等の血管から脈拍を検出する場合は、生体（指等）の通過光を受光する透過型に構成してもよい。
また、上述したように、本構成では、受光素子２２に外光Ｌ３を受光可能な素子を使用したため、太陽光や照明光等の外光Ｌ３を受光波長帯域に含まない受光素子を使用した従来の光学式脈拍計測装置に比して、汎用の受光素子を広く適用することができる。また、この受光素子２２には、汎用のフォトトランジスタに限らず、汎用のフォトダイオード等の広く流通する光電センサを適用することができる。
【００２２】
この脈拍計測装置１のリストバンド１１は、図２（Ｂ）に示すように、装置本体１０の裏側に延在し、装置本体１０の裏側部１１Ａ及びこの裏側部１１Ａに連続する連続部１１Ｂ、１１Ｃが、透明樹脂又は半透明樹脂等の透光部材で形成されている。このため、脈拍計測装置１をユーザの腕Ｘに装着した場合、上記透光部材からなる裏側部１１Ａ及び連続部１１Ｂ、１１Ｃが、受光素子２２と生体（腕Ｘ）との間に位置し、受光素子２２と生体（腕Ｘ）との間に外光Ｌ３を導入し易くすることができる。従って、図２（Ｂ）に示すように、ランニング時のユーザの腕振り等によってリストバンド１１（生体情報計測部２０）と腕Ｘとの間に所定の隙間Ｙが空いた状態では、太陽光や照明光等の外光Ｌ３が、上記隙間Ｙ、上記裏側部１１Ａ及び連続部１１Ｂ、１１Ｃを通過して、腕Ｘ等で反射した後に受光素子２２で受光されることとなる。但し、上記隙間Ｙが空いた場合でも、脈拍計測装置１の姿勢や外光Ｌ３の入射角度によっては外光Ｌ３が受光素子２２に殆ど受光されない場合がある。
一方、リストバンド１１（生体情報計測部２０）と腕Ｘとの間に殆ど隙間が空いていない状態では、外光Ｌ３がリストバンド１１や装置本体１０のケースで遮られて受光素子２２に殆ど受光されない。なお、外光Ｌ３の特に赤外光成分及び赤色光成分については、腕Ｘ等の生体で吸収されにくいために生体（腕Ｘ、生体表面（皮膚等））を通過し易く、かかる外光Ｌ３の生体通過光についても受光素子２２で受光される場合がある。
【００２３】
すなわち、これら透光部材で形成された裏側部１１Ａ及び連続部１１Ｂ、１１Ｃは、外光Ｌ３を受光素子２２の受光面側に導入させる外光導入部３０として機能する。以下の説明では、これらを外光導入部３０と表記する。
なお、外光導入部３０は、透光部材に限らず、リストバンド１１に形成された開口にしてもよく、要は、外光Ｌ３を、体動に応じて受光量が異なるように受光素子２２に導入可能であれば、任意に設計可能である。また、リストバンド１１は、外光導入部３０以外を着色部材にしても良いし、或いは、外光導入部３０以外も透光部材にしてもよい。
【００２４】
図３は脈拍計測装置１のブロック図であり、図４はその回路構成を示す図である。図３において、発光制御部５０は、発光素子２１の発光制御を行うものであり、図示せぬ発振回路の発振信号（例えば３２ｋＨｚ）を分周（１／１６分周）して得た２ｋＨｚの信号を入力し、この２ｋＨｚの信号がＨレベルのときに発光素子２１に駆動電力を供給し、Ｌレベルのときに駆動電力の供給を停止し、これによって、２ｋＨｚデューティ比５０％に対応する時間間隔で発光／消灯を繰り返させる。つまり、この脈拍計測装置１では、２ｋＨｚを発光の時間間隔を規定する発光周波数（間欠駆動周波数とも言う）としている。
【００２５】
電流電圧変換部５１は、受光素子２２から出力される受光信号Ｓ１を電流／電圧変換するものであり、図４に示すように、オペアンプ５１Ａの出力と入力−を抵抗５１Ｂを介して接続すると共にコンデンサ５１Ｃを並列接続して構成され、受光量に比例した電流を電圧に変換してオペアンプ５１Ａにより増幅して出力する。
【００２６】
この電流電圧変換部５１で変換された受光信号Ｓ２は、帯域通過フィルタ部（第１フィルタ）５２と低域通過フィルタ部（第２フィルタ）５３とに各々出力される。
帯域通過フィルタ部５２は、電流電圧変換部５１から出力される受光信号Ｓ２から、上記発光周波数（２ｋＨｚ）を含む第１周波数範囲の信号成分（第１信号成分）を抽出するものである。より具体的には、帯域通過フィルタ部５２は、上記発光周波数を中心周波数（ｆ０）とし、かつ、低域通過フィルタ部５３を通過する低周波数範囲より高い周波数範囲（本構成では、１ｋＨｚ以上）の信号成分を通過させる。これによって、生体を通った光（光Ｌ１の反射光Ｌ２）の信号成分を取り出すことができ、言い換えれば、血流の変動（脈動及び体動による変動）に応じて振幅レベルが変動する信号成分を取り出すことができる。
【００２７】
この帯域通過フィルタ部５２は、図４に示すように、オペアンプ５２Ａの出力と入力−を抵抗５２Ｂを介して接続すると共にコンデンサ５２Ｃ、５２Ｄを介して接続し、かつ、入力−への入力側に抵抗５２Ｅ、５２Ｆを接続した多重帰還型アクティブフィルタ回路で構成される。このため、パッシブフィルタ回路で構成する場合に比して、得られる周波数特性の自由度が高く、かつ、増幅を行うことができるといった利点がある。
ここで、帯域通過フィルタ部５２の通過周波数を１ｋＨｚ以上にした理由は、赤外線方式の家電リモコンで一般に使用される１ｋＨｚ以下の周波数信号（ノイズ成分に相当）を除くためである。
【００２８】
低域通過フィルタ部５３は、電流電圧変換部５１から出力される受光信号から、脈動の周波数より低い第２周波数範囲の信号成分を抽出するものであり、言い換えれば、体動に応じて受光素子２２に入射した外光Ｌ３による体動成分を抽出可能なフィルタに構成されている。
ここで、体動成分（歩行時やランニング時等の腕振り等の身体運動や手首等の関節を曲げたときに変化する皮膚表面の動きがある）は、通常０Ｈｚよりも大きく１０Ｈｚ以下の範囲となるため、この低域通過フィルタ部５３は、例えば、カットオフ周波数（ｆｃ）が１０Ｈｚのフィルタに構成される。この低域通過フィルタ部５３は、図４に示すように、抵抗５３Ａ、５３Ｂ、コンデンサ５３Ｃ、５３Ｄ及びオペアンプ５３Ｅを組み合わせたアクティブフィルタ回路で構成され、パッシブフィルタ回路で構成する場合に比して、得られる周波数特性の自由度が高く、かつ、増幅を行うことが可能である。
【００２９】
図３に示すように、帯域通過フィルタ部５２の後段（出力側）には、増幅率変更部５４によって増幅率を変更可能な増幅部（第１増幅部）５５と、検波部（第１検波部）５６とが順に接続され、低域通過フィルタ部５３の後段（出力側）には、増幅率変更部５７によって増幅率を変更可能な増幅部（第２増幅部）５８が接続される。
【００３０】
増幅部５５及び増幅率変更部５４は、帯域通過フィルタ部５２を通過した受光信号Ｓ３を演算処理部６０で演算処理（周波数解析処理）するのに適切な信号レベルに増幅するためのものである。この増幅部５５及び増幅率変更部５４は、図４に示すように、抵抗５５Ａ、５５Ｂを前後に配置したオペアンプ５５Ｃの出力と入力−をつなぐ帰還経路に、抵抗５４Ａ、５５Ｃ、５４Ｃを各々スイッチング素子５４Ｄ、５４Ｅ、５４Ｆを介して接続可能に構成した回路に構成され、演算処理部６０の制御の下、各スイッチング素子５４Ｄ〜５４Ｆのオン／オフが制御されることによって増幅率を変更する。
【００３１】
検波部５６は、増幅部５５で増幅された受光信号Ｓ３の包絡線を検波するものであり、図４に示すように、入力信号をショットキーダイオード５６Ａに通し、抵抗５６Ｂとコンデンサ５６Ｃの並列回路で受ける回路が適用される。この抵抗５６Ｂとコンデンサ５６Ｃの時定数を適切に設計することで、その出力波形を入力波形の包絡線に近い波形にすることができる。
【００３２】
なお、上記検波部５６に代えて、上記受光信号Ｓ３の極大値同士或いは極小値同士をピークホールドし、それらを結んで包絡線を検波する検波部を設けるようにしてもよい。この場合、この検波部は、上記受光信号Ｓ３をアナログデジタル変換した波形データを取得するデジタルメモリスコープを備え、このデジタルメモリスコープに取得された波形データの極大値同士或いは極小値同士をピークホールドして得た包絡線を取得するようにすればよい。このように、極大値同士或いは極小値同士を結んだ包絡線を得ることにより、ノイズ成分の影響を低減して包絡線を検波することができる。この場合、この包絡線には、脈動そのものの波形ではなく、脈動の１／２周波数の波形が含まれることになる。
【００３３】
増幅部５８及び増幅率変更部５７は、低域通過フィルタ部５３を通過した受光信号Ｓ４を演算処理部６０で演算処理（周波数解析処理）するのに適切な信号レベルに増幅するためのものであり、上記増幅部５５及び増幅率変更部５４と同様の回路が適用される。つまり、増幅部５８及び増幅率変更部５７は、図４に示すように、抵抗５８Ａ、５８Ｂを前後に配置したオペアンプ５８Ｃの出力と入力−をつなぐ帰還経路に、抵抗５７Ａ、５７Ｃ、５７Ｃを各々スイッチング素子５７Ｄ、５７Ｅ、５７Ｆを介して接続可能に構成した回路に構成され、演算処理部６０の制御の下、各スイッチング素子５７Ｄ〜５７Ｆのオン／オフが制御されることによって増幅率を変更する。
【００３４】
入力選択部６１は、演算処理部６０の制御の下、増幅部５５で増幅された受光信号Ｓ３、及び、増幅部５８で増幅された受光信号Ｓ４のいずれかをＡ／Ｄ変換部６２に切換出力するものであり、図３に示すように、増幅部５５及び５８の出力側にスイッチング素子６１Ａ、６１Ｂを各々配置し、一方のスイッチング素子６１Ａには、演算処理部６０からＡ／Ｄ入力選択信号ＳＳを出力してオンオフを切り換え、他方のスイッチング素子６１Ｂには、上記Ａ／Ｄ入力選択信号ＳＳの信号レベルをインバータ６１Ｃにより反転させて供給してオフオンを切り換える回路が適用される。
Ａ／Ｄ変換部６２は、受光信号Ｓ３の包絡線信号Ｓ３Ａ又は受光信号Ｓ４のいずれかを入力し、アナログデジタル変換して演算処理部６０に出力する。
【００３５】
演算処理部６０は、ＣＰＵ、ＲＯＭ及びＲＡＭ（本例では、ＦＦＴ用ＲＡＭ１と、ＦＦＴ用ＲＡＭ２とを有する）等を備えたコンピュータ構成を備え、ＣＰＵがＲＯＭに格納された制御プログラムを実行することにより、入力したデータの周波数解析処理（ＦＦＴ処理）を行う周波数解析部、周波数解析結果に基づき脈拍数を算出する脈拍数算出部、特定した脈拍数を表示部１５に表示させる処理等の脈拍計測装置１の各部を制御する制御部として機能する。
【００３６】
なお、本実施形態では、演算処理部６０がソフトウェア処理により周波数解析部、脈拍数算出部及び制御部として機能する場合を説明するが、これに限らず、これら周波数解析部、脈拍数算出部及び制御部を構成するハードウェア回路によって構成してもよい。また、この演算処理部６０には、図示せぬ発振回路の発振信号の分周信号（１Ｈｚ信号）に基づいて時刻を計時する時計回路も内蔵されている。時計回路の構成は公知の構成を適用すればよいため、詳細な説明は省略する。
【００３７】
次に、この脈拍計測装置１の脈拍計測時の動作を説明する。図５はこの場合の動作を示すフローチャートである。
ユーザが操作スイッチ１６Ａを操作して脈拍計測の開始を指示すると、まず、演算処理部６０内のＣＰＵ（以下、制御部という）は、図示せぬ発振回路の発振信号を分周した出力（本例では３２Ｈｚの信号）の割り込みを許可し（ステップＳＰ１）、この３２Ｈｚ割り込みをトリガーとしてＡ／Ｄ変換部６２の取り込みと、入力選択部６１による入力切り換えを行う。また、制御部は、変数Ｎを値０に設定し（ステップＳ２）、その後、発光制御部５０により２ｋＨｚ（発光周波数）のデューティ比５０％の時間間隔で発光素子２１の発光制御を開始させる。
【００３８】
このため、受光素子２２から受光信号Ｓ１が出力され、電流電圧変換部５１によって電流電圧変換されて、図６（Ａ）に示すように、２ｋＨｚの変調波である受光信号Ｓ２が出力される。そして、この受光信号Ｓ２は、２ｋＨｚを中心周波数とする帯域通過フィルタ部５２を通過することによって、図６（Ｂ）に示すように、２ｋＨｚの変調成分を示す受光信号Ｓ３が出力される。次に、受光信号Ｓ３に対し、検波部５６が検波処理を施すことによって、図６（Ｃ）に示すように、受光信号Ｓ３の包絡線を示す包絡線信号Ｓ３Ａが出力される。
この包絡線信号Ｓ３Ａの波形には、脈動成分（脈波形）と体動成分（体動波形）の情報が含まれている。
【００３９】
また、受光信号Ｓ２が、１０Ｈｚ以下を通過する低域通過フィルタ部５３を通過することによって、脈動成分が取り除かれた受光信号Ｓ４が取得され、この受光信号Ｓ４の波形には、体動成分とノイズ成分（体動に依存しない外乱光成分等）の情報が含まれる。
次に、上記ステップＳ２の処理実行後、制御部は、Ａ／Ｄ入力選択信号ＳＳをＨレベルに設定し、３２Ｈｚ割り込みフラグがオンになるまで待機する（ステップＳＰ４）。なお、この３２Ｈｚ割り込みフラグとは、上記３２Ｈｚの割り込み信号の１周期（１／３２秒）が経過したか否かを判定するためのフラグであり、言い換えれば、Ａ／Ｄ変換部６２の取り込みと入力切り換えとを行うタイミングを検出するフラグである。
【００４０】
Ａ／Ｄ入力選択信号ＳＳがＨレベルに設定された場合、入力選択部６１のスイッチング素子６１Ａが開状態に設定されると共に、スイッチング素子６１Ｂが閉状体に設定されるため、包絡線信号Ｓ３ＡがＡ／Ｄ変換部６２に出力され、ここでアナログデジタル変換されて演算処理部６０に出力される。
そして、制御部は、３２Ｈｚ割り込みフラグがオンになると（ステップＳＰ５：ＹＥＳ）、３２Ｈｚ割り込みフラグをクリアして（ステップＳＰ６）、デジタルデータ（包絡線の波形データ）をＦＦＴ用ＲＡＭ１に格納する（ステップＳＰ７）。
【００４１】
続いて、制御部は、Ａ／Ｄ入力選択信号ＳＳをＬレベルに切り換え、入力選択部６１のスイッチング素子６１Ａを閉状態に切り換えると共に、スイッチング素子６１Ｂを開状態に切り換える。このため、低域通過フィルタ部５３を通過した受光信号Ｓ４が、Ａ／Ｄ変換部６２に入力されてアナログデジタル変換されて演算処理部６０に出力される。そして、制御部は、３２Ｈｚ割り込みフラグがオンになると（ステップＳＰ８：ＹＥＳ）、３２Ｈｚ割り込みフラグをクリアして（ステップＳＰ９）、デジタルデータ（受光信号Ｓ４の波形データ）をＦＦＴ用ＲＡＭ２に格納する（ステップＳＰ１０）。
【００４２】
次に、制御部は、変数Ｎに値１を加算し（ステップＳＰ１１）、変数Ｎが値２５６に達したか否かを判定する（ステップＳＰ１２）。このとき、値２５６に達していない場合は（ステップＳＰ１２：ＮＯ）、ステップＳＰ３の処理へ移行して、包絡線の波形データの取り込みと、受光信号Ｓ４の波形データの取り込みとを再び開始させる。これによって、これらデータの取り込みは変数Ｎが値２５６に達するまで繰り返され、ＦＦＴ用ＲＡＭ１及びＦＦＴ用ＲＡＭ２に、各１６Ｈｚ２５６ポイントの波形データを格納するまで上記データの取り込みを繰り返す。
【００４３】
続いて、変数Ｎが値２５６に達すると（ステップＳＰ１２：ＹＥＳ）、制御部は、ＦＦＴ用ＲＡＭ１に格納された波形データ（受光信号Ｓ３の包絡線の波形データ）に対し、ＦＦＴ（高速フーリエ変換）処理（周波数解析）を実行することにより（ステップＳＰ１３）、図７（Ａ）に示すように、該波形に含まれる複数の周波数スペクトル（ローカルピーク）を特定可能な第１の周波数解析結果を得る。
次いで、制御部は、ＦＦＴ用ＲＡＭ２に格納された波形データ（受光信号Ｓ４の波形データ）に対し、同様のＦＦＴ（高速フーリエ変換）処理（周波数解析）を実行することにより（ステップＳＰ１４）、図７（Ｂ）に示すように、該波形に含まれる１又は複数の周波数スペクトル（ローカルピーク）を特定可能な第２の周波数解析結果を得る。
【００４４】
ここで、上述したように、ＦＦＴ用ＲＡＭ１に格納された波形データには、脈動成分（脈波形）と体動成分（体動波形）の情報が含まれることから、図７（Ａ）に示す２つの周波数スペクトルｆ１、ｆ２は、脈動成分と体動成分とに相当することが判る。しかし、この情報だけでは、どちらの周波数スペクトルｆ１、ｆ２が脈動成分であるかを特定することが困難である。
一方、ＦＦＴ用ＲＡＭ２に格納された波形データには、上述したように、体動成分とノイズ成分が含まれ、脈動成分（脈波形）が含まれないことから、図７（Ｂ）に示す大きな周波数スペクトルｆ３、ｆ４、ｆ５は、体動成分とノイズ成分のいずれかであることが判る。
【００４５】
このため、制御部は、上記第１の周波数解析結果と上記第２の周波数解析結果を比較することによって、両方の周波数解析結果に含まれる同一周波数の周波数スペクトルｆ２、ｆ３を体動成分（体動スペクトル）と特定し、第１の周波数解析結果に含まれるが、第２の周波数解析結果には含まれない周波数スペクトルｆ１を脈動成分（脈動スペクトル）と特定することができる（ステップＳＰ１５）。そして、制御部は、この周波数スペクトルｆ１の周波数を一分間当たりの振幅回数に換算して脈拍数を算出し、算出した脈拍数を表示部１５に表示させ（ステップＳ１６）、ステップＳＰ２の処理に移行する。以上説明した動作は、ユーザから脈拍計測の停止が指示されるまで継続的に繰り返され、リアルタイムの脈拍数を表示部１５に表示することができる。以上が、脈拍計測時の動作である。
【００４６】
なお、上述したように、周波数スペクトルｆ２、ｆ３を体動成分、周波数スペクトルｆ１を脈動成分と特定できた場合には、周波数スペクトルｆ４、ｆ５、ｆ６はノイズ成分であると特定できる。このため、これらノイズ成分を記憶して以降はノイズと推定する学習処理を行うことで、以降の検出精度向上に役立てることが可能である。
一方、ユーザから脈拍計測の停止が指示され、時刻表示の開始が指示された場合には、制御部は、時計回路で計時されている時刻を取得して表示部１５に表示させる。この場合、発光素子２１の駆動や受光素子２２からの電流電圧変換部５１への信号出力等を停止し、電力消費量の節約を図るようにしている。
【００４７】
以上説明したように、本実施形態の脈拍計測装置１は、所定の時間間隔（発光周波数に対応）で光Ｌ１を照射する単一の発光素子２１と、その光Ｌ１を生体を介して受光可能で、かつ、ユーザの体動に連動して外光Ｌ３の受光量が変化可能に設けられた単一の受光素子２２と、受光素子２２が出力する受光信号から、上記時間間隔に対応する発光周波数を中心周波数とする第１周波数範囲の第１信号成分を抽出する帯域通過フィルタ部５２と、上記受光信号から、脈動の周波数より低い第２周波数範囲の第２信号成分を抽出する低域通過フィルタ部５３とを備えるので、帯域通過フィルタ部５２によって、生体を通った光（光Ｌ１の反射光Ｌ２）だけの信号成分を取り出して血流の変動に表れる脈動成分及び体動成分を取り出すことができると共に、低域通過フィルタ部５３によって外光Ｌ３の受光量変化に表れる体動成分を取り出すことができる。
【００４８】
従って、帯域通過フィルタ部５２を通過した第１信号成分から脈動成分と体動成分とのいずれかに対応する周波数スペクトルｆ１、ｆ２を特定し、かつ、低域通過フィルタ部５３を通過した第２信号成分から体動成分の周波数スペクトルｆ３を特定する演算処理部６０を設け、この演算処理部６０が、第１信号成分に含まれるが、第２信号成分には含まれない周波数スペクトルｆ１を脈動スペクトルと特定してこの脈動スペクトルの周波数に基づいて脈拍数を算出するので、脈動成分ｆ１及び体動成分ｆ２（＝ｆ３）を精度良く検出することができ、脈拍検出精度及び体動検出精度を向上させることができる。
また、本構成では、発光素子２１及び受光素子２２を各々一個ずつ備えるだけでよいので、従来のものと比較して、脈拍検出精度が高いにも拘わらず、発光素子及び受光素子の数を低減することができる。
【００４９】
さらに、本構成では、リストバンド１１に、ユーザの体動によって受光素子２２への外光入射光量が変化するように外光Ｌ３を導入する外光導入部３０を設けたので、受光素子２２が体動に応じた入射光量で外光Ｌ３をより確実に受光することができる。
なお、この外光導入部３０を設けない構成にすることも可能である。外光導入部３０を設けない構成でも、腕時計型の脈拍計測装置１といった身体装着タイプでは、ユーザの腕振り等の運動によって生体情報計測部２０と腕Ｘとの間等に所定の隙間Ｙが空くため、この隙間Ｙを通過した外光Ｌ３が受光素子２２で受光される場合があるからである。また、外光Ｌ３が生体（腕Ｘ、生体表面（皮膚等））を通過して受光素子２２に届く場合もあるため、これによっても外光Ｌ３が受光素子２２で受光される場合が考えられる。但し、外光Ｌ３を確実に受光素子２２に受光させるには、外光導入部３０を設けた方が確実であることは明らかである。
【００５０】
また、本構成では、図２（Ａ）に示すように、受光素子２２を装置本体１０の中心を外した位置に設置することによって装置本体１０の重心位置を中心位置からずらすようにしている。これによって、ユーザの腕振り等の運動によって生体情報計測部２０と腕Ｘとの間等に隙間Ｙが空き易くなり、外光Ｌ３を受光素子２２で受光され易くして体動の検出精度を高めることができる。
【００５１】
さらに、本構成では、生体を通った光（光Ｌ１の反射光Ｌ２）だけの信号成分を取り出す第１フィルタを帯域通過フィルタ部５２に構成したので、ノイズ成分に相当する高周波成分を予め除去することが可能である。
また、帯域通過フィルタ部５２を通過した第１信号成分と、低域通過フィルタ部５３を通過した第２信号成分とを、Ａ／Ｄ変換部６２に切換出力する入力選択部６１を設けたので、これら信号成分の信号処理に使用するＡ／Ｄ変換部６２及び演算処理部６０（周波数解析部）を各々一つで共用でき、その分、構成部品を低減することができる。
なお、上記実施形態では、発光周波数を２ｋＨｚとする場合を例示したが、これに限らない。例えば、発振回路に時計用に一般に使用される水晶発振回路を使用した場合、その分周周波数である２０４８Ｈｚ或いは１０２４Ｈｚの整数倍の周波数を使用することが好ましい。
【００５２】
＜第２実施形態＞
図８は第２実施形態に係る脈拍計測装置１のブロック図である。この脈拍計測装置１は、帯域通過フィルタ部５２を通過した信号（包絡線信号Ｓ３Ａ）専用のＡ／Ｄ変換部６２Ａ及び第１周波数解析部（ＦＦＴ）７０Ａを備え、かつ、低域通過フィルタ部５３を通過した信号（受光信号Ｓ４）専用のＡ／Ｄ変換部６２Ｂ及び第２周波数解析部（ＦＦＴ）７０Ｂを備える点が、第１実施形態に係る脈拍計測装置１と異なる。これ以外の構成は第１実施形態と略同一の構成であるため、同一の符号を付して示し、重複する説明は省略する。
【００５３】
Ａ／Ｄ変換部６２Ａは、帯域通過フィルタ部５２を通過して検波部５６から出力される受光信号Ｓ３の包絡線信号Ｓ３Ａを入力し、これをアナログデジタル変換して第１周波数解析部７０Ａに出力する。そして、第１周波数解析部７０Ａは、入力した包絡線信号Ｓ３Ａのデータに対してＦＦＴ（高速フーリエ変換）処理を施し、この周波数解析結果を演算処理部６０に出力する。
【００５４】
また、Ａ／Ｄ変換部６２Ｂは、低域通過フィルタ部５３を通過した受光信号Ｓ４を入力し、これをアナログデジタル変換して第２周波数解析部７０Ｂに出力する。そして、第２周波数解析部７０Ｂは、入力した受光信号Ｓ４のデータに対してＦＦＴ（高速フーリエ変換）処理を施し、この周波数解析結果を演算処理部６０に出力する。
【００５５】
このように、本実施形態では、帯域通過フィルタ部５２を通過した信号専用のＡ／Ｄ変換部６２Ａ及び第１周波数解析部７０Ａを設けると共に、低域通過フィルタ部５３を通過した信号専用のＡ／Ｄ変換部６２Ａ及び第２周波数解析部７０Ｂを備えるので、計算量の多い周波数解析処理を、専用の周波数解析部７０Ａ、７０Ｂで行うことができ、演算処理部６０（ＣＰＵ）の負担を軽減することができる。これによって、演算処理部６０内のＣＰＵに演算能力が比較的低いものを採用することができ、また、汎用の周波数解析用のＩＣチップを周波数解析部７０Ａ、７０Ｂとして使用することができる。
【００５６】
＜第３実施形態＞
図９は第３実施形態に係る脈拍計測装置１のブロック図である。この脈拍計測装置１は、ユーザの指Ｘ２に向けて発光素子２１の光Ｌ１を照射し、その反射光Ｌ２を受光素子２２で受光する光学式の脈拍計測装置であり、外光導入部３０、増幅率変更部５４及び５７を備えない点が、上記第２実施形態に係る脈拍計測装置１と異なる。それ以外の構成は、これ以外の構成は第２実施形態と略同一の構成であるため、同一の符号を付して示し、重複する説明は省略する。
この脈拍計測装置１においては、増幅部５５及び５８の増幅率が予め定めた一定値に固定される。この固定値は、様々な環境において脈波を計測した場合に最も適切であった値に設定される。この構成によれば、増幅部５５及び５８の増幅率を変更する増幅率変更部等を備えない分、上記第１実施形態及び第２実施形態よりも構成部品を低減することができ、また、演算処理部６０（ＣＰＵ）の負担を軽減することができる。
【００５７】
なお、本発明は上述の実施形態に限定されるものではなく、本発明の目的を達成できる範囲での変形、改良等は本発明に含まれるものである。
例えば、上述の実施形態では、脈拍計測処理を実行するための制御プログラムを脈拍計測装置１内に予め記憶しておく場合について説明したが、この制御プログラムを、磁気記録媒体、光記録媒体、半導体記録媒体等のコンピュータが読み取り可能な記録媒体に格納し、コンピュータが記録媒体からこの制御プログラムを読み取って実行するようにしてもよい。また、この制御プログラムを通信ネットワーク上の配信サーバ等からダウンロードできるようにしてもよい。
また、上述の実施形態では、本発明を腕時計型の脈拍計測装置に適用する場合について説明したが、これに限らない。例えば、発光素子２１と受光素子２２とを備えた光学式センサユニット（生体情報検出部に相当）を人や動物等の生体の血流部位に取り付け、それ以外の構成は、別体とした脈拍計測装置等に広く適用が可能である。また、脈拍検出の血流部位は、手首、指に限らず、耳部等でもよい。
【図面の簡単な説明】
【００５８】
【図１】本発明の第１実施形態に係る脈拍計測装置を示す図である。
【図２】（Ａ）は脈拍計測装置の発光素子と受光素子のレイアウトを示す平面図であり、（Ｂ）はその側断面図である。
【図３】脈拍計測装置のブロック図である。
【図４】脈拍計測装置の回路構成を示す図である。
【図５】脈拍計測時の動作を示すフローチャートである。
【図６】（Ａ）〜（Ｄ）は脈拍計測装置における各部の信号波形を示す図である。
【図７】（Ａ）は帯域通過フィルタ部を通過した信号の周波数解析結果を示す図であり、（Ｂ）は低域通過フィルタ部を通過した信号の周波数解析結果を示す図である。
【図８】本発明の第２実施形態に係る脈拍計測装置のブロック図である。
【図９】本発明の第３実施形態に係る脈拍計測装置のブロック図である。
【符号の説明】
【００５９】
１…脈拍計測装置、１０…装置本体、１１…リストバンド、１５…表示部、２０…生体情報検出部、２１…発光素子、２２…受光素子、３０…外光導入部、５０…発光制御部、５１…電流電圧変換部、５２…帯域通過フィルタ部（第１フィルタ）、５３…低域通過フィルタ部（第２フィルタ）、５４、５７…増幅率変更部、５５、５８…増幅部、５６…検波部、６０…演算処理部（周波数解析部、脈拍数算出部、制御部）、６１…入力選択部、６２…Ａ／Ｄ変換部、ｆ１〜ｆ６…周波数スペクトル。
【図１】
2008220722.tif 000003
【図２】
2008220722.tif 000004
【図３】
2008220722.tif 000005
【図４】
2008220722.tif 000006
【図５】
2008220722.tif 000007
【図６】
2008220722.tif 000008
【図７】
2008220722.tif 000009
【図８】
2008220722.tif 000010
【図９】
2008220722.tif 000011