JP2017047193A - Biological measurement apparatus, biological measurement method, and detection method - Google Patents

Biological measurement apparatus, biological measurement method, and detection method Download PDF

Info

Publication number
JP2017047193A
JP2017047193A JP2016162687A JP2016162687A JP2017047193A JP 2017047193 A JP2017047193 A JP 2017047193A JP 2016162687 A JP2016162687 A JP 2016162687A JP 2016162687 A JP2016162687 A JP 2016162687A JP 2017047193 A JP2017047193 A JP 2017047193A
Authority
JP
Japan
Prior art keywords
living body
vibration
vibration wave
frequency
vibrator
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.)
Pending
Application number
JP2016162687A
Other languages
Japanese (ja)
Inventor
櫻井 直樹
Naoki Sakurai
直樹 櫻井
秀美 秋元
Hidemi Akimoto
秀美 秋元
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.)
Applied Vibro-Acoustics Inc
Original Assignee
Applied Vibro-Acoustics Inc
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 Applied Vibro-Acoustics Inc filed Critical Applied Vibro-Acoustics Inc
Publication of JP2017047193A publication Critical patent/JP2017047193A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

PROBLEM TO BE SOLVED: To measure information inside a living body by using a noninvasive and simple method without burdening a subject's body as much as possible.SOLUTION: An exciter 1 gives an oscillatory wave having a single frequency selected from an audible zone range to a living body. A geophone 2 receives the oscillatory wave propagating in the living body 10 at a position a prescribed distance apart from a position at which the exciter 1 has given the oscillatory wave to the living body. A search unit 3A, while varying a frequency of the oscillatory wave given to the living body 10 by the exciter 1, searches for the oscillatory wave's frequency at which a phase difference between the oscillatory wave given to the living body 10 by the exciter 1 and the oscillatory wave received by the geophone 2 is a specific phase difference. A calculation unit 3B, on the basis of the frequency having been searched by the search unit 3A, calculates the oscillatory wave's propagation speed in the living body 10.SELECTED DRAWING: Figure 1

Description

本発明は、生体計測装置、生体計測方法及び検出方法に関する。   The present invention relates to a biological measurement apparatus, a biological measurement method, and a detection method.

人体等の生体内部の情報を調べるには様々な方法があり、X線を用いたり、超音波を用いたり、生体組織の一部を取り出したりすることが行われている。例えば、超音波を用いて生体の硬さを計測する硬さ測定装置が開示されている(例えば、特許文献1参照)。   There are various methods for examining information inside a living body such as a human body, and X-rays, ultrasonic waves, and part of a living tissue are taken out. For example, a hardness measuring device that measures the hardness of a living body using ultrasonic waves is disclosed (for example, see Patent Document 1).

特許第4287200号公報Japanese Patent No. 4287200

生体内部の情報の計測は、被爆及び痛みを伴わず、被検者の身体にできるだけ負担をかけることなく、非侵襲かつ簡便な方法で行われるのが望ましい。また、例えば、超音波は、生体内を通りにくく、計測できる範囲に限界がある。   The measurement of information inside the living body is desirably performed in a non-invasive and simple manner without causing exposure and pain and without placing a burden on the subject's body as much as possible. In addition, for example, ultrasonic waves are difficult to pass through the living body, and there is a limit to the range that can be measured.

本発明は、被爆及び痛みを伴わず、被検者の身体にできるだけ負担をかけることなく、非侵襲かつ簡便な方法で生体内部の情報を計測することができる生体計測装置、生体計測方法及び検出方法を提供することを目的とする。   The present invention relates to a living body measuring apparatus, a living body measuring method, and a detection capable of measuring information inside a living body by a non-invasive and simple method without causing exposure and pain, and without burdening the body of the subject as much as possible. It aims to provide a method.

上記目的を達成するために、本発明の第1の観点に係る生体計測装置は、
可聴域の範囲から選択された単一の周波数の振動波を生体に与える加振部と、
前記加振部が前記生体に前記振動波を与える位置から所定距離離れた位置で、前記生体を伝わる振動波を受振する受振部と、
前記加振部によって前記生体に与えられる振動波の周波数を変更しながら、前記加振部により前記生体に与えられる振動波と前記受振部で受振される振動波との位相差が特定の位相差となるような前記振動波の周波数を探索する探索部と、
前記探索部で探索された周波数に基づいて、前記生体中における前記振動波の伝搬速度を算出する算出部と、
を備える。
In order to achieve the above object, a biological measurement apparatus according to the first aspect of the present invention provides:
An excitation unit that gives a living body a vibration wave having a single frequency selected from the range of the audible range;
A vibration receiving unit that receives vibration waves transmitted through the living body at a position away from a position where the vibration unit gives the vibration wave to the living body by a predetermined distance;
While changing the frequency of the vibration wave given to the living body by the vibration unit, the phase difference between the vibration wave given to the living body by the vibration unit and the vibration wave received by the vibration receiving unit is a specific phase difference. A search unit for searching for the frequency of the vibration wave such that
A calculation unit that calculates a propagation speed of the vibration wave in the living body based on the frequency searched by the search unit;
Is provided.

この場合、球面のプローブを前記生体に所定の深さ押し込んだ状態で計測される応力を計測する応力計測部を備え、
前記算出部は、
前記応力計測部で計測された応力に基づいて、前記生体のヤング率を算出する、
こととしてもよい。
In this case, a stress measurement unit that measures stress measured in a state where a spherical probe is pushed into the living body to a predetermined depth,
The calculation unit includes:
Based on the stress measured by the stress measurement unit, to calculate the Young's modulus of the living body,
It is good as well.

前記算出部は、
算出した前記伝搬速度と、前記生体のヤング率とに基づいて、前記生体の密度を算出する、
こととしてもよい。
The calculation unit includes:
Calculating the density of the living body based on the calculated propagation velocity and the Young's modulus of the living body;
It is good as well.

前記算出部は、
算出した前記生体の密度に基づいて、前記生体における筋肉組織及び/又は脂肪組織の比率に関する情報を算出する、
こととしてもよい。
The calculation unit includes:
Based on the calculated density of the living body, calculate information on the ratio of muscle tissue and / or adipose tissue in the living body,
It is good as well.

前記加振器は、
前記生体における異なる2つの加振点を前記各加振器でそれぞれ同位相又は逆位相で加振し、
前記受振器は、
前記2つの加振点の中央1点である受振点で受振する、
こととしてもよい。
The vibrator is
Two different excitation points in the living body are respectively excited in the same phase or in opposite phase by the respective vibrators,
The geophone is
Receiving at a receiving point which is one central point of the two excitation points;
It is good as well.

本発明の第2の観点に係る生体計測方法は、
コンピュータが、加振器によって生体に与えられる振動波の周波数を可聴域の範囲で変更しながら、前記加振器により前記生体に与えられる振動波と前記受振器で受振される振動波との位相差が特定の位相差となるような前記振動波の周波数を探索する探索ステップと、
コンピュータが、前記探索ステップで探索された周波数に基づいて、前記生体中における前記振動波の伝搬速度を前記生体の硬さに関する情報として算出する算出ステップと、
を含む。
The biological measurement method according to the second aspect of the present invention includes:
While the computer changes the frequency of the vibration wave applied to the living body by the vibrator in the range of the audible range, the vibration wave applied to the living body by the vibrator and the vibration wave received by the shaker A search step for searching for a frequency of the vibration wave such that the phase difference becomes a specific phase difference;
A calculation step in which the computer calculates the propagation speed of the vibration wave in the living body as information on the hardness of the living body based on the frequency searched in the searching step;
including.

この場合、球面型のプローブを前記生体の組織を所定の深さまで押し込んだ時に得られる応力を計測する応力ひずみ計を用いて、当該応力を計測する計測ステップを含み、
前記算出ステップでは、
前記コンピュータが、
前記計測ステップで計測された応力に基づいて、前記生体のヤング率を算出する、
こととしてもよい。
In this case, using a stress strain meter that measures the stress obtained when the spherical probe is pushed into the living tissue to a predetermined depth, the measurement step includes measuring the stress,
In the calculating step,
The computer is
Based on the stress measured in the measurement step, calculate the Young's modulus of the living body,
It is good as well.

前記算出ステップでは、
コンピュータが、前記伝搬速度と前記生体のヤング率とに基づいて、前記生体の密度を前記生体の硬さに関する情報として算出する、
こととしてもよい。
In the calculating step,
The computer calculates the density of the living body as information on the hardness of the living body based on the propagation speed and the Young's modulus of the living body.
It is good as well.

前記算出ステップでは、
コンピュータは、算出した前記生体の密度に基づいて、前記生体における筋肉組織及び/又は脂肪組織の比率に関する情報を前記生体の硬さに関する情報として算出する、
こととしてもよい。
In the calculating step,
Based on the calculated density of the living body, the computer calculates information on the ratio of muscle tissue and / or adipose tissue in the living body as information on the hardness of the living body,
It is good as well.

前記探索ステップでは、
前記加振器と前記受振器とを水平に配置して前記振動波の周波数を探索する、
こととしてもよい。
In the search step,
Search the frequency of the vibration wave by horizontally arranging the vibrator and the geophone,
It is good as well.

