JP4881574B2

JP4881574B2 - 経穴位置評価装置

Info

Publication number: JP4881574B2
Application number: JP2005124291A
Authority: JP
Inventors: 剛智福元; 恒秋山
Original assignee: Panasonic Corp; Matsushita Electric Industrial Co Ltd
Current assignee: Panasonic Corp; Panasonic Holdings Corp
Priority date: 2005-04-21
Filing date: 2005-04-21
Publication date: 2012-02-22
Anticipated expiration: 2025-04-21
Also published as: JP2006296809A; WO2006115072A1; EP1872718A4; US20090036798A1; CN101179988A; CN100581452C; US7818054B2; EP1872718A1

Description

本発明は、鍼治療などにおいて、鍼刺激の効果と経穴（ツボ）と呼ばれる刺激点の位置を非侵襲的かつ客観的に、しかも高精度に判定できる経穴位置評価装置に関する。

古くから皮膚表面上に周辺より皮膚抵抗が低い部位が小さな領域で存在し、その分布は鍼灸治療における刺激点である経穴の分布と比較的よく一致することが知られている。このような皮膚抵抗が周辺より低い部位（皮膚抵抗減弱点と言う）は、良導絡自律神経調整法や低周波治療器などにおいて治療に利用されており、これらの皮膚抵抗減弱点を探索する装置も市販され、実際に臨床での診断や治療に用いられている。

従来、これらの皮膚抵抗減弱点である経穴を探索（判定）する装置としては、経穴周辺の皮膚の電気的特性として皮膚の直流抵抗を測定し、その違いから経穴を判定する方法が多用されてきたが、これら直流抵抗の測定による経穴の判定では、測定結果が測定電極の接触状態に依存したり、皮膚電圧の自発的な変動との区別が明確にできないなどの問題があり、測定結果の信頼性・再現性に乏しい。

そこで、近年においては、これら測定結果の信頼性・再現性を向上させる手法として、経穴周辺の皮膚の電気的特性として、経穴周辺の皮膚への交流電流の印加により生成する電圧に基づいて皮膚インピーダンスの周波数特性の複素平面上での分布（皮膚インピーダンス軌跡）を測定し、皮膚インピーダンス軌跡の違いに基づいて経穴の判定を行うものがある。例えば皮膚インピーダンス軌跡の形状を特徴付ける４つのパラメータ（後述するZ₀、Z_∞、β、τ_m）のうち、測定電極の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映するパラメータとして Cole-Cole の円弧則における中心緩和時間τ_mに着目し、そのパラメータτ_mの違いに基づいて経穴の判定を行っている。このようにすることで、従来技術が有する問題点を回避できるようになっている（例えば下記の特許文献１参照）。
特開２００４−３３７３４９号公報（要約書、請求項２）

しかしながら、上記従来技術では、経穴とそうでない部位での電気的特性の差を評価する際に、皮膚インピーダンス軌跡の形状を特徴付けるパラメータとして、Cole-Cole の円弧則における中心緩和時間τ_mを算出しているが、以下のような問題点を有していた。
皮膚インピーダンス軌跡の形状を特徴付ける中心緩和時間τ_mを算出する際、皮膚インピーダンスの周波数特性において Cole-Cole の円弧則が成立している必要があるが、その成立の可否を確認する手段を有していない。また、測定された複数の周波数における皮膚インピーダンスから中心緩和時間τ_mを非線形最小２乗法などによって推定する必要があり、推定するための手順が煩雑で、これを実現するハードウエアの構成が複雑になり、装置の小型化には不向きであった。さらに、中心緩和時間τ_mを精度良く推定するには、図４に示すように、皮膚インピーダンス軌跡上に測定された皮膚インピーダンスデータが偏らず、まんべんなく分布している必要がある。

しかし、実際の皮膚で測定される皮膚インピーダンスデータの分布は、図５に示すように、インピーダンス軌跡上の一部分に偏って存在することが多く、又、その分布の仕方は測定点ごとに大きく異なっている場合が多い。このような場合、皮膚インピーダンス軌跡の形状を特徴付ける中心緩和時間τ_mを精度良く推定するには、極めて低い周波数も含めて広範囲の周波数範囲の皮膚インピーダンスを測定する必要があった。そのため、時間分解能が低下するほか、測定対象周波数が低くなればなるほど温度ドリフトや体動などによる外乱の影響を受けやすくなり、測定結果の信頼性・再現性を保証するのが困難であった。

