JP5633875B2 - Ventilation characteristics judgment device - Google Patents

Ventilation characteristics judgment device Download PDF

Info

Publication number
JP5633875B2
JP5633875B2 JP2010167481A JP2010167481A JP5633875B2 JP 5633875 B2 JP5633875 B2 JP 5633875B2 JP 2010167481 A JP2010167481 A JP 2010167481A JP 2010167481 A JP2010167481 A JP 2010167481A JP 5633875 B2 JP5633875 B2 JP 5633875B2
Authority
JP
Japan
Prior art keywords
value
bioelectrical impedance
breathing
cpu
centering
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.)
Active
Application number
JP2010167481A
Other languages
Japanese (ja)
Other versions
JP2012024413A (en
Inventor
増尾 善久
善久 増尾
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.)
Tanita Corp
Original Assignee
Tanita 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 Tanita Corp filed Critical Tanita Corp
Priority to JP2010167481A priority Critical patent/JP5633875B2/en
Priority to EP11173341A priority patent/EP2407100A1/en
Priority to US13/181,171 priority patent/US20120016254A1/en
Publication of JP2012024413A publication Critical patent/JP2012024413A/en
Application granted granted Critical
Publication of JP5633875B2 publication Critical patent/JP5633875B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Description

本発明は、被験者の換気機能の特性を判定可能な装置に関する。   The present invention relates to an apparatus capable of determining characteristics of a subject's ventilation function.

生体電気インピーダンスを測定し、測定結果に基づいて、生体の状態を推定する各種の装置が従来から知られている。そのような装置の一つとして、特許文献1には、体幹生体電気インピーダンスに基づいて、肺活量を推定する技術が開示されている。   Various devices that measure bioelectric impedance and estimate the state of a living body based on the measurement result have been known. As one of such devices, Patent Document 1 discloses a technique for estimating vital capacity based on trunk bioelectrical impedance.

特開2007−50127号公報(段落0020参照)Japanese Patent Laying-Open No. 2007-50127 (see paragraph 0020)

しかしながら、肺活量の値は、被験者の体格(身長、体重)や年齢などに応じて異なる値を示すので、肺活量の値が低いからといって、被験者の胸部のうち呼吸に寄与する部位(例えば胸部骨格筋、肺胞組織、胸部骨格筋および肺胞組織の両方を含む)の機能が正常であるとは限らない。すなわち、特許文献1に開示された技術では、被験者の肺活量を推定することはできても、被験者の換気機能の特性までを判定することはできない。
そこで、本発明は、被験者の換気機能の特性を判定可能な装置を提供することを課題とする。
However, the value of vital capacity shows different values depending on the physique (height, weight), age, etc. of the subject. Therefore, because the value of vital capacity is low, the part of the subject's chest that contributes to breathing (for example, the chest) The function of skeletal muscle, alveolar tissue, both thoracic skeletal muscle and alveolar tissue is not necessarily normal. That is, with the technique disclosed in Patent Document 1, even if the subject's vital capacity can be estimated, the characteristics of the subject's ventilation function cannot be determined.
Then, this invention makes it a subject to provide the apparatus which can determine the characteristic of a test subject's ventilation function.

上述した課題を解決するため、本発明に係る換気特性判定装置は、被験者の肺の上部を含み、腹部を含まない体幹上部の第1生体電気インピーダンスと、被験者の肺の中下部および腹部を含む体幹中部の第2生体電気インピーダンスとを測定する生体電気インピーダンス測定部(170、200)と、第1生体電気インピーダンスの測定値と、第2生体電気インピーダンスの測定値とに基づいて、被験者の換気機能の特性を判定する判定部(170)とを備えることを特徴とする。好適な態様として、判定部は、第1生体電気インピーダンスの測定値と、第2生体電気インピーダンスの測定値とに基づいて、被験者の胸部のうち呼吸に寄与する部位の機能が正常であるか否かを判定する。   In order to solve the above-described problem, a ventilation characteristic determination device according to the present invention includes a first bioelectric impedance of an upper trunk including an upper part of a subject's lung and not including an abdomen, and a lower part and an abdomen of the subject's lung. Based on the bioelectrical impedance measuring unit (170, 200) for measuring the second bioelectrical impedance of the middle trunk including the measured value of the first bioelectrical impedance and the measured value of the second bioelectrical impedance And a determination unit (170) for determining characteristics of the ventilation function. As a preferable aspect, the determination unit determines whether or not the function of the part contributing to breathing in the subject's chest is normal based on the measurement value of the first bioelectrical impedance and the measurement value of the second bioelectrical impedance. Determine whether.

例えば被験者が胸部に病歴を有しており、胸部のうち呼吸に寄与する部位(例えば胸部骨格筋、肺胞組織、胸部骨格筋および肺胞組織の両方を含む)の機能が低下していた場合は、呼吸に伴う胸部骨格筋の変化が健常者に比べて小さいので、換気量を確保するために、横隔膜の変動量(上昇量および下降量)が大きくなる。また、加齢などにより、胸部のうち呼吸に寄与する部位の機能が低下していた場合も同様である。ここで、肺の上部を含み、腹部を含まない体幹上部の第1生体電気インピーダンスの変化は、肺に出入りする空気量と連動している一方、肺の中下部および腹部を含む第2生体電気インピーダンスの変化は、横隔膜の動きと連動しており、横隔膜の変動量が大きいほど第2生体電気インピーダンスの変化も大きくなる。したがって、例えば胸部のうち呼吸に寄与する部位の機能が低下していた被験者と健常者との間で、1呼吸において肺に出入りする空気の量(換気量)が同じであり、当該1呼吸における第1生体電気インピーダンスの変化が同じであっても、胸部のうち呼吸に寄与する部位の機能が低下していた被験者は、健常者に比べて横隔膜の変動量が大きいので、当該1呼吸における第2生体電気インピーダンスの変化が大きくなる。上述の本発明の態様では、被験者の胸部のうち呼吸に寄与する部位の機能が低下している場合は、第2生体電気インピーダンスの変化が大きくなる点に着目し、第2生体電気インピーダンスの測定値と第1生体電気インピーダンスの測定値とに基づいて、被験者の胸部のうち呼吸に寄与する部位の機能が正常であるか否かを判定している。これにより、被験者の換気機能の特性(胸部のうち呼吸に寄与する部位の機能が低下している等)を判定できる。   For example, if the subject has a medical history in the chest and the function of the part of the chest that contributes to breathing (including chest skeletal muscle, alveolar tissue, both breast skeletal muscle and alveolar tissue) is reduced Since the change of the thoracic skeletal muscle accompanying breathing is smaller than that of a healthy person, the amount of fluctuation of the diaphragm (amount of rise and fall) is large in order to secure the ventilation. The same applies to the case where the function of the part of the chest that contributes to breathing is reduced due to aging or the like. Here, the change in the first bioelectric impedance of the upper trunk including the upper part of the lung and not including the abdomen is linked to the amount of air entering and exiting the lung, while the second living body including the middle lower part of the lung and the abdomen. The change in the electrical impedance is linked to the movement of the diaphragm, and the change in the second bioelectrical impedance increases as the amount of fluctuation of the diaphragm increases. Therefore, for example, the amount of air that enters and exits the lungs in one breath (ventilation amount) is the same between the subject who has reduced the function of the part of the chest that contributes to breathing and the healthy person, Even if the change in the first bioelectrical impedance is the same, the subject whose function of the portion of the chest that contributes to breathing has decreased has a larger amount of diaphragm fluctuation than the healthy subject. 2 Change in bioelectrical impedance increases. In the above-described aspect of the present invention, when the function of a part contributing to breathing in the subject's chest is reduced, the second bioelectrical impedance is measured by focusing on the fact that the change in the second bioelectrical impedance is increased. Based on the value and the measured value of the first bioelectrical impedance, it is determined whether or not the function of the part contributing to respiration in the subject's chest is normal. This makes it possible to determine the characteristics of the subject's ventilation function (for example, the function of the part of the chest that contributes to breathing is reduced).

ここで、胸式で胸郭を広げる呼吸法(胸式呼吸)では、内外肋間筋などの呼吸筋の伸縮変化と肺の伸縮(膨張収縮)変化が同じ方向に作用するので、肺の生体電気インピーダンスが増加すれば胸部骨格筋の生体電気インピーダンスも増加し、肺の生体電気インピーダンスが減少すれば胸部骨格筋の生体電気インピーダンスも減少する。一方、腹式呼吸は胸郭の変化がほとんど見られない呼吸法なので、胸部骨格筋の生体電気インピーダンスはほとんど変化せず、肺の生体電気インピーダンスが呼吸に伴って大きく変化する。被験者の呼吸が、胸式呼吸および腹式呼吸の何れの場合であっても、肺の上部を含み、腹部を含まない体幹上部の第1生体電気インピーダンスは、吸気では増加方向に変化する一方、呼気では減少方向に変化するという具合である。   Here, in the breathing method (thoracic breathing) that expands the rib cage with the chest type, the change in the expansion and contraction of the respiratory muscles such as the internal and external intercostal muscles and the change in the expansion and contraction (expansion and contraction) of the lungs act in the same direction. Increases the bioelectrical impedance of the thoracic skeletal muscle, and when the pulmonary bioelectrical impedance decreases, the bioelectrical impedance of the thoracic skeletal muscle also decreases. On the other hand, since abdominal breathing is a breathing method in which changes in the rib cage are hardly seen, the bioelectrical impedance of the thoracic skeletal muscle hardly changes, and the bioelectrical impedance of the lung changes greatly with respiration. Whether the subject's breathing is chest breathing or abdominal breathing, the first bioelectrical impedance of the upper trunk including the upper part of the lung and not including the abdomen changes in an increasing direction with inspiration. In exhalation, it changes in a decreasing direction.

一方、腹式呼吸の特徴は、腹部骨格筋の伸縮により内臓組織と伴に横隔膜を上下させる点にある。より具体的には、腹式呼吸の呼気時のみ、腹筋を緊張させて内臓組織と伴に横隔膜を押し上げ上昇させることで、内臓組織と腹部骨格筋との並列部の生体電気インピーダンスが上昇する。このとき、肺組織の生体電気インピーダンスは減少する。このため、横隔膜から上部の生体電気インピーダンスの減少を横隔膜から下部の生体電気インピーダンスの増加が打ち消すように作用する。胸式呼吸では、そのようなことはない。   On the other hand, the feature of abdominal breathing is that the diaphragm is moved up and down along with the visceral tissue by expansion and contraction of the abdominal skeletal muscle. More specifically, only during exhalation of abdominal breathing, the abdominal muscles are tensioned and the diaphragm is pushed up and lifted together with the visceral tissue, thereby increasing the bioelectrical impedance of the parallel portion of the visceral tissue and the abdominal skeletal muscle. At this time, the bioelectrical impedance of the lung tissue decreases. For this reason, it acts so that the increase in the bioelectric impedance in the upper part from the diaphragm cancels the increase in the bioelectric impedance in the lower part from the diaphragm. This is not the case with chest breathing.

腹式呼吸の呼気における第2生体電気インピーダンスの変化は、第1生体電気インピーダンスの変化とは異なる(変化を示す波形が異なる)のに対して、吸気における第2生体電気インピーダンスの変化は、胸式呼吸および腹式呼吸の何れの場合であっても、第1生体電気インピーダンスの変化と同じである(変化を示す波形が同じである)。すなわち、第2生体電気インピーダンスの経時的変化を示す波形には、呼気時の腹式呼吸に起因した波形歪みが含まれるので、吸気における第1生体電気インピーダンスの変化と第2生体電気インピーダンスの変化とに基づいて、被験者の胸部のうち呼吸に寄与する部位の機能が正常であるか否かを判定することが好ましい。   The change in the second bioelectrical impedance in the breath of abdominal breathing is different from the change in the first bioelectrical impedance (the waveform indicating the change is different), whereas the change in the second bioelectrical impedance in the inspiration It is the same as the change of the first bioelectrical impedance in both cases of the respiration and the abdominal respiration (the waveforms indicating the change are the same). That is, since the waveform indicating the change over time of the second bioelectrical impedance includes waveform distortion caused by abdominal breathing during expiration, the change in the first bioelectrical impedance and the change in the second bioelectrical impedance during inspiration Based on the above, it is preferable to determine whether or not the function of the portion of the subject's chest that contributes to breathing is normal.

この好適な態様として、第1生体電気インピーダンスの測定値の振幅基準レベルを示す第1センタリング値、および、第2生体電気インピーダンスの測定値の振幅基準レベルを示す第2センタリング値を生成するセンタリング値生成部(170)と、第1生体電気インピーダンスの測定値の第1センタリング値に対する相対値である第1相対値を求める第1相対値算出部(170)と、第2生体電気インピーダンスの測定値の第2センタリング値に対する相対値である第2相対値を求める第2相対値算出部(170)と、をさらに備え、判定部は、被験者の1呼吸ごとに、当該1呼吸の吸気における第1相対値と第2相対値との比に基づいて、被験者の胸部のうち呼吸に寄与する部位の機能が正常であるか否かを判定する。第1センタリング値とは、呼吸に伴う第1生体電気インピーダンスの経時的変化を示す波形から、呼吸に起因した情報を抽出するための振幅基準レベルである。また、第2センタリング値とは、呼吸に伴う第2生体電気インピーダンスの経時的変化を示す波形から、呼吸に起因した情報を抽出するための振幅基準レベルである。上述の態様では、呼吸に起因した情報を示す第1相対値および第2相対値に基づいた判定処理を行うことで、被験者の換気機能の特性をより正確に判定できるという利点がある。   As a preferred embodiment, a centering value for generating a first centering value indicating the amplitude reference level of the measurement value of the first bioelectrical impedance and a second centering value indicating the amplitude reference level of the measurement value of the second bioelectrical impedance. A generation unit (170), a first relative value calculation unit (170) that obtains a first relative value that is a relative value of the measured value of the first bioelectrical impedance to the first centering value, and a measured value of the second bioelectrical impedance A second relative value calculation unit (170) that obtains a second relative value that is a relative value to the second centering value, and the determination unit performs the first inspiration of the one breath for each breath of the subject. Based on the ratio between the relative value and the second relative value, it is determined whether or not the function of the part contributing to breathing in the subject's chest is normal. A 1st centering value is an amplitude reference level for extracting the information resulting from respiration from the waveform which shows the time-dependent change of the 1st bioelectrical impedance accompanying respiration. The second centering value is an amplitude reference level for extracting information resulting from respiration from a waveform indicating a change with time of the second bioelectric impedance accompanying respiration. In the above-described aspect, there is an advantage that the characteristics of the subject's ventilation function can be more accurately determined by performing the determination process based on the first relative value and the second relative value indicating information resulting from respiration.

より具体的な態様として、判定部は、第1相対値のピーク値に対する第2相対値のピーク値の割合が所定の閾値以上である場合は、被験者の胸部のうち呼吸に寄与する部位の機能は正常ではないと判定する一方、閾値を下回る場合は、被験者の胸部のうち呼吸に寄与する部位の機能は正常であると判定する。これにより、被験者の胸部のうち呼吸に寄与する部位の機能が正常であるか否かが正確に判定されるという具合である。   As a more specific aspect, when the ratio of the peak value of the second relative value to the peak value of the first relative value is equal to or greater than a predetermined threshold, the determination unit functions as a part of the subject's chest that contributes to breathing. On the other hand, if it is less than the threshold, it is determined that the function of the part contributing to breathing in the subject's chest is normal. Thus, it is accurately determined whether or not the function of the part contributing to breathing in the subject's chest is normal.

本発明に係る換気特性判定装置の態様として、第1生体電気インピーダンスの測定値と、第1センタリング値とが等しくなるゼロクロスタイミングを抽出するゼロクロスタイミング抽出部(170)をさらに備え、生体電気インピーダンス測定部は、所定の周期でサンプリングタイミングに到達するたびに、第1生体電気インピーダンスおよび第2生体電気インピーダンスを測定し、センタリング値生成部は、所定数のサンプリングタイミングの各々における第1生体電気インピーダンスの測定値に基づいて第1センタリング値を生成する一方、ゼロクロスタイミング抽出部で抽出されたゼロクロスタイミングにおける第2生体電気インピーダンスの測定値に基づいて第2センタリング値を生成する。   The aspect of the ventilation characteristic determination device according to the present invention further includes a zero cross timing extraction unit (170) that extracts a zero cross timing at which the measured value of the first bioelectrical impedance is equal to the first centering value, and the bioelectrical impedance measurement is performed. The unit measures the first bioelectrical impedance and the second bioelectrical impedance each time the sampling timing is reached at a predetermined cycle, and the centering value generating unit is configured to measure the first bioelectrical impedance at each of the predetermined number of sampling timings. The first centering value is generated based on the measurement value, and the second centering value is generated based on the measurement value of the second bioelectrical impedance at the zero cross timing extracted by the zero cross timing extraction unit.

被験者の呼吸が胸式呼吸であっても腹式呼吸であっても、第1生体電気インピーダンスの変化を示す波形は略正弦波状となる。体動などによる影響で、第1生体電気インピーダンスの変化を示す波形が乱れても、それに応じた第1センタリング値(第1生体電気インピーダンスの変化を示す波形から、呼吸に起因した情報を抽出するための振幅基準レベル)が得られるように、センタリング値生成部は、所定数のサンプリングタイミングの各々における第1生体電気インピーダンスの測定値に基づいて第1センタリング値を生成する。より具体的には、センタリング値生成部は、サンプリングタイミングごとに、当該サンプリングタイミングよりも所定の時間長だけ前の時点を始点とし、当該サンプリングタイミングを終点とするセンタリング期間内の複数のサンプリングタイミングの各々における第1生体電気インピーダンスの測定値を用いた移動平均処理を行い、その結果に基づいて、当該サンプリングタイミングにおける第1センタリング値を求める。これにより、体動などによる影響で、第1生体電気インピーダンスの変化を示す波形が乱れても、それに応じた第1センタリング値を精度良く生成できる。上記センタリング期間の時間長は、当該サンプリングタイミングにおける被験者の呼吸速度に応じて可変に設定される。なお、ここでの「移動平均処理」には、重み付けのない平均処理だけでなく、重み付けのある平均処理も含まれる。例えば各サンプリングタイミングにおける周波数の相違に応じた重み付けがされたうえで、平均処理が行われる態様であってもよい。   Whether the subject's breathing is chest breathing or abdominal breathing, the waveform indicating the change in the first bioelectrical impedance is substantially sinusoidal. Even if the waveform indicating the change in the first bioelectrical impedance is disturbed due to the influence of body movement or the like, the first centering value corresponding to that is extracted (information resulting from respiration from the waveform indicating the change in the first bioelectrical impedance). The centering value generation unit generates the first centering value based on the measured value of the first bioelectrical impedance at each of the predetermined number of sampling timings. More specifically, for each sampling timing, the centering value generator generates a plurality of sampling timings within a centering period starting from a time point that is a predetermined time length before the sampling timing and ending with the sampling timing. A moving average process using the measured value of the first bioelectric impedance in each is performed, and based on the result, the first centering value at the sampling timing is obtained. Thereby, even if the waveform indicating the change in the first bioelectrical impedance is disturbed due to the influence of body movement or the like, the first centering value corresponding to the waveform can be generated with high accuracy. The time length of the centering period is variably set according to the breathing rate of the subject at the sampling timing. Here, the “moving average process” includes not only an unweighted average process but also an average process with a weight. For example, the averaging process may be performed after weighting according to the difference in frequency at each sampling timing.

一方、腹式呼吸に伴う第2生体電気インピーダンスの変化は、第1生体電気インピーダンスの変化とは異なる態様(非正弦波状)となるので、第1センタリング値を求める場合と同様に、所定数のサンプリングタイミングの各々における第2生体電気インピーダンスの測定値を用いた移動平均処理を行っても、これらの振幅基準レベル、すなわち、第2生体電気インピーダンスの変化を示す波形から、呼吸に起因した情報を抽出するための振幅基準レベルを精度良く求めることは困難である。   On the other hand, the change in the second bioelectric impedance due to the abdominal breathing has a different form (non-sinusoidal) from the change in the first bioelectrical impedance. Even if the moving average processing using the measurement value of the second bioelectric impedance at each sampling timing is performed, information derived from respiration is obtained from the amplitude reference level, that is, the waveform indicating the change in the second bioelectric impedance. It is difficult to accurately obtain an amplitude reference level for extraction.

そこで、上記態様では、センタリング値生成部は、ゼロクロスタイミング抽出部で抽出されたゼロクロスタイミングにおける第2生体電気インピーダンスの測定値に基づいて、当該サンプリングタイミングにおける第2センタリング値を生成する。具体的には、センタリング値生成部は、サンプリングタイミングごとに、当該サンプリングタイミングがゼロクロスタイミングであるか否かを判定し、当該サンプリングタイミングがゼロクロスタイミングである場合は、当該サンプリングタイミングにおける第2生体電気インピーダンスの測定値に基づいて、当該サンプリングタイミングにおける第2センタリング値を生成する一方、当該サンプリングタイミングがゼロクロスタイミングでない場合は、当該サンプリングタイミングの直前のサンプリングタイミングで生成した第2センタリング値を、当該サンプリングタイミングにおける第2センタリング値として採用する。これにより、第2生体電気インピーダンスの測定値の振幅基準レベルを精度良く求めることができる。   Therefore, in the above aspect, the centering value generation unit generates the second centering value at the sampling timing based on the measurement value of the second bioelectrical impedance at the zero cross timing extracted by the zero cross timing extraction unit. Specifically, for each sampling timing, the centering value generation unit determines whether or not the sampling timing is a zero-cross timing. If the sampling timing is the zero-cross timing, the second bioelectric at the sampling timing is determined. Based on the measured impedance value, the second centering value at the sampling timing is generated. On the other hand, if the sampling timing is not zero cross timing, the second centering value generated at the sampling timing immediately before the sampling timing is The second centering value at the timing is adopted. Thereby, the amplitude reference level of the measured value of the second bioelectrical impedance can be obtained with high accuracy.

本発明に係る実施形態の生体測定装置の電気的構成について示すブロック図である。It is a block diagram shown about the electrical constitution of the biometric device of the embodiment concerning the present invention. 生体測定装置の外観例を示す斜視図である。It is a perspective view which shows the example of an external appearance of a biometric apparatus. 生体測定装置の電極配置を示す説明図である。It is explanatory drawing which shows electrode arrangement | positioning of a biometric apparatus. 電流電極の選択と電圧電極の選択とを説明するための説明図である。It is explanatory drawing for demonstrating selection of a current electrode and selection of a voltage electrode. 生体測定装置の動作内容を示すフローチャートである。It is a flowchart which shows the operation | movement content of a biometric apparatus. 体幹を構成する組織の概略を示す模式図である。It is a schematic diagram which shows the outline of the structure | tissue which comprises a trunk. 体幹の生体電気インピーダンスの等価回路を示す回路図である。It is a circuit diagram which shows the equivalent circuit of the bioelectrical impedance of a trunk. 呼吸と生体電気インピーダンスの変化の関係を説明する説明図である。It is explanatory drawing explaining the relationship between respiration and the change of bioelectrical impedance. 腹式呼吸における生体電気インピーダンスの変化を示す図である。It is a figure which shows the change of the bioelectrical impedance in abdominal respiration. 胸式呼吸における生体電気インピーダンスの変化を示す図である。It is a figure which shows the change of the bioelectrical impedance in chest type respiration. 腹式呼吸に伴う胸部および腹部の周囲径変化を示す図である。It is a figure which shows the circumference change of the chest and abdomen accompanying abdominal respiration. 胸式呼吸に伴う胸部および腹部の周囲径変化を示す図である。It is a figure which shows the circumference diameter change of the chest and abdomen accompanying chest type respiration. 胸部の周囲径変化と腹部の周囲径変化との比と、第1相対値および第2相対値との関係を示す相関図である。It is a correlation diagram showing the relationship between the ratio between the change in the circumference of the chest and the change in the circumference of the abdomen, and the first relative value and the second relative value. 呼吸解析処理の処理内容を示すフローチャートである。It is a flowchart which shows the processing content of a respiration analysis process. 第1センタリング処理の処理内容を示すフローチャートである。It is a flowchart which shows the processing content of a 1st centering process. 呼吸タイミング抽出処理の処理内容を示すフローチャートである。It is a flowchart which shows the processing content of a respiration timing extraction process. 呼吸タイミング抽出処理の処理内容を示すフローチャートである。It is a flowchart which shows the processing content of a respiration timing extraction process. 呼吸スピード判別フラッグ設定処理の処理内容を示すフローチャートである。It is a flowchart which shows the processing content of a breathing speed discrimination | determination flag setting process. 第1センタリング値抽出処理の処理内容を示すフローチャートである。It is a flowchart which shows the processing content of a 1st centering value extraction process. 第2センタリング値の生成方法を概念的に説明するための模式図である。It is a schematic diagram for demonstrating notionally the production | generation method of a 2nd centering value. 第2相対値算出処理の処理内容を示すフローチャートである。It is a flowchart which shows the processing content of a 2nd relative value calculation process. 腹式呼吸における第1相対値と第2相対値との経時的変化を示す図である。It is a figure which shows the time-dependent change of the 1st relative value and 2nd relative value in abdominal respiration. ΔRib/ΔAb推定演算処理の処理内容を示すフローチャートである。It is a flowchart which shows the processing content of (DELTA) Rib / (DELTA) Ab estimation calculation processing. ΔRib/ΔAb推定演算処理の処理内容を示すフローチャートである。It is a flowchart which shows the processing content of (DELTA) Rib / (DELTA) Ab estimation calculation processing. ΔRib/ΔAb推定演算処理の演算結果を示す図である。It is a figure which shows the calculation result of (DELTA) Rib / (DELTA) Ab estimation calculation processing. 呼吸レベル抽出処理の処理内容を示すフローチャートである。It is a flowchart which shows the processing content of a respiration level extraction process. 表示部での表示態様を示す図である。It is a figure which shows the display mode in a display part. 呼吸レベル表示処理の処理内容を示すフローチャートである。It is a flowchart which shows the processing content of a respiration level display process. 呼吸レベルと1回換気量との間の関係を示す相関図である。It is a correlation diagram which shows the relationship between a respiration level and a tidal volume. 換気特性判定/表示処理の処理内容を示すフローチャートである。It is a flowchart which shows the processing content of ventilation characteristic determination / display processing. 胸部に病歴が無い被験者の呼吸に伴う第1相対値および第2相対値の経時的変化を示す図である。It is a figure which shows the time-dependent change of the 1st relative value and 2nd relative value accompanying the respiration of the test subject who has no medical history in a chest. 胸部に病歴を有する被験者の呼吸に伴う第1相対値および第2相対値の経時的変化を示す図である。It is a figure which shows the time-dependent change of the 1st relative value and 2nd relative value accompanying the respiration of the test subject who has a medical history in the chest. 換気能バランスの表示態様を示す図である。It is a figure which shows the display mode of ventilation capacity balance.