前記探索ステップでは、
乳房の上側部分に前記加振器と前記受振器とを配置して前記探索ステップ及び前記算出ステップを行って、その部分における前記振動波の伝搬速度を計測し、
乳房の下側部分に前記加振器と前記受振器とを配置して前記探索ステップ及び前記算出ステップを行って、その部分における前記振動波の伝搬速度を計測する、
こととしてもよい。
In the search step,
Place the exciter and the geophone on the upper part of the breast, perform the search step and the calculation step, measure the propagation speed of the vibration wave in that part,
Placing the exciter and the geophone in the lower part of the breast and performing the search step and the calculation step to measure the propagation speed of the vibration wave in that part;
It is good as well.

前記探索ステップでは、
前記生体における異なる複数の加振点を前記加振器でそれぞれ同位相又は逆位相で加振し、
前記複数の加振点の中央1点を受振点として前記受振器で受振する、
こととしてもよい。
In the search step,
A plurality of different excitation points in the living body are respectively excited in the same phase or in opposite phase by the vibrator,
Receiving at the geophone as a receiving point at the center of the plurality of excitation points,
It is good as well.

本発明の第3の観点に係る検出方法は、
乳房を計測対象として、本発明の第2の観点に係る生体計測方法を行って、前記生体中における前記振動波の伝搬速度を求め、その伝搬速度の変化に基づいて、乳がんを検出する。
The detection method according to the third aspect of the present invention is:
Using the breast as a measurement target, the living body measuring method according to the second aspect of the present invention is performed to determine the propagation speed of the vibration wave in the living body, and breast cancer is detected based on the change in the propagation speed.

本発明によれば、可聴域における単一の周波数の振動波を生体に与える加振部と生体を伝わる振動波を受振する受振部との距離を予め定められた距離とする。そして、生体に与えられる振動波の周波数を変更しながら、加振部により生体に与えられる振動波と受振部で受振される振動波との位相差が特定の位相差となるような振動波の周波数を探索する。このようにすれば、生体を伝わる振動波の周波数と波長との関係がわかるため、生体における振動波の伝搬速度を求めることができる。生体における振動波の伝搬速度は、生体の硬さに応じて変化するため、生体内部の情報を得るための基礎的な情報となる。このように、本発明によれば、被爆及び痛みを伴わず、被検者の身体にできるだけ負担をかけることなく、非侵襲かつ簡便な方法で生体内部の情報を計測することができる。   According to the present invention, the distance between the excitation unit that applies a vibration wave having a single frequency in the audible range to the living body and the vibration reception unit that receives the vibration wave transmitted through the living body is set to a predetermined distance. Then, while changing the frequency of the vibration wave applied to the living body, the vibration wave such that the phase difference between the vibration wave applied to the living body by the vibration unit and the vibration wave received by the vibration receiving unit becomes a specific phase difference. Search for a frequency. In this way, since the relationship between the frequency and the wavelength of the vibration wave transmitted through the living body is known, the propagation speed of the vibration wave in the living body can be obtained. Since the propagation speed of the vibration wave in the living body changes according to the hardness of the living body, it becomes basic information for obtaining information inside the living body. As described above, according to the present invention, information inside the living body can be measured by a non-invasive and simple method without causing exposure and pain, and without burdening the body of the subject as much as possible.

本発明の実施の形態1に係る生体計測装置の構成を示すブロック図である。It is a block diagram which shows the structure of the biometric apparatus which concerns on Embodiment 1 of this invention. 図2(A)は、加振器モニタ信号と受振信号とで、振動波の位相差が180度ずれる場合を示す図であり、図2(B)は、加振器モニタ信号と受振信号とで、振動波の位相差が540度ずれる場合を示す図である。FIG. 2A is a diagram showing a case where the vibration wave phase difference is 180 degrees between the vibration monitor signal and the vibration signal. FIG. 2B is a diagram showing the vibration monitor signal and the vibration signal. FIG. 6 is a diagram illustrating a case where the phase difference of the vibration wave is shifted by 540 degrees. 振動波の周波数と、加振器モニタ信号と受振信号の振動波の位相差との関係を示す図である。It is a figure which shows the relationship between the frequency of a vibration wave, and the phase difference of the vibration wave of an exciter monitor signal and a vibration receiving signal. 本発明の実施の形態1に係る生体計測装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the biometric apparatus which concerns on Embodiment 1 of this invention. 本発明の実施の形態2に係る生体計測装置の構成を示すブロック図である。It is a block diagram which shows the structure of the biometric apparatus which concerns on Embodiment 2 of this invention. 本発明の実施の形態2に係る生体計測装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the biometric apparatus which concerns on Embodiment 2 of this invention. 本発明の実施の形態3に係る生体計測装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the bioinstrumentation apparatus which concerns on Embodiment 3 of this invention. 乳房における加振器と受振器との配置パターンを示す図である。It is a figure which shows the arrangement | positioning pattern of the vibrator and geophone in a breast. 図9(A)は、パターンAの配置(加振器を上にした場合)における位相グラフの一例である。図9(B)は、パターンAの配置(加振器を下にした場合)における位相グラフの一例である。FIG. 9A is an example of a phase graph in the pattern A arrangement (when the vibrator is placed on the top). FIG. 9B is an example of a phase graph in the pattern A arrangement (when the vibrator is placed below). 図10(A)は、パターンCの配置(加振器を上にした場合)における位相グラフの一例である。図10(B)は、パターンCの配置(加振器を下にした場合)における位相グラフの一例である。FIG. 10A is an example of a phase graph in the arrangement of pattern C (when the vibrator is placed on the top). FIG. 10B is an example of a phase graph in the pattern C arrangement (when the vibrator is placed on the bottom). 図11(A)は、パターンBの配置における位相グラフの一例である。図11(B)は、パターンDの配置における位相グラフの一例である。FIG. 11A is an example of a phase graph in the pattern B arrangement. FIG. 11B is an example of a phase graph in the pattern D arrangement. 図12(A)は、加振点から半波長分だけ離れた受振点での振動が測定される様子を示す模式図である。図12(B)は、加振点から1波長分だけ離れた受振点での振動が測定される様子を示す模式図である。FIG. 12A is a schematic diagram illustrating a state in which vibration is measured at a receiving point that is separated from the excitation point by a half wavelength. FIG. 12B is a schematic diagram illustrating a state in which vibration is measured at a receiving point that is one wavelength away from the excitation point. 1点加振法を行う場合の本発明の実施の形態5に係る生体計測装置の構成を示す図である。It is a figure which shows the structure of the biometric apparatus which concerns on Embodiment 5 of this invention in the case of performing 1 point | piece excitation method. 2点同位相加振法及び2点逆位相加振法を行う場合の本発明の実施の形態5に係る生体計測装置の構成を示す図である。It is a figure which shows the structure of the bioinstrumentation apparatus which concerns on Embodiment 5 of this invention in the case of performing a 2-point in-phase excitation method and a 2-point anti-phase excitation method. 1点加振法の測定結果を示す位相グラフである。It is a phase graph which shows the measurement result of a 1 point vibration method. 2点同位相加振法の測定結果を示す位相グラフである。It is a phase graph which shows the measurement result of a two-point in-phase excitation method. 2点逆位相加振法の測定結果を示す位相グラフである。It is a phase graph which shows the measurement result of a 2 point | piece antiphase excitation method.

以下、本発明の実施の形態について図面を参照して詳細に説明する。各図において同一の構成要素には同一の符号を付している。   Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings. In each figure, the same components are denoted by the same reference numerals.

実施の形態1.
図1に示すように、本発明の実施の形態1に係る生体計測装置100は、加振器1と、受振器2と、コントローラ3と、を備える。生体計測装置100は、生体10の内部の情報を計測する。
Embodiment 1 FIG.
As shown in FIG. 1, the biological measurement apparatus 100 according to Embodiment 1 of the present invention includes a vibrator 1, a shaker 2, and a controller 3. The living body measuring apparatus 100 measures information inside the living body 10.

加振器(スピーカ)1は、振動波(縦波)を生体10に与える。加振器1によって与えられる振動波は、コントローラ3により可聴域の範囲から選択された単一の周波数の振動波である。可聴域の範囲としては、0Hzより高く、かつ、20,000Hz以下の様々な範囲を設定することできるが、好ましくは、0Hzより高く、かつ、1000Hz以下である。加振器1は、コントローラ3から入力される加振信号に従って振動する。これにより、振動波(縦波)が生体10内を伝わっていく。   The vibrator (speaker) 1 gives a vibration wave (longitudinal wave) to the living body 10. The vibration wave given by the vibrator 1 is a vibration wave having a single frequency selected from the range of the audible range by the controller 3. As the range of the audible range, various ranges higher than 0 Hz and 20,000 Hz or lower can be set, but preferably higher than 0 Hz and 1000 Hz or lower. The vibrator 1 vibrates according to the vibration signal input from the controller 3. Thereby, a vibration wave (longitudinal wave) is transmitted through the living body 10.

加振器1には、加速度ピックアップ1Aが取り付けられている。加速度ピックアップ1Aは、加振器1から生体10に与えられる振動波の波形を検出する。加速度ピックアップ1Aは、検出した波形に対応する加振器モニタ信号をコントローラ3に出力する。   An acceleration pickup 1A is attached to the vibrator 1. The acceleration pickup 1 </ b> A detects a waveform of a vibration wave given from the vibrator 1 to the living body 10. The acceleration pickup 1 </ b> A outputs an exciter monitor signal corresponding to the detected waveform to the controller 3.

