Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP4688466B2 - Ultrasonic diagnostic equipment - Google Patents
[go: Go Back, main page]

JP4688466B2 - Ultrasonic diagnostic equipment - Google Patents

Ultrasonic diagnostic equipment Download PDF

Info

Publication number
JP4688466B2
JP4688466B2 JP2004295899A JP2004295899A JP4688466B2 JP 4688466 B2 JP4688466 B2 JP 4688466B2 JP 2004295899 A JP2004295899 A JP 2004295899A JP 2004295899 A JP2004295899 A JP 2004295899A JP 4688466 B2 JP4688466 B2 JP 4688466B2
Authority
JP
Japan
Prior art keywords
signal
wave
doppler
diagnostic apparatus
ultrasonic diagnostic
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.)
Expired - Fee Related
Application number
JP2004295899A
Other languages
Japanese (ja)
Other versions
JP2005253949A (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.)
Hitachi Ltd
Original Assignee
Aloka Co Ltd
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 Aloka Co Ltd filed Critical Aloka Co Ltd
Priority to JP2004295899A priority Critical patent/JP4688466B2/en
Publication of JP2005253949A publication Critical patent/JP2005253949A/en
Application granted granted Critical
Publication of JP4688466B2 publication Critical patent/JP4688466B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Description

本発明は、超音波診断装置に関し、特に連続波(CW)ドプラ法における送受信技術に関する。   The present invention relates to an ultrasonic diagnostic apparatus, and more particularly to a transmission / reception technique in a continuous wave (CW) Doppler method.

超音波診断装置において、連続波ドプラモードでは、数MHzの正弦波として構成される送信波が生体内へ連続的に放射され、生体内からの反射波が連続的に受波される。反射波には、生体内における運動体(例えば血流)によるドプラシフト情報が含まれる。そこで、その情報を取り出して周波数解析することによりドプラ波形が得られる。反射波には、生体内の強反射体(主に静止体)からの情報(クラッタ)も含まれる。そのクラッタをできる限り抑圧してドプラ情報だけを抽出することが望まれるが、特に低速の観測域においてはクラッタを十分に除去できないという問題がある。   In the ultrasonic diagnostic apparatus, in the continuous wave Doppler mode, a transmission wave configured as a sine wave of several MHz is continuously radiated into the living body, and a reflected wave from the living body is continuously received. The reflected wave includes Doppler shift information due to a moving body (for example, blood flow) in the living body. Therefore, a Doppler waveform can be obtained by extracting the information and performing frequency analysis. The reflected wave includes information (clutter) from a strong reflector (mainly a stationary body) in the living body. Although it is desired to suppress the clutter as much as possible and extract only Doppler information, there is a problem that the clutter cannot be sufficiently removed particularly in a low-speed observation area.

特開平3−48789号公報JP-A-3-48789

以上のように、ドプラ情報を観測する場合においてはクラッタが妨害となり、それをできる限り抑圧することが望まれている。   As described above, when observing Doppler information, clutter is an obstacle, and it is desired to suppress it as much as possible.

特許文献1には、周波数変調された鋸歯状の連続波を生体内へ送波し、反射波を受波して得られた受信信号と変調された送信波とからビート信号を発生させ、それを周波数解析する装置が開示されている。しかし、この文献に記載された装置は、変調波の高調波帯域に着目するものではない。   In Patent Document 1, a frequency-modulated sawtooth continuous wave is transmitted into a living body, and a beat signal is generated from a reception signal obtained by receiving a reflected wave and a modulated transmission wave. Has been disclosed. However, the device described in this document does not focus on the harmonic band of the modulated wave.

本発明の目的は、連続波ドプラ法の適用に当たって、クラッタをできる限り抑圧しつつドプラ情報を感度良く観測できるようにすることにある。   An object of the present invention is to enable observation of Doppler information with high sensitivity while suppressing clutter as much as possible when applying the continuous wave Doppler method.

上記目的を達成するために、本発明の好適な態様である超音波診断装置は、搬送波信号を変調波信号を用いて変調処理することにより、変調送信信号を生成する送信部と、前記変調送信信号の供給によって生体に対して超音波を送波し、生体からの反射波を受波して受信信号を出力する送受波器と、前記受信信号を前記変調送信信号を用いて復調し、これにより復調信号を得る復調手段と、前記復調信号から、前記変調波信号を基本波とした場合における基本波成分および第n(nは2以上の自然数)高調波成分のうちの少なくとも一つの成分の付近に存在するドプラ情報を抽出するドプラ情報抽出手段と、を有することを特徴とする。   In order to achieve the above object, an ultrasonic diagnostic apparatus according to a preferred aspect of the present invention includes a transmitter that generates a modulated transmission signal by modulating a carrier wave signal using a modulated wave signal, and the modulated transmission. A transmitter / receiver for transmitting an ultrasonic wave to a living body by supplying a signal, receiving a reflected wave from the living body and outputting a received signal; and demodulating the received signal using the modulated transmission signal; And a demodulating means for obtaining a demodulated signal from the demodulated signal and at least one of a fundamental wave component and an nth (n is a natural number of 2 or more) harmonic component when the modulated wave signal is a fundamental wave. And Doppler information extraction means for extracting Doppler information existing in the vicinity.

上記構成では、搬送波信号を変調波信号を用いて変調処理することにより変調送信信号を生成し、さらに、この変調送信信号を用いて受信信号を復調する。このため、無変調の搬送波信号の場合に搬送波信号の周波数付近に局在していた本来取り除くべき固定物からのエコー(クラッタ)が、変調波信号を基本波とした場合における基本波成分および高調波成分に分散される。分散された結果、基本波成分および各高調波成分付近でのクラッタは、無変調の場合のクラッタに比べて数十dB小さくなる。上記構成では、基本波成分および各高調波成分のうちの少なくとも一つの成分付近のドプラ情報を抽出するため、無変調の場合に局在していたクラッタの影響が低減(全く影響を受けないような構成がさらに望ましい)される。結果として、上記構成によればクラッタをできる限り抑圧しつつドプラ情報を感度良く観測できる。   In the above configuration, the modulated transmission signal is generated by modulating the carrier wave signal using the modulated wave signal, and the received signal is demodulated using the modulated transmission signal. For this reason, in the case of an unmodulated carrier signal, an echo (clutter) from a fixed object to be removed that was localized in the vicinity of the frequency of the carrier signal is a fundamental wave component and a harmonic when the modulated wave signal is a fundamental wave. Dispersed into wave components. As a result of the dispersion, the clutter in the vicinity of the fundamental wave component and each harmonic component becomes tens of dB smaller than the clutter in the case of no modulation. In the above configuration, Doppler information in the vicinity of at least one of the fundamental wave component and each harmonic component is extracted, so the influence of clutter localized in the case of no modulation is reduced (not affected at all). Is more desirable). As a result, according to the above configuration, Doppler information can be observed with high sensitivity while suppressing clutter as much as possible.

望ましくは、前記復調手段は、前記変調送信信号を参照信号として前記受信信号を検波処理する検波回路を含むことを特徴とする。なお、検波回路として、受信信号を直交検波処理する直交検波回路を利用してもよい。また望ましくは、前記ドプラ情報抽出手段は、前記少なくとも一つの成分およびその付近に存在するドプラ信号から成る側帯波を抽出するフィルタ部を含む、ことを特徴とする。さらに望ましくは、前記ドプラ情報抽出手段は、前記側帯波から前記ドプラ信号を抽出するドプラ信号抽出部を含む、ことを特徴とする。さらに望ましくは、前記ドプラ信号抽出部は、前記変調波信号の周波数をN倍(Nは自然数)した信号のうちの少なくとも一つを参照信号として、前記側帯波を検波処理する検波回路を含む、ことを特徴とする。   Preferably, the demodulation means includes a detection circuit for detecting the received signal using the modulated transmission signal as a reference signal. Note that a quadrature detection circuit that performs quadrature detection processing on the received signal may be used as the detection circuit. Preferably, the Doppler information extraction unit includes a filter unit that extracts a sideband composed of the at least one component and a Doppler signal existing in the vicinity thereof. More preferably, the Doppler information extraction unit includes a Doppler signal extraction unit that extracts the Doppler signal from the sideband. More preferably, the Doppler signal extraction unit includes a detection circuit that detects the sideband wave using at least one of signals obtained by multiplying the frequency of the modulated wave signal by N (N is a natural number) as a reference signal. It is characterized by that.

望ましくは、前記復調手段に与える前記変調波送信信号に対して遅延処理を施す遅延回路を含む、ことを特徴とする。この構成によれば、例えば、遅延時間を調整して各遅延時間ごとにドプラ信号/クラッタ電力の比を算出し、ドプラ電力の割合がなるべく大きくなるように調整することができる。また望ましくは、前記復調手段へ出力される受信信号に対して遅延処理を施す遅延回路を含む、ことを特徴とする。   Preferably, it includes a delay circuit that performs a delay process on the modulated wave transmission signal applied to the demodulation means. According to this configuration, for example, the ratio of Doppler signal / clutter power can be calculated for each delay time by adjusting the delay time, and the ratio of Doppler power can be adjusted to be as large as possible. Desirably, a delay circuit for delaying the received signal output to the demodulator is included.

望ましくは、前記ドプラ信号抽出部に与える参照信号に対して遅延処理を施す遅延回路を含む、ことを特徴とする。また望ましくは、前記ドプラ信号抽出部へ出力される側帯波に対して遅延処理を施す遅延回路を含む、ことを特徴とする。   Preferably, it includes a delay circuit that performs a delay process on a reference signal supplied to the Doppler signal extraction unit. Preferably, a delay circuit that performs a delay process on the sideband output to the Doppler signal extraction unit is included.

望ましくは、前記ドプラ情報抽出手段は、前記変調波信号を基本波とした場合における基本波成分および第n(nは2以上の自然数)高調波成分のうちの複数成分の各々からドプラ信号を抽出し、抽出した複数のドプラ信号を加算処理する、ことを特徴とする。望ましくは、前記変調波信号は正弦波、三角波または鋸歯状波である、ことを特徴とする。   Preferably, the Doppler information extraction unit extracts a Doppler signal from each of a plurality of components of a fundamental wave component and an nth (n is a natural number of 2 or more) harmonic component when the modulated wave signal is a fundamental wave. And adding a plurality of extracted Doppler signals. Preferably, the modulated wave signal is a sine wave, a triangular wave, or a sawtooth wave.

上記目的を達成するために、本発明の別の好適な態様である超音波診断装置は、周波数f0の搬送波信号を周波数fmの変調波信号を用いて周波数変調処理することにより、変調送信信号を生成する送信部と、前記変調送信信号の供給によって生体に対して超音波を送波し、生体からの反射波を受波して受信信号を出力する送受波器と、前記受信信号に対して、前記変調送信信号を参照信号として直交検波処理を行い、これにより復調信号を得る検波回路と、前記復調信号から、前記変調波信号を基本波とした場合における基本波成分および第n(nは2以上の自然数)高調波成分のうちの少なくとも一つの成分とその付近に存在するドプラ信号とから成る側帯波を抽出するバンドパスフィルタと、前記側帯波から前記ドプラ信号を抽出するドプラ信号抽出部と、前記ドプラ信号から、運動体に関するドプラ情報を取得するドプラ情報取得部と、を有することを特徴とする。 To achieve the above object, an ultrasonic diagnostic apparatus which is another preferred embodiment of the present invention, by performing frequency modulation processing by using the modulated wave signal of the frequency f m of the carrier signal of a frequency f 0, modulation transmission A transmitter that generates a signal, a transmitter / receiver that transmits ultrasonic waves to the living body by supplying the modulated transmission signal, receives a reflected wave from the living body, and outputs a received signal; and On the other hand, a quadrature detection process is performed using the modulated transmission signal as a reference signal, thereby obtaining a demodulated signal, and the fundamental wave component and the nth ( n is a natural number of 2 or more) a bandpass filter for extracting a sideband consisting of at least one of the harmonic components and a Doppler signal present in the vicinity thereof, and a Doppler signal for extracting the Doppler signal from the sideband. It has a plastic signal extraction part and a Doppler information acquisition part which acquires Doppler information about a moving body from the Doppler signal.

望ましくは、前記変調波信号として、周波数fm1の変調波信号をさらに周波数fm2の変調波信号を用いて周波数変調処理することにより得られる多重変調波信号を用いる、ことを特徴とする。周波数変調を多段で行うことにより、さらに低速の領域においてクラッタの影響を低減してドプラ信号を抽出することができる。 Preferably, as the modulation wave signal, a multiplex modulation wave signal obtained by further frequency-modulating the modulation wave signal having the frequency f m1 using the modulation wave signal having the frequency f m2 is used. By performing frequency modulation in multiple stages, the influence of clutter can be reduced and a Doppler signal can be extracted in a lower speed region.