本発明は上記従来技術における問題点を解決するもので、簡単なハードウエア構成で、信頼性・再現性のある経穴の判定を行うことができる経穴位置評価装置を提供することを目的とする。

上記従来の問題点を解決するために、本発明に係る経穴位置評価装置は、
少なくとも３つの周波数成分から構成される電流を発生する電流発生部と、
前記電流発生部からの出力電流を生体の皮膚に配置した複数の電極により複数の異なる皮膚上の測定点にほぼ同時に通電する電極系と、
前記複数の測定点にそれぞれ通電された電流を検出する電流検出器と、
前記電流検出器により検出された電流と、前記電極系の通電によって前記複数の測定点の皮膚に生じた電圧を測定する測定部と、
前記測定部により測定された電流と電圧の両方の周波数解析を行い、各測定点での皮膚インピーダンスの周波数特性を得て、前記各測定点での皮膚インピーダンスの周波数特性において Cole-Cole の円弧則が成立するか否かを確認する周波数解析部と、
前記周波数解析部により解析された前記複数の測定点での皮膚インピーダンスの周波数特性から、前記各測定点での皮膚インピーダンスの特徴量を生成する特徴量生成部と、
前記特徴量生成部により生成された前記各測定点での特徴量から経穴位置を判定する判定部と、
前記特徴量生成部により生成された前記各測定点での特徴量と前記判定部により判定された経穴位置とを表示する表示手段とを、
備え、前記周波数解析部により解析された前記複数の測定点での皮膚インピーダンスが測定電流密度に依存することなく前記Cole-Coleの円弧則が成立した場合に、前記特徴量生成部において生成される特徴量として通電電流に含まれる最低周波数での皮膚インピーダンスの実数部と虚数部との比を算出し、前記判定部で経穴位置を判定するよう構成されている。
この構成により、評価に用いる特徴量の算出手順が簡単になり、測定結果の信頼性・再現性を従来と比べ高くすることができるほか、評価に用いる特徴量は従来着目していた中心緩和時間τ_ｍと相関があるので、測定電極の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映する評価パラメータの違いに基づいて経穴の判定を行うことが可能になる。また、この構成により、測定された皮膚インピーダンスにおいて Cole-Cole の円弧則の成立可否を確認でき、測定結果の信頼性・再現性を従来と比べ高くすることができる。さらに、この構成により、パラメータの推定手順の煩雑さを回避することができるほか、経穴位置の判定を迅速に行える。

また、本発明に係る経穴位置評価装置は、
前記電流発生部において発生させる通電電流には、１つの最低周波数成分と、その整数倍の周波数成分が含まれることを特徴とする。
この構成により、離散フーリエ変換をした場合の隣接する周波数へのスペクトルの洩れは無く、時間分解能の下限は最低周波数のみに依存して決まるので、高時間分解能で皮膚インピーダンス軌跡を測定することができる。

また、本発明に係る経穴位置評価装置は、
前記電流発生部が、通電電流に含ませる周波数成分の大きさを、周波数ごと並びに測定位置ごとに変更可能とすることを特徴とする。
この構成により、各周波数、各測定点ごとに最適な信号対雑音比を確保でき、高い精度で皮膚インピーダンス軌跡を測定することができる。

また、本発明に係る経穴位置評価装置は、
前記特徴量生成部において、少なくとも１つの周波数における皮膚インピーダンスの実数部と虚数部との比のなかで、前記各測定点で差が顕著となる周波数での比を特徴量とする。
この構成により、パラメータの推定手順の煩雑さを回避することができるほか、判定部での経穴位置の判定の信頼性がより高くなる。

本発明によれば、測定された皮膚インピーダンスにおいて Cole-Cole の円弧則が成立するか否かを確認でき、評価に用いる特徴量の算出手順を簡単なハードウエア構成で実現できる。また、測定結果の信頼性・再現性を従来と比べ高くすることができるほか、評価に用いる特徴量は従来着目していた中心緩和時間τ_mと相関があるので、測定電極の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映する評価パラメータの違いに基づいて経穴の判定を行うことが可能になる。したがって、簡単なハードウエア構成でより信頼性・再現性の高い経穴位置装置を提供することができる。

