JP2010187993A - Blood pressure measuring device and program as well as recording medium - Google Patents

Blood pressure measuring device and program as well as recording medium Download PDF

Info

Publication number
JP2010187993A
JP2010187993A JP2009036789A JP2009036789A JP2010187993A JP 2010187993 A JP2010187993 A JP 2010187993A JP 2009036789 A JP2009036789 A JP 2009036789A JP 2009036789 A JP2009036789 A JP 2009036789A JP 2010187993 A JP2010187993 A JP 2010187993A
Authority
JP
Japan
Prior art keywords
blood pressure
learning
estimated
data
determination
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2009036789A
Other languages
Japanese (ja)
Other versions
JP5218139B2 (en
Inventor
Koki Futatsuyama
幸樹 二ツ山
Shinji Nanba
晋治 難波
Koji Oguri
宏次 小栗
Satomi Suzuki
里実 鈴木
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Denso Corp
Original Assignee
Denso Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Denso Corp filed Critical Denso Corp
Priority to JP2009036789A priority Critical patent/JP5218139B2/en
Publication of JP2010187993A publication Critical patent/JP2010187993A/en
Application granted granted Critical
Publication of JP5218139B2 publication Critical patent/JP5218139B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a blood pressure measuring device which can estimate blood pressure with high accuracy by carrying out good learning to blood pressure for which truth or falsehood is ambiguous, a program, and a recording medium. <P>SOLUTION: In a step 110, a pulse wave signal and an electrocardiogram signal are analyzed and a characteristic amount required for blood pressure measurement is computed. In a step 130, two or more binary discrimination are performed. Namely, the characteristic amount is used for a discrimination function established for each blood pressure for every 1 mmHg, being under or equal to above a blood pressure value for the determination subject is determined, and a binary series, which has the number of digits corresponding to the number of blood pressure for the determination subject, is created. In a step 140, the binary series and a code table, which are created in the step 130, are compared, Hamming distances are calculated, and a blood pressure value corresponding to the minimum Hamming distance is chosen. In a step 150, the determination of reliability is performed by the number of learning data. <P>COPYRIGHT: (C)2010,JPO&INPIT

Description

本発明は、例えば車両に搭載されて、カフ等を用いることなくドライバなどの血圧の測定が可能な血圧測定装置及びプログラム並びに記録媒体に関する。   The present invention relates to a blood pressure measurement device, a program, and a recording medium that are mounted on a vehicle and can measure blood pressure such as a driver without using a cuff or the like.

従来より、血圧を測定する場合には、聴診法においてもオシロメトリック法においても、カフによる締め付けが必要であり、装置が大型化したり、被測定者に苦痛があったり、連続的な血圧測定ができないなどの問題があった。   Conventionally, when measuring blood pressure, both the auscultation method and the oscillometric method require tightening with a cuff, the device becomes larger, the subject is suffering, and continuous blood pressure measurement is required. There were problems such as being unable to do so.

これに対して、近年では、脈波伝播速度や脈波特徴量などを用いて血圧を算出する測定法が検討されている(特許文献1〜10参照)。例えば、容積脈波であれば、光を用いて測定できるため、装置の小型化を実現できるととともに、被測定者の苦痛を解消でき、更には、連続的な血圧測定が可能になる。   On the other hand, in recent years, measurement methods for calculating blood pressure using a pulse wave velocity, a pulse wave feature amount, and the like have been studied (see Patent Documents 1 to 10). For example, volume pulse waves can be measured using light, so that the apparatus can be miniaturized, the pain of the measurement subject can be eliminated, and continuous blood pressure measurement can be performed.

また、血圧を推定する手法として、重回帰やニューラルネットワーク(NN)による学習なども提案されている(特許文献11、12参照)。   Further, as a technique for estimating blood pressure, learning using multiple regression, neural network (NN), and the like have been proposed (see Patent Documents 11 and 12).

特開平10−295656号公報JP-A-10-295656 特開平10−295657号公報Japanese Patent Laid-Open No. 10-295657 特開平11−318837号公報Japanese Patent Laid-Open No. 11-318837 特表2001−504362号公報JP-T-2001-504362 特表2003−555号公報Special table 2003-555 gazette 特表2003−527149号公報Special table 2003-527149 gazette 特開2006−263354号公報JP 2006-263354 A 特開2006−6897号公報JP 2006-6897 A 特開2007−7075号公報JP 2007-7075 A 特開2007−7077号公報JP 2007-7077 A 特開2000−126142号公報JP 2000-126142 A 特開2000−166885号公報JP 2000-166885 A

上述した特許文献1〜12の技術は、例えば脈波や心電の信号の特徴量(例えば脈波の2階微分、脈波伝播速度(PTT)、心拍周期(RR))などを用いて、或いは、重回帰やNNの手法を利用して、血圧を推定するものであるが、かならずしも十分ではない。   The techniques of Patent Documents 1 to 12 described above use, for example, the feature amount of a pulse wave or an electrocardiogram signal (for example, the second derivative of the pulse wave, the pulse wave velocity (PTT), the heart rate cycle (RR)), etc. Alternatively, blood pressure is estimated using multiple regression or NN techniques, but this is not always sufficient.

重回帰においては、推定を線形的に行うために、推定精度に限界があり、また、NNによる学習では計算時間が長く、かつ、学習のデータを真値としてしまうため、真値が曖昧な血圧推定に適用が難しいという課題があった。   In multiple regression, since estimation is performed linearly, estimation accuracy is limited, and learning by NN takes a long calculation time, and the learning data becomes a true value. There was a problem that it was difficult to apply to estimation.

つまり、血圧推定では、市販のカフ式の血圧計で計測した値を真値として学習する場合、カフ式の血圧計の測定精度(誤差の標準偏差)が4〜8mmHgばらつくため、誤ったデータを基に学習せざるを得ず、推定精度が低下してしまっていた。(参照:JIS T 1115 (非観血式電子血圧計)、http://www.healthcare.omron.co.jp/medical/product/technique/pdf/technique_nO1.pdf)
本発明は、上記点に鑑み、測定精度がばらつき易い血圧に対して高い精度で血圧を推定できる血圧測定装置及びプログラム並びに記録媒体を提供することを目的とする。
That is, in blood pressure estimation, when learning a value measured with a commercially available cuff sphygmomanometer as a true value, the measurement accuracy (standard deviation of error) of the cuff sphygmomanometer varies from 4 to 8 mmHg. Based on this, we had to learn, and the estimation accuracy had been reduced. (Reference: JIS T 1115 (Non-invasive electronic blood pressure monitor), http://www.healthcare.omron.co.jp/medical/product/technique/pdf/technique_nO1.pdf)
In view of the above points, an object of the present invention is to provide a blood pressure measurement device, a program, and a recording medium that can estimate blood pressure with high accuracy with respect to blood pressure whose measurement accuracy is likely to vary.

(1)請求項1の発明は、生体から検出された心電信号及び脈波信号のうち、少なくとも前記脈波信号から得られる生体情報を用いて、血圧を推定する血圧測定装置において、学習データである前記生体情報の特徴量と血圧との関係を予め機械学習することにより得られた識別手段であって、所定の血圧毎に前記生体情報の特徴量から血圧を2値化判別する学習識別手段と、前記血圧の推定を行う場合に、前記学習識別手段を用い、測定によって得られた前記生体情報の特徴量に基づいて、前記所定の血圧に対して2値化判別を行う2値化判別手段と、前記2値化判別手段によって、前記生体情報の特徴量に基づいて異なる前記所定の血圧に対してそれぞれ2値化判定を行い、該複数回の2値化判定の結果に基づいて前記血圧を推定する血圧推定手段と、を備えたことを特徴とする。   (1) The invention of claim 1 is a blood pressure measurement apparatus that estimates blood pressure using at least biological information obtained from a pulse wave signal among electrocardiogram signals and pulse wave signals detected from a living body. Learning means for obtaining a binarized determination of blood pressure from the feature value of the biological information for each predetermined blood pressure, which is obtained by previously learning the relationship between the feature value of the biological information and the blood pressure. And binarization for performing binarization determination on the predetermined blood pressure based on the feature amount of the biological information obtained by measurement using the learning identification unit when estimating the blood pressure. The determination unit and the binarization determination unit respectively perform binarization determination on the different predetermined blood pressures based on the feature amount of the biological information, and based on the results of the multiple binarization determinations Blood pressure to estimate the blood pressure Characterized by comprising a constant means.

本発明では、予め学習データである生体情報の特徴量と血圧との関係を機械学習することにより、所定の血圧毎に(例えば100mmHgや101mmHg等の1mmHg毎などに)生体情報の特徴量から血圧を2値化判別する学習識別手段(例えば識別関数:評価関数:判別関数)を生成する。次に、実際に血圧の推定を行う場合には、学習識別手段を用い、測定によって得られた生体情報の特徴量に基づいて、所定の血圧(例えば100mmHgや101mmHg等の1mmHg毎など)に対して、それぞれ2値化判別を行う。そして、この結果得られた複数回の判定データに基づいて血圧を推定する。   In the present invention, the blood pressure is calculated from the feature value of the biometric information for each predetermined blood pressure (for example, every 1 mmHg such as 100 mmHg or 101 mmHg) by machine learning the relationship between the feature value of the biometric information that is learning data and the blood pressure in advance. Learning discriminating means (for example, discriminant function: evaluation function: discriminant function) is generated. Next, when the blood pressure is actually estimated, the learning identification unit is used to determine a predetermined blood pressure (for example, every 1 mmHg such as 100 mmHg or 101 mmHg) based on the feature amount of the biological information obtained by the measurement. Then, binarization determination is performed respectively. Then, the blood pressure is estimated based on a plurality of determination data obtained as a result.

この様に、本発明では、所定の血圧毎に2値化判定を行い、その複数回のデータに基づいて血圧の推定を行うので、従来の手法に比較して、真値が曖昧な(真値が誤差を含む)血圧に対して、2値化することで、曖昧なデータの影響を低減し、高い精度で血圧を推定できるという効果を奏する。   In this way, in the present invention, binarization determination is performed for each predetermined blood pressure, and blood pressure is estimated based on the data of a plurality of times. Therefore, the true value is ambiguous compared to the conventional method (true By binarizing the blood pressure (which includes an error), the influence of ambiguous data is reduced, and the blood pressure can be estimated with high accuracy.

なお、ここで、特徴量とは、心電信号から得られる信号の状態を示す特徴的な値や脈波信号から得られる信号の状態を示す特徴的な値(従って生体の状態を示す指標となる特徴的な値)、或いは心電信号と脈波信号とから得られる生体の状態を示す指標となる特徴的な値であり、例えば下記の様な周知の特徴量の値を1又は複数採用できる。   Here, the feature amount is a characteristic value indicating a signal state obtained from an electrocardiogram signal or a characteristic value indicating a signal state obtained from a pulse wave signal (accordingly, an index indicating the state of a living body). Or a characteristic value that serves as an index indicating the state of a living body obtained from an electrocardiogram signal and a pulse wave signal. For example, one or a plurality of well-known feature values such as the following are used. it can.

心電信号と脈波信号との時間差(PTT:図3参照)、心拍数(HR)、心拍間隔(RRI)、脈波のAI(Augmentation Index:脈波における反射波(P2)と進行波(P1)の大きさの比(P2/P1):図4参照)、速度脈波(1階微分)の特徴量a1〜f1(図4参照)、加速度脈波(2階微分)の特徴量a〜f(時間軸に沿った波のピークを示すa波(収縮初期陽性波)、b波(収縮初期陰性波)、c波(収縮中期再上昇波)、d波(収縮後期再下降波)、e波(拡張初期陽性波)、f波(拡張初期陰性波):図4参照)など。   Time difference between electrocardiogram signal and pulse wave signal (PTT: see FIG. 3), heart rate (HR), heart rate interval (RRI), pulse wave AI (Augmentation Index: reflected wave (P2) in pulse wave and traveling wave ( P1) magnitude ratio (P2 / P1): see FIG. 4), velocity pulse wave (first-order derivative) feature quantities a1 to f1 (see FIG. 4), acceleration pulse wave (second-order derivative) feature quantity a ˜f (a wave indicating initial wave peak along time axis (contraction early positive wave), b wave (initial contraction negative wave), c wave (contraction middle re-rising wave), d wave (late systolic re-falling wave) E wave (extended initial positive wave), f wave (extended initial negative wave): see FIG.

なお、前記特徴量として、上述した特徴量に加えて、生体情報の特徴量である身体特徴量(身体のある特徴を示す値:例えば身長、体重、性別、年齢)や生体のインピーダンスなどを、前記血圧の推定に用いることもできる。   As the feature amount, in addition to the above-described feature amount, a body feature amount (a value indicating a feature of the body: for example, height, weight, sex, age), a biological impedance, etc. It can also be used to estimate the blood pressure.

(2)請求項2の発明は、前記学習識別手段は、アダブースト(adaboost)により、前記血圧の2値化のための識別関数を生成したことを特徴とする。   (2) The invention of claim 2 is characterized in that the learning identification means generates an identification function for binarizing the blood pressure by adaboost.

本発明では、ブースティングの一種である周知のアダブーストの手法を利用して、血圧の2値化判定を行うものである。   In the present invention, the blood pressure binarization determination is performed using a known Adaboost technique which is a kind of boosting.

このブースティングとは、弱い性能を持った学習機(識別機)を多数組み合わせて高性能の学習機を構成する手法のことであり、ブースティングの最も一般的なアルゴリズムであるアダブーストは、苦手な例題に指数関数的に更新される重みをかけることで学習し、逐次的に能力を高めていくアルゴリズムであり、非常に簡潔な手続きで高い汎用性能を得ることができるため、現在広く用いられている。   This boosting is a method of constructing a high-performance learning machine by combining a number of learning machines (identifiers) with weak performance, and AdaBoost, the most common boosting algorithm, is not good at it. It is an algorithm that learns by applying weights that are updated exponentially to examples, and gradually increases the ability, and because it can obtain high general-purpose performance with a very simple procedure, it is currently widely used Yes.

なお、アダブーストの手法については、特開2008−9843号公報、特開2007−307358号公報、特開2008−140230号公報、特開2008−117391号公報、数理解析研究所講究録 1439巻 2005年 111-127等に開示されている。   As for the method of Adaboost, JP 2008-9843 A, JP 2007-307358 A, JP 2008-140230 A, JP 2008-117391 A, Mathematical Analysis Laboratory, Vol. 1439, 2005. 111-127 and the like.

(3)請求項3の発明は、前記複数回の2値化判定の結果から、各血圧毎の判定結果を並べた2値数列を作成し、この2値数列と前もって準備された符号表における各血圧毎の数列とを対比して、前記血圧の推定を行うことを特徴とする。   (3) The invention of claim 3 creates a binary number sequence in which the determination results for each blood pressure are arranged from the result of the plurality of binarization determinations, and in the code table prepared in advance with this binary number sequence The blood pressure is estimated by comparing with a sequence for each blood pressure.

本発明は、血圧推定の手法(例えばECOC(Error Correcting Output Coding)法:T. G. Dietterich and G. Bakiri., "Solving multiclass learning problems via errorcorrecting output codes", Journal of Artificial Intelligence Research, 2:265-286,1995))参照)を例示したものである。本発明では、複数回の2値化判定の結果から得られた2値数列と符号表の各血圧毎の数列(例えば2値数列)とを対比して、その対応する各桁の一致の程度などから、血圧の推定を行うことができる。なお、両数列の各桁の対比を行うために、各数列の桁数は同じとする。   The present invention relates to a blood pressure estimation method (for example, an ECCC (Error Correcting Output Coding) method: TG Dietterich and G. Bakiri., “Solving multiclass learning problems via error correcting output codes”, Journal of Artificial Intelligence Research, 2: 265-286, 1995))). In the present invention, a binary number sequence obtained from the result of multiple binarization determinations is compared with a number sequence for each blood pressure (for example, a binary number sequence) in the code table, and the degree of coincidence of the corresponding digits. From this, blood pressure can be estimated. In addition, in order to compare each digit of both number sequences, the number of digits of each number sequence is the same.

この符号表とは、例えば図7(a)に示す様に、各血圧毎に設定された数列の所定の桁の数値が、血圧の増加(又は減少)するに従って、徐々に変化(例えば一桁ずつ0から1に増加)するものである。なお、符号表の数列としては、「0」、「1」の数列が考えられるが、重み付けのために、例えば図7(b)に示す様に、それ以外の数値の数列としてもよい。   For example, as shown in FIG. 7A, the code table is a numerical value of a predetermined digit of a numerical sequence set for each blood pressure gradually changes (for example, one digit) as the blood pressure increases (or decreases). Increments from 0 to 1). As the number sequence of the code table, the number sequences of “0” and “1” are conceivable. However, for weighting, for example, as shown in FIG.