受振器2は、加振器1が生体10と接触する位置から、予め定められた距離Lだけ離れた位置で生体10と接触している。受振器2は、加振器1によって与えられ、生体10を伝達される振動を受振する加速度ピックアップである。   The geophone 2 is in contact with the living body 10 at a position separated by a predetermined distance L from the position where the vibrator 1 is in contact with the living body 10. The geophone 2 is an acceleration pickup that receives vibrations given by the vibrator 1 and transmitted through the living body 10.

受振器2で受振される生体10を伝わる振動波の位相は、加振器1から生体10に与えられる振動波の位相に対して遅れる。図2(A)に示すように、加振器モニタ信号と受振信号とで、位相差が180度ずれる場合、加振器1と受振器2との距離Lは、生体10を伝わる振動波の半波長(L=λ/2)となる。また、図2(B)に示すように、加振器モニタ信号と受振信号とで、位相差が540度ずれる場合、加振器1と受振器2との間の距離は、生体10を伝わる振動波の半波長(L=3λ/2)となる。このように、加振器1と受振器2との間の距離Lと、加振器モニタ信号と受振信号との位相差を求めれば、生体10を伝わる振動波の波長を求めることができる。   The phase of the vibration wave transmitted through the living body 10 that is received by the vibration receiver 2 is delayed with respect to the phase of the vibration wave that is applied from the vibrator 1 to the living body 10. As shown in FIG. 2A, when the phase difference is 180 degrees between the vibrator monitor signal and the vibration receiving signal, the distance L between the vibration vibrator 1 and the vibration receiving body 2 is the vibration wave transmitted through the living body 10. Half wavelength (L = λ / 2). In addition, as shown in FIG. 2B, when the phase difference between the vibrator monitor signal and the vibration receiving signal is shifted by 540 degrees, the distance between the vibrator 1 and the vibration receiver 2 is transmitted through the living body 10. The half wave of the vibration wave (L = 3λ / 2). Thus, if the distance L between the vibrator 1 and the shaker 2 and the phase difference between the shaker monitor signal and the received signal are obtained, the wavelength of the vibration wave transmitted through the living body 10 can be obtained.

加振器1と受振器2との間の距離が同じであっても、加振器モニタ信号と受振信号との位相差は、振動波の周波数に応じて変化する。図3に示すように、振動波の周波数が上がっていくにつれ、受振信号は、加振器モニタ信号に対して遅れていく。周波数がf1、f2、f3になると、位相差は、180度のm(mは自然数)倍となる。   Even if the distance between the shaker 1 and the shaker 2 is the same, the phase difference between the shaker monitor signal and the received vibration signal changes according to the frequency of the vibration wave. As shown in FIG. 3, as the frequency of the vibration wave increases, the vibration receiving signal is delayed with respect to the vibrator monitor signal. When the frequencies are f1, f2, and f3, the phase difference is m of 180 degrees (m is a natural number).

本実施の形態に係る生体計測装置100は、加振器1と受振器2との距離Lと、生体10を伝わる振動波の周波数と、加振器モニタ信号と受振信号との位相差とに基づいて、生体10を伝わる振動波の伝搬速度Vを、生体内部の情報として計測する。例えば、加振器モニタ信号と受振信号との位相差が180度となる振動波の周波数f1を求める。そして、そのときの振動波の波長λ=2Lに基づいて、振動波の伝搬速度Vが求まる。   The living body measurement apparatus 100 according to the present embodiment includes a distance L between the vibrator 1 and the shaker 2, a frequency of a vibration wave transmitted through the living body 10, and a phase difference between the shaker monitor signal and the received vibration signal. Based on this, the propagation velocity V of the vibration wave transmitted through the living body 10 is measured as information inside the living body. For example, the frequency f1 of the vibration wave at which the phase difference between the exciter monitor signal and the vibration receiving signal is 180 degrees is obtained. Then, based on the wavelength λ = 2L of the vibration wave at that time, the propagation velocity V of the vibration wave is obtained.

図1に戻り、コントローラ3は、加振器1と受振器2とを制御して、上記計測を行う。具体的には、コントローラ3は、加振器1に、選択した周波数の加振信号を出力するとともに、加速度ピックアップ1Aから出力される加振器モニタ信号を入力し、受振器2から出力される受振信号を入力する。コントローラ3は、加振器1に出力する加振信号を調整しながら、加速度ピックアップ1Aから出力される加振器モニタ信号と受振器2で受振される受振信号とに対して信号処理を行うことにより、生体10内部の情報を計測する。   Returning to FIG. 1, the controller 3 controls the vibrator 1 and the geophone 2 to perform the measurement. Specifically, the controller 3 outputs an excitation signal having a selected frequency to the vibrator 1, inputs a vibrator monitor signal output from the acceleration pickup 1 </ b> A, and outputs it from the receiver 2. Input vibration receiving signal. The controller 3 performs signal processing on the vibration monitor signal output from the acceleration pickup 1 </ b> A and the vibration reception signal received by the vibration receiver 2 while adjusting the vibration signal output to the vibration generator 1. Thus, information inside the living body 10 is measured.

コントローラ3は、CPU(Central Processing Unit)、記憶装置及び入出力装置を有するコンピュータである。コントローラ3は、記憶装置に格納されたプログラムを実行することにより、加振信号の調整と、加振器モニタ信号と受振信号とに対する信号処理を行う。   The controller 3 is a computer having a CPU (Central Processing Unit), a storage device, and an input / output device. The controller 3 executes the program stored in the storage device, thereby adjusting the vibration signal and performing signal processing on the vibration monitor signal and the vibration signal.

図1に示すように、コントローラ3は、加振信号の調整と、加振器モニタ信号と受振信号とに対する信号処理を行うプログラムを実行することにより実現される機能の構成として、探索部3Aと、算出部3Bと、を備えている。   As shown in FIG. 1, the controller 3 has a search unit 3A as a configuration of functions realized by executing a program for adjusting an excitation signal and performing signal processing on the shaker monitor signal and the received signal. And a calculation unit 3B.

探索部3Aは、加振器1に加振信号を出力する。この加振信号は単一の周波数を有する正弦波状の信号である。加振器1は、加振信号に従って振動し、加振信号の周波数に合致する振動波(縦波)を生体10に伝達する。探索部3Aは、可聴域の範囲で、加振信号の周波数を変更可能である。加振信号の周波数を変更することにより、加振器1によって生体10に与えられる振動波の周波数を変更する。   The search unit 3 </ b> A outputs an excitation signal to the vibrator 1. This excitation signal is a sinusoidal signal having a single frequency. The vibrator 1 vibrates according to the vibration signal, and transmits a vibration wave (longitudinal wave) that matches the frequency of the vibration signal to the living body 10. Search unit 3A can change the frequency of the excitation signal within the audible range. By changing the frequency of the vibration signal, the frequency of the vibration wave applied to the living body 10 by the vibrator 1 is changed.

探索部3Aは、加振信号の周波数を変更しながら、加振器モニタ信号と受振信号を受信する。探索部3Aは、加振器モニタ信号と受振信号との位相差を算出する。探索部3Aは、算出した位相差が特定の位相差となるまで、加振信号の周波数を変更する。例えば、探索部3Aは、加振器1によって生体10に与えられる振動波の周波数を変更しながら、加振器1により生体10に与えられる振動波と受振器2で受振される振動波との位相差が特定の位相差(−180度)となるような振動波の周波数f1を探索する。   The search unit 3A receives the vibration monitor signal and the vibration receiving signal while changing the frequency of the vibration signal. Search unit 3A calculates the phase difference between the vibrator monitor signal and the vibration receiving signal. Search unit 3A changes the frequency of the excitation signal until the calculated phase difference becomes a specific phase difference. For example, the search unit 3A changes the frequency of the vibration wave given to the living body 10 by the vibrator 1 and the vibration wave given to the living body 10 by the shaker 1 and the vibration wave received by the shaker 2. The frequency f1 of the vibration wave is searched for such that the phase difference becomes a specific phase difference (−180 degrees).

算出部3Bは、探索部3Aで探索された周波数に基づいて、生体10中における振動波の伝搬速度Vを算出する。
V=f1×2L…(1)
The calculation unit 3B calculates the vibration wave propagation velocity V in the living body 10 based on the frequency searched by the search unit 3A.
V = f1 × 2L (1)

次に、生体計測装置100の計測動作について説明する。   Next, the measurement operation of the biological measurement apparatus 100 will be described.

図4に示すように、加振器1と受振器2とのセッティングを行う(ステップS1)。具体的には、加振器1と受振器2を所定の間隔L開いた状態でサンプル(生体10)に接触させる。   As shown in FIG. 4, the setting of the vibrator 1 and the geophone 2 is performed (step S1). Specifically, the exciter 1 and the geophone 2 are brought into contact with the sample (living body 10) with a predetermined distance L open.

続いて、加振器1を振動波の周波数を0Hzから1000Hzまでスイープさせる(ステップS2)。   Subsequently, the vibrator 1 is swept from 0 Hz to 1000 Hz of the vibration wave frequency (step S2).

続いて、探索部3Aは、それぞれの周波数について加振器モニタ信号の位相と受振信号の位相とを検出し、両者の位相差を検出する(ステップS3)。   Subsequently, the search unit 3A detects the phase of the exciter monitor signal and the phase of the received signal for each frequency, and detects the phase difference between the two (step S3).

続いて、探索部3Aは、位相差が−180度となる振動波の周波数f1を決定する(ステップS4)。   Subsequently, the search unit 3A determines the frequency f1 of the vibration wave at which the phase difference is −180 degrees (step S4).