以下、図面を参照して本発明の実施の形態について説明する。図１は本発明の経穴位置評価装置の第１の実施の形態の概略構成を示すブロック図である。電流発生部１から発生された交流電流Ｉ₁〜Ｉ_nはそれぞれ電流検出器２₁〜２_n、電流印加電極３₁〜３_nを介して、被検体の皮膚３０の各測定点＃１〜＃ｎへ印加され、接地電極４へと流れていく。この通電によって生じた電流印加電極３₁〜３_nとマイナス電極５（−）との間の皮膚に生じた電圧Ｖ₁〜Ｖ_nの降下を接地電極４の電位を基準にして測定部６により測定する。このような電極系により測定する手法は３電極法と呼ばれ、電流印加電極３₁〜３_nの直下、つまり、測定点＃１〜＃ｎの直下における皮膚インピーダンスを測定する際によく用いられる。

図２に電流発生部１並びに測定部６のより詳細な構成例を示すブロック図を示す。電流発生部１は少なくとも１つの電流源１₁〜１_nから構成されており、電流源１₁〜１_nはそれぞれ各測定点＃１〜＃ｎへ通電する交流電流Ｉ₁〜Ｉ_nを発生する。ここで、各電流源１₁〜１_nから発生される電流Ｉ₁〜Ｉ_nに含ませる周波数成分の数、周波数、大きさは、制御部２０からの制御信号２１０により設定可能になっている。測定部６は、通電によって各測定点＃１〜＃ｎ直下の皮膚に生じた電圧Ｖ₁〜Ｖ_nを測定するための少なくとも１つの差動増幅器６１₁〜６１_nと、少なくとも１つのハイパスフィルタ（ＨＰＦ）６２₁〜６２_n（並びに６７₁〜６７_n）と、少なくとも１つのプログラマブルゲインアンプ（ＰＧＡ）６３₁〜６３_n（並びに６８₁〜６８_n）と、少なくとも１つのローパスフィルタ（ＬＰＦ）６４₁〜６４_n（並びに６９₁〜６９_n）と、少なくとも１つのＡ／Ｄ変換器６５₁〜６５_n（並びに７０₁〜７０_n）とから構成されている。

通電によって測定点＃１〜＃ｎの皮膚に生じた電圧Ｖ₁〜Ｖ_nはそれぞれ差動増幅器６１₁〜６１_nにより測定されるが、測定された電圧Ｖ₁〜Ｖ_nには通電に起因しない不要な信号成分が含まれている。これは周波数解析部７での正確な周波数解析の妨げになるため、ハイパスフィルタ６２₁〜６２_nによりそれぞれ差動増幅器６１₁〜６１_nの出力から不要な低周波成分を除去する。不要な低周波成分が除去された測定点＃１〜＃ｎの電圧Ｖ₁〜Ｖ_nはそれぞれ、必要に応じてプログラマブルゲインアンプ６３₁〜６３_nにより増幅され、次いでローパスフィルタ６４₁〜６４_nにより不要な高周波成分が除去される。

また、測定点＃１〜＃ｎの皮膚に通電された電流Ｉ₁〜Ｉ_nは、それぞれ電流検出器２₁〜２_nにより測定されるが、測定点＃１〜＃ｎでの電圧Ｖ₁〜Ｖ_nと同じ信号処理を、測定点＃１〜＃ｎに通電した電流Ｉ₁〜Ｉ_nに対して行なわなければ、図１に示す周波数解析部７での周波数解析の手順が繁雑になり、解析結果も不正確になってしまう。そのため、測定点＃１〜＃ｎの電圧Ｖ₁〜Ｖ_nと同じ信号処理を行なうためにハイパスフィルタ６７₁〜６７_n、プログラマブルゲインアンプ６８₁〜６８_n、ローパスフィルタ６９₁〜６９_nが必要となる。ここで、プログラマブルゲインアンプ６３₁〜６３_n、６８₁〜６８_nの増幅率は制御部２０からの制御信号２１１、２１２によって制御可能にしておく。