また、上記目的を達成するために、本発明の別の好適な態様である超音波診断装置は、周波数f0の搬送波信号を周波数fmの変調波信号を用いて周波数変調処理することにより、変調送信信号を生成する送信部と、前記変調送信信号の供給によって生体に対して超音波を送波し、生体からの反射波を受波して受信信号を出力する送受波器と、前記受信信号に対して、前記変調送信信号を参照信号として直交検波処理を行い、これにより同相成分復調信号および直交成分復調信号を得る検波回路と、前記同相成分復調信号および直交成分復調信号の各々から、前記変調波信号を基本波とした場合における基本波成分および第n(nは2以上の自然数)高調波成分のうちの少なくとも一つの成分とその付近に存在するドプラ信号とから成る側帯波を抽出するバンドパスフィルタと、前記同相成分復調信号の側帯波および前記直交成分復調信号の側帯波から、前記ドプラ信号を抽出するドプラ信号抽出部と、前記ドプラ信号から、運動体に関するドプラ情報を取得するドプラ情報取得部と、を有し、前記ドプラ信号抽出部は、前記変調波信号の周波数をN倍(Nは自然数)した信号のうちの少なくとも一つを参照信号として、前記同相成分復調信号から抽出された側帯波を直交検波処理する同相成分検波回路および前記直交成分復調信号から抽出された側帯波を直交検波処理する直交成分検波回路を含む、ことを特徴とする。 In order to achieve the above object, an ultrasonic diagnostic apparatus which is another preferred embodiment of the present invention, by performing frequency modulation processing by using the modulated wave signal of the frequency f m of the carrier signal of a frequency f 0, A transmitter for generating a modulated transmission signal; a transducer for transmitting an ultrasonic wave to a living body by supplying the modulated transmission signal; receiving a reflected wave from the living body; and outputting a received signal; The signal is subjected to quadrature detection processing using the modulated transmission signal as a reference signal, thereby obtaining an in-phase component demodulated signal and a quadrature component demodulated signal, and each of the in-phase component demodulated signal and the quadrature component demodulated signal, A sideband consisting of at least one of a fundamental wave component and an nth (n is a natural number of 2 or more) harmonic component and a Doppler signal present in the vicinity thereof when the modulated wave signal is a fundamental wave. A Doppler signal extraction unit for extracting the Doppler signal from the sideband wave of the in-phase component demodulated signal and the sideband wave of the quadrature component demodulated signal, and the Doppler signal from the Doppler signal. A Doppler information acquisition unit for acquiring, wherein the Doppler signal extraction unit is configured to demodulate the in-phase component using at least one of signals obtained by multiplying the frequency of the modulated wave signal by N times (N is a natural number) as a reference signal. An in-phase component detection circuit for performing quadrature detection processing on the sideband wave extracted from the signal and a quadrature component detection circuit for performing quadrature detection processing on the sideband wave extracted from the quadrature component demodulation signal are included.

望ましくは、前記ドプラ信号抽出部は、前記同相成分検波回路から出力される二つの信号に対する二乗和の平方根として得られるドプラ信号同相成分、および、前記直交成分検波回路から出力される二つの信号に対する二乗和の平方根として得られるドプラ信号直交成分、の二つのドプラ信号成分を抽出する、ことを特徴とする。   Preferably, the Doppler signal extraction unit applies the Doppler signal in-phase component obtained as the square root of the sum of squares of the two signals output from the in-phase component detection circuit and the two signals output from the quadrature component detection circuit. Two Doppler signal components of the Doppler signal orthogonal component obtained as the square root of the sum of squares are extracted.

本発明により、連続波ドプラ法の適用に当たって、クラッタをできる限り抑圧しつつドプラ情報を感度良く観測できる。   According to the present invention, when applying the continuous wave Doppler method, Doppler information can be observed with high sensitivity while suppressing clutter as much as possible.

以下、本発明の好適な実施形態を図面に基づいて説明する。   DESCRIPTION OF EXEMPLARY EMBODIMENTS Hereinafter, preferred embodiments of the invention will be described with reference to the drawings.

図1には、本発明に係る超音波診断装置の好適な実施形態が示されており、図1はその全体構成を示すブロック図である。送信用振動子10は生体内へ送信波を連続的に送波し、また、受信用振動子12は生体内からの反射波を連続的に受波する。このように、送信および受信がそれぞれ異なる振動子で行われて、いわゆる連続波ドプラ法による送受信が実行される。   FIG. 1 shows a preferred embodiment of an ultrasonic diagnostic apparatus according to the present invention, and FIG. 1 is a block diagram showing the overall configuration thereof. The transmitting vibrator 10 continuously transmits a transmission wave into the living body, and the receiving vibrator 12 continuously receives a reflected wave from the living body. In this way, transmission and reception are performed by different vibrators, and transmission / reception is performed by a so-called continuous wave Doppler method.

送信制御部14は、送信用振動子10を制御して超音波の送信制御を行う。送信制御部14には、FM変調処理が施されたFM連続波(FMCW波)が入力され、このFM連続波に対応する送信波が送信用振動子10から送波される。FM変調器20は、FM連続波を送信制御部14に出力する。FM変調器20は、RF波発振器22から供給されるRF波、および、FM変調波発振器24から供給される変調波に基づいてFM連続波を発生する。このFM連続波の波形については後の原理説明で詳述する。   The transmission control unit 14 controls the transmission transducer 10 to perform ultrasonic transmission control. An FM continuous wave (FMCW wave) that has been subjected to FM modulation processing is input to the transmission control unit 14, and a transmission wave corresponding to the FM continuous wave is transmitted from the transmission transducer 10. The FM modulator 20 outputs an FM continuous wave to the transmission control unit 14. The FM modulator 20 generates an FM continuous wave based on the RF wave supplied from the RF wave oscillator 22 and the modulated wave supplied from the FM modulated wave oscillator 24. The waveform of this FM continuous wave will be described in detail later in the explanation of the principle.

受信制御部16は、受信用振動子12から供給される受波信号に対して増幅処理等の受信処理を施し、受信RF信号を形成して受信ミキサ30へ出力する。受信ミキサ30は、受信RF信号に対して直交検波を施して複素信号を生成する回路であり、2つのミキサ32及び34で構成される。各ミキサは受信RF信号を所定の参照信号と混合する回路である。受信ミキサ30の各ミキサに供給される参照信号は、FM変調器20から出力される。つまり、ミキサ32にはFM変調器20が出力するFM連続波が直接供給され、一方、ミキサ34にはFM連続波がπ/2シフト回路26を経由して供給される。π/2シフト回路26はFM連続波の位相をπ/2だけずらす回路である。   The reception control unit 16 performs reception processing such as amplification processing on the reception signal supplied from the reception transducer 12, forms a reception RF signal, and outputs the reception RF signal to the reception mixer 30. The reception mixer 30 is a circuit that performs quadrature detection on the received RF signal to generate a complex signal, and includes two mixers 32 and 34. Each mixer is a circuit that mixes the received RF signal with a predetermined reference signal. The reference signal supplied to each mixer of the reception mixer 30 is output from the FM modulator 20. That is, the FM continuous wave output from the FM modulator 20 is directly supplied to the mixer 32, while the FM continuous wave is supplied to the mixer 34 via the π / 2 shift circuit 26. The π / 2 shift circuit 26 is a circuit that shifts the phase of the FM continuous wave by π / 2.

後の原理説明で詳述するが、各ミキサで実行される受信RF信号と参照信号との混合処理の結果である受信ミキサ出力信号には、FM変調波発振器24から供給される変調波の変調波周波数fmに関する基本波成分および各高調波成分が含まれている。そして、受信ミキサ30の後段に設けられるBPF(バンドパスフィルタ)36,38によって、これら基本波成分および各高調波成分の中から所望の周波数成分が抽出されて2つのミキサ28,29の入力となる。ミキサ28,29では、N倍回路25で基本波fmをN倍(Nは自然数)して得られる変調波fmの基本波成分あるいは高調波成分を参照信号として、BPF36,38からの各入力を直交検波してベースバンド信号に変換する。ベースバンド信号に変換された2つの信号はLPF37,39を介してドプラ演算処理部40に出力される。この結果、通常の連続波CWに比べて本実施形態のFM連続波(FMCW)ではクラッタの影響が低減される。 As will be described in detail later in the explanation of the principle, the received mixer output signal, which is the result of the mixing process of the received RF signal and the reference signal executed by each mixer, is used to modulate the modulated wave supplied from the FM modulated wave oscillator 24. It contains the fundamental wave component and the harmonic component relating wave frequency f m. A desired frequency component is extracted from these fundamental wave components and each harmonic component by BPFs (band pass filters) 36 and 38 provided at the subsequent stage of the receiving mixer 30, and the inputs of the two mixers 28 and 29. Become. In the mixers 28 and 29, the fundamental wave component or harmonic component of the modulated wave f m obtained by multiplying the fundamental wave f m by N times the N frequency circuit 25 (N is a natural number) is used as a reference signal. The input is quadrature detected and converted to a baseband signal. The two signals converted into the baseband signals are output to the Doppler arithmetic processing unit 40 via the LPFs 37 and 39. As a result, the influence of clutter is reduced in the FM continuous wave (FMCW) of the present embodiment compared to the normal continuous wave CW.

ドプラ演算処理部40は、LPF37,39で抽出された信号に基づいてドプラ演算を実行する。つまり、LPF37,39の出力信号から、ドプラ信号を抽出して例えば血流速度を演算する。もちろん、血流以外の計測対象についての速度演算も可能である。表示処理部42は、速度演算の結果に基づいて例えばドプラ波形を形成し、形成したドプラ波形を表示部44に表示する。   The Doppler computation processing unit 40 performs Doppler computation based on the signals extracted by the LPFs 37 and 39. That is, the Doppler signal is extracted from the output signals of the LPFs 37 and 39 and, for example, the blood flow velocity is calculated. Of course, it is possible to calculate the speed of a measurement target other than the blood flow. The display processing unit 42 forms, for example, a Doppler waveform based on the speed calculation result, and displays the formed Doppler waveform on the display unit 44.

以上、概説したように、本実施形態では、連続波(CW)を正弦波でFM変調した超音波(FMCW波)を送受波し、変調波の基本波あるいは高調波の近傍に発生するドプラ信号の周波数を計測することによりクラッタの影響を低減することができる。そこで、次にその原理について詳述する。   As described above, in this embodiment, a Doppler signal generated in the vicinity of a fundamental wave or a harmonic wave of a modulated wave by transmitting and receiving an ultrasonic wave (FMCW wave) obtained by FM-modulating a continuous wave (CW) with a sine wave. By measuring the frequency, the influence of clutter can be reduced. Then, the principle is explained in full detail next.

周波数f0のCWに変調周波数fmのFM変調を施したFMCW送信波は次式のように表現できる。

Figure 0004688466
数1において、Δfは周波数変動幅の0−P値(ゼロピーク値)である。また、ドプラシ
フトを伴わない場合のFMCW受信波は生体による減衰を無視すると次式で表現できる。
Figure 0004688466
FMCW送信波の周波数スペクトラムは、数1を展開することで得られる。数1に示すFMCW送信波は次式のように展開できる。
Figure 0004688466
数3において、J0(β),J2n(β),J2n+1(β)は、第1種ベッセル関数である。各項の振幅は、変調指数βおよびそれに対応するベッセル関数によって決定される。ちなみに、数3に示されるFMCW送信波の周波数スペクトルを図示すると図2のようになる。 An FMCW transmission wave obtained by subjecting CW of frequency f 0 to FM modulation of modulation frequency f m can be expressed as the following equation.
Figure 0004688466
In Equation 1, Δf is a 0-P value (zero peak value) of the frequency fluctuation range. In addition, the FMCW received wave without Doppler shift can be expressed by the following equation if attenuation by the living body is ignored.
Figure 0004688466
The frequency spectrum of the FMCW transmission wave can be obtained by expanding Equation 1. The FMCW transmission wave shown in Equation 1 can be developed as follows.
Figure 0004688466
In Equation 3, J 0 (β), J 2n (β), and J 2n + 1 (β) are Bessel functions of the first kind. The amplitude of each term is determined by the modulation index β and the corresponding Bessel function. Incidentally, the frequency spectrum of the FMCW transmission wave expressed by Equation 3 is illustrated in FIG.

また、ドプラシフトを伴わない場合の受信波vR(t)の周波数スペクトラムは、数2を展開することで得られる。数2に示すFMCW受信波は次式のように展開できる。

Figure 0004688466
数4に示されるように、受信波の周波数スペクトラムは送信波と同じ周波数成分を持っている。しかし受信波の各周波数成分の振幅は、位相差φ0とφmに応じて変化している。 Further, the frequency spectrum of the received wave v R (t) without Doppler shift can be obtained by developing Equation 2. The FMCW received wave shown in Equation 2 can be developed as follows.
Figure 0004688466
As shown in Equation 4, the frequency spectrum of the received wave has the same frequency component as the transmitted wave. However, the amplitude of each frequency component of the received wave changes according to the phase differences φ 0 and φ m .

さらに、ドプラシフトを伴う場合、数2のvR(t)は以下のように書き換えられる。

Figure 0004688466
なお、数5においてfmに対するドプラシフトは、f0のシフト分fdに比較して小さいので無視している。ちなみに、数5に示されるFMCW受信波の周波数スペクトルを図示すると図3のようになる。図3において、数5に示されるFMCW受信波の周波数スペクトルは破線で示されている。実線で示されるスペクトルはドプラシフトを伴わない場合のFMCW受信波の周波数スペクトルを示している。 Furthermore, when accompanied by Doppler shift, v R (t) in Equation 2 is rewritten as follows.
Figure 0004688466
In Formula 5, the Doppler shift with respect to f m is neglected because it is smaller than the shift amount f d of f 0 . Incidentally, the frequency spectrum of the received FMCW wave shown in Equation 5 is shown in FIG. In FIG. 3, the frequency spectrum of the FMCW received wave shown in Equation 5 is indicated by a broken line. The spectrum indicated by the solid line indicates the frequency spectrum of the FMCW received wave without Doppler shift.