続いて、算出部3Bは、上記式(1)を用いて、振動波の伝搬速度Vを算出する(ステップS5)。振動波の伝搬速度Vは、生体10の硬さを示す情報となる。伝播速度Vが大きければ大きいほど、生体10の計測部位は硬いということになる。   Subsequently, the calculation unit 3B calculates the propagation velocity V of the vibration wave using the above formula (1) (step S5). The propagation speed V of the vibration wave is information indicating the hardness of the living body 10. The greater the propagation velocity V, the harder the measurement site of the living body 10 is.

以上詳細に説明したように、本実施の形態によれば、可聴域における単一の周波数の振動波を生体10に与える加振器1と生体10を伝わる振動波を受振する受振器2との距離を予め定められた距離Lとする。そして、生体10に与えられる振動波の周波数を変更しながら、加振器1により生体に与えられる振動波と受振器2で受振される振動波との位相差が特定の位相差となるような振動波の周波数を探索する。   As described above in detail, according to the present embodiment, the vibration exciter 1 that applies a vibration wave having a single frequency in the audible range to the living body 10 and the vibration receiver 2 that receives the vibration wave transmitted through the living body 10 are provided. The distance is set to a predetermined distance L. Then, while changing the frequency of the vibration wave applied to the living body 10, the phase difference between the vibration wave applied to the living body by the vibrator 1 and the vibration wave received by the vibration receiver 2 becomes a specific phase difference. Search for vibration wave frequency.

このようにすれば、加振器1と受振器2との距離L、振動波の位相差180°及び周波数f1の関係に基づいて、生体10における振動波の伝搬速度Vを求めることができる。生体10における振動波の伝搬速度Vは、生体10の硬さに応じて変化する、生体10の内部の情報を得るための基礎的な情報となる。このように、本実施の形態に係る生体計測装置100によれば、被爆及び痛みを伴わず、被検者の身体にできるだけ負担をかけることなく、非侵襲かつ簡便な方法で生体10内部の情報を計測することができる。   In this way, the propagation speed V of the vibration wave in the living body 10 can be obtained based on the relationship between the distance L between the vibration exciter 1 and the vibration receiver 2, the phase difference of vibration wave 180 °, and the frequency f1. The propagation speed V of the vibration wave in the living body 10 is basic information for obtaining information inside the living body 10 that changes according to the hardness of the living body 10. As described above, according to the living body measuring apparatus 100 according to the present embodiment, the information inside the living body 10 can be obtained in a non-invasive and simple manner without causing exposure and pain, and without burdening the subject's body as much as possible. Can be measured.

実施の形態2.
次に、本発明の実施の形態2について説明する。本実施の形態では、生体10内部の情報の計測に、生体計測装置100の他に、図5に示す応力ひずみ計101が用いられる。
Embodiment 2. FIG.
Next, a second embodiment of the present invention will be described. In the present embodiment, a stress strain meter 101 shown in FIG. 5 is used in addition to the living body measuring apparatus 100 for measuring information inside the living body 10.

図5に示すように、応力ひずみ計101は、棒状体20と、応力計測部21と、押し当て部22と、を備える。棒状体20は、応力計測部21から外部に突き出しており、突き出した先端には、原点Oを中心とする半径Rの球の球面の一部が形成されている。応力計測部21は、棒状体20の先端に物を押し当てた場合の応力Fを計測する。押し当て部22は、生体10に棒状体20の先端を押し当てた場合に、生体10表面に押し当てられる。この押し当て部22により、棒状体20の先端が生体10に入り込む深さはhとなる。   As shown in FIG. 5, the stress strain gauge 101 includes a rod-shaped body 20, a stress measurement unit 21, and a pressing unit 22. The rod-shaped body 20 protrudes outward from the stress measurement unit 21, and a part of a spherical surface of a sphere having a radius R centered at the origin O is formed at the protruding tip. The stress measurement unit 21 measures the stress F when an object is pressed against the tip of the rod-shaped body 20. The pressing unit 22 is pressed against the surface of the living body 10 when the tip of the rod-shaped body 20 is pressed against the living body 10. The depth at which the tip of the rod-shaped body 20 enters the living body 10 by the pressing portion 22 is h.

応力ひずみ計101を、図5に示すように生体10に押し付けた場合に求められる応力Fは、コントローラ3に送られる。   The stress F obtained when the stress strain meter 101 is pressed against the living body 10 as shown in FIG. 5 is sent to the controller 3.