(4)請求項4の発明は、前記複数回の2値化判定の結果から得られた2値数列と前記符号表の2値数列とから、両数列における最小ハミング距離を求め、該最小ハミング距離に基づいて、前記血圧の推定を行うことを特徴とする。   (4) According to the invention of claim 4, the minimum Hamming distance in both number sequences is obtained from the binary number sequence obtained from the result of the plurality of binarization determinations and the binary number sequence of the code table, and the minimum Hamming The blood pressure is estimated based on the distance.

本発明は、符号表を用いた具体的な手法を例示したものである。例えば両数列において対応する桁の数値が同じであれば0とし異なっていれば1として、各桁の判定結果を加算してハミング距離を求め、その最小値を最小ハミング距離として求めることができる。   The present invention exemplifies a specific method using a code table. For example, if the numerical values of the corresponding digits in both number sequences are the same, it is 0, and if they are different, it is 1 and the determination result of each digit is added to obtain the Hamming distance, and the minimum value can be obtained as the minimum Hamming distance.

この最小ハミング距離に対応した血圧は、その推定精度が高いと考えられるので、本発明では、最小ハミング距離に基づいて血圧の推定を行う。   Since it is considered that the blood pressure corresponding to the minimum hamming distance has high estimation accuracy, the present invention estimates the blood pressure based on the minimum hamming distance.

(5)請求項5の発明は、前記符号表の数列が、複数個ずつの幅で数値が異なるように設定されていることを特徴とする(実施例3参照)。   (5) The invention of claim 5 is characterized in that the numerical sequences of the code table are set so that the numerical values are different with a plurality of widths (see Example 3).

例えば図11に示す様に、符号表の各数列において各桁の数値が各桁毎に増加する場合に比べて、複数桁(例えば2桁)毎に増加する場合の方が、血圧の推定の際の演算の負担が軽減される(計算時間が短縮される)。   For example, as shown in FIG. 11, in the case where the numerical value of each digit in each numerical sequence of the code table is increased for each digit, the blood pressure estimation is more likely to occur when the digit number is increased by a plurality of digits (for example, two digits). The calculation burden is reduced (calculation time is shortened).

(6)請求項6の発明は、前記符号表として、被験者の状態に適応して設定された数列の符号表を用いることを特徴とする(実施例4参照)。   (6) The invention of claim 6 is characterized in that a code table of a numerical sequence set in accordance with the condition of the subject is used as the code table (see Example 4).

例えば被験者が高齢の場合には、低い血圧である可能性は低いので、その場合は、例えば図12に示す様に、低い血圧を判定する場合の数列のない符号表を用いることができる。これにより、高い血圧の判定しか実施しないので、計算時間が短縮される。   For example, when the subject is elderly, there is a low possibility that the blood pressure is low. In this case, as shown in FIG. 12, for example, a code table without a sequence for determining low blood pressure can be used. As a result, only high blood pressure is determined, so that the calculation time is shortened.

(7)請求項7の発明は、前記血圧の2値化の判定を、所定の血圧範囲の群ごとに行うことを特徴とする(実施例5参照)。   (7) The invention of claim 7 is characterized in that the determination of binarization of the blood pressure is performed for each group of a predetermined blood pressure range (see Example 5).

学習において真値となる市販のカフ式血圧計には、通常4〜8mmHgの誤差があると言われている。従って、例えば1mmHgの様な細かい血圧毎に2値化判定を行うのではなく、例えば4〜8mmHgの群(ステップ)毎に2値化判定を行ってもよい。これによっても、ある程度の測定精度が得られ、且つ、計算時間が短縮される。   It is said that a commercially available cuff sphygmomanometer that is a true value in learning usually has an error of 4 to 8 mmHg. Therefore, for example, instead of performing binarization determination for each fine blood pressure such as 1 mmHg, binarization determination may be performed for each group (step) of 4 to 8 mmHg, for example. This also provides a certain degree of measurement accuracy and shortens the calculation time.

(8)請求項8の発明は、前記所定の血圧の範囲では前記群の幅を狭くし、前記所定の血圧の範囲外では前記群の幅を広くすることを特徴とする(実施例6参照)。   (8) The invention of claim 8 is characterized in that the width of the group is narrowed within the range of the predetermined blood pressure, and the width of the group is widened outside the range of the predetermined blood pressure (see Example 6). ).

学習データは、平均的な血圧の範囲(データの中央部分:例えば100〜140mmHg)に集中しており、この部分のデータ数は多く、それ以外の部分のデータ数は少ない。   The learning data is concentrated in an average blood pressure range (central portion of data: for example, 100 to 140 mmHg), and the number of data in this portion is large, and the number of data in other portions is small.

従って、学習データの数の多い領域の血圧判定の幅を狭くし(例えば1mmHg)、学習データの数の少ない領域の血圧判定の幅を広くする(例えば4〜8mmHg)。   Therefore, the range of blood pressure determination in a region with a large amount of learning data is narrowed (for example, 1 mmHg), and the range of blood pressure determination in a region with a small amount of learning data is widened (for example, 4 to 8 mmHg).

これにより、計算時間の短縮とリファレンス血圧(学習データ)の曖昧性を低減することができる。   Thereby, shortening of calculation time and ambiguity of reference blood pressure (learning data) can be reduced.

なお、所定の血圧の範囲を、100〜140mmHgに設定することが、測定精度を確保しつつ、計算時間を短縮できるので好ましい。また、群の幅を狭くする際には、2mmHg以上としてもよいが、最小単位の1mmHgとしてもよい。   In addition, it is preferable to set the predetermined blood pressure range to 100 to 140 mmHg because the calculation time can be shortened while ensuring the measurement accuracy. Further, when the width of the group is narrowed, it may be 2 mmHg or more, but may be 1 mmHg as a minimum unit.

(9)請求項9の発明は、前記学習識別手段における事前学習の際に、所定の血圧の範囲の学習データを学習に用いないことを特徴とする(実施例7参照)。   (9) The invention of claim 9 is characterized in that learning data in a predetermined blood pressure range is not used for learning in the prior learning in the learning identification means (see Example 7).

市販のカフ血圧計には、通常4〜8mmHgの誤差があるので、その点を考慮して、二値化の学習データの血圧値から、所定幅以内(例えば±4mmHg以内など)のデータを用いないようにしたものである。これにより、血圧の推定精度が向上する。   Since a commercially available cuff sphygmomanometer usually has an error of 4 to 8 mmHg, in consideration of this, data within a predetermined range (for example, within ± 4 mmHg) is used from the blood pressure value of binarized learning data. It is not to be. This improves blood pressure estimation accuracy.

(10)請求項10の発明は、前記学習識別手段における事前学習の際に、異なる血圧の範囲における二値化のデータ数の比を、所定の範囲内に設定することを特徴とする(実施例8参照)。   (10) The invention of claim 10 is characterized in that the ratio of the number of binarized data in a range of different blood pressures is set within a predetermined range at the time of pre-learning in the learning discriminating means. (See Example 8).

各血圧層における学習データ数が大きく異なる場合には、その脈波が本来属する血圧に同じような特徴量を持つ学習データがあったとしても、その脈波が本来属しない血圧にも同じような特徴量を持つ学習データが存在する確率が高くなり、結果として、血圧の推定精度が低下することがある。   If the number of learning data in each blood pressure layer is significantly different, even if there is learning data with the same feature amount in the blood pressure to which the pulse wave originally belongs, the same applies to blood pressure to which the pulse wave does not originally belong. The probability that there exists learning data having a feature amount increases, and as a result, blood pressure estimation accuracy may decrease.

従って、本発明では、2値化判定のために所定の血圧値で2つに区分された各血圧層(例えば140mmHg未満/以上に区分された各血圧層)における学習データのデータ数の比を、所定範囲内(例えば2倍以内)に設定する。これにより、本来属しない血圧に同じような特徴量を持つデータが存在する確率が低くなるため、血圧の測定精度が向上する。   Therefore, in the present invention, the ratio of the number of learning data in each blood pressure layer (for example, each blood pressure layer divided into less than 140 mmHg) divided into two by a predetermined blood pressure value for binarization determination is calculated. , Set within a predetermined range (for example, within 2 times). As a result, the probability that there is data having a similar feature amount in blood pressure that does not originally belong is reduced, and blood pressure measurement accuracy is improved.

(11)請求項11の発明は、前記血圧を推定する際に、データの信頼度を算出し、血圧推定結果とともに、前記データの信頼度を表示することを特徴とする。   (11) The invention of claim 11 is characterized in that when the blood pressure is estimated, the reliability of the data is calculated, and the reliability of the data is displayed together with the blood pressure estimation result.

学習データの数などによって、血圧の推定値の信頼度(信頼性の程度)が変化する。そのため、信頼度を併せて結果を表示することで、計測者や被計測者に正確な情報を伝えることができる。   The reliability (degree of reliability) of the estimated value of blood pressure varies depending on the number of learning data and the like. Therefore, accurate information can be conveyed to the measurer and the person to be measured by displaying the result together with the reliability.

(12)請求項12の発明は、前記データの信頼度を算出し、この信頼度が所定の判定値より低い場合には、より信頼度が高いと予想される異なる手法を採用して血圧の推定を行うことを特徴とする。   (12) The invention of claim 12 calculates the reliability of the data, and when the reliability is lower than a predetermined determination value, adopts a different method that is expected to have higher reliability, It is characterized by performing estimation.

例えば学習データの数に応じて信頼度を設定する。例えば学習データの数が所定値以上の場合には信頼度高と設定し、所定値未満の場合には信頼度低と設定する。   For example, the reliability is set according to the number of learning data. For example, when the number of learning data is greater than or equal to a predetermined value, the reliability is set high, and when the number is less than the predetermined value, the reliability is set low.

従って、この信頼度に応じて、異なる手法を採用するか否かを判定することができる。   Therefore, it is possible to determine whether or not to adopt a different method according to the reliability.

(13)請求項13の発明は、前記データの信頼度は、学習データの数又は最小ハミング距離に基づいて判断することを特徴とする。   (13) The invention of claim 13 is characterized in that the reliability of the data is determined based on the number of learning data or the minimum Hamming distance.

本発明は、データの信頼度の根拠を例示したものである。例えば学習データが少ない場合には、学習データが不十分なため血圧の推定値の信頼度が低く、また、最小ハミング距離が小さいほど、符号表との一致度が高く、血圧の信頼度は高いと考えられる。従って、学習データの数や最小ハミング距離に基づいて、血圧の信頼度を判断することができる。   The present invention exemplifies the basis of data reliability. For example, when the learning data is small, the reliability of the estimated blood pressure value is low because the learning data is insufficient, and the smaller the minimum Hamming distance, the higher the degree of coincidence with the code table and the higher the blood pressure reliability. it is conceivable that. Therefore, the reliability of blood pressure can be determined based on the number of learning data and the minimum Hamming distance.

(14)請求項14の発明は、前記異なる手法は、重回帰による推定の手法であることを特徴とする。   (14) The invention of claim 14 is characterized in that the different method is a method of estimation by multiple regression.

本発明は、異なる手法を例示したものである。学習データの数が少ない場合の血圧の推定値は、重回帰による血圧の推定値より精度が低いと考えられる。つまり、重回帰では、すべてのデータを用いて、ある値を外挿し推定するため、データ数が少ない領域でもある程度の精度が確保できるが、学習では、その1点の値を推定するため、データ数が少ない領域では精度が低下する。そのため、学習データの数が所定値より少ない場合は、重回帰による手法を採用することにより、血圧推定の精度が向上する。   The present invention illustrates different approaches. The estimated value of blood pressure when the number of learning data is small is considered to be less accurate than the estimated value of blood pressure by multiple regression. In other words, in multiple regression, a certain value is extrapolated and estimated using all data, so that a certain degree of accuracy can be ensured even in an area where the number of data is small. The accuracy decreases in a region where the number is small. Therefore, when the number of learning data is less than a predetermined value, the accuracy of blood pressure estimation is improved by adopting a technique based on multiple regression.

(15)請求項15の発明は、生体から検出された心電信号及び脈波信号のうち、少なくとも前記脈波信号から得られる生体情報を用いて、血圧を推定する血圧測定装置において、前記生体情報の特徴量に基づいて、重回帰によって血圧を推定する第1血圧推定手段と、前記生体情報の特徴量と血圧との関係の学習によって得られた判別手段を用い、前記特徴量から血圧を推定する第2血圧推定手段と、前記第1血圧推定手段によって推定された第1推定血圧と、前記第2血圧推定手段によって推定された第2推定血圧とを、所定の信頼度の高低を判定する判定条件に基づいて比較して、前記信頼度が高い推定血圧を決定する血圧決定手段と、を備えたことを特徴とする。   (15) The invention of claim 15 is the blood pressure measurement apparatus for estimating blood pressure using at least biological information obtained from the pulse wave signal among the electrocardiogram signal and the pulse wave signal detected from the living body. Based on the feature quantity of the information, the first blood pressure estimation means for estimating the blood pressure by multiple regression and the discrimination means obtained by learning the relationship between the feature quantity of the biological information and the blood pressure, the blood pressure is calculated from the feature quantity. The second blood pressure estimating means to be estimated, the first estimated blood pressure estimated by the first blood pressure estimating means, and the second estimated blood pressure estimated by the second blood pressure estimating means are used to determine whether the predetermined reliability is high or low. Blood pressure determining means for determining an estimated blood pressure having a high degree of reliability by comparison based on the determination condition.

本発明では、生体情報の特徴量に基づいて、重回帰によって血圧を推定する第1血圧推定手段と、生体情報の特徴量と血圧との関係の学習によって得られた判定手段(例えば前記学習識別手段)を用い、特徴量から血圧を推定する第2血圧推定手段とを備えているので、信頼度の高い方を判定により採用することにより、精度の高い血圧推定が可能となる。   In the present invention, a first blood pressure estimation unit that estimates blood pressure by multiple regression based on a feature amount of biological information, and a determination unit obtained by learning a relationship between the feature amount of biological information and blood pressure (for example, the learning identification) And the second blood pressure estimating means for estimating the blood pressure from the feature amount, the highly reliable blood pressure can be estimated by adopting the one with higher reliability by the determination.