上述の数2や数5で表される受信波形は、超音波振動子が受信する信号波形(受信RF信号)である。超音波診断装置は受信RF信号に対して復調処理を実行する。FMCWの受信RF信号を復調する場合、復調系では参照信号としてFMCW波を用いて受信波と乗算を行う。復調系における受信ミキサ出力は、vT(t)とvR(t)を乗算した結果として次式のように算出される。

Figure 0004688466
ここで、ベッセル関数に関する次の公式を利用する。
Figure 0004688466
数7の公式を用いると、数6はさらに次式のように計算される。なお、数8では数6における係数1/2を省略する。
Figure 0004688466
数8で表現される受信ミキサ出力の周波数スペクトルを図示すると図4のようになる。図4において線スペクトルで示されるドプラエコー46がドプラ信号に相当する。つまり、ドプラシフトを伴う場合、受信ミキサ出力の基本波および各高調波成分は、FM変調周波数が抑圧されたDSB−SC(Double Sideband−Suppressed Carrier)のスペクトルの形をしている。つまり、ドプラエコーには、FM変調周波数(fm)およびその高調波成分(2fm,3fm,・・・)のスペクトラムは存在しない。なお、図4においてFM変調周波数およびその高調波成分に存在する固定物エコー48のスペクトルについては次のように表現できる。 The reception waveform represented by the above-described formulas 2 and 5 is a signal waveform (received RF signal) received by the ultrasonic transducer. The ultrasonic diagnostic apparatus performs demodulation processing on the received RF signal. When demodulating an FMCW received RF signal, the demodulation system performs multiplication with the received wave using an FMCW wave as a reference signal. The reception mixer output in the demodulation system is calculated by the following equation as a result of multiplying v T (t) and v R (t).
Figure 0004688466
Here, the following formula for the Bessel function is used.
Figure 0004688466
If the formula of Formula 7 is used, Formula 6 is further calculated as follows: In Equation 8, the coefficient 1/2 in Equation 6 is omitted.
Figure 0004688466
The frequency spectrum of the reception mixer output expressed by Equation 8 is as shown in FIG. A Doppler echo 46 indicated by a line spectrum in FIG. 4 corresponds to a Doppler signal. That is, when Doppler shift is involved, the fundamental wave and each harmonic component of the reception mixer output are in the form of a DSB-SC (Double Sideband-Suppressed Carrier) spectrum in which the FM modulation frequency is suppressed. That is, the spectrum of the FM modulation frequency (f m ) and its harmonic components (2f m , 3f m ,...) Does not exist in the Doppler echo. In FIG. 4, the spectrum of the fixed object echo 48 existing in the FM modulation frequency and its harmonic component can be expressed as follows.

ドプラシフトを伴わない場合の受信ミキサ出力は、数8に示されるvD(t)についてfd=0と置いて、以下のように求めることができる。この場合の目標は、移動体と同じ位置にあり、同じ反射率を持った固定目標を想定している。

Figure 0004688466
数9に示すように、固定目標からのエコー(クラッタ)は変調周波数fmの整数倍のところにスペクトラムを生じる。つまり、図4に示した周波数スペクトルにおいて、FM変調周波数およびその高調波成分に固定物エコー48のスペクトラムが発生する。 The reception mixer output without Doppler shift can be obtained as follows by setting f d = 0 for v D (t) shown in Equation 8. The target in this case is assumed to be a fixed target at the same position as the moving body and having the same reflectance.
Figure 0004688466
As shown in Formula 9, echoes from stationary targets (clutter) produces a spectrum at an integral multiple of the modulation frequency f m. That is, in the frequency spectrum shown in FIG. 4, the spectrum of the fixed object echo 48 is generated at the FM modulation frequency and its harmonic component.

上述した原理を利用して、以下にクラッタ電力とドプラ信号電力について説明する。   The clutter power and the Doppler signal power will be described below using the above principle.

まず受信RF信号について、通常のCWとFMCWに関するクラッタ電力を考察する。図5は、受信RF信号における通常のCWでのクラッタの電力スペクトルを示す図であり、生体内深さに応じた各クラッタ電力50とその総和52が示されている。クラッタは固定目標からのエコーの集合体と見なせるので、クラッタ全体の電力スペクトラムの総和52は、各深さにおける固定目標からのエコーを固定目標の数だけ加えた形となる。クラッタ電力は生体内深さに応じて、深いものほど減衰の結果小さくなる。クラッタの電力スペクトラムは、統計的には固定目標からのエコーと相似形となるから、クラッタの電力スペクトラムの総和は、各固定目標からのクラッタと似たような形となる。   First, regarding the received RF signal, the clutter power related to normal CW and FMCW is considered. FIG. 5 is a diagram showing the power spectrum of clutter in a normal CW in the received RF signal, and shows each clutter power 50 and its sum 52 according to the in-vivo depth. Since the clutter can be regarded as a collection of echoes from the fixed target, the total power spectrum 52 of the entire clutter has a form obtained by adding the echoes from the fixed target at each depth by the number of fixed targets. Depending on the in-vivo depth of the clutter power, the deeper the clutter power, the smaller the attenuation. The clutter power spectrum is statistically similar to echoes from fixed targets, so the sum of the clutter power spectrum is similar to the clutter from each fixed target.

一方、図6は、受信RF信号におけるFMCWでのクラッタの電力スペクトルを示す図であり、各周波数ごとに、生体内深さに応じた各クラッタ電力60とその総和62が示されている。FMCWのクラッタは、総電力、つまり、各周波数ごとのFMCWのクラッタ総和を全ての周波数に亘って加算した総和は通常のCWと同じであるが、スペクトラムの形は、送信波形が側帯波を伴っている分、通常のCWとは異なっている。つまり、FMCWでは、FMの変調度が大きくなるにしたがって、f0に集中していた電力が、fmの基本波(f0+fm)および高調波(f0+2fm,f0+3fm,f0+4fm,・・・)に分散していく。FMCWでもクラッタ電力は、深いものほど減衰するが、通常のCWでもFMCWでもこの影響は同様である。 On the other hand, FIG. 6 is a diagram showing the power spectrum of the clutter in the FMCW in the received RF signal, and shows each clutter power 60 according to the in-vivo depth and its sum 62 for each frequency. The FMCW clutter has the same total power, that is, the sum of the FMCW clutter sum for each frequency over all frequencies, but it is the same as the normal CW. Therefore, it is different from normal CW. That is, in the FMCW, according FM modulation degree is increased, the electric power were concentrated in f 0 is the fundamental wave of f m (f 0 + f m ) and harmonic (f 0 + 2f m, f 0 + 3f m, f 0 + 4f m ,...). Even in FMCW, the clutter power attenuates as the depth increases, but this effect is the same in both normal CW and FMCW.

FMCWにおいて、fmの基本波および高調波成分の電圧(一般的には、側帯波の基本波および高調波成分と表現する場合もある)は第1次ベッセル関数により表現できる。第1次ベッセル関数の概略をグラフに示すと図7のようになる。図7には、横軸にβの値が示され縦軸にベッセル関数の値が示されている。この第1次ベッセル関数の曲線を用いて、つまり、FMCWにおける周波数f0に対応するJ0(β)、fmの基本波に対応するJ1(β)、第2高調波に対応するJ2(β)、第3高調波に対応するJ3(β)の値をクラッタの電圧として、通常のCWとFMCW間のクラッタ電力比を推定する。変調度としてはβ=1を仮定すると、図7から、以下の値を得る。

Figure 0004688466
例えば、fmの基本波近傍に注目すると、FMCWでは、通常のCWに比べ、クラッタ電力が7.0dB低減されている。但し、受信RF信号は復調処理が実行される前の何らの信号処理も施していない信号である。つまり、受信RF信号としては、ドプラ情報を必要とする目標からのエコーも、変調度βに応じてクラッタと同様の割合で電力が低減する。ゆえに、受信RF信号における目標からのエコーのクラッタに対する電力比は、FMCWにおいても通常のCWと同じである。 In FMCW, (generally, it may be expressed as fundamental and harmonic components of sidebands) voltage of the fundamental wave and the harmonic component of f m can be represented by a first-order Bessel function. The outline of the first-order Bessel function is shown in FIG. In FIG. 7, the value of β is shown on the horizontal axis, and the value of the Bessel function is shown on the vertical axis. Using this first-order Bessel function curve, that is, J 0 (β) corresponding to the frequency f 0 in FMCW, J 1 (β) corresponding to the fundamental wave of f m , and J corresponding to the second harmonic. 2 (β), the value of J 3 (β) corresponding to the third harmonic is used as the clutter voltage, and the clutter power ratio between normal CW and FMCW is estimated. Assuming that β = 1 as the degree of modulation, the following values are obtained from FIG.
Figure 0004688466
For example, focusing on the fundamental wave near the f m, the FMCW, compared to conventional CW, clutter power is reduced 7.0 dB. However, the received RF signal is a signal that has not been subjected to any signal processing before the demodulation processing is executed. That is, as the received RF signal, the power of the echo from the target that requires Doppler information is also reduced at the same rate as the clutter according to the modulation degree β. Therefore, the power ratio of the echo from the target to the clutter in the received RF signal is the same as the normal CW in the FMCW.

次に、受信ミキサ出力について、通常のCWとFMCWに関するクラッタ電力を考察する。通常のCWでは、受信ミキサによりIQチャンネルに分離し、ベースバンド帯域で信号処理している。この場合、ミキサの機能は単純な周波数シフトであるから、RF帯におけるドプラ信号とクラッタのレベル関係は、そのままベースバンド帯域でも維持される。FMCW波を送信した場合でも、受信ミキサのレファレンスにCW信号を用いれば、ドプラ信号とクラッタのレベル関係は、そのままベースバンド帯域でも維持される。したがって、ドプラ信号とクラッタのレベル関係は改善されない。ところが、本実施形態で用いるFMCWでは、受信ミキサ(図1の符号30)の参照信号(レファレンス)として、送信波であるFMCW波を用いる。このため、ドプラ信号とクラッタのレベル比改善ができる。この場合、ドプラ信号レベルは、数8で表現されるように、変調度βと距離に依存するパラメータkの影響を受ける。すなわち、距離に依存してエコー電力が増減する。クラッタも同様の効果を受ける。ただし、ドプラ信号の場合と異なるのは、クラッタは深さ方向に渡って分布しているので、各場所におけるクラッタの距離依存性を全て考慮しなければならない点である。この様子を定量的に解析する。   Next, the clutter power for normal CW and FMCW is considered for the receiving mixer output. In normal CW, the signal is separated into IQ channels by a receiving mixer and signal processing is performed in the baseband. In this case, since the function of the mixer is a simple frequency shift, the level relationship between the Doppler signal and the clutter in the RF band is maintained as it is in the baseband band. Even when an FMCW wave is transmitted, if a CW signal is used for the reference of the receiving mixer, the level relationship between the Doppler signal and the clutter is maintained as it is in the baseband band. Therefore, the level relationship between the Doppler signal and the clutter is not improved. However, in the FMCW used in the present embodiment, an FMCW wave that is a transmission wave is used as a reference signal (reference) of the reception mixer (reference numeral 30 in FIG. 1). For this reason, the level ratio of the Doppler signal and the clutter can be improved. In this case, the Doppler signal level is affected by a modulation factor β and a parameter k that depends on the distance, as expressed by Equation (8). That is, the echo power increases or decreases depending on the distance. Clutter has the same effect. However, the difference from the case of the Doppler signal is that the clutter is distributed in the depth direction, so that all the distance dependency of the clutter at each location must be considered. This situation is analyzed quantitatively.

解析に先立ち、数8における各成分の振幅に比例する係数であるベッセル関数の距離依存性を求める。代表例として、以下の条件を仮定する。