なお、本発明は、測定された電流と電圧に対して行う信号処理の順序や手段を限定するものではなく、周波数解析部７において正確な周波数解析が行える処理であればどのような順序や手段を用いてもよい。このようにして測定された、測定点＃１〜＃ｎにそれぞれ印加した電流Ｉ₁〜Ｉ_n、各測定点での電圧Ｖ₁〜Ｖ_nは、Ａ／Ｄ変換器７０₁〜７０_n、６５₁〜６５_nによってディジタル信号へと変換され、周波数解析部７へと送出される。

図１に示す周波数解析部７では、測定点＃１〜＃ｎの皮膚にそれぞれ通電された電流Ｉ₁〜Ｉ_nと、測定点＃１〜＃ｎの皮膚の電圧Ｖ₁〜Ｖ_nとに対し周波数解析を行い、その結果を基に各測定点＃１〜＃ｎ直下の皮膚インピーダンスの周波数特性を算出する。以下ではその原理について述べる。今、皮膚インピーダンスが通電電流密度に依存せず、皮膚が線形システムとみなせると仮定する。被検体の皮膚に印加された電流が次式（１）のような電流i(t)である場合、皮膚が線形システムと仮定できる場合、通電によって皮膚に生じた電圧は同様に次式（２）のような電圧v(t)になる。

ここで、ω=２πf（fは測定周波数）、|In|、|Vn|はそれぞれ各周波数成分の電流と電圧の実効値であり、θ_in、θ_vnはそれぞれ各周波数成分fnの電流及び電圧の位相である。式（１）、（２）におけるi(t)、v(t)のフーリエ変換をI(jω)、V(jω)とすると、

であるから、皮膚インピーダンスZ(jω)は

と表すことができる。ただし、

とする。

つまり、必要な周波数fnをi(t)に複数含ませておけば、i(t)、ｖ(t)をフーリエ変換し、それぞれの実数部と虚数部から各周波数におけるインピーダンスが求まる。また、電流波形に含ませる周波数を最低周波数fminの整数倍にすることで、離散フーリエ変換をした場合の隣接する周波数へのスペクトルの洩れは無く、時間分解能の下限は最低周波数のみに依存して決まる。よって、高時間分解能で皮膚インピーダンス軌跡が測定できる。

ところで、皮膚インピーダンスZ(jω)=Rz(ω)-jXz（ω）は次式(9)で表される Cole-Cole の円弧則を満たし、周波数特性の複素平面上での軌跡は、図３に示すような円弧上に
分布することが知られている。

ただし、ω=2πf（fは測定周波数）、τ_mは中心緩和時間、βは円弧の中心角をあらわすパラメータであり、

とする。一般に皮膚インピーダンスの場合、高周波領域のインピーダンスは非常に小さいことからZ_∞≒0とみなせる。式(9)を角周波数ωではなく周波数fで表現すると、

となり、fm = 1/2πτ_mを特性周波数と呼ぶ。式(9)と(12)は等価であるが、以下の説明では式(12)を用いる。

ここで、測定された皮膚インピーダンスの周波数特性の複素平面上での軌跡が円弧上に分布していても、式 (9)、(12)で示されるCole-Cole の円弧則を満たさない場合がある。その場合も含めた一般的な皮膚インピーダンスを表現すると、上式(12)にもう１つのパラメータαを加えた次式(13)(14)にて表現される。

ここで、皮膚インピーダンスが測定電流密度に依存し、皮膚が線形システムとみなせない場合が Cole-Cole の円弧則を満たさない例として挙げられる。従来例においては、測定電極の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映するパラメータに着目している。具体的には、測定された皮膚インピーダンスが Cole-Cole の円弧則を満足すると仮定し、式(12)中の４つのパラメータ（Z₀、Z_∞、β、τ_m=2π/fm）のうち、Cole-Cole の円弧則における中心緩和時間τ_mに着目し、その違いに基づいて経穴の判定を行っている。しかし、前述したように以下の問題点があった。
・測定された皮膚インピーダンスが Cole-Cole の円弧則の成立を確認する手段を設けていない。
・中心緩和時間τ_mの推定手順が煩雑である。
・測定データが皮膚インピーダンス軌跡上に一様に分布しない。