本実施の形態では、算出部3Bは、応力Fに基づいて、生体10のヤング率を算出する。生体10のヤング率Eは、次式を用いて求められる。
E=(3/4)×(F/(R1/2×h2/3)×(1−σ)…(2)
ここで、σはポアソン比である。σは、計測対象の部位が筋肉組織である場合には、例えば0.490を使用することができ、部位が脂肪組織である場合には、例えば0.495を使用することができる。
In the present embodiment, the calculation unit 3B calculates the Young's modulus of the living body 10 based on the stress F. The Young's modulus E of the living body 10 is obtained using the following equation.
E = (3/4) × (F / (R 1/2 × h 2/3 ) × (1-σ 2 ) (2)
Here, σ is a Poisson's ratio. For example, 0.490 can be used when the site to be measured is a muscle tissue, and 0.495 can be used when the site is a fat tissue.

なお、上記式(2)は、以下の公知の固体接触式に基づいて、被接触側の半径R’を無限大と置くことにより得られる式である。
h=F2/3[D(1/R+1/R’)]1/3 …(3)
ここで、D=3/4((1−σ)/E+(1−σ’)/E’)である。σ’、E’は、棒状体20のポアソン比及びヤング率である。
In addition, said Formula (2) is a type | formula obtained by setting radius R 'of the to-be-contacted side to infinity based on the following well-known solid contact type | formulas.
h = F 2/3 [D 2 (1 / R + 1 / R ′)] 1/3 (3)
Here, D = 3/4 ((1-σ) 2 / E + (1-σ ′) 2 / E ′). σ ′ and E ′ are the Poisson's ratio and Young's modulus of the rod-shaped body 20.

さらに、算出部3Bは、上記実施の形態1と同様に、算出した伝搬速度Vと、生体10のヤング率Eとに基づいて、生体10の密度ρを、次式を用いて算出する。
ρ=(a×E)/V … (4)
上記式(4)のaは、生体10の形状、振動波の振動数と測定場所によって決まる定数である。
Furthermore, the calculation unit 3B calculates the density ρ of the living body 10 using the following equation based on the calculated propagation velocity V and the Young's modulus E of the living body 10 as in the first embodiment.
ρ = (a 2 × E) / V 2 (4)
In the above formula (4), a is a constant determined by the shape of the living body 10, the frequency of the vibration wave, and the measurement location.

図6に示すように、探索部3A及び算出部3Bは、上記実施の形態1と同様に、ステップS1〜S5を行う。続いて、応力ひずみ計101を用いて、応力Fの計測を行う(ステップS6)。計測された応力Fは、算出部3Bに送られる。算出部3Bは、応力ひずみ計101で計測される生体10の応力Fに基づいて、生体10のヤング率Eを、上記式(2)を用いて算出する(ステップS7)。続いて、算出部3Bは、生体10の密度ρを、上記式(4)を用いて算出する(ステップS8)。   As shown in FIG. 6, the search unit 3A and the calculation unit 3B perform steps S1 to S5 as in the first embodiment. Subsequently, the stress F is measured using the stress strain meter 101 (step S6). The measured stress F is sent to the calculation unit 3B. The calculation unit 3B calculates the Young's modulus E of the living body 10 using the above equation (2) based on the stress F of the living body 10 measured by the stress strain meter 101 (step S7). Subsequently, the calculation unit 3B calculates the density ρ of the living body 10 using the above formula (4) (step S8).

以上説明したように、本実施の形態によれば、被爆及び痛みを伴わず、被検者の身体にできるだけ負担をかけることなく、非侵襲かつ簡便な方法で、生体10の密度ρを求めることができる。   As described above, according to the present embodiment, the density ρ of the living body 10 is obtained by a non-invasive and simple method without causing exposure and pain, and placing as little burden on the subject's body as possible. Can do.

なお、本実施の形態では、応力ひずみ計101で計測された生体10の応力Fを自動的に算出部3Bに送り、算出部3Bが応力ひずみ計101から送られた応力Fを用いてヤング率Eを算出したが、本発明はこれには限られない。例えば、応力ひずみ計101で計測された応力Fを、計測者が、入出力装置を操作してコントローラ3に応力Fを入力し、算出部3Bが、入力された応力Fを用いてヤング率Eを算出するようにしてもよい。   In the present embodiment, the stress F of the living body 10 measured by the stress strain meter 101 is automatically sent to the calculation unit 3B, and the calculation unit 3B uses the stress F sent from the stress strain meter 101 to perform the Young's modulus. Although E is calculated, the present invention is not limited to this. For example, the measurer inputs the stress F measured by the stress strain meter 101 to the controller 3 by operating the input / output device, and the calculation unit 3B uses the input stress F to make the Young's modulus E. May be calculated.

実施の形態3.
次に、本発明の実施の形態3について説明する。本実施の形態では、生体10の脂肪組織の比率を測定する。
Embodiment 3 FIG.
Next, a third embodiment of the present invention will be described. In the present embodiment, the ratio of the adipose tissue of the living body 10 is measured.

算出部3Bは、上記実施の形態2と同様に、算出した生体10の密度ρに基づいて、生体10における筋肉組織及び/又は脂肪組織の比率に関する情報を算出する。ここで、筋肉組織の比率をxとし、脂肪組織の比率をyとする。一般に、筋肉組織の密度は1.06に近く、脂肪組織の密度は0.90に近いので、生体の密度ρは、次式で表される。
ρ=1.06x+0.90y … (5)
x+y=1.0…(6)
上記式(5)、(6)より、脂肪組織の比率yは、次式のようになる。
y=(1.06−ρ)/0.16 … (7)
The calculation unit 3B calculates information related to the ratio of muscle tissue and / or adipose tissue in the living body 10 based on the calculated density ρ of the living body 10 as in the second embodiment. Here, the ratio of muscle tissue is x, and the ratio of fat tissue is y. In general, since the density of muscle tissue is close to 1.06 and the density of adipose tissue is close to 0.90, the density ρ of the living body is expressed by the following equation.
ρ = 1.06x + 0.90y (5)
x + y = 1.0 (6)
From the above formulas (5) and (6), the fat tissue ratio y is expressed by the following formula.
y = (1.06-ρ) /0.16 (7)

図7に示すように、探索部3A及び算出部3B等は、上記実施の形態2と同様に、ステップS1〜S8を実行する。ステップS8実行後、算出部3Bは、算出した密度ρに基づいて、式(6)を用いて脂肪組織の比率yを算出する(ステップS9)。   As illustrated in FIG. 7, the search unit 3A, the calculation unit 3B, and the like execute steps S1 to S8 as in the second embodiment. After execution of step S8, the calculation unit 3B calculates the fat tissue ratio y using the equation (6) based on the calculated density ρ (step S9).

以上説明したように、本実施の形態によれば、被爆及び痛みを伴わず、被検者の身体にできるだけ負担をかけることなく、非侵襲かつ簡便な方法で、生体10の脂肪組織の比率yを求めることができる。   As described above, according to the present embodiment, the ratio y of the adipose tissue of the living body 10 is a non-invasive and simple method that does not involve exposure and pain, and places as little burden on the subject's body as possible. Can be requested.

実施の形態4.
次に、本発明の実施の形態4について説明する。本実施の形態では、乳房の硬さを計測する。
Embodiment 4 FIG.
Next, a fourth embodiment of the present invention will be described. In the present embodiment, the hardness of the breast is measured.

本実施の形態では、図8に示すように、以下に示す幾つかの配置パターンで加振器1及び受振器2を配置してそれぞれ計測を行った。なお、加振器1と受振器2との間隔Lは5cmとしている。
パターンA:乳房30の左側で加振器1と受振器2とを鉛直方向に配置する配置パターンである。この場合、振動波の伝達方向は上下方向となる。
パターンB:乳房30の上側で加振器1と受振器2とを水平方向に配置する配置パターンである。この場合、振動波の伝達方向は水平方向となる。
パターンC:乳房30の右側で加振器1と受振器2とを鉛直方向に配置する配置パターンである。この場合、振動波の伝達方向は上下方向となる。
パターンD:乳房30の下側で加振器1と受振器2とを水平方向に配置する配置パターンである。この場合、振動波の伝達方向は水平方向となる。
In the present embodiment, as shown in FIG. 8, the vibration exciter 1 and the vibration receiver 2 are arranged in the following several arrangement patterns, and measurement is performed. In addition, the space | interval L of the vibrator 1 and the geophone 2 is 5 cm.
Pattern A: An arrangement pattern in which the vibrator 1 and the geophone 2 are arranged in the vertical direction on the left side of the breast 30. In this case, the transmission direction of the vibration wave is the vertical direction.
Pattern B: An arrangement pattern in which the vibrator 1 and the geophone 2 are arranged in the horizontal direction above the breast 30. In this case, the transmission direction of the vibration wave is the horizontal direction.
Pattern C: An arrangement pattern in which the vibrator 1 and the geophone 2 are arranged in the vertical direction on the right side of the breast 30. In this case, the transmission direction of the vibration wave is the vertical direction.
Pattern D: An arrangement pattern in which the vibrator 1 and the geophone 2 are arranged in the horizontal direction below the breast 30. In this case, the transmission direction of the vibration wave is the horizontal direction.

図9(A)には、被検者の姿勢を立居とした状態で、パターンAで加振器1を上側とした場合の位相グラフが示されている。また、図9(B)には、パターンAで加振器1を下側とした場合の位相グラフが示されている。図9(A)と、図9(B)を比較するとわかるように、パターンAでは、加振器1及び受振器2を上下反転しただけで、周波数f1の値が大幅に異なるようになった。   FIG. 9A shows a phase graph in the case where the vibrator 1 is on the upper side in the pattern A in a state where the posture of the subject is standing. FIG. 9B shows a phase graph when the vibrator 1 is placed on the lower side in the pattern A. As can be seen by comparing FIG. 9 (A) and FIG. 9 (B), in pattern A, the value of the frequency f1 is significantly different just by turning the vibrator 1 and the receiver 2 upside down. .

図10(A)には、被検者の姿勢を立居とした状態で、パターンCでかつ加振器1を上側とした場合の位相グラフが示されている。また、図10(B)には、パターンCでかつ加振器1を下側とした場合の位相グラフが示されている。図10(A)と、図10(B)を比較するとわかるように、パターンCでは、加振器1及び受振器2を上下反転しただけで、周波数f1の値が大幅に異なるようになった。   FIG. 10A shows a phase graph in the case where the posture of the subject is standing and the pattern C and the vibrator 1 are on the upper side. FIG. 10B shows a phase graph in the case where the pattern C and the vibrator 1 are on the lower side. As can be seen by comparing FIG. 10 (A) and FIG. 10 (B), in the pattern C, the value of the frequency f1 is significantly different just by turning the vibrator 1 and the shaker 2 upside down. .

図11(A)には、被検者の姿勢を立居とした状態で、パターンBでの位相グラフが示されている。パターンBでは、加振器1を左側にしても右側にしても周波数f1の値に変化はなく、その値はほぼ同じで44.4Hzであった。このとき、伝播速度Vは以下のように算出される。
V=0.05×2×44.4=4.4m/s
FIG. 11A shows a phase graph in the pattern B in a state where the posture of the subject is standing. In the pattern B, the value of the frequency f1 was not changed regardless of whether the vibrator 1 was on the left side or the right side, and the value thereof was substantially the same and was 44.4 Hz. At this time, the propagation velocity V is calculated as follows.
V = 0.05 × 2 × 44.4 = 4.4 m / s

図11(B)には、被検者の姿勢を立居とした状態で、パターンDでの位相グラフが示されている。パターンDでは、加振器1を左側にしても右側にしても周波数f1の値に変化はなく、その値はほぼ同じで372.8Hzであった。このとき、伝搬速度Vは以下のように算出される。
V=0.05×2×372.8=37.2m/s
FIG. 11B shows a phase graph of the pattern D in a state where the posture of the subject is standing. In the pattern D, the value of the frequency f1 did not change regardless of whether the vibrator 1 was on the left side or the right side, and the value was almost the same, 372.8 Hz. At this time, the propagation velocity V is calculated as follows.
V = 0.05 × 2 × 372.8 = 37.2 m / s

このように、本実施の形態では、加振器1と受振器2とを水平に配置した方が、精度良く計測を行うことができることが明らかとなった。   Thus, in this Embodiment, it became clear that the direction which arrange | positioned the vibrator 1 and the geophone 2 horizontally can perform a measurement with sufficient precision.

ところで、同じ被検者が寝た状態でパターンBでの計測を行った場合には、周波数f1は、74.0Hzとなり、振動波の伝播速度Vは7.4m/sとなった。一方、被検者が寝た状態でパターンDの計測を行った場合には、周波数f1は、134.6Hzとなり、伝搬速度Vは13.5m/sとなった。   By the way, when measurement was performed with the pattern B while the same subject was sleeping, the frequency f1 was 74.0 Hz, and the propagation velocity V of the vibration wave was 7.4 m / s. On the other hand, when the pattern D was measured with the subject sleeping, the frequency f1 was 134.6 Hz, and the propagation velocity V was 13.5 m / s.

以上の結果から同じ被検者で同じ箇所を計測しても、立った状態と寝た状態とで、それぞれ振動波の伝搬速度Vが異なることがわかる。具体的には、パターンBの計測位置、すなわち乳房30の上半分は、立った方が寝た姿勢よりも柔らかくなり、パターンDの計測位置、すなわち乳房30の下側では、寝た姿勢の方が、立った状態よりも柔らかくなった。被検者の姿勢に応じた乳房30の各部位の硬さの変化等は、その乳房30に装着される下着の形状などを決める際に有用な情報となる。例えば、立ったとき、座ったときの乳房30の各部位の硬さを測定し、それぞれの姿勢に応じて、ほどよい締め付け力を有する下着をオーダーメイドでつくることも可能となる。   From the above results, it can be seen that even if the same location is measured by the same subject, the propagation speed V of the vibration wave differs between the standing state and the sleeping state. Specifically, the measurement position of the pattern B, that is, the upper half of the breast 30, becomes softer than the posture of sleeping, and the measurement position of the pattern D, that is, the lower side of the breast 30, the sleeping posture. However, it became softer than standing. Changes in the hardness of each part of the breast 30 according to the posture of the subject are useful information for determining the shape of the underwear to be worn on the breast 30. For example, it is possible to measure the hardness of each part of the breast 30 when standing or sitting, and to make an underwear having an appropriate tightening force according to each posture.

この他、本実施の形態に係る生体計測装置100は、乳がんの検出にも適用することができる。乳がんの腫瘍は、正常な組織に比べ、10倍以上硬いことが知られている。したがって、生体計測装置100を用いて、乳房30の同じ部位における振動波の伝搬速度Vの変化を観察する。例えば、その部位で、ある時点で、振動波の伝搬速度Vが速くなった場合には、その部位の細胞組織が硬くなったことを示しているので、乳がんの発生を疑うことができる。このような計測を行うことにより、マンモグラフィ検査等を行うことなく、乳がんを早期発見することも可能となる。   In addition, the biological measurement apparatus 100 according to the present embodiment can be applied to detection of breast cancer. Breast cancer tumors are known to be ten times harder than normal tissues. Therefore, the living body measuring apparatus 100 is used to observe a change in the propagation speed V of the vibration wave in the same part of the breast 30. For example, if the propagation velocity V of the vibration wave increases at a certain point in that region, it indicates that the cell tissue at that region has become hard, and thus the occurrence of breast cancer can be suspected. By performing such measurement, it becomes possible to detect breast cancer early without performing mammography examination or the like.

実施の形態5.
次に、本発明の実施の形態5について説明する。ここでは、上記各実施の形態に係る生体計測装置100を用いて非破壊で生体10の硬さを測定する方法の原理について改めてまとめるとともに、加振器1及び受振器2を用いた振動の伝達方法のバリエーションについて説明する。本実施の形態に係る生体計測装置100は、加振器1を2つ備える。なお、本実施の形態では、加振器1は、縦波だけでなく横波も生体10に伝える。
Embodiment 5. FIG.
Next, a fifth embodiment of the present invention will be described. Here, the principle of the method for measuring the hardness of the living body 10 in a non-destructive manner using the living body measuring apparatus 100 according to each of the above embodiments is summarized, and vibration transmission using the vibrator 1 and the shaker 2 is performed. A method variation will be described. The biological measurement apparatus 100 according to the present embodiment includes two vibrators 1. In the present embodiment, the vibrator 1 transmits not only a longitudinal wave but also a transverse wave to the living body 10.

一般に、波は媒質の硬いものでは速く伝播し、軟らかいものでは遅く伝播する性質がある。これは、波の伝播速度V[m/s]と、媒質の硬さ(ヤング率)E、密度ρの間で、以下の式(1)’のような関係があるからである。
V=C×√(E/ρ) …(1)’
ここで、係数Cは、縦波と横波の場合でポアソン比によってのみ表される式である。ポアソン比は、生体10であれば0.3から0.5の間で一定の値をとるので、係数Cはほぼ定数になる。また、密度ρについても、生体10であれば水の密度に近い定数をとる。このように係数Cと密度ρは定数であるので、式(1)’から、伝播速度Vを測定すれば、硬さEを定量的に評価することができる。
In general, a wave propagates fast when the medium is hard, and propagates slowly when the medium is soft. This is because there is a relationship represented by the following equation (1) ′ between the wave propagation velocity V [m / s], the medium hardness (Young's modulus) E, and the density ρ.
V = C × √ (E / ρ) (1) ′
Here, the coefficient C is an expression expressed only by the Poisson's ratio in the case of a longitudinal wave and a transverse wave. Since the Poisson's ratio takes a constant value between 0.3 and 0.5 in the case of the living body 10, the coefficient C is almost constant. Further, the density ρ is a constant close to the density of water in the case of the living body 10. Since the coefficient C and the density ρ are constants as described above, the hardness E can be quantitatively evaluated by measuring the propagation velocity V from the equation (1) ′.

振動の伝播速度Vを評価する方法には、振動の伝播の距離lと時間tを直接測定して評価する方法と、加振点と受振点の位相差θを用いて評価する方法の2つがある。伝播の距離lと時間tを直接測る方法では、加振点でパルス状の振動を一定の周期で与え、加振点から距離lの位置に置かれた受振点までパルス波が到達する時間tを測定する。しかし、この方法では、超音波などを用いたパルス波発生装置や、コンピュータによる振動制御部と受振部の製作費用が高額になる。   There are two methods for evaluating the vibration propagation velocity V: a method of directly measuring and evaluating the vibration propagation distance l and time t, and a method of evaluating using the phase difference θ between the excitation point and the receiving point. is there. In the method of directly measuring the propagation distance l and the time t, a pulsed vibration is given at a certain period at the excitation point, and the time t when the pulse wave reaches the receiving point located at a distance l from the excitation point. Measure. However, in this method, the manufacturing cost of a pulse wave generator using ultrasonic waves or the like, or a vibration control unit and a vibration receiving unit by a computer becomes high.

一方、位相差θを用いて評価する方法では、次の式(2)’を用いて伝播速度V[m/s]を求める。
V=360×f×l/θ …(2)’
ここで、f[Hz]は振動数であり、l[m]は加振点と受振点の間の距離である。また、θ[°]は、加振器1の加振点での振動の位相から受振器2の受振点の振動の位相を引いた位相差である。ここで、例えば、加振点と受振点の位相差が180度である場合には、図12(A)に示すように加振点から半波長分だけ離れた受振点での振動が測定され、加振点と受振点の位相差が360度の場合には、図12(B)に示すように、加振点から1波長分だけ離れた受振点での振動が測定される。
On the other hand, in the evaluation method using the phase difference θ, the propagation velocity V [m / s] is obtained using the following equation (2) ′.
V = 360 × f × l / θ (2) ′
Here, f [Hz] is the frequency, and l [m] is the distance between the excitation point and the receiving point. Θ [°] is a phase difference obtained by subtracting the vibration phase at the receiving point of the shaker 2 from the vibration phase at the excitation point of the shaker 1. Here, for example, when the phase difference between the excitation point and the receiving point is 180 degrees, the vibration at the receiving point separated by a half wavelength from the excitation point is measured as shown in FIG. When the phase difference between the excitation point and the receiving point is 360 degrees, as shown in FIG. 12B, the vibration at the receiving point that is separated from the excitation point by one wavelength is measured.

このように、加振点での振動と受振点で受振される振動との位相差θを用いる方法では、可聴域内(20〜7000Hz)の正弦波の機械振動を生体10に与えるだけよいため、装置の構造を簡素にすることができるので、装置の故障を発生しにくくすることができるうえ、装置の製造費用を抑えることができる。また、加振器1の加振点から継続的に正弦波を与えるため、フーリエ変換(FFT)を安定的に実行でき、位相差θの検出精度を向上することができる。そこで、本実施の形態に係る生体計測装置100は、一定の振動数fと距離lの下で、位相差θを測定する。   Thus, in the method using the phase difference θ between the vibration at the excitation point and the vibration received at the receiving point, it is only necessary to give the living body 10 a sine wave mechanical vibration within the audible range (20 to 7000 Hz). Since the structure of the apparatus can be simplified, it is possible to make it difficult for the apparatus to fail and to reduce the manufacturing cost of the apparatus. Further, since the sine wave is continuously given from the excitation point of the vibrator 1, the Fourier transform (FFT) can be stably executed, and the detection accuracy of the phase difference θ can be improved. Therefore, living body measuring apparatus 100 according to the present embodiment measures phase difference θ under a constant frequency f and distance l.

受振感度の向上のため、本実施の形態に係る生体計測装置100では、生体10の特徴に応じて以下に示す3つの方法のいずれかを用いて、加受振間の振動の伝播速度を測定し硬さを評価する。なお、生体計測装置100は、2台の加振器1を備えている。   In order to improve the vibration receiving sensitivity, the living body measurement apparatus 100 according to the present embodiment measures the propagation speed of vibration between vibrations using one of the following three methods according to the characteristics of the living body 10. Evaluate hardness. The biological measurement apparatus 100 includes two vibrators 1.

(1)1点加振法
図13に示すように、加振器1を1台用いて1点で加振し受振器2により1点で受振する方法である。この方法は、生体10の硬さが位置によって変化する場合、すなわち硬さが不均一な場合に、測定箇所ごとの硬さを測定できる。
(1) One-Point Excitation Method As shown in FIG. 13, this is a method in which one exciter 1 is used to excite at one point and the receiver 2 receives the signal at one point. This method can measure the hardness of each measurement location when the hardness of the living body 10 changes depending on the position, that is, when the hardness is non-uniform.

(2)2点同位相加振法
図14に示すように、加振器1を2台用いて異なる2つの加振点を同位相で加振し、2つの加振器1の中央1点を受振点として受振器により受振する方法である。2つの加振点を同位相で加振すると、加振点で発生した横波は、波動の重ね合わせの原理により受振点での横波の振幅が大きくなり受振が容易になる。ただし、加振点で発生した縦波は受振点で相殺されるので、この方法で縦波を利用した伝播速度の測定には向かない。また、2つの加振点間の硬さが均一な場合では受振が容易であるが、不均一な場合には加振点から異なる伝播速度で受振点に振動波が到達するために、受振点での位相がそろわず伝播速度Vの測定精度が低下する。
(2) Two-point in-phase excitation method As shown in FIG. 14, two different excitation points are excited in the same phase using two excitation units 1, and one center point of the two excitation units 1. Is received by the geophone using the geophone as a receiving point. When the two excitation points are excited in the same phase, the transverse wave generated at the excitation point becomes easy to receive due to the increased amplitude of the transverse wave at the receiving point due to the principle of wave superposition. However, since the longitudinal wave generated at the excitation point is canceled out at the receiving point, this method is not suitable for measuring the propagation velocity using the longitudinal wave. In addition, when the hardness between the two excitation points is uniform, the vibration is easily received. When the hardness is not uniform, the vibration waves reach the reception point at different propagation velocities from the excitation point. As a result, the measurement accuracy of the propagation velocity V decreases.

(3)2点逆位相加振法
図14に示すように、加振器1を2台用いて異なる2つの加振点を逆位相で加振し、2つの加振器1の中央1点を受振点として受振器により受振する受振する方法である。2つの加振点を逆位相で加振すると、加振点で発生した縦波は、受振点での縦波の振幅が大きくなり受振が容易になる。一方、加振点で発生した横波は受振点で相殺される。受振が容易な縦波を用いると、生体10の内部で硬さの異なる組織が受振器の直下付近に局所的に存在する場合、振動の一部が硬さの異なる組織で反射し伝播速度Vが変化する。したがって、加振点及び受振点の位置を変えながら伝播速度Vの変化を見ることにより、硬さの異なる組織を見つけることができる。また、2つの加振点間の硬さが不均一な場合、(2)と同様の理由で伝播速度の測定精度が悪化する。
(3) Two-point anti-phase excitation method As shown in FIG. 14, two different excitation points are excited in opposite phases using two vibrators 1, and the center 1 of the two vibrators 1 is In this method, the point is received by a receiving device with the receiving point as a receiving point. When the two excitation points are excited in opposite phases, the longitudinal wave generated at the excitation point increases the amplitude of the longitudinal wave at the reception point and is easy to receive. On the other hand, the transverse wave generated at the excitation point is canceled at the receiving point. When longitudinal waves that are easy to receive vibrations are used, when a tissue having different hardness is locally present in the living body 10 in the vicinity of directly under the vibration receiver, a part of the vibration is reflected by the tissue having different hardness and the propagation velocity V Changes. Therefore, it is possible to find tissues having different hardnesses by observing changes in the propagation velocity V while changing the positions of the excitation point and the receiving point. Further, when the hardness between the two excitation points is not uniform, the measurement accuracy of the propagation speed is deteriorated for the same reason as in (2).

以上、(1)、(2)、(3)の方法についての特徴をまとめると、以下の表に示すようになる。このように、1点加振法、2点同位相加振法、2点逆位相加振法は、生体10の状態に応じて使い分けるのが望ましい。生体10の状態が不明な場合には、1点加振法、2点同位相加振法、2点逆位相加振法全てで測定を行うのが望ましい。

Figure 2017047193
The characteristics of the methods (1), (2), and (3) are summarized as shown in the following table. As described above, it is desirable that the one-point excitation method, the two-point in-phase excitation method, and the two-point anti-phase excitation method are properly used according to the state of the living body 10. When the state of the living body 10 is unknown, it is desirable to perform measurement by all of the one-point excitation method, the two-point in-phase excitation method, and the two-point anti-phase excitation method.
Figure 2017047193

生体10の組織に見立てた寒天の伝播速度Vを、上記の3つの方法を用いて振動の測定をした。寒天は、縦7cm、横16cm、高さ4cmの直方体として作成した。   The vibration speed of the agar propagation velocity V assumed to be the tissue of the living body 10 was measured using the above three methods. The agar was prepared as a rectangular parallelepiped 7 cm long, 16 cm wide and 4 cm high.

なお、1点加振法では、加振点と受振点とを、直方体上面の中央から端にそれぞれ5cm離れた点に置いた。また、2点同位相加振法と2点逆位相加振法では、受振点を直方体上面の中央に置き、2つの加振点をそれぞれ中央から2.5cm離れた点に置いた。そして、1000Hzまでの振動を加振点で寒天に与えた。1点加振法、2点同位相加振法、2点逆位相加振法の3方法により、寒天に振動を与えて測定した結果を図15、図16、図17に示す。   In the one-point excitation method, the excitation point and the reception point were placed at points 5 cm away from the center to the end of the upper surface of the rectangular parallelepiped. In the two-point in-phase excitation method and the two-point anti-phase excitation method, the receiving point was placed at the center of the rectangular parallelepiped upper surface, and the two excitation points were placed at points 2.5 cm away from the center. And vibration to 1000 Hz was given to the agar at the excitation point. 15, 16, and 17 show the results obtained by applying vibration to the agar by three methods of the one-point excitation method, the two-point in-phase excitation method, and the two-point anti-phase excitation method.

2点同位相・逆位相加振法のときの加振と受振の距離2.5cmは、1点加振法の時の半分なので、式(2)より位相差も半分になることが予想される。実際、0[Hz]から1000[Hz]までの間で測定された2点同位相・逆位相加振法の位相差(図16と図17参照)は、1点加振法の位相差(図15参照)のほぼ半分になっていることがわかった。また、2点同位相加振法では横波を、2点逆位相加振法では縦波を受振点で測定する。両者とも加振と受振の距離が同じであれば、式(2)より位相差は同じであることが予想される。実際、測定された2点同位相加振法の位相差(図16参照)は2点逆位相加振法の位相差(図17参照)と同じであることがわかった。   The distance of 2.5 cm between excitation and reception in the two-point in-phase / anti-phase excitation method is half that in the one-point excitation method, so the phase difference is expected to be halved from equation (2). Is done. Actually, the phase difference of the two-point in-phase / anti-phase excitation method (see FIGS. 16 and 17) measured between 0 [Hz] and 1000 [Hz] is the phase difference of the one-point excitation method. It was found that it was almost half of that (see FIG. 15). In the two-point in-phase excitation method, a transverse wave is measured at the receiving point, and in the two-point anti-phase excitation method, a longitudinal wave is measured at the receiving point. In both cases, if the distance between excitation and reception is the same, the phase difference is expected to be the same from equation (2). Actually, it was found that the measured phase difference of the two-point in-phase excitation method (see FIG. 16) was the same as the phase difference of the two-point anti-phase excitation method (see FIG. 17).

以上説明したように、本実施の形態によれば、被爆及び痛みを伴わず、被検者に身体にできるだけ負担をかけることなく、非侵襲かつ簡便な方法で、乳房30の各部位の硬さを求めることができる。   As described above, according to the present embodiment, the hardness of each part of the breast 30 can be obtained in a non-invasive and simple manner without causing exposure and pain, and placing as little burden on the subject as possible. Can be requested.

なお、上記各実施の形態では、コンピュータが、位相差が180度になるような周波数を探索したが、本発明はこれには限られない。例えば、位相差が90度、270度になるような周波数を探索するようにしてもよい。   In each of the above embodiments, the computer searches for a frequency at which the phase difference is 180 degrees, but the present invention is not limited to this. For example, a frequency where the phase difference is 90 degrees or 270 degrees may be searched.

この発明は、この発明の広義の精神と範囲を逸脱することなく、様々な実施の形態及び変形が可能とされるものである。また、上述した実施の形態は、この発明を説明するためのものであり、この発明の範囲を限定するものではない。すなわち、この発明の範囲は、実施の形態ではなく、特許請求の範囲によって示される。そして、特許請求の範囲内及びそれと同等の発明の意義の範囲内で施される様々な変形が、この発明の範囲内とみなされる。   Various embodiments and modifications can be made to the present invention without departing from the broad spirit and scope of the present invention. The above-described embodiments are for explaining the present invention and do not limit the scope of the present invention. In other words, the scope of the present invention is shown not by the embodiments but by the claims. Various modifications within the scope of the claims and within the scope of the equivalent invention are considered to be within the scope of the present invention.

本発明は、生体内の情報を計測するのに適用することができる。   The present invention can be applied to measuring in vivo information.

1 加振器、1A 加速度ピックアップ、2 受振器、3 コントローラ、3A 探索部、3B 算出部、10 生体、20 棒状体、21 応力計測部、22 押し当て部、30 乳房、100 生体計測装置、101 応力ひずみ計。   DESCRIPTION OF SYMBOLS 1 Exciter, 1A Acceleration pick-up, 2 Geophone, 3 controller, 3A search part, 3B calculation part, 10 biological body, 20 stick-shaped body, 21 stress measurement part, 22 pressing part, 30 breast, 100 biological measurement apparatus, 101 Stress strain gauge.

Claims (13)

可聴域の範囲から選択された単一の周波数の振動波を生体に与える加振部と、
前記加振部が前記生体に前記振動波を与える位置から所定距離離れた位置で、前記生体を伝わる振動波を受振する受振部と、
前記加振部によって前記生体に与えられる振動波の周波数を変更しながら、前記加振部により前記生体に与えられる振動波と前記受振部で受振される振動波との位相差が特定の位相差となるような前記振動波の周波数を探索する探索部と、
前記探索部で探索された周波数に基づいて、前記生体中における前記振動波の伝搬速度を算出する算出部と、
を備える生体計測装置。
An excitation unit that gives a living body a vibration wave having a single frequency selected from the range of the audible range;
A vibration receiving unit that receives vibration waves transmitted through the living body at a position away from a position where the vibration unit gives the vibration wave to the living body by a predetermined distance;
While changing the frequency of the vibration wave given to the living body by the vibration unit, the phase difference between the vibration wave given to the living body by the vibration unit and the vibration wave received by the vibration receiving unit is a specific phase difference. A search unit for searching for the frequency of the vibration wave such that
A calculation unit that calculates a propagation speed of the vibration wave in the living body based on the frequency searched by the search unit;
A biological measurement apparatus comprising:
球面のプローブを前記生体に所定の深さ押し込んだ状態で計測される応力を計測する応力計測部を備え、
前記算出部は、
前記応力計測部で計測された応力に基づいて、前記生体のヤング率を算出する、
請求項1に記載の生体計測装置。
A stress measuring unit that measures stress measured in a state where a spherical probe is pushed into the living body to a predetermined depth;
The calculation unit includes:
Based on the stress measured by the stress measurement unit, to calculate the Young's modulus of the living body,
The biological measurement apparatus according to claim 1.
前記算出部は、
算出した前記伝搬速度と、前記生体のヤング率とに基づいて、前記生体の密度を算出する、
請求項2に記載の生体計測装置。
The calculation unit includes:
Calculating the density of the living body based on the calculated propagation velocity and the Young's modulus of the living body;
The biological measurement apparatus according to claim 2.
前記算出部は、
算出した前記生体の密度に基づいて、前記生体における筋肉組織及び/又は脂肪組織の比率に関する情報を算出する、
請求項3に記載の生体計測装置。
The calculation unit includes:
Based on the calculated density of the living body, calculate information on the ratio of muscle tissue and / or adipose tissue in the living body,
The biological measurement apparatus according to claim 3.
前記加振器は、
前記生体における異なる2つの加振点を前記各加振器でそれぞれ同位相又は逆位相で加振し、
前記受振器は、
前記2つの加振点の中央1点である受振点で受振する、
請求項1に記載の生体計測装置。
The vibrator is
Two different excitation points in the living body are respectively excited in the same phase or in opposite phase by the respective vibrators,
The geophone is
Receiving at a receiving point which is one central point of the two excitation points;
The biological measurement apparatus according to claim 1.
コンピュータが、加振器によって生体に与えられる振動波の周波数を可聴域の範囲で変更しながら、前記加振器により前記生体に与えられる振動波と前記受振器で受振される振動波との位相差が特定の位相差となるような前記振動波の周波数を探索する探索ステップと、
コンピュータが、前記探索ステップで探索された周波数に基づいて、前記生体中における前記振動波の伝搬速度を前記生体の硬さに関する情報として算出する算出ステップと、
を含む生体計測方法。
While the computer changes the frequency of the vibration wave applied to the living body by the vibrator in the range of the audible range, the vibration wave applied to the living body by the vibrator and the vibration wave received by the shaker A search step for searching for a frequency of the vibration wave such that the phase difference becomes a specific phase difference;
A calculation step in which the computer calculates the propagation speed of the vibration wave in the living body as information on the hardness of the living body based on the frequency searched in the searching step;
A biological measurement method including:
球面型のプローブを前記生体の組織を所定の深さまで押し込んだ時に得られる応力を計測する応力ひずみ計を用いて、当該応力を計測する計測ステップを含み、
前記算出ステップでは、
前記コンピュータが、
前記計測ステップで計測された応力に基づいて、前記生体のヤング率を算出する、
請求項6に記載の生体計測方法。
A measurement step of measuring the stress using a stress strain meter that measures the stress obtained when the spherical tissue probe is pushed into the living tissue to a predetermined depth;
In the calculating step,
The computer is
Based on the stress measured in the measurement step, calculate the Young's modulus of the living body,
The biological measurement method according to claim 6.
前記算出ステップでは、
コンピュータが、前記伝搬速度と前記生体のヤング率とに基づいて、前記生体の密度を前記生体の硬さに関する情報として算出する、
請求項7に記載の生体計測方法。
In the calculating step,
The computer calculates the density of the living body as information on the hardness of the living body based on the propagation speed and the Young's modulus of the living body.
The biological measurement method according to claim 7.
前記算出ステップでは、
コンピュータは、算出した前記生体の密度に基づいて、前記生体における筋肉組織及び/又は脂肪組織の比率に関する情報を前記生体の硬さに関する情報として算出する、
請求項8に記載の生体計測方法。
In the calculating step,
Based on the calculated density of the living body, the computer calculates information on the ratio of muscle tissue and / or adipose tissue in the living body as information on the hardness of the living body,
The biological measurement method according to claim 8.
前記探索ステップでは、
前記加振器と前記受振器とを水平に配置して前記振動波の周波数を探索する、
請求項6から9のいずれか一項に記載の生体計測方法。
In the search step,
Search the frequency of the vibration wave by horizontally arranging the vibrator and the geophone,
The biological measurement method according to any one of claims 6 to 9.
前記探索ステップでは、
乳房の上側部分に前記加振器と前記受振器とを配置して前記探索ステップ及び前記算出ステップを行って、その部分における前記振動波の伝搬速度を計測し、
乳房の下側部分に前記加振器と前記受振器とを配置して前記探索ステップ及び前記算出ステップを行って、その部分における前記振動波の伝搬速度を計測する、
請求項10に記載の生体計測方法。
In the search step,
Place the exciter and the geophone on the upper part of the breast, perform the search step and the calculation step, measure the propagation speed of the vibration wave in that part,
Placing the exciter and the geophone in the lower part of the breast and performing the search step and the calculation step to measure the propagation speed of the vibration wave in that part;
The biological measurement method according to claim 10.
前記探索ステップでは、
前記生体における異なる複数の加振点を前記加振器でそれぞれ同位相又は逆位相で加振し、
前記複数の加振点の中央1点を受振点として前記受振器で受振する、
請求項6に記載の生体計測方法。
In the search step,
A plurality of different excitation points in the living body are respectively excited in the same phase or in opposite phase by the vibrator,
Receiving at the geophone as a receiving point at the center of the plurality of excitation points,
The biological measurement method according to claim 6.
乳房を計測対象として、請求項6から12のいずれか一項に記載の生体計測方法を行って、前記生体中における前記振動波の伝搬速度を求め、その伝搬速度の変化に基づいて、乳がんを検出する、
検出方法。
The biological measurement method according to any one of claims 6 to 12 is performed on a breast as a measurement target, a propagation speed of the vibration wave in the living body is obtained, and breast cancer is determined based on a change in the propagation speed. To detect,
Detection method.
JP2016162687A 2015-09-02 2016-08-23 Biological measurement apparatus, biological measurement method, and detection method Pending JP2017047193A (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2015172707 2015-09-02
JP2015172707 2015-09-02

Publications (1)

Publication Number Publication Date
JP2017047193A true JP2017047193A (en) 2017-03-09

Family

ID=58280592

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016162687A Pending JP2017047193A (en) 2015-09-02 2016-08-23 Biological measurement apparatus, biological measurement method, and detection method

Country Status (1)

Country Link
JP (1) JP2017047193A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017169332A1 (en) * 2016-03-28 2017-10-05 株式会社日立製作所 Ultrasonic imaging device and ultrasonic reception signal processing method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017169332A1 (en) * 2016-03-28 2017-10-05 株式会社日立製作所 Ultrasonic imaging device and ultrasonic reception signal processing method
JPWO2017169332A1 (en) * 2016-03-28 2018-06-28 株式会社日立製作所 Ultrasonic imaging apparatus and ultrasonic reception signal processing method

Similar Documents

Publication Publication Date Title
KR101914021B1 (en) Method and Apparatus for Detecting Viscoelastic Parameters of Viscoelastic Media
Catheline et al. A solution to diffraction biases in sonoelasticity: the acoustic impulse technique
EP2310829B1 (en) A system and method for detection, characterization and imaging of heterogeneity using shear wave induced resonance
JP6437525B2 (en) Measurement method and configuration using electromagnetic waves
US9585631B2 (en) Devices, methods, and systems for measuring elastic properties of biological tissues using acoustic force
RU2688294C1 (en) Method and device for determining elasticity
US20130018240A1 (en) Body measurement and imaging with a mobile device
CN106233134B (en) Apparatus for ultrasonic examination and ultrasonic inspection method
CN106232013B (en) Ultrasonic method and apparatus for characterizing a weakly anisotropic soft medium, and ultrasonic probe assembly for use in the characterizing apparatus
BRPI1010214B1 (en) shear wave dispersion vibrometry (sduv), non-transitory computer-readable media and ultrasound device configured for shear wave dispersion vibrometry (sduv)
WO2015077086A2 (en) System and method for analyzing tissue using shear waves
JP2018531138A6 (en) Elasticity detection method and apparatus
CN111631687B (en) Intelligent spine detection system, detector and detection method
US20170055836A1 (en) Apparatus for Dynamic Stress Measurement
Laugier Quantitative ultrasound instrumentation for bone in vivo characterization
Chen et al. Shear wave dispersion ultrasound vibrometry based on a different mechanical model for soft tissue characterization
JP2017047193A (en) Biological measurement apparatus, biological measurement method, and detection method
JP2021013575A (en) Ultrasonic diagnostic apparatus and control method of ultrasonic diagnostic apparatus
WO2021212578A1 (en) Viscoelasticity measuring method, elasticity measuring method and ultrasonic measuring system
Lin et al. Viscoelastic properties of normal rat liver measured by ultrasound elastography: Comparison with oscillatory rheometry
JP6633854B2 (en) Ultrasonic measuring device
Vogl et al. Reliability of phase velocity measurements of flexural acoustic waves in the human tibia in-vivo
JP6065601B2 (en) Ultrasonic diagnostic equipment
JPH1176176A (en) Organic tissue viscoelasticity characteristic measuring apparatus and method and medium storing organic tissue viscoelasticity characteristic calculating program
JP2006230861A (en) Muscle fatigue grade measurement device