Figure 0004688466
ここで、dは振動子から目標までの距離である。数10に示される各パラメータを用いて、J0(kβ), J1(kβ),J2(kβ),J3(kβ)をd=30cmまで計算し、横軸を診断の深さ(生体内の深さ)d、縦軸をベッセル関数値としてプロットすると図8のようになる。図8に示される物理的意味は以下のとおりである。 Prior to the analysis, the distance dependency of the Bessel function, which is a coefficient proportional to the amplitude of each component in Equation 8, is obtained. As a representative example, the following conditions are assumed.
Figure 0004688466
Here, d is the distance from the transducer to the target. Using each parameter shown in Equation 10, J 0 (kβ), J 1 (kβ), J 2 (kβ), and J 3 (kβ) are calculated up to d = 30 cm, and the horizontal axis is the depth of diagnosis ( Plotting the in vivo depth (d) and the vertical axis as a Bessel function value is as shown in FIG. The physical meanings shown in FIG. 8 are as follows.

(1)d=0では、送信FMCW波を参照波としているために、FMCWの高調波成分はすべて直流付近(RF帯ではf0相当)に集まってくる。つまり、fmの基本波および高調波成分は出現しない。(2)dの値がだんだん大きくなるにつれ、kβも大きくなり、fmの基本波および高調波成分が出現してくる。(3)したがって、深部からのエコーのミキサ出力ほどfmの基本波および高調波成分の電力が増大する。この例では、0〜15cmの範囲で増大し、逆に15〜30cmにおけるエコーでは減少している。(4)kβは30cmの周期で、0と±2の間の値をとる。(5)したがって、ベッセル関数のkβ値は、0と±2の間の値を周期的に往復する。(6)図8のグラフを参照すると、0と2の間で支配的なベッセル関数の次数は0次、1次、2次および3次のみである。しかし、この効果は、βが大きくなるにつれて支配的な次数も増加する傾向にある。(7)したがって、FMCWにおけるクラッタの高調波成分は、β=1の場合、3次まで考慮すれば充分であるが、βをより大きい値に設定すると3次以上でも改善効果がある。 (1) Since the transmission FMCW wave is used as a reference wave at d = 0, all the harmonic components of the FMCW gather near the direct current (corresponding to f 0 in the RF band). In other words, the fundamental wave and the harmonic component of f m does not appear. (2) as the value of d is gradually increased, Keibeta also increases, fundamental and harmonic components of f m comes to appear. (3) Therefore, the power of the fundamental and harmonic components of f m as echo mixer output of the deep is increased. In this example, it increases in the range of 0 to 15 cm, and conversely decreases in the echo at 15 to 30 cm. (4) kβ takes a value between 0 and ± 2 with a period of 30 cm. (5) Therefore, the kβ value of the Bessel function periodically reciprocates between 0 and ± 2. (6) Referring to the graph of FIG. 8, the order of the Bessel function dominant between 0 and 2 is only the 0th order, the 1st order, the 2nd order and the 3rd order. However, this effect tends to increase the dominant order as β increases. (7) Therefore, the harmonic component of the clutter in FMCW is sufficient if the third order is considered when β = 1, but if β is set to a larger value, there is an improvement effect even at the third or higher order.

図8は、FMCW方式の受信ミキサ出力におけるfmの基本波および高調波信号成分の電圧が深さによってどのように変化するかに相当する。すなわち、生体内に分布する組織からの各クラッタは、その位置に応じて、基本波および高調波信号成分ごとに振幅が変動する。その様子を各周波数成分ごとに表現しなおしたのが図9である。 8 corresponds to how varies with the fundamental wave and the voltage is the depth of the harmonic signal components of f m in the receiving mixer output of FMCW system. That is, the amplitude of each clutter from the tissue distributed in the living body varies for each fundamental wave and harmonic signal component according to the position. FIG. 9 shows the state again for each frequency component.

図9は、基本波および高調波クラッタ電圧の深さ依存性を説明するための図であり、図9には、各周波数ごとに、生体内深さに応じた各クラッタ電力90とその総和92が示されている。なお、図9は、生体内の減衰が存在しないと仮定した場合、各組織からのクラッタ電力を示したものである。ミキサ出力は、各周波数において、各深さからのクラッタの総和であるから、総和の電力スペクトラムは通常のCWの場合とは異なった形となる。   FIG. 9 is a diagram for explaining the depth dependency of the fundamental wave and the harmonic clutter voltage. FIG. 9 shows each clutter power 90 corresponding to the in-vivo depth and the total sum 92 for each frequency. It is shown. FIG. 9 shows the clutter power from each tissue when it is assumed that there is no attenuation in the living body. Since the mixer output is the sum of the clutter from each depth at each frequency, the power spectrum of the sum is different from that in the case of normal CW.

図8および図9から、β=1の場合、クラッタが支配的な周波数帯は、DC、基本波(fm)、その2倍(2fm)および3倍(3fm)の高調波帯域であることがわかる。しかし、βが大きくなるにつれて次数も高次となる。また、図9は、クラッタの電力が、見かけ上深さに対して、周期的に変動していることを示している。これら周期的な変動は図8に示したベッセル関数の挙動に対応している。実際のエコーは深さ方向のすべてのエコーからの総和である。各組織からのエコーは、その位置に応じて、基本波および高調波信号成分の振幅変動を受けると同時に、生体における減衰を受ける。生体内の減衰特性を定性的に表現したのが図9に示す生体内の減衰特性である。生体内の減衰特性は、組織の深さをd(cm)とすると次式で与えられる。

Figure 0004688466
数11において、0.6に掛けられた2は、生体内での超音波の往復を意味している。 From FIG. 8 and FIG. 9, when β = 1, the frequency band in which the clutter is dominant is DC, the fundamental wave (f m ), the harmonic band of twice (2f m ) and three times (3f m ) thereof. I know that there is. However, the order becomes higher as β increases. FIG. 9 shows that the power of the clutter fluctuates periodically with respect to the apparent depth. These periodic fluctuations correspond to the behavior of the Bessel function shown in FIG. The actual echo is the sum of all echoes in the depth direction. Echoes from each tissue are subject to attenuation in the living body at the same time that they undergo amplitude fluctuations of the fundamental wave and harmonic signal components according to their positions. The attenuation characteristic in the living body shown in FIG. 9 expresses the attenuation characteristic in the living body qualitatively. The attenuation characteristic in the living body is given by the following equation when the depth of the tissue is d (cm).
Figure 0004688466
In Equation 11, 2 multiplied by 0.6 means the reciprocation of the ultrasonic wave in the living body.

図9には数11の減衰特性も表現されている。図9における減衰特性(Attenuation)は、数11において、f0=3MHzとした場合の計算結果をプロットしたものである。クラッタの総和は、図9において、減衰特性と各周波数成分ごとの深さ方向の成分を、各深さごとに求めてそれを積算すればよい。つまり、実際のミキサ出力は、図9において、各深さにおけるベッセル関数値と減衰量を掛け合せた値(電圧)に比例する。クラッタ電力を算出するには、各周波数について、深さ方向からのすべてのエコーを加算すればよい。つまり、各ベッセル関数と減衰曲線と掛け合わせた曲線について、その曲線の面積を求めることによって、クラッタ電力を求めることができる。ただし、組織内には、反射係数の等しい固定目標が均一に分布しているとする。こうして算出された電力の結果を表2に示す。

Figure 0004688466
なお、表2では、各成分の相対電圧積分値をそれぞれ求め、相対電圧の合計が“1”となるように規格化した。 In FIG. 9, the attenuation characteristic of Formula 11 is also expressed. The attenuation characteristic (Attenuation) in FIG. 9 is obtained by plotting the calculation result when f 0 = 3 MHz in Equation 11. The sum of clutter may be obtained by obtaining the attenuation characteristics and the components in the depth direction for each frequency component for each depth in FIG. That is, the actual mixer output is proportional to the value (voltage) obtained by multiplying the Bessel function value at each depth by the attenuation amount in FIG. In order to calculate the clutter power, all echoes from the depth direction may be added for each frequency. That is, the clutter power can be obtained by obtaining the area of the curve obtained by multiplying each Bessel function and the attenuation curve. However, it is assumed that fixed targets having the same reflection coefficient are uniformly distributed in the tissue. Table 2 shows the power results thus calculated.
Figure 0004688466
In Table 2, the relative voltage integral value of each component was obtained and normalized so that the sum of the relative voltages was “1”.

ここで、受信RF信号に関する表1と受信ミキサ出力に関する表2のそれぞれに示される結果をまとめると表3のようになる。

Figure 0004688466
表3では、いずれの方式でも、送信CW電力が一定であれば、クラッタ電力の総和は不変であるということを前提としている。この表のなかで、相対電力とは、通常のCWクラッタ電力の総和を0dBとした場合の相対値である。例えば、FMCWのfo 成分が−1.9dBということは、受信RF信号で、無変調の場合にくらべ、FM変調により、fo の成分の電力が1.9dBだけ高調波に移動し、その分だけfo の電力が減少しているということを意味している。 Here, the results shown in Table 1 for the received RF signal and Table 2 for the received mixer output are summarized in Table 3.
Figure 0004688466
In Table 3, it is assumed that the sum of the clutter power remains unchanged if the transmission CW power is constant in any method. In this table, the relative power is a relative value when the sum of normal CW clutter power is 0 dB. For example, if the F o component of FMCW is −1.9 dB, the power of the component of f o is shifted to a harmonic by 1.9 dB by FM modulation compared to the case of no reception modulation in the received RF signal. This means that the power of f o is reduced by the amount.

表3から、次のことがわかる。(1)foの成分、あるいはDC成分は、受信RF信号でも、受信ミキサ出力でも、電力に大差は無い。これは、受信ミキサ出力では、受信RF信号とほとんど時間差の無い参照信号を乗算することにより発生するクラッタDC成分が支配的となっているからである。深部からのクラッタ電力は見かけ上、変調度が大きくなるように見え、DC成分の寄与が少なくなるばかりか、減衰が大きいので、DC成分に比べれば、その電力は非常に小さい。(2)受信ミキサ出力のクラッタfm基本波成分は、通常のCWにくらべ、25.7dB(約1/370)低い。これは、受信ミキサ出力では、参照信号との乗算により、fmの基本波成分および高調波成分が極端に抑圧されるからである。この抑圧効果は、クラッタが体表面に近いほど効果が大きい。深部になるにしたがい、この抑圧効果は少なくなる。すなわち、参照信号とクラッタの時間差のために、fmの基本波成分および高調波成分が増大する。しかし、これら深部領域からのクラッタは、数11の減衰項が支配的になるので、ほとんど無視できる。(3)fmの基本波成分および高調波のクラッタ成分がfoのクラッタ成分に比較して数十dBも抑圧されるので、fmの基本波成分あるいは高調波近傍のドプラ信号を抽出すれば、クラッタの除去が非常に容易となる。 Table 3 shows the following. (1) There is no great difference in power between the component of f o or the DC component, whether it is a received RF signal or a received mixer output. This is because the clutter DC component generated by multiplying the received RF signal by the reference signal having almost no time difference is dominant in the receiving mixer output. The clutter power from the deep part appears to increase the modulation degree, and the contribution of the DC component is not only reduced, but also the attenuation is large, so that the power is very small compared to the DC component. (2) clutter f m fundamental component of the receiving mixer output is compared with the normal CW, 25.7DB (about 1/370) low. This is because the receiving mixer output by multiplying the reference signal, because the fundamental wave component and harmonic components of f m is extremely suppressed. This suppression effect is more effective as the clutter is closer to the body surface. As the depth increases, this suppression effect decreases. That is, because of the time difference between the reference signal and the clutter, the fundamental wave component and harmonic components of f m increases. However, the clutter from these deep regions is almost negligible because the attenuation term of Formula 11 becomes dominant. (3) Since the clutter components of the fundamental wave component and harmonic of f m is tens dB compared to the clutter components of f o is suppressed, by extracting the Doppler signal of the fundamental wave component or harmonic vicinity of f m In this case, the clutter can be removed very easily.

次に、受信ミキサ出力について、通常のCWとFMCWに関するドプラ信号を考察する。例として、信号レベルとしては、深さ15cmからのエコーを想定する。通常のCWクラッタを基準とし、信号レベルとして、基本波fmの近傍のドプラ信号抽出を仮定する。数8によれば、基本波信号は次式で表現される。

Figure 0004688466
深さによるkβの影響は、図8における15cmでのJ1(kβ)の値として、0.58と求めることができる。0.58の意味は、FMCW方式において、受信ミキサのレファレンスとして送信FMCW波を用いると、通常のCWにくらべ、往復で信号が0.58(電力比で、0.34倍、−4.7dB)となるということである。信号レベルとして、第2高調波2fmの近傍のドプラ信号抽出を行う場合は、図8における15cmでのJ2(kβ)の値として、0.37(電力比で、0.137倍、−8.6dB)を得ることができる。 Next, consider the Doppler signals for normal CW and FMCW for the receiving mixer output. As an example, an echo from a depth of 15 cm is assumed as the signal level. With respect to the normal CW clutter, as a signal level, it is assumed Doppler signal extracted in the vicinity of the fundamental wave f m. According to Equation 8, the fundamental wave signal is expressed by the following equation.
Figure 0004688466
The influence of kβ due to the depth can be obtained as 0.58 as the value of J 1 (kβ) at 15 cm in FIG. The meaning of 0.58 means that when a transmission FMCW wave is used as a reference of the reception mixer in the FMCW system, the signal is 0.58 (round and round by power ratio, 0.34 times, -4.7 dB) compared to normal CW. ). As a signal level, when performing Doppler signal extraction in the vicinity of the second harmonic 2f m is the value of J 2 (kβ) at 15cm in FIG 8, 0.37 (power ratio, 0.137 times, - 8.6 dB) can be obtained.

一方、通常のCWにおけるドプラ信号は、FMCWのようなFM変調に伴う距離依存性が無い。したがって、生体による減衰を両者共通とすると、FMCWのドプラ信号は、通常のCWにくらべ、基本波の場合4.7dB、第2高調波の場合、8.6dB低いことになる。これらの結果を、直流および第3高調波3fmの近傍のドプラ信号電力とともにまとめると表4のようになる。

Figure 0004688466
表4において、相対電圧は図8における15cmでの各ベッセル関数の値であり、相対電圧の2乗として相対電力が示されている。 On the other hand, a Doppler signal in a normal CW has no distance dependency associated with FM modulation like FMCW. Therefore, assuming that attenuation by the living body is common, the FMCW Doppler signal is 4.7 dB lower for the fundamental wave and 8.6 dB lower for the second harmonic wave than the normal CW. The results are summarized together with the Doppler signal power in the vicinity of the direct current and the third harmonic 3f m is as in Table 4.
Figure 0004688466
In Table 4, the relative voltage is the value of each Bessel function at 15 cm in FIG. 8, and the relative power is shown as the square of the relative voltage.