なお、第2血圧推定手段としては、アダブーストによる学習識別手段以外に、前記特許文献11、12等に記載のニューラルネットワーク(NN:Neural Network)や、SVM(Support Vector Machine)等の他のマシンラーンイングの手法により血圧を推定する手段を採用することができる。   As the second blood pressure estimation means, in addition to the learning identification means by Adaboost, other machine learners such as a neural network (NN: Neural Network) described in Patent Documents 11 and 12, the SVM (Support Vector Machine), etc. A means for estimating blood pressure by the Ing technique can be employed.

また、所定の信頼度の高低を判定する判定条件としては、学習データの数、血圧の範囲、脈波特徴量の範囲、(請求項1の手法による)最小ハミング距離などを採用できる。   As the determination conditions for determining the level of the predetermined reliability, the number of learning data, the blood pressure range, the pulse wave feature amount range, the minimum Hamming distance (by the method of claim 1), and the like can be adopted.

例えば学習データが少ない場合には、推定血圧の信頼度が低いと考えられるので、学習データの数が所定値より少ない場合には、重回帰によって求めた血圧を採用することができる。   For example, when the learning data is small, it is considered that the reliability of the estimated blood pressure is low. Therefore, when the number of learning data is less than a predetermined value, the blood pressure obtained by multiple regression can be employed.

なお、前記特徴量として、上述した心電信号や脈波信号から得られる特徴量に加えて、生体情報の特徴量である身体特徴量(身体のある特徴を示す値:例えば身長、体重、性別、年齢)や生体のインピーダンスなどを、前記血圧の推定に用いることもできる。   As the feature amount, in addition to the feature amount obtained from the above-described electrocardiogram signal or pulse wave signal, a body feature amount that is a feature amount of biological information (a value indicating a feature of the body: for example, height, weight, sex, etc. , Age), biological impedance, and the like can also be used for the estimation of the blood pressure.

(16)請求項16の発明は、前記第2血圧推定手段は、学習データである前記生体情報の特徴量と血圧との関係を予め機械学習することにより得られた識別手段であって、所定の血圧毎に前記生体情報の特徴量から血圧を2値化判別する学習識別手段と、前記血圧の推定を行う場合に、前記学習識別手段を用い、測定によって得られた前記生体情報の特徴量に基づいて、前記所定の血圧に対して2値化判別を行う2値化判別手段と、前記2値化判別手段によって、前記生体情報の特徴量に基づいて異なる前記所定の血圧に対してそれぞれ2値化判定を行い、該複数回の2値化判定の結果に基づいて前記血圧を推定する血圧推定手段と、を備えたことを特徴とする。   (16) In the invention of claim 16, the second blood pressure estimation means is an identification means obtained by machine learning in advance of the relationship between the characteristic amount of the biological information, which is learning data, and blood pressure, Learning identification means for binarizing and discriminating blood pressure from the feature amount of the biological information for each blood pressure, and the feature amount of the biological information obtained by measurement using the learning identification means when estimating the blood pressure And binarization determination means for performing binarization determination on the predetermined blood pressure, and the binarization determination means for each of the predetermined blood pressures that differ based on the feature amount of the biological information. Blood pressure estimation means for performing binarization determination and estimating the blood pressure based on the result of the plurality of binarization determinations.

本発明は、前記第2血圧推定手段を例示したものである。従って、第2血圧推定手段として、前記請求項1〜13の発明を採用することができる。   The present invention exemplifies the second blood pressure estimating means. Therefore, the invention of claims 1 to 13 can be employed as the second blood pressure estimating means.

(17)請求項17の発明は、前記重回帰により求めた第1推定血圧が、100〜140mmHgの場合には、血圧の推定値として、前記学習による第2推定血圧を採用することを特徴とする。   (17) The invention of claim 17 is characterized in that, when the first estimated blood pressure obtained by the multiple regression is 100 to 140 mmHg, the second estimated blood pressure by the learning is adopted as an estimated value of blood pressure. To do.

100〜140mmHgの血圧は、通常の血圧の範囲であり、学習データ数も多いと考えられるので、この場合は学習による血圧の推定値を採用するものである。   The blood pressure of 100 to 140 mmHg is a range of normal blood pressure, and it is considered that the number of learning data is large. In this case, an estimated value of blood pressure by learning is adopted.

(18)請求項18の発明は、前記学習データが、1mmHg当たり5個以上の場合には、血圧の推定値として、前記学習による第2推定血圧を採用することを特徴とする。   (18) The invention of claim 18 is characterized in that when the learning data is 5 or more per 1 mmHg, the second estimated blood pressure by the learning is adopted as an estimated value of blood pressure.

本発明は、学習データの数の好ましい範囲を例示したものである。   The present invention exemplifies a preferable range of the number of learning data.

(19)請求項19の発明は、前記請求項1〜18のいずれかに記載の血圧測定装置の機能を実現させるためのプログラムである。   (19) The invention of claim 19 is a program for realizing the function of the blood pressure measurement device according to any one of claims 1 to 18.

よって、マイクロコンピュータ等でこのプログラムを実施することにより、上述した生体情報などを用いて、血圧を推定することが可能である。   Therefore, by executing this program with a microcomputer or the like, it is possible to estimate the blood pressure using the biological information described above.

(20)請求項20の発明は、前記請求項19に記載のプログラムを記録したコンピュータ読取可能な記録媒体である。   (20) A twentieth aspect of the invention is a computer-readable recording medium on which the program according to the nineteenth aspect is recorded.

よって、この記録媒体に記録したプログラムを用いて、上述したように、血圧を推定することができる。   Therefore, using the program recorded on this recording medium, the blood pressure can be estimated as described above.

車両内における血圧測定装置のシステム構成を示す説明図である。It is explanatory drawing which shows the system configuration | structure of the blood-pressure measurement apparatus in a vehicle. 血圧測定装置等の機能を示すブロック図であるIt is a block diagram which shows functions, such as a blood pressure measuring device 心電信号と脈波信号を示すグラフである。It is a graph which shows an electrocardiogram signal and a pulse wave signal. 脈波信号とそれを微分した信号を示すグラフである。It is a graph which shows a pulse wave signal and the signal which differentiated it. 実施例1における血圧推定の手順を示す説明図である。It is explanatory drawing which shows the procedure of the blood-pressure estimation in Example 1. FIG. 血圧の2値化判別の手法を示す説明図である。It is explanatory drawing which shows the method of the binarization discrimination | determination of a blood pressure. (a)は重み付け無しの符号表を示す説明図、(b)は重み付け有りの符号表を示す説明図である。(A) is explanatory drawing which shows a code table without weighting, (b) is explanatory drawing which shows a code table with weighting. 実施例1における血圧推定の処理を示すフローチャートである。3 is a flowchart illustrating a blood pressure estimation process according to the first embodiment. 実施例2における血圧推定の手順を示す説明図である。It is explanatory drawing which shows the procedure of the blood-pressure estimation in Example 2. FIG. 実施例2における血圧推定の処理を示すフローチャートである。10 is a flowchart illustrating a blood pressure estimation process in the second embodiment. 実施例3において用いる符号表等を示す説明図である。It is explanatory drawing which shows the code | cord | chord table etc. which are used in Example 3. FIG. 実施例4において用いる符号表等を示す説明図である。It is explanatory drawing which shows the code | cord | chord table etc. which are used in Example 4. FIG. 実験例1(アダブーストと重回帰の比較)の実験結果を示す説明図である。It is explanatory drawing which shows the experimental result of Experimental example 1 (comparison of AdaBoost and multiple regression). 実験例2(各血圧層の推定精度の比較)の実験結果を示す説明図である。It is explanatory drawing which shows the experimental result of Experimental example 2 (comparison of the estimation precision of each blood-pressure layer). 実験例3(群による血圧推定)の実験結果を示す説明図である。It is explanatory drawing which shows the experimental result of Experimental example 3 (blood pressure estimation by a group). 実験例4(最小ハミング距離と誤差)の実験結果を示す説明図である。It is explanatory drawing which shows the experimental result of Experimental example 4 (minimum hamming distance and error). 実験例5(学習データの選別の効果)の実験結果を示す説明図である。It is explanatory drawing which shows the experimental result of Experimental example 5 (effect of selection of learning data). 実験例6(学習データ比の影響)の実験結果を示す説明図である。It is explanatory drawing which shows the experimental result of Experimental example 6 (effect of learning data ratio).

次に、本発明の血圧測定装置の実施例について説明する。   Next, an embodiment of the blood pressure measurement device of the present invention will be described.

本実施例では、自動車に搭載されてドライバの血圧を測定する血圧測定装置について説明する。   In the present embodiment, a blood pressure measurement device that is mounted on an automobile and measures the blood pressure of a driver will be described.

a)まず、本実施例の血圧測定装置のシステム構成を説明する。   a) First, the system configuration of the blood pressure measurement device according to the present embodiment will be described.

図1に示す様に、本実施例では、自動車に搭載された血圧測定装置1と、ドライバ等に情報の報知を行う報知装置3と、手動にてデータの入力を行うマニュアル入力部5と、ステアリング7に取り付けられた脈波センサ9と、ステアリング7に取り付けられた一対の電極11、13とを備えている。   As shown in FIG. 1, in this embodiment, a blood pressure measurement device 1 mounted on an automobile, a notification device 3 that notifies a driver or the like, a manual input unit 5 that manually inputs data, A pulse wave sensor 9 attached to the steering 7 and a pair of electrodes 11 and 13 attached to the steering 7 are provided.

前記血圧測定装置1は、周知のマイクロコンピュータを中心とした電子制御装置であり、脈波センサ9や電極11、13等からの信号に基づいて、血圧の算出や報知装置3の制御等を行う。   The blood pressure measurement device 1 is an electronic control device centered on a well-known microcomputer, and calculates blood pressure, controls the notification device 3, and the like based on signals from the pulse wave sensor 9, the electrodes 11, 13, and the like. .

前記報知装置3は、血圧などを表示する液晶等のディスプレイ15や、その内容を音声等で出力するスピーカ17により構成される。   The notification device 3 includes a display 15 such as a liquid crystal that displays blood pressure and the like, and a speaker 17 that outputs the contents thereof by voice or the like.

前記マニュアル入力部5は、マニュアルにて、体重、身長、年齢、性別等の個人別の身体特徴情報が入力可能な例えばキーボードやテンキーやリモコンなどの入力装置である。なお、ディスプレイ15の表示画面をタッチパネルとして、データを入力するようにしてもよい。   The manual input unit 5 is an input device such as a keyboard, a numeric keypad, or a remote controller capable of inputting individual body characteristic information such as weight, height, age, and sex manually. Note that data may be input using the display screen of the display 15 as a touch panel.

前記脈波センサ9は、周知の発光素子(LED)や受光素子(PD)を備えた光学式のセンサであり、例えばドライバの指先等に光を照射し、その反射波を利用して脈波(容積脈波)を検出することができる。従って、後述する血圧の算出に用いる脈波信号は、この脈波センサ9から得ることができる。   The pulse wave sensor 9 is an optical sensor provided with a known light emitting element (LED) or light receiving element (PD). For example, the pulse wave sensor 9 irradiates light on a fingertip of a driver and uses the reflected wave to make a pulse wave. (Volume pulse wave) can be detected. Therefore, a pulse wave signal used for blood pressure calculation described later can be obtained from the pulse wave sensor 9.

前記一対の電極11、13は、ドライバの左右の手にそれぞれ接触するように、ステアリング9の左右に配置されている。この電極11、13は、心電信号を得るために電圧を印加する心電計の電極として用いられる。ここでは、両電極11、13と血圧測定装置1が、心電計の機能を果たす。従って、後述する血圧の算出に用いる心電信号は、この電極11、13を利用して得ることができる。   The pair of electrodes 11 and 13 are arranged on the left and right of the steering wheel 9 so as to come into contact with the left and right hands of the driver, respectively. These electrodes 11 and 13 are used as electrodes of an electrocardiograph that applies a voltage to obtain an electrocardiographic signal. Here, both the electrodes 11 and 13 and the blood pressure measuring device 1 serve as an electrocardiograph. Therefore, an electrocardiographic signal used for blood pressure calculation described later can be obtained using the electrodes 11 and 13.

b)次に、本実施例の血圧測定装置1等の機能をブロック図にて更に詳しく説明する。   b) Next, the function of the blood pressure measurement device 1 and the like of this embodiment will be described in more detail with reference to a block diagram.

図2に示す様に、本実施例の血圧測定装置1には、心電信号取得部25と、脈波信号取得部27と、心電信号解析部31と、心電及び脈波信号解析部33と、脈波信号解析部35と、血圧演算部37とを備えている。   As shown in FIG. 2, the blood pressure measurement device 1 of the present embodiment includes an electrocardiogram signal acquisition unit 25, a pulse wave signal acquisition unit 27, an electrocardiogram signal analysis unit 31, and an electrocardiogram and pulse wave signal analysis unit. 33, a pulse wave signal analysis unit 35, and a blood pressure calculation unit 37.

このうち、心電信号取得部25は、心臓の活動に伴う電気的興奮を、電極11、13間の電位差(心電信号)として計測するものである。   Among these, the electrocardiogram signal acquisition unit 25 measures electrical excitement accompanying the activity of the heart as a potential difference (electrocardiogram signal) between the electrodes 11 and 13.

脈波信号取得部27は脈波センサ9を駆動して脈波信号を取得するものである。   The pulse wave signal acquisition unit 27 drives the pulse wave sensor 9 to acquire a pulse wave signal.

心電信号解析部31は、心電信号を解析して、例えば心拍間隔(HR)等を算出するものである。   The electrocardiogram signal analysis unit 31 analyzes the electrocardiogram signal and calculates, for example, a heartbeat interval (HR).

心電及び脈波信号解析部33は、図3に示す様に、心電信号と脈波信号を用いて、心電信号に対する脈波信号の遅れ時間である脈波伝播時間(PTT)を求めるものである。   As shown in FIG. 3, the electrocardiogram and pulse wave signal analysis unit 33 uses the electrocardiogram signal and the pulse wave signal to obtain a pulse wave propagation time (PTT) that is a delay time of the pulse wave signal with respect to the electrocardiogram signal. Is.

脈波信号解析部35は、図4に示す様に、脈波信号を解析して、1階微分(速度脈波)、2階微分(加速度脈波)、3階微分、4階微分を行うとともに、各微分における特徴量(例えば速度脈波のa1〜k1、加速度脈波のa〜f等)の算出、AI(容積AI)の算出などを行うものである。   As shown in FIG. 4, the pulse wave signal analysis unit 35 analyzes the pulse wave signal and performs first-order differentiation (velocity pulse wave), second-order differentiation (acceleration pulse wave), third-order differentiation, and fourth-order differentiation. At the same time, feature amounts (for example, velocity pulse waves a1 to k1, acceleration pulse waves a to f, etc.) and AI (volume AI) are calculated.

血圧演算部37は、生体情報である脈波信号及び心電信号のうち少なくとも脈波信号から得られる特徴量、例えばPTT、HR、速度脈波や加速度脈波などの特徴量、容積AI等を用い、後述する学習や重回帰などによって血圧を推定(算出)するものである。   The blood pressure calculation unit 37 obtains at least a feature amount obtained from the pulse wave signal among the pulse wave signal and the electrocardiogram signal as biological information, for example, a feature amount such as PTT, HR, velocity pulse wave, acceleration pulse wave, volume AI, and the like. The blood pressure is estimated (calculated) by learning, multiple regression, etc., which will be described later.