<1.実施形態>
<1−1:生体測定装置の構成>
図1は、本発明の実施形態に係る生体測定装置1の構成を示すブロック図である。この生体測定装置1は、生体の状態を測定するものであるが、その機能の一部は、被験者の換気機能の特性を判定する換気特性判定装置としての役割を担う。
生体測定装置1は、体重を測定するとともに装置全体の動作を管理する管理部100と、被験者の各部位の生体電気インピーダンスを測定する生体電気インピーダンス測定部200とを備える。管理部100は、体重計110、第1記憶部120、第2記憶部130、音声処理部140、スピーカ145、入力部150、並びに表示部160を備える。これらの構成要素は、バスを介してCPU(Central Processing Unit)170と接続されている。CPU170は、装置全体を制御する制御中枢として機能する。なお、CPU170は図示せぬクロック信号発生回路からクロック信号の供給を受けて動作する。また、各構成要素には図示せぬ電源スイッチがオン状態になると、電源回路から電源が供給される。
<1. Embodiment>
<1-1: Configuration of Biometric Apparatus>
FIG. 1 is a block diagram showing a configuration of a biometric apparatus 1 according to an embodiment of the present invention. This living body measuring apparatus 1 measures the state of a living body, but a part of its function plays a role as a ventilation characteristic determining apparatus that determines the characteristics of a subject's ventilation function.
The biometric apparatus 1 includes a management unit 100 that measures body weight and manages the operation of the entire apparatus, and a bioelectrical impedance measurement unit 200 that measures bioelectrical impedance of each part of the subject. The management unit 100 includes a weight scale 110, a first storage unit 120, a second storage unit 130, a sound processing unit 140, a speaker 145, an input unit 150, and a display unit 160. These components are connected to a CPU (Central Processing Unit) 170 via a bus. The CPU 170 functions as a control center that controls the entire apparatus. The CPU 170 operates by receiving a clock signal from a clock signal generation circuit (not shown). Further, when a power switch (not shown) is turned on to each component, power is supplied from the power supply circuit.

体重計110は、被験者の体重を測定し、その測定した体重データを、バスを介してCPU170に出力する。第1記憶部120は、不揮発性のメモリであって、例えばROM(Read Only Memory)で構成される。第1記憶部120には、装置全体を制御する制御プログラムが記憶されている。CPU170は、制御プログラムにしたがって所定の演算を実行する。   The weight scale 110 measures the weight of the subject, and outputs the measured weight data to the CPU 170 via the bus. The first storage unit 120 is a non-volatile memory, and is composed of, for example, a ROM (Read Only Memory). The first storage unit 120 stores a control program for controlling the entire apparatus. CPU 170 executes a predetermined calculation according to the control program.

第2記憶部130は、揮発性のメモリであり、例えば、DRAM(Dynamic Random Access Memory)等によって構成される。第2記憶部130はCPU170の作業領域として機能し、CPU170が所定の演算を実行する際にデータを記憶する。また、音声処理部140は、CPU170の制御の下、音声データをDA変換して得た音声信号を増幅してスピーカ145に出力する。スピーカ145は増幅した音声信号を振動に変換して放音する。これによって、呼吸方法の指導などのアドバイス情報を音によって被験者に報知することができる。   The second storage unit 130 is a volatile memory, and includes, for example, a DRAM (Dynamic Random Access Memory) or the like. The second storage unit 130 functions as a work area of the CPU 170, and stores data when the CPU 170 executes a predetermined calculation. In addition, the audio processing unit 140 amplifies an audio signal obtained by DA converting audio data under the control of the CPU 170 and outputs the amplified audio signal to the speaker 145. The speaker 145 converts the amplified audio signal into vibration and emits the sound. Thus, advice information such as guidance on the breathing method can be notified to the subject by sound.

入力部150は、各種のスイッチから構成され、被験者がスイッチを操作すると、身長、年齢、及び性別といった情報が入力される。表示部160は、体重や呼吸の種別といった測定結果や、腹式呼吸に導くための呼気と吸気のリズムやパターンなどのアドバイス情報を知らせる機能、あるいは被験者に各種の情報の入力を促すメッセージを表示する機能を有する。表示部160は、例えば、液晶表示装置などで構成される。   The input unit 150 includes various switches, and when a subject operates the switches, information such as height, age, and sex is input. The display unit 160 displays a measurement result such as weight and type of breath, a function for notifying advice information such as a rhythm and a pattern of exhalation and inspiration for leading to abdominal breathing, or a message prompting the subject to input various information. It has the function to do. The display unit 160 is configured by, for example, a liquid crystal display device.

次に、生体電気インピーダンス測定部200は、被験者(人体)の生体電気インピーダンスを測定する。生体電気インピーダンス測定部200は、交流電流出力回路210、基準電流検出回路220、電位差検出回路230、A/D変換器240、電極切換回路251及び252を備える。
交流電流出力回路210は基準電流Irefを生成する手段である。交流電流出力回路210は、基準電流Irefの実効値が予め定められた値となるように、当該基準電流Irefを生成する。基準電流検出回路220は、被測定対象に流れる基準電流Irefの大きさを検出して電流データDiとしてCPU170に出力するとともに、被験者に基準電流Irefを通電する。この場合、電極切換回路252は、電流電極X1〜X4の中から2つを選択して電流を供給する。
さらに、電位差検出回路230は、電圧電極Y1〜Y4の中から選択された2つの電圧電極の間の電位差を検出して電位差信号ΔVを生成する。A/D変換器240は電位差信号ΔVをアナログ信号からデジタル信号に変換し電圧データDvとしてCPU170に出力する。CPU170は電圧データDvと電流データDiとに基づいて生体電気インピーダンスZ(=Dv/Di)を計算する。
Next, the bioelectrical impedance measuring unit 200 measures the bioelectrical impedance of the subject (human body). The bioelectrical impedance measurement unit 200 includes an alternating current output circuit 210, a reference current detection circuit 220, a potential difference detection circuit 230, an A / D converter 240, and electrode switching circuits 251 and 252.
The alternating current output circuit 210 is a means for generating a reference current Iref. The alternating current output circuit 210 generates the reference current Iref so that the effective value of the reference current Iref becomes a predetermined value. The reference current detection circuit 220 detects the magnitude of the reference current Iref flowing through the object to be measured and outputs it as current data Di to the CPU 170 and energizes the subject with the reference current Iref. In this case, the electrode switching circuit 252 selects two of the current electrodes X1 to X4 and supplies current.
Further, the potential difference detection circuit 230 detects a potential difference between two voltage electrodes selected from the voltage electrodes Y1 to Y4 and generates a potential difference signal ΔV. The A / D converter 240 converts the potential difference signal ΔV from an analog signal to a digital signal and outputs it to the CPU 170 as voltage data Dv. The CPU 170 calculates a bioelectric impedance Z (= Dv / Di) based on the voltage data Dv and the current data Di.

第1記憶部120は、各種データを予め記憶することができる。たとえば、各部位の生体電気インピーダンスを変数として体脂脂肪率や筋肉量を算出するための相関式又は相関テーブルが記憶されている。
CPU170は、体重、被験者の各種の部位の生体電気インピーダンス(例えば、上肢生体電気インピーダンス、下肢生体電気インピーダンス、体幹生体電気インピーダンス)、を演算し、かつ、各種の入出力、測定、演算等について制御する。なお、生体電気インピーダンスなどに基づいて、内臓脂肪/皮下脂肪、内臓脂肪量、皮下脂肪率、皮下脂肪量、全身の脂肪率、身体の各部位の脂肪率(上肢脂肪率、下肢脂肪率、体幹脂肪率など)を演算することもできる。
The first storage unit 120 can store various data in advance. For example, a correlation equation or a correlation table for calculating body fat fat percentage and muscle mass using the bioelectrical impedance of each part as a variable is stored.
The CPU 170 calculates body weight, bioelectrical impedance (for example, upper limb bioelectrical impedance, lower limb bioelectrical impedance, trunk bioelectrical impedance) of various parts of the subject, and various input / output, measurement, calculation, and the like. Control. Based on bioelectrical impedance, etc., visceral fat / subcutaneous fat, visceral fat mass, subcutaneous fat percentage, subcutaneous fat mass, whole body fat percentage, fat percentage of each part of body (upper limb fat percentage, lower limb fat percentage, body Stem fat percentage etc.) can also be calculated.

図2に、生体測定装置1の外観例を示す。生体測定装置1は、L字型の形状をしており、台座部20の上に柱状の筐体部30を備える。台座部20には、左足用の電流電極X1及び電圧電極Y1と、右足用の電流電極X2及び電圧電極Y2が設けられている。また、筐体部30の上部には、表示部160が設けられている。この表示部160は、タッチパネルで構成されており、入力部150としても機能する。さらに、筐体部30の左右の側面には、左手用の電極部30Lと右手用の電極部30Rが設けられている。   In FIG. 2, the example of an external appearance of the biometric apparatus 1 is shown. The biometric device 1 has an L-shape and includes a columnar casing 30 on the pedestal 20. The pedestal 20 is provided with a current electrode X1 and a voltage electrode Y1 for the left foot, and a current electrode X2 and a voltage electrode Y2 for the right foot. A display unit 160 is provided on the upper portion of the housing unit 30. The display unit 160 is configured with a touch panel and also functions as the input unit 150. Further, left and right electrode portions 30 </ b> L and right hand electrode portions 30 </ b> R are provided on the left and right side surfaces of the housing portion 30.

図3は筐体部30の上部を拡大した拡大図である。この図に示すように、左手用の電極部30Lは電流電極X3及び電圧電極Y3を備え、右手用の電極部30Rは電流電極X4及び電圧電極Y4を備える。被験者は、台座30の上に立ち、左右の手を下げた状態で電極部30L及び電極部30Rを握ることによって、測定を行う。   FIG. 3 is an enlarged view in which the upper part of the housing part 30 is enlarged. As shown in this figure, the left-hand electrode portion 30L includes a current electrode X3 and a voltage electrode Y3, and the right-hand electrode portion 30R includes a current electrode X4 and a voltage electrode Y4. The subject performs measurement by standing on the pedestal 30 and holding the electrode part 30L and the electrode part 30R with the left and right hands lowered.

電極切換回路251及び252は、CPU170の制御の下、両手及び両足に装着される8個の電極を選択する。この8個の電極を適宜選択することによって、人体の所定の部位における生体電気インピーダンスZを計測することが可能となる。例えば、図4(A)に示すように基準電流Irefを左足用の電流電極X1と左手用の電流電極X3との間に供給し、左足用の電圧電極Y1と左手用の電圧電極Y3との間で電位差を計測すれば、全身の生体電気インピーダンスを計測することができる。なお、基準電流Irefを電流電極X2及びX4の間に流し、電圧電極Y2及びY4の間の電位差を計測しても全身の生体電気インピーダンスを計測することができる。さらに、図4(K)に示すように両掌を短絡させ、両足を短絡させ、両掌から両足までの生体電気インピーダンスを全身の生体電気インピーダンスとして測定してもよい。   The electrode switching circuits 251 and 252 select eight electrodes to be worn on both hands and both feet under the control of the CPU 170. By appropriately selecting these eight electrodes, it is possible to measure the bioelectrical impedance Z at a predetermined part of the human body. For example, as shown in FIG. 4A, the reference current Iref is supplied between the current electrode X1 for the left foot and the current electrode X3 for the left hand, and the voltage electrode Y1 for the left foot and the voltage electrode Y3 for the left hand If the potential difference is measured between them, the bioelectrical impedance of the whole body can be measured. Note that the bioelectric impedance of the whole body can be measured even if the reference current Iref is passed between the current electrodes X2 and X4 and the potential difference between the voltage electrodes Y2 and Y4 is measured. Furthermore, as shown in FIG. 4 (K), both palms may be short-circuited, both feet may be short-circuited, and the bioelectric impedance from both palms to both feet may be measured as the bioelectric impedance of the whole body.

また、図4(B)に示すように基準電流Irefを右足用の電流電極X2と右手用の電流電極X4との間に供給し、右足用の電圧電極Y2と左足用の電圧電極Y1との間で電位差を計測すれば、右下肢の生体電気インピーダンスZを計測することができる。また、図4(C)に示すように基準電流Irefを左足用の電流電極X1と左手用の電流電極X3との間に供給し、左足用の電圧電極Y1と右足用の電圧電極Y2との間で電位差を計測すれば、左下肢の生体電気インピーダンスZを計測することができる。   Further, as shown in FIG. 4B, the reference current Iref is supplied between the right foot current electrode X2 and the right hand current electrode X4, and the right foot voltage electrode Y2 and the left foot voltage electrode Y1 are connected. If the potential difference is measured between them, the bioelectrical impedance Z of the right lower limb can be measured. Further, as shown in FIG. 4C, the reference current Iref is supplied between the left foot current electrode X1 and the left hand current electrode X3, and the left foot voltage electrode Y1 and the right foot voltage electrode Y2 are connected. If the potential difference is measured between them, the bioelectrical impedance Z of the left lower limb can be measured.

また、図4(D)に示すように基準電流Irefを右手用の電流電極X4と右足用の電流電極X2との間に供給し、右手用の電圧電極Y4と左手用の電圧電極Y3との間の電位差を計測すれば、右上肢(体幹上部)の生体電気インピーダンスZを計測することができる。ただし、これに限らず、基準電流Irefを左手用の電流電極X3と右手用の電流電極X4との間に供給し、右手用の電圧電極Y4と右足用の電圧電極Y2との間の電位差を計測することでも、右上肢(体幹上部)の生体電気インピーダンスZを計測することができる。   Further, as shown in FIG. 4D, the reference current Iref is supplied between the right-hand current electrode X4 and the right-foot current electrode X2, and the right-hand voltage electrode Y4 and the left-hand voltage electrode Y3 are connected. By measuring the potential difference between them, the bioelectrical impedance Z of the upper right limb (upper trunk) can be measured. However, the present invention is not limited thereto, and the reference current Iref is supplied between the left-hand current electrode X3 and the right-hand current electrode X4, and the potential difference between the right-hand voltage electrode Y4 and the right-foot voltage electrode Y2 is determined. The bioelectric impedance Z of the upper right limb (upper trunk) can also be measured by measuring.

また、図4(E)に示すように基準電流Irefを左手用の電流電極X3と左足用の電流電極X1との間に供給し、右手用の電圧電極Y4と左手用の電圧電極Y3との間で電位差を計測すれば、左上肢(体幹上部)の生体電気インピーダンスZを計測することができる。ただし、これに限らず、基準電流Irefを左手用の電流電極X3と右手用の電流電極X4との間に供給し、左手用の電圧電極Y3と左足用の電圧電極Y1との間の電位差を計測することでも、左上肢(体幹上部)の生体電気インピーダンスZを計測することができる。   Further, as shown in FIG. 4E, the reference current Iref is supplied between the left-hand current electrode X3 and the left-foot current electrode X1, and the right-hand voltage electrode Y4 and the left-hand voltage electrode Y3 are connected. If the potential difference is measured between them, the bioelectrical impedance Z of the left upper limb (upper trunk) can be measured. However, the reference current Iref is supplied between the left-hand current electrode X3 and the right-hand current electrode X4, and the potential difference between the left-hand voltage electrode Y3 and the left-foot voltage electrode Y1 is not limited thereto. The bioelectrical impedance Z of the left upper limb (upper trunk) can also be measured by measuring.

また、図4(F)に示すように基準電流Irefを右手用の電流電極X4と左手用の電流電極X3との間に供給し、右手用の電圧電極Y4と左手用の電圧電極Y3との間で電位差を計測すれば、両掌間の生体電気インピーダンスZを計測することができる。ここで、体幹を体幹上部と体幹中部に分けた場合、左上肢、右上肢、及び掌間の生体電気インピーダンスは、いずれも体幹上部が含まれる。このため、左上肢、右上肢、及び掌間の生体電気インピーダンスを体幹上部の生体電気インピーダンスとして取り扱うことも可能である。   Further, as shown in FIG. 4F, the reference current Iref is supplied between the right-hand current electrode X4 and the left-hand current electrode X3, and the right-hand voltage electrode Y4 and the left-hand voltage electrode Y3 are connected. If the potential difference is measured between the two, the bioelectric impedance Z between both palms can be measured. Here, when the trunk is divided into the upper trunk and the middle trunk, the bioelectrical impedance between the left upper limb, the upper right limb, and the palm includes the upper trunk. For this reason, it is also possible to handle the bioelectric impedance between the left upper limb, the upper right limb, and the palm as the bioelectric impedance of the upper trunk.

さらに、図4(G)に示すように基準電流Irefを左手用の電流電極X3と左足用の電流電極X1との間に供給し、右手用の電圧電極X4と右足用の電圧電極Y2との間で電位差を計測すれば、体幹中部の生体電気インピーダンスを計測することができる。また、図4(H)に示すように基準電流Irefを右手用の電流電極X4と右足用の電流電極X2との間に供給し、左手用の電圧電極Y3と左足用の電圧電極Y1との間で電位差を計測すれば、体幹中部の生体電気インピーダンスを計測することができる。また、図4(I)に示すように基準電流Irefを右手用の電流電極X4と左足用の電流電極X1との間に供給し、左手用の電圧電極Y3と右足用の電圧電極Y2との間で電位差を計測すれば、体幹中部を斜めに横切る生体電気インピーダンスを計測することができる。図4(J)に示すように基準電流Irefを右足用の電流電極X2と左手用の電流電極X3との間に供給し、右手用の電圧電極Y4と左足用の電圧電極Y1との間で電位差を計測すれば、体幹中部を斜めに横切る生体電気インピーダンスを計測することができる。   Further, as shown in FIG. 4G, the reference current Iref is supplied between the left hand current electrode X3 and the left foot current electrode X1, and the right hand voltage electrode X4 and the right foot voltage electrode Y2 are connected. If the potential difference is measured between the two, the bioelectrical impedance in the middle of the trunk can be measured. Further, as shown in FIG. 4H, the reference current Iref is supplied between the right-hand current electrode X4 and the right-foot current electrode X2, and the left-hand voltage electrode Y3 and the left-foot voltage electrode Y1 are connected. If the potential difference is measured between the two, the bioelectrical impedance in the middle of the trunk can be measured. Further, as shown in FIG. 4 (I), the reference current Iref is supplied between the current electrode X4 for the right hand and the current electrode X1 for the left foot, and the voltage electrode Y3 for the left hand and the voltage electrode Y2 for the right foot By measuring the potential difference between the two, it is possible to measure the bioelectrical impedance that crosses the middle part of the trunk diagonally. As shown in FIG. 4J, the reference current Iref is supplied between the current electrode X2 for the right foot and the current electrode X3 for the left hand, and between the voltage electrode Y4 for the right hand and the voltage electrode Y1 for the left foot. By measuring the potential difference, bioelectrical impedance that crosses the middle of the trunk diagonally can be measured.

なお、体幹部の生体電気インピーダンスZxの測定方法は、上述した方法に限定されるものではなく、両手両足の電極のうち、基準電流Irefを供給する電極と電位差を検出する電極とを適宜選択することによって、手、足、あるいは全身といった人体の各部位の生体電気インピーダンスZを各々測定し、測定結果を加減算して体幹中部の生体電気インピーダンスZを算出すればよい。さらに、四肢以外に頭部の耳たぶなどに四肢のいずれかの代用として使用しても、体幹部の生体電気インピーダンスZxの測定は可能である。くわえて、体幹に接触電極を設ける場合には言うに及ばない。   The method for measuring the bioelectrical impedance Zx of the trunk is not limited to the above-described method, and an electrode for supplying a reference current Iref and an electrode for detecting a potential difference are appropriately selected from the electrodes of both hands and feet. Thus, the bioelectrical impedance Z of each part of the human body such as the hand, foot, or whole body is measured, and the bioelectrical impedance Z of the middle trunk is calculated by adding and subtracting the measurement results. Furthermore, the bioelectrical impedance Zx of the trunk can be measured even if it is used as a substitute for any of the four limbs in the earlobe of the head in addition to the four limbs. In addition, it goes without saying that a contact electrode is provided on the trunk.

<1−2:生体測定装置の動作>
図5は、生体測定装置1の動作を示すフローチャートである。まず、入力部150における電源スイッチ(図示省略)がオンされると、図示せぬ電力供給部から電気系統各部に電力を供給し、表示部160により身長を含む身体特定情報(身長、性別、年齢など)を入力するための画面を表示する(ステップS1)。
続いて、入力部150から身長、性別、年齢等が入力されると、体重計110により体重が測定され、CPU170は体重を取得する(ステップS2)。
<1-2: Operation of Biometric Device>
FIG. 5 is a flowchart showing the operation of the biometric apparatus 1. First, when a power switch (not shown) in the input unit 150 is turned on, power is supplied to each part of the electric system from a power supply unit (not shown), and body identification information including height (height, gender, age) is displayed on the display unit 160. Etc.) is displayed (step S1).
Subsequently, when height, sex, age, and the like are input from the input unit 150, the weight is measured by the weight scale 110, and the CPU 170 acquires the weight (step S2).

ステップS2の後、CPU170は、呼吸解析処理を実行する(ステップS3)。この処理では、CPU170は、被験者の呼吸が腹式呼吸なのか胸式呼吸なのかを判別可能な判別情報を求める。ステップS3の後、CPU170は、被験者の1呼吸ごとに、当該1呼吸における腹式呼吸および胸式呼吸の各々の大きさと余裕度とを表示する呼吸レベル表示処理を実行する(ステップS4)。また、CPU170は、被験者の1呼吸ごとに、当該1呼吸における被験者の換気機能の特性(換気特性)を判定するとともに、その判定結果を表示する換気特性判定/表示処理を実行する(ステップS5)。これらの詳細な内容については後述する。   After step S2, CPU 170 executes a breath analysis process (step S3). In this processing, the CPU 170 obtains discrimination information that can discriminate whether the subject's breathing is abdominal breathing or chest breathing. After step S3, the CPU 170 executes a respiration level display process for displaying the magnitude and margin of each of abdominal respiration and chest respiration in one respiration for each respiration of the subject (step S4). Further, for each breath of the subject, the CPU 170 determines the characteristics of the subject's ventilation function (ventilation characteristics) in one breath, and executes a ventilation characteristic determination / display process for displaying the determination result (step S5). . Details of these will be described later.

<1−3:呼吸解析の原理>
次に、呼吸解析の原理について説明する。図6は、体幹部の組織の概略を示す模式図である。図6に示すように、体幹部の組織は、横隔膜によって上下に分けられている。上部には、肺と、内外肋間筋などの胸部骨格筋とが形成されている。一方、下部には、内臓組織と、内外腹斜筋・腹横筋や腹直筋などからなる腹部骨格筋とが形成されている。
腹式呼吸および胸式呼吸のいずれの場合であっても、呼気時に横隔膜は上昇して肺が圧縮され、吸気時に横隔膜は下降して肺は伸長拡大する。胸式呼吸に無い腹式呼吸の特徴は、腹直筋や内外腹斜筋・腹横筋などの腹部呼吸筋の伸縮により内臓組織と供に横隔膜を上下させる点にある。
<1-3: Principle of respiratory analysis>
Next, the principle of respiratory analysis will be described. FIG. 6 is a schematic diagram showing an outline of the tissue of the trunk. As shown in FIG. 6, the trunk tissue is divided into upper and lower portions by a diaphragm. In the upper part, lungs and thoracic skeletal muscles such as internal and external intercostals are formed. On the other hand, a visceral tissue and abdominal skeletal muscles including internal and external oblique / lateral abdominal muscles and rectus abdominis muscles are formed in the lower part.
In either case of abdominal breathing or chest breathing, the diaphragm rises and the lungs are compressed during expiration, and the diaphragm descends and the lungs expand and expand during inspiration. A feature of abdominal breathing that does not exist in chest breathing is that the diaphragm is moved up and down along with the visceral tissues by expansion and contraction of abdominal respiratory muscles such as the rectus abdominis, internal and external oblique and transverse abdominal muscles.

ここで、体幹上部の第1生体電気インピーダンスZaと体幹中部の第2生体電気インピーダンスZbとは、図7に示す等価回路で表すことができる。図7に示すように、第1生体電気インピーダンスZaは、胸部骨格筋の生体電気インピーダンスZ1、および、肺の生体電気インピーダンスZ2の並列インピーダンスと、上肢骨格筋の生体電気インピーダンスZ3とが直列に接続されたものとなる。ここで、Z1およびZ2の並列インピーダンスは、肺の上葉部の生体インピーダンスに相当する。つまり、第1生体電気インピーダンスZaは、肺の上部を含み、腹部を含まない体幹上部の生体電気インピーダンスである。   Here, the first bioelectric impedance Za of the upper trunk and the second bioelectric impedance Zb of the middle trunk can be represented by an equivalent circuit shown in FIG. As shown in FIG. 7, the first bioelectrical impedance Za is a serial connection of the bioelectrical impedance Z1 of the thoracic skeletal muscle and the parallel impedance of the bioelectrical impedance Z2 of the lung and the bioelectrical impedance Z3 of the upper limb skeletal muscle. Will be. Here, the parallel impedance of Z1 and Z2 corresponds to the bioimpedance of the upper lobe of the lung. That is, the first bioelectrical impedance Za is the bioelectrical impedance of the upper trunk including the upper part of the lung and not including the abdomen.

また、図7に示すように、第2生体電気インピーダンスZbは、胸部骨格筋の生体電気インピーダンスZ4、および、肺の生体電気インピーダンスZ5の並列インピーダンスと、腹部骨格筋の生体インピーダンスZ6、および、内臓組織の生体インピーダンスZ7の並列インピーダンスとが直列に接続されたものとなる。ここで、Z4およびZ5の並列インピーダンスは、肺の中下葉部の生体インピーダンスに相当する。また、横隔膜の生体電気インピーダンスは、内臓組織に代表される生体電気インピーダンスZ7に含ませて考えることができる。つまり、第2生体電気インピーダンスZbは、肺の中下部および腹部を含む体幹中部の生体電気インピーダンスである。   Further, as shown in FIG. 7, the second bioelectric impedance Zb includes the bioelectric impedance Z4 of the thoracic skeletal muscle, the parallel impedance of the bioelectric impedance Z5 of the lung, the bioimpedance Z6 of the abdominal skeletal muscle, and the internal organs. The parallel impedance of the biological impedance Z7 of the tissue is connected in series. Here, the parallel impedance of Z4 and Z5 corresponds to the bioimpedance of the middle and lower lobe of the lung. Further, the bioelectric impedance of the diaphragm can be considered to be included in the bioelectric impedance Z7 typified by visceral tissue. That is, the second bioelectrical impedance Zb is the bioelectrical impedance of the middle trunk including the middle lower part and the abdomen of the lung.

次に、図8を参照して、呼吸と生体電気インピーダンスの変化との関係を説明する。呼吸に連動した第1生体電気インピーダンスZaの変化は、肺に絶縁性の高い空気が出入りすることによる電気的特質(電気導電性、1/体積抵抗率)の変化が主な原因であると考えられる。つまり、呼気(呼息)では肺組織中に含まれる空気量が減るため肺の生体電気インピーダンスZ2は減少方向に変化する(ΔZlu<0)。一方、吸気(吸息)では空気量が増加するため、肺の生体電気インピーダンスZ2は増加方向に変化する(ΔZlu>0)。   Next, with reference to FIG. 8, the relationship between respiration and changes in bioelectrical impedance will be described. Changes in the first bioelectric impedance Za linked to respiration are thought to be mainly due to changes in electrical characteristics (electrical conductivity, 1 / volume resistivity) due to the entry and exit of highly insulating air into and from the lungs. It is done. That is, in exhalation (exhalation), the amount of air contained in the lung tissue decreases, so that the bioelectric impedance Z2 of the lung changes in a decreasing direction (ΔZlu <0). On the other hand, in the inspiration (inhalation), the amount of air increases, so the bioelectric impedance Z2 of the lung changes in the increasing direction (ΔZlu> 0).