そこで、周波数解析部７で算出された各測定点＃１〜＃ｎでの皮膚インピーダンスが測定電流密度に依存しないことをあらかじめ確認してから、判定部９での経穴位置の判定を行なうようにする。こうすることで、測定された皮膚インピーダンスにおいて Cole-Cole の円弧則が成立するか否かを間接的に確認できるので、特徴量生成部８で生成される特徴量の信頼性・再現性を従来例より高くすることができる。

また、特徴量生成部８においては、周波数解析部７から送出された皮膚インピーダンスZ(jf)の周波数特性から、皮膚インピーダンスZ(jf)の実数部R_Z(f)と虚数部X_Z(f)の比の周波数特性K(ｆ)を算出し、制御部２０より設定される固定周波数f0でのK(f0)を各測定点ごとに生成する。このK(f0)を、従来例におけるパラメータτ_mの代わりに評価パラメータとして採用する。このような構成にすることで、パラメータの推定手順の煩雑さが回避できるほか、前記K(ｆ)と従来例において着目しているパラメータτ_mの大小と正の相関があるので、従来例と同様、測定電極３、４、５の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映する評価パラメータの違いに基づいて経穴の判定を行える。以下では、その理由について述べる。まず、式(14)を変形すると、次式(15)(16)が得られる。

ここで、皮膚インピーダンスが Cole-Cole の円弧則を満たし、高周波領域での皮膚インピーダンスが０とみなせる場合、つまり、α=βとZ_∞=0が成り立つ場合、上式(15)(16)は次式(17)(18)のように表される。

さらに、式(17)(18)の実数部＝R_Z(f) と虚数部＝X_Z(f) の比K(f)＝X_Z(f)／R_Z(f) は次式(19)のようになる。

よって、βが一定であれば、同一測定周波数fでのK(f)の大小はτ_mの大小と正の相関を持ち、fmの大小と負の相関を持つ。

ここで、上記説明が成立するのは、βが一定で、かつ、Z_∞=0でなければならない。一般に皮膚インピーダンスの場合、βは0.65〜0.9と極めて限定された範囲の値をとることが知られており、さらに、測定電極３、４、５の面積が同一で、測定点＃１〜＃ｎが比較的隣接していれば、ほぼ一定とみなすことができる。また、高周波領域における皮膚インピーダンスの大きさは低周波領域と比較すると小さく、特に、上記３電極法によって測定された皮膚インピーダンスの場合、極めて小さいため、Z_∞=0とみなしてよい。

よって、特徴量生成部８において、例えば電流発生部１から発生された電流Ｉには１つの周波数f0の周波数成分を含ませておき、皮膚インピーダンスの実数部と虚数部の比K(f0)を各測定点＃１〜＃ｎごとに算出し、判定部９において、K(f0)が他の測定点よりも低い値になった測定点が経穴であると判定することで、従来例と比べ簡単なハードウエア構成で、かつ、測定電極３、４、５の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映する評価パラメータの違いに基づいて経穴の判定を行えることになる。また、上記K(f0)はある固定周波数f0のみの値であり、従来例のように直接、パラメータτ_mを精度良く求める際のように、極めて低い周波数も含む広範囲の皮膚インピーダンスを測定する必要がない。

なお、前述の判定手法においては、電流Ｉに含ませる周波数はf0１つのみとして説明したが、本発明はそれを限定するものではなく、K(ｆ0)が経穴とそうでない部位での電気的特性の差を明確に反映する周波数であればその数を限定しない。例えばあらかじめ電流発生部１において複数の周波数成分を含む電流Ｉを発生させておいて、その中から各測定点＃１〜＃ｎでのK(f)の値の差が顕著に見られる周波数を特徴量生成部８にて選択し、その周波数でのK(f)を判定部９に送出してもよい。こうすることで、簡単なハードウエア構成で、判定部９での経穴位置の判定の信頼性をより高くすることができるほか、前述したβやZ_∞に関する仮定の成立可否を確認できる。

なお、式(19)中のｆm=τ_m/2πは虚数部X_Zの絶対値が最大値をとる周波数である。よって特徴量生成部８において生成する特徴量を、周波数解析部７において算出される皮膚インピーダンスの周波数特性の虚数部X_Zの絶対値が最大値をとる周波数ｆmとしてもよい。その場合、判定部９では、ｆmが他の測定点よりも大きくなった点を経穴と判定するようにすればよい。こうすることで、従来例と比べ簡単なハードウエア構成で、かつ、測定電極３、４、５の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映する評価パラメータの違いに基づいて経穴の判定を行えることになる。以上のようにして特徴量生成部８で生成された各測定点＃１〜＃ｎでの特徴量や、判定部９での判定結果は表示部１０へと送出され、モニタなどの表示手段により適宜表示される。