なお、この血圧演算部37には、後述する様に、予め学習データによって(所定血圧毎に)血圧の2値化判別を行うための識別関数が生成されて記憶されている。   As will be described later, an identification function for performing binarization determination of the blood pressure is generated and stored in advance in the blood pressure calculation unit 37 (for each predetermined blood pressure).

b)次に、本実施例における血圧測定の手順について説明する。   b) Next, the procedure of blood pressure measurement in the present embodiment will be described.

本実施例では、ブースティングの一種である周知のアダブーストの手法を利用して、血圧の2値化判定を行うものである。   In the present embodiment, the blood pressure binarization is determined using a well-known Adaboost technique which is a kind of boosting.

(1)まず、アダブーストにおける事前学習モードについて説明する。   (1) First, the pre-learning mode in AdaBoost will be described.

ここでは、アダブーストの手法により血圧の2値化判別を行う識別関数を作成する。つまり、血圧が既知である場合の特徴量のデータを予め学習することにより、(実際の血圧の測定に際して得られた特徴量が所定の血圧未満か以上かを判別するための)識別関数を作成する。   Here, a discriminant function for performing binarization determination of blood pressure is created by the Adaboost method. In other words, an identification function is created (to determine whether the feature value obtained during actual blood pressure measurement is less than or equal to the predetermined blood pressure) by learning in advance the feature value data when the blood pressure is known To do.

このアダブーストによって識別関数を生成する手順としては、以下に述べる様に、例えば、数理解析研究所講究録 1439巻 2005年 111-127等に開示されている。   As described below, the procedure for generating the discriminant function by Adaboost is disclosed in, for example, Research Institute of Mathematical Analysis Vol. 1439 2005 111-127.

以下、識別関数(F)を生成する方法について詳細に説明する。   Hereinafter, a method for generating the discriminant function (F) will be described in detail.

まず、血圧が既知であるサンプルデータから、予めポジティブサンプル又はネガティブサンプルであることが正解付けされた学習サンプル(x1,y1),・・・,(xN,yN)を用意する。   First, learning samples (x1, y1),..., (XN, yN) that are correctly identified as positive samples or negative samples from sample data with known blood pressure are prepared.

例えばある血圧(例えば100mmHg)を境界として、100未満と100以上とを区分し、100未満の場合を誤とし100以上を正とした場合を考えると、予めポジティブサンプル(正解)とネガティブサンプル(不正解)であることが正解付けされた学習サンプル(x1、y2)・・・(xN、yN)を用意し、このサンプルデータを用いて学習を行う。   For example, if a blood pressure (for example, 100 mmHg) is used as a boundary and less than 100 is divided into 100 or more, and a case where the blood pressure is less than 100 is mistaken and 100 or more is considered positive, a positive sample (correct answer) and a negative sample (incorrect Learning samples (x1, y2)... (XN, yN) with correct answers are prepared and learning is performed using the sample data.

なお、xiは学習サンプルのデータ、すなわち、上述した複数の特徴量を示す。また、yiは学習サンプルのラベルであり、yi=1は、学習サンプルがポジティブサンプル(血圧が所定値以上であることを示すサンプル:正)としてラベリングされている場合を示し、yi=−1は、学習サンプルがネガティブサンプル(例えば血圧が所定値以上ではない(即ち所定値未満):誤)であることを示すサンプルとしてラベリングされていることを示す。   Note that xi represents learning sample data, that is, the plurality of feature amounts described above. Further, yi is a label of the learning sample, yi = 1 indicates a case where the learning sample is labeled as a positive sample (a sample indicating that the blood pressure is equal to or higher than a predetermined value: positive), and yi = −1 , Indicates that the learning sample is labeled as a sample indicating that it is a negative sample (for example, blood pressure is not greater than or equal to a predetermined value (ie, less than a predetermined value): false).

ここで、定義関数を下記式(1)とした場合、アダブーストのアルゴリズムは、以下1〜3の様に書ける。   Here, when the defining function is the following formula (1), the Adaboost algorithm can be written as follows:

1.重みの初期値を、w1(i)=1/N、識別関数をF0(x)=0とする。 1. The initial value of the weight is w 1 (i) = 1 / N, and the discriminant function is F 0 (x) = 0.

2.For t=1・・・T
(a)重み付き誤り率を下記式(2)とした場合、重み付き誤り率を最小にする弱各 学習機(下記式(3))を選択する。
2. For t = 1 ・ ・ ・ T
(A) When the weighted error rate is expressed by the following equation (2), each weak learning machine (the following equation (3)) that minimizes the weighted error rate is selected.

(b)選択した弱学習機に対する係数(下記式(4))を計算する。 (B) A coefficient (the following formula (4)) for the selected weak learning machine is calculated.

この学習機を、Ft=Ft+1+αttで更新する。 This learning machine is updated with F t = F t + 1 + α t f t.

(c)重みを下記式(5)の様に更新する。     (C) The weight is updated as shown in the following formula (5).

3.下記式(6)の識別関数の符号で判別する。 3. It discriminate | determines with the code | symbol of the discriminant function of following formula (6).

従って、前記識別関数が正か負かを判断することで、ポジティブサンプルであるかどうかを評価することができる。 Therefore, by determining whether the discrimination function is positive or negative, it is possible to evaluate whether it is a positive sample.

なお、ここでは、アダブースティングアルゴリズムを用いて識別関数Fを取得する場合について例示しているが、それに代えてニューラルネットワーク(Neural network)、SVM(Support Vector Machine)等他のマシンラーニングの手法を用いることもできる。   Here, the case where the discriminant function F is acquired by using the Adaboosting algorithm is illustrated, but instead, other machine learning methods such as a neural network and SVM (Support Vector Machine) are used. It can also be used.

(2)次に、実際に血圧の測定を行う検知モードにおける手順について説明する。   (2) Next, the procedure in the detection mode for actually measuring blood pressure will be described.

・図5及び図6に示す様に、脈波センサ9及び電極11、13からの信号から、心電信号及び脈波信号の特徴量を求める。   As shown in FIG. 5 and FIG. 6, feature quantities of the electrocardiogram signal and the pulse wave signal are obtained from the signals from the pulse wave sensor 9 and the electrodes 11 and 13.

具体的には、血圧の算出に用いるために、上述した特徴量(PTT、HR、速度脈波の特徴量a1〜k1、加速度脈波の特徴量a〜f、AI(容積AI)等)を算出する。   Specifically, the above-described feature amounts (PTT, HR, velocity pulse wave feature amounts a1 to k1, acceleration pulse wave feature amounts a to f, AI (volume AI), etc.) are used to calculate blood pressure. calculate.

・次に、アダブーストの検知モードにおける2値化判定を複数回実施する。   Next, the binarization determination in the AdaBoost detection mode is performed a plurality of times.

具体的には、事前学習モードで生成した識別関数、即ち、所定の血圧(例えば80mmHg)未満か以上かを2値化判別するため識別関数を用い、その識別関数のパラメータである特徴量から、血圧が80mmHg未満か以上かを2値化判別する。その結果、例えば80mmHg未満であれば判定結果を0と設定し、80mmHg以上であれば判定結果を1と設定する。   Specifically, the discriminant function generated in the pre-learning mode, that is, the discriminant function is used to binarize and discriminate whether the blood pressure is less than or equal to a predetermined blood pressure (for example, 80 mmHg), and from the feature quantity that is the parameter of the discriminant function, Whether the blood pressure is less than 80 mmHg or not is binarized. As a result, for example, if it is less than 80 mmHg, the determination result is set to 0, and if it is 80 mmHg or more, the determination result is set to 1.

同様な判定を、事前学習モードによって生成された(1mmHg毎に2値化判定を行う)各識別関数を用いて、1mmHg毎に判定血圧を徐々に上げて行う。なお、ここでは、血圧は80〜160mmHgの範囲で判定を行う。   Similar determination is performed by gradually increasing the determination blood pressure for each 1 mmHg using each discriminant function generated in the pre-learning mode (performing binarization determination for each 1 mmHg). Here, the blood pressure is determined in the range of 80 to 160 mmHg.

そして、各識別関数による判定結果を、血圧の低い順に配列して2値数列(例えば(0・・011・・1))を作成する。   Then, the determination result by each discrimination function is arranged in order of increasing blood pressure to create a binary number sequence (for example, (0 ·· 011 ·· 1)).

即ち、2値数列の左端の桁が80mmHgの判定結果を示し、右端の桁が160mmHgの判定結果を示す様に、右の桁にゆくほど判定血圧が高くなるようにして、2値数列を作成する。   In other words, a binary number sequence is created so that the determined blood pressure increases as it goes to the right digit so that the leftmost digit of the binary number sequence indicates the determination result of 80 mmHg and the rightmost digit indicates the determination result of 160 mmHg. To do.

・次に、ECOC法により、血圧を推定する。   Next, blood pressure is estimated by the ECOC method.

具体的には、まず、前記数列を、予め1mmHg毎に作成された符号表と対比する。通常は、図6(b)に示す様に、1回判定する(1mmHg毎の81回の判定を1回の判定とする)。   Specifically, first, the numerical sequence is compared with a code table created in advance for each 1 mmHg. Usually, as shown in FIG.6 (b), it determines once (81 determinations for every 1 mmHg are set as 1 determination).

この符号表は、前記判定によって得られる2値数列と同じ桁数の2値数列が、各判定血圧毎に順次配置されたものであり、その2値数列は、前記判定血圧に合わせて例えば80〜160の範囲で1mmHg毎に設定されている。   In this code table, a binary number sequence having the same number of digits as the binary number sequence obtained by the determination is sequentially arranged for each determination blood pressure, and the binary number sequence is, for example, 80 according to the determination blood pressure. It is set for every 1 mmHg in the range of ~ 160.

詳しくは、符号表では、例えば血圧が160から80に行くに従って、図7(a)に示す様に(尚、一部桁数を省略してある)、(0・・0001)、(0・・0011)、(0・・0111)と、2値数列の数字の1を示す桁数が、右側から順に増加するように設定されている。   Specifically, in the code table, for example, as the blood pressure goes from 160 to 80, as shown in FIG. 7A (note that some digits are omitted), (0... 0001), (0. (0011), (0... 0111), and the number of digits indicating the number 1 in the binary number sequence is set to increase sequentially from the right side.

・次に、各血圧毎に、前記判定によって求めた2値数列と符号表の2値数列の各桁の数字を比較する。具体的には、各桁の数値が一致している場合には0とし、不一致の場合は1として、各数字を加算し、いわゆるハミング距離を求める。   Next, for each blood pressure, the binary number sequence obtained by the determination is compared with the numbers in each digit of the binary number sequence of the code table. Specifically, 0 is set when the numerical values of the respective digits match, and 1 is set when the numbers do not match, thereby obtaining a so-called Hamming distance.

そして、各血圧毎にハミング距離を求め、その最小の値(最小ハミング距離)を求め、この最小ハミング距離となる血圧を選択する。ここで、最小ハミング距離となる血圧を選択する理由は、血圧推定の誤差が少ないと考えられるからである。   Then, a hamming distance is obtained for each blood pressure, a minimum value (minimum hamming distance) is obtained, and a blood pressure that becomes the minimum hamming distance is selected. Here, the reason for selecting the blood pressure that is the minimum Hamming distance is that it is considered that there is little error in blood pressure estimation.

つまり、ここでは、最小ハミング距離となる血圧を、最も確からしい血圧として選択するのである。   That is, here, the blood pressure that provides the minimum hamming distance is selected as the most probable blood pressure.

・次に、信頼度を計算する。   -Next, the reliability is calculated.

例えば最小ハミング距離と学習データ数をパラメータとし、例えば選択された血圧に関する学習データが所定値以上十分にある場合には、信頼度が高いとする。例えば学習データ数が多いほど信頼度が高いとし、また、学習データが所定の判定値以下の場合には、信頼度0として、選択しないようにしてもよい。   For example, the minimum hamming distance and the number of learning data are used as parameters. For example, when the learning data regarding the selected blood pressure is sufficiently larger than a predetermined value, the reliability is high. For example, the greater the number of learning data, the higher the reliability, and when the learning data is equal to or less than a predetermined determination value, the reliability may be set to 0 and not selected.

・次に、信頼度を考慮して判定された血圧を表示する。なお、この際には、信頼度に関するデータも表示してもよい。例えば信頼度を学習データの数に応じて、低、中、高に分けて表示することができる。なお、学習データが所定の判定値より少ない場合には、その旨(例えば「測定の信頼性が低いので、血圧は表示しません」等)を表示して、血圧を表示しないことが好ましい。   -Next, the blood pressure determined in consideration of the reliability is displayed. At this time, data relating to reliability may also be displayed. For example, the reliability can be displayed divided into low, medium, and high according to the number of learning data. When the learning data is less than the predetermined determination value, it is preferable to display that fact (for example, “the blood pressure is not displayed because measurement reliability is low”) and the blood pressure is not displayed.

c)次に、上述した手順に従って血圧演算部37にて行われる血圧測定処理を説明する。   c) Next, blood pressure measurement processing performed in the blood pressure calculation unit 37 according to the above-described procedure will be described.

図8に示す様に、ステップ100(S)では、脈波センサ9及び電極11、13から、脈波信号及び心電信号を取得する。   As shown in FIG. 8, in step 100 (S), a pulse wave signal and an electrocardiogram signal are acquired from the pulse wave sensor 9 and the electrodes 11 and 13.

続くステップ110では、脈波信号及び心電信号を解析して、上述した血圧測定に必要な特徴量(ここでは、例えばPTT、HR、AI、加速度脈波a〜fなど)を算出する。   In the subsequent step 110, the pulse wave signal and the electrocardiogram signal are analyzed, and the characteristic quantities (for example, PTT, HR, AI, acceleration pulse waves a to f, etc.) necessary for the blood pressure measurement described above are calculated.

続くステップ120では、その様にして求めた特徴量が適正な値であるか否か(即ちあり得ないような異常値ではないか)を判定し、ここで肯定判断されるとステップ130に進み、一方、否定判断されるとステップ100に戻る。   In the following step 120, it is determined whether or not the feature quantity obtained in this way is an appropriate value (that is, an abnormal value that is not possible). If an affirmative determination is made here, the process proceeds to step 130. On the other hand, if a negative determination is made, the process returns to step 100.

ステップ130では、複数回の2値判別を行う。即ち、1mmHg毎の各血圧毎に設定された識別関数に前記特徴量を用い、その判定対象の血圧値未満か以上か(例えば80mmHg未満か以上か)を判定し、判定対象の血圧の数に対応した桁数を有する2値数列を作成する。   In step 130, binary discrimination is performed a plurality of times. That is, the feature amount is used for the discrimination function set for each blood pressure for each 1 mmHg, and it is determined whether the blood pressure value is less than or equal to the determination target (for example, less than 80 mmHg). Create a binary sequence with the corresponding number of digits.

続くステップ140では、前記ステップ130にて作成した2値数列と前記図6(b)に示す符号表とを比較し、ハミング距離を求め、その最小ハミング距離に対応した血圧値を選択する。   In the following step 140, the binary sequence created in step 130 is compared with the code table shown in FIG. 6B to obtain a Hamming distance, and a blood pressure value corresponding to the minimum Hamming distance is selected.