受信ミキサ出力のクラッタ電力についての結果を示す表2と、受信ミキサ出力のドプラ信号電力についての結果を示す表4とを、さらにまとめたものが表5である。

Figure 0004688466
表5では、通常のCWドプラ信号電力のクラッタ電力に対する比をαdBとした。また、クラッタは生体内のすべてのエコーからの総電力和、ドプラ信号は深さ15cmにおける移動目標からの信号を想定した。さらに、表5に示されるドプラ信号電力とクラッタ電力との電力比関係を図10に示す。図10は、通常のCWとFMCWのそれぞれについて、横軸を周波数、縦軸を電力値としてドプラ信号電力100とクラッタ電力102のそれぞれの波形が示されている。 Table 5 shows a further summary of Table 2 showing the result of the clutter power of the reception mixer output and Table 4 showing the result of the Doppler signal power of the reception mixer output.
Figure 0004688466
In Table 5, the ratio of normal CW Doppler signal power to clutter power is α dB. The clutter was assumed to be the sum of the total power from all echoes in the living body, and the Doppler signal was assumed to be a signal from a moving target at a depth of 15 cm. Furthermore, the power ratio relationship between the Doppler signal power and the clutter power shown in Table 5 is shown in FIG. FIG. 10 shows the respective waveforms of the Doppler signal power 100 and the clutter power 102 for normal CW and FMCW, with the horizontal axis representing the frequency and the vertical axis representing the power value.

表5および図10から、FMCWは通常のCWに比較して、深さ15cm、β=1の場合、クラッタを以下のとおり抑圧できることがわかる。(1)基本波近傍のドプラを抽出する場合:21.0dB(約1/126)、(2)第2高調波近傍のドプラを抽出する場合:30.8dB(約1/1200)、(3)第3高調波近傍のドプラを抽出する場合:45.6dB(約1/36000)。このように、FMCW波を利用することにより、つまり、図1に示す超音波診断装置により、通常のCWに比べてクラッタの影響を低減することができる。   From Table 5 and FIG. 10, it can be seen that the FMCW can suppress the clutter as follows when the depth is 15 cm and β = 1 as compared with the normal CW. (1) When extracting the Doppler near the fundamental wave: 21.0 dB (about 1/126), (2) When extracting the Doppler near the second harmonic: 30.8 dB (about 1/1200), (3 ) When extracting Doppler in the vicinity of the third harmonic: 45.6 dB (about 1/36000). As described above, by using the FMCW wave, that is, by the ultrasonic diagnostic apparatus shown in FIG. 1, the influence of clutter can be reduced as compared with normal CW.

図11は、図1に示す超音波診断装置を改良して、クラッタの影響をさらに低減した超音波診断装置を示している。図1と異なる点は、図11に示す超音波診断装置では、可変遅延回路110が、参照波として用いている送信FMCW変調波の遅延時間調整を行ってドプラ信号電力/クラッタ電力の最適化を行うことである。   FIG. 11 shows an ultrasonic diagnostic apparatus in which the influence of clutter is further reduced by improving the ultrasonic diagnostic apparatus shown in FIG. The difference from FIG. 1 is that in the ultrasound diagnostic apparatus shown in FIG. 11, the variable delay circuit 110 adjusts the delay time of the transmission FMCW modulated wave used as the reference wave to optimize the Doppler signal power / clutter power. Is to do.

生体内の深さに応じて、受信ミキサ出力における積算されたクラッタ電力とドプラ信号電力と相対的な大きさの関係は異なってくる。表5においては生体内の深さを15cmとして計算した結果を示したが、深さに応じて、ドプラ信号対クラッタ電力の比が変化する。診断深さはターゲットの位置に応じて決定されるため容易に変更することはできない。しかし、受信ミキサに供給される参照波に遅延時間を設けて調整することにより、擬似的にターゲットの位置、つまり診断深さを変更させることが可能になる。そこで、図11に示す超音波診断装置では、受信ミキサの前段に可変遅延回路110を設けて、参照波の遅延時間調整を行えるようにした。この結果、遅延時間を調整して、各遅延時間ごとにドプラ信号電力/クラッタ電力の比を算出し、ドプラ信号電力の割合がなるべく大きくなるように調整することができる。一方、変調周波数および変調波の振幅を変化させる(βを変化させる)ことにより、生体内の深さに応じて受信ミキサ出力における積算クラッタ電力とドプラ信号電力の相対比を改善することができる。   Depending on the depth in the living body, the relationship between the relative magnitudes of the accumulated clutter power and Doppler signal power at the output of the receiving mixer differs. Table 5 shows the results of calculation assuming that the depth in the living body is 15 cm, but the ratio of Doppler signal to clutter power varies depending on the depth. Since the diagnostic depth is determined according to the position of the target, it cannot be easily changed. However, by adjusting the reference wave supplied to the receiving mixer by providing a delay time, it becomes possible to change the position of the target, that is, the diagnostic depth in a pseudo manner. Therefore, in the ultrasonic diagnostic apparatus shown in FIG. 11, the variable delay circuit 110 is provided in the preceding stage of the reception mixer so that the delay time of the reference wave can be adjusted. As a result, it is possible to adjust the delay time, calculate the ratio of Doppler signal power / clutter power for each delay time, and adjust the Doppler signal power ratio to be as large as possible. On the other hand, by changing the modulation frequency and the amplitude of the modulation wave (by changing β), it is possible to improve the relative ratio of the accumulated clutter power and the Doppler signal power at the output of the receiving mixer according to the depth in the living body.

なお、遅延時間の調整は、受信ミキサ30へ供給される受信RF信号に対して行われてもよい。つまり、受信制御部16の後段に可変遅延回路を設けて、受信ミキサ30に供給される受信RF信号の遅延時間を調整することにより、参照波の遅延時間調整と同様な効果を得ることができる。   The delay time may be adjusted with respect to the reception RF signal supplied to the reception mixer 30. That is, by providing a variable delay circuit in the subsequent stage of the reception control unit 16 and adjusting the delay time of the reception RF signal supplied to the reception mixer 30, the same effect as the delay time adjustment of the reference wave can be obtained. .

図12は、図11に示す超音波診断装置をさらに改良した超音波診断装置を示している。図11と異なる点は、図12に示す超音波診断装置では、可変遅延回路112が、N倍回路25の後段に設けられている点である。   FIG. 12 shows an ultrasonic diagnostic apparatus in which the ultrasonic diagnostic apparatus shown in FIG. 11 is further improved. The difference from FIG. 11 is that the variable delay circuit 112 is provided in the subsequent stage of the N-times circuit 25 in the ultrasonic diagnostic apparatus shown in FIG.

ミキサ28,29は、N倍回路25で基本波fmをN倍(Nは自然数)して得られる変調波fmの基本波成分あるいは高調波成分を参照信号として、BPF36,38から出力される各信号を直交検波してベースバンド信号に変換する。受信ミキサ30から出力される受信ミキサ出力信号には、FM変調波発振器24から供給される変調波の変調波周波数fmに関する基本波成分および各高調波成分が含まれている。受信ミキサ出力信号の一例を数式で表現したのが前述の数8である。BPF36,38は、数8に示されるfmに関する基本波成分および各高調波成分の中から所望の周波数成分を抽出して、2つのミキサ28,29に出力する。例えば、BPF36によって、数8に示されるfmに関する基本波成分が抽出されたとすると、BPF36からミキサ28に供給される信号は次式で表現される。

Figure 0004688466
したがって、ミキサ28の参照信号としては、cos(2πfmt+φm/2)を用いるのがSNR(信号対ノイズ比)の良い乗算結果となる。このため、可変遅延回路112は、N倍回路25から出力される参照信号の位相を調整してミキサ28にcos(2πfmt+φm/2)の参照信号を供給する。この結果、SNRの良い乗算結果を得ることができる。 The mixer 28 and 29, the fundamental wave f m as N times (N is a natural number) the reference signal fundamental wave component or a harmonic component of the modulation wave f m obtained by N times circuit 25, is output from the BPF36,38 Each signal is orthogonally detected and converted to a baseband signal. The receiving mixer output signal output from the receiving mixer 30 includes a fundamental wave component and the harmonic component about the modulation wave frequency f m of the modulation wave supplied from the FM modulation wave oscillator 24. An example of the reception mixer output signal is expressed by the mathematical formula 8 described above. The BPFs 36 and 38 extract a desired frequency component from the fundamental wave component and the harmonic components related to f m shown in Equation 8 and output them to the two mixers 28 and 29. For example, the BPF 36, when the fundamental wave component relating f m as shown in Equation 8 is extracted, the signal supplied from the BPF 36 to the mixer 28 is expressed by the following equation.
Figure 0004688466
Therefore, using cos (2πf m t + φ m / 2) as a reference signal for the mixer 28 is a multiplication result with a good SNR (signal to noise ratio). Therefore, the variable delay circuit 112 adjusts the phase of the reference signal output from the N-times circuit 25 and supplies the reference signal of cos (2πf m t + φ m / 2) to the mixer 28. As a result, a multiplication result having a good SNR can be obtained.

なお、遅延時間の調整は、BPF36,38の出力に対して行われてもよい。つまり、BPF36,38の各々の後段に可変遅延回路を設けて、ミキサ28,29に供給される信号の遅延時間を調整することにより、参照信号を遅延時間調整する場合と同様な効果を得ることができる。   The delay time may be adjusted for the outputs of the BPFs 36 and 38. In other words, by providing a variable delay circuit in the subsequent stage of each of the BPFs 36 and 38 and adjusting the delay time of the signals supplied to the mixers 28 and 29, the same effect as that of adjusting the delay time of the reference signal can be obtained. Can do.

ここで、ドプラ信号を検出する際の遅延時間の調整について検討する。   Here, the adjustment of the delay time when detecting the Doppler signal will be considered.

図12における受信ミキサ30(以降、必要に応じて「第1ミキサ」と称する)により低周波帯域に復調されたベースバンド信号は、ミキサ28,29(以降、必要に応じて「第2ミキサ」と称する)により通常のドプラ信号と同様、直流付近のドプラ信号に変換される。図13に示すように、FMCWでは、第1ミキサ出力におけるドプラ信号(ドプラエコー46)は、FM変調周波数(fm)およびその高調波成分(2fm,3fm,・・・)の固体物エコー(クラッタ信号)48の近傍にDSB−SC(Double Sideband−Suppressed Carrier)の形で出現するから、この信号を直流付近の帯域に変換するには、AMの同期検波に相当する復調プロセスが必要となる。 The baseband signal demodulated to the low frequency band by the receiving mixer 30 (hereinafter referred to as “first mixer” as necessary) in FIG. 12 is mixed with the mixers 28 and 29 (hereinafter referred to as “second mixer” as necessary). Is converted into a Doppler signal in the vicinity of DC. As shown in FIG. 13, in the FMCW, the Doppler signal (Doppler echo 46) at the first mixer output is a solid object echo of the FM modulation frequency (f m ) and its harmonic components (2f m , 3f m ,...). Since it appears in the form of DSB-SC (Double Sideband-Suppressed Carrier) in the vicinity of (clutter signal) 48, a demodulation process corresponding to AM synchronous detection is required to convert this signal to a band near DC. Become.

例えば、高調波成分2fmの側帯波として出現しているドプラ信号を抽出する場合は、第1ミキサ出力を、第2ミキサにおいて周波数2fmの第2参照波により同期検波する。図12において、N倍回路25の次数Nを2に設定すればこのプロセスを実現できる。第1ミキサ後段のBPF36,38がローパスフィルタとして機能することにより、第2高調波成分を選択したドプラ信号に対する同期検波過程は、以下のように表現できる。 For example, to extract Doppler signals have emerged as a side band harmonic component 2f m is the first mixer output to synchronous detection by the second reference wave frequency 2f m in the second mixer. In FIG. 12, this process can be realized by setting the order N of the N multiplication circuit 25 to 2. Since the BPFs 36 and 38 subsequent to the first mixer function as a low-pass filter, the synchronous detection process for the Doppler signal in which the second harmonic component is selected can be expressed as follows.

まず、第1ミキサ出力(この場合はIチャンネルのミキサ32の出力を仮定)は前述の数8で表現できる。BPF36により第2高調波成分を選択すると次式の信号が抽出される。

Figure 0004688466
また、周波数2fmの第2参照波を次式のように表現する。
Figure 0004688466
数14および数15から、第1ミキサ後段のBPF36により第2高調波成分を選択したドプラ信号に対する同期検波過程は、次式によって表現される。
Figure 0004688466
数16は、第2ミキサ出力をLPF37によりドプラ周波数帯域以上の高域成分を除去した結果に相当する。 First, the first mixer output (in this case, the output of the I-channel mixer 32 is assumed) can be expressed by the above equation (8). When the second harmonic component is selected by the BPF 36, the following signal is extracted.
Figure 0004688466
Also, to represent the second reference wave frequency 2f m as follows.
Figure 0004688466
From the equations 14 and 15, the synchronous detection process for the Doppler signal for which the second harmonic component is selected by the BPF 36 subsequent to the first mixer is expressed by the following equation.
Figure 0004688466
Equation 16 corresponds to the result of removing the high-frequency component higher than the Doppler frequency band from the second mixer output by the LPF 37.