胸式で胸郭を広げる呼吸法(胸式呼吸)では、内外肋間筋などの呼吸骨格筋の伸縮変化と肺の伸縮変化が同じ方向に作用するので、肺の生体電気インピーダンスZ2が増加すれば胸部骨格筋の生体電気インピーダンスZ1も増加し、肺の生体電気インピーダンスZ2が減少すれば胸部骨格筋の生体電気インピーダンスZ1も減少する。一方、腹式呼吸は胸郭の変化がほとんど見られない呼吸法なので、胸部骨格筋の生体電気インピーダンスZ1はほとんど変化せず、肺の生体電気インピーダンスZ2が呼吸に伴って大きく変化する。なお、第1生体電気インピーダンスZaには、上肢骨格筋の生体電気インピーダンスZ3が含まれているが、上肢骨格筋は、呼吸に直接的に寄与する筋肉ではない。本実施形態では、被験者は、図2に示す測定装置の台座部20の上に立ち、左右の腕を下げた状態で30L及び30Rを握り計測を行うので、計測中に上肢骨格筋(Z3)が動くことは殆ど無い。図9および図10に示すように、被験者の呼吸が、胸式呼吸および腹式呼吸の何れの場合であっても、吸気では第1生体電気インピーダンスZaは増加方向に変化し、呼気では第1生体電気インピーダンスZaは減少方向に変化するという具合である。   In the breathing method (thoracic breathing) that expands the thorax with the chest type, the expansion and contraction of the respiratory skeletal muscles such as the internal and external intercostal muscles and the expansion and contraction of the lungs act in the same direction, so if the lung bioelectrical impedance Z2 increases, the chest The bioelectrical impedance Z1 of the skeletal muscle also increases, and if the bioelectrical impedance Z2 of the lung decreases, the bioelectrical impedance Z1 of the thoracic skeletal muscle also decreases. On the other hand, since abdominal breathing is a breathing method in which changes in the rib cage are hardly seen, the bioelectrical impedance Z1 of the thoracic skeletal muscle hardly changes, and the bioelectrical impedance Z2 of the lung changes greatly with respiration. The first bioelectric impedance Za includes the bioelectrical impedance Z3 of the upper limb skeletal muscle, but the upper limb skeletal muscle is not a muscle that directly contributes to respiration. In the present embodiment, the subject stands on the pedestal portion 20 of the measuring apparatus shown in FIG. 2 and grasps 30L and 30R in a state where the left and right arms are lowered, so the upper limb skeletal muscle (Z3) is measured during the measurement. Rarely moves. As shown in FIGS. 9 and 10, the first bioelectrical impedance Za changes in an increasing direction in inspiration and the first in exhalation, regardless of whether the subject's breathing is chest breathing or abdominal breathing. The bioelectrical impedance Za changes in a decreasing direction.

一方、呼吸に連動した第2生体電気インピーダンスZbの変化は、横隔膜の動きと連動している。上述したように、横隔膜は呼気時上昇し、吸気時下降するが、腹式呼吸の呼気時のみ、腹筋が緊張して内臓組織と伴に横隔膜が押し上げられ(横隔膜上昇)、内臓組織と腹部骨格筋との並列部の生体電気インピーダンスが上昇する(ΔZst>0)。このとき、肺組織の生体電気インピーダンスは減少する(ΔZlu<0)。このため、横隔膜から上部の生体電気インピーダンスの減少を横隔膜から下部の生体電気インピーダンスの増加が打ち消すように作用する。このように、胸式呼吸と腹式呼吸とでは、横隔膜から下部にある腹部骨格筋と内臓組織の動きが異なる。   On the other hand, the change of the second bioelectric impedance Zb linked to respiration is linked to the movement of the diaphragm. As described above, the diaphragm rises during exhalation and descends during inhalation, but only during abdominal breathing, the abdominal muscles are tense and the diaphragm is pushed up with the visceral tissue (diaphragm rise), and the visceral tissue and abdominal skeleton The bioelectrical impedance of the parallel part with the muscle increases (ΔZst> 0). At this time, the bioelectrical impedance of the lung tissue decreases (ΔZlu <0). For this reason, it acts so that the increase in the bioelectric impedance in the upper part from the diaphragm cancels the increase in the bioelectric impedance in the lower part from the diaphragm. Thus, the movements of the abdominal skeletal muscle and the visceral tissue below the diaphragm differ between the chest breathing and the abdominal breathing.

図9に示すように、被験者の呼吸が腹式呼吸の場合、吸気では第2生体電気インピーダンスZbは増加方向に変化する一方、呼気では、横隔膜から上部の生体電気インピーダンスの減少を横隔膜から下部の生体電気インピーダンスの増加が打ち消すように作用するので、第2生体電気インピーダンスZbは増加方向に変化する。また、図10に示すように、被験者の呼吸が胸式呼吸の場合は、上述した第1生体電気インピーダンスZaの変化と同様に、吸気では第2生体電気インピーダンスZbは増加方向に変化する一方、呼気では第2生体電気インピーダンスZbは減少方向に変化するという具合である。   As shown in FIG. 9, when the subject's breathing is abdominal breathing, the second bioelectrical impedance Zb changes in an increasing direction during inhalation, while in exhalation, a decrease in the bioelectrical impedance from the diaphragm to the lower part is performed from the diaphragm to the lower part. Since the increase in bioelectrical impedance acts to cancel out, the second bioelectrical impedance Zb changes in the increasing direction. Also, as shown in FIG. 10, when the subject's breathing is chest breathing, the second bioelectrical impedance Zb changes in an increasing direction during inspiration, as in the change in the first bioelectrical impedance Za described above, In exhalation, the second bioelectrical impedance Zb changes in a decreasing direction.

次に、図11および図12を参照して、被験者の呼吸と、被験者の胸部の周囲径Ribおよび腹部の周囲径Abとの関係について説明する。まず、被験者の呼吸が腹式呼吸である場合を想定する。図11は、レスピトレース(米国A.M.I社製)で、被験者の腹式呼吸に連動した胸部および腹部の各々の周囲径変化を時系列的に捕捉した結果を示す図である。図11からも理解されるように、被験者の呼吸が腹式呼吸の場合は、その呼吸に応じて腹部の周囲径Abが変化する一方、胸部の周囲径Ribは殆ど変化しない。したがって、腹式呼吸の場合は、被験者の胸部の周囲径Ribの変化(Ribの測定値の基準レベルを示すRib基準値に対するRib(測定値)の相対値)ΔRibと、被験者の腹部の周囲径Abの変化(Abの測定値の基準レベルを示すAb基準値に対するAb(測定値)の相対値)ΔAbとの比を示すΔRib/ΔAbは、「1」を下回るという具合である。なお、レスピトレースの情報は、測定値の基準値に対する相対値のピーク値またはボトム値の絶対値(0-P)、および、ピーク値とボトム値との絶対値の和(P-P)のうちの何れかで検出される。ここでは、被験者の1呼吸ごとに、呼吸の種別の判定が行われるので、レスピトレースの情報は、P-Pの形で検出される。   Next, with reference to FIG. 11 and FIG. 12, the relationship between the breathing of the subject and the peripheral diameter Rib of the subject's chest and the peripheral diameter Ab of the abdomen will be described. First, it is assumed that the subject's breathing is abdominal breathing. FIG. 11 is a diagram showing a result of capturing time-series changes in the peripheral diameter of each of the chest and abdomen in conjunction with the subject's abdominal breathing with Respitrace (manufactured by A.M.I., USA). As can be understood from FIG. 11, when the subject's breathing is abdominal breathing, the peripheral diameter Ab of the abdomen changes according to the breathing, while the peripheral diameter Rib of the chest hardly changes. Therefore, in the case of abdominal breathing, the change in the circumference of the subject's chest circumference Rib (the relative value of Rib (measurement value) relative to the Rib reference value indicating the reference level of the Rib measurement value) ΔRib and the circumference of the subject's abdomen ΔRib / ΔAb indicating the ratio of Ab change (relative value of Ab (measured value) to Ab reference value indicating the reference level of the measured value of Ab) ΔAb is less than “1”. The response trace information includes the peak value or absolute value of the bottom value relative to the reference value of the measurement value (0-P), and the sum of the absolute value of the peak value and the bottom value (PP). It is detected by either. Here, since the determination of the type of breathing is performed for each breath of the subject, the information on the response trace is detected in the form of PP.

次に、被験者の呼吸が胸式呼吸である場合を想定する。図12は、レスピトレース(米国A.M.I社製)で、被験者の胸式呼吸に連動した胸部および腹部の各々の周囲径変化を時系列的に捕捉した結果を示す図である。被験者の呼吸が胸式呼吸の場合は、その呼吸に応じた胸部の周囲径Ribの変化は、腹部の周囲径Abの変化よりも大きいので、上述のΔRib/ΔAbは、「1」を上回るという具合である。ここでは、ΔRib/ΔAbは、被験者の呼吸が腹式呼吸なのか胸式呼吸なのかを判別可能な判別情報であると捉えることができる。   Next, it is assumed that the subject's breathing is chest breathing. FIG. 12 is a diagram showing a result of capturing time-series changes in the peripheral diameter of each of the chest and abdomen in conjunction with the subject's chest breathing with Respitrace (manufactured by A.M.I., USA). When the subject's breathing is thoracic breathing, the change in the peripheral diameter Rib of the chest corresponding to the respiration is larger than the change in the peripheral diameter Ab of the abdomen, and thus the above-described ΔRib / ΔAb exceeds “1”. Condition. Here, ΔRib / ΔAb can be regarded as discrimination information that can discriminate whether the subject's breathing is abdominal breathing or chest breathing.

本実施形態では、被験者の胸部の周囲径の変化ΔRibと腹部の周囲径の変化ΔAbとの比(=ΔRib/ΔAb)と、第1生体電気インピーダンスZaの測定値の第1センタリング値Za0に対する相対値である第1相対値ΔZa、および、第2生体電気インピーダンスZbの測定値の第2センタリング値Zb0に対する第2相対値ΔZbとの間には相関関係があることを見出し、その相関関係を表す式を用いて、第1相対値ΔZaおよび第2相対値ΔZbに対応するΔRib/ΔAbの値を求める。そして、その求めたΔRib/ΔAbの値から、被験者の呼吸の種別(胸式呼吸なのか腹式呼吸なのか)を推定できるという具合である。この詳細な内容は後述するが、「第1センタリング値Za0」とは、第1生体電気インピーダンスZaの変化を示す波形から、呼吸に起因した情報を抽出するための振幅基準レベルを示すものである。また、「第2センタリング値Zb0」とは、第2生体電気インピーダンスZbの変化を示す波形から、呼吸に起因した情報を抽出するための振幅基準レベルを示すものである。   In the present embodiment, the ratio (= ΔRib / ΔAb) of the change ΔRib in the circumference of the chest of the subject and the change ΔAb in the circumference of the abdomen, and the relative value of the measured value of the first bioelectrical impedance Za to the first centering value Za0. It is found that there is a correlation between the first relative value ΔZa that is a value and the second relative value ΔZb of the measured value of the second bioelectrical impedance Zb with respect to the second centering value Zb0, and represents the correlation Using the equation, the value of ΔRib / ΔAb corresponding to the first relative value ΔZa and the second relative value ΔZb is obtained. The type of breathing (whether chest breathing or abdominal breathing) of the subject can be estimated from the obtained ΔRib / ΔAb value. Although the detailed contents will be described later, the “first centering value Za0” indicates an amplitude reference level for extracting information resulting from respiration from a waveform indicating a change in the first bioelectric impedance Za. . The “second centering value Zb0” indicates an amplitude reference level for extracting information resulting from respiration from a waveform indicating a change in the second bioelectric impedance Zb.

図13は、複数の被験者の測定データから得られた、ΔRib/ΔAbと、ΔZb/ΔZaとの関係を示す相関図である。図13からも理解されるように、ΔRib/ΔAbと、ΔZb/ΔZaとの間には相関係数R=0.651,P<0.01という高い相関が得られ、以下の回帰式(1)が成立する。
正の相関がみられ、以下の回帰式(1)が成立する。
ΔRib/ΔAb=a0×ΔZb/ΔZa+b0 ・・・(1)
a0:回帰係数,b0:定数。
FIG. 13 is a correlation diagram showing the relationship between ΔRib / ΔAb and ΔZb / ΔZa obtained from the measurement data of a plurality of subjects. As understood from FIG. 13, a high correlation of correlation coefficient R = 0.651 and P <0.01 is obtained between ΔRib / ΔAb and ΔZb / ΔZa, and the following regression equation (1 ) Holds.
A positive correlation is observed, and the following regression equation (1) is established.
ΔRib / ΔAb = a0 × ΔZb / ΔZa + b0 (1)
a0: regression coefficient, b0: constant.

また、上記回帰式(1)は以下のように変形できる。
ΔRib/ΔAb=(a0×ΔZb−ΔZa)/ΔZa+b1 ・・・(2)
b1:定数(=b0+1)。
The regression equation (1) can be modified as follows.
ΔRib / ΔAb = (a0 × ΔZb−ΔZa) / ΔZa + b1 (2)
b1: Constant (= b0 + 1).

ここで、呼吸に伴う第1生体電気インピーダンスZaの変化は、肺の上葉部の生体電気インピーダンス(Z1およびZ2の並列インピーダンス)の変化であると捉えることができる。一方、第2生体電気インピーダンスZbの変化は、肺の中下葉部の生体電気インピーダンス(Z4およびZ5の並列インピーダンス)の変化と、腹部の生体電気インピーダンス(Z6およびZ7の並列インピーダンス)の変化との和であると捉えることができる。肺の上葉部の生体電気インピーダンスの変化、および、肺の中下葉部の生体電気インピーダンスの変化は、同じ部位(胸部)の生体電気インピーダンスの変化であるとみなせば、第2生体電気インピーダンスZbの変化と第1生体電気インピーダンスZaの変化との差分は、腹部の生体電気インピーダンスの変化に相当する。そうすると、上記式(2)は、腹部の生体電気インピーダンスの変化と胸部の生体電気インピーダンスの変化との比と、ΔRib/ΔAbとの関係を表す式であると捉えることもできる。上記式(2)のa0は、肺の上葉部と中下葉部との測定感度の相違を補正するための補正係数であるとみなすことができる。   Here, the change in the first bioelectric impedance Za accompanying respiration can be regarded as a change in the bioelectric impedance (parallel impedance of Z1 and Z2) of the upper lobe of the lung. On the other hand, changes in the second bioelectrical impedance Zb include changes in the bioelectrical impedance (parallel impedance of Z4 and Z5) of the middle and lower lobes of the lung, and changes in the bioelectrical impedance of the abdomen (parallel impedance of Z6 and Z7). Can be seen as the sum of If the change of the bioelectric impedance of the upper lobe of the lung and the change of the bioelectric impedance of the middle and lower lobe of the lung are considered to be changes of the bioelectric impedance of the same part (chest), the second bioelectric impedance The difference between the change in Zb and the change in the first bioelectric impedance Za corresponds to the change in the bioelectric impedance of the abdomen. Then, the above equation (2) can be regarded as an equation representing the relationship between ΔRib / ΔAb and the ratio between the change in the bioelectrical impedance of the abdomen and the change in the bioelectrical impedance of the chest. In the above equation (2), a0 can be regarded as a correction coefficient for correcting a difference in measurement sensitivity between the upper lobe portion and the middle lower lobe portion of the lung.

<1−4:呼吸解析処理>
次に、CPU170が実行する呼吸解析処理について説明する。図14は、呼吸解析処理の具体的な内容を説明するためのフローチャートである。本実施形態では、通常の1呼吸(=1回の吸気+1回の呼気)につき、10回の呼吸解析処理を実行するように設定される。ここでは、通常の1呼吸に要する時間を4秒とみなし、CPU170は、0.4秒ごとに、呼吸解析処理を実行するという具合である。以下では、呼吸解析処理を実行するタイミング(0.4秒ごとのタイミング)をサンプリングタイミングと呼ぶ。なお、これは一例であり、呼吸解析処理を実行するタイミングは任意に設定可能である。
<1-4: Respiration analysis processing>
Next, the breath analysis process executed by the CPU 170 will be described. FIG. 14 is a flowchart for explaining specific contents of the breath analysis process. In the present embodiment, it is set so that 10 breath analysis processes are executed for one normal breath (= one inspiration + 1 exhalation). Here, the time required for one normal breath is regarded as 4 seconds, and the CPU 170 executes the breath analysis process every 0.4 seconds. Hereinafter, the timing for executing the breath analysis process (timing every 0.4 seconds) is referred to as sampling timing. This is only an example, and the timing for executing the breath analysis process can be arbitrarily set.

図14に示すように、まず、CPU170は、サンプリングタイミングに到達したか否かを判定し(ステップS10)、ステップS10の結果が肯定である場合はステップS20に進む。ここでは、第n番目(n≧1)のサンプリングタイミングに到達した場合を想定して、ステップS20以下の各ステップの具体的な内容を説明する。ステップS20以下の各ステップの具体的な説明に先立ち、まずは、各ステップの内容の概略を簡単に説明する。ステップS10の後のステップS20において、CPU170は、肺の上部を含み、腹部を含まない体幹上部の第1生体電気インピーダンスZaを測定する。ステップS20の後のステップS30において、CPU170は、肺の中下部および腹部を含む体幹中部の第2生体電気インピーダンスZbを測定する。ステップS30の後のステップS40において、CPU170は、ステップS20で測定した第1生体電気インピーダンスZaおよびステップS30で測定した第2生体電気インピーダンスZbの各々について、スムージング処理を実行する。ステップS40の後のステップS50において、CPU170は、第1生体電気インピーダンスZaの測定値の振幅基準レベルを示す第1センタリング値Za0を生成する。ステップS50の後のステップS60において、CPU170は、第1生体電気インピーダンスZaの測定値の第1センタリング値Za0に対する相対値である第1相対値ΔZaを算出する。ステップS60の後のステップS70において、CPU170は、第2生体電気インピーダンスZbの測定値の振幅基準レベルを示す第2センタリング値Zb0を生成し、その生成した第2センタリング値Zb0を用いて、第2生体電気インピーダンスZbの測定値の第2センタリング値Zb0に対する相対値である第2相対値ΔZbを算出する。ステップS70の後のステップS80において、CPU170は、上記回帰式(2)にしたがって演算処理を実行することで、第1相対値ΔZaおよび第2相対値ΔZbに対応するΔRib/ΔAbの値を求める。ステップS80の後のステップS90において、CPU170は、被験者の呼吸の深さを示す呼吸レベルを抽出する呼吸レベル抽出処理を実行する。以下、各ステップの具体的な内容を順番に説明していく。   As shown in FIG. 14, first, the CPU 170 determines whether or not the sampling timing has been reached (step S10), and proceeds to step S20 if the result of step S10 is affirmative. Here, assuming the case where the nth (n ≧ 1) sampling timing has been reached, the specific contents of each step after step S20 will be described. Prior to specific description of each step after step S20, first, an outline of the contents of each step will be briefly described. In step S20 after step S10, the CPU 170 measures the first bioelectric impedance Za of the upper trunk including the upper part of the lung and not including the abdomen. In step S30 after step S20, the CPU 170 measures the second bioelectric impedance Zb of the middle trunk including the middle lower part and the abdomen of the lung. In step S40 after step S30, the CPU 170 performs a smoothing process on each of the first bioelectric impedance Za measured in step S20 and the second bioelectric impedance Zb measured in step S30. In step S50 after step S40, the CPU 170 generates a first centering value Za0 indicating the amplitude reference level of the measured value of the first bioelectrical impedance Za. In step S60 after step S50, the CPU 170 calculates a first relative value ΔZa that is a relative value of the measured value of the first bioelectrical impedance Za with respect to the first centering value Za0. In step S70 after step S60, the CPU 170 generates a second centering value Zb0 indicating the amplitude reference level of the measurement value of the second bioelectrical impedance Zb, and uses the generated second centering value Zb0 to generate the second centering value Zb0. A second relative value ΔZb that is a relative value of the measured value of the bioelectrical impedance Zb with respect to the second centering value Zb0 is calculated. In step S80 after step S70, the CPU 170 calculates the value of ΔRib / ΔAb corresponding to the first relative value ΔZa and the second relative value ΔZb by executing arithmetic processing according to the regression equation (2). In step S90 after step S80, the CPU 170 executes a respiration level extraction process for extracting a respiration level indicating the respiration depth of the subject. Hereinafter, the specific contents of each step will be described in order.

図14に示すように、ステップS20において、CPU170は、体幹上部の第1生体電気インピーダンスZaを測定する。より具体的には、CPU170は、左手用の電流電極X3と右手用の電流電極X4とを選択するように電極切換回路252を制御するとともに、右足用の電圧電極Y2と右手用の電圧電極Y4とを選択するように電極切替回路251を制御する。CPU170は、右手と左手との間に流れる基準電流Irefの大きさを示す電流データDiと、右足用の電圧電極Y2と右手用の電圧電極Y4との間の電位差を示す電圧データDvとから、右上肢(体幹上部)の第1生体電気インピーダンスZaを測定する。ここでは、第n番目(n≧1)のサンプリングタイミングにおける第1生体電気インピーダンスの実測値をZa(n)’と表記する。   As shown in FIG. 14, in step S20, the CPU 170 measures the first bioelectric impedance Za of the upper trunk. More specifically, the CPU 170 controls the electrode switching circuit 252 so as to select the current electrode X3 for the left hand and the current electrode X4 for the right hand, and the voltage electrode Y2 for the right foot and the voltage electrode Y4 for the right hand. The electrode switching circuit 251 is controlled so as to select. From the current data Di indicating the magnitude of the reference current Iref flowing between the right hand and the left hand, and the voltage data Dv indicating the potential difference between the right foot voltage electrode Y2 and the right hand voltage electrode Y4, the CPU 170 The first bioelectric impedance Za of the upper right limb (upper trunk) is measured. Here, the measured value of the first bioelectrical impedance at the nth (n ≧ 1) sampling timing is expressed as Za (n) ′.

ステップS20の後、CPU170は、体幹中部の第2生体電気インピーダンスZbを測定する(ステップS30)。より具体的には、CPU170は、左足用の電流電極X1と右手用の電流電極X4とを選択するように電極切換回路252を制御するとともに、左手用の電圧電極Y3と右足用の電圧電極Y2とを選択するように電極切換回路251を制御する。CPU170は、左足と右手との間に流れる基準電流Irefの大きさを示す電流データDiと、左手用の電圧電極Y3と右足用の電圧電極Y2との間の電位差を示す電圧データDvとから、体幹中部の第2生体電気インピーダンスZbを測定する。ここでは、第n番目のサンプリングタイミングにおける第2生体電気インピーダンスの実測値をZb(n)’と表記する。   After step S20, the CPU 170 measures the second bioelectric impedance Zb of the middle trunk (step S30). More specifically, the CPU 170 controls the electrode switching circuit 252 so as to select the current electrode X1 for the left foot and the current electrode X4 for the right hand, and the voltage electrode Y3 for the left hand and the voltage electrode Y2 for the right foot. The electrode switching circuit 251 is controlled so as to select. From the current data Di indicating the magnitude of the reference current Iref flowing between the left foot and the right hand, and the voltage data Dv indicating the potential difference between the voltage electrode Y3 for the left hand and the voltage electrode Y2 for the right foot, The second bioelectric impedance Zb in the middle of the trunk is measured. Here, the measured value of the second bioelectrical impedance at the nth sampling timing is expressed as Zb (n) ′.

ステップS30の後、CPU170は、ステップS20で測定した第1生体電気インピーダンスZa(n)’およびステップ30で測定した第2生体電気インピーダンスZb(n)’の各々について、スムージング処理を実行する(ステップS40)。まず、第1生体電気インピーダンスZa(n)’のスムージング処理について具体的に説明する。CPU170は、第n−2番目のサンプリングタイミングにおける第1生体電気インピーダンスの実測値Za(n−2)’と、第n−1番目のサンプリングタイミングにおける第1生体電気インピーダンスの実測値Za(n−1)’と、第n番目のサンプリングタイミングにおける第1生体電気インピーダンスの実測値Za(n)’とを用いた移動平均処理を行う。そして、その処理結果を、第n番目のサンプリングタイミングにおける第1生体電気インピーダンスの測定値として採用する(スムージング処理)。ここでは、スムージング処理が行われた後の、第n番目のサンプリングタイミングにおける第1生体電気インピーダンスの測定値を、Za(n)と表記する。   After step S30, the CPU 170 executes a smoothing process for each of the first bioelectrical impedance Za (n) ′ measured in step S20 and the second bioelectrical impedance Zb (n) ′ measured in step 30 (step S30). S40). First, the smoothing process of the first bioelectric impedance Za (n) ′ will be specifically described. The CPU 170 measures the measured value Za (n-2) ′ of the first bioelectrical impedance at the (n−2) th sampling timing and the measured value Za (n−) of the first bioelectrical impedance at the (n−1) th sampling timing. 1) ′ and a moving average process using the measured value Za (n) ′ of the first bioelectrical impedance at the nth sampling timing. And the process result is employ | adopted as a measured value of the 1st bioelectrical impedance in the nth sampling timing (smoothing process). Here, the measured value of the first bioelectrical impedance at the nth sampling timing after the smoothing process is performed is denoted as Za (n).

次に、第2生体電気インピーダンスZb(n)’のスムージング処理について具体的に説明する。CPU170は、第n−2番目のサンプリングタイミングにおける第2生体電気インピーダンスの実測値Zb(n−2)’と、第n−1番目のサンプリングタイミングにおける第2生体電気インピーダンスの実測値Zb(n−1)’と、第n番目のサンプリングタイミングにおける第2生体電気インピーダンスの実測値Zb(n)’とを用いた移動平均処理を行う。そして、その処理結果を、第n番目のサンプリングタイミングにおける第2生体電気インピーダンスの測定値として採用する(スムージング処理)。ここでは、スムージング処理が行われた後の、第n番目のサンプリングタイミングにおける第2生体電気インピーダンスの測定値を、Zb(n)と表記する。   Next, the smoothing process of the second bioelectrical impedance Zb (n) ′ will be specifically described. The CPU 170 measures the second measured value Zb (n-2) ′ of the second bioelectrical impedance at the (n−2) th sampling timing and the measured value Zb (n−) of the second bioelectrical impedance at the (n−1) th sampling timing. 1) ′ and a moving average process using the measured value Zb (n) ′ of the second bioelectrical impedance at the n-th sampling timing are performed. And the process result is employ | adopted as a measured value of the 2nd bioelectrical impedance in the nth sampling timing (smoothing process). Here, the measured value of the second bioelectrical impedance at the nth sampling timing after the smoothing process is performed is denoted as Zb (n).

ステップS40の後、CPU170は、第1生体電気インピーダンスZaの測定値の振幅基準レベルを示す第1センタリング値Za0を生成する第1センタリング処理を実行する(ステップS50)。ここでは、第n番目のサンプリングタイミングにおける第1センタリング値をZa0(n)と表記する。本実施形態では、CPU170は、第n番目のサンプリングタイミングよりも所定の時間長だけ前の時点を始点とし、第n番目のサンプリングタイミングを終点とするセンタリング期間内の複数のサンプリングタイミングの各々における第1生体電気インピーダンスZaの測定値に基づいて、第n番目のサンプリングタイミングにおける第1センタリング値Za0(n)を生成する。センタリング期間の時間長は、第n番目のサンプリングタイミングにおける被験者の呼吸の速度に応じて可変に設定される。以下、その具体的な内容について詳細に説明する。   After step S40, the CPU 170 executes a first centering process for generating a first centering value Za0 indicating the amplitude reference level of the measured value of the first bioelectrical impedance Za (step S50). Here, the first centering value at the nth sampling timing is expressed as Za0 (n). In the present embodiment, the CPU 170 starts the time point a predetermined time before the n-th sampling timing as the start point and the n-th sampling timing as the end point in each of the plurality of sampling timings in the centering period. Based on the measured value of the 1 bioelectrical impedance Za, the first centering value Za0 (n) at the nth sampling timing is generated. The time length of the centering period is variably set according to the breathing speed of the subject at the nth sampling timing. The specific contents will be described in detail below.

図15は、第1センタリング処理の具体的な内容を示すフローチャートである。図15に示すように、まず、CPU170は、10個のサンプリングタイミングの各々における第1生体電気インピーダンスZaの測定値の振幅基準レベルを示すMA10を抽出するMA10抽出処理を実行する(ステップS51)。より具体的には、CPU170は、第n−9番目〜第n番目の10個のサンプリングタイミングの各々における第1生体電気インピーダンスZaの測定値(Za(n−9)〜Za(n))を用いて移動平均処理を行い、その処理結果を、第n番目のサンプリングタイミングにおけるMA10(n)として抽出する([Za(n−9)+Za(n−8)+・・・+Za(n)]/10→MA10(n))。   FIG. 15 is a flowchart showing specific contents of the first centering process. As shown in FIG. 15, first, the CPU 170 executes an MA10 extraction process for extracting MA10 indicating the amplitude reference level of the measurement value of the first bioelectrical impedance Za at each of the ten sampling timings (step S51). More specifically, the CPU 170 obtains the measured values (Za (n-9) to Za (n)) of the first bioelectric impedance Za at each of the n-9th to nth sampling timings. The moving average processing is performed, and the processing result is extracted as MA10 (n) at the nth sampling timing ([Za (n−9) + Za (n−8) +... + Za (n)] / 10 → MA10 (n)).

ステップS51の後、CPU170は、20個のサンプリングタイミングの各々における第1生体電気インピーダンスZaの測定値の振幅基準レベルを示すMA20を抽出するMA20抽出処理を実行する(ステップS52)。より具体的には、CPU170は、第n−19番目〜第n番目の20個のサンプリングタイミングの各々における第1生体電気インピーダンスZaの測定値(Za(n−19)〜Za(n))を用いて移動平均処理を行い、その処理結果を、第n番目のサンプリングタイミングにおけるMA20(n)として抽出する([Za(n−19)+Za(n−18)+・・・+Za(n)]/20→MA20(n))。   After step S51, the CPU 170 executes MA20 extraction processing for extracting MA20 indicating the amplitude reference level of the measurement value of the first bioelectrical impedance Za at each of the 20 sampling timings (step S52). More specifically, the CPU 170 obtains measured values (Za (n-19) to Za (n)) of the first bioelectric impedance Za at each of the n-19th to nth 20 sampling timings. The moving average processing is performed, and the processing result is extracted as MA20 (n) at the nth sampling timing ([Za (n−19) + Za (n−18) +... + Za (n)] / 20 → MA20 (n)).

ステップS52の後、CPU170は、10個のサンプリングタイミングの各々における第1生体電気インピーダンスZaの測定値のうち最大の値をMAX10として抽出するMAX10抽出処理を実行する(ステップS53)。より具体的には、CPU170は、第n−9番目〜第n番目の10個のサンプリングタイミングの各々における第1生体電気インピーダンスZaの測定値のうち最大の値を、第n番目のサンプリングタイミングにおけるMAX10(n)として抽出するという具合である。   After step S52, the CPU 170 executes MAX10 extraction processing for extracting the maximum value as MAX10 among the measured values of the first bioelectrical impedance Za at each of the ten sampling timings (step S53). More specifically, the CPU 170 calculates the maximum value among the measured values of the first bioelectrical impedance Za at each of the n-9th to nth sampling timings at the nth sampling timing. That is, it is extracted as MAX10 (n).

ステップS53の後、CPU170は、10個のサンプリングタイミングの各々における第1生体電気インピーダンスZaの測定値のうち最小の値をMIN10として抽出するMIN10抽出処理を実行する(ステップS54)。より具体的には、CPU170は、第n−9番目〜第n番目の10個のサンプリングタイミングの各々における第1生体電気インピーダンスZaの測定値のうち最小の値を、第n番目のサンプリングタイミングにおけるMIN10(n)として抽出するという具合である。   After step S53, the CPU 170 executes a MIN10 extraction process for extracting the minimum value of the measured values of the first bioelectric impedance Za at each of the ten sampling timings as MIN10 (step S54). More specifically, the CPU 170 obtains the minimum value among the measured values of the first bioelectrical impedance Za at each of the n-9th to nth sampling timings at the nth sampling timing. That is, it is extracted as MIN10 (n).

ステップS54の後、CPU170は、20個のサンプリングタイミングの各々におけるMAX10とMIN10との平均値(第n番目のサンプリングタイミングにおける平均値をAV10(n)と表記)について移動平均処理を行い、その処理結果を、中央値として算出する中央値算出処理を実行する(ステップS55)。より具体的には、CPU170は、第n−19番目〜第n番目の20個のサンプリングタイミングの各々における平均値(AV10(n−19)〜AV10(n))について移動平均処理を行い、その処理結果を、第n番目のサンプリングタイミングにおける中央値CNT20(n)として抽出する([AV10(n−19)+AV10(n−18)+・・・+AV10(n)]/20→CNT20(n))。ここでは、説明を省略するが、中央値CNT20(n)は、体動などに起因するアーチファクト(データ波形の歪み)等による処理に適さない異常波形の抽出に用いられる。   After step S54, the CPU 170 performs a moving average process on the average value of MAX10 and MIN10 at each of the 20 sampling timings (the average value at the nth sampling timing is expressed as AV10 (n)). A median calculation process for calculating the result as a median is executed (step S55). More specifically, the CPU 170 performs a moving average process on average values (AV10 (n-19) to AV10 (n)) at each of the n-19th to nth 20 sampling timings, The processing result is extracted as the median value CNT20 (n) at the nth sampling timing ([AV10 (n−19) + AV10 (n−18) +... + AV10 (n)] / 20 → CNT20 (n) ). Although the description is omitted here, the median value CNT20 (n) is used to extract an abnormal waveform that is not suitable for processing due to artifacts (distortion of data waveform) caused by body movements or the like.

ステップS55の後、CPU170は、第n番目のサンプリングタイミングにおける被験者の呼吸のタイミングを抽出する呼吸タイミング抽出処理を実行する(ステップS56)。以下では、図16および図17を参照しながら、呼吸タイミング抽出処理の具体的な内容を説明する。図16および図17は、呼吸タイミング抽出処理の具体的な内容を示すフローチャートである。図16に示すように、まずCPU170は、第n番目のサンプリングタイミングにおける第1生体電気インピーダンスZaの微分係数dZa(n)を抽出する微分係数抽出処理を実行する(ステップS201)。より具体的には、CPU170は、以下の式(3)にしたがって演算処理を実行することで、微分係数dZa(n)を抽出する。
[Za(n)−Za(n−2)]/0.8=dZa(n) ・・・(3)
After step S55, the CPU 170 executes a breathing timing extraction process for extracting the breathing timing of the subject at the nth sampling timing (step S56). Hereinafter, specific contents of the breathing timing extraction process will be described with reference to FIGS. 16 and 17. 16 and 17 are flowcharts showing specific contents of the breathing timing extraction process. As shown in FIG. 16, first, the CPU 170 executes a differential coefficient extraction process for extracting the differential coefficient dZa (n) of the first bioelectric impedance Za at the nth sampling timing (step S201). More specifically, the CPU 170 extracts the differential coefficient dZa (n) by executing arithmetic processing according to the following equation (3).
[Za (n) −Za (n−2)] / 0.8 = dZa (n) (3)

次に、CPU170は、ステップS201で抽出した微分係数dZa(n)の絶対値が0.1より小さいか否かを判定する(ステップS202)。ステップS202の結果が肯定である場合、CPU170は、微分係数dZa(n)の極性判別フラッグF0(n)を「0」に設定してステップS204に進む。極性判別フラッグF0(n)が「0」であるとは、第n番目のサンプリングタイミングにおいて、第1生体電気インピーダンスZaの値は極大値(ピーク値)または極小値(ボトム値)であることを意味する。   Next, the CPU 170 determines whether or not the absolute value of the differential coefficient dZa (n) extracted in step S201 is smaller than 0.1 (step S202). If the result of step S202 is affirmative, the CPU 170 sets the polarity determination flag F0 (n) of the differential coefficient dZa (n) to “0” and proceeds to step S204. The polarity determination flag F0 (n) being “0” means that the value of the first bioelectrical impedance Za is a maximum value (peak value) or a minimum value (bottom value) at the n-th sampling timing. means.

一方、ステップS202の結果が否定である場合、CPU170は、微分係数dZa(n)の値が0より大きいか否かを判定する(ステップS203)。ステップS203の結果が肯定である場合、CPU170は、極性判別フラッグF0(n)を「+1」に設定してステップS204へ進む。極性判別フラッグF0(n)が「+1」であるとは、第n番目のサンプリングタイミングにおいて、第1生体電気インピーダンスZaの変化の方向は正側であることを意味する。ステップS203の結果が否定である場合、CPU170は、極性判別フラッグF0(n)を「-1」に設定してステップS204へ進む。極性判別フラッグF0(n)が「-1」であるとは、第n番目のサンプリングタイミングにおいて、第1生体電気インピーダンスZaの変化の方向は負側であることを意味する。   On the other hand, when the result of step S202 is negative, the CPU 170 determines whether or not the value of the differential coefficient dZa (n) is greater than 0 (step S203). If the result of step S203 is affirmative, the CPU 170 sets the polarity determination flag F0 (n) to “+1” and proceeds to step S204. The polarity determination flag F0 (n) being “+1” means that the change direction of the first bioelectrical impedance Za is on the positive side at the nth sampling timing. If the result of step S203 is negative, the CPU 170 sets the polarity determination flag F0 (n) to “−1” and proceeds to step S204. The polarity determination flag F0 (n) being “−1” means that the direction of change of the first bioelectrical impedance Za is negative on the nth sampling timing.

ステップS204において、CPU170は、第n番目のサンプリングタイミングにおける極性判別フラッグF0(n)の絶対値と、第n−1番目のサンプリングタイミングにおける極性判別フラッグF0(n−1)の絶対値とが等しく、かつ、F0(n−1)の値とF0(n)の値とが等しくないか否かを判定する。ステップS204の結果が肯定である場合、CPU170は、F0(n)を「0」に設定して、次のステップS206(図17参照)へ進む。ステップS204の結果が否定である場合、CPU170は、ステップS204の直前で設定したF0(n)の値を維持したまま、次のステップS206へ進む。   In step S204, the CPU 170 makes the absolute value of the polarity determination flag F0 (n) at the nth sampling timing equal to the absolute value of the polarity determination flag F0 (n-1) at the (n-1) th sampling timing. In addition, it is determined whether the value of F0 (n-1) is not equal to the value of F0 (n). If the result of step S204 is affirmative, the CPU 170 sets F0 (n) to “0” and proceeds to the next step S206 (see FIG. 17). If the result of step S204 is negative, the CPU 170 proceeds to the next step S206 while maintaining the value of F0 (n) set immediately before step S204.

図17を参照しながら、呼吸タイミング抽出処理の具体的な内容の説明を続ける。CPU170は、第n番目のサンプリングタイミングにおける極性判別フラッグF0(n)が「0」であるか否かを判定する(ステップS206)。ステップS206の結果が否定である場合、CPU170は、ピーク/ボトム判別フラッグF1(n)を「0」に設定して、ステップS209へ進む。ピーク/ボトム判別フラッグF1(n)が「0」であるとは、第n番目のサンプリングタイミングにおける第1生体電位インピーダンスの測定値Za(n)は、ピーク値またはボトム値ではないことを意味する。
一方、ステップS206の結果が肯定である場合、CPU170は、3個のサンプリングタイミングの各々における極性判別フラッグF0の和が「+1」よりも大きいか否かを判定する(ステップS207)。より具体的には、CPU170は、第n−2番目のサンプリングタイミング〜第n番目のサンプリングタイミングの各々における極性判別フラッグ(F0(n−2)〜F0(n))の和が「+1」よりも大きいか否かを判定する。
The description of the specific content of the breathing timing extraction process will be continued with reference to FIG. The CPU 170 determines whether or not the polarity determination flag F0 (n) at the nth sampling timing is “0” (step S206). If the result of step S206 is negative, the CPU 170 sets the peak / bottom determination flag F1 (n) to “0” and proceeds to step S209. The peak / bottom determination flag F1 (n) being “0” means that the measured value Za (n) of the first biopotential impedance at the nth sampling timing is not a peak value or a bottom value. .
On the other hand, when the result of step S206 is affirmative, the CPU 170 determines whether or not the sum of the polarity determination flags F0 at each of the three sampling timings is greater than “+1” (step S207). More specifically, the CPU 170 determines that the sum of the polarity determination flags (F0 (n-2) to F0 (n)) at each of the (n-2) th sampling timing to the nth sampling timing is "+1". It is judged whether it is larger than.

ステップS207の結果が肯定である場合、CPU170は、ピーク/ボトム判別フラッグF1(n)を「+1」に設定して、ステップS209へ進む。ピーク/ボトム判別フラッグF1(n)が「+1」であるとは、第n番目のサンプリングタイミングにおける第1生体電気インピーダンスの測定値Za(n)はピーク値(最大値)であることを意味する。   If the result of step S207 is affirmative, the CPU 170 sets the peak / bottom determination flag F1 (n) to “+1”, and proceeds to step S209. The peak / bottom discrimination flag F1 (n) being “+1” means that the measured value Za (n) of the first bioelectrical impedance at the nth sampling timing is a peak value (maximum value). To do.

ステップS207の結果が否定である場合、CPU170は、3個のサンプリングタイミングの各々における極性判別フラッグF0の和が「-1」よりも小さいか否かを判定する(ステップS208)。より具体的には、CPU170は、第n−2番目のサンプリングタイミング〜第n番目のサンプリングタイミングの各々における極性判別フラッグ(F0(n−2)〜F0(n))の和が「-1」よりも小さいか否かを判定する。   If the result of step S207 is negative, the CPU 170 determines whether or not the sum of the polarity determination flags F0 at each of the three sampling timings is smaller than “−1” (step S208). More specifically, the CPU 170 determines that the sum of the polarity determination flags (F0 (n-2) to F0 (n)) at each of the (n-2) th sampling timing to the nth sampling timing is "-1". Or less.

ステップS208の結果が肯定である場合、CPU170は、ピーク/ボトム判別フラッグF1(n)を「-1」に設定して、ステップS209へ進む。ピーク/ボトム判別フラッグF1(n)が「-1」であるとは、第n番目のサンプリングタイミングにおける第1生体電気インピーダンスの測定値Za(n)はボトム値(最小値)であることを意味する。一方、ステップS208の結果が否定である場合、CPU170は、ピーク/ボトム判別フラッグF1(n)を「0」に設定して、ステップS209へ進むという具合である。   If the result of step S208 is affirmative, the CPU 170 sets the peak / bottom determination flag F1 (n) to “−1”, and proceeds to step S209. The peak / bottom discrimination flag F1 (n) being “−1” means that the measured value Za (n) of the first bioelectrical impedance at the nth sampling timing is a bottom value (minimum value). To do. On the other hand, if the result of step S208 is negative, the CPU 170 sets the peak / bottom discrimination flag F1 (n) to “0” and proceeds to step S209.

ステップS209において、CPU170は、ピーク/ボトム判別フラッグF1(n)が「+1」であるか否かを判定する。ステップS209の結果が否定である場合、CPU170は、直前の第n−1番目のサンプリングタイミングにおけるサンプリングカウンタ値N(n−1)に1を加算する一方(ステップS210)、ステップS209の結果が肯定である場合、CPU170は、サンプリングカウンタ値Nを初期化する(ステップS211)。ここで、図9および図10からも理解されるように、被験者の呼吸が胸式呼吸および腹式呼吸の何れの場合であっても、呼吸に伴う第1生体電気インピーダンスZaの変化を示す波形は略正弦波状であるところ、サンプリングカウンタ値Nは、第1生体電気インピーダンスZaの測定値がピーク値に到達するたびに初期化(サンプリングカウンタ値=0)され、次のピーク値に到達するまでのサンプリングタイミングの回数が順次にカウントされていくという具合である。以上で、図15のステップS56における呼吸タイミング抽出処理が終了する。   In step S209, the CPU 170 determines whether or not the peak / bottom determination flag F1 (n) is “+1”. When the result of step S209 is negative, the CPU 170 adds 1 to the sampling counter value N (n−1) at the immediately preceding n−1th sampling timing (step S210), while the result of step S209 is positive. If so, the CPU 170 initializes the sampling counter value N (step S211). Here, as can be understood from FIGS. 9 and 10, a waveform showing a change in the first bioelectrical impedance Za accompanying respiration regardless of whether the subject's respiration is chest respiration or abdominal respiration. Is substantially sinusoidal, the sampling counter value N is initialized (sampling counter value = 0) every time the measured value of the first bioelectrical impedance Za reaches the peak value, and until the next peak value is reached. The number of sampling timings is sequentially counted. This completes the breathing timing extraction process in step S56 of FIG.

再び図15に戻って説明を続ける。上述の呼吸タイミング抽出処理が終了すると、CPU170は、被験者の呼吸が速めの呼吸なのか遅めの呼吸なのかを判別する呼吸スピード判別フラッグを設定する(ステップS57)。以下、図18を参照しながら、ステップS57でCPU170が実行する呼吸スピード判別フラッグ設定処理の具体的な内容を説明する。図18は、呼吸スピード判別フラッグ設定処理の具体的な内容を示すフローチャートである。図18に示すように、まずCPU170は、極性判別フラッグF0(n)が「0」であるか否かを判定する(ステップS301)。ステップS301の結果が否定である場合、CPU170は、第n番目のサンプリングタイミングにおける呼吸スピード判別フラッグFma(n)は、直前の第n−1番目のサンプリングタイミングにおける呼吸スピード判別フラッグFma(n−1)に等しいとみなして処理を終了する。   Returning to FIG. 15 again, the description will be continued. When the above-described respiration timing extraction process ends, the CPU 170 sets a respiration speed determination flag for determining whether the subject's respiration is fast or slow respiration (step S57). Hereinafter, specific contents of the breathing speed determination flag setting process executed by the CPU 170 in step S57 will be described with reference to FIG. FIG. 18 is a flowchart showing specific contents of the breathing speed determination flag setting process. As shown in FIG. 18, first, the CPU 170 determines whether or not the polarity determination flag F0 (n) is “0” (step S301). When the result of step S301 is negative, the CPU 170 determines that the breathing speed determination flag Fma (n) at the nth sampling timing is the breathing speed determination flag Fma (n−1) at the immediately preceding n−1th sampling timing. ) And the process ends.

ステップS301の結果が肯定である場合、CPU170は、ピーク/ボトム判別フラッグF1(n)が「+1」であるか否かを判定する(ステップS302)。ステップS302の結果が肯定である場合、CPU170は、直前の第n−1番目のサンプリングタイミングにおけるサンプリングカウンタ値N(n−1)が「10」よりも大きいか否かを判定する(ステップS303)。ここで、被験者の呼吸のスピードが遅ければ、第1生体電気インピーダンスZaの測定値がピーク値に到達してから、次のピーク値に到達するまでの時間長は長くなり、次のピーク値に到達する直前のサンプリングカウンタ値Nも大きくなる。本実施形態では、CPU170は、第n番目のサンプリングタイミングにおいて第1生体電気インピーダンスZaの測定値がピーク値に到達したと判断した場合は、その直前の第n−1番目のサンプリングタイミングにおけるサンプリングカウンタ値が「10」よりも大きいか否かを判定し、当該サンプリングカウンタ値が「10」よりも大きいと判定した場合は、被験者の呼吸は遅めの呼吸であると判断する。具体的には、ステップS303の結果が肯定である場合、CPU170は、呼吸スピード判別フラッグFma(n)を「20」に設定して処理を終了する。呼吸スピード判別フラッグFma(n)が「20」であるとは、被験者の呼吸が遅めの呼吸であることを意味する。また、ステップS303の結果が否定である場合、CPU170は、被験者の呼吸は速めの呼吸であると判断して呼吸スピード判別フラッグFma(n)を「10」に設定して処理を終了する。呼吸スピード判別フラッグFma(n)が「10」であるとは、被験者の呼吸が速めの呼吸であることを意味する。   If the result of step S301 is affirmative, the CPU 170 determines whether or not the peak / bottom determination flag F1 (n) is “+1” (step S302). If the result of step S302 is affirmative, the CPU 170 determines whether or not the sampling counter value N (n−1) at the immediately preceding n−1th sampling timing is greater than “10” (step S303). . Here, if the breathing speed of the subject is slow, the time length from when the measured value of the first bioelectrical impedance Za reaches the peak value until it reaches the next peak value becomes longer, and the next peak value is reached. The sampling counter value N immediately before reaching also increases. In the present embodiment, when the CPU 170 determines that the measured value of the first bioelectrical impedance Za has reached the peak value at the nth sampling timing, the sampling counter at the n−1th sampling timing immediately before the measured value. It is determined whether or not the value is greater than “10”, and if it is determined that the sampling counter value is greater than “10”, it is determined that the subject's breathing is slow breathing. Specifically, if the result of step S303 is affirmative, the CPU 170 sets the breathing speed determination flag Fma (n) to “20” and ends the process. The breathing speed determination flag Fma (n) being “20” means that the subject's breathing is slow breathing. If the result of step S303 is negative, the CPU 170 determines that the subject's breathing is fast breathing, sets the breathing speed discrimination flag Fma (n) to “10”, and ends the process. The breathing speed determination flag Fma (n) being “10” means that the subject's breathing is fast breathing.

一方、ステップS302の結果が否定である場合、CPU170は、ピーク/ボトム判別フラッグF1(n)が「-1」であるか否かを判定する(ステップS304)。ステップS304の結果が否定である場合、CPU170は、第n番目のサンプリングタイミングにおける呼吸スピード判別フラッグFma(n)は、直前の第n−1番目のサンプリングタイミングにおける呼吸スピード判別フラッグFma(n−1)に等しいとみなして処理を終了する。ステップS304の結果が肯定である場合、CPU170は、直前の第n−1番目のサンプリングタイミングにおけるサンプリングカウンタ値N(n−1)が「5」よりも大きいか否かを判定する(ステップS305)。本実施形態では、CPU170は、ボトム値に到達する直前のサンプリングカウンタ値N(n−1)が「5」よりも大きい場合は、被験者の呼吸は遅めの呼吸であると判断する一方、ボトム値に到達する直前のサンプリングカウンタ値N(n−1)が「5」よりも小さい場合は、被験者の呼吸は速めの呼吸であると判断する。具体的には、CPU170は、ステップS305の結果が肯定である場合は、呼吸スピード判別フラッグFma(n)を「20」に設定して処理を終了する一方、ステップS305の結果が否定である場合は、呼吸スピード判別フラッグFma(n)を「10」に設定して処理を終了するという具合である。以上で、図15のステップS57における呼吸スピード判別フラッグ設定処理が終了する。   On the other hand, when the result of step S302 is negative, the CPU 170 determines whether or not the peak / bottom determination flag F1 (n) is “−1” (step S304). When the result of step S304 is negative, the CPU 170 determines that the respiration speed determination flag Fma (n) at the nth sampling timing is the respiration speed determination flag Fma (n−1) at the immediately preceding n−1th sampling timing. ) And the process ends. If the result of step S304 is affirmative, the CPU 170 determines whether or not the sampling counter value N (n−1) at the immediately preceding (n−1) th sampling timing is greater than “5” (step S305). . In the present embodiment, when the sampling counter value N (n−1) immediately before reaching the bottom value is larger than “5”, the CPU 170 determines that the subject's breathing is slow breathing, while the bottom When the sampling counter value N (n−1) immediately before reaching the value is smaller than “5”, it is determined that the subject's breathing is fast breathing. Specifically, if the result of step S305 is affirmative, CPU 170 sets the respiration speed determination flag Fma (n) to “20” and ends the process, while the result of step S305 is negative. Indicates that the breathing speed determination flag Fma (n) is set to “10” and the process is terminated. This completes the breathing speed determination flag setting process in step S57 of FIG.

再び図15に戻って説明を続ける。上述の呼吸スピード判別フラッグ設定処理が終了すると、CPU170は、第n番目のサンプリングタイミングにおける第1センタリング値Za0(n)を抽出する(ステップS58)。以下、図19を参照しながら、ステップS58でCPU170が実行する第1センタリング値抽出処理の具体的な内容を説明する。図19は、第1センタリング値抽出処理の具体的な内容を示すフローチャートである。図19に示すように、まずCPU170は、呼吸スピード判別フラッグFma(n)が「10」であるか否かを判定する(ステップS401)。言い換えれば、CPU170は、第n番目のサンプリングタイミングにおける被験者の呼吸が速めの呼吸であるか否かを判定するという具合である。   Returning to FIG. 15 again, the description will be continued. When the above-described breathing speed determination flag setting process ends, the CPU 170 extracts the first centering value Za0 (n) at the nth sampling timing (step S58). Hereinafter, specific contents of the first centering value extraction process executed by the CPU 170 in step S58 will be described with reference to FIG. FIG. 19 is a flowchart showing specific contents of the first centering value extraction process. As shown in FIG. 19, first, the CPU 170 determines whether or not the respiration speed determination flag Fma (n) is “10” (step S401). In other words, the CPU 170 determines whether or not the subject's breathing at the n-th sampling timing is fast breathing.

本実施形態では、被験者の呼吸速度に応じて可変に設定されるセンタリング期間内の複数のサンプリングタイミングの各々における第1生体電気インピーダンスZaの測定値に基づいて、第n番目のサンプリングタイミングにおける第1センタリング値Za0(n)が生成される。上記ステップS401の結果が肯定である場合、つまりは第n番目のサンプリングタイミングにおける被験者の呼吸が速めの呼吸である場合は、その速めの1呼吸に要する時間長(ここでは約4.0秒)がセンタリング期間として設定される。すなわち、被験者の呼吸が速めの呼吸である場合は、第n−9番目のサンプリングを始点とし、第n番目のサンプリングタイミングを終点とする期間がセンタリング期間として設定され、第n−9番目〜第n番目のサンプリングタイミングの各々における第1生体電気インピーダンスの測定値を用いた移動平均処理の結果に基づいて第1センタリング値Za0(n)が生成される。より具体的には、ステップS401の結果が肯定である場合、CPU170は、図15のステップS51で求めたMA10(n)に基づいて、第n番目のサンプリングタイミングにおける第1センタリング値Za0(n)を生成する(ステップS402)。さらに詳述すると、CPU170は、第n−2番目のサンプリングタイミングにおける第1センタリング値Za0(n−2)と、第n−1番目のサンプリングタイミングにおける第1センタリング値Za0(n−1)と、MA10(n)とを用いて移動平均処理を行い、その処理結果を、第n番目のサンプリングタイミングにおける第1センタリング値Za0(n)として採用する([Za0(n−2)+Za0(n−1)+MA10(n)]/3→Za0(n))。   In this embodiment, based on the measured value of the first bioelectrical impedance Za at each of a plurality of sampling timings within a centering period that is variably set according to the breathing rate of the subject, the first at the nth sampling timing. A centering value Za0 (n) is generated. If the result of step S401 is affirmative, that is, if the subject's breathing at the n-th sampling timing is a fast breath, the time length required for one fast breath (here, approximately 4.0 seconds) is centered. Set as period. That is, when the subject's breathing is fast breathing, a period from the n-9th sampling as the start point and the nth sampling timing as the end point is set as the centering period, and the n-9th to A first centering value Za0 (n) is generated based on the result of the moving average process using the measurement value of the first bioelectrical impedance at each of the nth sampling timings. More specifically, when the result of step S401 is affirmative, the CPU 170 determines the first centering value Za0 (n) at the nth sampling timing based on MA10 (n) obtained in step S51 of FIG. Is generated (step S402). More specifically, the CPU 170 has a first centering value Za0 (n-2) at the (n-2) th sampling timing, a first centering value Za0 (n-1) at the (n-1) th sampling timing, The moving average processing is performed using MA10 (n), and the processing result is adopted as the first centering value Za0 (n) at the nth sampling timing ([Za0 (n−2) + Za0 (n−1). ) + MA10 (n)] / 3 → Za0 (n)).