続くステップ150では、例えば学習データの数により、信頼度の判定を行う。例えば血圧値が100mmHgと判断された場合には、99mmHg、100mmHg、101mmHgにおける学習データ数が各10以上あるか否かを判定する。ここで、学習データ数が所定値以上であり、信頼度が高いと判定された場合には、ステップ170に進み、一方、信頼度が低いと判定された場合には、ステップ160に進む。   In the subsequent step 150, the reliability is determined based on the number of learning data, for example. For example, when it is determined that the blood pressure value is 100 mmHg, it is determined whether there are 10 or more learning data numbers at 99 mmHg, 100 mmHg, and 101 mmHg. If it is determined that the number of learning data is equal to or greater than the predetermined value and the reliability is high, the process proceeds to step 170. On the other hand, if the reliability is determined to be low, the process proceeds to step 160.

ステップ170では、高い信頼度の血圧が得られたので、この血圧を測定値として決定する。   In step 170, since a highly reliable blood pressure is obtained, this blood pressure is determined as a measurement value.

続くステップ180では、前記ステップ170で決定された血圧を、報知装置3のディスプレイに表示し、一旦本処理を終了する。   In the subsequent step 180, the blood pressure determined in step 170 is displayed on the display of the notification device 3, and the process is temporarily terminated.

一方、前記ステップ170で、信頼度が低いと判断されて進む、ステップ160では、特徴量に基づいて予め生体情報の特徴量と血圧との関係により得られた重回帰式によって血圧を算出し、続くステップ170で、その血圧を今回の測定として決定し、ステップ180にて、その血圧をディスプレイに表示し、一旦本処理を終了する。   On the other hand, in step 160, it is determined that the reliability is low, and in step 160, the blood pressure is calculated by a multiple regression equation obtained in advance from the relationship between the feature amount of the biological information and the blood pressure based on the feature amount, In step 170, the blood pressure is determined as the current measurement. In step 180, the blood pressure is displayed on the display, and the process is temporarily terminated.

なお、重回帰により血圧を算出する手法については、下記実施例2において説明するが、信頼度が低い場合には、重回帰による血圧推定を行うことなく、そのまま、血圧測定できない旨を表示して、一旦本処理を終了してもよい。   The technique for calculating blood pressure by multiple regression will be described in Example 2 below, but if the reliability is low, it is displayed that blood pressure cannot be measured without performing blood pressure estimation by multiple regression. This processing may be temporarily terminated.

d)この様に、本実施例では、アダブーストを利用した2値化判定により、各血圧値に対応した複数回の2値化判定を実施し、それによって得られた2値数列を符号表の2値数列と比較し、最小ハミング距離に対応した血圧を今回の測定値として決定するので、真値の曖昧な測定対象であっても、学習精度を向上でき、結果、精度良く血圧測定を行うことができる。   d) As described above, in this embodiment, binarization determination using Adaboost is performed a plurality of times of binarization determination corresponding to each blood pressure value, and the binary sequence obtained thereby is represented in the code table. Compared with the binary sequence, the blood pressure corresponding to the minimum Hamming distance is determined as the current measurement value, so that even if the measurement value is ambiguous, the learning accuracy can be improved. As a result, the blood pressure measurement is performed with high accuracy. be able to.

また、2値化判定によって得られた血圧値の信頼性の判定を行っているので、その点からも精度の高い測定値が得られるという利点がある。   In addition, since the reliability of the blood pressure value obtained by the binarization determination is determined, there is an advantage that a highly accurate measurement value can be obtained from this point.

次に、実施例2を説明するが、前記実施例1と同様な内容の説明は省略する。   Next, the second embodiment will be described, but the description of the same contents as the first embodiment will be omitted.

なお、本実施例は、ハード構成は前記実施例1と同様であるので、血圧測定の手順などについて説明する。   Since the hardware configuration of this embodiment is the same as that of the first embodiment, the procedure for measuring blood pressure will be described.

本実施例は、学習によって血圧を推定する処理と、重回帰によって血圧を推定する処理の両方を行って、より高い精度で血圧を推定するものである。   In this embodiment, the blood pressure is estimated with higher accuracy by performing both the process of estimating the blood pressure by learning and the process of estimating the blood pressure by multiple regression.

a)ここでは、学習によって血圧を推定する手法は、前記実施例1と同様であるので、血圧の算出に用いる重回帰の演算式について説明する。   a) Here, since the technique for estimating blood pressure by learning is the same as that in the first embodiment, an arithmetic expression for multiple regression used for calculating blood pressure will be described.

この重回帰による血圧の推定については、例えば既に本発明者により提案された手法(特願2007−154086号参照)を採用できる。以下簡単に説明する。   For estimation of blood pressure by this multiple regression, for example, a method already proposed by the present inventor (see Japanese Patent Application No. 2007-154086) can be employed. This will be briefly described below.

ここでは、中枢血管特徴量(PTT、脈波面積)と末梢血管特徴量(容積AI)とを用いた演算式、例えば下記式(7)〜(11)のいずれかの演算式を用いて、血圧(推定血圧:EBP)を求める。   Here, using an arithmetic expression using the central vascular feature (PTT, pulse wave area) and the peripheral vascular feature (volume AI), for example, any one of the following expressions (7) to (11), The blood pressure (estimated blood pressure: EBP) is obtained.

・PTTを使用する場合
EBP=α・PTT+β・容積AI+γ ・・(7)
EBP=α・PTT+β・容積AI+γ・(脈波特徴量)+・・+ ・・(8)
EBP=α・PTT+β・容積AI+γ・d+δ・e+ε・f
+ζ・HR+・・+ ・・(9)
・PTTを使用しない場合
EBP=α・脈波面積+β・容積AI+γ ・・(10)
・より精度の良い推定式
EBP=α・PTT+β・容積AI+γ・(脈波特徴量)+δ・(身体特徴量)
+ε・(インピーダンス)+ζ ・・(11)
前記各式で、α、β、γ、δ、ε、ζは係数である。脈波特徴量は、例えばa1、b1など(図4参照)である。HRは脈拍であり、脈拍に代えて脈拍間隔(RRI)を用いることができる。脈波面積とは、脈波の基線をゼロとした場合の波形に囲まれる部分の面積である。身体特徴量は、身長、体重、性別、年齢などである。
・ When using PTT EBP = α ・ PTT + β ・ Volume AI + γ ・ ・ (7)
EBP = α · PTT + β · Volume AI + γ · (Pulse wave feature) + · · + · · · (8)
EBP = α · PTT + β · Volume AI + γ · d + δ · e + ε · f
+ Ζ · HR + · · · · · · (9)
・ When PTT is not used EBP = α ・ Pulse wave area + β ・ Volume AI + γ (10)
・ Estimated expression with higher accuracy EBP = α ・ PTT + β ・ Volume AI + γ ・ (Pulse wave feature) + δ ・ (Body feature)
+ Ε · (Impedance) + ζ (11)
In the above equations, α, β, γ, δ, ε, and ζ are coefficients. The pulse wave feature amounts are, for example, a1 and b1 (see FIG. 4). HR is a pulse, and a pulse interval (RRI) can be used instead of a pulse. The pulse wave area is the area of the portion surrounded by the waveform when the baseline of the pulse wave is zero. The body feature amount includes height, weight, sex, age, and the like.

以下、各式について説明する。   Hereinafter, each formula will be described.

前記式(7)では、推定血圧の算出に、PTTと容積AIを用いる。このうち、PTTは中枢血管の血管弾性を示すパラメータであり、容積AIは末梢血管の血管弾性を示すパラメータであるので、両パラメータを用いることにより、中枢、末梢系の情報を取り入れられているため、精度良く血圧を推定することができる。   In the formula (7), PTT and volume AI are used for calculating the estimated blood pressure. Among these, PTT is a parameter indicating the vascular elasticity of the central blood vessel, and the volume AI is a parameter indicating the vascular elasticity of the peripheral blood vessel. Therefore, by using both parameters, information on the central and peripheral systems is taken in. The blood pressure can be estimated with high accuracy.

前記式(8)では、推定血圧の算出に、PTTと容積AIに加えて、脈波特徴量を用いる。よって、前記式(7)よりも精度良く血圧を推定することができる。ここで、脈波特徴量を加味すると推定精度が向上するのは、特徴量の非線形性の効果、及び心臓の拍出の影響等を取り入れられているからである。   In the equation (8), the pulse wave feature amount is used in addition to the PTT and the volume AI for calculation of the estimated blood pressure. Therefore, the blood pressure can be estimated with higher accuracy than the equation (7). The reason why the estimation accuracy is improved when the pulse wave feature amount is taken into account is that the effect of nonlinearity of the feature amount, the influence of the cardiac output, and the like are incorporated.

前記式(9)では、推定血圧の算出に、PTTと容積AIに加えて、加速度脈波のd、e、f、及びHRを用いる。よって、前記式(7)よりも精度良く血圧を推定することができる。ここで、加速度脈波のd、e、fとHRを加味すると推定精度が向上するのは、後述の容積AIを圧AIに補正する効果、及び心臓の拍出の影響の取り入れがあるからである。   In the formula (9), the estimated blood pressure is calculated using d, e, f, and HR of the acceleration pulse wave in addition to the PTT and the volume AI. Therefore, the blood pressure can be estimated with higher accuracy than the equation (7). Here, when the acceleration pulse wave d, e, f and HR are taken into account, the estimation accuracy is improved because of the effect of correcting the volume AI described later to the pressure AI and the influence of the cardiac output. is there.

前記式(10)では、推定血圧の算出に、脈波面積(S)と容積AIを用いている。このうち、Sは中枢血管の血管弾性を示す特徴量であると考えられ(実験により推定)、容積AIは末梢血管の血管弾性を示す特徴量であるので、両特徴量を用いることにより、マニュアルにて体重を入力しなくても、精度良く血圧を推定することができる。   In the equation (10), the pulse wave area (S) and the volume AI are used to calculate the estimated blood pressure. Of these, S is considered to be a feature value indicating the vascular elasticity of the central blood vessel (estimated by experiment), and the volume AI is a feature value indicating the vascular elasticity of the peripheral blood vessel. The blood pressure can be accurately estimated without inputting the body weight.

前記式(11)では、推定血圧の算出に、PTTと容積AIに加えて、脈波特徴量、身体特徴量、インピーダンス(例えば両手間で測定したインピーダンス)を用いている。よって、前記式(8)よりも精度良く血圧を推定することができる。なお、身長や体重などの身体特徴量は、血流量と相関関係があり、また、インピーダンスは体重や体脂肪と相関関係がある(よって血流量とも相関関係がある)ので、これらを加味することにより、血圧の推定精度が向上すると考えられる。   In the equation (11), in addition to the PTT and the volume AI, a pulse wave feature value, a body feature value, and impedance (for example, impedance measured between both hands) are used for calculating the estimated blood pressure. Therefore, the blood pressure can be estimated with higher accuracy than the equation (8). Body features such as height and weight have a correlation with blood flow, and impedance has a correlation with body weight and body fat (and thus also with blood flow). Thus, it is considered that blood pressure estimation accuracy is improved.

b)次に、本実施例における血圧測定の手順について説明する。   b) Next, the procedure of blood pressure measurement in the present embodiment will be described.

本実施例では、前記実施例1のアダブースト及びECOCの手法を利用して、血圧の推定を行うとともに、重回帰により血圧の推定を行うものである。   In this embodiment, blood pressure is estimated using the Adaboost and ECOC methods of the first embodiment, and blood pressure is estimated by multiple regression.

・図9に示す様に、脈波センサ9及び電極11、13からの信号から、心電信号及び脈波信号の特徴量を求める。   As shown in FIG. 9, the feature quantity of the electrocardiogram signal and the pulse wave signal is obtained from the signals from the pulse wave sensor 9 and the electrodes 11 and 13.

具体的には、血圧の算出に用いるために、上述した特徴量(PTT、HR、速度脈波の特徴量a1〜k1、加速度脈波の特徴量a〜f、AI(容積AI)等)を算出する。   Specifically, the above-described feature amounts (PTT, HR, velocity pulse wave feature amounts a1 to k1, acceleration pulse wave feature amounts a to f, AI (volume AI), etc.) are used to calculate blood pressure. calculate.

・次に、アダブーストの検知モードにおける2値化判定を複数回実施して血圧を算出する。   Next, blood pressure is calculated by performing binarization determination in the AdaBoost detection mode a plurality of times.

即ち、前記実施例1と同様にして、事前学習モードで生成した識別関数に、今回算出された特徴量(例えばPTTやHRなど)を適用し、各血圧における2値化判定を行い、その各識別関数の判定結果の2値数列を作成する。そして、この2値数列と符号表とを対比して、最小ハミング距離を求め、最小ハミング距離となる血圧を選択する。   That is, as in the first embodiment, the feature amount calculated this time (for example, PTT or HR) is applied to the discriminant function generated in the pre-learning mode, and binarization determination is performed for each blood pressure. A binary number sequence of discrimination function determination results is created. Then, the binary number sequence and the code table are compared to determine the minimum hamming distance, and the blood pressure that is the minimum hamming distance is selected.

・また、前記重回帰の算出式を利用して、特徴量(例えばPTTやHRなど)から、血圧を算出する。   Further, the blood pressure is calculated from the characteristic amount (for example, PTT, HR, etc.) using the calculation formula for the multiple regression.

・次に、学習によって求めた血圧と重回帰によって求めた血圧とを比較して、血圧を決定する。   Next, the blood pressure obtained by learning is compared with the blood pressure obtained by multiple regression to determine the blood pressure.

例えば重回帰によって求めた血圧が115mmHg、学習によって求めた血圧が112mmHgとした場合を考えると、重回帰による血圧が、100mmHgと140mmHgの範囲であれば、学習による血圧を信頼度を高いデータとして採用する。これは、100mmHg〜140mmHgの範囲の学習のデータ数が多いからである。   For example, when the blood pressure obtained by multiple regression is 115 mmHg and the blood pressure obtained by learning is 112 mmHg, if the blood pressure obtained by multiple regression is in the range of 100 mmHg and 140 mmHg, the blood pressure obtained by learning is used as data with high reliability. To do. This is because the number of learning data in the range of 100 mmHg to 140 mmHg is large.

また、実際に学習による血圧のデータ数が10個の場合は、学習による血圧を採用し、3個以下ならば重回帰によるデータを採用する。これは、学習データが少ないほど、算出された血圧の信頼性が低いと考えられるからである。   Further, when the number of blood pressure data actually obtained by learning is 10, the blood pressure obtained by learning is adopted, and when the number is 3 or less, data obtained by multiple regression is adopted. This is because it is considered that the smaller the learning data, the lower the reliability of the calculated blood pressure.

c)次に、上述した手順に従って行われる血圧測定処理を説明する。   c) Next, blood pressure measurement processing performed according to the above-described procedure will be described.

図10に示す様に、ステップ200(S)では、脈波センサ9及び電極11、13から、脈波信号及び心電信号を取得する。   As shown in FIG. 10, in step 200 (S), a pulse wave signal and an electrocardiogram signal are acquired from the pulse wave sensor 9 and the electrodes 11 and 13.