数16によれば、ドプラ信号vD2(t)の振幅は、J2(kβ),φmおよびφrに依存する。J2(kβ)とφmはターゲットの深さにより決定されるが、φrは信号処理系(例えば図12における可変遅延回路112)で設定されるパラメータで、vD2(t)はφrにより周期的に変化する。 According to Equation 16, the amplitude of the Doppler signal v D2 (t) depends on J 2 (kβ), φ m and φ r . J 2 (kβ) and φ m are determined by the depth of the target, φ r is a parameter set in the signal processing system (for example, variable delay circuit 112 in FIG. 12), and v D2 (t) is φ r It changes periodically.

ドプラターゲットの深さを固定し、第2参照波vr2(t)の位相φrを変化させた場合の第2ミキサ出力vD2(t)の電力を測定した結果を図14に示す。 FIG. 14 shows the result of measuring the power of the second mixer output v D2 (t) when the depth of the Doppler target is fixed and the phase φ r of the second reference wave v r2 (t) is changed.

図14において、(a)は基本波(FM変調周波数fm)成分、(b)は二次高調波成分、(c)は三次高調波成分に関する測定結果であり、各々、横軸は遅延時間(位相φr)を表し縦軸は第2ミキサ出力電力を示している。また、図14において、白丸はIチャンネルの第2ミキサ出力測定値、黒丸はQチャンネルの第2ミキサ出力測定値、点線はIチャンネル測定値の外挿曲線、2点鎖線はQチャンネル測定値の外挿曲線、実線は計算値を示している。 In FIG. 14, (a) is the fundamental wave (FM modulation frequency f m ) component, (b) is the second harmonic component, and (c) is the measurement result regarding the third harmonic component, and the horizontal axis represents the delay time. (Phase φ r ) is represented, and the vertical axis represents the second mixer output power. In FIG. 14, the white circle is the I-channel second mixer output measurement value, the black circle is the Q-channel second mixer output measurement value, the dotted line is the extrapolation curve of the I-channel measurement value, and the two-dot chain line is the Q-channel measurement value. Extrapolated curves and solid lines show calculated values.

図14から、位相φrの変化に応じてドプラ信号が周期的に変化する様子がわかる。また、数16を考慮すると、位相φrを固定してターゲットの深さを変化させてもドプラ信号が周期的に変化することが理解できる。ドプラ信号がターゲットの深さにより周期的に変化するという現象は、表現を代えれば、ターゲットの深さによりドプラ信号の受信感度が変動するということであり、超音波診断装置にこれを応用する場合、受信感度の変動(ドプラ信号の深さ依存性)への対策が望まれる。 From Figure 14, the Doppler signal it can be seen that changes periodically in accordance with the change of the phase phi r. Further, considering Equation 16, it can be understood that the Doppler signal periodically changes even if the depth of the target is changed with the phase φ r fixed. The phenomenon that the Doppler signal periodically changes depending on the depth of the target is that if the expression is changed, the reception sensitivity of the Doppler signal varies depending on the depth of the target, and this is applied to an ultrasonic diagnostic apparatus. Measures against fluctuations in reception sensitivity (depth dependence of Doppler signals) are desired.

図15は、図12に示す超音波診断装置を改良して、ドプラ信号の深さ依存性を回避した超音波診断装置を示している。図12と異なる点は、図15に示す超音波診断装置では、BPF36の出力に対して直交検波方式による検波処理を施し、さらに、BPF38の出力に対して直交検波方式による検波処理を施すことである。以下、図12と異なる部分を中心に図15の超音波診断装置について説明する。   FIG. 15 shows an ultrasonic diagnostic apparatus in which the ultrasonic diagnostic apparatus shown in FIG. 12 is improved to avoid the depth dependence of the Doppler signal. The ultrasonic diagnostic apparatus shown in FIG. 15 is different from FIG. 12 in that the output of the BPF 36 is subjected to detection processing by the quadrature detection method, and further, the output of the BPF 38 is subjected to detection processing by the quadrature detection method. is there. Hereinafter, the ultrasonic diagnostic apparatus of FIG. 15 will be described with a focus on differences from FIG.

図15において、BPF36,38は、数8に示されるfmに関する基本波成分および各高調波成分の中から所望の周波数成分を抽出する。つまり、基本波成分および第n(nは2以上の自然数)高調波成分のうちの少なくとも一つの成分とその付近に存在するドプラ信号とから成る側帯波を抽出する。 In FIG. 15, the BPFs 36 and 38 extract a desired frequency component from the fundamental wave component and each harmonic component related to f m shown in Equation 8. That is, a sideband consisting of at least one of the fundamental wave component and the nth (n is a natural number of 2 or more) harmonic component and a Doppler signal existing in the vicinity thereof is extracted.

そして、BPF36から出力される信号(側帯波信号)に対して、二つのミキサ150,152および二つのローパスフィルタ(LPF)160,162によって同期検波および直交検波を行い、さらに、BPF38から出力される信号(側帯波信号)に対して、二つのミキサ154,156および二つのローパスフィルタ(LPF)164,166によって同期検波および直交検波を行う。   The signal (sideband signal) output from the BPF 36 is subjected to synchronous detection and quadrature detection by the two mixers 150 and 152 and the two low-pass filters (LPF) 160 and 162, and further output from the BPF 38. The signal (sideband signal) is subjected to synchronous detection and quadrature detection by two mixers 154 and 156 and two low-pass filters (LPF) 164 and 166.

ミキサ150,154は、N倍回路25で基本波fmをN倍(Nは自然数)して得られる変調波fmの基本波成分あるいは高調波成分を参照信号とする。一方、ミキサ152,156は、N倍回路25で基本波fmをN倍(Nは自然数)した後、π/2シフト回路180を経由して供給される信号を参照信号とする。π/2シフト回路180は、入力信号の位相をπ/2だけずらす回路である。 Mixer 150, 154 are, is N times the fundamental wave f m N times circuit 25 (N is a natural number) and the reference signal fundamental wave component or a harmonic component of the modulation wave f m obtained by. On the other hand, the mixers 152 and 156 use the signal supplied via the π / 2 shift circuit 180 after the fundamental wave f m is multiplied N times (N is a natural number) by the N multiplication circuit 25 as a reference signal. The π / 2 shift circuit 180 is a circuit that shifts the phase of the input signal by π / 2.

この構成により、BPF36の出力に対して直交検波方式による検波処理が実行され、また、BPF38の出力に対して直交検波方式による検波処理が実行される。なお、ミキサ150,152,154,156を利用する検波処理、および、ミキサ32,34を利用する検波処理は、アナログ回路による検波処理の他、デジタル処理(例えば、サンプリング処理などを含む)による検波処理でもよい。   With this configuration, detection processing by the quadrature detection method is executed on the output of the BPF 36, and detection processing by the quadrature detection method is executed on the output of the BPF 38. Note that the detection processing using the mixers 150, 152, 154, and 156 and the detection processing using the mixers 32 and 34 are detection by digital processing (for example, including sampling processing) in addition to detection processing by an analog circuit. Processing is also acceptable.

