CN101179988A - 经穴位置评价装置 - Google Patents

Abstract

本发明的经穴位置评价装置公开了用简单的硬件结构进行有可靠性的经穴判定的技术，根据该技术，电流发生部(1)发生的电流通过电流施加电极(3₁、3₂～3_n)施加于被检测体的皮肤30的检测点(#1～#n)，由测定部(6)测定通电的电流和通电电流在皮肤上产生的电压，由频率分析部(7)对所测定的电流和电压的频率进行分析，计算各检测点的皮肤的阻抗Z(jf)，由特征量生成部(8)根据在检测点#1～#n测得的皮肤阻抗Z(jf)的频率特性计算出皮肤阻抗Z(jf)的实数部RZ(f)与虚数部XZ(f)之比的频率特性K(f)＝XZ(f)/RZ(f)，由判定部(9)根据其频率特性K(f)的不同，从各检测点中判定经穴位置。

Description

经穴位置评价装置

技术领域

本发明涉及在针疗等治疗中能够非侵袭、客观且高精度地判定针刺的效果和称为经穴的刺激点的位置的经穴位置评价装置。

背景技术

很久以来就知道在皮肤表面上存在小区域的皮肤电阻比周边小的部位，其分布与针灸治疗的刺激点即经穴的分布比较一致。这样的皮肤电阻比周边的低的部位(也称为“皮肤电阻减弱点”)在良导络自律神经调整法和低频治疗器等中被用于治疗，这些探索皮肤电阻减弱点的装置也在市场上有销售，在实际临床中使用于诊断和治疗。

以往，作为对这些皮肤电阻减弱点、即经穴进行探索(判定)的装置，测定皮肤的直流电阻作为经穴周边的皮肤的电学特性，根据其差异判定经穴，这种方法已经得到广泛使用，但是，在通过测定这些直流电阻进行经穴判定的方法中，其测定结果依赖于测定电极的接触状态，而且还存在与皮肤电压的自发变动不能够进行明确区别等的问题，测定结果的可靠性和重现性差。

因此，近年来作为提高这些测定结果的可靠性、重现性的方法，有下述方法，根据对经穴周边的皮肤施加交流电流而生成的电压，测定皮肤阻抗的频率特性在复数平面上的分布(皮肤阻抗轨迹)，作为经穴周边的皮肤的电学特性，根据皮肤阻抗轨迹的差异进行经穴的判定。例如作为表现皮肤阻抗轨迹的形状这样的特征的4个参数(下述Z₀、Z_∞、β、τ_m)中与测定电极的接触状态没有太大关系的、明确反映经穴与非经穴部位的电学特性的差异的参数，着眼于Cole-Cole的圆弧法则的中心张驰时间τ_m，根据该参数τ_m的不同进行经穴的判定，通过这样做能够避免已有技术存在的问题(参照例如下述专利文献1)。

专利文献1：日本特开2004-337349号公报(说明书摘要、权利要求2)

但是，在上述已有技术中，对经穴和非经穴部位的电学特性的差异进行评价时，作为带有皮肤阻抗轨迹的形状的特征的参数，计算出Cole-Cole的圆弧法则中的中心张驰时间τ_m，但是存在如下所述问题。

在计算出表示皮肤阻抗轨迹的形状这样的特征的中心张驰时间τ_m时，在皮肤阻抗的频率特性中Cole-Cole的圆弧法则成立是必要的，但是没有确认它是否能够成立的手段。而且有必要根据所测定的多个频率的皮肤阻抗利用非线性最小二乘法等推定中心张驰时间τ_m，由于推定用的步骤比较复杂，实现它的硬件的结构复杂，不利于装置的小型化。而且为了高精度推定中心张驰时间τ_m，必须如图4所示，在皮肤阻抗轨迹上测定的皮肤阻抗数据没有偏颇地普遍分布。

但是，实际的皮肤上测定的皮肤阻抗数据的分布如图5所示，往往偏向阻抗轨迹上的一部分存在，而且其分布的方式往往对于各检测点有很大的不同。在这样的情况下，为了高精度地推定表示皮肤阻抗轨迹的形状这样的特征的中心张驰时间τ_m，有必要测定也包含极低的频率范围的很宽的频率范围的皮肤阻抗。因此，除了时间分辨率低外，测定对象频率越低则越容易受到温度漂移和身体运动造成的外部干扰的影响，难于保证测定结果的可靠性和重现性。

发明内容

本发明是为了解决上述已有技术存在的问题而作出的，其目的在于，提供一种能够用简单的硬件结构进行具有可靠性和重现性的经穴判定的经穴位置评价装置。

为了解决上述已有技术存在的问题，本发明的经穴位置评价装置，具备发生至少由一个频率分量构成的电流的电流发生部、

利用配置于生物体的皮肤上的多个电极将所述电流发生部的输出电流大致同时地对多个不同的皮肤上的检测点通电的电极系统、

检测对所述多个检测点的分别通电的电流的电流检测器、

测定所述电流检测器检测出的电流和所述电极系统的通电在所述多个检测点的皮肤上产生的电压的测定部、

对所述测定部测定的电流和电压两者进行频率分析，取得各检测点的皮肤的阻抗的频率特性的频率分析部、

根据所述频率分析部分析的所述多个检测点的皮肤的阻抗的频率特性，生成所述各检测点的皮肤阻抗的特征量的特征量生成部、

根据所述特征量生成部生成的所述各检测点的特征量判定经穴位置的判定部、以及

显示所述特征量生成部生成的所述各检测点的特征量与所述判定部判定的经穴位置的显示装置。

利用这样的结构，用于评价的特征量的计算步骤将变得简单，能够使测定结果的可靠性和重现性比以往高，除此以外，用于评价的特征量由于与一向注重的中心张驰时间τ_m有关，所以能够与测定电极的接触状态不大相干地，根据明确反映经穴与非经穴部位的电学特性之差的评价参数的不同进行经穴的判定。

又，本发明的经穴位置评价装置，其特征在于，所述频率分析部确认在所述各检测点的皮肤阻抗的频率特性中Cole-Cole的圆弧法则是否成立。

利用这样的结构，能够确认在所测定的皮肤阻抗中Cole-Cole的圆弧法则可否成立，能够使得测定结果的可靠性和重现性比以往高。

又，本发明的经穴位置评价装置，其特征在于，所述电流发生部发生的通电电流包含一个最低频率分量和其整数倍的频率分量。

利用这样的结构，由于在进行离散傅利叶变换情况下频谱没有向相邻的频率浊漏，时间分辨率的下限只由最低频率决定，因此能够以高时间分辨率测定皮肤阻抗轨迹。

又，本发明的经穴位置评价装置，其特征在于，

所述电流发生部能够对每一频率、以及每一测定位置改变通电电流中包含的频率分量的大小。

利用这样的结构，能够对各频率、各检测点确保最合适的信噪比，能够以高精度测定皮肤阻抗轨迹。

又，本发明的经穴位置评价装置，其特征在于，

在所述特征量生成部生成的特征量是至少一个频率的皮肤阻抗的实数部与虚数部之比。

利用这样的结构，能够避免参数的推定步骤的复杂性。

又，本发明的经穴位置评价装置，其特征在于，

在所述特征量生成部以至少一个频率的皮肤阻抗的实数部与虚数部之比中在所述各检测点差值显著的频率上的比为特征值。

利用这样的结构，除了能够避免参数的推定步骤的复杂性外，能够使在判定部进行的经穴位置判定有更高的可靠性。

又，本发明的经穴位置评价装置，其特征在于，在所述特征量生成部生成的特征量是通电电流中包含的最低的频率下的皮肤阻抗的实数部与虚数部之比。

利用这样的结构，除了能够避免参数的推定步骤的复杂性外，能够迅速进行经穴位置的判定。

如果采用本发明，则能够确认在所测定的皮肤阻抗中Cole-Cole的圆弧法则是否成立，能够用简单的硬件实现使用于评价的特征量的计算步骤。而且除了能够使测定结果的可靠性和重现性比以往更高外，由于使用于评价的特征量与向来重视的中心张驰时间τ_m相关，所以能够与测定电极的接触状态不大相干地，根据明确反映经穴与非经穴部位的电学特性之差的评价参数的不同进行经穴的判定。从而，能够以简单的硬件结构提供可靠性和重现性更高的经穴位置装置。

附图说明

图1是表示本发明的经穴位置评价装置的第1实施形态的大概结构的方框图。

图2是表示图1的电极系统与测定部的更详细的结构例的方框图。

图3是说明Cole-Cole的圆弧法则的说明图。

图4是说明测定的数据在阻抗轨迹上普遍分布的说明图。

图5是说明测定的数据偏向于阻抗轨迹上的一部分分布的说明图。

具体实施方式

下面参照附图对本发明的实施形态进行说明。图1是表示本发明的经穴位置评价装置的第1实施形态的概要结构的方框图。电流发生部1发生的交流电流I₁～I_n分别通过电流检测器2₁～2_n、电流施加电极3₁～3_n施加于被检测物体的皮肤30的各检测点#1～#n，流向接地电极4。以接地电极4的电位为基准，利用测定部6测定因该通电产生的、在电流施加电极3₁～3_n与负电极5(-)之间的皮肤上生成的电压V₁～V_n的下降。利用这样的电极系统测定的方法被称为3电极法，在测定电流施加电极3₁～3_n的正下方、也就是检测点#1～#n的正下方的皮肤阻抗时经常被采用。

图2是表示电流发生部1以及测定部6的更详细的结构例的方框图。电流发生部1由至少1个电流源I₁～I_n构成，电流源I₁～I_n分别发生向各检测点#1～#n通电的交流电流I₁～I_n。在这里，各电流源I₁～I_n发生的电流I₁～I_n中包含的频率分量的数目、频率、大小可利用来自控制部20的控制信号210进行设定。测定部6由用于测定因通电在检测点#1～#n正下方的皮肤上产生的电压V₁～V_n的至少一个差动放大部61₁～61_n、至少一个高通滤波器(HPF)62₁～62_n(以及67₁～67_n)、至少一个可编程增益放大器(PGA)63₁～63_n(以及68₁～68_n)、至少一个低通滤波器(LPF)64₁～64_n(以及69₁～69_n)、以及至少一个A/D变换器65₁～65_n(以及70₁～70_n)构成。

分别利用差动放大器61₁～61_n测定因通电在检测点#1～#n的皮肤上产生的电压V₁～V_n，测得的电压V₁～V_n中包含不是因通电引起的不需要的信号分量。这妨碍频率分析部7进行正确的频率分析，因此，利用高通滤波器62₁～62_n分别从差动放大器61₁～61_n的输出中去除不需要的低频分量。不需要的低频分量被去除后的检测点#1～#n的电压V₁～V_n分别根据需要由可编程增益放大器63₁～63_n放大，接着利用低通滤波器64₁～64_n去除不需要的高频分量。

又，分别由电流检测器2₁～2_n测定对检测点#1～#n的皮肤通电的电流I₁～I_n，如果对向检测点#1～#n的皮肤通电的电流I₁～I_n进行与在检测点#1～#n的电压V₁～V_n相同的信号处理，则如图1所示的频率分析部7进行的频率分析的步骤将变得复杂，分析结果也不正确。因此为了进行与在检测点#1～#n的电压V₁～V_n相同的信号处理，需要高通滤波器67₁～67_n、可编程增益放大器68₁～68_n、低通滤波器69₁～69_n。在这里，可编程增益放大器63₁～63_n、68₁～68_n的放大率可预先利用来自控制部20的控制信号211、212进行控制。

还有，本发明不限定对所测定的电流和电压进行的信号处理的步骤和手段，只要是在频率分析部7能够进行正确的频率分析的处理，无论采用哪一种步骤和手段都可以。利用A/D变换器70₁～70_n、65₁～65_n将这样测定到的分别施加于检测点#1～#n的电流I₁～I_n、各测定点的电压V₁～V_n变换为数字信号，发送到频率分析部7。

在图1所示的频率分析部7，对检测点#1～#n的皮肤分别通电的电流I₁～I_n以及检测点#1～#n的皮肤的电压V₁～V_n进行频率分析，根据其结果计算出各检测点#1～#n正下方的皮肤阻抗的频率特性。以下对其原理进行叙述。现在假定皮肤阻抗与通电电流密度无关，皮肤可以看作线性系统。施加于被检测体的皮肤上的电流为下式(1)那样的电流i(t)的情况下，皮肤可以假定为线性系统的情况下，因通电在皮肤上产生的电压同样为下式(2)的电压v(t)。

(数学式1)

$i\left(t\right)=\underset{n}{\mathrm{\Σ}}\sqrt{2}\left|I_n\right|\sin ({\mathrm{\ω}}_{n}t+{\mathrm{\θ}}_{\mathrm{in}})---\left(1\right)$

$v\left(t\right)=\underset{n}{\mathrm{\Σ}}\sqrt{2}\left|V_n\right|\sin ({\mathrm{\ω}}_{n}t+{\mathrm{\θ}}_{v_n})---\left(2\right)$

在这里，ω＝2πf(f为测定频率)，|In|、|Vn|分别为各频率分量的电流和电压的有效值，θ_in、θ_vn分别为各频率分量fn的电流和电压的相位。将式(1)、(2)中的i(t)、v(t)的傅里叶变换作为I(jω)、V(jω)时，由于

(数学式2)

$I\left(\mathrm{j\ω}\right)=\underset{n}{\mathrm{\Σ}}\sqrt{2}\mathrm{\π}\left|I_n\right|(\cos {\mathrm{\θ}}_{i_n}-j\sin {\mathrm{\θ}}_{i_n})\mathrm{\δ}(\mathrm{\ω}-{\mathrm{\ω}}_n)$

$=\underset{n}{\mathrm{\Σ}}(R_{I_n}\left(\mathrm{\ω}\right)-j{X}_{J_n}\left(\mathrm{\ω}\right))\·\mathrm{\δ}(\mathrm{\ω}-{\mathrm{\ω}}_n)---\left(3\right)$

$V\left(\mathrm{j\ω}\right)=\underset{n}{\mathrm{\Σ}}\sqrt{2}\mathrm{\π}\left|V_n\right|(\cos {\mathrm{\θ}}_{v_n}-j\sin {\mathrm{\θ}}_{v_n})\mathrm{\δ}(\mathrm{\ω}-{\mathrm{\ω}}_n)$

$=\underset{n}{\mathrm{\Σ}}(R_{V_n}\left(\mathrm{\ω}\right)-j{X}_{V_n}\left(\mathrm{\ω}\right))\·\mathrm{\δ}(\mathrm{\ω}-{\mathrm{\ω}}_n)---\left(4\right)$

δ(ω)：狄拉克δ函数

所以皮肤阻抗Z(jω)可以表示为

(数学式3)

$Z\left(\mathrm{j\ω}\right)=\underset{\text{n}}{\mathrm{\Σ}}(R_{Z_n}-j{X}_{Z_n})\·\mathrm{\δ}(\mathrm{\ω}-{\mathrm{\ω}}_n)---\left(5\right)$

并且假定

(数学式4)

$R_{Z_n}=\frac{R_{I_n}\left(\mathrm{\ω}\right)R_{V_n}\left(\mathrm{\ω}\right)+X_{I_n}\left(\mathrm{\ω}\right)X_{V_n}\left(\mathrm{\ω}\right)}{R_{I_n}^2\left(\mathrm{\ω}\right)+X_{I_n}^2\left(\mathrm{\ω}\right)}---\left(6\right)$

$X_{Z_n}=\frac{R_{I_n}\left(\mathrm{\ω}\right)X_{V_n}\left(\mathrm{\ω}\right)-X_{I_n}{R}_{V_n}\left(\mathrm{\ω}\right)}{R_{I_n}^2\left(\mathrm{\ω}\right)+X_{I_n}^2\left(\mathrm{\ω}\right)}---\left(7\right)$

也就是说，如果使i(t)中包含多个所需要的频率fn，则对i(t)、v(t)进行傅里叶变换，能够从各实数部和虚数部求出各频率的阻抗。又，通过使电流波形中包含的频率为最低频率f_min的整数倍，可以消除进行离散傅里叶变换的情况下光谱向相邻的频率浊漏，时间分辨率的下限只取决于最低频率。因此能够用高时间分辨率测定皮肤阻抗轨迹。

但是，皮肤阻抗Z(jω)＝Rz(ω)-jxz(ω)满足下式(9)表示的Cole-Cole的圆弧法则，频率特性的复数平面上的轨迹，已知如图3所示分布在圆弧上。

(数学式5)

$Z\left(\mathrm{j\ω}\right)=R_Z\left(\mathrm{\ω}\right)-j{X}_Z\left(\mathrm{\ω}\right)---\left(8\right)$

$=Z_{\∞}+\frac{Z_0-Z_{\∞}}{1+{\left(\mathrm{j\ω}{\mathrm{\τ}}_m\right)}^{\mathrm{\β}}}---\left(9\right)$

其中，ω＝2πf(f为测定频率)、τ_m是中心张弛时间、β为表示圆弧的中心角的参数，假定

(数学式6)

$Z_0=\underset{\mathrm{\ω}\→0}{\lim }Z\left(\mathrm{\ω}\right),Z_{\∞}=\underset{\mathrm{\ω}\→\∞}{\lim }Z\left(\mathrm{\ω}\right)$

通常在皮肤阻抗的情况下，高频区域的阻抗非常小，因此可以看作Z_∞0。式(9)用频率f而不用角频率ω表示时，

(数学式7)

$Z\left(\mathrm{jf}\right)=R_Z\left(f\right)-j{X}_Z\left(f\right)---\left(10\right)$

$=Z_{\∞}+\frac{Z_0-Z_{\∞}}{1+{\left(j2\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{\mathrm{\β}}}---\left(11\right)$

$=Z_{\∞}+\frac{Z_0-Z_{\∞}}{1+{\left(j\frac{f}{\mathrm{fm}}\right)}^{\mathrm{\β}}}---\left(12\right)$

fm＝1/2πτ_m被称为特性频率。式(9)和(12)是等效的，但是在下面说明中采用式(12)。

在这里，即使所测定的皮肤阻抗的频率特性在复数平面上的轨迹分布在圆弧上，也有Cole-Cole的圆弧法则得不到满足的情况。若表达也包含这种情况的一般的皮肤阻抗，则用在上式(12)中加上又一个参数α的下式(13)、(14)表达。

(数学式8)

$Z\left(\mathrm{jf}\right)=R_Z\left(f\right)-j{X}_Z\left(f\right)---\left(13\right)$

$=Z_{\∞}+\frac{Z_0-Z_{\∞}}{1+j^{\mathrm{\β}}{\left(2\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{\mathrm{\α}}}=Z_{\∞}+\frac{Z_0-Z_{\∞}}{1+j^{\mathrm{\β}}{\left(\frac{f}{\mathrm{fm}}\right)}^{\mathrm{\α}}}---\left(14\right)$

在这里，皮肤阻抗取决于测定电流密度，作为Cole-Cole的圆弧法则得不到满足例子可以举出皮肤不能看作线性系统的情况。在以往的例子中，着眼于与测定电极的接触状态没有太大关系，明确反映经穴与非经穴部位的电气特性的差别的参数。具体地说，假定所测定的皮肤阻抗满足Cole-Cole的圆弧法则，式(12)中的四个参数(Z₀、Z_∞、β、τ_m＝2π/fm)中，着眼于Cole-Cole的圆弧法则中的中心张弛时间τ_m，根据该差异进行经穴的判定。但是如上所述存在以下所述问题。

·没有设置确认所测定的皮肤阻抗是否使Cole-Cole的圆弧法则成立的手段。

·中心张弛时间τ_m的推定步骤繁杂。

·测定数据在皮肤阻抗轨迹上分布不相同。

因此，在预先确认频率分析部7计算出的各检测点#1～#n的皮肤阻抗与测定电流密度无关之后在判定部9进行经穴位置的判定。通过这样做，能够间接确认所测定的皮肤阻抗中Cole-Cole的圆弧法则是否成立，因此能够使得特征量生成部8生成的特征量的可靠性、重现性比以往的例子高。

又，在特征量生成部8中，根据频率分析部7发送出的皮肤阻抗Z(jf)的频率特性计算出皮肤阻抗Z(jf)的实数部RZ(f)与虚数部XZ(f)之比的频率特性K(f)，对每一个测定点生成控制部20所设定的固定频率f0下的K(f0)。采用该K(f0)取代已有例中的参数τ_m作为评价参数。通过形成这样的结构，除了可以避免参数的推定步骤的繁杂外，由于所述K(f)与已有例中着眼的参数τ_m的大小有正的相关，因此与已有例一样，与测定电极即电流施加电极3₁～3_n、接地电极4、负电极5的接触状态没有很大的关系，可以根据明确反映经穴与非经穴的部位的电气特性的差异的评价参数的不同进行经穴的判定。以下对其理由进行叙述。首先，使式(14)变形，就得到下式(15)、(16)。

(数学式9)

$\frac{R_Z\left(f\right)-Z_{\∞}}{Z_0-Z_{\∞}}=\frac{1+{\left(2\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{\mathrm{\α}}\cos \frac{\mathrm{\β\π}}{2}}{1+2{\left(1\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{\mathrm{\β}}\cos \frac{\mathrm{\β\π}}{2}+{\left(\mathrm{\ω}{\mathrm{\τ}}_m\right)}^{2\mathrm{\β}}}=\frac{1+{(f/f_m)}^{\mathrm{\α}}\cos \frac{\mathrm{\β\π}}{2}}{1+2{(f/f_m)}^{\mathrm{\β}}\cos \frac{\mathrm{\β\π}}{2}+{(f/f_m)}^{2\mathrm{\β}}}---\left(15\right)$

$\frac{X_Z\left(f\right)}{Z_0-Z_{\∞}}=\frac{{\left(2\mathrm{\πf}{\mathrm{\γ}}_m\right)}^{\mathrm{\α}}\sin \frac{\mathrm{\β\π}}{2}}{1+2{\left(2\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{\mathrm{\α}}\cos \frac{\mathrm{\β\π}}{2}+{\left({\mathrm{\ω\τ}}_m\right)}^{2\mathrm{\β}}}=\frac{{(f/f_m)}^{\mathrm{\α}}\sin \frac{\mathrm{\β\π}}{2}}{1+2{(f/f_m)}^{\mathrm{\α}}\cos \frac{\mathrm{\β\π}}{2}+{(f/f_m)}^{2\mathrm{\β}}}---\left(16\right)$

在这里，皮肤阻抗满足Cole-Cole的圆弧法则，高频区域的皮肤阻抗看作0的情况下，也就是α＝β和Z_∞＝0成立的情况下，上式(15)、(16)如下式(17)、(18)所示表达。

(数学式10)

$\frac{R_Z\left(f\right)}{Z_0}=\frac{1+{\left(2\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{\mathrm{\β}}\cos \frac{\mathrm{\β\π}}{2}}{1+2{\left(2\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{\mathrm{\β}}\cos \frac{\mathrm{\β\π}}{2}+{\left(2\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{2\mathrm{\β}}}=\frac{1+{(f/f_m)}^{\mathrm{\β}}\cos \frac{\mathrm{\β\π}}{2}}{1+2{(f/f_m)}^{\mathrm{\β}}\cos \frac{\mathrm{\β\π}}{2}+{(f/f_m)}^{2\mathrm{\β}}}---\left(17\right)$

$\frac{X_Z\left(f\right)}{Z_0}=\frac{{\left(2\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{\mathrm{\β}}\sin \frac{\mathrm{\β\π}}{2}}{1+2{\left(2\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{\mathrm{\β}}\cos \frac{\mathrm{\β\π}}{2}+{\left(2\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{2\mathrm{\β}}}=\frac{{(f/f_m)}^{\mathrm{\β}}\sin \frac{\mathrm{\β\π}}{2}}{1+2{(f/f_m)}^{\mathrm{\β}}\cos \frac{\mathrm{\β\π}}{2}+{(f/f_m)}^{2\mathrm{\β}}}---\left(18\right)$

而且式(17)、(18)的实数部＝RZ(f)与虚数部＝XZ(f)之比K(f)＝XZ(f)/RZ(f)如下式(19)所示。

(数学式11)

$K\left(f\right)=\frac{X_Z\left(f\right)}{R_Z\left(f\right)}=\frac{{\left(2\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{\mathrm{\β}}\sin \frac{\mathrm{\β\π}}{2}}{1+{\left(2\mathrm{\πf}{\mathrm{\τ}}_m\right)}^{\mathrm{\β}}\cos \frac{\mathrm{\β\π}}{2}}=\frac{{(f/f_m)}^{\mathrm{\β}}\sin \frac{\mathrm{\β\π}}{2}}{1+{(f/f_m)}^{\mathrm{\β}}\cos \frac{\mathrm{\β\π}}{2}}---\left(19\right)$

因此，如果β为一定，则在同一测定频率f下的K(f)的大小与τ_m有正的相关关系，与fm的大小有负的相关关系。

在这里，上述说明成立必须是β为一定，而且Z_∞＝0。通常在皮肤阻抗的情况下，已知β取0.65～0.9这样极其有限的范围的数值，而且测定电极3、4、5的面积相同，如果检测点#1～#n相互比较靠近，则可以看作大致一定。高频区域的皮肤阻抗的大小与低频区域相比比较小，特别是在利用上述3电极法测定的皮肤阻抗的情况下，数值极小，因此可以看作Z_∞＝0。

因此，在特征量生成部8中，使例如电流发生部1发生的电流I中包含一个频率f0的频率分量，对每一个检测点#1～#n计算出皮肤阻抗的实数部与虚数部之比K(fD)，在判定部9中，将K(f0)比其他检测点数值低的检测点判定为经穴，由此，能够以比已有技术例简单的硬件结构，而且根据与测定电极3、4、5的接触状态关系不大的、明确反映经穴与非经穴部位的电气特性之差的评价参数的不同进行经穴的判定。而且上述K(f0)只是某固定频率f0的值，不需要像已有技术的例子那样直接地像高精度求参数τ_m时那样测定也包含极低频率的宽范围的皮肤阻抗。

还有，上述判定方法中，对于电流I中包含的频率只有一个f0的情况进行了说明，但是本发明并不限于此，只要是K(f0)明确反映经穴与非经穴的部位的电气特性的差的频率，不限定其个数。也可以例如预先使电流发生部1中发生包含多个频率分量的电流I，由特征量生成部8从其中选择在格检测点#1～#n的K(f)的值的差别明显能够看得出的频率，将该频率下的K(f)发送到判定部9。通过这样做，能够以简单的硬件结构进一步提高判定部9进行的经穴位置判定的可靠性，除此以外还能够确认关于上述β和Z_∞的假定是否能够成立。

还有，式(19)中的fm＝τ_m/2π是虚数部XZ的绝对值取最大值的频率。因此也可以将特征量生成部8生成的特征量作为频率分析部7计算出的皮肤阻抗的频率特性的虚数部XZ的绝对值取最大值的频率fm。在这种情况下，判定部9只要将fm比其他检测点大的点判定为经穴即可。通过这样做，能够以比已有技术的例子简单的硬件结构、而且能够根据和与测定电极3、4、5的接触状态关系不大的、明确反映经穴与非经穴部位的电气特性的差异的评价参数的不同，进行经穴的判定。这样在特征量生成部8生成的各检测点#1～#n的特征量和判定部9进行判定的结果被发送到显示部10，利用监视器等显示装置适当进行显示。

工业应用性

如上所述，本发明的经穴位置评价装置能够确认在所测得的皮肤阻抗中Cole-Cole的圆弧法则是否能够成立，评价中使用的特征量的计算步骤简单。而且测定结果的可靠性、重现性比以往高，除此以外，评价中使用的特征量，作为以往关注的参数与Cole-Cole的圆弧法则中的中心张弛时间τ_m存在相关性，因此能够根据和与测定电极的接触状态关系不大的、明确反映经穴与非经穴部位的电气特性的差异的评价参数的不同，进行经穴的判定。从而，与已有技术的例子相比，能够以更简单的硬件结构进行更可靠的重现性高的经穴位置判定，因此在针刺治疗等方法中，在针刺刺激的效果和被称为经穴的刺激点的位置的探索以及良导电点的探索等医疗领域中，以非侵袭方式且客观地评价皮肤的电气特性的差异时，是有用的。

Claims (7)

1.一种经穴位置评价装置，其特征在于，具备

电流发生部，发生至少由一个频率分量构成的电流、

电极系统，利用配置于生物体的皮肤上的多个电极将来自所述电流发生部的输出电流大致同时地对多个不同的皮肤上的检测点通电、

电流检测器，检测分别对所述多个检测点的通电的电流、

测定部，测定所述电流检测器检测出的电流和由所述电极系统的通电在所述多个检测点的皮肤上产生的电压、

频率分析部，进行对所述测定部测定的电流和电压两者的频率分析，取得各检测点的皮肤的阻抗的频率特性、

特征量生成部，根据所述频率分析部分析的所述多个检测点的皮肤的阻抗的频率特性，生成所述各检测点的皮肤阻抗的特征量、

判定部，根据所述特征量生成部生成的所述各检测点的特征量判定经穴位置、以及

显示装置，显示所述特征量生成部生成的所述各检测点的特征量与所述判定部判定的经穴位置。

2.根据权利要求1所述的经穴位置评价装置，其特征在于，所述频率分析部确认在所述各检测点的皮肤阻抗的频率特性中Cole-Cole的圆弧法则是否成立。

3.根据权利要求1所述的经穴位置评价装置，其特征在于，所述电流发生部发生的通电电流包含一个最低频率分量和其整数倍的频率分量。

4.根据权利要求1所述的经穴位置评价装置，其特征在于，所述电流发生部能够对每一频率、以及每一测定位置改变通电电流中包含的频率分量的大小。

5.根据权利要求1所述的经穴位置评价装置，其特征在于，在所述特征量生成部生成的特征量是至少一个频率的皮肤阻抗的实数部与虚数部之比。

6.根据权利要求1所述的经穴位置评价装置，其特征在于，在所述特征量生成部以至少一个频率的皮肤阻抗的实数部与虚数部之比中在所述各检测点差值显著的频率上的比为特征值。

7.根据权利要求1所述的经穴位置评价装置，其特征在于，在所述特征量生成部生成的特征量是通电电流中包含的最低的频率下的皮肤阻抗的实数部与虚数部之比。

Images