以上のように、本発明にかかる経穴位置評価装置は、測定された皮膚インピーダンスにおいて Cole-Cole の円弧則の成立可否を確認でき、評価に用いる特徴量の算出手順が簡単である。また測定結果の信頼性・再現性を従来と比べ高くすることができるほか、評価に用いる特徴量は、従来着目していたパラメータとしてCole-Cole の円弧則における中心緩和時間τ_mと相関があるので、測定電極の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映する評価パラメータの違いに基づいて経穴の判定を行うことが可能になる。よって、従来例に比べ、より簡単なハードウエア構成で、より信頼性・再現性の高い経穴位置判定を行えることから、鍼治療などにおいて、鍼刺激の効果と経穴（ツボ）と呼ばれる刺激点の位置探索や良導点探索など、医療分野にて皮膚の電気的特性の差を非侵襲的かつ客観的に評価する際に有用である。

本発明の経穴位置評価装置の第１の実施の形態における概略構成を示すブロック図図１の電極系と測定部のより詳細な構成例を示すブロック図 Cole-Cole の円弧則を説明する説明図測定データがインピーダンス軌跡上にまんべんなく分布することを説明する図測定データがインピーダンス軌跡上の一部分に偏って分布することを説明する図

符号の説明

１電流発生部
１₁、１₂〜１_n 電流源
２₁、２₂〜２_n 電流検出器
３₁、３₂〜３_n 電流印加電極
４接地電極
５マイナス電極
６測定部
７周波数解析部
８特徴量生成部
９判定部
１０表示部
２０制御部
６１₁〜６１_n 差動増幅器
６２₁〜６２_n、６７₁〜６７_n ハイパスフィルタ（ＨＰＦ）
６３₁〜６３_n、６８₁〜６８_n プログラマブルゲインアンプ（ＰＧＡ）
６４₁〜６４_n、６９₁〜６９_n ローパスフィルタ（ＬＰＦ）
６５₁〜６５_n、７０₁〜７０_n Ａ／Ｄ変換器
２１１、２１２制御信号

Claims

少なくとも３つの周波数成分から構成される電流を発生する電流発生部と、
前記電流発生部からの出力電流を生体の皮膚に配置した複数の電極により複数の異なる皮膚上の測定点にほぼ同時に通電する電極系と、
前記複数の測定点にそれぞれ通電された電流を検出する電流検出器と、
前記電流検出器により検出された電流と、前記電極系の通電によって前記複数の測定点の皮膚に生じた電圧を測定する測定部と、
前記測定部により測定された電流と電圧の両方の周波数解析を行い、各測定点での皮膚インピーダンスの周波数特性を得る周波数解析部と、
前記周波数解析部により解析された前記複数の測定点での皮膚インピーダンスの周波数特性から、前記各測定点での皮膚インピーダンスの特徴量を生成する特徴量生成部と、
前記特徴量生成部により生成された前記各測定点での特徴量から経穴位置を判定する判定部と、
前記特徴量生成部により生成された前記各測定点での特徴量と前記判定部により判定された経穴位置とを表示する表示手段とを、
備え、前記周波数解析部により解析された前記複数の測定点での皮膚インピーダンスが測定電流密度に依存することなくCole-Coleの円弧則が成立した場合に、前記特徴量生成部において生成される特徴量として通電電流に含まれる最低周波数での皮膚インピーダンスの実数部と虚数部との比を算出し、前記判定部で経穴位置を判定するよう構成された経穴位置評価装置。
前記電流発生部において発生させる通電電流は、１つの最低周波数成分と、その整数倍の周波数成分を含むことを特徴とする請求項１に記載の経穴位置評価装置。
前記電流発生部は、通電電流に含ませる周波数成分の大きさを、周波数ごと並びに測定位置ごとに変更可能とすることを特徴とする請求項１又は２に記載の経穴位置評価装置。