そして、2乗和平方根演算部170においてLPF160の出力(I成分)およびLPF162の出力(Q成分)から、これらの2乗和平方根、つまり(I2+Q21/2が演算される。同様に、2乗和平方根演算部172においてLPF164の出力(I成分)およびLPF166の出力(Q成分)から2乗和平方根が演算される。 Then, the square sum square root calculation unit 170 calculates the square sum square, that is, (I 2 + Q 2 ) 1/2 from the output (I component) of the LPF 160 and the output (Q component) of the LPF 162. Similarly, the square sum square root calculation unit 172 calculates the square sum square root from the output (I component) of the LPF 164 and the output (Q component) of the LPF 166.

そして、2乗和平方根演算部170における演算結果がドプラ信号同相成分として出力され、また、2乗和平方根演算部172における演算結果がドプラ信号直交成分として出力される。ドプラ演算処理部40は、ドプラ信号同相成分およびドプラ信号直交成分に基づいてドプラ演算を実行する。   The calculation result in the square sum square root calculation unit 170 is output as a Doppler signal in-phase component, and the calculation result in the square sum square root calculation unit 172 is output as a Doppler signal orthogonal component. The Doppler calculation processing unit 40 performs Doppler calculation based on the Doppler signal in-phase component and the Doppler signal quadrature component.

このように、図15に示す超音波診断装置では、BPF36の出力およびBPF38の出力の各々に対して直交検波方式による検波処理を施すことにより、ドプラ信号の深さ依存性を回避している。その原理は以下のとおりである。   As described above, the ultrasonic diagnostic apparatus shown in FIG. 15 avoids the depth dependence of the Doppler signal by performing the detection processing by the orthogonal detection method on each of the output of the BPF 36 and the output of the BPF 38. The principle is as follows.

まず、BPF36により第2高調波成分が選択されたとすると数14の信号が抽出される。そして、ミキサ150の参照信号を数17のように表現する。

Figure 0004688466
その結果、LPF160の出力は次のように表現される。
Figure 0004688466
First, assuming that the second harmonic component is selected by the BPF 36, the signal of Expression 14 is extracted. Then, the reference signal of the mixer 150 is expressed as in Expression 17.
Figure 0004688466
As a result, the output of the LPF 160 is expressed as follows.
Figure 0004688466

また、ミキサ152の参照信号を数19のように表現する。

Figure 0004688466
その結果、LPF162の出力は次のように表現される。
Figure 0004688466
Further, the reference signal of the mixer 152 is expressed as in Expression 19.
Figure 0004688466
As a result, the output of the LPF 162 is expressed as follows.
Figure 0004688466

さらに、2乗和平方根演算部170において、数18および数20によって得られる二つの信号の2乗和の平方根が以下のように演算される。

Figure 0004688466
数21におけるvD3I(t)には、数16のvD2(t)に含まれていたφmおよびφrが存在しない。したがって、vD3I(t)には、深さ依存性が存在しないことになる。 Further, in the square sum square root calculation unit 170, the square root of the square sum of the two signals obtained by Expression 18 and Expression 20 is calculated as follows.
Figure 0004688466
In v D3I (t) in Equation 21, φ m and φ r included in v D2 (t) in Equation 16 do not exist. Therefore, there is no depth dependence in v D3I (t).

2乗和平方根演算部172の出力も、2乗和平方根演算部170の出力と同様な計算プロセスにより次式のように演算される。

Figure 0004688466
数22におけるvD3Q(t)にもφmおよびφrが含まれず、したがって、深さ依存性が存在しないことになる。 The output of the square sum square root calculation unit 172 is also calculated by the same calculation process as the output of the square sum square root calculation unit 170 as follows.
Figure 0004688466
V D3Q (t) in Equation 22 also does not include φ m and φ r , and therefore there is no depth dependency.

数21および数22はドプラ信号の直交検波出力であり、これら出力に基づいて、ドプラ演算処理部40において、従来から知られている通常の信号処理によりターゲットの速度成分などが演算される。   Equations (21) and (22) are quadrature detection outputs of the Doppler signal, and based on these outputs, the Doppler computation processing unit 40 computes the velocity component of the target and the like by conventional signal processing conventionally known.

以上説明したように、図15の超音波診断装置により、ドプラ信号の深さ依存性を回避することができる。すなわち、ドプラターゲットがどのような深さにあろうとも、ドプラ信号はその深さに依存しない信号として抽出でき、後段の処理において、通常のドプラ信号処理と同様に、FFTおよび表示系により画像として観測ができる。   As described above, the depth dependence of the Doppler signal can be avoided by the ultrasonic diagnostic apparatus of FIG. That is, regardless of the depth of the Doppler target, the Doppler signal can be extracted as a signal that does not depend on the depth, and in the subsequent processing, as in normal Doppler signal processing, an image is obtained by FFT and a display system. Can observe.

図16は、図11に示す超音波診断装置をさらに改良した超音波診断装置を示している。図16に示す超音波診断装置では、図11のN倍回路25に換えて、1倍回路25a、2倍回路25bおよび3倍回路25cの3つの回路が設けられ、また、これら3つの回路の各々に対応して、受信ミキサ30の後段が3系統に分離されている。   FIG. 16 shows an ultrasonic diagnostic apparatus in which the ultrasonic diagnostic apparatus shown in FIG. 11 is further improved. In the ultrasonic diagnostic apparatus shown in FIG. 16, three circuits of a 1 × circuit 25a, a 2 × circuit 25b, and a 3 × circuit 25c are provided instead of the N × circuit 25 of FIG. Corresponding to each, the subsequent stage of the receiving mixer 30 is separated into three systems.

BPF36a,38aは、それぞれ、受信ミキサ30の出力信号から変調波周波数fmに関する基本波成分を抽出して、ミキサ28a,29aへ出力する。ミキサ28a,29aは、それぞれ、1倍回路25aから得られる変調波fmの基本波成分を参照信号として、BPF36a,38aから出力される信号を検波してベースバンド信号に変換する。また、BPF36b,38bは、それぞれ、受信ミキサ30の出力信号から変調波周波数fmに関する2次高調波成分を抽出して、ミキサ28b,29bへ出力する。ミキサ28b,29bは、それぞれ、2倍回路25bから得られる変調波fmの2次高調波成分を参照信号として、BPF36b,38bから出力される信号を検波してベースバンド信号に変換する。さらに、BPF36c,38cは、それぞれ、受信ミキサ30の出力信号から変調波周波数fmに関する3次高調波成分を抽出して、ミキサ28c,29cへ出力する。ミキサ28c,29cは、それぞれ、3倍回路25cから得られる変調波fmの3次高調波成分を参照信号として、BPF36c,38cから出力される信号を検波してベースバンド信号に変換する。 BPF36a, 38a, respectively, to extract the fundamental component about the modulation wave frequency f m from the output signal of the receiving mixer 30, and outputs the mixers 28a, to 29a. Mixers 28a, 29a, respectively, as a reference signal to the fundamental wave component of the modulated wave f m obtained from 1-fold circuit 25a, BPF36a, by detecting a signal outputted from 38a into a baseband signal. Further, BPF36b, 38b, respectively, to extract the second harmonic component about the modulation wave frequency f m from the output signal of the receiving mixer 30, and outputs the mixer 28b, the 29b. Mixers 28b, 29 b, respectively, as a reference signal to the second harmonic component of the modulation wave f m derived from doubled circuits 25b, BPF36b, by detecting a signal output from 38b converts the baseband signal. Furthermore, BPF36c, 38c, respectively, to extract the third-order harmonic component about the modulation wave frequency f m from the output signal of the reception mixer 30, the mixer 28c, and outputs it to 29c. Mixer 28c, 29c, respectively, as a reference signal to the third harmonic component of the modulation wave f m obtained from 3-fold circuit 25c, BPF36c, converted into a baseband signal by detecting the signal output from 38c.

そして、ミキサ28a,28b,28cから出力される信号は、加算器120において加算処理されてドプラ演算処理部40に出力され、また、ミキサ29a,29b,29cから出力される信号は、加算器122において加算処理されてドプラ演算処理部40に出力される。加算器120,122は、各ミキサから出力される信号に対して、重み付けなどを行って加算処理する。   The signals output from the mixers 28a, 28b, and 28c are added in the adder 120 and output to the Doppler arithmetic processing unit 40, and the signals output from the mixers 29a, 29b, and 29c are added to the adder 122. Is added to and output to the Doppler computation processing unit 40. Adders 120 and 122 perform addition processing by weighting the signals output from the mixers.

こうして、図16に示す超音波診断装置では、変調波周波数fmの基本波、2倍の高調波および3倍の高調波からそれぞれ抽出した信号(ドプラ信号)を加算処理してドプラ情報を抽出する。このため、図11に示す超音波診断装置のように、変調波fmの基本波成分および高調波成分のいずれか一つに対応するドプラ信号を利用する場合に比べてSNRを改善することができる。 Thus, the ultrasonic diagnostic apparatus illustrated in FIG. 16, extracts the fundamental wave, Doppler information by adding treated twice the harmonic and 3-fold, respectively extracted signal from harmonic (Doppler signal) of the modulation wave frequency f m To do. Therefore, as in the ultrasonic diagnostic apparatus shown in FIG. 11, to improve the SNR as compared with the case of using the Doppler signal corresponding to one of the fundamental wave component and harmonic components of the modulated wave f m it can.

なお、必要に応じて、図16の1倍回路25a、2倍回路25bおよび3倍回路25cのそれぞれの後段に、位相や電圧調整を行う回路を追加して1倍回路25a、2倍回路25bおよび3倍回路25cのそれぞれから出力される参照信号に対して、位相や電圧の調整を行ってもよい。参照信号の位相や電圧を調整して、最適な検波(例えばSNRの良い検波)が実行されるように設定することができる。   If necessary, a circuit for adjusting the phase and voltage is added to the subsequent stage of each of the 1 × circuit 25a, 2 × circuit 25b, and 3 × circuit 25c in FIG. 16 to add 1 × circuit 25a, 2 × circuit 25b. The phase and voltage may be adjusted for the reference signal output from each of the triple circuit 25c. The phase and voltage of the reference signal can be adjusted so that optimum detection (for example, detection with good SNR) is performed.

また、1倍回路25a、2倍回路25bおよび3倍回路25cに加えて、4倍回路、5倍回路を設けるなど、さらに高次の参照信号を利用してもよい。さらに高次の参照信号を利用する場合、各次元の参照信号に対応して、受信ミキサ30の後段も4系統あるいは5系統などに拡張される。   Further, in addition to the 1 × circuit 25a, the 2 × circuit 25b, and the 3 × circuit 25c, a higher order reference signal such as a 4 × circuit or a 5 × circuit may be used. Further, when a higher-order reference signal is used, the subsequent stage of the reception mixer 30 is also expanded to four systems or five systems in accordance with each dimension of the reference signal.

図17は、図1に示す超音波診断装置を簡略化した超音波診断装置を示している。図17に示す超音波診断装置では、図1のミキサ32,34に換えて一つのミキサ33が設けられ、また、ミキサ28,29に換えて一つのミキサ27が設けられている。つまり、図1の超音波診断装置において、ミキサ32,34およびミキサ28,29で実行されていた直交検波が、図17では単純な検波に簡略化されている。単純な検波に簡略化したことに伴い、図17では、図1で必要とされていたπ/2シフト回路26、BPF38およびLPF39が省略されている。   FIG. 17 shows an ultrasonic diagnostic apparatus in which the ultrasonic diagnostic apparatus shown in FIG. 1 is simplified. In the ultrasonic diagnostic apparatus shown in FIG. 17, one mixer 33 is provided in place of the mixers 32 and 34 in FIG. 1, and one mixer 27 is provided in place of the mixers 28 and 29. That is, the quadrature detection performed in the mixers 32 and 34 and the mixers 28 and 29 in the ultrasonic diagnostic apparatus of FIG. 1 is simplified to simple detection in FIG. With the simplification of the simple detection, the π / 2 shift circuit 26, the BPF 38, and the LPF 39 that are required in FIG. 1 are omitted in FIG.

図17に示す超音波診断装置では、直交検波に換えて単純な検波としたことにより、血流の方向に関する情報が失われるものの、血流の速度計測は直交検波の場合と同程度の計測が可能である。しかも、図1の超音波診断装置に比べて装置の構成を簡略化することができる。なお、図17では、LPF37の後段にウォールモーションフィルタ(WMF)130が設けられ、このWMF130においてクラッタが除去される。   In the ultrasonic diagnostic apparatus shown in FIG. 17, information relating to the direction of blood flow is lost because simple detection is performed instead of quadrature detection, but blood flow velocity measurement is approximately the same as in quadrature detection. Is possible. In addition, the configuration of the apparatus can be simplified as compared with the ultrasonic diagnostic apparatus of FIG. In FIG. 17, a wall motion filter (WMF) 130 is provided after the LPF 37, and clutter is removed in the WMF 130.

以上、本発明の好適な実施形態を説明したが、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。例えば、上述の実施形態では周波数fmの正弦波によるFM変調の例を示したが、周波数fm1の変調波にさらに周波数fm2の変調波でFM変調を施したFMCW送信波を形成してもよい。FM変調を多段で行うことにより、さらに低速の領域においてクラッタの影響を低減してドプラ信号を抽出することができる。 As mentioned above, although preferred embodiment of this invention was described, embodiment mentioned above is only a mere illustration in all the points, and does not limit the scope of the present invention. For example, an example of the FM-modulated by a sine wave of a frequency f m in the above embodiment, by forming the FMCW transmission wave subjected to FM modulation by the modulation wave of further frequency f m2 to the modulation wave frequency f m1 Also good. By performing FM modulation in multiple stages, the influence of clutter can be reduced and a Doppler signal can be extracted in a lower speed region.

また、上述の実施形態では周波数fmの正弦波によるFM変調の例を示したが、正弦波に換えて、例えば図18に示す周期fmの三角波140を利用してもよい。もちろん、三角波として鋸歯状の波を利用してもよい。 Also, although an example of the FM-modulated by a sine wave of a frequency f m in the above embodiment, instead of the sine wave may be utilized triangular wave 140 of the period f m shown in FIG. 18 for example. Of course, a sawtooth wave may be used as the triangular wave.

本発明に係る超音波診断装置の全体構成を示すブロック図である。1 is a block diagram showing an overall configuration of an ultrasonic diagnostic apparatus according to the present invention. FMCW送信波の周波数スペクトルを示す図である。It is a figure which shows the frequency spectrum of a FMCW transmission wave. FMCW受信波の周波数スペクトルを示す図である。It is a figure which shows the frequency spectrum of a FMCW received wave. 受信ミキサ出力の周波数スペクトルを示す図である。It is a figure which shows the frequency spectrum of a receiving mixer output. 受信RF信号における通常のCWでのクラッタの電力スペクトルを示す図である。It is a figure which shows the power spectrum of the clutter in normal CW in a received RF signal. 受信RF信号におけるFMCWでのクラッタの電力スペクトルを示す図である。It is a figure which shows the power spectrum of the clutter in FMCW in a received RF signal. 第1次ベッセル関数を示す図である。It is a figure which shows a primary Bessel function. 横軸を診断の深さ、縦軸をベッセル関数値としてプロットした図である。It is the figure which plotted the horizontal axis as the depth of diagnosis and the vertical axis as the Bessel function value. 基本波および高調波クラッタ電圧の深さ依存性を説明するための図である。It is a figure for demonstrating the depth dependence of a fundamental wave and a harmonic clutter voltage. ドプラ信号電力とクラッタ電力との電力比関係を説明するための図である。It is a figure for demonstrating the power ratio relationship of Doppler signal electric power and clutter electric power. 本発明に係る別の超音波診断装置の全体構成を示す図である。It is a figure which shows the whole structure of another ultrasonic diagnostic apparatus which concerns on this invention. 本発明に係るさらに別の超音波診断装置の全体構成を示す図である。It is a figure which shows the whole structure of another ultrasonic diagnostic apparatus which concerns on this invention. 第1ミキサ出力におけるドプラ信号を説明するための図である。It is a figure for demonstrating the Doppler signal in a 1st mixer output. 第2ミキサ出力電力の参照波遅延時間依存性を説明するための図である。It is a figure for demonstrating the reference wave delay time dependence of 2nd mixer output electric power. ドプラ信号の深さ依存性を回避した超音波診断装置を説明するための図である。It is a figure for demonstrating the ultrasonic diagnosing device which avoided the depth dependence of the Doppler signal. 本発明に係るさらに別の超音波診断装置の全体構成を示す図である。It is a figure which shows the whole structure of another ultrasonic diagnostic apparatus which concerns on this invention. 本発明に係るさらに別の超音波診断装置の全体構成を示す図である。It is a figure which shows the whole structure of another ultrasonic diagnostic apparatus which concerns on this invention. 三角波を説明するための図である。It is a figure for demonstrating a triangular wave.

符号の説明Explanation of symbols

20 FM変調器、22 RF波発振器、24 FM変調波発振器、26 π/2シフト回路、30 受信ミキサ、40 ドプラ演算処理部。   20 FM modulator, 22 RF wave oscillator, 24 FM modulated wave oscillator, 26 π / 2 shift circuit, 30 receiving mixer, 40 Doppler arithmetic processing unit.

Claims (17)

搬送波信号を変調波信号を用いて変調処理することにより、周波数を周期的に変化させた連続波の変調送信信号を生成する送信部と、
前記変調送信信号の供給によって生体に対して超音波を送波し、生体からの超音波を受波して受信信号を出力する送受波器と、
前記受信信号を前記変調送信信号を用いて復調し、これにより復調信号を得る復調手段と、
前記復調信号から、前記変調波信号を基本波とした場合における基本波成分および第n(nは2以上の自然数)高調波成分のうちの少なくとも一つの成分の付近に存在するドプラ情報を抽出するドプラ情報抽出手段と、
を有することを特徴とする超音波診断装置。
A modulation unit that modulates a carrier wave signal using a modulated wave signal, thereby generating a continuous wave modulated transmission signal whose frequency is periodically changed ; and
A transmitter / receiver for transmitting an ultrasonic wave to a living body by supplying the modulated transmission signal, receiving an ultrasonic wave from the living body and outputting a reception signal;
Demodulating means for demodulating the received signal using the modulated transmission signal, thereby obtaining a demodulated signal;
Doppler information existing in the vicinity of at least one of the fundamental wave component and the n-th (n is a natural number of 2 or more) harmonic component when the modulated wave signal is a fundamental wave is extracted from the demodulated signal. Doppler information extraction means;
An ultrasonic diagnostic apparatus comprising:
請求項1に記載の超音波診断装置において、
前記復調手段は、前記変調送信信号を参照信号として前記受信信号を検波処理する検波回路を含む、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 1,
The demodulating means includes a detection circuit that detects the received signal using the modulated transmission signal as a reference signal.
An ultrasonic diagnostic apparatus.
請求項1に記載の超音波診断装置において、
前記復調手段は、前記変調送信信号を参照信号として前記受信信号を直交検波処理する直交検波回路を含む、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 1,
The demodulating means includes a quadrature detection circuit that performs quadrature detection processing on the received signal using the modulated transmission signal as a reference signal.
An ultrasonic diagnostic apparatus.
請求項1に記載の超音波診断装置において、
前記ドプラ情報抽出手段は、前記少なくとも一つの成分およびその付近に存在するドプラ信号から成る側帯波を抽出するフィルタ部を含む、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 1,
The Doppler information extraction unit includes a filter unit that extracts a sideband composed of the at least one component and a Doppler signal existing in the vicinity thereof.
An ultrasonic diagnostic apparatus.
請求項4に記載の超音波診断装置において、
前記ドプラ情報抽出手段は、前記側帯波から前記ドプラ信号を抽出するドプラ信号抽出部を含む、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 4,
The Doppler information extraction unit includes a Doppler signal extraction unit that extracts the Doppler signal from the sideband.
An ultrasonic diagnostic apparatus.
請求項5に記載の超音波診断装置において、
前記ドプラ信号抽出部は、前記変調波信号の周波数をN倍(Nは自然数)した信号のうちの少なくとも一つを参照信号として、前記側帯波を検波処理する検波回路を含む、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 5,
The Doppler signal extraction unit includes a detection circuit that detects the sideband using at least one of signals obtained by multiplying the frequency of the modulated wave signal by N (N is a natural number) as a reference signal.
An ultrasonic diagnostic apparatus.
請求項2に記載の超音波診断装置において、
前記復調手段に与える参照信号に対して遅延処理を施す遅延回路を含む、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 2,
Including a delay circuit that performs a delay process on a reference signal to be supplied to the demodulation means
An ultrasonic diagnostic apparatus.
請求項2に記載の超音波診断装置において、
前記復調手段へ出力される受信信号に対して遅延処理を施す遅延回路を含む、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 2,
A delay circuit that performs a delay process on the received signal output to the demodulating means;
An ultrasonic diagnostic apparatus.
請求項6に記載の超音波診断装置において、
前記ドプラ信号抽出部に与える参照信号に対して遅延処理を施す遅延回路を含む、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 6,
A delay circuit that performs a delay process on a reference signal to be provided to the Doppler signal extraction unit;
An ultrasonic diagnostic apparatus.
請求項6に記載の超音波診断装置において、
前記ドプラ信号抽出部へ出力される側帯波に対して遅延処理を施す遅延回路を含む、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 6,
Including a delay circuit that performs a delay process on the sideband output to the Doppler signal extraction unit,
An ultrasonic diagnostic apparatus.
請求項1に記載の超音波診断装置において、
前記ドプラ情報抽出手段は、前記変調波信号を基本波とした場合における基本波成分および第n(nは2以上の自然数)高調波成分のうちの複数成分の各々からドプラ信号を抽出し、抽出した複数のドプラ信号を加算処理する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 1,
The Doppler information extraction unit extracts and extracts a Doppler signal from each of a plurality of components of a fundamental wave component and an nth (n is a natural number of 2 or more) harmonic component when the modulated wave signal is a fundamental wave. Add multiple Doppler signals
An ultrasonic diagnostic apparatus.
請求項1から11のいずれか1項に記載の超音波診断装置において、
前記変調波信号は正弦波である、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to any one of claims 1 to 11,
The modulated wave signal is a sine wave;
An ultrasonic diagnostic apparatus.
請求項1から11のいずれか1項に記載の超音波診断装置において、
前記変調波信号は三角波または鋸歯状波である、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to any one of claims 1 to 11,
The modulated wave signal is a triangular wave or a sawtooth wave,
An ultrasonic diagnostic apparatus.
周波数fの搬送波信号を周波数fの変調波信号を用いて周波数変調処理することにより、変調送信信号を生成する送信部と、
前記変調送信信号の供給によって生体に対して超音波を送波し、生体からの反射波を受波して受信信号を出力する送受波器と、
前記受信信号に対して、前記変調送信信号を参照信号として直交検波処理を行い、これにより復調信号を得る検波回路と、
前記復調信号から、前記変調波信号を基本波とした場合における基本波成分および第n(nは2以上の自然数)高調波成分のうちの少なくとも一つの成分とその付近に存在するドプラ信号とから成る側帯波を抽出するバンドパスフィルタと、
前記側帯波から前記ドプラ信号を抽出するドプラ信号抽出部と、
前記ドプラ信号から、運動体に関するドプラ情報を取得するドプラ情報取得部と、
を有することを特徴とする超音波診断装置。
By frequency modulation processing by using the modulated wave signal of the frequency f m of the carrier signal of a frequency f 0, and a transmission unit for generating a modulated transmission signal,
A transmitter / receiver for transmitting ultrasonic waves to a living body by supplying the modulated transmission signal, receiving a reflected wave from the living body and outputting a reception signal;
A detection circuit that performs quadrature detection processing on the received signal using the modulated transmission signal as a reference signal, thereby obtaining a demodulated signal, and
From the demodulated signal, from at least one of the fundamental wave component and the n-th (n is a natural number of 2 or more) harmonic component when the modulated wave signal is a fundamental wave, and a Doppler signal present in the vicinity thereof A bandpass filter for extracting sidebands consisting of:
A Doppler signal extraction unit for extracting the Doppler signal from the sideband;
A Doppler information acquisition unit that acquires Doppler information about the moving body from the Doppler signal;
An ultrasonic diagnostic apparatus comprising:
請求項14に記載の超音波診断装置において、
前記変調波信号として、周波数fm1の変調波信号をさらに周波数fm2の変調波信号を用いて周波数変調処理することにより得られる多重変調波信号を用いる、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 14,
As the modulation wave signal, a multiplex modulation wave signal obtained by further frequency-modulating the modulation wave signal of frequency f m1 using the modulation wave signal of frequency f m2 is used.
An ultrasonic diagnostic apparatus.
周波数fの搬送波信号を周波数fの変調波信号を用いて周波数変調処理することにより、変調送信信号を生成する送信部と、
前記変調送信信号の供給によって生体に対して超音波を送波し、生体からの反射波を受波して受信信号を出力する送受波器と、
前記受信信号に対して、前記変調送信信号を参照信号として直交検波処理を行い、これにより同相成分復調信号および直交成分復調信号を得る検波回路と、
前記同相成分復調信号および直交成分復調信号の各々から、前記変調波信号を基本波とした場合における基本波成分および第n(nは2以上の自然数)高調波成分のうちの少なくとも一つの成分とその付近に存在するドプラ信号とから成る側帯波を抽出するバンドパスフィルタと、
前記同相成分復調信号の側帯波および前記直交成分復調信号の側帯波から、前記ドプラ信号を抽出するドプラ信号抽出部と、
前記ドプラ信号から、運動体に関するドプラ情報を取得するドプラ情報取得部と、
を有し、
前記ドプラ信号抽出部は、前記変調波信号の周波数をN倍(Nは自然数)した信号のうちの少なくとも一つを参照信号として、前記同相成分復調信号から抽出された側帯波を直交検波処理する同相成分検波回路および前記直交成分復調信号から抽出された側帯波を直交検波処理する直交成分検波回路を含む、
ことを特徴とする超音波診断装置。
By frequency modulation processing by using the modulated wave signal of the frequency f m of the carrier signal of a frequency f 0, and a transmission unit for generating a modulated transmission signal,
A transmitter / receiver for transmitting ultrasonic waves to a living body by supplying the modulated transmission signal, receiving a reflected wave from the living body and outputting a reception signal;
A detection circuit that performs quadrature detection processing on the received signal using the modulated transmission signal as a reference signal, thereby obtaining an in-phase component demodulated signal and a quadrature component demodulated signal;
From each of the in-phase component demodulated signal and the quadrature component demodulated signal, at least one of a fundamental wave component and an nth (n is a natural number of 2 or more) harmonic component when the modulated wave signal is a fundamental wave; A bandpass filter that extracts sidebands composed of Doppler signals present in the vicinity thereof;
A Doppler signal extraction unit that extracts the Doppler signal from the sideband wave of the in-phase component demodulated signal and the sideband wave of the quadrature component demodulated signal;
A Doppler information acquisition unit that acquires Doppler information about the moving body from the Doppler signal;
Have
The Doppler signal extraction unit performs quadrature detection processing on the sideband extracted from the in-phase component demodulated signal, using at least one of the signals obtained by multiplying the frequency of the modulated wave signal by N (N is a natural number) as a reference signal. Including an in-phase component detection circuit and a quadrature component detection circuit for performing quadrature detection processing on the sideband wave extracted from the quadrature component demodulated signal,
An ultrasonic diagnostic apparatus.
請求項16に記載の超音波診断装置において、
前記ドプラ信号抽出部は、前記同相成分検波回路から出力される二つの信号に対する二乗和の平方根として得られるドプラ信号同相成分、および、前記直交成分検波回路から出力される二つの信号に対する二乗和の平方根として得られるドプラ信号直交成分、の二つのドプラ信号成分を抽出する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 16, wherein
The Doppler signal extraction unit obtains the Doppler signal in-phase component obtained as the square root of the sum of squares of the two signals output from the in-phase component detection circuit, and the square sum of the two signals output from the quadrature component detection circuit. Extract the two Doppler signal components of the Doppler signal orthogonal component obtained as the square root,
An ultrasonic diagnostic apparatus.
JP2004295899A 2003-10-10 2004-10-08 Ultrasonic diagnostic equipment Expired - Fee Related JP4688466B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2004295899A JP4688466B2 (en) 2003-10-10 2004-10-08 Ultrasonic diagnostic equipment

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
JP2003352636 2003-10-10
JP2003352636 2003-10-10
JP2004034828 2004-02-12
JP2004034828 2004-02-12
JP2004295899A JP4688466B2 (en) 2003-10-10 2004-10-08 Ultrasonic diagnostic equipment

Publications (2)

Publication Number Publication Date
JP2005253949A JP2005253949A (en) 2005-09-22
JP4688466B2 true JP4688466B2 (en) 2011-05-25

Family

ID=35080244

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004295899A Expired - Fee Related JP4688466B2 (en) 2003-10-10 2004-10-08 Ultrasonic diagnostic equipment

Country Status (1)

Country Link
JP (1) JP4688466B2 (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4787683B2 (en) * 2006-06-15 2011-10-05 日立アロカメディカル株式会社 Ultrasonic diagnostic equipment
JP5459963B2 (en) * 2007-04-27 2014-04-02 日立アロカメディカル株式会社 Ultrasonic diagnostic equipment
US8100832B2 (en) 2007-04-27 2012-01-24 Hitachi Aloka Medical, Ltd. Ultrasound diagnostic apparatus
JP5297082B2 (en) * 2008-05-12 2013-09-25 日立アロカメディカル株式会社 Ultrasonic diagnostic equipment
JP5629052B2 (en) 2008-06-03 2014-11-19 日立アロカメディカル株式会社 Ultrasonic diagnostic equipment
JP5654198B2 (en) 2008-06-16 2015-01-14 日立アロカメディカル株式会社 Ultrasonic diagnostic equipment
JP5235110B2 (en) * 2008-07-16 2013-07-10 日立アロカメディカル株式会社 Ultrasonic diagnostic equipment
JP5325503B2 (en) * 2008-08-27 2013-10-23 日立アロカメディカル株式会社 Ultrasonic diagnostic equipment
JP5416499B2 (en) * 2008-09-03 2014-02-12 富士フイルム株式会社 Ultrasonic diagnostic equipment
JP2010075289A (en) * 2008-09-24 2010-04-08 Aloka Co Ltd Ultrasonic diagnostic apparatus
JP7025886B2 (en) * 2017-10-19 2022-02-25 フクダ電子株式会社 Oscillator drive circuit and ultrasonic diagnostic equipment

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0348789A (en) * 1989-07-17 1991-03-01 Yokogawa Medical Syst Ltd Cw doppler device

Also Published As

Publication number Publication date
JP2005253949A (en) 2005-09-22

Similar Documents

Publication Publication Date Title
US8323200B2 (en) Ultrasound diagnostic apparatus
JP5459963B2 (en) Ultrasonic diagnostic equipment
US11123044B2 (en) Signal processing device, ultrasonic diagnostic apparatus, and method
US8343055B2 (en) Ultrasound diagnostic apparatus
US8100832B2 (en) Ultrasound diagnostic apparatus
JP4688466B2 (en) Ultrasonic diagnostic equipment
JP4931910B2 (en) Ultrasonic imaging device
JP2007330541A (en) Ultrasonic diagnostic equipment
JP2006288974A (en) Ultrasonic diagnostic equipment
JP4787683B2 (en) Ultrasonic diagnostic equipment
JP5308748B2 (en) Ultrasonic diagnostic equipment
JP2011240006A (en) Ultrasonograph
JP2009261749A (en) Ultrasonic diagnostic apparatus
JP2011234846A (en) Ultrasonograph
JP2009297444A (en) Ultrasonic diagnostic apparatus
JP2000237187A (en) Ultrasonic diagnostic device
JP5202138B2 (en) Ultrasonic diagnostic equipment
JP5297082B2 (en) Ultrasonic diagnostic equipment
JP2011036599A (en) Ultrasonograph
JP5325503B2 (en) Ultrasonic diagnostic equipment
JP2011036506A (en) Ultrasonograph
JP2010194109A (en) Ultrasonic diagnostic apparatus
JP2009247365A (en) Ultrasonic diagnostic system
JP2009261827A (en) Ultrasonic diagnostic apparatus
JP2011045603A (en) Ultrasonic diagnostic apparatus

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070925

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20101109

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110105

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

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110215

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20140225

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees