(57)【要約】      （修正有）【課題】心拍変異度の周波数領域を精密に分析する心拍変異度分析装置及び心拍変異度検出方法を提供する。【解決手段】コンバータ、Ｒ−ピーク検出部、周波数変換部、及び心拍変異度分析部を含む。コンバータは、アナログ心電図信号を受信してデジタル心電図信号に変換する。Ｒ−ピーク検出部はデジタル心電図信号を受信してＲ−ピーク値を検出し、Ｒ−ピーク値に基づいてＲＲ間隔時系列信号を生成する。周波数変換部は周波数移動部、低域通過フィルタ、及び高速フーリエ変換部を含む。周波数移動部は、ＲＲ間隔時系列信号から対象周波数帯域の信号を抽出して対象周波数信号を生成する。心拍変異度分析部は、対象周波数信号に基づいて心拍変異度信号を生成する。【選択図】図３
2017225808.tif 000002
【特許請求の範囲】
【請求項１】
ＱＲＳ波を含むアナログ心電図信号を受信してデジタル心電図信号に変換するコンバータと、
前記デジタル心電図信号を受信してＲ−ピーク値を検出し、前記Ｒ−ピーク値に基づいてＲＲ間隔時系列信号を生成するＲ−ピーク検出部と、
前記ＲＲ間隔時系列信号から対象周波数帯域の信号を抽出して対象周波数信号を生成する周波数変換部と、
前記対象周波数信号に基づいて前記対象周波数帯域に対する心拍変異度信号を生成する心拍変異度分析部と、を含み、
前記周波数変換部は、
前記ＲＲ間隔時系列信号の対象周波数帯域を基底帯域に移動してシフティング検出信号を生成する周波数移動部と、
前記シフティング検出信号を受信し、前記基底帯域のシフティング検出信号を通過させる低域通過フィルタと、
前記低域通過フィルタを通過したシフティング検出信号に基づいて対象周波数信号を生成する高速フーリエ変換部と、を含む心拍変異度分析装置。
【請求項２】
前記周波数移動部は、
Ｉ（Ｉｎ−ｐｈａｓｅ）チャネル信号とＱ（ｑｕａｄｒａｔｕｒｅ−ｐｈａｓｅ）チャネル信号を含む前記シフティング検出信号を生成する請求項１に記載の心拍変異度分析装置。
【請求項３】
前記周波数移動部は、
前記ＲＲ間隔時系列信号に第１正弦信号をかけてＩチャネル信号を生成し、前記ＲＲ間隔時系列信号に第２正弦信号をかけてＱチャネル信号を生成し、
前記第１正弦信号及び前記第２正弦信号は９０度の位相差を有し、前記第１正弦信号及び前記第２正弦信号の周波数は前記対象周波数帯域の中心周波数に基づく請求項２に記載の心拍変異度分析装置。
【請求項４】
前記コンバータは、
サンプリング周波数でサンプリングされた前記デジタル心電図信号を生成する請求項１に記載の心拍変異度分析装置。
【請求項５】
前記周波数変換部は、
前記低域通過フィルタを通過した前記シフティング検出信号のサンプリングレートを減少させるデシメータを更に含む請求項４に記載の心拍変異度分析装置。
【請求項６】
前記周波数変換部は、
前記低域通過フィルタを通過した前記シフティング検出信号にゼロデータを挿入するデータ挿入部を更に含む請求項４に記載の心拍変異度分析装置。
【請求項７】
前記周波数変換部に対象周波数帯域モードに対応する対象周波数帯域情報を提供する帯域選択部を更に含む請求項１に記載の心拍変異度分析装置。
【請求項８】
前記周波数帯域モードは、
高周波数帯域モード、低周波数帯域モード、及び超低周波数帯域モードを含み、
前記帯域選択部は、
前記高周波数帯域モードで前記高周波数帯域情報を前記周波数変換部に提供し、前記低周波数帯域モードで前記低周波数帯域情報を前記周波数変換部に提供し、前記超低周波数帯域モードで前記超低周波数帯域情報を前記周波数変換部に提供する請求項７に記載の心拍変異度分析装置。
【請求項９】
前記対象周波数帯域は高周波数帯域を含み、
前記高周波数帯域は０．１５乃至０．４Ｈｚの間の帯域を含む請求項８に記載の心拍変異度分析装置。
【請求項１０】
前記対象周波数帯域は低周波数帯域を含み、
前記低周波数帯域は０．０４乃至０．１５Ｈｚの間の帯域を含む請求項８に記載の心拍変異度分析装置。
【請求項１１】
前記対象周波数帯域は超低周波数帯域を含み、
前記超低周波数帯域は０．００３乃至０．０４Ｈｚの間の帯域を含む請求項８に記載の心拍変異度分析装置。
【請求項１２】
前記心拍変異度信号は、前記対象周波数帯域のパワースペクトル信号を含む請求項１に記載の心拍変異度分析装置。
【請求項１３】
前記対象周波数帯域は高周波数帯域及び低周波数帯域を含み、
前記心拍変異度分析部は、前記パワースペクトル信号に基づいて前記高周波帯域のパワー値、前記低周波帯域のパワー値、及び前記高周波数帯域に対する前記低周波数帯域のパワー比のうち少なくとも一つを算出する請求項１２に記載の心拍変異度分析装置。
【請求項１４】
心電図信号を受信するステップと、
前記心電図信号のＲ−ピーク値を検出するステップと、
前記Ｒ−ピーク値に基づいて生成されたＲＲ間隔時系列信号の対象周波数帯域を抽出して対象周波数信号を生成する周波数変換ステップと、
前記対象周波数信号に基づいて前記対象周波数帯域に対する心拍変異度信号を計算及び分析するステップと、を含み、
周波数変換ステップは、
前記ＲＲ間隔時系列信号の前記対象周波数帯域を基底帯域に移動させるステップと、
前記基底帯域に移動したＲＲ間隔時系列信号を低域通過フィルタでフィルタリングするステップと、
前記フィルタリングされたＲＲ間隔時系列信号を高速フーリエ変換ステップと、を含む心拍変異度検出方法。
【請求項１５】
前記周波数変換ステップは、
前記ＲＲ間隔時系列信号をＩ（Ｉｎ−ｐｈａｓｅ）チャネル信号とＱ（ｑｕａｄｒａｔｕｒｅ−ｐｈａｓｅ）チャネル信号に分流するステップを更に含む請求項１４に記載の心拍変異度検出方法。
【請求項１６】
前記周波数変換ステップは、
前記フィルタリングされたＲＲ間隔時系列信号のサンプリングレートを減少させるステップと、
前記フィルタリングされたＲＲ間隔時系列信号にゼロデータを挿入するステップと、を更に含む請求項１５に記載の心拍変異度検出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、心拍変異度の検出及び分析に関するものであって、より詳しくは、心拍変異度分析装置及びそれを利用した心拍変異度検出方法に関するものである。
【背景技術】
【０００２】
自律神経系は臓器の機能と物質の代謝を調節し、体内外の環境要因の変化に対して適切な均衡を図ることで生命維持活動及び体内恒常性の維持に関わる。自律神経系は様々な精神疾患、身体疾患、及びストレス性疾患に関わる。よって、自律神経系の状態を分析することで環境的なストレスを評価することができる。心拍変異度の分析は、交感神経及び副交感神経の状態を定量的に評価する。
【０００３】
心拍変異度の分析方法は、時間領域分析方法及び周波数領域分析方法を含む。時間領域分析方法は、ＱＲＳ波の時間間隔に基づいて心拍変異度を分析する方法である。例えば、時間領域分析方法はＳＤＮ、ＮＳＤＡＮ、ＮＲＭＳＳ、ＤＮＮ５０、ｐＮＮ５０などがある。時間領域分析方法は計算が容易であるが、交感神経系及び副交感神経系の活性度を区分するか、自律神経系の均衡性を定量化して示すことができない短所がある。
【０００４】
一方、周波数領域分析方法は心拍変異度の周波数帯域別の特性を検出し、該当帯域のパワー値を評価する方法である。周波数領域分析方法は、交感神経及び副交感神経に関わる自律神経系の均衡性を定量化して示すことができる。よって、周波数領域を分析するための心拍変異度分析装置及びその方法に対する要求が持続的に提起されている。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明は、心拍変異度の周波数領域を精密に分析することができる心拍変異度分析装置及びそれを利用した心拍変異度検出方法を提供する。
【課題を解決するための手段】
【０００６】
本発明の実施例による心拍変異度分析装置は、コンバータ、Ｒ−ピーク検出部、周波数変換部、及び心拍変異度分析部を含む。
【０００７】
コンバータは、ＱＲＳ波を含むアナログ心電図信号を受信してデジタル心電図信号に変換する。コンバータは、サンプリング周波数でサンプリングされたデジタル心電図信号を生成する。
【０００８】
Ｒ−ピーク検出部はデジタル心電図信号を受信してＲ−ピーク値を検出し、Ｒ−ピーク値に基づいてＲＲ間隔時系列信号を生成する。
【０００９】
周波数変換部は、ＲＲ間隔時系列信号から対象周波数帯域の信号を抽出して対象周波数信号を生成する。周波数変換部は周波数移動部、低域通過フィルタ、デシメータ、データ挿入部、及び高速フーリエ変換部を含む。
【００１０】
周波数移動部は、ＲＲ間隔時系列信号の対象周波数帯域を基底帯域に移動してシフティング検出信号を生成する。周波数移動部はＲＲ間隔時系列信号に第１正弦信号をかけてＩチャネル信号を生成し、ＲＲ間隔時系列信号に第２正弦信号をかけてＱチャネル信号を生成する。第１正弦信号及び第２正弦信号は９０度の位相差を有し、第１正弦信号及び第２正弦信号の周波数は対象周波数帯域の中心周波数に基づく。低域通過フィルタはシフティング検出信号を受信し、基底帯域のシフティング検出信号を通過させる。デシメータはシフティング検出信号のサンプリングレートを減少する。データ挿入部はシフティング検出信号にゼロデータを挿入する。
【００１１】
心拍変異度分析部は、対象周波数信号に基づいて対象周波数帯域に対する心拍変異度を生成する。心拍変異度信号は対象周波数帯域のパワースペクトル信号を含む。心拍変異度分析部はパワースペクトル信号に基づいて高周波数帯域のパワー値、低周波数帯域のパワー値、超低周波数帯域のパワー値、及び高周波数帯域に対する低周波数帯域のパワー比のうち少なくとも一つを算出する。
【００１２】
心拍変異度分析装置は、対象周波数帯域モードに対応する対象周波数情報を周波数変換部に提供する帯域選択部を更に含む。対象周波数帯域モードは高周波数帯域モード、低周波数帯域モード、及び超低周波数帯域モードを含む。
【００１３】
本発明の実施例による心拍変異度検出方法は、心電図信号を受信するステップと、Ｒ−ピーク値を検出するステップと、周波数変換ステップと、心拍変異度信号を計算及び分析するステップと、を含む。周波数変換ステップは、ＲＲ間隔時系列信号の対象周波数帯域を基底帯域に移動するステップと、低域通過フィルタでフィルタリングするステップと、ＲＲ間隔時系列信号のサンプリングレートを減少するステップと、ＲＲ間隔時系列信号にゼロデータを挿入するステップと、高速フーリエ変換ステップと、を含む。
【発明の効果】
【００１４】
本発明の実施例による心拍変異度分析装置及びそれを利用した心拍変異度検出方法は、対象周波数帯域の心拍変異度を検出して精密に分析することができる。
【図面の簡単な説明】
【００１５】
【図１】心電図信号の波形を示すグラフである。
【図２】心拍変異度信号の周波数パワースペクトルを示すグラフである。
【図３】本発明の一実施例による心拍変異度分析装置のブロック図である。
【図４】本発明の一実施例による周波数変換部のブロック図である。
【図５】図４の周波数変換部が対象周波数帯域の信号を検出する過程を示すグラフである。
【図６】図４の周波数変換部が対象周波数帯域の信号を検出する過程を示すグラフである。
【図７】図４の周波数変換部が対象周波数帯域の信号を検出する過程を示すグラフである。
【図８】図４の周波数変換部が対象周波数帯域の信号を検出する過程を示すグラフである。
【図９】本発明の他の実施例による心拍変異度分析装置のブロック図である。
【図１０】本発明の一実施例による心拍変異度検出方法を示す順序図である。
【図１１】本発明の一実施例による周波数を変換するステップを示す順序図である。
【発明を実施するための形態】
【００１６】
以下、本発明の技術分野における通常の知識を有する者が本発明を容易に実施し得るように、本発明の実施例を明確に詳細に記載する。
【００１７】
図１は、心電図信号の波形を示すグラフである。
【００１８】
心電図（Ｅｌｅｃｔｒｏｃａｒｄｉｏｇｒａｍ、ＥＣＧ）は、心臓の拍動によって心筋から発生する電気的活動をグラフで示すものである。図１を参照すると、心電図信号はＰ波、Ｑ波、Ｒ波、Ｓ波、及びＴ波を含む。図１は１回の心臓拍動によるＰ波、Ｑ波、Ｒ波、Ｓ波及びＴ波を示したが、このような波形は繰り返し形成される。
【００１９】
Ｐ波は心房の収縮に起因する。Ｐ波とＱ波との間隔は心房の電気的刺激が房室結節に伝導される時間に対応する。一般に、Ｐ波とＱ波との間隔は０．２秒以内に形成される。Ｑ波は心臓の電気的活動電流が減少する下向波である。Ｒ波はＱ波の後、心臓の電気的活動電流が増加する上向波である。Ｒ波によって心電図は最大値を有し、このような最大値はＲ−ピーク値と定義される。Ｓ波はＲ波の後、心臓の電気的活動電流が減少する下向波である。Ｑ波、Ｒ波、及びＳ波は心室筋の脱分極過程に起因する。Ｔ波は脱分極による心室の収縮の後、回復され弛緩される再分極過程に起因する。Ｑ波、Ｒ波、Ｓ波はＱＲＳ波と定義され、ＱＲＳ波は洞房結節から始まった電気刺激が房室結節を経て心室に伝導されて心室を収縮させる脱分極によって発生する。一般に、ＱＲＳ波は０．０６秒乃至０．１秒の間に形成される。
【００２０】
図１の心電図信号の波形は、時間に対する心臓の電気的活動を示すものである。本発明の心拍変異度分析装置は、心電図信号の波形のＲ−ピーク値を検出し、Ｒ−ピーク値の間の間隔を連続的に示す心拍数タコグラム（ｈｅａｒｔ  ｒａｔｅ  ｔａｃｈｏｇｒａｍ）による波形を周波数変換して、心拍変異度の周波数領域を分析する。具体的な心電図信号の周波数変換及び心拍変異度信号の生成過程は後述する。
【００２１】
図２は、心拍変異度信号の周波数パワースペクトルを示すグラフである。
【００２２】
心拍変異度（Ｈｅａｒｔ  ｒａｔｅ  ｖａｒｉａｂｉｌｔｙ、ＨＲＶ）は心拍数の拍動間変動を意味し、心拍の瞬間周期を測定して周期の変化量を検出する。心拍変異度は交感神経及び副交感神経の相互作用を反映する。
【００２３】
図２を参照すると、心拍変異度信号の波形として周波数領域に対するパワースペクトル密度（Ｐｏｗｅｒ  ｓｐｅｃｔｒａｌ  ｄｅｎｓｉｔｙ、ＰＳＤ）が図示されている。Ｒ−ピーク値間の間隔のように心拍動間の間隔が時間単位に測定され、このような時間間隔は周波数変換のための信号を生成する。パワースペクトルを計算する際、信号の単位は自乗される。心拍変異度の分析は、時間の測定から得られるＲ−ピーク値間の間隔に基づいてパワースペクトルを計算するため、縦軸の単位はｍｓｅｃ２単位と定義される。
【００２４】
心拍変異度信号の波形は高周波数帯域（Ｈｉｇｈ  ｆｒｅｑｕｅｎｃｙ  ｂａｎｄ、ＨＦ）、低周波数帯域（Ｌｏｗ  ｆｒｅｑｕｅｎｃｙ  ｂａｎｄ、ＬＦ）、及び超低周波数帯域（Ｖｅｒｙ  ｌｏｗ  ｆｒｅｑｕｅｎｃｙ  ｂａｎｄ、ＶＬＦ）に区分される。高周波数帯域、低周波数帯域、及び超低周波数帯域それぞれのパワー値は、交感神経系または副交感神経系の活動性を判断するために利用される。
【００２５】
高周波数帯域（ＨＦ）は０．１５乃至０．４Ｈｚの間の周波数帯域と定義される。高周波数帯域は副交感神経系の活動性に対する指標として活用される。高周波数帯域は呼吸周期に関わる。高周波数帯域は心臓の電気的安定度に関わる。よって、高周波数帯域のパワー値を分析することで心肺の活動性及び心臓の電気的安定度などを判断することができる。
【００２６】
低周波数帯域（ＬＦ）は０．０４乃至０．１５Ｈｚの間の周波数帯域と定義される。低周波数帯域は圧受容体の反射または血圧の調節などによる心拍数の変化に関わる。低周波数帯域は交感神経系の活動性に対する指標として活用される。例えば、偏頭痛患者のように交感神経系の活動性が過剰に示される場合、低周波数帯域の変動幅が大きく示される。但し、これに限らず、低周波数帯域ＬＦは交感神経系及び副交感神経系の活動性に対する指標として活用される。
【００２７】
超低周波数帯域（ＶＬＦ）は０．００３乃至０．０４Ｈｚの間の帯域と定義される。超低周波数帯域は交感神経系の活動性に対する指標として活用される。例えば、睡眠時無呼吸症候群の患者は無呼吸の間に交感神経が刺激されて超低周波数帯域のパワー値が増加する。また、超低周波数帯域は体温の調節、ホルモン、及び血管運動などに対する指標として活用される。
【００２８】
図３は、本発明の一実施例による心拍変異度分析装置のブロック図である。
【００２９】
図３を参照すると、心拍変異度分析装置１０００はコンバータ１００、Ｒピーク検出部２００、周波数変換部３００、心拍変異度分析部４００、帯域選択部５００、貯蔵部６００、及び表示部７００を含む。
【００３０】
コンバータ１００はアナログ信号を受信してデジタル信号に変換する。コンバータ１００はアナログ心電図信号（ＥＣＧＴ）を受信してデジタル心電図信号（ＥＣＧｎ）に変換する。コンバータ１００は心電図センサ（図示せず）と連結されるが、ユーザの心電図を感知した心電図センサはコンバータ１００にアナログ心電図信号を提供する。アナログ心電図信号は特定された時間の間にコンバータ１００を利用してデジタル心電図信号に変換される。デジタル心電図信号は離散信号である。
【００３１】
コンバータ１００はサンプリング信号（ｆｓ）を受信する。サンプリング信号は心電図信号のサンプリング周波数に関する情報を含む。サンプリング周波数値は２００乃至５００Ｈｚの間に形成されるが、これに限ることはない。コンバータ１００はサンプリング信号に基づいて変換されたデジタル心電図信号のデータ量を決定する。サンプリング信号の周波数値が大きければ高い安定性を有するデジタル心電図信号が生成され、サンプリング信号の周波数値が小さければ少ないデータを有するデジタル心電図信号が生成される。
【００３２】
Ｒ−ピーク検出部２００はコンバータ１００からデジタル心電図信号を受信する。Ｒ−ピーク検出部２００はデジタル心電図信号からＲ−ピーク値を検出する。Ｒ−ピーク検出部２００はＲ−ピーク値を検出するためのＱＲＳ波検出アルゴリズムを含む。例えば、Ｒ−ピーク検出部２００はＰａｎ＆ＴｏｍｐｋｉｎｓアルゴリズムまたはＨａｍｉｌｔｏｎ＆Ｔｏｍｐｋｉｎｓアルゴリズムを含む。Ｒ−ピーク検出部２００はＱＲＳ波検出アルゴリズムを介して測定時間の間に存在する複数のＲ波を検出し、各Ｒ波のピーク値であるＲ−ピーク値を順次に検出する。
【００３３】
Ｒ−ピーク検出部２００は順次に検出されたＲ−ピーク値に基づいてＲＲ間隔時系列信号（ＲＰｎ）を生成する。Ｒ−ピーク検出部２００は時間の経過によって順次に検出されたＲ−ピーク値間の間隔であるＲＲ間隔（Ｐ−ｐｅａｋ  ｔｏ  Ｒ−ｐｅａｋＩｎｔｅｒｖａｌ、ＲＲＩ）を計算する。Ｒ−ピーク検出部２００はＲＲ間隔（ＲＲＩ）を補間（ｉｎｔｅｒｐｏｌａｔｉｏｎ）してＲＲＩタコグラム（Ｔａｃｏｇｒａｍ）に対応するＲＲ間隔時系列信号を生成する。例えば、Ｒ−ピーク検出部２００はＲＲ間隔の連続的な情報に基づき、ＲＲ間隔の持続時間の標本密度分布または隣接したＲＲ間隔差の標本密度分布のような幾何学的パターンに変換するための補間を介して、ＲＲ間隔時系列信号を生成する。
【００３４】
周波数変換部３００はＲ−ピーク検出部２００からＲＲ感覚時系列信号を受信する。周波数変換部３００はＲＲ間隔時系列信号の対象周波数帯域（ｆｔ）を抽出して対象周波数信号（Ｔｆ）を生成する。周波数変換部３００は高速フーリエ変換（Ｆａｓｔ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ、ＦＦＴ）を利用してＲＲ間隔時系列信号を周波数系列の信号に変換する。
【００３５】
周波数変換部３００はＲＲ間隔時系列信号の周波数変換を行う前に、ＲＲ間隔時系列信号の対象周波数帯域に対応する信号値を抽出する。対象周波数帯域に対応する信号値を抽出するために、周波数変換部３００は後述する帯域選択部５００から周波数帯域情報（ｆｔｉ）を入力される。
【００３６】
対象周波数帯域とは、ＲＲ間隔時系列信号で分析の対象となる周波数帯域を意味する。例えば、分析の対象が高周波数帯域であれば、対象周波数帯域は０．１５乃至０．４Ｈｚの間の帯域を含む。分析の対象が低周波数帯域であれば、対象周波数帯域は０．０４乃至０．１５Ｈｚの間の帯域を含む。分析の対象が超低周波数帯域であれば、対象周波数帯域は０．００３乃至０．０４Ｈｚの間の帯域を含む。周波数変換部３００は対象周波数帯域情報に基づいて対象周波数信号を生成する。対象周波数信号は周波数領域に対する心拍変異度信号で対象周波数帯域に対応する信号である。具体的な周波数変換部３００の構成及び対象周波数信号の生成過程は後述する。
【００３７】
心拍変異度分析部４００は周波数変換部３００から対象周波数信号を受信する。心拍変異度分析部４００は対象周波数信号に基づいて心拍変異度を分析する。心拍変異度分析部４００は対象周波数信号に基づいて対象周波数帯域に対する心拍変異度信号（ＨＲＶ）を生成する。心拍変異度信号は対象周波数帯域のパワースペクトル信号（ＰＳＤ）を含む。心拍変異度信号は対象周波数信号に対する強度を計算し、パワー値などを算出して生成される。
【００３８】
例えば、心拍変異度信号は対象周波数帯域に応じて、高周波帯域のパワー値、低周波帯域のパワー値、または超低周波帯域のパワー値のうちいずれか一つである。上述したように、高周波帯域のパワー値、低周波帯域のパワー値、または超低周波帯域のパワー値は自律神経系の活動性に対する指標として利用される。
【００３９】
心拍変異度信号は高周波帯域に対する低周波帯域のパワー比（ＬＦ／ＨＦ）を含む。パワー比が高ければ交感神経が活性化されているか副交感神経の活性が抑制されていることを意味する。パワー比は交感神経系及び副交感神経系の間の均衡性を定量化して分析するようにする。心拍変異度信号はこれに限らず、極超低周波数帯域（Ｕｌｔｒａ−ｌｏｗ  ｆｒｅｑｕｅｎｃｙ  ｂａｎｄ）のパワー値など多様な分析信号を含む。このような多様な心拍変異度信号に基づいて、心拍変異度の分析が行われる。
【００４０】
帯域選択部５００は周波数変換部３００に対象周波数帯域情報を提供する。対象周波数帯域情報は対象周波数帯域または対象周波数帯域の中心周波数（ｆｃ)に関する情報を含む。帯域選択部５００は対象周波数帯域モードに対応する対象周波数帯域情報を周波数変換部３００に提供し、心拍変異度分析部４００が対象周波数帯域の心拍変異度信号を生成するようにする。
【００４１】
対象周波数帯域モードは複数の周波数帯域モードを含む。対象周波数帯域モードは高周波数帯域モード、低周波数帯域モード、及び超低周波数帯域モードを含み、更に多い帯域モードを含んでもよい。それぞれの対象周波数帯域モードにおける対象周波数帯域は互いに重畳されない。それぞれの対象周波数帯域モードにおいて、帯域選択部５００はそれぞれの対象周波数帯域に対応する対象周波数帯域情報を周波数変換部３００に提供する。
【００４２】
高周波数帯域モードにおいて、帯域選択部５００は高周波数帯域情報を周波数変換部３００に提供する。高周波数帯域情報を受信した周波数変換部３００は高周波数帯域に対する対象周波数信号を生成する。心拍変異度分析部４００は対象周波数信号に基づいて高周波数帯域に対する心拍変異度信号を生成する。
【００４３】
低周波数帯域モードにおいて、帯域選択部５００は低周波数帯域情報を周波数変換部３００に提供する。低周波数帯域情報を受信した周波数変換部３００は低周波数帯域に対する対象周波数信号を生成する。心拍変異度分析部４００は対象周波数信号に基づいて低周波数帯域に対する心拍変異度信号を生成する。
【００４４】
超低周波数帯域モードにおいて、帯域選択部５００は超低周波数帯域情報を周波数変換部３００に提供する。超低周波数帯域情報を受信した周波数変換部３００は超低周波数帯域に対する対象周波数信号を生成する。心拍変異度分析部４００は対象周波数信号に基づいて超低周波数帯域に対する心拍変異度信号を生成する。
【００４５】
対象周波数帯域モードは時間に応じて変更される。例えば、帯域選択部５００は一定時間の間に高周波数帯域モードに動作して高周波数帯域に関する対象周波数帯域情報を提供し、次に対象周波数帯域モードの変更によって低周波数帯域に対する対象周波数帯域情報を提供する。よって、帯域選択部５００のモード変更によって周波数変換部３００は多様な周波数帯域に対する対象周波数信号を生成する。また、これに基づいて心拍変異度分析部４００は多様な心拍変異度信号を生成し、心拍変異度に対する多角的な分析が可能になる。
【００４６】
貯蔵部６００は心拍変異度分析部４００に対する記憶装置として使用される。貯蔵部６００は心拍変異度分析部４００から心拍変異度信号、対象周波数信号、及び対象周波数帯域情報などを受信して貯蔵する。
【００４７】
貯蔵部６００はＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、ＰＲＯＭ（Ｐｒｏｇｒａｍｍａｂｌｅ  ＲＯＭ）、ＥＰＲＯＭ（Ｅｌｅｃｔｒｉｃａｌｌｙ  Ｐｒｏｇｒａｍｍａｂｌｅ  ＲＯＭ）、ＥＥＰＲＯＭ（Ｅｌｅｃｔｒｉｃａｌｌｙ  Ｅｒａｓａｂｌｅ  ａｎｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  ＲＯＭ）、フラッシュメモリ装置、ＰＲＡＭ（Ｐｈａｓｅ−ｃｈａｎｇｅ  ＲＡＭ）、ＭＲＡＭ（Ｍａｇｎｅｔｉｃ  ＲＡＭ）、ＲＲＡＭ（Ｒｅｓｉｓｔｉｖｅ  ＲＡＭ）またはＦＲＡＭ（Ｆｅｒｒｏｅｌｅｃｔｒｉｃ  ＲＡＭ）のような不揮発性メモリ装置のうち少なくとも一つを含み、ＳＲＡＭ（Ｓｔａｔｉｃ  ＲＡＭ）、ＤＲＡＭ（Ｄｙｎａｍｉｃ  ＲＡＭ）、またはＳＤＲＡＭ（Ｓｙｎｃｈｒｏｎｏｕｓ  ＤＲＡＭ）のような揮発性メモリ装置のうち少なくとも一つを含む。
【００４８】
表示部７００は心拍変異度分析部４００によって生成された心拍変異度信号、周波数変換部３００によって生成された対象周波数信号、または貯蔵部６００に貯蔵されたデータに基づいたイメージ信号（ＩＨＲＶ）を受信して表示する。
【００４９】
表示部７００はＬＣＤ（Ｌｉｑｕｉｄ  Ｃｒｙｓｔａｌ  Ｄｉｓｐｌａｙ）、ＯＬＥＤ（Ｏｒｇａｎｉｃ  Ｌｉｇｈｔ  Ｅｍｉｔｔｉｎｇ  Ｄｉｏｄｅ）、ＡＭＯＬＥＤ（Ａｃｔｉｖｅ  Ｍａｔｒｉｘ  ＯＬＥＤ）、フレキシブルディスプレイ、または電子インクのうち少なくとも一つを含む。
【００５０】
図４は、本発明の一実施例による周波数変換部３００に対するブロック図である。図５乃至図８は、図４の周波数変換部３００が対象周波数帯域の信号を検出する過程を説明するためのグラフである。
【００５１】
図４を参照すると、周波数変換部３００は周波数移動部３１１、３１２、低域通過フィルタ３１２、３２２、デシメータ３３１、３３２、データ挿入部３４１、３４２、及び高速フーリエ変換部３５０を含む。但し、図４とは異なって、周波数変換部３００は前記構成要素のうち一部を含まなくてもよい。例えば、周波数変換部３００はデシメータ３３１、３３２を含まないか、データ挿入部３４１、３４２を含まなくてもよい。
【００５２】
周波数移動部３１１、３１２は第１周波数移動部３１１及び第２周波数移動部３１２を含む。第１周波数移動部３１１及び第２周波数移動部３１２はＲ−ピーク検出部２００からＲＲ間隔時系列信号を受信する。
【００５３】
第１周波数移動部３１１はＲＲ間隔時系列信号を受信してＩ（Ｉｎ−ｐｈａｓｅ）チャネル信号を生成する。第１周波数移動部３１１はＩチャネル信号を生成するために、第１正弦信号（Ｉｆｔ）をＲＲ間隔時系列信号にかける。ＲＲ間隔時系列信号がｎに対する離散信号であれば、第１正弦信号はｃｏｓ（２＊ｐｉ＊ｆｃ＊ｎ＊ｔｓ）と定義される。第１正弦信号のｆｃは分析しようとする対象周波数帯域の中心周波数と定義される。第１正弦信号のｔｓはサンプリング周波数の逆数に対応する。
【００５４】
第２周波数移動部３１２はＲＲ間隔時系列信号を受信してＱ（ｑｕａｄｒａｔｕｒｅ−ｐｈａｓｅ）チャネル信号を生成する。第２周波数移動部３１２はＱチャネル信号を生成するために、第２正弦信号（Ｑｆｔ）をＲＲ間隔時系列信号にかける。ＲＲ間隔時系列信号がｎに対する離散信号であれば、第２正弦信号は−ｓｉｎ（２＊ｐｉ＊ｆｃ＊ｎ＊ｔｓ）と定義される。第１周波数移動部３１１及び第２周波数移動部３１２はＲＲ間隔時系列信号をＩチャネル信号及びＱチャネル信号に分流する。第１周波数移動部３１１はＲＲ間隔時系列信号の実数部を抽出し、第２周波数移動部３１２はＲＲ間隔時系列信号の虚数部を抽出する。
【００５５】
周波数変換部３００は第１正弦信号及び第２正弦信号を生成してＲＲ時系列信号を分流するための分流器（図示せず）を更に含む。分流器は第１正弦信号及び第２正弦信号を生成する。第１正弦信号及び第２正弦信号の周波数は対象周波数帯域及びサンプリング周波数に基づくため、分流器は帯域選択部５００から対象周波数帯域情報を受信し、サンプリング信号を受信する。これとは異なって、分流器は帯域選択部５００に提供されて第１周波数移動部３１１及び第２周波数移動部３１２に第１正弦信号及び第２正弦信号を提供する。または、分流器は心拍変異度分析装置１０００に別途の構成要素として備えられてもよい。
【００５６】
第１正弦信号または第２正弦信号がＲＲ間隔時系列信号にかけられると、周波数領域においてＲＲ間隔時系列信号は中心周波数だけ移動する。つまり、対象周波数帯域は基底帯域（ｂａｓｅｂａｎｄ）に移動する。第１周波数移動部３１１及び第２周波数移動部３１２はＲＲ間隔時系列信号をＩチャネルとＱチャネル信号に分流し、対象周波数帯域を基底帯域にシフティングするためにそれぞれ第１正弦信号と第２正弦信号をかけて、シフティング検出信号を生成する。
【００５７】
第１低域通過フィルタ３２１及び第２低域通過フィルタ３２２はシフティング検出信号を受信する。第１低域通過フィルタ３２１はＩチャネル信号を受信し、第２低域通過フィルタ３２２はＱチャネル信号を受信する。第１低域通過フィルタ３２１及び第２低域通過フィルタ３２２は基底帯域のシフティング検出信号を通過し、高い周波数帯域のシフティング検出信号を遮断する。よって、ＲＲ間隔時系列信号から対象周波数帯域に対応する信号を選択的に抽出することができる。また、分析の対象となる周波数帯域の信号を抽出し、残りの帯域の信号を遮断するため、不必要なデータを除去してデータの量を減らすことができる。
【００５８】
第１低域通過フィルタ３２１及び第２低域通過フィルタ３２２は、帯域選択部５００の対象周波数帯域モードに応じて遮断周波数を変換する。例えば、高周波数帯域モードにおいて高周波数の帯域幅は約０．２５Ｈｚであり、低周波数帯域モードにおいて低周波数の帯域幅は約０．１１Ｈｚであって、超低周波数帯域モードにおいて超低周波数の帯域幅は約０．０４Ｈｚである。よって、高周波数帯域モードでは約０．１２５Ｈｚ程度の遮断周波数が要求され、低周波数帯域モードでは約０．０５５Ｈｚ程度の遮断周波数が要求され、超低周波数の帯域幅は約０．０２Ｈｚの遮断周波数が要求される。よって、第１低域通過フィルタ３２１及び第２低域通過フィルタ３２２は帯域選択部５００から対象周波数帯域情報を受信して遮断周波数を変換する。対象周波数帯域情報は該当対象周波数帯域モードに対する帯域幅情報を含む。
【００５９】
第１デシメータ３３１は第１低域通過フィルタ３２１を通過したシフティング検出信号を受信し、第２デシメータ３３２は第２低域通過フィルタ３２２を通過したシフティング検出信号を受信する。第１デシメータ３３１及び第２デシメータ３３２は信号のサンプリングレート（ｓａｍｐｌｉｎｇ  ｒａｔｅ）を減少させる。例えば、第１デシメータ３３１及び第２デシメータ３３２が４次デシメーションを行う場合、従来のサンプリング周波数の１／４だけサンプリング周波数が減少する。この場合、シフティング検出信号のデータは１／４だけ減少する。また、基底帯域のシフティング検出信号は、周波数領域においてサンプリング周波数を基準に４倍拡張された波形で示される。第１デシメータ３３１及び第２デシメータ３３２がＭ次デシメーションを行う場合、Ｍ値に基づいて周波数分解能が決定される。Ｍ次デシメーションを行う際、周波数分解能は（サンプリング周波数／Ｍ／（ＦＦＴポイントの個数））と定義される。
【００６０】
第１デシメータ３３１及び第２デシメータ３３２は倍率信号（ＭＭ）を受信してデシメーションの次数を決定する。図示していないが、周波数変換部３００は倍率信号を生成して第１デシメータ３３１及び第２デシメータ３３２に倍率信号を提供する倍率器を更に含む。または、倍率器は心拍変異度分析装置１０００に別途の構成要素として備えられてもよい。
【００６１】
第１データ挿入部３４１及び第２データ挿入部３４２はデシメーションが行われた信号に０のデータを挿入（ｚｅｒｏ  ｐａｄｄｉｎｇ）する。０のデータはシフティング検出信号それぞれのデータの間に挿入される。第１データ挿入部３４１及び第２データ挿入部３４２は、第１及び第２デシメータ３３１、３３２から出力される信号が直接高速フーリエ変換部３５０に入力される場合のデータの数が基準ＦＦＴポイントの個数より小さければ、そのデータの個数の差だけ０のデータを挿入することで、周波数分解能を改善し、信号の精密度を向上させる。
【００６２】
第１デシメータ３３１及び第２デシメータ３３２はパディング信号（ＰＰ）を受信してデータの挿入個数を決定する。図示していないが、周波数変換部３００はパディング信号を生成して第１データ挿入部３４１及び第２データ挿入部３４２にパディング信号を提供する分解能調整器を更に含む。または、分解能調整器は心拍変異度分析装置１０００に別途の構成要素として備えられてもよい。
【００６３】
第１デシメータ３３１及び第２デシメータ３３２は対象周波数帯域に対するＲＲ間隔時系列信号を集中して観察するようにし、データの量を減らす。第１データ挿入部３４１及び第２データ挿入部３４２は対象周波数帯域に対するＲＲ間隔時系列信号の周波数分解能を改善し、精密度を向上するためにデータを挿入する。よって、データの効率性及び処理速度の観点で第１データ挿入部３４１及び第２データ挿入部３４２は除去される。
【００６４】
高速フーリエ変換部３５０は第１データ挿入部３４１及び第２データ挿入部３４２から０のデータ挿入された信号を受信する。高速フーリエ変換部３５０は第１及び第２低域通過フィルタ３２１、３２２、第１及び第２デシメータ３３１、３３２、第１及び第２データ挿入部３４１、３４２を介して伝達されたＲＲ間隔時系列信号を高速フーリエ変換する。高速フーリエ変換部３５０は離散フーリエ変換（Ｄｉｓｃｒｅｔｅ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）の複雑な計算を減らす高速フーリエ変換アルゴリズムを含み、高速フーリエ変換アルゴリズムに基づいて対象周波数信号を生成する。
【００６５】
図５は、周波数変換部３００に入力される信号を例示的に示している。図５の波形は周波数変換部３００に入力される信号を周波数変換して示したグラフである。ユーザは対象周波数帯域を観測するために周波数変換部３００に対象周波数帯域情報を入力する。
【００６６】
図６は、図５の信号が第１周波数移動部５１１及び第２周波数移動部５１２によって変換された信号の大きさを示している。第１周波数移動部５１１または第２周波数移動部５１２は図５の信号に第１正弦信号または第２正弦信号をかけて対象周波数帯域の信号を基底帯域を移動する。周波数領域において、図５ａの信号は中心周波数だけ負の方向に移動する。
【００６７】
図７は、図６の信号が第１及び第２低域通過フィルタ３２１、３２２によって変換された信号の大きさを示している。第１及び第２低域通過フィルタ３２１、３２２は対象周波数帯域／２の遮断周波数値を有する。第１及び第２低域通過フィルタ３２１、３２２によって、図６の信号は基底帯域の波形だけ残り、残りの波形は除去される。
【００６８】
図８は、図７の信号が第１及び第２デシメータ３３１、３３２によって変換された信号の大きさを示している。第１及び第２デシメータ３３１、３３２はＭ次デシメーションを行う。デシメーションによってサンプリング周波数は１／Ｍだけ減少する。よって、第１及び第２デシメータ３３１、３３２は対象周波数帯域の信号を図８の信号のように集中して観察するようにする。
【００６９】
図９は、本発明の他の実施例による心拍変異度分析装置のブロック図である。
【００７０】
図９を参照すると、心拍変異度分析装置２０００はコンバータ２１００、Ｒ−ピーク検出部２２００、第１周波数変換部２３００ａ、第２周波数変換部２３００ｂ、第３周波数変換部２３００ｃ、及び心拍変異度分析部２４００を含む。図９のコンバータ２１００、Ｒ−ピーク検出部２２００、及び心拍変異度分析部２４００の構成はそれぞれ図３のコンバータ１００、Ｒ−ピーク検出部２００、及び心拍変異度分析部４００に対応するため、具体的な説明は省略する。
【００７１】
第１周波数変換部２３００ａはＲーピーク検出部２２００からＲＲ間隔時系列信号を受信して第１対象周波数信号（Ｔｆ１）を出力する。第１周波数変換部２３００ａはＲＲ間隔時系列信号のうち第１対象周波数帯域（ｆｔ１）に対する信号を抽出する。例えば、第１対象周波数帯域は高周波数帯域である。
【００７２】
第２周波数変換部２３００ｂはＲーピーク検出部２２００からＲＲ間隔時系列信号を受信して第２対象周波数信号（Ｔｆ２）を出力する。第２周波数変換部２３００ｂはＲＲ間隔時系列信号のうち第２対象周波数帯域（ｆｔ２）に対する信号を抽出する。例えば、第２対象周波数帯域は低周波数帯域である。
【００７３】
第３周波数変換部２３００ｃはＲーピーク検出部２２００からＲＲ間隔時系列信号を受信して第３対象周波数信号（Ｔｆ３）を出力する。第３周波数変換部２３００ｃはＲＲ間隔時系列信号のうち第３対象周波数帯域（ｆｔ３）に対する信号を抽出する。例えば、第３対象周波数帯域は超低周波数帯域である。
【００７４】
図３の心拍変異度分析装置１０００に比べ、図９の心拍変異度分析装置２０００は複数個の周波数変換部を含む。よって、多様な周波数帯域の心拍変異度信号を同時に分析することができる。例えば、心拍変異度分析装置２０００は高周波数帯域、低周波数帯域、及び超低周波数帯域の心拍変異度を同時に算出し分析する。図９の心拍変異度分析装置２０００に比べ、図３の心拍変異度分析装置１０００は相対的に少ないデータ容量で心拍変異度を分析する。
【００７５】
図１０は、本発明の他の実施例による心拍変異度検出方法Ｓ１０００を示す順序図である。
【００７６】
図１０を参照すると、心拍変異度検出方法Ｓ１０００は心電図信号を受信するステップＳ１００と、Ｒ−ピーク値を検出するステップＳ２００と、周波数を変換するステップＳ３００と、心拍変異度信号を計算及び分析するステップＳ４００と、を含む。
【００７７】
心電図信号を受信するステップＳ１００において、心拍変異度検出装置は外部からアナログ心電図信号を受信する。心電図信号を受信するステップＳ１００はアナログ心電図信号をデジタル心電図信号に変換するステップを含む。心電図信号を受信するステップＳ１００は、アナログ−デジタルコンバータ１００、２１００によって行われる。この場合、アナログ心電図信号はアナログ−デジタルコンバータ１００、２１００が受信したサンプリング信号に含まれたサンプリング周波数に基づいてデジタル心電図信号に変換される。
【００７８】
Ｒ−ピーク値を検出するステップＳ２００において、心拍変異度検出装置は心電図信号のＲ−ピーク値を検出する。Ｒ−ピーク値を検出するステップＳ２００はＲ−ピーク検出部２００、２２００で行われる。Ｒ−ピーク値を検出するステップＳ２００は、Ｒ−ピーク値に基づいてＲＲ間隔時系列信号を生成するステップを含む。ＲＲ間隔時系列信号は順次に検出されたＲ−ピーク値に基づいてＲＲ間隔を計算して生成される。
【００７９】
周波数を変換数ステップＳ３００において、心拍変異度検出装置はＲＲ時系列信号の対象周波数帯域を抽出して対象周波数信号に変換する。周波数変換するステップＳ３００は、図３の周波数変換部３００または図９の第１乃至第３周波数変換部２３００ａ、２３００ｂ、２３００ｃで行われる。具体的な内容は後述する。
【００８０】
心拍変異度信号を計算及び分析するステップＳ４００において、心拍変異度検出装置は対象周波数信号に基づいて対象周波数帯域に対する心拍変異度信号を生成する。心拍変異度信号を計算及び分析するステップＳ４００は、心拍変異度分析部４００、２４００で行われる。心拍変異度信号を計算及び分析するステップＳ４００において、心拍変異度分析部４００、２４００は高周波数帯域のパワー値、低周波数帯域のパワー値、または超低周波数帯域のパワー値、パワー比などを計算し分析して自律神経系などの活動性を判断し、精神健康などを診断する。
【００８１】
図１１は、図１０の周波数を変換するステップＳ３００を示す順序図である。
【００８２】
図１１を参照すると、周波数を変換するステップＳ３００はＱチャネルとＩチャネルを分流するステップＳ３１０と、基底帯域に移動させるステップＳ３２０と、低域通過フィルタでフィルタリングするステップＳ３３０と、デシメーションするステップＳ３４０と、ゼロパディング（ｚｅｒｏ  ｐａｄｄｉｎｇ）するステップＳ３５０と、高速フーリエ変換するステップＳ３６０と、を含む。図８の各ステップは、図３の周波数変換部３００または図９の第１乃至図３周波数変換部２３００ａ、２３００ｂ、２３００ｃで行われる。
【００８３】
ＱチャネルとＩチャネルを分流するステップＳ３１０において、心拍変異度分析装置はＲＲ間隔時系列信号を実数部に対応するＩチャネルと虚数部に対応するＱチャネルに分流する。ＱチャネルとＩチャネルを分流するステップＳ３１０は、第１正弦信号をＲＲ間隔時系列信号にかけてＩチャネルを生成するステップ、及び第２正弦信号をＲＲ間隔時系列信号にかけてＱチャネルを生成するステップを含む。Ｉチャネルを生成するステップは第１周波数の移動部３１１で行われ、Ｑチャネルを生成するステップは第２周波数の移動部３１２によって行われる。
【００８４】
基底帯域に移動させるステップＳ３２０において、心拍変異度分析装置はＲＲ間隔時系列信号の対象周波数帯域を基底帯域に移動させてシフティング検出信号を生成する。基底帯域に移動させるステップＳ３２０は第１周波数移動部３１１及び第２周波数移動部３１２によって行われる。基底帯域に移動させるステップＳ３２０において、ＲＲ間隔時系列信号は対象周波数帯域の中心周波数だけ移動する。基底帯域に移動するために第１正弦信号及び第２正弦信号は正弦波関数であり、正弦波関数の周波数は中心周波数に比例する。ＱチャネルとＩチャネルを分流するステップＳ３１０と基底帯域に移動させるステップＳ３２０は、第１周波数移動部３１１及び第２周波数移動部３１２で同時に行われる。
【００８５】
低域通過フィルタでフィルタリングするステップＳ３３０において、心拍変異度分析装置は基底帯域のシフティング検出信号を通過させ、残りの周波数帯域のシフティング検出信号は除去する。低域通過フィルタでフィルタリングするステップＳ３３０は、第１及び第２低域通過フィルタ３１２、３２２で行われる。
【００８６】
デシメーションするステップＳ３４０において、心拍変異度分析装置はフィルタリングされたシフティング検出信号のサンプリングレートを減少させる。デシメーションするステップＳ３４０は第１及び第２デシメータ３３１、３３２で行われる。デシメーションするステップＳ３４０において、データの量は減り、サンプリング周波数はデシメーション次数だけ減少する。
【００８７】
ゼロパディングするステップＳ３５０において、心拍変異度分析装置はフィルタリングされたシフティング検出信号にゼロデータを挿入する。ゼロパディングするステップＳ３５０は第１及び第２データ挿入部３４１、３４２で行われる。ゼロパディングするステップＳ３５０において、デシメーションするステップＳ３４０によってデータが減少した際、基準ＦＦＴポイントの個数より高速フーリエ変換部３５０に入力されるデータの個数が足りなければ、第１及び第２データ挿入部３４１、３４２はそのデータの個数の差だけゼロデータを挿入する。ゼロパディングするステップ３５０において、データの量が増え、周波数分解能が改善される。
【００８８】
高速フーリエ変換するステップＳ３６０において、心拍変異度分析装置はフィルタリング、デシメーション、及びゼロパディングされたシフティング検出信号を時間領域から周波数領域に高速フーリエ変換する。高速フーリエ変換するステップＳ３６０は高速フーリエ変換部３５０で行われる。
【００８９】
上述した内容は本発明を実施するための具体的な例である。本発明には上述した実施例だけでなく、単純な設計変更や容易に変更し得る実施例も含まれるはずである。また、本発明には上述した実施例を利用してこれから容易に変形して実施し得る技術も含まれるはずである。
【符号の説明】
【００９０】
１０００、２０００：心拍変異度分析装置
１００、２１００：コンバータ
２００、２２００：Ｒ−ピーク検出部
３００、２３００ａ、２３００ｂ、２３００ｃ：周波数変換部
４００、２４００：心拍変異度分析部
５００：帯域選択部
ｆｓ：サンプリング信号
ｆｔ：対象周波数帯域
ｆｃ：中心周波数
ＥＣＧＴ：アナログ心電図信号
ＥＣＧｎ：デジタル心電図信号
ＲＲｎ：ＲＲ間隔時系列信号
Ｔｆ：対象周波数信号
ＨＲＶ：心拍変異度信号
【図１】
2017225808.tif 000003
【図２】
2017225808.tif 000004
【図３】
2017225808.tif 000005
【図４】
2017225808.tif 000006
【図５】
2017225808.tif 000007
【図６】
2017225808.tif 000008
【図７】
2017225808.tif 000009
【図８】
2017225808.tif 000010
【図９】
2017225808.tif 000011
【図１０】
2017225808.tif 000012
【図１１】
2017225808.tif 000013