一方、上記ステップS401の結果が否定である場合、つまりは第n番目のサンプリングタイミングにおける被験者の呼吸が遅めの呼吸である場合は、その遅めの1呼吸に要する時間長(ここでは約8.0秒)がセンタリング期間として設定される。すなわち、被験者の呼吸が遅めの呼吸である場合は、第n−19番目のサンプリングを始点とし、第n番目のサンプリングタイミングを終点とする期間がセンタリング期間として設定され、第n−19番目〜第n番目のサンプリングタイミングの各々における第1生体電気インピーダンスの測定値を用いた移動平均処理の結果に基づいて第1センタリング値Za0(n)が生成される。より具体的には、ステップS401の結果が否定である場合、CPU170は、図15のステップS52で求めたMA20(n)に基づいて、第n番目のサンプリングタイミングにおける第1センタリング値Za0(n)を生成する(ステップS403)。さらに詳述すると、CPU170は、Za0(n−2)と、Za0(n−1)と、MA20(n)とを用いて移動平均処理を行い、その処理結果を、第n番目のサンプリングタイミングにおける第1センタリング値Za0(n)として採用する([Za0(n−2)+Za0(n−1)+MA20(n)]/3→Za0(n))。以上で、図14のステップS50における第1センタリング処理が終了する。   On the other hand, if the result of step S401 is negative, that is, if the subject's breathing at the n-th sampling timing is a slow breath, the length of time required for that slow breath (here, approximately 8.0) Second) is set as the centering period. That is, when the subject's breathing is slow breathing, the period from the n-19th sampling as the start point and the nth sampling timing as the end point is set as the centering period. A first centering value Za0 (n) is generated based on the result of the moving average process using the measurement value of the first bioelectrical impedance at each of the nth sampling timings. More specifically, if the result of step S401 is negative, the CPU 170 determines the first centering value Za0 (n) at the nth sampling timing based on MA20 (n) obtained in step S52 of FIG. Is generated (step S403). More specifically, the CPU 170 performs moving average processing using Za0 (n-2), Za0 (n-1), and MA20 (n), and the processing result is obtained at the nth sampling timing. Adopted as the first centering value Za0 (n) ([Za0 (n−2) + Za0 (n−1) + MA20 (n)] / 3 → Za0 (n)). This completes the first centering process in step S50 of FIG.

前述したように、被験者の呼吸が胸式呼吸であっても腹式呼吸であっても、第1生体電気インピーダンスZaの変化を示す波形は略正弦波状となる。CPU170は、体動などによる影響で、第1生体電気インピーダンスZaの変化を示す波形が乱れても、それに応じた第1センタリング値Za0が得られるように、所定数のサンプリングタイミングの各々における第1生体電気インピーダンスZaの測定値に基づいて第1センタリング値Za0を生成する。より具体的には、CPU170は、サンプリングタイミングごとに、当該サンプリングタイミングよりも所定の時間長だけ前の時点を始点とし、当該サンプリングタイミングを終点とするセンタリング期間内の複数のサンプリングタイミングの各々における第1生体電気インピーダンスの測定値を用いた移動平均処理を行い、その結果に基づいて、当該サンプリングタイミングにおける第1センタリング値Za0を求めるので、体動などによる影響で、第1生体電気インピーダンスZaの変化を示す波形が乱れても、それに応じた第1センタリング値Za0を精度良く生成できる。そして、各サンプリングタイミングに対応するセンタリング期間の時間長は、当該サンプリングタイミングにおける被験者の呼吸速度に応じて可変に設定されるという具合である。   As described above, the waveform indicating the change in the first bioelectrical impedance Za is substantially sinusoidal regardless of whether the subject's breathing is chest breathing or abdominal breathing. Even if the waveform indicating the change in the first bioelectric impedance Za is disturbed due to the influence of body movement or the like, the CPU 170 obtains the first centering value Za0 according to the first waveform at each of the predetermined number of sampling timings. A first centering value Za0 is generated based on the measured value of the bioelectric impedance Za. More specifically, for each sampling timing, the CPU 170 starts at a time point that is a predetermined time length before the sampling timing and starts at the sampling timing in each of a plurality of sampling timings within the centering period with the sampling timing as an end point. Since the moving average process using the measurement value of one bioelectrical impedance is performed and the first centering value Za0 at the sampling timing is obtained based on the result, the change in the first bioelectrical impedance Za due to the influence of body movement or the like Even if the waveform indicating is disturbed, the first centering value Za0 corresponding to the waveform can be generated with high accuracy. The time length of the centering period corresponding to each sampling timing is variably set according to the breathing rate of the subject at the sampling timing.

図14に戻って説明を続ける。図14に示すように、ステップS50の第1センタリング処理が終了すると、CPU170は、第1生体電気インピーダンスの測定値Za(n)の第1センタリング値Za0(n)に対する相対値である第1相対値ΔZa(n)を算出する第1相対値算出処理を実行する(ステップS60)。より具体的には、CPU170は、ステップS40で求めた第1生体電気インピーダンスの測定値Za(n)と、ステップS50で求めた第1センタリング値Za0(n)との差分を求め、その求めた差分値を、第n番目のサンプリングタイミングにおける第1相対値ΔZa(n)として採用するという具合である。   Returning to FIG. 14, the description will be continued. As illustrated in FIG. 14, when the first centering process in step S50 is completed, the CPU 170 performs a first relative that is a relative value of the first bioelectrical impedance measurement value Za (n) with respect to the first centering value Za0 (n). A first relative value calculation process for calculating the value ΔZa (n) is executed (step S60). More specifically, the CPU 170 obtains the difference between the first bioelectrical impedance measurement value Za (n) obtained in step S40 and the first centering value Za0 (n) obtained in step S50. The difference value is employed as the first relative value ΔZa (n) at the nth sampling timing.

ステップS60の後、CPU170は、第2生体電気インピーダンスZbの測定値の振幅基準レベルを示す第2センタリング値Zb0を生成し、その生成した第2センタリング値Zb0を用いて、第2生体電気インピーダンスZbの測定値の第2センタリング値Zb0に対する相対値である第2相対値ΔZbを算出する(ステップS70)。   After step S60, the CPU 170 generates a second centering value Zb0 indicating the amplitude reference level of the measurement value of the second bioelectrical impedance Zb, and uses the generated second centering value Zb0 to generate the second bioelectrical impedance Zb. A second relative value ΔZb, which is a relative value of the measured value to the second centering value Zb0, is calculated (step S70).

ここで、前述したように、腹式呼吸の呼気における第2生体電気インピーダンスZbの変化は、第1生体電気インピーダンスZaの変化とは異なる態様を示すので、第1センタリング値Za0を求める場合と同様に、所定数のサンプリングタイミングの各々における第2生体電気インピーダンスZbの測定値を用いた移動平均処理を行っても、これらの振幅基準レベル、すなわち、第2生体電気インピーダンスZbの変化を示す波形から、呼吸に起因した情報を抽出するための振幅基準レベル(第2センタリング値Zb0)を精度良く求めることは困難である。   Here, as described above, the change in the second bioelectrical impedance Zb in the exhalation of abdominal breathing is different from the change in the first bioelectrical impedance Za, and thus is similar to the case of obtaining the first centering value Za0. In addition, even if the moving average process using the measurement value of the second bioelectric impedance Zb at each of the predetermined number of sampling timings is performed, the amplitude reference level, that is, the waveform indicating the change of the second bioelectric impedance Zb is used. It is difficult to accurately obtain an amplitude reference level (second centering value Zb0) for extracting information resulting from respiration.

そこで、本実施形態では、図20に示すように、第1生体電気インピーダンスZaの測定値と、第1センタリング値Za0とが等しくなるゼロクロスタイミングを抽出し、当該ゼロクロスタイミングにおける第2生体電気インピーダンスZbの測定値に基づいて第2センタリング値Zb0を生成している。これにより、第2生体電気インピーダンスZbの測定値の振幅基準レベルを精度良く抽出できる。図20は、第2センタリング値Zb0の生成方法を概念的に説明するための図である。以下では、図21を参照しながら、ステップS70でCPU170が実行する第2相対値算出処理の具体的な内容を説明する。   Therefore, in the present embodiment, as shown in FIG. 20, the zero cross timing at which the measured value of the first bioelectric impedance Za is equal to the first centering value Za0 is extracted, and the second bioelectric impedance Zb at the zero cross timing is extracted. The second centering value Zb0 is generated based on the measured value. Thereby, the amplitude reference level of the measured value of the second bioelectrical impedance Zb can be extracted with high accuracy. FIG. 20 is a diagram for conceptually explaining a method of generating the second centering value Zb0. Hereinafter, specific contents of the second relative value calculation process executed by the CPU 170 in step S70 will be described with reference to FIG.

図21は、第2相対値算出処理の具体的な内容を示すフローチャートである。図21に示すように、CPU170は、5個のサンプリングタイミングの各々における第1相対値ΔZaのうち最小の値を抽出する(ステップS71)。より具体的には、CPU170は、第n−4番目〜第n番目のサンプリングタイミングの各々における第1相対値ΔZaの絶対値|ΔZa|(|ΔZa(n−4)|〜|ΔZa(n)|)のうち最小の値を、第n番目のサンプリングタイミングにおけるクロスポイント判定値ΔMIN5(n)として抽出するという具合である。   FIG. 21 is a flowchart showing specific contents of the second relative value calculation process. As shown in FIG. 21, the CPU 170 extracts the minimum value from the first relative values ΔZa at each of the five sampling timings (step S71). More specifically, the CPU 170 determines the absolute value | ΔZa | (| ΔZa (n−4) | to | ΔZa (n) of the first relative value ΔZa at each of the (n−4) th to n-th sampling timings. The minimum value of |) is extracted as the cross-point determination value ΔMIN5 (n) at the n-th sampling timing.

ステップS71の後、CPU170は、直前のサンプリングタイミングにおけるクロスポイント判定値と、ステップS71で抽出したクロスポイント判定値とが等しく、且つ、直前のサンプリングタイミングがゼロクロスタイミングであるか否かを判定する(ステップS72)。より具体的には、CPU170は、第n−1番目のサンプリングタイミングにおけるクロスポイント判定値ΔMIN5(n−1)とΔMIN5(n)とが等しく、且つ、第n−1番目のサンプリングタイミングにおけるクロスポイント判定フラッグF2(n−1)が「+1」に設定されているか否かを判定する。クロスポイント判定フラッグF2(n−1)が「+1」に設定されている場合は、第n−1番目のサンプリングタイミングがゼロクロスタイミングであるとみなされる。なお、クロスポイント判定フラッグF2の初期値(デフォルト値)、つまりは第1番目のサンプリングタイミングにおけるクロスポイント判定フラッグF2(1)の値は「0」に設定されている。   After step S71, the CPU 170 determines whether or not the cross point determination value at the immediately preceding sampling timing is equal to the cross point determination value extracted at step S71 and the immediately preceding sampling timing is zero cross timing ( Step S72). More specifically, the CPU 170 determines that the cross point determination values ΔMIN5 (n−1) and ΔMIN5 (n) at the (n−1) th sampling timing are equal and the crosspoint at the (n−1) th sampling timing. It is determined whether or not the determination flag F2 (n−1) is set to “+1”. When the cross point determination flag F2 (n-1) is set to "+1", the (n-1) th sampling timing is regarded as the zero cross timing. The initial value (default value) of the cross point determination flag F2, that is, the value of the cross point determination flag F2 (1) at the first sampling timing is set to “0”.

ステップS72の結果が肯定の場合、CPU170は、第n番目のサンプリングタイミングはゼロクロスタイミングではないと判定し、クロスポイント判定フラッグF2(n)を「0」に設定してステップS74へ進む。ステップS72の結果が否定の場合、CPU170は、第n番目のサンプリングタイミングにおける第1相対値ΔZa(n)の絶対値が0.3以下であるか否かを判定する(ステップS73)。ステップS73の結果が否定の場合、CPU170は、第n番目のサンプリングタイミングはゼロクロスタイミングではないと判定し、クロスポイント判定フラッグF2(n)を「0」に設定してステップS74へ進む。ステップS73の結果が肯定の場合、CPU170は、第n番目のサンプリングタイミングはゼロクロスタイミングであると判定し、クロスポイント判定フラッグF2(n)を「+1」に設定してステップS74へ進む。   If the result of step S72 is affirmative, the CPU 170 determines that the nth sampling timing is not the zero cross timing, sets the cross point determination flag F2 (n) to “0”, and proceeds to step S74. If the result of step S72 is negative, the CPU 170 determines whether or not the absolute value of the first relative value ΔZa (n) at the nth sampling timing is 0.3 or less (step S73). If the result of step S73 is negative, the CPU 170 determines that the nth sampling timing is not the zero cross timing, sets the cross point determination flag F2 (n) to “0”, and proceeds to step S74. If the result of step S73 is affirmative, the CPU 170 determines that the nth sampling timing is zero cross timing, sets the cross point determination flag F2 (n) to “+1”, and proceeds to step S74.

次に、CPU170は、第n番目のサンプリングタイミングにおけるクロスポイント判定フラッグF2(n)が「+1」であるか否かを判定する(ステップS74)。ステップS74の結果が肯定の場合、CPU170は、第2センタリング値Zb0を抽出する(ステップS75)。より具体的には、CPU170は、第n番目のサンプリングタイミングにおける第2生体電気インピーダンスの測定値Zb(n)を、第2センタリング値Zb0(n)として抽出する。そして、その抽出した第2センタリング値Zb0(n)と、第n−2番目のサンプリングタイミングにおける第2センタリング値Zb0(n−2)と、第n−1番目のサンプリングタイミングにおける第2センタリング値Zb0(n−1)とを用いて移動平均処理を行い、その処理結果を、正規の第2センタリング値Zb0(n)として生成する([Zb0(n−2)+Zb0(n−1)+Zb0(n)]/3→Zb0(n))。一方、ステップS74の結果が否定の場合、CPU170は、直前の第n−1番目のサンプリングタイミングにおける第2センタリング値Zb0(n−1)を、第n番目のサンプリングタイミングにおける第2センタリング値Zb0(n)として採用する(Zb0(n−1)→Zb0(n))。   Next, the CPU 170 determines whether or not the cross point determination flag F2 (n) at the nth sampling timing is “+1” (step S74). If the result of step S74 is affirmative, the CPU 170 extracts the second centering value Zb0 (step S75). More specifically, the CPU 170 extracts the measured value Zb (n) of the second bioelectrical impedance at the nth sampling timing as the second centering value Zb0 (n). Then, the extracted second centering value Zb0 (n), the second centering value Zb0 (n-2) at the (n-2) th sampling timing, and the second centering value Zb0 at the (n-1) th sampling timing. (N-1) and a moving average process is performed, and the process result is generated as a normal second centering value Zb0 (n) ([Zb0 (n-2) + Zb0 (n-1) + Zb0 (n )] / 3 → Zb0 (n)). On the other hand, if the result of step S74 is negative, the CPU 170 changes the second centering value Zb0 (n−1) at the immediately preceding n−1th sampling timing to the second centering value Zb0 ( n) (Zb0 (n−1) → Zb0 (n)).

そして、CPU170は、第n番目のサンプリングタイミングにおける第2相対値ΔZb(n)を算出する(ステップS76)。より具体的には、CPU170は、第2生体電気インピーダンスの測定値Zb(n)と、第2センタリング値Zb0(n)との差分を、第n番目のサンプリングタイミングにおける第2相対値ΔZb(n)として採用するという具合である。以上で、図14のステップS70における第2相対値算出処理が終了する。   Then, the CPU 170 calculates a second relative value ΔZb (n) at the nth sampling timing (step S76). More specifically, the CPU 170 calculates the difference between the second bioelectrical impedance measurement value Zb (n) and the second centering value Zb0 (n) as the second relative value ΔZb (n at the nth sampling timing. ) Is adopted. This completes the second relative value calculation process in step S70 of FIG.

例えば被験者の呼吸が腹式呼吸であって、第1生体電気インピーダンスZaおよび第2生体電気インピーダンスZbが図9のように変化する場合、第1相対値ΔZaおよび第2相対値ΔZbの経時的変化を示す波形は、図22のようになる。第1相対値ΔZaの経時的変化を示す波形は、第1生体電気インピーダンスZaの測定値の振幅基準レベルである第1センタリング値Za0をゼロ基準とするものであり、第2相対値ΔZbの経時的変化を示す波形は、第2生体電気インピーダンスZbの測定値の振幅基準レベルである第2センタリング値Zb0をゼロ基準とするものである。両波形を重ねることで、吸気における第1相対値ΔZaおよび第2相対値ΔZbの各々の波形と、呼気における第1相対値ΔZaおよび第2相対値ΔZbの各々の波形とを判別することができる。本実施形態では、CPU170は、第1相対値ΔZaの値に基づいて、被験者の呼吸が吸気であるか呼気であるかを判定する。より具体的には、CPU170は、第1相対値ΔZaが正の値である場合は吸気であると判定し、負の値である場合は呼気であると判定するという具合である。また、吸気における第1相対値ΔZaおよび第2相対値ΔZbの各々の波形の振幅および積分値の違いに基づいて、前述の回帰式(2)の係数を補正してもよい。これにより、測定精度の向上が図られる。   For example, when the subject's breathing is abdominal breathing and the first bioelectrical impedance Za and the second bioelectrical impedance Zb change as shown in FIG. 9, the first time relative value ΔZa and the second relative value ΔZb change over time. The waveform showing is as shown in FIG. The waveform showing the temporal change of the first relative value ΔZa is based on the first centering value Za0, which is the amplitude reference level of the measured value of the first bioelectrical impedance Za, as the zero reference, and the second relative value ΔZb is elapsed with time. The waveform indicating the change in the state is based on the second centering value Zb0, which is the amplitude reference level of the measurement value of the second bioelectrical impedance Zb, as a zero reference. By superimposing both waveforms, it is possible to determine the waveforms of the first relative value ΔZa and the second relative value ΔZb in inhalation and the waveforms of the first relative value ΔZa and the second relative value ΔZb in expiration. . In the present embodiment, the CPU 170 determines whether the subject's breathing is inspiration or expiration based on the value of the first relative value ΔZa. More specifically, the CPU 170 determines that it is inspiration when the first relative value ΔZa is a positive value, and determines that it is expiration when the first relative value ΔZa is a negative value. Further, the coefficient of the regression equation (2) described above may be corrected based on the difference in the amplitude and integral value of each waveform of the first relative value ΔZa and the second relative value ΔZb in the intake air. Thereby, the measurement accuracy is improved.

再び図14に戻って説明を続ける。図14に示すように、ステップ70の第2相対値算出処理が終了すると、CPU170は、第1相対値ΔZa(n)および第2相対値ΔZb(n)に対応するΔRib/ΔAbを推定するΔRib/ΔAb推定演算処理を実行する(ステップS80)。以下、図23および図24を参照しながら、ステップS80でCPU170が実行するΔRib/ΔAb推定演算処理の具体的な内容を説明する。図23および図24は、ΔRib/ΔAb推定演算処理の具体的な内容を示すフローチャートである。図23に示すように、まずCPU170は、第n番目のサンプリングタイミングにおける第1相対値ΔZa(n)の値が「0」以上であるか否かを判定する(ステップS81)。   Returning to FIG. 14 again, the description will be continued. As shown in FIG. 14, when the second relative value calculation process in step 70 is completed, the CPU 170 estimates ΔRib / ΔAb corresponding to the first relative value ΔZa (n) and the second relative value ΔZb (n). / ΔAb estimation calculation processing is executed (step S80). The specific contents of the ΔRib / ΔAb estimation calculation process executed by the CPU 170 in step S80 will be described below with reference to FIGS. 23 and 24 are flowcharts showing specific contents of the ΔRib / ΔAb estimation calculation process. As shown in FIG. 23, first, the CPU 170 determines whether or not the value of the first relative value ΔZa (n) at the n-th sampling timing is “0” or more (step S81).

ステップS81の結果が否定の場合、CPU170は、被験者の呼吸が呼気であると判断し、第n番目のサンプリングタイミングにおけるΔRib/ΔAb(n)の推定演算を行う(ステップS82)。より具体的には、CPU170は、上述の回帰式(2)にしたがって演算処理を実行することで、第1相対値ΔZa(n)および第2相対値ΔZb(n)に対応するΔRib/ΔAb(n)を求める。一方、ステップS81の結果が肯定の場合、CPU170は、被験者の呼吸が吸気であると判断し、第n番目のサンプリングタイミングにおけるΔRib/ΔAb(n)の値を初期値に設定する。本実施形態では、ステップS81の結果が肯定の場合、CPU170は、第n番目のサンプリングタイミングにおけるΔRib/ΔAb(n)の値を、初期値である「1.0」に設定する。   If the result of step S81 is negative, the CPU 170 determines that the subject's breath is exhalation, and performs an estimation calculation of ΔRib / ΔAb (n) at the nth sampling timing (step S82). More specifically, the CPU 170 executes arithmetic processing according to the above-described regression equation (2), whereby ΔRib / ΔAb (corresponding to the first relative value ΔZa (n) and the second relative value ΔZb (n). n). On the other hand, if the result of step S81 is affirmative, the CPU 170 determines that the subject's breathing is inspiration, and sets the value of ΔRib / ΔAb (n) at the nth sampling timing to an initial value. In this embodiment, when the result of step S81 is positive, the CPU 170 sets the value of ΔRib / ΔAb (n) at the n-th sampling timing to “1.0” which is an initial value.

次に、CPU170は、ΔRib/ΔAb(n)の値が、−2.5以上であり、且つ4.5以下であるか否かを判定する(ステップS83)。ステップS83の結果が否定の場合、CPU170は、ΔRib/ΔAbの値を、初期値である「1.0」に設定してステップS84へ進む。ステップS83の結果が肯定の場合、CPU170は、そのままステップS84へ進む。   Next, the CPU 170 determines whether or not the value of ΔRib / ΔAb (n) is −2.5 or more and 4.5 or less (step S83). If the result of step S83 is negative, the CPU 170 sets the value of ΔRib / ΔAb to the initial value “1.0” and proceeds to step S84. If the result of step S83 is affirmative, the CPU 170 proceeds to step S84 as it is.

ステップS84において、CPU170は、ΔRib/ΔAb(n)の絶対値と、直前の第n−1番目のサンプリングタイミングにおけるΔRib/ΔAb(n−1)の絶対値との差分(|ΔRib/ΔAb(n)|−|ΔRib/ΔAb(n−1)|)が0.3よりも大きいか否かを判定する。ステップS84の結果が否定である場合、CPU170は、直前の第n−1番目のサンプリングタイミングにおけるΔRib/ΔAb(n−1)と、ΔRib/ΔAbとの平均を求め、その求めた平均値を、第n番目のサンプリングタイミングにおけるΔRib/ΔAb(n)として採用([ΔRib/ΔAb(n−1)+ΔRib/ΔAb(n)]/2→ΔRib/ΔAb(n))して次のステップS85(図24参照)へ進む。一方、ステップS84の結果が肯定である場合、CPU170は、第n番目のサンプリングタイミングにおけるΔRib/ΔAb(n)の値を、初期値である「1.0」に設定して次のステップS85へ進む。   In step S84, the CPU 170 determines the difference (| ΔRib / ΔAb (n) between the absolute value of ΔRib / ΔAb (n) and the absolute value of ΔRib / ΔAb (n−1) at the immediately preceding (n−1) th sampling timing. ) | − | ΔRib / ΔAb (n−1) |) is determined as to whether it is greater than 0.3. When the result of step S84 is negative, the CPU 170 obtains the average of ΔRib / ΔAb (n−1) and ΔRib / ΔAb at the immediately preceding (n−1) -th sampling timing, and calculates the obtained average value. Adopted as ΔRib / ΔAb (n) at the n-th sampling timing ([ΔRib / ΔAb (n−1) + ΔRib / ΔAb (n)] / 2 → ΔRib / ΔAb (n)) and the next step S85 (FIG. 24). On the other hand, if the result of step S84 is affirmative, the CPU 170 sets the value of ΔRib / ΔAb (n) at the n-th sampling timing to the initial value “1.0” and proceeds to the next step S85.

図24を参照しながら、ΔRib/ΔAb推定演算処理の具体的な内容の説明を続ける。図24に示すように、ステップS85において、CPU170は、第n番目のサンプリングタイミングにおける第1相対値ΔZa(n)が「0」よりも小さいか否かを判定する。ステップS85の結果が肯定の場合、CPU170は、被験者の呼吸が呼気であると判断し、直前の積分回数のカウント値Niに1を加算する。より具体的には、CPU170は、直前の第n−1番目のサンプリングタイミングにおける積分回数のカウント値Ni(n−1)に1を加算した値を、第n番目のサンプリングタイミングにおける積分回数のカウント値Ni(n)として採用する。   With reference to FIG. 24, the description of the specific contents of the ΔRib / ΔAb estimation calculation process will be continued. As shown in FIG. 24, in step S85, the CPU 170 determines whether or not the first relative value ΔZa (n) at the nth sampling timing is smaller than “0”. If the result of step S85 is affirmative, the CPU 170 determines that the subject's breathing is exhalation, and adds 1 to the count value Ni of the previous integration count. More specifically, the CPU 170 counts the number of integrations at the nth sampling timing by adding a value obtained by adding 1 to the count value Ni (n−1) of the number of integrations at the immediately preceding n−1th sampling timing. The value Ni (n) is adopted.

一方、ステップS85の結果が否定の場合、CPU170は、被験者の呼吸が吸気であると判断し、積分回数のカウント値Niの値を「0」に初期化する。すなわち、この場合、第n番目のサンプリングタイミングにおける積分回数のカウント値Ni(n)は「0」に設定されるという具合である。   On the other hand, if the result of step S85 is negative, the CPU 170 determines that the subject's breathing is inspiration, and initializes the count value Ni of the number of integrations to “0”. That is, in this case, the count value Ni (n) of the number of integrations at the nth sampling timing is set to “0”.

続いて、図24に示すように、CPU170は、第1相対値ΔZa(n)が「0」よりも小さいか否かを再び判定する(ステップS86)。ステップS86の結果が肯定の場合、CPU170は、直前の第n−1番目のサンプリングタイミングにおけるΔRib/ΔAbの積分値(ΣΔRib/ΔAb(n−1))と、ΔRib/ΔAb(n)との和を求めることで、第n番目のサンプリングタイミングにおけるΔRib/ΔAbの積分値(ΣΔRib/ΔAb(n))を求めて次のステップS87へ進む。一方、ステップS86の結果が否定の場合、つまりは、被験者の呼吸状態が吸気であると判断した場合は、CPU170は、ピーク/ボトム判別フラッグF1(n)を「+1」に設定して次のステップS87へ進む。   Subsequently, as illustrated in FIG. 24, the CPU 170 determines again whether or not the first relative value ΔZa (n) is smaller than “0” (step S86). When the result of step S86 is affirmative, the CPU 170 calculates the sum of ΔRib / ΔAb integral value (ΣΔRib / ΔAb (n−1)) and ΔRib / ΔAb (n) at the immediately preceding n−1th sampling timing. Is obtained, and the integrated value (ΣΔRib / ΔAb (n)) of ΔRib / ΔAb at the n-th sampling timing is obtained, and the process proceeds to the next step S87. On the other hand, if the result of step S86 is negative, that is, if it is determined that the subject's breathing state is inhalation, the CPU 170 sets the peak / bottom determination flag F1 (n) to “+1” and then continues. The process proceeds to step S87.

次に、CPU170は、第n番目のサンプリングタイミングにおける積分回数のカウント値Ni(n)がゼロであるか否かを判定する(ステップS87)。ステップS87の結果が否定の場合、CPU170は、第n番目のサンプリングタイミングにおけるΔRib/ΔAbの積分値[ΣΔRib/ΔAb(n)]を、第n番目のサンプリングタイミングにおける積分回数のカウント値Ni(n)で割ることで、ΔRib/ΔAbの平均値を求める。一方、ステップS87の結果が肯定の場合、CPU170は、直前の第n−1番目のサンプリングタイミングにおけるΔRib/ΔAbの平均値([ΣΔRib/ΔAb(n−1)]/Ni(n−1))を、第n番目のサンプリングタイミングにおけるΔRib/ΔAbの平均値として採用する。以上で、図14のステップS80におけるΔRib/ΔAb推定演算処理が終了する。   Next, the CPU 170 determines whether or not the count value Ni (n) of the number of integrations at the nth sampling timing is zero (step S87). If the result of step S87 is negative, the CPU 170 obtains the integrated value [ΣΔRib / ΔAb (n)] of ΔRib / ΔAb at the nth sampling timing, and the count value Ni (n of the number of integrations at the nth sampling timing). ) To obtain an average value of ΔRib / ΔAb. On the other hand, when the result of step S87 is affirmative, the CPU 170 determines the average value of [Delta] Rib / [Delta] Ab at the immediately preceding (n-1) th sampling timing ([[Sigma] [Delta] Rib / [Delta] Ab (n-1)] / Ni (n-1)). Is adopted as the average value of ΔRib / ΔAb at the nth sampling timing. This completes the ΔRib / ΔAb estimation calculation process in step S80 of FIG.

例えば被験者の呼吸が腹式呼吸であって、第1相対値ΔZaおよび第2相対値ΔZbが図22のように変化する場合、上述のΔRib/ΔAbの平均値([ΣΔRib/ΔAb]/Ni)は、図25のように変化する。図25において、呼気におけるΔRib/ΔAbの値(平均値)が1.0よりも小さい場合は、被験者の呼吸は腹式呼吸であると推定される一方、1.0を上回る場合は、被験者の呼吸は胸式呼吸であると推定されるという具合である。   For example, when the subject's breathing is abdominal breathing and the first relative value ΔZa and the second relative value ΔZb change as shown in FIG. 22, the average value of the above-described ΔRib / ΔAb ([ΣΔRib / ΔAb] / Ni) Changes as shown in FIG. In FIG. 25, when the value (average value) of ΔRib / ΔAb in exhalation is smaller than 1.0, the subject's breathing is estimated to be abdominal breathing. Respiration is estimated to be chest respiration.

再び図14に戻って説明を続ける。図14に示すように、ステップS80のΔRib/ΔAb推定演算処理が終了すると、CPU170は、被験者の呼吸の深さを示す呼吸レベルを抽出する呼吸レベル抽出処理を実行する(ステップS90)。以下、図26を参照しながら、ステップS90でCPU170が実行する呼吸レベル抽出処理の具体的な内容を説明する。図26は、呼吸レベル抽出処理の具体的な内容を示すフローチャートである。図26に示すように、まずCPU170は、サンプリングタイミングに到達したか否かを判定し(ステップS501)、ステップS501の結果が肯定である場合は次のステップS502へ進む。ステップS502において、CPU170は、被験者の呼吸が吸気であるか否かを判定する。具体的には、CPU170は、当該サンプリングタイミングにおける第1相対値ΔZaの値が正の値である場合には吸気であると判定する一方、負の値である場合には呼気であると判定するという具合である。CPU170は、吸気であると判定した場合は、吸気判別フラッグF3を「1」に設定する一方、呼気であると判定した場合は、吸気判別フラッグF3を「0」に設定する。なお、初期状態においては、吸気判別フラッグF3は「1」に設定される(つまりはF3のデフォルト値は1に設定される)。   Returning to FIG. 14 again, the description will be continued. As shown in FIG. 14, when the ΔRib / ΔAb estimation calculation process in step S80 ends, the CPU 170 executes a breathing level extraction process for extracting a breathing level indicating the breathing depth of the subject (step S90). Hereinafter, the specific content of the respiration level extraction process executed by the CPU 170 in step S90 will be described with reference to FIG. FIG. 26 is a flowchart showing specific contents of the respiration level extraction process. As shown in FIG. 26, first, the CPU 170 determines whether or not the sampling timing has been reached (step S501), and if the result of step S501 is affirmative, the process proceeds to the next step S502. In step S502, the CPU 170 determines whether or not the subject's breathing is inspiration. Specifically, the CPU 170 determines that it is inspiration when the value of the first relative value ΔZa at the sampling timing is a positive value, and determines that it is expiration when the value is a negative value. That's it. If the CPU 170 determines that it is inhalation, it sets the inhalation determination flag F3 to “1”, whereas if it determines that it is expiration, it sets the inspiration determination flag F3 to “0”. In the initial state, the intake determination flag F3 is set to “1” (that is, the default value of F3 is set to 1).

ステップS502の結果が肯定の場合、CPU170は、ピークホールド処理を実行する(ステップS503)。より具体的には、CPU170は、吸気での複数のサンプリングタイミングの各々における第1相対値ΔZaのうち最大の値をピーク値ΔZa(MAX)として保持する。一方、ステップS502の結果が否定の場合、CPU170は、ボトムホールド処理を実行する(ステップS504)。より具体的には、CPU170は、呼気での複数のサンプリングタイミングの各々における第1相対値ΔZaのうち最小の値をボトム値ΔZa(MIN)として保持する。   If the result of step S502 is positive, the CPU 170 executes a peak hold process (step S503). More specifically, the CPU 170 holds the maximum value as the peak value ΔZa (MAX) among the first relative values ΔZa at each of a plurality of sampling timings during intake. On the other hand, if the result of step S502 is negative, the CPU 170 executes bottom hold processing (step S504). More specifically, the CPU 170 holds the minimum value among the first relative values ΔZa at each of a plurality of sampling timings during expiration as the bottom value ΔZa (MIN).

次に、CPU170は、当該サンプリングタイミングがゼロクロスタイミングであるか否かを判定する(ステップS505)。より具体的には、CPU170は、当該サンプリングタイミングにおけるクロスポイント判定フラッグF2が「+1」であるか否かを判定することで、当該サンプリングタイミングがゼロクロスタイミングであるか否かを判定するという具合である。ステップS505の結果が否定の場合、当該サンプリングタイミングにおける呼吸レベル抽出処理は終了する。一方、ステップS505の結果が肯定の場合、CPU170は、当該サンプリングタイミングにおける第1生体電気インピーダンスZaの微分係数dZaが正極性(>0)であるか否かを判定する(ステップS506)。言い換えれば、CPU170は、当該サンプリングタイミングが、呼気から吸気へと変化するタイミングであるか否かを判定する。   Next, the CPU 170 determines whether or not the sampling timing is a zero cross timing (step S505). More specifically, the CPU 170 determines whether or not the sampling timing is the zero cross timing by determining whether or not the cross point determination flag F2 at the sampling timing is “+1”. It is. If the result of step S505 is negative, the respiration level extraction process at the sampling timing ends. On the other hand, when the result of step S505 is positive, the CPU 170 determines whether or not the differential coefficient dZa of the first bioelectric impedance Za at the sampling timing is positive (> 0) (step S506). In other words, the CPU 170 determines whether or not the sampling timing is a timing for changing from expiration to inspiration.

ステップS506の結果が肯定の場合、CPU170は、当該サンプリングタイミングが、呼気から吸気へと変化するタイミングであると判断して、そのときホールドされているピーク値ΔZa(MAX)とボトム値ΔZa(MIN)との絶対値の和を、直前の1呼吸における呼吸レベルΔZap−pとして抽出する(ステップS507)。その後、CPU170は、吸気フラッグ設定処理を実行する(ステップS508)。より具体的には、CPU170は、吸気判別フラッグF3を「1」に設定する。そして、CPU170は、ピークホールド処理を初期化する(ステップS509)。より具体的には、CPU170は、ステップS503で保持していたピーク値ΔZa(MAX)を、「0」に設定(初期化)して、当該サンプリングタイミングにおける呼吸レベル抽出処理を終了する。   If the result of step S506 is affirmative, the CPU 170 determines that the sampling timing is a timing at which the expiration changes from inspiration to inspiration, and the peak value ΔZa (MAX) and bottom value ΔZa (MIN) held at that time are determined. ) And the absolute value is extracted as the respiration level ΔZap-p in the immediately preceding breath (step S507). Thereafter, the CPU 170 executes an intake flag setting process (step S508). More specifically, the CPU 170 sets the intake determination flag F3 to “1”. Then, the CPU 170 initializes the peak hold process (step S509). More specifically, the CPU 170 sets (initializes) the peak value ΔZa (MAX) held in step S503 to “0”, and ends the respiration level extraction process at the sampling timing.

一方、ステップS506の結果が否定の場合、CPU170は、当該サンプリングタイミングが、吸気から呼気へと変化するタイミングであると判断して、呼気フラッグ設定処理を実行する(ステップS510)。より具体的には、CPU170は、吸気判別フラッグF3を「0」に設定する。そして、CPU170は、ボトムホールド処理を初期化する(ステップS511)。より具体的には、CPU170は、ステップS504で保持していたボトム値ΔZa(MIN)を、「0」に設定(初期化)して、当該サンプリングタイミングにおける呼吸レベル抽出処理を終了する。以上で、呼吸レベル抽出処理が終了し、第n番目のサンプリングタイミングにおける呼吸解析処理が終了する。   On the other hand, if the result of step S506 is negative, the CPU 170 determines that the sampling timing is a timing at which the sampling changes from inspiration to expiration, and executes an expiration flag setting process (step S510). More specifically, the CPU 170 sets the intake determination flag F3 to “0”. Then, the CPU 170 initializes the bottom hold process (step S511). More specifically, the CPU 170 sets (initializes) the bottom value ΔZa (MIN) held in step S504 to “0”, and ends the respiration level extraction process at the sampling timing. Thus, the respiration level extraction process ends, and the respiration analysis process at the nth sampling timing ends.

以上に説明したように、CPU170が呼吸解析処理を実行することで、サンプリングタイミングごとに、当該サンプリングタイミングにおける第1相対値ΔZaおよび第2相対値ΔZbに対応するΔRib/ΔAbが求められるので、被験者の呼吸の種別(腹式呼吸なのか胸式呼吸なのか)をリアルタイムで正確に判別できる。   As described above, since the CPU 170 executes the breath analysis process, ΔRib / ΔAb corresponding to the first relative value ΔZa and the second relative value ΔZb at the sampling timing is obtained for each sampling timing. Type of breathing (whether abdominal breathing or chest breathing) can be accurately determined in real time.

<1−5:呼吸レベル表示処理>
次に、CPU170が実行する呼吸レベル表示処理について説明する。本実施形態では、CPU170は、被験者の1呼吸ごとに、当該1呼吸における腹式呼吸および胸式呼吸の各々の大きさと余裕度とを表示(報知)する呼吸レベル表示処理を実行する。より具体的には、CPU170は、被験者の1呼吸ごとに、当該1呼吸における呼吸の深さを示す呼吸レベルと、当該1呼吸における呼吸解析処理の結果とから、当該1呼吸における胸式呼吸および腹式呼吸の各々の大きさと余裕度とを表示部160に表示するように制御する。図27に示すように、本実施形態では、1呼吸における胸式呼吸および腹式呼吸の各々の大きさと余裕度のほか、当該1呼吸に占める腹式呼吸の割合を示す腹式レベルを表示部160に表示するように制御する。図27に示す第1バーグラフBG1は、胸式呼吸および腹式呼吸の各々の大きさと余裕度とを表示するためのものである。一方、図27に示す第2バーグラフBG2は、被験者の腹式レベルを表示するためのものである。これらの詳細な内容については後述する。
<1-5: Respiration level display process>
Next, a respiration level display process executed by the CPU 170 will be described. In the present embodiment, the CPU 170 executes a breathing level display process for displaying (notifying) the size and margin of each abdominal breathing and chest breathing in one breath for each breath of the subject. More specifically, for each breath of the subject, the CPU 170 calculates the chest breathing in the one breath and the respiratory level indicating the depth of the breath in the one breath and the result of the breath analysis process in the one breath. Control is performed so that the magnitude and margin of each abdominal breath are displayed on the display unit 160. As shown in FIG. 27, in this embodiment, in addition to the magnitude and margin of each of the chest breathing and abdominal breathing in one breath, the abdominal level indicating the proportion of the abdominal breathing in one breath is displayed on the display unit. The display is controlled to be displayed at 160. The first bar graph BG1 shown in FIG. 27 is for displaying the magnitude and margin of each of the chest breathing and the abdominal breathing. On the other hand, the second bar graph BG2 shown in FIG. 27 is for displaying the abdominal level of the subject. Details of these will be described later.

図28を参照しながら、CPU170が実行する呼吸レベル表示処理の具体的な内容を説明する。図28は、呼吸レベル表示処理の具体的な内容を示すフローチャートである。図28に示すように、CPU170は、直前の1呼吸における呼吸レベルを正規化する(ステップS601)。ここで、「呼吸レベルを正規化する」とは、前述の呼吸レベル抽出処理で抽出された直前の1呼吸における呼吸レベルの値を、被験者間における体格等の個人差の影響を取り除いた値に補正することを意味する。具体的には、ステップS601において、CPU170は、直前の1呼吸における呼吸レベルΔZap−pを、当該1呼吸における第2センタリング値Zb0(詳述すれば、当該1呼吸における複数のサンプリングタイミングの各々における第2センタリング値Zb0の平均値)で割った後に100を乗じて求めた値を、呼吸レベルの正規化値%ΔZa(=(ΔZap−p/Zb0)×100)として採用する。   The specific contents of the respiration level display process executed by the CPU 170 will be described with reference to FIG. FIG. 28 is a flowchart showing specific contents of the respiration level display process. As shown in FIG. 28, the CPU 170 normalizes the respiratory level in the immediately preceding one breath (step S601). Here, “normalize the respiratory level” is a value obtained by removing the influence of individual differences such as physique among subjects from the value of the respiratory level immediately before one extracted by the above-described respiratory level extraction processing. It means to correct. Specifically, in step S601, the CPU 170 sets the respiration level ΔZap-p in the immediately preceding one breath to the second centering value Zb0 in the one breath (more specifically, in each of a plurality of sampling timings in the one breath. A value obtained by dividing by 100 after being divided by the average value of the second centering value Zb0 is employed as a normalized value% ΔZa (= (ΔZap−p / Zb0) × 100) of the respiration level.

次に、CPU170は、ステップS601で求めた呼吸レベルの正規化値%ΔZaの値に応じて、着色表示すべき第1バーグラフBG1の段階数(以下、「第1表示段階数」と呼ぶ)を決定する(ステップS602)。より具体的には、CPU170は、ステップS601で求めた呼吸レベルの正規化値%ΔZaに対応する1回換気量の正規化値%ΔTVを求め、その求めた%ΔTVの値に応じて、第1表示段階数を決定する。ここで、「1回換気量」とは、被験者の1呼吸で肺に出入りする空気量を示すものであり、ΔTVと表記する。また、「1回換気量を正規化する」とは、スパイロなどで計測された1回換気量ΔTVの値(実測値)を、被験者間における体格等の個人差の影響を取り除いた値に補正することを意味する。本実施形態では、1回換気量ΔTVの正規化を行う際には、スパイロなどで計測された被験者の肺活量の値(実測値)VCを、標準肺活量を示す標準VCで割り算した形で表される係数%VC(=実測VC/標準VC)が用いられる。具体的には、計測された1回換気量ΔTVの値(実測値)を上述の係数%VCで割ることで、1回換気量の正規化値%ΔTV(=ΔTV/%VC)が求められる。なお、例えば男性の標準肺活量VCmale(ml)は、(27.63−0.112×年齢)×身長(cm)で表され、女性の標準肺活量VCfemale(ml)は、(21.78−0.101×年齢)×身長(cm)で表される。   Next, the CPU 170 determines the number of stages of the first bar graph BG1 to be colored and displayed (hereinafter referred to as “first display stage number”) according to the value of the normalized value% ΔZa of the respiratory level obtained in step S601. Is determined (step S602). More specifically, the CPU 170 obtains a normalized value% ΔTV of the tidal volume corresponding to the normalized value% ΔZa of the respiratory level obtained in step S601, and according to the obtained value of% ΔTV, The number of display stages is determined. Here, the “tidal volume” indicates the amount of air that enters and exits the lungs in one breath of the subject, and is expressed as ΔTV. “Normalizing the tidal volume” means correcting the tidal volume value ΔTV (measured value) measured by a spiro etc. to a value that removes the influence of individual differences such as physique between subjects. It means to do. In the present embodiment, when normalizing the tidal volume ΔTV, the value (actually measured value) VC of the subject's vital capacity measured by a spiro or the like is divided by the standard VC indicating the standard vital capacity. Coefficient% VC (= actual VC / standard VC) is used. Specifically, the normalized value% ΔTV (= ΔTV /% VC) of the tidal volume is obtained by dividing the measured value (actually measured value) of the tidal volume ΔTV by the above-described coefficient% VC. . For example, male standard vital capacity VCmale (ml) is expressed by (27.63-0.112 × age) × height (cm), and female standard vital capacity VCfemale (ml) is (21.78-0. 101 × age) × height (cm).

図29は、男女合わせて20人で、各々の1回換気量ΔTVを3回ずつ(小程度の換気、中程度の換気、大程度の換気を各1回)測定したときの、1回換気量の正規化値%ΔTVと呼吸レベルの正規化値%ΔZaとの関係を示す相関図である。図29に示すように、1回換気量の正規化値%ΔTVと呼吸レベルの正規化値%ΔZaとは相関係数r=0.75,P<0.01という高い相関が得られ、以下の回帰式(4)が成立する。
%ΔZa=c0×%ΔTV ・・・(4)
C0:回帰係数。
FIG. 29 shows tidal ventilation when 20 people in both men and women measure each tidal volume ΔTV three times (small ventilation, medium ventilation, and large ventilation once). It is a correlation diagram showing the relationship between the normalized value% ΔTV of the quantity and the normalized value% ΔZa of the respiratory level. As shown in FIG. 29, the normalized value% ΔTV of the tidal volume and the normalized value% ΔZa of the respiration level have a high correlation coefficient of r = 0.75 and P <0.01. The regression equation (4) is established.
% ΔZa = c0 ×% ΔTV (4)
C0: regression coefficient.

CPU170は、上述の回帰式(4)にしたがって演算処理を実行することで、ステップS601で求めた呼吸レベルの正規化値%ΔZaに対応する1回換気量の正規化値%ΔTVを求める。そして、CPU170は、その求めた1回換気量の正規化値%ΔTVに応じて、第1表示段階数を決定する。本実施形態では、第1表示段階数の最大値は、胸式5段階と腹式5段階とを合わせた「10」に設定されている(図27参照)。   The CPU 170 obtains a normalized value% ΔTV of the tidal volume corresponding to the normalized value% ΔZa of the respiratory level obtained in step S601 by executing a calculation process according to the regression equation (4) described above. Then, the CPU 170 determines the first display step number according to the obtained normalized value% ΔTV of the tidal volume. In the present embodiment, the maximum value of the number of first display steps is set to “10”, which is a combination of the five chest steps and the five abdominal steps (see FIG. 27).

図29に示すように、本実施形態では、1回換気量の正規化値%ΔTVの値が第1所定値α1以上である場合は、被験者の呼吸レベルは「最大」とみなされる。この場合、CPU170は、第1表示段階数を最大値である「10」に決定する。また、1回換気量の正規化値%ΔTVの値が第2所定値α2(<α1)以上であって、かつ第1所定値α1を下回る場合は、被験者の呼吸レベルは、「大程度」とみなされる(図29参照)。この場合、CPU170は、第1表示段階数を「8」に決定する。さらに、1回換気量の正規化値%ΔTVが第3所定値α3(<α2)以上であって、かつ第2所定値α2を下回る場合は、被験者の呼吸レベルは「中程度」とみなされる(図29参照)。この場合、CPU170は、第1表示段階数を「6」に決定する。また、1回換気量の正規化値%ΔTVの値が第4所定値α4(<α3)以上であって、かつ第3所定値α3を下回る場合は、被験者の呼吸レベルは「小程度」とみなされる(図29参照)。この場合、CPU170は、第1表示段階数を「4」に決定する。また、1回換気量の正規化値%ΔTVの値が第4所定値α4を下回る場合は、被験者の呼吸レベルは、安静時における必須レベル程度とみなされる(図29参照)。この場合、CPU170は、第1表示段階数を「2」に決定するという具合である。   As shown in FIG. 29, in this embodiment, when the value of the normalized value% ΔTV of the tidal volume is equal to or greater than the first predetermined value α1, the breathing level of the subject is regarded as “maximum”. In this case, the CPU 170 determines the first display step number to be “10” which is the maximum value. Further, when the value of the normalized tidal volume value% ΔTV is equal to or larger than the second predetermined value α2 (<α1) and is lower than the first predetermined value α1, the respiration level of the subject is “large”. (See FIG. 29). In this case, the CPU 170 determines the first display step number as “8”. Further, when the normalized value% ΔTV of the tidal volume is equal to or larger than the third predetermined value α3 (<α2) and is lower than the second predetermined value α2, the subject's respiratory level is regarded as “medium”. (See FIG. 29). In this case, the CPU 170 determines the first display step number as “6”. When the value of the tidal volume normalized value% ΔTV is equal to or larger than the fourth predetermined value α4 (<α3) and lower than the third predetermined value α3, the subject's breathing level is “small”. Considered (see FIG. 29). In this case, the CPU 170 determines the first display step number as “4”. In addition, when the value of the normalized value% ΔTV of the tidal volume is lower than the fourth predetermined value α4, the subject's breathing level is regarded as an essential level at rest (see FIG. 29). In this case, the CPU 170 determines the number of first display steps as “2”.

再び図28に戻って説明を続ける。図28に示すように、ステップS602の段階数決定処理が終了すると、CPU170は、直前の1呼吸における腹式レベルを決定する(ステップS603)。具体的には、CPU170は、直前の1呼吸におけるΔRib/ΔAbの値(平均値)に応じて、腹式レベルの数値を決定する。腹式レベルの数値は0〜100までの値で表され、その値が0に近いほど腹式呼吸の程度は小さく(被験者の呼吸に占める腹式呼吸の割合は小さく)、100に近いほど腹式呼吸の程度は大きい(被験者の呼吸に占める腹式呼吸の割合は大きい)。   Returning to FIG. 28 again, the description will be continued. As shown in FIG. 28, when the stage number determination process in step S602 ends, the CPU 170 determines the abdominal level in the immediately preceding one breath (step S603). Specifically, CPU 170 determines the numerical value of the abdominal level according to the value (average value) of ΔRib / ΔAb in the immediately preceding one breath. The numerical value of the abdominal level is represented by a value from 0 to 100. The closer the value is to 0, the smaller the level of abdominal breathing (the proportion of the abdominal breathing in the subject's breathing is smaller). The degree of breathing is large (the proportion of abdominal breathing in the subject's breathing is large).

図28に示すように、ステップS603の後、CPU170は、ステップS602で決定した第1表示段階数を、腹式と胸式とに分割する段階数分割処理を実行する(ステップS604)。より具体的には、CPU170は、ステップS602で決定した第1表示段階数と、ステップS603で決定した腹式レベルとに基づいて、段階数分割処理を実行する。ここでは、全体の呼吸レベルを示す第1表示段階数が「6」である一方、腹式レベルは「70」である場合を想定して説明を続ける。腹式レベルが「70」であるとは、腹式呼吸と胸式呼吸との割合は、7:3であるとみなされるので、第1表示段階数「6」は、胸式で2段階、腹式で4段階に分割される。つまり、図27のように、胸式呼吸の大きさを示す第1バーグラフのBG1の表示段階数は「2」、腹式呼吸の大きさを示す第1バーグラフBG1の表示段階数は「4」になるように設定される。   As shown in FIG. 28, after step S603, the CPU 170 executes a stage number division process for dividing the first display stage number determined in step S602 into an abdominal type and a chest type (step S604). More specifically, the CPU 170 executes the stage number division process based on the first display stage number determined in step S602 and the abdominal level determined in step S603. Here, the description will be continued assuming that the first display stage number indicating the entire respiration level is “6” while the abdominal level is “70”. If the abdominal level is “70”, the ratio between the abdominal breathing and the chest breathing is considered to be 7: 3, so the first display stage number “6” Divided into 4 stages in the abdomen. That is, as shown in FIG. 27, the number of display stages of BG1 of the first bar graph indicating the magnitude of chest breathing is “2”, and the number of display stages of the first bar graph BG1 indicating the magnitude of abdominal breathing is “ 4 ".

図28に示すように、ステップS604の後、CPU170は、被験者の呼吸速度に応じたマージンレベルを決定する(ステップS605)。いま、直前の1呼吸が安静呼吸であって、当該1呼吸に要する時間長が4秒以上であった場合を想定する。本実施形態では、被験者の呼吸速度に応じた必須の呼吸レベルに対応する第1バーグラフBG1の段階数が予め定められており、安静呼吸時の呼吸速度に応じた必須の呼吸レベルに対応する第1バーグラフBG1の段階数は「2」に設定されている。CPU170は、その段階数「2」を、胸式と腹式に1段階ずつ割り振る。上述のステップS604で説明したように、ここでは、胸式呼吸の大きさを示す第1バーグラフBG1の表示段階数は「2」、腹式呼吸の大きさを示す第1バーグラフBG1の表示段階数は「4」となるように設定されているので、胸式呼吸については1段階分の余裕があり、腹式呼吸については3段階分の余裕があるという具合である。   As shown in FIG. 28, after step S604, the CPU 170 determines a margin level according to the breathing rate of the subject (step S605). Now, it is assumed that the last breath is a rest breath and the time length required for the breath is 4 seconds or more. In the present embodiment, the number of steps of the first bar graph BG1 corresponding to the essential respiratory level according to the breathing rate of the subject is determined in advance, and corresponds to the essential respiratory level according to the respiratory rate during rest breathing. The number of steps of the first bar graph BG1 is set to “2”. The CPU 170 assigns the stage number “2” to the breast type and the abdominal type one by one. As described in step S604 above, here, the number of display stages of the first bar graph BG1 indicating the magnitude of the chest breathing is “2”, and the display of the first bar graph BG1 indicating the magnitude of the abdominal breathing is displayed. Since the number of stages is set to “4”, there is a margin for one stage for chest breathing and a margin for three stages for abdominal breathing.

この場合は、腹式呼吸の大きさを示す第1バーグラフBG1の表示段階数が「1」以上であれば、被験者の腹式呼吸の大きさは必須レベルを満たしていることを意味し、その表示段階数が「2」であれば、被験者の腹式呼吸の大きさは必須レベルを上回る「小程度」のレベルを満たしていることを意味し、その表示段階数が「3」であれば、さらに上の「中程度」のレベルを満たしていることを意味し、その表示段階数が「4」であれば、さらに上の「大程度」のレベルを満たしていることを意味し、その表示段階数が「5」であれば、被験者の腹式呼吸の大きさは「最大」レベルを満たしていることを意味するという具合である。前述したように、ここでは、腹式呼吸の大きさを示す第1バーグラフBG1の表示段階数数は「4」に設定されているので、被験者の腹式呼吸の大きさは「大程度」のレベルを満たしており、必須レベル(表示段階数「1」)に対して十分な余裕があることが分かる。胸式呼吸についても同様であり、ここでは、胸式呼吸の大きさを示す第1バーグラフBG1の表示段階数数は「2」に設定されているので、被験者の胸式呼吸の大きさは「小程度」のレベルを満たしており、必須レベル(表示段階数「1」)を上回っていることが分かる。   In this case, if the number of display steps of the first bar graph BG1 indicating the magnitude of abdominal breathing is “1” or more, it means that the subject's abdominal breathing magnitude meets the required level, If the number of display steps is “2”, it means that the subject's abdominal breathing level satisfies the “small” level exceeding the required level, and if the display step number is “3”. Means that the “medium” level above is met, and if the number of display stages is “4”, it means that the “large” level is met, If the number of display steps is “5”, it means that the magnitude of the subject's abdominal breathing satisfies the “maximum” level. As described above, since the number of display steps of the first bar graph BG1 indicating the magnitude of abdominal breathing is set to “4”, the magnitude of the subject's abdominal breathing is “large”. It can be seen that there is a sufficient margin with respect to the required level (the number of display stages “1”). The same applies to thoracic breathing. Here, since the number of display steps of the first bar graph BG1 indicating the magnitude of thoracic breathing is set to “2”, the magnitude of thoracic breathing of the subject is It can be seen that the level of “small” is satisfied and exceeds the required level (number of display steps “1”).

なお、被験者の呼吸速度が速いほど、当該呼吸速度に応じた必須の呼吸レベルに対応する第1バーグラフBG1の段階数は増加するので、結果として、胸式呼吸および腹式呼吸の必須レベルに対応する段階数も増加する。例えば1呼吸に要する時間長が3秒以上であって、かつ4秒未満である場合は、必須の呼吸レベルに対応する第1バーグラフBG1の段階数は「4」に設定され、その段階数「4」は、胸式と腹式に2段階ずつ割り振られる。これに応じて、胸式呼吸および腹式呼吸の各々のマージンレベルが変化するという具合である。   As the subject's respiratory rate increases, the number of stages of the first bar graph BG1 corresponding to the essential respiratory level corresponding to the respiratory rate increases. As a result, the required level of chest and abdominal breathing is achieved. The corresponding number of stages also increases. For example, if the time length required for one breath is 3 seconds or more and less than 4 seconds, the number of stages of the first bar graph BG1 corresponding to the required breathing level is set to “4”, and the number of stages “4” is assigned in two stages for the chest type and the abdominal type. Correspondingly, the respective margin levels of chest and abdominal breathing change.

ステップS605の後、CPU170は、被験者の胸式呼吸および腹式呼吸の各々の大きさと余裕度とを第1バーグラフBG1で表示する一方、被験者の腹式レベルを第2バーグラフBG2で表示する(ステップS606)。図26に示すように、CPU170は、着色表示すべき第1バーグラフBG1の段階数(胸式2段階、腹式4段階)のうち腹式呼吸および胸式呼吸の各々の必須レベルに対応する段階数と、マージン分に対応する段階数とを異なる色で表示するので、被験者は、自分の腹式呼吸および胸式呼吸の各々の大きさと余裕度とを容易に把握することができる。   After step S605, the CPU 170 displays the magnitude and margin of each of the subject's chest breathing and abdominal breathing as a first bar graph BG1, while displaying the subject's abdominal level as a second bar graph BG2. (Step S606). As shown in FIG. 26, the CPU 170 corresponds to the essential levels of the abdominal breathing and the chest breathing among the number of stages of the first bar graph BG1 to be colored and displayed (two breast stages, four abdominal stages). Since the number of steps and the number of steps corresponding to the margin are displayed in different colors, the subject can easily grasp the size and margin of each of his abdominal breathing and chest breathing.

また、図27に示すように、第2バーグラフBG2で表示される腹式レベルの段階数は「5」に設定され、上述のステップS603で求められた腹式レベルの数値に応じて、5つの段階のうちの何れかが選択的に着色表示される。具体的には、CPU170は、腹式レベルの数値が0〜20の場合は1段階目のみを着色表示するように制御し、腹式レベルの数値が21〜40の場合は2段階目のみを着色表示するように制御し、腹式レベルの数値が41〜60の場合は3段階目のみを着色表示するように制御し、腹式レベルの数値が61〜80の場合は4段階目のみを着色表示するように制御し、腹式レベルの数値が81〜100の場合は5段階目のみを着色表示するように制御する。前述したように、ここでは、腹式レベルは「70」である場合を想定しているので、図27に示すように、CPU170は、第2バーグラフBG2の4段階目のみを着色表示するように制御するという具合である。   Further, as shown in FIG. 27, the number of steps of the abdominal level displayed in the second bar graph BG2 is set to “5”, and 5 is set according to the numerical value of the abdominal level determined in the above step S603. Any one of the two stages is selectively colored. Specifically, the CPU 170 performs control so that only the first step is colored when the abdominal level value is 0 to 20, and only the second step is displayed when the abdominal level value is 21 to 40. If the numerical value of the abdominal level is 41 to 60, control is performed so that only the third step is displayed in color. If the numerical value of the abdominal level is 61 to 80, only the fourth step is controlled. Control is performed so that the display is colored, and when the numerical value of the abdominal level is 81 to 100, control is performed so that only the fifth stage is displayed in color. As described above, since it is assumed here that the abdomen level is “70”, as shown in FIG. 27, the CPU 170 color-displays only the fourth stage of the second bar graph BG2. It is a condition to control to.

以上に説明したように、CPU170は、被験者の1呼吸ごとに、当該1呼吸における腹式呼吸および胸式呼吸の大きさと余裕度とを表示部160に表示するように制御するので、被験者は、自分の胸部呼吸筋および腹部呼吸筋の活用の強みと弱みとを認識できる。これにより、被験者に対して、自分の強みを自覚させつつも、自分の弱みを活性化させるための腹式呼吸等の呼吸筋トレーニングへのモチベーションを確保させることができる。また、本実施形態によれば、スパイロなどのように、被験者が最大の呼吸を行わなくとも、当該被験者の呼吸能力の余裕を把握できるので、被験者の安全を確保するという観点からも好ましいという利点がある。   As described above, the CPU 170 controls the display unit 160 to display the magnitude and margin of abdominal breathing and chest breathing in one breath for each breath of the subject. Recognize the strengths and weaknesses of using your chest and abdominal respiratory muscles. This makes it possible to ensure motivation for respiratory muscle training such as abdominal breathing for activating one's weakness while making the subject aware of his strength. In addition, according to the present embodiment, it is preferable from the viewpoint of ensuring the safety of the subject because the subject can grasp the margin of the breathing ability of the subject without performing the maximum breathing, such as a spiro. There is.

<1−6:換気特性判定/表示処理>
次に、CPU170が実行する換気特性判定/表示処理について説明する。本実施形態では、CPU170は、被験者の1呼吸ごとに、当該1呼吸における換気機能の特性(換気特性)を判定するとともに、その判定結果を表示(報知)する換気特性判定/表示処理を実行する。以下、図30を参照しながら、CPU170が実行する換気特性判定/表示処理の具体的な内容を説明する。図30は、換気特性判定/表示処理の具体的な内容を示すフローチャートである。
<1-6: Ventilation characteristic determination / display processing>
Next, ventilation characteristic determination / display processing executed by the CPU 170 will be described. In the present embodiment, the CPU 170 determines a ventilation function characteristic (ventilation characteristic) for each breath of the subject and executes a ventilation characteristic determination / display process for displaying (notifying) the determination result. . Hereinafter, specific contents of the ventilation characteristic determination / display process executed by the CPU 170 will be described with reference to FIG. FIG. 30 is a flowchart showing specific contents of the ventilation characteristic determination / display process.

図30に示すように、まずCPU170は、直前の1呼吸における被験者の呼吸レベルが、安静時必須レベル以上であるか否かを判断する(ステップS700)。より具体的には、前述の呼吸レベル表示処理で説明した内容(図28のステップS601、ステップS602)と同様に、CPU170は、直前の1呼吸における呼吸レベルの正規化値%ΔZaの値に対応する1回換気量の正規化値%ΔTVを求める。そして、その求めた1回換気量の正規化値%ΔTVに基づいて、被験者の呼吸レベルが、安静時必須レベルを上回っているか否かを判断する。CPU170は、1回換気量の正規化値%ΔTVの値が第4所定値α4以上である場合は、被験者の呼吸レベルは、安静時における必須レベル以上であると判断するという具合である(図29参照)。   As shown in FIG. 30, first, the CPU 170 determines whether or not the breathing level of the subject in the immediately preceding one breath is equal to or higher than the essential level at rest (step S700). More specifically, similar to the content described in the above-described respiration level display process (steps S601 and S602 in FIG. 28), the CPU 170 corresponds to the value of the normalized value% ΔZa of the respiration level in the immediately preceding one respiration. Obtain the normalized value% ΔTV of the tidal volume. Then, based on the obtained normalized value% ΔTV of the tidal volume, it is determined whether or not the breathing level of the subject exceeds the essential level at rest. When the normalized value% ΔTV of the tidal volume is equal to or greater than the fourth predetermined value α4, the CPU 170 determines that the subject's breathing level is equal to or higher than the essential level at rest (FIG. 29).

ステップS700の結果が否定である場合、CPU170は、被験者の呼吸レベルが安静時必須レベル以上となるように、指導情報を報知する(ステップS701)。呼吸を指導するためのメッセージが表示部160に表示される態様であってもよいし、音声による報知が行われる態様であってもよい。また、両者を組み合わせてもよい。   If the result of step S700 is negative, the CPU 170 notifies the guidance information so that the subject's breathing level is equal to or higher than the essential level at rest (step S701). A mode in which a message for instructing breathing may be displayed on the display unit 160, or a mode in which voice notification is performed may be used. Moreover, you may combine both.

一方、ステップS700の結果が肯定である場合、あるいは上述のステップS701の後、CPU170は、直前の1呼吸における被験者の換気機能の特性を判定する(ステップS702)。より具体的には、CPU170は、当該1呼吸における第1生体電気インピーダンスZaの測定値および第2生体電気インピーダンスZbの測定値に基づいて、当該1呼吸における被験者の換気機能の特性を判定するという具合である。以下、具体的な内容について説明する。   On the other hand, if the result of step S700 is affirmative or after step S701 described above, the CPU 170 determines the characteristics of the subject's ventilation function in the immediately preceding one breath (step S702). More specifically, the CPU 170 determines the characteristics of the ventilation function of the subject in the one breath based on the measured value of the first bioelectric impedance Za and the measured value of the second bioelectric impedance Zb in the one breath. Condition. Hereinafter, specific contents will be described.

例えば被験者が胸部に病歴を有しており、胸部のうち呼吸に寄与する部位(例えば胸部骨格筋、肺胞組織、胸部骨格筋および肺胞組織の両方)の機能が低下していた場合は、呼吸に伴う胸部骨格筋の変化が健常者に比べて小さい、または不足するので、換気量を確保するために、横隔膜の変動量(上昇量および下降量)が大きくなる。また、加齢などにより、胸部のうち呼吸に寄与する部位の機能が低下していた場合も同様である。ここで、肺の上部を含み、腹部を含まない体幹上部の第1生体電気インピーダンスZaの変化は、肺に出入りする空気量と連動している一方、肺の中下部および腹部を含む第2生体電気インピーダンスZbの変化は、横隔膜の動きと連動しており、横隔膜の変動量が大きいほど第2生体電気インピーダンスZbの変化も大きくなる。したがって、例えば胸部のうち呼吸に寄与する部位の機能が低下していた被験者と健常者との間で、1呼吸において肺に出入りする空気の量(換気量)が同じであり、当該1呼吸における第1生体電気インピーダンスZaの変化が同じであっても、胸部のうち呼吸に寄与する部位の機能が低下していた被験者は、健常者に比べて横隔膜の変動量が大きいので、当該1呼吸における第2生体電気インピーダンスZbの変化が大きくなる。   For example, if the subject has a medical history in the chest and the function of the chest part that contributes to breathing (for example, thoracic skeletal muscle, alveolar tissue, both thoracic skeletal muscle and alveolar tissue) Since the change of the thoracic skeletal muscle accompanying respiration is smaller or insufficient than that of a healthy person, the amount of diaphragm fluctuation (increase and decrease) is increased in order to ensure ventilation. The same applies to the case where the function of the part of the chest that contributes to breathing is reduced due to aging or the like. Here, the change in the first bioelectrical impedance Za including the upper part of the lung and not including the abdomen is linked to the amount of air entering and exiting the lung, while the second including the middle lower part of the lung and the abdomen. The change in the bioelectrical impedance Zb is linked to the movement of the diaphragm, and the change in the second bioelectrical impedance Zb increases as the amount of fluctuation of the diaphragm increases. Therefore, for example, the amount of air that enters and exits the lungs in one breath (ventilation amount) is the same between the subject who has reduced the function of the part of the chest that contributes to breathing and the healthy person, Even if the change of the first bioelectrical impedance Za is the same, the subject whose function of the portion of the chest that contributes to respiration has decreased has a larger amount of fluctuation of the diaphragm than that of a healthy person. The change in the second bioelectrical impedance Zb becomes large.

図31は、被験者の呼吸に伴う第1相対値ΔZaおよび第2相対値ΔZbの経時的変化を示す図であって、この場合の被験者は、胸部に病歴の無い男性(健常者)である。一方、図32は、被験者の呼吸に伴う第1相対値ΔZaおよび第2相対値ΔZbの経時的変化を示す図であって、この場合の被験者は、胸部に病歴を有する男性である。図31および図32からも理解されるように、被験者の胸部のうち呼吸に寄与する部位の機能が低下している場合は、第2相対値ΔZbの振幅値は大きくなることが分かる。   FIG. 31 is a diagram showing temporal changes in the first relative value ΔZa and the second relative value ΔZb associated with the subject's breathing. In this case, the subject is a male (healthy person) who has no medical history in the chest. On the other hand, FIG. 32 is a diagram showing temporal changes in the first relative value ΔZa and the second relative value ΔZb accompanying the breathing of the subject, and the subject in this case is a man who has a medical history in the chest. As can be understood from FIGS. 31 and 32, it can be seen that the amplitude value of the second relative value ΔZb increases when the function of the portion of the subject's chest that contributes to breathing is reduced.

本実施形態では、被験者の胸部のうち呼吸に寄与する部位の機能が低下している場合は、肺の中下部および腹部を含む体幹中部の第2生体電気インピーダンスZbの変化が大きくなる点に着目し、第2生体電気インピーダンスZbの測定値と第1生体電気インピーダンスZaの測定値とに基づいて、被験者の胸部のうち呼吸に寄与する部位の機能が正常であるか否かを判定している(換気機能の特性を判定している)。ここで、前述したように、腹式呼吸の呼気における第2生体電気インピーダンスZbの変化は、第1生体電気インピーダンスZaの変化とは異なる(変化を示す波形が異なる)のに対して、吸気における第2生体電気インピーダンスZbの変化は、胸式呼吸および腹式呼吸の何れの場合であっても、第1生体電気インピーダンスZaの変化と同じである(変化を示す波形が同じである)。すなわち、被験者の呼吸に伴う第2生体電気インピーダンスZbの経時的変化を示す波形には、呼気時の腹式呼吸に起因した波形歪みが含まれるので、吸気における第1生体電気インピーダンスZaの変化と第2生体電気インピーダンスZbの変化とに基づいて、被験者の胸部のうち呼吸に寄与する部位の機能が正常であるか否かを判定することが好ましい。   In the present embodiment, when the function of a part contributing to breathing in the subject's chest is reduced, the change in the second bioelectric impedance Zb in the middle trunk including the middle and lower abdomen of the lung is increased. Pay attention and determine whether the function of the part contributing to breathing in the chest of the subject is normal based on the measured value of the second bioelectrical impedance Zb and the measured value of the first bioelectrical impedance Za Yes (judging the characteristics of ventilation function). Here, as described above, the change in the second bioelectrical impedance Zb in the expiration of abdominal breathing is different from the change in the first bioelectrical impedance Za (the waveform indicating the change is different), whereas in the inspiration The change in the second bioelectrical impedance Zb is the same as the change in the first bioelectrical impedance Za in both cases of chest respiration and abdominal respiration (the waveforms indicating the change are the same). That is, since the waveform showing the change over time of the second bioelectric impedance Zb accompanying the breathing of the subject includes waveform distortion caused by abdominal breathing during expiration, the change in the first bioelectric impedance Za during inspiration It is preferable to determine whether or not the function of the part contributing to breathing in the subject's chest is normal based on the change in the second bioelectrical impedance Zb.

そこで、本実施形態では、吸気における第1相対値ΔZaおよび第2相対値ΔZbに基づいて、被験者の胸部のうち呼吸に寄与する部位の機能が正常であるか否かを判定している。より具体的には以下のとおりである。ステップS703において、CPU170は、被験者の直前の1呼吸の吸気における第1相対値ΔZaのピーク値ΔZa(MAX)、および、第2相対値ΔZbのピーク値ΔZb(MAX)を抽出する。そして、CPU170は、第1相対値ΔZaのピーク値ΔZa(MAX)に対する第2相対値ΔZbのピーク値ΔZb(MAX)の割合(=ΔZb(MAX)/ΔZa(MAX))を、被験者の換気機能の特性を示す換気能バランスBPとして算出する。CPU170は、その算出した換気能バランスBPの値が所定の閾値以上である場合は、被験者の胸部のうち呼吸に寄与する部位の機能は正常ではないと判定し、所定の閾値を下回る場合は、被験者の胸部のうち呼吸に寄与する部位の機能は正常であると判定するという具合である。   Therefore, in the present embodiment, it is determined based on the first relative value ΔZa and the second relative value ΔZb in inspiration whether or not the function of the part contributing to breathing in the subject's chest is normal. More specifically, it is as follows. In step S703, the CPU 170 extracts the peak value ΔZa (MAX) of the first relative value ΔZa and the peak value ΔZb (MAX) of the second relative value ΔZb in the inspiration of one breath immediately before the subject. Then, the CPU 170 determines the ratio of the peak value ΔZb (MAX) of the second relative value ΔZb to the peak value ΔZa (MAX) of the first relative value ΔZa (= ΔZb (MAX) / ΔZa (MAX)), and the ventilation function of the subject. Calculated as a ventilation capacity balance BP indicating the characteristics of When the calculated value of the ventilation capacity balance BP is equal to or greater than a predetermined threshold, the CPU 170 determines that the function of the part that contributes to breathing in the subject's chest is not normal, and if below the predetermined threshold, In other words, the function of the part contributing to breathing in the subject's chest is determined to be normal.

また、女性は胸式優位の呼吸を行い、男性は腹式優位の呼吸を行うのが一般的である。すなわち、換気量が同じであっても、呼吸に伴う横隔膜の変動量は男性の方が大きいので、上述の換気能バランスBPの値は、男性の方が高い値を示す。また、内臓脂肪の蓄積量が多い被験者においては、腹部骨格筋の活動が制限される。女性は、妊娠などによって腹部骨格筋の活動が制限されても、それを補って換気能力を確保するための代償能力を有しているが、男性は、そのような代償能力を有していないので、内臓脂肪の蓄積量が多い場合はクロージングボリュームが低下し、上述の換気能バランスBPの値は高い値を示す。このように、上述の換気能バランスBPから、被験者の換気機能の特性を判定できるという具合である。   Moreover, it is common for a woman to take a breath of the chest type and a man to take a breath of the abdominal type. That is, even if the ventilation volume is the same, the amount of fluctuation of the diaphragm accompanying breathing is larger in men, so the value of the above-described ventilation capacity balance BP shows a higher value in men. In addition, the activity of abdominal skeletal muscle is restricted in subjects with a large amount of visceral fat accumulation. Women have the compensatory ability to compensate for the limited activity of abdominal skeletal muscles, such as pregnancy, while men have no such compensatory ability. Therefore, when the accumulation amount of visceral fat is large, the closing volume is reduced, and the value of the ventilation capacity balance BP is high. Thus, the characteristic of the subject's ventilation function can be determined from the above-described ventilation capacity balance BP.

図30に示すように、ステップS702の後、CPU170は判定結果を表示部160に表示するように制御する(ステップS703)。本実施形態では、図33に示すように、ステップS702の判定結果は第3バーグラフBG3で表示される。図33に示すように、第3バーグラフBG3で表示される換気能バランスBPの段階数は「5」に設定され、ステップS702で得られた換気能バランスBPの値に応じて、5つの段階のうちの何れかが選択的に着色表示される。具体的には、CPU170は、換気能バランスBPの値が0〜0.4の場合は1段階目のみを着色表示するように制御し、換気能バランスBPの値が0.4〜0.8の場合は2段階目のみを着色表示するように制御し、換気能バランスBPの値が0.8〜1.2の場合は3段階目のみを着色表示するように制御し、換気能バランスBPの値が1.2〜1.6の場合は4段階目のみを着色表示するように制御し、換気能バランスBPの値が1.6を超える場合は5段階目のみを着色表示するように制御する。例えば換気能バランスBPの値が0.5の場合は、図33に示すように、CPU170は、第3バーグラフBG3の2段階目のみを着色表示するように制御する。   As shown in FIG. 30, after step S702, the CPU 170 controls to display the determination result on the display unit 160 (step S703). In the present embodiment, as shown in FIG. 33, the determination result in step S702 is displayed as a third bar graph BG3. As shown in FIG. 33, the number of steps of the ventilation capacity balance BP displayed in the third bar graph BG3 is set to “5”, and five stages are set according to the value of the ventilation capacity balance BP obtained in step S702. Any one of these is selectively colored and displayed. Specifically, when the value of the ventilation capacity balance BP is 0 to 0.4, the CPU 170 performs control so that only the first stage is colored, and when the value of the ventilation capacity balance BP is 0.4 to 0.8, the second stage. Only when the value of ventilation capacity balance BP is 0.8 to 1.2, control is performed so that only the third stage is colored and displayed, and when the value of ventilation capacity balance BP is 1.2 to 1.6, 4 is controlled. Control is performed so that only the stage is colored, and when the value of the ventilation capacity balance BP exceeds 1.6, control is performed so that only the fifth stage is colored. For example, when the value of the ventilation capacity balance BP is 0.5, as shown in FIG. 33, the CPU 170 performs control so that only the second stage of the third bar graph BG3 is displayed in color.

例えば、被験者の呼吸が胸式優位(典型的な女性の呼吸)である場合は、第3バーグラフBG3の1段階目が着色表示され、被験者の呼吸が腹式優位(典型的な男性の呼吸)である場合は、第3バーグラフBG3の2段階目が着色表示され、被験者の換気機能が低下している場合(胸部のうち呼吸に寄与する部位の機能が正常でない場合、内蔵脂肪蓄積量が過大な場合等)は、第3バーグラフBG3の3段階目以降が着色表示されるという具合である。被験者は、第3バーグラフBG3を見ることで、自分の換気機能の特性(例えば胸部のうち呼吸に寄与する部位の機能が低下している等)を正確かつ容易に認識できるという具合である。以上より、本実施形態によれば、被験者の換気機能の特性を正確かつ容易に判定できる。   For example, if the subject's breathing is chest dominant (typical female breathing), the first stage of the third bar graph BG3 is colored and the subject's breathing is abdominal dominant (typical male breathing). ), The second stage of the third bar graph BG3 is colored and displayed, and the subject's ventilation function is reduced (if the function of the part of the chest that contributes to breathing is not normal, the amount of built-in fat accumulation Is excessive), the third and subsequent stages of the third bar graph BG3 are colored and displayed. The test subject can recognize the characteristic of his / her ventilation function (for example, the function of a part contributing to breathing in the chest is reduced) accurately and easily by looking at the third bar graph BG3. As mentioned above, according to this embodiment, the characteristic of a subject's ventilation function can be determined correctly and easily.

<2.変形例>
本発明は上述した実施形態に限定されるものではなく、例えば、以下の変形が可能である。また、以下に示す変形例のうちの2以上の変形例を組み合わせることもできる。
<2. Modification>
The present invention is not limited to the above-described embodiments, and for example, the following modifications are possible. Also, two or more of the modifications shown below can be combined.

(1)変形例1
上述の実施形態では、被験者の肺全体の換気機能の特性を示すバランス指標BPを求めて表示しているが、これに限らず、例えば右肺の換気機能の特性を示す指標(右換気能バランスBPR)と左肺の換気機能の特性を示す指標(左換気能バランスBPL)とを独立に求めて表示する態様であってもよい。
(1) Modification 1
In the above-described embodiment, the balance index BP indicating the characteristics of the ventilation function of the entire lungs of the subject is obtained and displayed. However, the present invention is not limited thereto, and for example, an index indicating the characteristics of the ventilation function of the right lung (right ventilation capacity balance) BP R ) and an index indicating the characteristics of the ventilation function of the left lung (left ventilation capacity balance BP L ) may be independently obtained and displayed.

右換気能バランスBPRは、右肺の上部を含み、腹部を含まない体幹上部右側の第3生体電気インピーダンスZaRの測定値と、右肺の中下部および腹部を含む体幹中部右側の第4生体電気インピーダンスZbRの測定値とに基づいて求められる。より具体的には以下のとおりである。
第3生体電気インピーダンスZaRは、右手と左手との間に流れる基準電流Irefの大きさを示す電流データDiと、右足用の電圧電極Y2と右手用の電圧電極Y4との間の電位差を示す電圧データDvとから求められる。なお、被験者の呼吸に伴う第3生体電気インピーダンスZaRの変化の仕方は、第1生体電気インピーダンスZaの変化の仕方と同様である。
第4生体電気インピーダンスZbRは、右手と右足との間に流れる基準電流Irefの大きさを示す電流データDiと、左手と左足との間の電位差を示す電圧データDvとから求められる。なお、呼吸に伴う第4生体電気インピーダンスZbRの変化の仕方は、第2生体電気インピーダンスZbの変化の仕方と同様である。
また、第3生体電気インピーダンスZaRの変化を示す波形から、呼吸に起因した情報を抽出するための振幅基準レベルを第3センタリング値ZaR0と呼ぶ。第3センタリング値ZaR0の生成方法は、上述の第1センタリング値Za0の生成方法と同様である。また、第4生体電気インピーダンスZbRの変化を示す波形から、呼吸に起因した情報を抽出するための振幅基準レベルを第4センタリング値ZbR0と呼ぶ。第4センタリング値ZbR0の生成方法は、上述の第2センタリング値Zb0の生成方法と同様である。
さらに、第3生体電気インピーダンスZaRの測定値の第3センタリング値ZaR0に対する相対値を第3相対値ΔZaR、第4生体電気インピーダンスZbRの第4センタリング値ZbR0に対する相対値を第4相対値ΔZbRとする。CPU170は、被験者の1呼吸ごとに、当該1呼吸の吸気における第3相対値のピーク値ΔZaR(MAX)に対する第4相対値のピーク値ΔZbR(MAX)の割合(=ΔZbR(MAX)/ΔZaR(MAX))を右換気能バランスBPRとして算出し、その算出結果を表示するように表示部160を制御するという具合である。
Right ventilation capacity balance BP R includes the upper right lung, and the measurement value of the third bioelectrical impedance ZaR the upper trunk right without the abdomen, the middle trunk right including the lower and middle abdomen of the right lung the 4 It is calculated | required based on the measured value of bioelectrical impedance ZbR. More specifically, it is as follows.
The third bioelectrical impedance ZaR is the current data Di indicating the magnitude of the reference current Iref flowing between the right hand and the left hand, and the voltage indicating the potential difference between the right foot voltage electrode Y2 and the right hand voltage electrode Y4. It is obtained from the data Dv. The method of changing the third bioelectrical impedance ZaR accompanying the breathing of the subject is the same as the method of changing the first bioelectrical impedance Za.
The fourth bioelectrical impedance ZbR is obtained from current data Di indicating the magnitude of the reference current Iref flowing between the right hand and the right foot, and voltage data Dv indicating the potential difference between the left hand and the left foot. Note that the manner of changing the fourth bioelectrical impedance ZbR accompanying respiration is the same as the way of changing the second bioelectrical impedance Zb.
An amplitude reference level for extracting information resulting from respiration from a waveform indicating a change in the third bioelectrical impedance ZaR is referred to as a third centering value ZaR0. The method for generating the third centering value ZaR0 is the same as the method for generating the first centering value Za0. An amplitude reference level for extracting information resulting from respiration from a waveform indicating a change in the fourth bioelectrical impedance ZbR is referred to as a fourth centering value ZbR0. The method for generating the fourth centering value ZbR0 is the same as the method for generating the second centering value Zb0.
Further, the relative value of the measured value of the third bioelectrical impedance ZaR with respect to the third centering value ZaR0 is the third relative value ΔZaR, and the relative value of the fourth bioelectrical impedance ZbR with respect to the fourth centering value ZbR0 is the fourth relative value ΔZbR. . The CPU 170 takes the ratio of the peak value ΔZbR (MAX) of the fourth relative value to the peak value ΔZaR (MAX) of the third relative value in the inspiration of the one breath for each breath of the subject (= ΔZbR (MAX) / ΔZaR ( MAX)) is calculated as a right ventilation capacity balance BP R, and so that controls the display unit 160 to display the calculation result.

一方、左換気能バランスBPLは、左肺の上部を含み、腹部を含まない体幹上部左側の第5生体電気インピーダンスZaLの測定値と、左肺の中下部および腹部を含む体幹中部左側の第6生体電気インピーダンスZbLの測定値とに基づいて求められる。より具体的には以下のとおりである。
第5生体電気インピーダンスZaLは、右手と左手との間に流れる基準電流Irefの大きさを示す電流データDiと、左足用の電圧電極Y1と左手用の電圧電極Y3との間の電位差を示す電圧データDvとから求められる。なお、被験者の呼吸に伴う第5生体電気インピーダンスZaLの変化の仕方は、第1生体電気インピーダンスZaの変化の仕方と同様である。
第6生体電気インピーダンスZbLは、左手と左足との間に流れる基準電流Irefの大きさを示す電流データDiと、右手と右足との間の電位差を示す電圧データDvとから求められる。なお、呼吸に伴う第6生体電気インピーダンスZbLの変化の仕方は、第2生体電気インピーダンスZbの変化の仕方と同様である。
また、第5生体電気インピーダンスZaLの変化を示す波形から、呼吸に起因した情報を抽出するための振幅基準レベルを第5センタリング値ZaL0と呼ぶ。第5センタリング値ZaL0の生成方法は、上述の第1センタリング値Za0の生成方法と同様である。また、第6生体電気インピーダンスZbLの変化を示す波形から、呼吸に起因した情報を抽出するための振幅基準レベルを第6センタリング値ZbL0と呼ぶ。第6センタリング値ZbL0の生成方法は、上述の第2センタリング値Zb0の生成方法と同様である。
さらに、第5生体電気インピーダンスZaLの測定値の第5センタリング値ZaL0に対する相対値を第5相対値ΔZaL、第6生体電気インピーダンスZbLの測定値の第6センタリング値ZbL0に対する相対値を第6相対値ΔZbLとする。CPU170は、被験者の1呼吸ごとに、当該1呼吸の吸気における第5相対値のピーク値ΔZaL(MAX)に対する第6相対値のピーク値ΔZbL(MAX)の割合(=ΔZbL(MAX)/ΔZaL(MAX))を左換気能バランスBPLとして算出し、その算出結果を表示するように表示部160を制御するという具合である。以上の態様によれば、右肺の換気特性と左肺の換気特性とを別々に判定できる。
On the other hand, the left ventilation capacity balance BP L includes the upper left lung, the measured value of the fifth bioelectrical impedance ZaL the upper trunk left without the abdomen, middle trunk left including the lower and middle abdomen of the left lung And the measured value of the sixth bioelectric impedance ZbL. More specifically, it is as follows.
The fifth bioelectrical impedance ZaL is current data Di that indicates the magnitude of the reference current Iref that flows between the right hand and the left hand, and a voltage that indicates a potential difference between the voltage electrode Y1 for the left foot and the voltage electrode Y3 for the left hand. It is obtained from the data Dv. Note that the way of changing the fifth bioelectrical impedance ZaL accompanying the breathing of the subject is the same as the way of changing the first bioelectrical impedance Za.
The sixth bioelectrical impedance ZbL is obtained from current data Di indicating the magnitude of the reference current Iref flowing between the left hand and the left foot and voltage data Dv indicating the potential difference between the right hand and the right foot. Note that the manner of changing the sixth bioelectrical impedance ZbL accompanying respiration is the same as the manner of changing the second bioelectrical impedance Zb.
An amplitude reference level for extracting information resulting from respiration from a waveform indicating a change in the fifth bioelectrical impedance ZaL is referred to as a fifth centering value ZaL0. The method for generating the fifth centering value ZaL0 is the same as the method for generating the first centering value Za0. An amplitude reference level for extracting information resulting from respiration from a waveform indicating a change in the sixth bioelectrical impedance ZbL is referred to as a sixth centering value ZbL0. The method for generating the sixth centering value ZbL0 is the same as the method for generating the second centering value Zb0.
Furthermore, the relative value of the measured value of the fifth bioelectrical impedance ZaL with respect to the fifth centering value ZaL0 is the fifth relative value ΔZaL, and the relative value of the measured value of the sixth bioelectrical impedance ZbL with respect to the sixth centering value ZbL0 is the sixth relative value. Let ΔZbL. The CPU 170 performs, for each breath of the subject, the ratio of the peak value ΔZbL (MAX) of the sixth relative value to the peak value ΔZaL (MAX) of the fifth relative value in the inspiration of the one breath (= ΔZbL (MAX) / ΔZaL ( MAX)) is calculated as the left ventilation capacity balance BP L, and so that controls the display unit 160 to display the calculation result. According to the above aspect, the ventilation characteristics of the right lung and the ventilation characteristics of the left lung can be determined separately.

(2)変形例2
上述の実施形態では、被験者の1呼吸の吸気における第1相対値のピーク値ΔZa(MAX)に対する第2相対値ΔZb(MAX)の割合(=ΔZb(MAX)/ΔZa(MAX))を、当該1呼吸における被験者の換気機能の特性を示す換気能バランスBPとして採用しているが、これに限らず、例えば第2相対値のピーク値ΔZb(MAX)に対する第1相対値のピーク値ΔZa(MAX)の割合(=ΔZa(MAX)/ΔZb(MAX))を換気能バランスBPとして採用することもできる。また、例えば被験者の1呼吸の吸気における第1相対値ΔZaの積分値に対する第2相対値ΔZbの積分値の割合(=第2相対値ΔZbの積分値/第1相対値ΔZaの積分値)を換気能バランスBPとして採用することもできるし、1呼吸の吸気における第2相対値ΔZbの積分値に対する第1相対値ΔZaの積分値の割合(=第1相対値ΔZaの積分値/第2相対値ΔZbの積分値)を換気能バランスBPとして採用することもできる。要するに、1呼吸の吸気における第1相対値ΔZaと第2相対値ΔZbとの比を換気能バランスBPとして採用することができる。
(2) Modification 2
In the above embodiment, the ratio (= ΔZb (MAX) / ΔZa (MAX)) of the second relative value ΔZb (MAX) to the peak value ΔZa (MAX) of the first relative value in inhalation of one breath of the subject is calculated as Although employed as a ventilation capacity balance BP indicating the characteristics of the ventilation function of the subject in one breath, the present invention is not limited to this. For example, the peak value ΔZa (MAX of the first relative value relative to the peak value ΔZb (MAX) of the second relative value is used. ) Ratio (= ΔZa (MAX) / ΔZb (MAX)) can also be employed as the ventilation capacity balance BP. Further, for example, the ratio of the integrated value of the second relative value ΔZb to the integrated value of the first relative value ΔZa in one breath of the subject (= the integrated value of the second relative value ΔZb / the integrated value of the first relative value ΔZa). It can also be adopted as a ventilation capacity balance BP, and the ratio of the integrated value of the first relative value ΔZa to the integrated value of the second relative value ΔZb in inspiration of one breath (= the integrated value of the first relative value ΔZa / the second relative value). The integrated value of the value ΔZb) can also be adopted as the ventilation capacity balance BP. In short, the ratio between the first relative value ΔZa and the second relative value ΔZb in inhalation of one breath can be adopted as the ventilation capacity balance BP.

(3)変形例3
上述の実施形態では、肺の上部の換気機能を示す指標(ΔZa(MAX))と肺の中下部の換気機能を示す指標(ΔZb(MAX))との比を換気能バランスBPとして採用しているが、これに限らず、例えば肺の上部と中下部とを合わせた肺全体の換気機能を示す指標(ΔZa(MAX)+ΔZb(MAX))と肺の上部の換気機能を示す指標(ΔZa(MAX))との比、あるいは、肺全体の換気機能を示す指標(ΔZa(MAX)+ΔZb(MAX))と肺の中下部の換気機能を示す指標(ΔZb(MAX))との比を換気能バランスBPとして採用することもできる。
(3) Modification 3
In the above-described embodiment, the ratio of the index indicating the ventilation function of the upper part of the lung (ΔZa (MAX)) and the index indicating the ventilation function of the lower part of the lung (ΔZb (MAX)) is adopted as the ventilation capacity balance BP. However, the present invention is not limited to this. For example, an index indicating the ventilation function of the entire lung (ΔZa (MAX) + ΔZb (MAX)) combined with the upper and middle parts of the lung and an index indicating the ventilation function of the upper part of the lung (ΔZa ( MAX)), or the ratio of the index indicating the ventilation function of the entire lung (ΔZa (MAX) + ΔZb (MAX)) and the index indicating the ventilation function of the middle and lower part of the lung (ΔZb (MAX)). It can also be adopted as a balance BP.

(4)変形例4
上述の実施形態において、電流電極および電圧電極の一例として、両手両足を電極の接点とする四肢誘導八電極法を一例として説明したが、本発明はこれに限定されるものではない。例えば、耳電極との四肢誘導法とを組み合わせて、体幹上部の生体電気インピーダンスを測定してもよい。
この場合には、耳電極を用いることによって、体幹上部の生体電気インピーダンスの測定について両腕計測ではなく片腕計測が可能となる。なお、耳電極を用いる場合には、イヤホンやヘッドホンに耳電極を組み込むことによって、音声等の音報知・音刺激との組み合わせが効果的である。
(4) Modification 4
In the above-described embodiment, as an example of the current electrode and the voltage electrode, the limb induction eight electrode method using both hands and both feet as electrode contacts has been described as an example, but the present invention is not limited thereto. For example, the bioelectrical impedance of the upper trunk may be measured in combination with the limb induction method with the ear electrode.
In this case, by using the ear electrode, the measurement of the bioelectrical impedance of the upper trunk can be performed by one-arm measurement rather than by both-arm measurement. In the case where an ear electrode is used, a combination with sound notification / stimulation such as voice is effective by incorporating the ear electrode into an earphone or a headphone.

また、上述した実施形態では、立位での計測であったが、本発明はこれに限定されるものではなく、便座での生体電気インピーダンスを計測するようにしてもよい。この場合には、便座や手摺に電極確保することができる。さらに、ポケッタブルやウェアラブルでのリラクゼーション姿勢(椅子座位等)で生体電気インピーダンスを計測するようにしてもよい。この場合には、マッサージチェアー等の手摺と足置き等に電極確保することができる。   In the above-described embodiment, the measurement is performed in the standing position. However, the present invention is not limited to this, and the bioelectric impedance in the toilet seat may be measured. In this case, electrodes can be secured on the toilet seat and handrail. Furthermore, the bioelectrical impedance may be measured in a pocketable or wearable relaxation posture (chair sitting position, etc.). In this case, an electrode can be secured for a handrail such as a massage chair and a footrest.

さらに、入浴中の呼吸計測も可能である。この場合には、浴槽手摺部と浴槽底の御尻や足裏接触側面部に電極を設ける。浴槽内のお湯よりも、体幹の方が、生理食塩水でできているので電流通電が支配的になる。よって、入浴中にリラックスした状態で呼吸法のトレーニングを行うことができる。
くわえて、上述した実施形態の生体測定装置1に、血圧計の腕帯と手で握る等で接触させる血圧計を付加し、血圧計に電極配置して呼吸変化や腕の筋の緊張具合を血圧測定時の補正情報として活用してもよい。
Furthermore, respiration measurement during bathing is also possible. In this case, electrodes are provided on the bathtub handrail part, the bottom of the bathtub bottom, and the sole contact side part. Current flow becomes dominant because the trunk is made of physiological saline rather than hot water in the bathtub. Therefore, breathing training can be performed in a relaxed state during bathing.
In addition, a sphygmomanometer that is brought into contact with the sphygmomanometer arm band by hand or the like is added to the biometric device 1 of the above-described embodiment, and electrodes are arranged on the sphygmomanometer to control respiratory changes and tension of arm muscles. You may utilize as correction information at the time of blood-pressure measurement.

(5)変形例5
上述の実施形態では、CPU170は、右手と左手との間に流れる基準電流Irefの大きさを示す電流データDiと、右足用の電圧電極Y2と右手用の電圧電極Y4との間の電位差を示す電圧データDvとから、右上肢の生体電気インピーダンスを測定しているが、これに限らず、肺の上部を含み、腹部を含まない体幹上部の第1生体電気インピーダンスZaの測定方法は任意である。例えば右手と左手との間に流れる基準電流Irefの大きさを示す電流データDiと、左足用の電圧電極Y1と左手用の電圧電極Y3との間の電位差を示す電圧データDvとから、左上肢の生体電気インピーダンスを測定し、それを第1生体電気インピーダンスZaとして採用してもよい。
(5) Modification 5
In the above-described embodiment, the CPU 170 indicates the current data Di indicating the magnitude of the reference current Iref flowing between the right hand and the left hand, and the potential difference between the voltage electrode Y2 for the right foot and the voltage electrode Y4 for the right hand. The bioelectrical impedance of the upper right limb is measured from the voltage data Dv. However, the measurement method of the first bioelectrical impedance Za of the upper trunk including the upper part of the lung and not including the abdomen is arbitrary. is there. For example, from the current data Di indicating the magnitude of the reference current Iref flowing between the right hand and the left hand and the voltage data Dv indicating the potential difference between the voltage electrode Y1 for the left foot and the voltage electrode Y3 for the left hand, Alternatively, the bioelectric impedance may be measured and used as the first bioelectric impedance Za.

また、上述の実施形態では、CPU170は、左足と右手との間に流れる基準電流Irefの大きさを示す電流データDiと、左手用の電圧電極Y3と右足用の電圧電極Y2との間の電位差を示す電圧データDvとから、体幹中部の第2生体電気インピーダンスZbを測定しているが、これに限らず、肺の中下部および腹部を含む体幹中部の第2生体電気インピーダンスZbの測定方法は任意である。例えば右足と左手との間に流れる基準電流Irefの大きさを示す電流データDiと、右手用の電圧電極Y4と左足用の電圧電極Y1との間で電位差を示す電圧データDvとから、体幹中部を斜めに横切る生体電気インピーダンスを測定し、それを第2生体電気インピーダンスZbとして採用してもよい。   In the embodiment described above, the CPU 170 determines the potential difference between the current data Di indicating the magnitude of the reference current Iref flowing between the left foot and the right hand, and the voltage electrode Y3 for the left hand and the voltage electrode Y2 for the right foot. The second bioelectrical impedance Zb in the middle trunk is measured from the voltage data Dv indicating, but not limited to this, the second bioelectrical impedance Zb in the middle trunk including the middle lower part and the abdomen of the lung is measured. The method is arbitrary. For example, from the current data Di indicating the magnitude of the reference current Iref flowing between the right foot and the left hand and the voltage data Dv indicating the potential difference between the voltage electrode Y4 for the right hand and the voltage electrode Y1 for the left foot, the trunk The bioelectrical impedance that crosses the middle part diagonally may be measured and used as the second bioelectrical impedance Zb.

1 生体測定装置
120 第1記憶部
150 入力部
170 CPU
200 生体電気インピーダンス測定部
1 Biological Measuring Device 120 First Storage Unit 150 Input Unit 170 CPU
200 Bioelectrical impedance measurement unit

Claims (6)

被験者の肺の上部を含み、腹部を含まない体幹上部の第1生体電気インピーダンスと、前記被験者の肺の中下部および腹部を含む体幹中部の第2生体電気インピーダンスとを測定する生体電気インピーダンス測定部と、
前記第1生体電気インピーダンスの測定値の振幅基準レベルを示す第1センタリング値、および、前記第2生体電気インピーダンスの測定値の振幅基準レベルを示す第2センタリング値を生成するセンタリング値生成部と、
前記第1生体電気インピーダンスの測定値の前記第1センタリング値に対する相対値である第1相対値を求める第1相対値算出部と、
前記第2生体電気インピーダンスの測定値の前記第2センタリング値に対する相対値である第2相対値を求める第2相対値算出部と
被験者の1呼吸ごとに、当該1呼吸の吸気における前記第1相対値と前記第2相対値との比に基づいて、前記被験者の胸部のうち呼吸に寄与する部位の機能が正常であるか否かを判定する判定部と、
を備えることを特徴とする換気特性判定装置。
A bioelectrical impedance for measuring the first bioelectrical impedance of the upper trunk including the upper part of the lung of the subject and not including the abdomen, and the second bioelectrical impedance of the middle part of the trunk including the middle lower part and the abdomen of the subject's lung A measuring section;
A centering value generator for generating a first centering value indicating an amplitude reference level of the measurement value of the first bioelectrical impedance and a second centering value indicating an amplitude reference level of the measurement value of the second bioelectrical impedance;
A first relative value calculation unit for obtaining a first relative value that is a relative value of the measured value of the first bioelectrical impedance with respect to the first centering value;
A second relative value calculation unit for obtaining a second relative value that is a relative value of the measured value of the second bioelectrical impedance to the second centering value ;
Whether or not the function of the portion of the subject's chest that contributes to breathing is normal based on the ratio of the first relative value and the second relative value in the inspiration of the one breath for each breath of the subject A determination unit for determining whether or not
Ventilation characteristics determining device you comprising: a.
前記判定部は、前記第1相対値のピーク値に対する前記第2相対値のピーク値の割合が所定の閾値以上である場合は、前記被験者の胸部のうち呼吸に寄与する部位の機能は正常ではないと判定する一方、前記閾値を下回る場合は、当該部位の機能は正常であると判定する、
ことを特徴とする請求項に記載の換気特性判定装置。
When the ratio of the peak value of the second relative value to the peak value of the first relative value is equal to or greater than a predetermined threshold, the determination unit is normal in a function of a part contributing to breathing in the chest of the subject On the other hand, if it falls below the threshold, it is determined that the function of the part is normal.
The ventilation characteristic determination apparatus according to claim 1 .
前記第1生体電気インピーダンスの測定値と、前記第1センタリング値とが等しくなるゼロクロスタイミングを抽出するゼロクロスタイミング抽出部をさらに備え、
前記生体電気インピーダンス測定部は、
所定の周期でサンプリングタイミングに到達するたびに、前記第1生体電気インピーダンスおよび前記第2生体電気インピーダンスを測定し、
前記センタリング値生成部は、
所定数の前記サンプリングタイミングの各々における前記第1生体電気インピーダンスの測定値に基づいて前記第1センタリング値を生成する一方、前記ゼロクロスタイミング抽出部で抽出された前記ゼロクロスタイミングにおける前記第2生体電気インピーダンスの測定値に基づいて前記第2センタリング値を生成する、
ことを特徴とする請求項または請求項に記載の換気特性判定装置。
A zero cross timing extraction unit for extracting a zero cross timing at which the measured value of the first bioelectrical impedance is equal to the first centering value;
The bioelectrical impedance measuring unit is
Measuring the first bioelectrical impedance and the second bioelectrical impedance each time the sampling timing is reached at a predetermined period;
The centering value generation unit
While generating the first centering value based on the measured value of the first bioelectrical impedance at each of the predetermined number of sampling timings, the second bioelectrical impedance at the zero crossing timing extracted by the zero crossing timing extracting unit Generating the second centering value based on the measured value of
The ventilation characteristic determination apparatus according to claim 1 or claim 2 , wherein
前記センタリング値生成部は、
前記サンプリングタイミングごとに、当該サンプリングタイミングよりも所定の時間長だけ前の時点を始点とし、当該サンプリングタイミングを終点とするセンタリング期間内の複数の前記サンプリングタイミングの各々における前記第1生体電気インピーダンスの測定値を用いて移動平均処理を行い、その結果に基づいて、当該サンプリングタイミングにおける前記第1センタリング値を生成する、
ことを特徴とする請求項に記載の換気特性判定装置。
The centering value generation unit
For each of the sampling timings, the measurement of the first bioelectrical impedance at each of the plurality of sampling timings within a centering period starting from a time point a predetermined time before the sampling timing and ending at the sampling timing A moving average process is performed using the value, and the first centering value at the sampling timing is generated based on the result.
The ventilation characteristic judging device according to claim 3 characterized by things.
前記センタリング期間の時間長は、被験者の呼吸速度に応じて可変に設定される、
ことを特徴とする請求項に記載の換気特性判定装置。
The time length of the centering period is variably set according to the breathing rate of the subject.
The ventilation characteristic judging device according to claim 4 characterized by things.
前記センタリング値生成部は、
前記サンプリングタイミングごとに、当該サンプリングタイミングが前記ゼロクロスタイミングであるか否かを判定し、当該サンプリングタイミングが前記ゼロクロスタイミングである場合は、当該サンプリングタイミングにおける前記第2生体電気インピーダンスの測定値に基づいて、当該サンプリングタイミングにおける前記第2センタリング値を生成する一方、当該サンプリングタイミングがゼロクロスタイミングでない場合は、当該サンプリングタイミングの直前の前記サンプリングタイミングで生成した前記第2センタリング値を、当該サンプリングタイミングにおける前記第2センタリング値として採用する、
ことを特徴とする請求項から請求項の何れかに記載の換気特性判定装置。
The centering value generation unit
For each sampling timing, it is determined whether or not the sampling timing is the zero cross timing, and when the sampling timing is the zero cross timing, based on the measurement value of the second bioelectrical impedance at the sampling timing And generating the second centering value at the sampling timing, if the sampling timing is not zero-cross timing, the second centering value generated at the sampling timing immediately before the sampling timing is the second centering value at the sampling timing. 2 Adopted as a centering value
The ventilation characteristic determination device according to any one of claims 3 to 5 , wherein
JP2010167481A 2010-07-15 2010-07-26 Ventilation characteristics judgment device Active JP5633875B2 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2010167481A JP5633875B2 (en) 2010-07-26 2010-07-26 Ventilation characteristics judgment device
EP11173341A EP2407100A1 (en) 2010-07-15 2011-07-08 Respiration characteristic analysis
US13/181,171 US20120016254A1 (en) 2010-07-15 2011-07-12 Respiration characteristic analysis apparatus and respiration characteristic analysis system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010167481A JP5633875B2 (en) 2010-07-26 2010-07-26 Ventilation characteristics judgment device

Publications (2)

Publication Number Publication Date
JP2012024413A JP2012024413A (en) 2012-02-09
JP5633875B2 true JP5633875B2 (en) 2014-12-03

Family

ID=45778093

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010167481A Active JP5633875B2 (en) 2010-07-15 2010-07-26 Ventilation characteristics judgment device

Country Status (1)

Country Link
JP (1) JP5633875B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6273823B2 (en) * 2013-12-19 2018-02-07 株式会社デンソー State estimation device and state estimation program
CN115300873B (en) * 2022-09-29 2022-12-30 江苏南通鼎顺网络科技有限责任公司 Data processing method for respiratory rehabilitation training

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS54120997A (en) * 1978-03-14 1979-09-19 Tokyo Shibaura Electric Co Respiration function monitor device
JP3072724B2 (en) * 1998-02-19 2000-08-07 日本電気株式会社 Respiration measurement device
EP1292224B2 (en) * 2000-06-09 2009-12-30 Böhm, Stephan, Dr. Method and apparatus for displaying information obtained by electrical impedance tomography data
US8808193B2 (en) * 2007-09-11 2014-08-19 Carefusion 207, Inc. Regional oxygen uptake/perfusion measuring device and method

Also Published As

Publication number Publication date
JP2012024413A (en) 2012-02-09

Similar Documents

Publication Publication Date Title
JP4732527B2 (en) Respiratory judgment device
JP5841716B2 (en) Method and system for interpreting and analyzing physiological information, performance information and situation information
US20120016255A1 (en) Respiration characteristic analysis apparatus and respiration characteristic analysis system
US20150342518A1 (en) System and method to monitor, guide, and evaluate breathing, utilizing posture and diaphragm sensor signals
EP2407100A1 (en) Respiration characteristic analysis
CN105496377A (en) Heart rate variability biofeedback exercise systematic method and apparatus
JP6637175B2 (en) A method and apparatus for estimating the transit time of a pulse in an artery from measurements of the tips of both hands or feet.
JP2002253534A (en) Balance detector capable of evaluating exercise capacity
US20230381590A1 (en) Method and equipment for monitoring the effectiveness of physical activities, especially sports activities
JP5610159B2 (en) Breathing training apparatus and breathing training system
JP6307457B2 (en) Exercise state feedback method, system, and program
JP5633875B2 (en) Ventilation characteristics judgment device
JP6330009B2 (en) Feedback method, system, and program for exercise state and psychological state
JP5622202B2 (en) Breathing training apparatus and breathing training system
CA3141123C (en) Portable dehydration monitoring system
JPS62277968A (en) Biofeed bag apparatus
JP4965205B2 (en) Heart rate fluctuation detection device
TWI580400B (en) Method and apparatus for calculating heart health index by using resting heart rate
JP6431785B2 (en) Psychological state feedback method, system, and program
WO2023120232A1 (en) Breathing instruction device
CN220778334U (en) Wearable equipment
CN118044819A (en) Wearable equipment
JP6093787B2 (en) Feedback method, system, and program for exercise state and psychological state
EP4322842A1 (en) System, method and computer program for monitoring health of a person
GB2621222A (en) Method and system for estimating cardiovascular fitness and maximum heart rate for a user

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20130517

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20131225

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140107

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140306

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: 20140917

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20141006

R150 Certificate of patent or registration of utility model

Ref document number: 5633875

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250