続くステップ210では、脈波信号及び心電信号を解析して、上述した血圧測定に必要な特徴量(ここでは、例えばPTT、HR、AI、加速度脈波a〜fなど)を算出する。   In the subsequent step 210, the pulse wave signal and the electrocardiogram signal are analyzed, and the characteristic quantities (for example, PTT, HR, AI, acceleration pulse waves a to f, etc.) necessary for the blood pressure measurement described above are calculated.

続くステップ220では、その様にして求めた特徴量が適正な値であるか否か(即ちあり得ないような異常値ではないか)を判定し、ここで肯定判断されるとステップ230に進み、一方、否定判断されるとステップ200に戻る。   In the following step 220, it is determined whether or not the feature quantity obtained in this way is an appropriate value (that is, an abnormal value that is impossible). If an affirmative determination is made here, the process proceeds to step 230. On the other hand, if a negative determination is made, the process returns to step 200.

ステップ230では、重回帰による血圧の推定を行う。   In step 230, blood pressure is estimated by multiple regression.

続くステップ240では、前記アダブーストを利用した複数回の2値判別を行い、学習による血圧の推定を行う。   In the following step 240, binary discrimination is performed a plurality of times using the AdaBoost, and blood pressure is estimated by learning.

続きステップ250では、重回帰による血圧の推定値と学習による血圧の推定値を比較し、例えば上述した学習データの数に基づいて信頼性を判定し、信頼性の高い血圧値を採用する。   In subsequent step 250, the estimated value of blood pressure by multiple regression is compared with the estimated value of blood pressure by learning, for example, reliability is determined based on the number of learning data described above, and a highly reliable blood pressure value is adopted.

続くステップ260では、前記ステップ250で決定された血圧を、報知装置3のディスプレイに表示し、一旦本処理を終了する。   In the subsequent step 260, the blood pressure determined in the step 250 is displayed on the display of the notification device 3, and the process is temporarily terminated.

d)この様に、本実施例では、重回帰による血圧の推定と学習による血圧の推定を行って、信頼性の高い血圧を採用するので、常に精度の高い血圧測定を行うことができるという効果がある。   d) As described above, in this embodiment, blood pressure is estimated by multiple regression and blood pressure is estimated by learning, and a highly reliable blood pressure is adopted. Therefore, it is possible to always perform highly accurate blood pressure measurement. There is.

なお、本実施例では、前記実施例1の様に、識別関数をアダブーストにより生成したが、それに代えて、上述したニューラルネットワーク、SVM等の他のマシンラーンイングの手法を用いることができる。   In this embodiment, the discriminant function is generated by Adaboost as in the first embodiment, but other machine learning methods such as the above-described neural network and SVM can be used instead.

次に、実施例3を説明するが、前記実施例1と同様な内容の説明は省略する。   Next, the third embodiment will be described, but the description of the same contents as the first embodiment will be omitted.

本実施例では、幅を持つ符号表を使用する。つまり、例えば2飛ばしの符号表を用いる。   In this embodiment, a code table having a width is used. That is, for example, a code table with two skips is used.

前記実施例1では、図11(a)に示す様に、符号表の各2値数列は、血圧が1mmHgずつ減少するにつれて、右の桁から1ずつ増加するように設定されている。   In the first embodiment, as shown in FIG. 11A, each binary number sequence in the code table is set to increase by 1 from the right digit as the blood pressure decreases by 1 mmHg.

しかし、本実施例では、図11(b)に示す様に、符号表の各2値数列は、血圧が2mmHgずつ減少するにつれて、右側から2桁ずつ1が増加するように設定されている。よって、例えば100mmHgと101mmHgとが同じ数列で評価される。   However, in this embodiment, as shown in FIG. 11B, each binary number sequence in the code table is set so that 1 increases from the right by 2 digits as the blood pressure decreases by 2 mmHg. Therefore, for example, 100 mmHg and 101 mmHg are evaluated in the same number sequence.

例えば理想的な100mmHgの数列が入ったとき(0・・00111・・1)は、ハミング距離(符号表とのズレ)は、98で2、100で1、102で2となり、単純に100mmHgと判定が可能である。   For example, when an ideal sequence of 100 mmHg is entered (0 ··· 00111 ··· 1), the Hamming distance (deviation from the code table) is 2 at 98, 1 at 100, and 2 at 102, and is simply 100 mmHg. Judgment is possible.

また、例えば理想的な101mmHgの数列が入ったとき(0・・00011・・1)は、ハミング距離、98で3、100で1、102で1となり、100mmHgと102mmHgでハミング距離が最小となり、101mmHgと判定が可能である。   For example, when an ideal sequence of 101 mmHg is entered (0..00011..1), the Hamming distance is 3, 98 at 1, 100 at 1, 102 at 1, and the Hamming distance is minimum at 100 mmHg and 102 mmHg. It can be determined as 101 mmHg.

これにより、符号表への照会が、図11(a)の場合に比べて半分で済む(約半分の41回)ため、計算時間を短縮できるという利点がある。   Thereby, since the reference to the code table is half as compared with the case of FIG. 11A (about 41 times, which is about half), there is an advantage that the calculation time can be shortened.

次に、実施例4を説明するが、前記実施例1と同様な内容の説明は省略する。   Next, the fourth embodiment will be described, but the description of the same contents as the first embodiment will be omitted.

本実施例では、被験者に適合した符号表を用いる。   In this embodiment, a code table suitable for the subject is used.

例えば高齢者の場合、血圧が100mmHgである可能性は低いので、例えば図12に示す様に、枠内の2値数列の無い符号表を用いる。つまり、例えば年齢に応じて(2値化判別する血圧である可能性が少ない場合は、その血圧に対応した数列を除いた)数列の少ない符号表を用いることができる。   For example, in the case of an elderly person, there is a low possibility that the blood pressure is 100 mmHg. Therefore, for example, as shown in FIG. 12, a code table without a binary sequence within the frame is used. In other words, for example, a code table with a small number of sequences can be used according to age (excluding the number sequence corresponding to the blood pressure when there is little possibility of binarized blood pressure).

これにより、計算時間を短縮することができる。   Thereby, calculation time can be shortened.

次に、実施例5を説明するが、前記実施例1と同様な内容の説明は省略する。   Next, the fifth embodiment will be described, but the description of the same contents as the first embodiment will be omitted.

本実施例では、血圧推定を所定の血圧範囲の群として行う。   In this embodiment, blood pressure estimation is performed as a group of predetermined blood pressure ranges.

つまり、下記の様に、例えば4mmHg毎のステップ毎に血圧推定を行う。   That is, as described below, for example, blood pressure is estimated for each step of 4 mmHg.

・・・
100mmHg未満/104mmHg以上
104mmHg未満/108mmHg以上
・・・
これにより、計算時間を短縮することができる。
...
Less than 100mmHg / 104mmHg or more Less than 104mmHg / 108mmHg or more
...
Thereby, calculation time can be shortened.

次に、実施例6を説明するが、前記実施例1と同様な内容の説明は省略する。   Next, the sixth embodiment will be described, but the description of the same contents as the first embodiment will be omitted.

本実施例では、血圧推定を所定の血圧範囲の群として行うが、ある範囲では群の幅を狭く、ある範囲では群の幅を広くする。   In this embodiment, blood pressure estimation is performed as a group of a predetermined blood pressure range, but the width of the group is narrow in a certain range, and the width of the group is widened in a certain range.

つまり、例えば100〜140mmHgの場合には、群の幅を1mmHgと狭く、100mmHg未満や140〜180mmHgでは、群の幅を4mmHgと中程度に、180mmHg以上では、群の幅を8mmHg以上と広くする。   That is, for example, in the case of 100 to 140 mmHg, the width of the group is as narrow as 1 mmHg, and in the case of less than 100 mmHg or in the range of 140 to 180 mmHg, the width of the group is as moderate as 4 mmHg. .

このように設定した理由は、学習データが少ない場合にはデータの信頼性が低くなるということと、血圧のデータは平均値の周辺の学習データは多く、それを外れるに従って学習データが少なくなるという現象があるからである。   The reason for setting in this way is that the reliability of the data is low when the learning data is small, and that the blood pressure data has a lot of learning data around the average value, and the learning data decreases as it goes away. This is because there is a phenomenon.

これにより、計算時間を短縮することができるとともに、学習データの数による曖昧性を低減することができるという利点がある。   Thereby, there is an advantage that calculation time can be shortened and ambiguity due to the number of learning data can be reduced.

次に、実施例7を説明するが、前記実施例1と同様な内容の説明は省略する。   Next, although Example 7 is demonstrated, description of the content similar to the said Example 1 is abbreviate | omitted.

一般に学習の真値とするカフ血圧計には、通常4〜8mmHgの誤差がある。従って、本実施例では、その点を考慮して、2値化の学習データの血圧値から、所定幅以内(例えば4mmHg幅以内、8mmHg幅以内など)のデータを用いないようにした。   Generally, a cuff sphygmomanometer that is a true value of learning usually has an error of 4 to 8 mmHg. Therefore, in this embodiment, in consideration of this point, data within a predetermined range (for example, within 4 mmHg width, within 8 mmHg width, etc.) is not used from the blood pressure value of the binarized learning data.

これにより、正しい真値のデータを学習に使用でき、血圧の推定精度が向上する。   Thereby, correct true value data can be used for learning, and blood pressure estimation accuracy is improved.

次に、実施例8を説明するが、前記実施例1と同様な内容の説明は省略する。   Next, although Example 8 will be described, description of the same contents as in Example 1 will be omitted.

本実施例では、2値化のための学習データ数をある範囲内に納めるものである。   In this embodiment, the number of learning data for binarization is within a certain range.

各血圧層における学習データ数が大きく異なる場合には、血圧の推定精度が低下することがある。例えば100mmHg未満/以上の学習データ数が1:1の様に、学習データ数の割合が同じ程度の場合には、血圧値の推定精度が高いが、例えば140mmHg未満/以上のデータ数が8:1の様に、学習データ数の割合が大きく異なる場合には、血圧値の推定精度が低い傾向が見られる。   When the number of learning data in each blood pressure layer is greatly different, blood pressure estimation accuracy may decrease. For example, when the ratio of the number of learning data is the same, such as the number of learning data of less than 100 mmHg is 1: 1, the accuracy of blood pressure estimation is high, but the number of data of less than 140 mmHg is 8: As shown in FIG. 1, when the ratio of the number of learning data is greatly different, the blood pressure value estimation accuracy tends to be low.

この理由は、その脈波が本来属する血圧に同じような特徴量を持つ学習データがあったとしても、その脈波が本来属しない血圧にも同じような特徴量を持つ学習データが存在する確率が高くなり、結果として、血圧の推定精度が低下するからである。   The reason for this is that even if there is learning data with the same feature amount in the blood pressure to which the pulse wave originally belongs, the probability that there is learning data with the same feature amount in the blood pressure to which the pulse wave does not originally belong This is because the accuracy of blood pressure estimation decreases as a result.

従って、本実施例では、2値化判定のために所定の血圧値で区分した各血圧層(例えば140mmHg未満/以上と2つに区分された各血圧層)における学習データのデータ数の比が、所定範囲内(例えば2倍以内)に設定している。   Therefore, in this embodiment, the ratio of the number of learning data in each blood pressure layer (for example, each blood pressure layer divided into two, less than 140 mmHg) divided by a predetermined blood pressure value for binarization determination is , Within a predetermined range (for example, within 2 times).

これにより、データ数に引きずられる効果は少なくなるので、血圧の測定精度が向上する。   Thereby, since the effect dragged by the number of data is reduced, blood pressure measurement accuracy is improved.

つまり、学習データを利用する場合には、予め所定の血圧範囲のデータをどの割合で使用するかを設定することが望ましい。なお、好ましいデータ数の割合として、例えば5倍以内に納まる程度が挙げられる。
<<実験例>>
次に、本発明の効果を確認するために行った実験例について説明する。
That is, when learning data is used, it is desirable to set in advance at what ratio data in a predetermined blood pressure range is used. In addition, as a preferable ratio of the number of data, for example, a degree within 5 times can be mentioned.
<< Experimental example >>
Next, experimental examples conducted for confirming the effects of the present invention will be described.

本実験は、インターネットにより公開されているMatlabのAdaboost Toolbox(http://research.graphicon.ru/machine-learnining/gml-adaboost-matlab-toolbox.html)を使用して行ったものである。   This experiment was conducted using Matlab's Adaboost Toolbox (http://research.graphicon.ru/machine-learnining/gml-adaboost-matlab-toolbox.html) published on the Internet.

以下の実験では、学習に用いる特徴量としては、PTT、HR、AI、加速度脈波a〜f、及び、身体特徴量(身長、体重、性別、年齢)、インピーダンスを採用した。また、重回帰に用いる特徴量としては、PTT、AI、脈波特徴量(加速度脈波a〜f)、身体特徴量(身長、体重、性別、年齢)、インピーダンスを採用し、重回帰に用いる式としては、前記式(11)を採用した。   In the following experiments, PTT, HR, AI, acceleration pulse waves a to f, body feature amounts (height, weight, sex, age), and impedance were adopted as feature amounts used for learning. Further, as feature quantities used for multiple regression, PTT, AI, pulse wave feature quantities (acceleration pulse waves a to f), body feature quantities (height, weight, sex, age), and impedance are adopted and used for multiple regression. As the formula, the formula (11) was adopted.

また、弱学習機は、CART(Classification And Regression Trees(http://research.graphicon.ru/machine-learning/index.php)を用い、個数は20個とした。   Moreover, the weak learning machine used CART (Classification And Regression Trees (http://research.graphicon.ru/machine-learning/index.php), and the number was 20 pieces.

<実験例1>アダブーストと重回帰の比較
本実験例は、全データに対して、単純にアダブースト、重回帰によって血圧を推定したものである。なお、アダブーストは、1mmHgステップで実施した。実験結果を図13に示す。
<Experimental Example 1> Comparison between AdaBoost and Multiple Regression In this experimental example, blood pressure is simply estimated by Adaboost and multiple regression for all data. Note that AdaBoost was implemented in 1 mmHg steps. The experimental results are shown in FIG.

図13(a)は、事前学習のための学習データ(個数N=345)における推定血圧とカフ血圧のグラフを示し、図13(b)は、検知モードの評価データ(個数N=86)における推定血圧とカフ血圧のグラフを示したものである。   FIG. 13A shows a graph of estimated blood pressure and cuff blood pressure in learning data (number N = 345) for pre-learning, and FIG. 13B shows detection data in evaluation mode (number N = 86). A graph of estimated blood pressure and cuff blood pressure is shown.

図13(b)から明かな様に、血圧の中心部(標準的な範囲:100〜140mmHg)では、アダブーストは重回帰と同程度もしくは重回帰より推定精度が良いことが分かる。   As is clear from FIG. 13B, it can be seen that AdaBoost has the same degree of accuracy as the multiple regression or better estimation accuracy than the multiple regression in the central part of blood pressure (standard range: 100 to 140 mmHg).

従って、血圧の中心部では、重回帰よりもアダブーストによって得られた血圧値を採用することが望ましい。   Therefore, it is desirable to adopt the blood pressure value obtained by Adaboost rather than multiple regression at the center of blood pressure.

<実験例2>各血圧層の推定精度
本実験例は、前記実験例1と同様に、各血圧層における重回帰とアダブーストの血圧推定精度を調べたものである。実験結果を図14に示す。
<Experimental Example 2> Estimated Accuracy of Each Blood Pressure Layer In this experimental example, similar to Experimental Example 1, the blood pressure estimation accuracy of multiple regression and Adaboost in each blood pressure layer was examined. The experimental results are shown in FIG.

図14(a)は、学習データにおける各血圧範囲と標準偏差のグラフを示し、図14(b)は、評価データにおける各血圧範囲と標準偏差のグラフを示したものである。   FIG. 14 (a) shows a graph of each blood pressure range and standard deviation in the learning data, and FIG. 14 (b) shows a graph of each blood pressure range and standard deviation in the evaluation data.

図14から、110〜130mmHgの血圧範囲、学習データ数Nが50個以上(1mmHg当たり5個以上)では、重回帰よりアダブーストの方が精度が良いことが分かる。   From FIG. 14, it can be seen that in the blood pressure range of 110 to 130 mmHg and the learning data number N is 50 or more (5 or more per 1 mmHg), Adaboost has better accuracy than multiple regression.

従って、上記条件を満たす場合には、アダブーストによって推定した血圧を採用することが望ましい。   Therefore, when the above conditions are satisfied, it is desirable to adopt the blood pressure estimated by Adaboost.

<実験例3>群(ステップ)による血圧推定
本実験例は、アダブーストによる推定を行う場合に、学習データとして、1mmHg毎ではなく、ある幅(例えば4mmHg等のステップ幅)を持ったデータを用いて、群として血圧の推定を行うものである。実験結果を図15に示す。
<Experimental example 3> Blood pressure estimation by group (step) This experimental example uses data having a certain width (for example, a step width of 4 mmHg, etc.) instead of every 1 mmHg as the learning data when estimation is performed by Adaboost. Thus, blood pressure is estimated as a group. The experimental results are shown in FIG.

図15(a)は、推定血圧とカフ血圧との関係を示すグラフであり、図15(b)は、誤差の標準偏差(誤差の平均値を0として換算)とステップ毎の血圧との関係を示すグラフである。   FIG. 15A is a graph showing the relationship between the estimated blood pressure and the cuff blood pressure, and FIG. 15B is the relationship between the standard deviation of the error (converted with the average value of the error being 0) and the blood pressure at each step. It is a graph which shows.

図15から、血圧の中央部はステップ幅が小さく方が、周辺部はステップ幅が大きい方が推定精度が良いことが分かる。   From FIG. 15, it can be seen that the estimation accuracy is better when the step width is smaller in the central portion of the blood pressure and the step width is larger in the peripheral portion.

従って、血圧の中央部、例えば100〜140mmHgの範囲では血圧推定のステップ幅を小さく、それより外側の領域ではステップ幅を大きくすることが望ましい。   Therefore, it is desirable to decrease the blood pressure estimation step width in the central portion of the blood pressure, for example, in the range of 100 to 140 mmHg, and increase the step width in the outer region.

<実験例4>ハミング距離と誤差
本実験例は、ハミング距離と推定血圧の誤差(カフ血圧と推定血圧との差)との関係を調べたものである。実験結果を図16に示す。
<Experimental example 4> Hamming distance and error In this experimental example, the relationship between the Hamming distance and the error of the estimated blood pressure (the difference between the cuff blood pressure and the estimated blood pressure) was examined. The experimental results are shown in FIG.

図16(a)は、最小ハミング距離と誤差(mmHg)との関係を調べたグラフを示し、図16(b)は、最小ハミング距離が10未満と10以上の割合と誤差との関係を示すグラフである。   FIG. 16A shows a graph in which the relationship between the minimum hamming distance and the error (mmHg) is examined, and FIG. 16B shows the relationship between the error when the minimum hamming distance is less than 10 and 10 or more. It is a graph.

図16から、最小ハミング距離が大きい方が誤差が大きくなる傾向がある。   From FIG. 16, the error tends to increase as the minimum Hamming distance increases.

従って、最小ハミング距離が小さいものを採用することが好ましい。   Accordingly, it is preferable to employ a small hamming distance.

<実験例5>学習データの選別の効果
この学習データの選別とは、間違っている可能性のあるデータを除くことにより、推定精度を高めるものである。
<Experimental example 5> Effect of selection of learning data This selection of learning data is to improve estimation accuracy by excluding data that may be wrong.

本実験例は、カフ血圧計には、通常4〜8mmHgの誤差があることを考慮して、2値化の学習データの血圧値から、所定幅以内のデータを用いないようにしたものである。   In this experimental example, in consideration of the fact that the cuff sphygmomanometer normally has an error of 4 to 8 mmHg, data within a predetermined range is not used from the blood pressure value of binarized learning data. .

つまり、学習の真値となるリファレンス血圧自体に誤差があるので、誤ったデータを学習に用いると、誤った結果になることがあるため、例えば100mmHg未満/以上で2値化を行う場合に、誤っている可能性の高い96〜104mmHg(±4mmHg)のデータを、学習データから外すことで、誤っている可能性の高いデータを除去し、正しい可能性の高いデータを用いて学習させるというものである。実験結果を図17に示す。   In other words, since there is an error in the reference blood pressure itself that is the true value of learning, using incorrect data for learning may result in an incorrect result. For example, when binarization is performed at less than or equal to 100 mmHg, By removing 96-104 mmHg (± 4 mmHg) data that is highly likely to be incorrect from the learning data, data that is likely to be incorrect is removed, and learning is performed using data that is likely to be correct. It is. The experimental results are shown in FIG.

図17(a)は、データの除去に対応した推定血圧とカフ血圧との関係を示すグラフであり、図17(b)は、データの除去に対応した誤差と血圧との関係を示すグラフであり、図17(c)は、各血圧層に分けたデータ数とデータ除去幅(データを除去する幅)に対応した推定血圧の誤差の標準偏差(誤差の平均値を0として換算)を示す表である。   FIG. 17A is a graph showing the relationship between estimated blood pressure and cuff blood pressure corresponding to data removal, and FIG. 17B is a graph showing the relationship between error and blood pressure corresponding to data removal. Yes, FIG. 17 (c) shows the standard deviation of error of estimated blood pressure corresponding to the number of data divided into each blood pressure layer and the data removal width (data removal width) (converted with the average value of the error as 0). It is a table.

図17から、4〜8mmHg幅より内側のデータを学習に用いずに、2値化学習させた方が、精度が1mmHgほど向上することが分かる。   FIG. 17 shows that the accuracy is improved by 1 mmHg when binarized learning is performed without using the data inside the 4-8 mmHg width for learning.

<実験例6>学習データ比の影響
本実験例は、学習データを所定の血圧で区分する各血圧層において、学習データ数の多寡が血圧推定に及ぼす影響を調べたものである。
<Experimental Example 6> Influence of Learning Data Ratio In this experimental example, the influence of the number of learning data on blood pressure estimation is examined in each blood pressure layer that classifies learning data by a predetermined blood pressure.

この実験では、各血圧層において2倍以上のデータ数の違いがある場合には、多い方のデータ数を半分にして実験を行った。実験結果を図18に示す。   In this experiment, when there was a difference in the number of data more than twice in each blood pressure layer, the experiment was performed with the larger number of data halved. The experimental results are shown in FIG.

詳しくは、例えば学習データの数が、99mmHg以下が39個、100〜109mmHgが60個、110〜119mmHgが69個、120〜129mmHgが79個、130〜139mmHgが55個、140〜149mmHgが24個、150mmHg以上が19個の場合、140mmHg未満は302個、140mmHg以上は43個となる。従って、この場合は、140mmHg未満の血圧層のデータを半分とする(例えば血圧の低いものから奇数又は偶数番を選択する)。   Specifically, for example, the number of learning data is 39 for 99 mmHg or less, 60 for 100 to 109 mmHg, 69 for 110 to 119 mmHg, 79 for 120 to 129 mmHg, 55 for 130 to 139 mmHg, 24 for 140 to 149 mmHg When there are 19 pieces of 150 mmHg or more, 302 pieces are below 140 mmHg and 43 pieces are 140 mmHg or more. Therefore, in this case, the data of the blood pressure layer of less than 140 mmHg is halved (for example, an odd or even number is selected from those having a low blood pressure).

図18(a)は、データ数を変化させた場合における推定血圧とカフ血圧との関係を示すグラフであり、図18(b)は、各血圧層に分けたデータ数を変化させた場合の推定血圧の誤差の標準偏差(誤差の平均値を0として換算)を示す表である。   FIG. 18A is a graph showing the relationship between the estimated blood pressure and the cuff blood pressure when the number of data is changed, and FIG. 18B is a graph when the number of data divided into each blood pressure layer is changed. It is a table | surface which shows the standard deviation (The average value of an error is converted into 0) of the error of estimated blood pressure.

図18から明かな様に、各血圧層において2倍以上のデータ数の差がある場合、データ数を調整しない場合(そのまま)に比べて、データ数の比(データ比)を補正した方(比の調整)が、精度が向上したことが確認できた。   As is clear from FIG. 18, when there is a difference in the number of data more than twice in each blood pressure layer, the ratio of the number of data (data ratio) is corrected compared to the case where the number of data is not adjusted (as it is) ( It was confirmed that the accuracy was improved.

従って、各血圧層のデータ数を近づけることが望ましいことが分かる。   Therefore, it can be seen that it is desirable to make the number of data of each blood pressure layer closer.

尚、本発明は前記実施例になんら限定されるものではなく、本発明を逸脱しない範囲において種々の態様で実施しうることはいうまでもない。   Needless to say, the present invention is not limited to the above-described embodiments, and can be implemented in various modes without departing from the scope of the present invention.

(1)例えば図7(b)に示す様に、符号表として、重み付けした数列を用いることができる。この場合は、対応する桁数が異なる場合には、異なる符号表の数値を加算してよく。これによって、大きく外れているデータのハミング距離が増大し、外れたデータを除去し、精度を向上できるという利点がある。   (1) For example, as shown in FIG. 7B, a weighted number sequence can be used as the code table. In this case, if the corresponding number of digits is different, the numerical values of different code tables may be added. This has the advantage that the hamming distance of the greatly deviated data is increased, the deviated data is removed, and the accuracy can be improved.

(2)また、例えば、上述した血圧測定装置の機能は、コンピュータのプログラムにより実行される処理により実現することができ、このプログラムは、記録媒体に記録できる。   (2) Further, for example, the function of the blood pressure measurement device described above can be realized by processing executed by a computer program, and this program can be recorded on a recording medium.

つまり、上述した様なプログラムをコンピュータシステムにて実現する機能は、例えば、コンピュータシステム側で起動するプログラムとして備えることができる。このようなプログラムの場合、例えば、フレキシブルディスク、光磁気ディスク、CD−ROM、DVD、ハードディスク等のコンピュータ読み取り可能な記録媒体に記録し、必要に応じてコンピュータシステムにロードして起動することにより用いることができる。この他、ROMやバックアップRAM等をコンピュータ読み取り可能な記録媒体として前記プログラムを記録しておき、このROMあるいはバックアップRAM等をコンピュータシステムに組み込んで用いても良い。   That is, the function for realizing the above-described program in the computer system can be provided as a program that is activated on the computer system side, for example. In the case of such a program, for example, the program is recorded on a computer-readable recording medium such as a flexible disk, a magneto-optical disk, a CD-ROM, a DVD, and a hard disk, and is used by being loaded into a computer system and started as necessary. be able to. In addition, the program may be recorded as a computer-readable recording medium such as a ROM or a backup RAM, and the ROM or the backup RAM may be incorporated into a computer system.

1…血圧測定装置
3…報知装置
5…マニュアル入力部
7…ステアリング
9…脈波センサ
11、13…電極
DESCRIPTION OF SYMBOLS 1 ... Blood pressure measuring device 3 ... Notification apparatus 5 ... Manual input part 7 ... Steering 9 ... Pulse wave sensor 11, 13 ... Electrode

Claims (20)

生体から検出された心電信号及び脈波信号のうち、少なくとも前記脈波信号から得られる生体情報を用いて、血圧を推定する血圧測定装置であって、
事前学習によって得られた前記生体情報の特徴量と血圧との関係に基づいて、異なる血圧毎に前記生体情報の特徴量から血圧を2値化判別する識別手段と、
血圧の推定を行う場合に、前記識別手段を用い、測定によって得られた前記生体情報の特徴量に基づいて、所定の血圧に対して2値化判別を行う2値化判別手段と、
前記2値化判別手段によって、前記生体情報の特徴量に基づいて異なる血圧に対してそれぞれ2値化判定を行い、該複数回の2値化判定の結果に基づいて前記血圧を推定する血圧推定手段と、
を備えたことを特徴とする血圧測定装置。
A blood pressure measurement device that estimates blood pressure using at least biological information obtained from a pulse wave signal among an electrocardiogram signal and a pulse wave signal detected from a living body,
Identification means for binarizing and discriminating blood pressure from the feature amount of the biological information for each different blood pressure based on the relationship between the feature amount of the biological information and blood pressure obtained by prior learning;
A binarization determination unit that performs binarization determination on a predetermined blood pressure based on a feature amount of the biological information obtained by measurement using the identification unit when estimating blood pressure;
Blood pressure estimation in which the binarization determination unit performs binarization determination on different blood pressures based on the feature amount of the biological information, and estimates the blood pressure based on the results of the multiple binarization determinations Means,
A blood pressure measurement device comprising:
前記学習識別手段は、アダブーストにより、前記血圧の2値化のための識別関数を生成したことを特徴とする請求項1に記載の血圧測定装置。   The blood pressure measurement device according to claim 1, wherein the learning identification unit generates an identification function for binarizing the blood pressure by Adaboost. 前記複数回の2値化判定の結果から、各血圧毎の判定結果を並べた2値数列を作成し、この2値数列と前もって準備された符号表における各血圧毎の数列とを対比して、前記血圧の推定を行うことを特徴とする請求項1又は2に記載の血圧測定装置。   A binary number sequence in which the determination results for each blood pressure are arranged from the result of the binarization determination of the plurality of times is created, and this binary number sequence is compared with the number sequence for each blood pressure in the code table prepared in advance. The blood pressure measurement apparatus according to claim 1, wherein the blood pressure is estimated. 前記複数回の2値化判定の結果から得られた2値数列と前記符号表の2値数列とから、両数列における最小ハミング距離を求め、該最小ハミング距離に基づいて、前記血圧の推定を行うことを特徴とする請求項3に記載の血圧測定装置。   A minimum hamming distance in both sequences is obtained from the binary sequence obtained from the result of the binary determination and the binary sequence in the code table, and the blood pressure is estimated based on the minimum hamming distance. The blood pressure measurement device according to claim 3, wherein the blood pressure measurement device is performed. 前記符号表の数列が、複数個ずつの幅で数値が異なるように設定されていることを特徴とする請求項3又は4に記載の血圧測定装置。   5. The blood pressure measurement device according to claim 3, wherein the numerical sequence of the code table is set so that a numerical value is different in a plurality of widths. 前記符号表として、被験者の状態に適応して設定された数列を有する符号表を用いることを特徴とする請求項3又は4に記載の血圧測定装置。   The blood pressure measurement device according to claim 3 or 4, wherein a code table having a numerical sequence set in conformity with the condition of the subject is used as the code table. 前記血圧の2値化の判定を、所定の血圧範囲の群ごとに行うことを特徴とする請求項1〜6のいずれかに記載の血圧測定装置。   The blood pressure measurement device according to claim 1, wherein the determination of binarization of the blood pressure is performed for each group of a predetermined blood pressure range. 所定の血圧の範囲では前記群の幅を狭くし、前記所定の血圧の範囲外では前記群の幅を広くして、前記血圧の推定を行うことを特徴とする請求項7に記載の血圧測定装置。   8. The blood pressure measurement according to claim 7, wherein the blood pressure is estimated by reducing the width of the group within a predetermined blood pressure range and increasing the width of the group outside the predetermined blood pressure range. apparatus. 前記学習識別手段における事前学習の際に、所定の血圧の範囲のデータを学習に用いないことを特徴とする請求項1〜8のいずれかに記載の血圧測定装置。   The blood pressure measurement device according to any one of claims 1 to 8, wherein data in a predetermined blood pressure range is not used for learning in advance learning in the learning identification unit. 前記学習識別手段における事前学習の際に、異なる血圧の範囲における二値化のデータ数の比を、所定の範囲内に設定することを特徴とする請求項1〜9のいずれかに記載の血圧測定装置。   The blood pressure according to any one of claims 1 to 9, wherein a ratio of the number of binarized data in different blood pressure ranges is set within a predetermined range at the time of prior learning in the learning identification means. measuring device. 前記血圧を推定する際に、データの信頼度を算出し、血圧推定結果とともに、前記データの信頼度を表示することを特徴とする請求項1〜10のいずれかに記載の血圧測定装置。   The blood pressure measurement device according to any one of claims 1 to 10, wherein when estimating the blood pressure, the reliability of data is calculated, and the reliability of the data is displayed together with a blood pressure estimation result. 前記データの信頼度を算出し、この信頼度が所定の判定値より低い場合には、より信頼度が高いと予想される異なる手法を採用して血圧の推定を行うことを特徴とする請求項11に記載の血圧測定装置。   The blood pressure is estimated by calculating a reliability of the data, and when the reliability is lower than a predetermined determination value, the blood pressure is estimated using a different method that is expected to have a higher reliability. 11. The blood pressure measurement device according to 11. 前記データの信頼度は、学習データの数又は最小ハミング距離に基づいて判断することを特徴とする請求項11又は12に記載の血圧測定装置。   The blood pressure measurement device according to claim 11 or 12, wherein the reliability of the data is determined based on the number of learning data or a minimum Hamming distance. 前記異なる手法は、重回帰による推定の手法であることを特徴とする請求項12又は13に記載の血圧測定装置。   The blood pressure measurement device according to claim 12 or 13, wherein the different method is an estimation method by multiple regression. 生体から検出された心電信号及び脈波信号のうち、少なくとも前記脈波信号から得られる生体情報を用いて、血圧を推定する血圧測定装置において、
前記生体情報の特徴量に基づいて、重回帰によって血圧を推定する第1血圧推定手段と、
前記生体情報の特徴量と血圧との関係の学習によって得られた判別手段を用い、前記特徴量から血圧を推定する第2血圧推定手段と、
前記第1血圧推定手段によって推定された第1推定血圧と、前記第2血圧推定手段によって推定された第2推定血圧とを、所定の信頼度の高低を判定する判定条件に基づいて比較して、前記信頼度が高い推定血圧を決定する血圧決定手段と、
を備えたことを特徴とする血圧測定装置。
In a blood pressure measurement device that estimates blood pressure using at least biological information obtained from the pulse wave signal among an electrocardiogram signal and a pulse wave signal detected from a living body,
First blood pressure estimating means for estimating blood pressure by multiple regression based on the feature amount of the biological information;
Second blood pressure estimating means for estimating blood pressure from the feature quantity using a discrimination means obtained by learning the relationship between the feature quantity of the biological information and blood pressure;
The first estimated blood pressure estimated by the first blood pressure estimating means and the second estimated blood pressure estimated by the second blood pressure estimating means are compared based on a determination condition for determining whether the predetermined reliability is high or low. Blood pressure determining means for determining the estimated blood pressure with high reliability,
A blood pressure measurement device comprising:
前記第2血圧推定手段は、
学習データである前記生体情報の特徴量と血圧との関係を予め機械学習することにより得られた識別手段であって、所定の血圧毎に前記生体情報の特徴量から血圧を2値化判別する学習識別手段と、
前記血圧の推定を行う場合に、前記学習識別手段を用い、測定によって得られた前記生体情報の特徴量に基づいて、前記所定の血圧に対して2値化判別を行う2値化判別手段と、
前記2値化判別手段によって、前記生体情報の特徴量に基づいて異なる前記所定の血圧に対してそれぞれ2値化判定を行い、該複数回の2値化判定の結果に基づいて前記血圧を推定する血圧推定手段と、
を備えたことを特徴とする請求項15に記載の血圧測定装置。
The second blood pressure estimating means includes
An identification means obtained by previously learning the relationship between the feature value of the biological information, which is learning data, and the blood pressure, and binarizing and determining the blood pressure from the feature value of the biological information for each predetermined blood pressure Learning identification means;
Binarization discriminating means for performing binarization discriminating on the predetermined blood pressure based on the feature quantity of the biological information obtained by measurement using the learning discriminating means when estimating the blood pressure; ,
The binarization determination means performs binarization determination on each of the different predetermined blood pressures based on the feature amount of the biological information, and estimates the blood pressure based on the results of the multiple binarization determinations Blood pressure estimating means for
The blood pressure measurement device according to claim 15, comprising:
前記重回帰により求めた第1推定血圧が、100〜140mmHgの場合には、血圧の推定値として、前記学習による第2推定血圧を採用することを特徴とする請求項15又は16に記載の血圧測定装置。   17. The blood pressure according to claim 15, wherein when the first estimated blood pressure obtained by the multiple regression is 100 to 140 mmHg, the second estimated blood pressure obtained by the learning is adopted as an estimated value of blood pressure. measuring device. 学習データが、1mmHg当たり5個以上の場合には、血圧の推定値として、前記学習による第2推定血圧を採用することを特徴とする請求項15〜17のいずれかに記載の血圧測定装置。   The blood pressure measuring device according to any one of claims 15 to 17, wherein when the learning data is 5 or more per 1 mmHg, the second estimated blood pressure obtained by the learning is adopted as an estimated value of blood pressure. 前記請求項1〜18のいずれかに記載の血圧測定装置の機能を実現させるためのプログラム。   The program for implement | achieving the function of the blood pressure measurement apparatus in any one of the said Claims 1-18. 前記請求項19に記載のプログラムを記録したコンピュータ読取可能な記録媒体。   A computer-readable recording medium on which the program according to claim 19 is recorded.
JP2009036789A 2009-02-19 2009-02-19 Blood pressure measuring device, program, and recording medium Expired - Fee Related JP5218139B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009036789A JP5218139B2 (en) 2009-02-19 2009-02-19 Blood pressure measuring device, program, and recording medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009036789A JP5218139B2 (en) 2009-02-19 2009-02-19 Blood pressure measuring device, program, and recording medium

Publications (2)

Publication Number Publication Date
JP2010187993A true JP2010187993A (en) 2010-09-02
JP5218139B2 JP5218139B2 (en) 2013-06-26

Family

ID=42814701

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009036789A Expired - Fee Related JP5218139B2 (en) 2009-02-19 2009-02-19 Blood pressure measuring device, program, and recording medium

Country Status (1)

Country Link
JP (1) JP5218139B2 (en)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012211847A (en) * 2011-03-31 2012-11-01 Oki Electric Ind Co Ltd Microvibration feature amount calculating device, microvibration feature amount calculation method and program
WO2015098977A1 (en) * 2013-12-25 2015-07-02 旭化成株式会社 Cardiac pulse waveform measurement device, portable device, medical device system, and vital sign information communication system
JP2016016295A (en) * 2014-07-11 2016-02-01 株式会社デンソー Blood pressure estimation apparatus
JP2016057701A (en) * 2014-09-05 2016-04-21 オムロン株式会社 Identification apparatus and control method for the same
JP2017067489A (en) * 2015-09-28 2017-04-06 富士フイルムRiファーマ株式会社 Diagnosis assistance device, method, and computer program
JP2020092738A (en) * 2018-12-10 2020-06-18 株式会社アルム Blood pressure estimation device, blood pressure estimation system, and blood pressure estimation program
CN112597789A (en) * 2019-09-03 2021-04-02 株式会社东书产业 Blood pressure estimation system, blood pressure estimation method, learning device, learning method, and program
JP2021078538A (en) * 2019-11-14 2021-05-27 日本光電工業株式会社 Learned model generation method, learned model generation system, inference device, and computer program
JP2021516387A (en) * 2018-03-05 2021-07-01 オムロン株式会社 Methods and devices for controlling the operation of mechanical devices and determining the reliability of data
JP2022042920A (en) * 2020-09-03 2022-03-15 Ssst株式会社 Biological information calculation system
JP7083194B1 (en) 2021-03-02 2022-06-10 Ssst株式会社 Biological information calculation system
JP7083195B1 (en) 2021-03-02 2022-06-10 Ssst株式会社 Biological information calculation system
JP7083193B1 (en) 2021-03-02 2022-06-10 Ssst株式会社 Biological information calculation system
KR20220094032A (en) * 2020-12-28 2022-07-05 이광춘 Health examination method

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102486700B1 (en) 2015-08-11 2023-01-11 삼성전자주식회사 Apparatus and method for estimating blood pressure

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000237153A (en) * 1999-02-19 2000-09-05 Nippon Colin Co Ltd Noninvasive continuous blood pressure estimating device and noninvasive continuous blood pressure predicting device
JP2006263354A (en) * 2005-03-25 2006-10-05 Denso Corp Apparatus for acquiring individual information of human body, sphygmomanometer, and pulse wave analyzer

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000237153A (en) * 1999-02-19 2000-09-05 Nippon Colin Co Ltd Noninvasive continuous blood pressure estimating device and noninvasive continuous blood pressure predicting device
JP2006263354A (en) * 2005-03-25 2006-10-05 Denso Corp Apparatus for acquiring individual information of human body, sphygmomanometer, and pulse wave analyzer

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012211847A (en) * 2011-03-31 2012-11-01 Oki Electric Ind Co Ltd Microvibration feature amount calculating device, microvibration feature amount calculation method and program
WO2015098977A1 (en) * 2013-12-25 2015-07-02 旭化成株式会社 Cardiac pulse waveform measurement device, portable device, medical device system, and vital sign information communication system
JPWO2015098977A1 (en) * 2013-12-25 2017-03-23 旭化成株式会社 Pulse wave measuring device, portable device, medical device system, and biological information communication system
US10624586B2 (en) 2013-12-25 2020-04-21 Asahi Kasei Kabushiki Kaisha Pulse wave measuring device, mobile device, medical equipment system and biological information communication system
JP2016016295A (en) * 2014-07-11 2016-02-01 株式会社デンソー Blood pressure estimation apparatus
JP2016057701A (en) * 2014-09-05 2016-04-21 オムロン株式会社 Identification apparatus and control method for the same
JP2017067489A (en) * 2015-09-28 2017-04-06 富士フイルムRiファーマ株式会社 Diagnosis assistance device, method, and computer program
JP2021516387A (en) * 2018-03-05 2021-07-01 オムロン株式会社 Methods and devices for controlling the operation of mechanical devices and determining the reliability of data
JP7099531B2 (en) 2018-03-05 2022-07-12 オムロン株式会社 Methods and devices for controlling the operation of mechanical devices and determining the reliability of data
JP2020092738A (en) * 2018-12-10 2020-06-18 株式会社アルム Blood pressure estimation device, blood pressure estimation system, and blood pressure estimation program
CN112597789A (en) * 2019-09-03 2021-04-02 株式会社东书产业 Blood pressure estimation system, blood pressure estimation method, learning device, learning method, and program
JP2021078538A (en) * 2019-11-14 2021-05-27 日本光電工業株式会社 Learned model generation method, learned model generation system, inference device, and computer program
JP7381301B2 (en) 2019-11-14 2023-11-15 日本光電工業株式会社 Trained model generation method, trained model generation system, inference device, and computer program
JP2022042920A (en) * 2020-09-03 2022-03-15 Ssst株式会社 Biological information calculation system
JP7083185B2 (en) 2020-09-03 2022-06-10 Ssst株式会社 Biological information calculation system
KR20220094032A (en) * 2020-12-28 2022-07-05 이광춘 Health examination method
KR102546779B1 (en) * 2020-12-28 2023-06-22 이광춘 Health examination method
JP7083193B1 (en) 2021-03-02 2022-06-10 Ssst株式会社 Biological information calculation system
JP7083195B1 (en) 2021-03-02 2022-06-10 Ssst株式会社 Biological information calculation system
JP2022133918A (en) * 2021-03-02 2022-09-14 Ssst株式会社 Biological information arithmetic system
JP2022133920A (en) * 2021-03-02 2022-09-14 Ssst株式会社 Biological information calculation system
JP2022133919A (en) * 2021-03-02 2022-09-14 Ssst株式会社 Biological information arithmetic system
JP7083194B1 (en) 2021-03-02 2022-06-10 Ssst株式会社 Biological information calculation system

Also Published As

Publication number Publication date
JP5218139B2 (en) 2013-06-26

Similar Documents

Publication Publication Date Title
JP5218139B2 (en) Blood pressure measuring device, program, and recording medium
JP4971041B2 (en) Blood pressure measuring device, program, and recording medium
EP3366203B1 (en) Method and system for cuffless blood pressure estimation using photoplethysmogram features and pulse transit time
US8668649B2 (en) System for cardiac status determination
CN101061950B (en) Atrial fibrilation detection by SP02
Tjahjadi et al. Noninvasive classification of blood pressure based on photoplethysmography signals using bidirectional long short-term memory and time-frequency analysis
US6503206B1 (en) Apparatus having redundant sensors for continuous monitoring of vital signs and related methods
WO2014132713A1 (en) Pulse wave propagation time calculation device
CN102551699B (en) The system that non-invasive cardiac output is determined
Mohebbian et al. Blind, cuff-less, calibration-free and continuous blood pressure estimation using optimized inductive group method of data handling
JP2006263354A (en) Apparatus for acquiring individual information of human body, sphygmomanometer, and pulse wave analyzer
US8382679B2 (en) Autonomic nerve activity measuring apparatus and autonomic nerve activity measuring method
JP2009089829A (en) Biological state estimating device, program, and recording medium
US10357164B2 (en) Method and device for non-invasive blood pressure measurement
US20130218038A1 (en) System for Non-invasive Cardiac Output Determination
KR20070056925A (en) Blood pressure measuring method and apparatus
Lee et al. Deep learning ensemble with asymptotic techniques for oscillometric blood pressure estimation
KR102154652B1 (en) Method for estimating continuous blood pressure using recurrent neural network and apparatus thereof
Argha et al. Artificial intelligence based blood pressure estimation from auscultatory and oscillometric waveforms: a methodological review
JP5141716B2 (en) Pulse wave determination device and blood pressure estimation device
Bote et al. Evaluation of blood pressure estimation models based on pulse arrival time
Jain et al. A sparse regression based approach for cuff-less blood pressure measurement
US9320445B2 (en) System for cardiac condition detection responsive to blood pressure analysis
Yang et al. Analysis of the radial pulse wave and its clinical applications: a survey
CN114340484A (en) Unsupervised real-time classification of arterial blood pressure signals

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110413

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20121101

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20121120

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130115

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20130205

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130218

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20160315

Year of fee payment: 3

R151 Written notification of patent or utility model registration

Ref document number: 5218139

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20160315

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees