JPH05340952A - ドップラ流速計 - Google Patents

Abstract

(57)【要約】 （修正有） 【目的】超音波を反射する粒子を含む流体の流速を測定
するためのドップラ流速計を提供する。 【構成】送受信器２０２に接続された超音波変換器２０
１は、送信パルスを送り出す。送信パルスとその発生装
置とは、２つの分離したチャンネルを有し、受信された
反射信号の処理を行なう。回路２０３によって構成され
る上側チャンネルでは、反射信号は送信パルスの２つの
周波数帯域の中で処理される。このチャンネルの出力の
ところに、流体路の中の解析される点での流速の平均値
に対応する出力信号がえられる。下側チャンネルは、２
個の直角位相復調器２０４，２０５と、低域フイルタ２
０６，２０７と、積分器および保持回路２０８，２０９
と、一定符号消去器２１１，２１２とを有する。コンピ
ュータ装置２１３の出力に、信号対のシーケンスが伝送
され、コンピュータ装置２１４の入力の信号対の各信号
は音響的に再生する装置２１５，２１６に送られる。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は超音波を反射する粒子を
含む流体の流速を測定するためのドップラ測定装置に関
する。

【０００２】

【従来の技術】この装置は下記の部品を有する。（イ）
予め定められたパルス繰り返し速度でそれに供給される
送信パルスに対応して応答する超音波パルスで流体を照
射し、かつ、流体内の粒子によって反射された反射信号
を受け取り、かつ、対応する反射信号を発生する、超音
波変換器。（ロ）送信パルスを発生し、かつ、この送信
パルスによって超音波変換器が励振されて超音波パルス
が放射される、超音波変換器に接続された送信器。
（ハ）送信された波の第１パルスおよび第２パルスに応
答して流体路内の設定点の中にある粒子によって反射さ
れた少なくとも２つの異なる反射波に対応する反射信号
を受け取りかつ処理するために超音波変換器に接続さ
れ、かつ、反射信号の隣接しているがしかし分離した周
波数帯域が別々の信号処理チヤンネル内でそれぞれ処理
される、受信器。（ニ）受信器の出力に接続され、か
つ、受信器でえられたドップラ情報に基づく流速に対応
する出力信号がえられる、評価装置。

【０００３】パルス・ドップラ法において、超音波を反
射する粒子を含む流体、例えば血液、の流速は、予め設
定された体積要素の中で測定される。この目的のため
に、一連の超音波パルスが変換器によって送信される。
パルスを送信してから時間τが経過する毎に、短い時間
の間、受信ゲートが開かれる。τは下記の式で定められ
る。

【０００４】

【数１】 τ＝２ｄ／ｃ， （１）

【０００５】ここで、ｄは体積要素までの距離であり、
そしてｃは音速である。このようにして、短い受信信号
の連続体がえられる。これらの信号が適切な電子装置に
よって処理され、連続したドップラ信号がえられる。こ
のドップラ信号の周波数は、よく知られた下記の公式に
より、流速ｖを与える。

【０００６】

【数２】 ｖ＝−（ｆ_D ．ｃ）／（２ｆ₁ cos θ） （２）

【０００７】ここで、ｆ_D はドップラ周波数、ｃは音
速、ｆ₁ は送信された超音波信号の周波数、θは音波の
伝播方向と速度の方向との間の角度である。この式にお
いて、ドップラ周波数ｆ_D は受信された信号の周波数と
送信された周波数との間の差であることとに注意しなけ
ればならない。

【０００８】全体的にいえば、ドップラ信号の高速フー
リエ変換が実行される。このようにして、流速の分布が
時間の関数として決定される。

【０００９】パルス・ドップラ法において、反射信号が
周期的に解析される、すなわち、周期的な静止点でのみ
測定される。したがって、ナイキストの定理により、一
義的に測定可能な最大ドップラ周波数が限定され、した
がって、一義的に測定可能な最大速度が限定される。

【００１０】一義的に則定可能な最大流速の値は、下記
の方程式で決定することができる。

【００１１】

【数３】 ｜ｖ_max ｜＝（ｃ² ）／（８ｆ₁d coｓθ） （３）

【００１２】ここで、ｃは音速、ｆ₁ は超音波周波数、
ｄは深さ、θは音波の伝播方向と速度の方向との間の角
度である。この方程式を用いて、例えば、ｄ＝１５ｃ
ｍ、 θ＝０、ｆ₁ ＝３ＭＨ、ｃ＝１５４０ｍ／秒の場
合、値ｖ_max ＝０．６６ｍ／秒を計算することができ
る。実際には、さらに大きな速度のえられることがしば
しばある。これらのさらに大きな速度は曲げられる
（「偽」ひずみ）、すなわち、これらはより小さな速度
として、または反対方向の速度として、表される。これ
は、この方法の重大な欠点である。

【００１３】パルス・ドップラ法で一義的に測定するこ
とができる点流速の測定領域を拡大するために、下記の
測定量が採用される。けれども、これらの測定量はすべ
て一定の欠点を有する。

【００１４】１． ｆ₁ より小さな放射周波数の採用 方程式（３）に示されているように、ｖ_max はより大き
くなる。しかし、この方法の空間分解能と感度は小さく
なる。それは、より小さな周波数では、より少ないエネ
ルギが再分散されるからである。したがって、実際に
は、用いられる低い周波数の限界は約２ＭHzである。

【００１５】２． 基線の移動 通常の測定領域は−ｖ_max と＋ｖ_max の間にある。評価
装置を単に変更することにより、この測定領域を移動さ
せることができる。例えば、０と２ｖ_max との間で測定
することができる。もし正の速度のみがえられる場合に
は、このことは有用である。この方法では、移動のみが
えられ、測定領域を拡大することはえられない。

【００１６】３． パルス継起周波数の増大 パルス系列の周波数の増大により、したがって、解析速
度が正規のレベルを超えることにより、測定可能最大速
度の増大がえられる。この方法の欠点は、要求されない
体積要素からの反射信号、例えば要求された体積要素の
外側からの反射信号、がえられることである。この要求
されない体積要素は要求された体積要素よりも変換器の
近くにあることがあり、従って、その反射信号の方が強
くて、要求された反射信号の受信が妨げられる。

【００１７】上記の要点１から要点３に示された方法は
すべて実際に用いられているが、これらのいずれも完全
に満足な解決法を表していないことが示されている。

【００１８】その後、第４の方法が知られるようになっ
たが、この方法は実際には用いられていない。 ４． トラッキング・ドップラ法 この方法では、速度がいつナイキスト限界を越えたかが
決定され、そしてそれから、方程式（３）によって定め
られる速度ｖ_max の倍数が加算される。この方法の主な
問題点は、ナイキスト領域の１つが内部に存在すること
である。したがって、この方法は実際には用いられてい
ない。

【００１９】流速の平均値ｖ_a を決定するために、いわ
ゆる２周波数法が知られている。この方法により、流速
を測定するのにパルス・ドップラ法を用いた場合、測定
領域を大きくすることができる。（米国特許明細書第
４，５３４，３５７号を参照。）このよく知られた２周
波数法では、隣接しているがしかし分離した反射信号周
波数帯域が、分離した信号処理チヤンネルの中で処理さ
れる。この方法によって定めることができる最大の値ｖ
_{a max} は下記の式によって与えられる。

【００２０】

【数４】 ｖ_{a max} ＝（ｃ² ）／（８（ｆ₂ −ｆ₁ ）ｄ cos θ） （４）

【００２１】ここで、ｆ₁ およびｆ₂ は周波数帯域の平
均周波数である。その他の記号は方程式（３）と同じで
ある。

【００２２】（ｆ₂ −ｆ₁ ）＜ｆ₀ であるから、２周波
数法の測定可能な最大速度ｖ_{a max}は、方程式（３）に
よって定められるｖ_max より大きい。

【００２３】すぐ前記で説明した２周波数法と前記で説
明したまた別のパルス・ドップラ法との間に、下記の相
異が存在することを確認することは重要である。

【００２４】前記２周波数法を用いる時、平均流速値は
より少数個の測定点、典型的には４個ないし８個の測定
点、から決定される。したがって、この方法は、カラー
・フロー・マップをトレースするのに特に適切である。

【００２５】前記で説明したまた別のパルス・ドップラ
法では、スペクトルは比較的多数個の測定点、典型的に
は６４個ないし２５６個の測定点、で決定される。ここ
ではしたがって、速度分布に関する大幅に詳細な情報が
時間の関数としてえられるが、より多くの測定点が必要
であり、そのために、より多くの時間が必要である。し
たがって、これらの方法は単一体積要素の中の流速を測
定するのに用いられる。

【００２６】

【発明が解決しようとする課題】本発明の目的は前記方
式のドップラ計をうることであり、このドップラ計によ
り、流速の時間分布の測定を実行することができ、およ
び流速が一義的に測定される測定領域が著しく広がる。

【００２７】

【課題を解決するための手段】この問題点を解決するた
めに、本発明による前記方式のドップラ計は、 （ホ） 受信器が（ホ１） 信号処理チャンネルのそれ
ぞれにそなえられ、かつ、各周波数帯域中で超音波変換
器から受け取った反射信号の直角位相の復調を実行し各
処理チャンネルの出力に、測定された複素値（Ｐ
（０）、Ｐ（Δｔ）、Ｐ（２Δｔ）、…）を定める１対
の直角位相信号がその都度えられる装置と、（ホ２）
両方の信号チャンネルから来る直角位相信号を処理し、
かつ、流体路の中の解析される点における平均流速値に
対応する第１出力信号をうることができる装置と、（ホ
３） 第１出力信号と信号処理チャンネルのうち１つの
チャンネルの出力の直角位相信号とを論理的に相関さ
せ、かつ、決定された平均流速値および測定された測定
値と一致する計算された複素値（Ｐ（Δｔ／２）、Ｐ
（３Δｔ／２）、…）を定める信号をこの論理的相関に
よりうることができる装置と、（ホ４） 測定された複
素値に対応する信号および計算された複素値に対応する
信号（Ｐ（０）、Ｐ（Δｔ／２）、Ｐ（Δｔ）、Ｐ（３
Δｔ／２）、Ｐ（２Δｔ）、…）から成る信号シーケン
スを評価し、かつ、この評価によって瞬間流速値に関す
る情報がえられる装置と、を有することを特徴とする。

【００２８】

【作用】前記種類の流速計を用いた本発明による解決法
により、流速の測定領域を増大するために今日まで知ら
れている方法の前記欠点をなくすることができる。

【００２９】１つの好ましい実施例では、送信パルスの
周期的シーケンスのシーケンス・スペクトルが２つの隣
接しているがしかし分離した周波数帯域から成るよう
に、送信器が構成される。

【００３０】信号シーケンスの評価のための装置は、測
定された複素値に対応する信号と計算された複素値に対
応する信号（Ｐ（０）、Ｐ（Δｔ／２）、Ｐ（Δｔ）、
Ｐ（３Δｔ／２）、Ｐ（２Δｔ）、…）から成る信号シ
ーケンスの解析を実行し、かつ、信号シーケンスの中の
信号によって定められる指標の位相変化に対応し、した
がって流速の瞬間値に対応する振幅を有する第２出力信
号をこれらの解析により生ずる装置を、有することが好
ましい。

【００３１】

【実施例】本発明のいくつかの実施例を添付図面を参照
して説明する。

【００３２】図１は、本発明によるドップラ速度計の原
理構造を示した非常に簡単化されたブロック線図であ
る。送受信器２０２に接続された超音波変換器２０１
は、送信パルスを送り出す。この送信パルスのスペクト
ルは、隣接しているがしかし分離した２個の周波数帯域
で構成されている。これらの送信パルスとその発生装置
とは、図７から図１７までの図面で詳細に説明される。
図１の装置は２つの分離したチャンネルを有し、これら
のチャンネルのおのおのは受信された反射信号の処理を
行なう。回路２０３によって構成される上側チャンネル
では、反射信号は送信パルスの２つの周波数帯域の中で
処理される。このチャンネルの出力のところに、流体路
の中の解析される点での流速の平均値に対応する出力信
号がえられる。この平均値は方程式（４）で与えられた
上限値ｖ_{a max} を有する。したがって、回路２０３を用
いて、従来のドップラ法に比べて、より大きな測定領域
において流速の平均値を測定することができる。流体路
の中の解析される点での流速の平均値ｖ_a を決定するた
めに、４個ないし８個（これは十分過ぎる数である）の
送信パルスの反射信号が通常の方法で処理される。

【００３３】図１の下側チャンネルは、周波数ｆ₁ の従
来のパルス・ドップラ・チャンネルであって、２個の直
角位相復調器２０４および２０５と、逐次低域フイルタ
２０６，２０７と、積分器および保持回路２０８，２０
９と、逐次一定符号消去器２１１，２１２とを有する。
２個の直角位相復調器２０４および２０５は、信号をsi
n ω₁t および cos ω₁t の形の信号と混合する。積
分器および保持回路２０８，２０９は、要求された時刻
において、信号を解析する。逐次一定符号消去器２１
１，２１２は、定常反射信号およびほぼ定常な反射信号
を消去する。一定符号消去器２１１の出力のところに、
直角位相信号 a(t) が伝送され、そして直角位相信号 b
(t) が一定符号消去器２１２の出力から伝送される。信
号 a(t) および b(t) は連続した信号である。けれど
も、下記の解析のために、信号を決定しなければならな
い時刻におけるこれらの信号の値のみを、考察しなけれ
ばならない。これらの時刻は、０、Δｔ、２Δｔであ
り、したがって、信号は次の値を有する。

【００３４】

【数５】 ａ（０），ａ（Δｔ），ａ（２Δｔ），ａ（３Δｔ）…… （５ａ） ｂ（０），ｂ（Δｔ），ｂ（２Δｔ），ｂ（３Δｔ）…… （５ｂ）

【００３５】ここで、Δｔは送信されたパルスの間の時
間間隔であり、ａ（ｋΔｔ）およびｂ（ｋΔｔ）は測定
された値である。これらの値は、複素信号の実数部およ
び虚数部であると考えることができる。

【００３６】

【数６】Ｐ（ｋΔｔ）＝ａ（ｋΔｔ）＋ｊｂ（ｋΔｔ）

【００３７】このように、各時刻における逐次信号の対
ａ（ｋΔｔ）、ｂ（ｋΔｔ）は、測定された複素数値Ｐ
（０）、Ｐ（Δｔ）、Ｐ（２Δｔ）、Ｐ（３Δｔ）を定
める。

【００３８】速度の値は、位相変化Ｐからうることがで
きる。この位相変化Ｐは、方程式（３）により、時間と
共により大きな正規速度限界ｖ_max を有する。

【００３９】本発明により、コンピュータ装置２０３の
出力信号と、信号ａ（ｋΔｔ）およびｂ（ｋΔｔ）が供
給される一定符号消去器２１１および２１２の出力が、
コンピュータ装置２１３の入力に供給される。これらの
信号はコンピュータ装置２１３の中で論理的に相関が行
なわれ、補間法により信号対の計算値ａ′（ｋΔｔ／
２）、ｂ′ (ｋΔｔ／２）がえられる。これらの信号対
の計算値は複素数値の計算値Ｐ（Δｔ／２）、Ｐ（３Δ
ｔ／２）、Ｐ（５Δｔ／２）、を定める。これらの値に
対し、次の方程式が成り立つ。

【００４０】

【数７】 Ｐ（ｋΔｔ／２）＝ａ′（ｋΔｔ／２）＋ｊｂ′（ｋΔｔ／２）

【００４１】値ａ′（ｋΔｔ／２）、ｂ′（ｋΔｔ／
２）は、コンピュータ装置２０３の出力信号によって決
定された流速の平均値ｖ_a と一致し、かつ、一定符号消
去器２１１，２１２の出力における信号の測定値ａ（ｋ
Δｔ）、ｂ（ｋΔｔ）と一致する。測定された値ａ（ｋ
Δｔ）、ｂ（ｋΔｔ）が対応する信号の引き続く対の間
で、例えばコンピュータ装置２１３の中で、補間法によ
り計算された値に対応する１対の信号ａ′（ｋΔｔ／
２）、ｂ′（ｋΔｔ／２）をうることができる。コンピ
ュータ装置２１３の中で、２個またはさらに多数個のこ
れらの補間値をまた、逐次の測定値の間で計算すること
ができる。もし２個の補間値が計算されるならば、その
時には信号の対応する対はａ′（ｋΔｔ／３）、ｂ′
（ｋΔｔ／３）である。もしｎ個の補間値が計算される
ならば、その場合には信号の対応する対はａ′（ｋΔｔ
／〔ｎ＋１〕）、ｂ′（ｋΔｔ／〔ｎ＋１〕）である。
ここで、ｋおよびｎは１，２，３，などの値である。

【００４２】コンピュータ装置２１３の出力に、信号対
のシーケンスが伝送される。この信号対は、測定された
信号ａ（ｋΔｔ），ｂ（ｋΔｔ）の対と、計算された信
号ａ′（ｋΔｔ／２）、ｂ′（ｋΔｔ／２）の対とを有
する。コンピュータ装置２１３の出力の信号シーケンス
の中の信号対のおのおのは、コンピュータ装置２１４に
供給される。コンピュータ装置２１４の中で、この信号
シーケンスの解析、例えばフーリエ変換による解析、が
実行される。この解析の結果えられる出力信号の振幅
は、信号シーケンスの信号によって定められる指標の位
相変化に対応する、したがって、流速の瞬間値に対応す
る。

【００４３】コンピュータ装置２１４の中で解析が実行
される時、６４個から１２８個までの反射信号の中に含
まれている情報が考慮される。反射信号のこの数はま
た、コンピュータ装置２０３の中で流速の平均値ｖ_a を
計算する際に情報量として考慮された反射信号の数（４
個ないし８個）よりも大幅に大きい。

【００４４】コンピュータ装置２１４の入力の信号対の
各信号はまた、デイジタル／アナログ変換器によって、
信号を音響的に再生する装置２１５または２１６に送ら
れる。

【００４５】測定された値ａ（ｋΔｔ）、ｂ（ｋΔｔ）
と計算された値ａ′（ｋΔｔ／２）、ｂ′（ｋΔｔ／
２）とに対応する信号を有する信号シーケンスの中の信
号は、測定された値ａ（ｋΔｔ）、ｂ（ｋΔｔ）に対応
する信号よりも短い時間間隔をそれらの間に有する。

【００４６】図２および図３の簡単な例を使って、コン
ピュータ装置２１３の機能を説明する。図２には、信号
Ｐ（ｋΔｔ）が示されており、これらの信号は測定点の
間にΔφ＝＋２７０°の位相変化を有する。けれども、
解析の周波数は非常に小さいから、これはΔφ＝−９０
°の位相変化として測定される。これは偽信号歪みの典
型的な例である。すなわち、測定された位相の大きさと
符号の両方が偽であり、したがって、測定された速度が
偽である。コンピュータ装置２０３の出力信号がコンピ
ュータ装置２１３に送られる。この信号は流速の平均値
ｖに対応する。この平均値から、解析値の間の位相変化
Δφは＋２７０°であることを決定することができる。
したがって、信号の測定された値の間で、計算値を補間
することができる。

【００４７】図３の点線で示された値Ｐ（Δｔ／２）、
Ｐ（３Δｔ／２）、Ｐ（５Δｔ／２）は補間された値で
ある。２周波数法からの平均速度情報を用いて、測定さ
れた値の間で値を補間することが、ここで説明される本
発明の本質的な特徴である。

【００４８】複素信号の補間が極座標表示で行なわれ
る。すなわち、大きさと位相が補間される。図３から、
信号はこの場合正しく再生されていることがわかる。す
なわち、単一値の間では位相変化は＋１３５°であり、
測定された値の間では位相変化は２７０°であることが
わかる。全体として利用しうる複素数値の数は２倍にな
り、そして測定範囲もまた２倍になる。

【００４９】図３に示された複素測定値と複素計算値の
信号シーケンスが、コンピュータ装置２１４に送られ
る。コンピュータ装置２１４はフーリエ変換を実行す
る。

【００５０】測定された値ａ（ｋΔｔ）、ｂ（ｋΔｔ）
の間だけでなく、また複数個の計算された値ａ′（ｋΔ
ｔ／（ｎ＋１））、ｂ′（ｋΔｔ／ｎ＋１））の間でも
補間を行なうことができる。解析定理から、一義的で測
定可能な最大周波数ｆ_max は下記の式で与えられること
がわかっている。

【００５１】

【数８】 ｆ_max ＝±ｆ_s ／２ ， （６）

【００５２】ここで、ｆ_s は解析の周波数を示す。もし
ｎ個の値が測定された値の間で補間されるならば、その
場合には解析の周波数は（ｎ＋１）で乗算される。した
がって、新しい周波数限界ｆ_max は下記の式で表され
る。

【００５３】

【数９】 ｆ_max ＝±（ｎ＋１）ｆ_s ／２＝（ｎ＋１）ｆ_max （７）

【００５４】ｆ_max に対応する新しい上部速度限界ｖ
_max は方程式（３）からうることができ、その結果は下
記の通りである。

【００５５】

【数１０】 ｖ_max ＝（ｎ＋１）ｖ_max ｖ_max ＝（ｎ＋１）ｃ² ／８ｆ₁ ｄ cosθ （８）

【００５６】方程式（４）において、２周波数法の平均
速度の決定のために、上部限界ｖ_{a max} が存在する。ｖ
_max ＝ｖ_{a max} であることが好ましく、したがって、下
記の式がえられる。

【００５７】

【数１１】 （ｎ＋１）ｃ² ／８ｆ₁ ｄ cosθ＝ｃ² ／８（ｆ₂ −ｆ₁ ）ｄ cosθ （ｎ＋１）／ｆ₁ ＝１／（ｆ₂ −ｆ₁ ） ｎ＝（ｆ₁ ／（ｆ₂ −ｆ₁ ））−１ （９）

【００５８】方程式（９）は、測定された値の間で補間
されることが好ましい測定された値の総数を与える。

【００５９】従来の知識に従い、２周波数法で測定領域
を拡大することは、因子５だけ可能である。このこと
は、測定された２個の値毎の間で、４個の値までに補間
を行なうことである。

【００６０】第１実施例の説明 図４は、図１に示された本発明によるドップラ流速計の
第１実施例の図面である。図４において、送受信器２０
２に接続された超音波変換器２０１は送信パルスを送信
する。この送信パルスの周波数スペクトルは、２個の隣
接しているが分離している周波数帯域を有する。これら
の送信パルスとこれらの送信パルスを発生する装置は、
下記において、図７から図１６までの図面を参照して詳
細に説明される。図４による回路において、超音波変換
器２０１と送受信器２０２とからえられる受信信号は、
２個の直角位相チャンネルに送られる。これらのチャン
ネルのうちの１つのチャンネルは、平均周波数ｆ₂ を有
する周波数帯域の中の反射信号を処理する。第２直角位
相チャンネルは、平均周波数ｆ₁ を有する周波数帯域の
中の反射信号を処理する。これらの直角位相チャンネル
の構造は、パルス・ドップラ技術において、それ自身は
よく知られている。これらの直角位相チャンネルのおの
おのは帯域フイルタ２２１または２２２を有し、そして
その後に同期復調器を有する。この同期復調器は、乗算
器２２４，２２５または２２６，２２７を有する。これ
らの乗算器の出力信号は、低域フイルタ２３１，２３
２，２３３，２３４によってフイルタ作用を受け、そし
て積分保持回路２３５，２３６，２３７，２３８によっ
て解析される。さらに、直列に接続された一定符号消去
器２４１，２４２，２４３，２４４は、一定の反射信号
またはほぼ一定の反射信号を抑制する役割を果たす。こ
れらの一定符号消去器の出力信号Ｉ_{1 ,} Ｑ_{1 ,} Ｉ _{2 ,}
Ｑ₂ は、コンピュータ装置２０３に送られる。ここでま
ず、これらの信号がアナログからデイジタルに変換さ
れ、そして処理されることにより、コンピュータ装置２
０３の出力に平均速度ｖ_a に対応する出力信号がえられ
る。一定符号消去器２４３，２４４の出力信号Ｉ_{1 ,} Ｑ
₁ と、コンピュータ装置２０３の出力信号とが、コンピ
ュータ装置２１３の中で処理される。コンピュータ装置
の出力信号として、信号対Ｉ，Ｑのシーケンスが生ず
る。これらの信号対のシーケンスは、出力信号Ｉ₁ およ
びＱ₂ と、周波数ｆ₁ またはｆ₂ のうちの１つの周波数
に対し補間によってえられた信号対Ｉ′およびＱ′とで
構成される。

【００６１】コンピュータ装置２１３の出力の信号シー
ケンスは、コンピュータ装置２１４に送られる。コンピ
ュータ装置２１４の中で、よく知られているアルゴリズ
ムによる高速フーリエ変換により、スペクトルが計算さ
れる。コンピュータ装置２０３，２１３，２１４は、デ
イジタル信号処理装置、例えばアナログ・デバイスのＡ
Ｄ ２１００処理装置、で実現されることが好ましい。
コンピュータ装置２０３および２１３の中で使用される
アルゴリズムを下記で説明する。コンピュータ装置２１
３の出力信号は、デイジタル／アナログ変換器２１５，
２１６によって変換され、そしてそれらの出力が２個の
拡声器２１７，２１８に送られる。流速の拡大された測
定領域に対応するこの出力信号を用いて、音響的再生が
行なわれる。

【００６２】第２実施例の説明 図５は、図１に示された本発明によるドップラ流速計の
第２実施例の図面である。図５による回路において、超
音波変換器２０１および送受信器２０２は、図１および
図４で既に説明された。帯域フイルタ２５１を用いて、
反射信号は、送信信号の２つの周波数帯域の中で、すな
わち、信号処理のために透過することが許されている２
つの周波数帯域の中で、フイルタ作用を受ける。

【００６３】図５による回路によって、２重復調法が実
行される。帯域フイルタ２５１を透過した反射信号は、
cos ω₀t および−sin ω₀t を用いて、乗算器２５２
および２５３と、低域フイルタ２５４および２５５とに
より、まず復調される。ここで、ω₀ は下記の式で与え
られる。

【００６４】

【数１２】 ω₀ ＝（ω₁ ＋ω₂ ）／２ （１０）

【００６５】その後、cos Ωｔおよび−sin Ωｔを用い
て、乗算器２５６，２５７，２５８，２５９と低域フイ
ルタ２６１，２６２，２６３，２６４とにより、復調さ
れる。ここで、Ωは下記の式で与えられる。

【００６６】

【数１３】 Ω＝（ω₂ ＋ω₁ ）／２ （１１）

【００６７】したがって、要求された直角位相信号Ｉ
_{1 ,} Ｑ_{1 ,} Ｉ_{2 ,} Ｑ₂ が、加算器２６５，２６６，２６
７，２６８により適切な加算および減算を行って、これ
らの信号からえられる。これらの直角位相信号Ｉ_{1 ,} Ｑ
_{1 ,} Ｉ_{2 ,} Ｑ₂ をうるさいに、図５に示されているよう
に、積分および保持回路２７１，２７２，２７３，２７
４、および一定符号消去器２７５，２７６，２７７，２
７８が、図１および図４による回路と同じように用いら
れる。図５による回路は、図４の回路よりも複雑であ
る。けれども、図５の回路は、第１直角位相復調におい
て、送信信号の周波数ω₀ が用いられるという利点を有
する。この送信信号の周波数ω₀ は、高い安定度をもっ
て利用することができる。２つの直角位相チャンネルの
出力信号をうるための信号の処理を、下記で説明する。

【００６８】受け取られた信号のスペクトルは線スペク
トルである。単一線の間の間隔は、よく知られているよ
うに、パルスの繰り返し周波数である。ここで、下記の
ドップラ・シフトを有する１つの線ｓ(t) を考える。

【００６９】

【数１４】 s(t)＝cos （ω₀ ＋Ω＋ω_d2）ｔ＋cos ( ω₀ −Ω＋ω_d1）ｔ （１２）

【００７０】ここで、ω_d1およびω_d2はω₀ −Ωおよび
（ω₀ ＋Ω）におけるドップラ周波数である。

【００７１】この第１変調により、信号Ｉ(t) および信
号Ｑ(t) がえられる。これらの信号は下記のようにして
定めることができる。

【００７２】乗算器２５２の出力に、下記の信号がえら
れる。

【００７３】

【数１５】cos ω₀t cos （ω₀ ＋Ω＋ω_d2）ｔ ＋cos ω₀t cos （ω₀ −Ω＋ω_d1）ｔ.

【００７４】低域フイルタ２５４の出力に、周波数差を
有する要素のみが存在する。したがって、下記の式がえ
られる。

【００７５】

【数１６】 I(t)＝（１／２）cos ( Ω＋ω_d2)t＋（１／２）cos ( ω_d1−Ω)t （１３）

【００７６】乗算器２５３の出力における直角位相チャ
ンネルに対し、下記の信号がえられる。

【００７７】

【数１７】−sin ω₀t cos( ω₀ ＋Ω＋ω_d2）t −sin ω₀t cos( ω₀ −Ω＋ω_d1）t

【００７８】低域フイルタ２５５の出力に、周波数差を
有する要素のみが存在する。したがって、下記の式がえ
られる。

【００７９】

【数１８】 Ｑ(t) ＝（１／２）sin ( Ω＋ω_d2)t−（１／２）sin ( Ω−ω_d1)t (14)

【００８０】第２変調により、復調器の出力に下記の信
号がその度にえられる。乗算器２５６の出力に下記の信
号がえられる。

【００８１】

【数１９】 〔（１／２）cos ( Ω＋ω_d2)t＋（１／２）cos ( Ω−ω_d1)t〕ｃosΩｔ

【００８２】低域フイルタ２６１の出力に、下記の信号
がえられる

【００８３】

【数２０】 ａ(t) ＝（１／４）cos ω_d2ｔ＋（１／４）cos ω_d1ｔ （１５）

【００８４】乗算器２５７の出力に、下記の信号がえら
れる。

【００８５】

【数２１】 〔(1／２）cos ( Ω＋ω_d2)t＋(1／２）cos ( Ω−ω_d1)t〕( −sin Ωｔ）．

【００８６】低域フイルタ２６２の出力に、下記の信号
がえられる。

【００８７】

【数２２】 ｂ(t) ＝（１／４）sin ω_d2t −（１／４）sin ω_d1ｔ （１６）

【００８８】乗算器２５８の出力に、下記の出力がえら
れる。

【００８９】

【数２３】 〔(1／2) sin( Ω＋ω_d2)t−(1／２）sin ( Ω−ω_d1)t〕( −sin Ωｔ）

【００９０】低域フイルタ２６３の出力に、下記の信号
がえられる。

【００９１】

【数２４】 c(t)＝−（１／４）cos ω_d2t ＋（１／４）cos ω_d1ｔ （１７）

【００９２】乗算器２５９の出力に、下記の信号がえら
れる。

【００９３】

【数２５】 〔(1／2) sin( Ω＋ω_d2)t− 1／２ sin (Ω−ω_d1) t 〕 cosΩｔ

【００９４】低域フイルタ２６４の出力に、下記の信号
がえられる。

【００９５】

【数２６】 d(t)＝（１／４）sin ω_d2t ＋（１／４）sin ω_d1ｔ （１８）

【００９６】したがって、加算器２６５，２６８と、直
列接続された回路とを用いて、図５に示された下記で定
められる出力信号Ｉ_{2 ,} Ｑ_{2 ,} Ｉ_{1 ,} Ｑ₁ がえられる。

【００９７】

【数２７】 Ｉ₂(t)＝ａ(t) −ｃ(t) ＝（１／２）cos ω_d2t （１９） Ｑ₂(t)＝ｂ(t) ＋ｄ(t) ＝（１／２）sin ω_d2t （２０） Ｉ₁(t)＝ａ(t) ＋ｃ(t) ＝（１／２）cos ω_d1ｔ （２１） Ｑ₁(t)＝ｄ(t) −ｂ(t) ＝（１／２）sin ω_d1ｔ （２２）

【００９８】これらの方程式から、図５の回路により、
要求された直角位相信号をうることができることがわか
る。ここで注目しなければならないことは、ω_d1は送信
周波数ω₁ ＝ω₀ −Ωに付随した周波数であり、そして
ω_d2は送信周波数ω₂ ＝ω₀＋Ωの周波数ω₂ に付随す
る周波数であることである。したがって、下記のドップ
ラ効果の公式により、対応する速度の値を計算すること
ができる。

【００９９】

【数２８】 ｖ₁ ＝−ω_d1 ^.ｃ／２ ^.（ω₀ −Ω）cos θ （２３） ｖ₂ ＝−ω_d2 ^.ｃ／２ ^.（ω₀ ＋Ω）cos θ （２４）

【０１００】信号Ｉ₁ ，Ｑ₁ ，Ｉ₂ ，Ｑ₂ をさらに処理
するための回路は、図４のところで既に説明した。した
がって、これらの回路は図５には示されていない。

【０１０１】本発明のアウトラインにおいて使用される
アルゴリズムの説明 図４のコンピュータ装置２０３において、平均速度が２
周波数法により計算される。コンピュータ装置２１３に
おいてさらに、ドップラ路の信号の値は、図２および図
３のところで説明したように補間される。これらの計算
のためのアルゴリズムをここで説明しなくてはならな
い。コンピュータ装置２１４は高速フーリエ変換を計算
する。この高速フーリエ変換の計算のアルゴリズムはよ
く知られており、したがって、ここでは説明しない。

【０１０２】２周波数法による平均速度の計算のための
アルゴリズム 方程式（１９）〜（２２）のドップラ信号は、複素信号
Ｐ（ｔ）の形で下記のように表すことができる。

【０１０３】

【数２９】 P₁(t) ＝I₁(t) ＋jQ₁(t)＝(1/2)(cos ω_d1ｔ＋jsinω_d1ｔ） ＝(1/2)exp(jω_d1ｔ） （２５） P₂(t) ＝I₂(t) ＋jQ₂(t)＝(1/2)(cos ω_d2ｔ＋jsinω_d2ｔ） ＝(1/2)exp(jω_d2ｔ） （２６）

【０１０４】方程式（２５）および（２６）は、２つの
複素指標を定める。この複素指標は、角ドップラ速度ω
_d1またはω_d2で回転する。これらの角ドップラ速度のう
ちの１つ、例えばω_d1、を決定するために、引き続く２
つの時刻と一致してこの指標の位相差を測定し、そして
この位相差から角ドップラ速度の値が下記の式で計算さ
れる。

【０１０５】

【数３０】 ω_d1＝φ₁ ／Δｔ （２７）

【０１０６】ここで、φ₁ はＰ₁ （ｔ＋Δｔ）とＰ
₁ （ｔ）との間の位相差である。この位相差は複素乗数
因子によって定めることができる。

【０１０７】

【数３１】 φ₁ ＝arg 〔Ｐ₁ （ｔ＋Δｔ）Ｐ^* ₁ （ｔ）〕 （２８）

【０１０８】ここで、arg 〔 〕は複素数の偏角
を意味し、そしてＰ₁ ^* （ｔ）はＰ ₁ （ｔ）の複素共役
値を意味する。

【０１０９】測定されるべき速度に対応する信号の中の
雑音と振動のために、方程式（２８）によってのφ₁ の
計算値は必らずしも信頼できるものではない。この値の
信頼性は平均をとることによって改善することができ
る。下記のベクトルを計算することが有効であることが
わかっている。

【０１１０】

【数３２】

【０１１１】ここで、ΣはＮ個の値についての和を意味
する。方程式（２９）から、Ｐ_1Nおよびしたがってφ₁
を下記のように決定することができる。

【０１１２】

【数３３】 φ₁ ＝arg （Ｐ_1N） （３０）

【０１１３】指標Ｐ_1Nは規格化される。すなわち、｜Ｐ
_1N｜は、ゼロに等しいかまたはゼロより大きく、および
１に等しいかまたは１より小さい範囲内にある。すなわ
ち

【０１１４】

【数３４】 ０≦｜Ｐ_1N｜≦１ （３１）

【０１１５】指標Ｐ_1Nの長さは信号の品質を表す。Ｐ_1N
の計算で考察された指標は図６に示されている。もし加
算されるすべての指標が同じ方向にあるならば、｜Ｐ_1N
｜＜＜１である。

【０１１６】同様の方式で、第２周波数に対して下記の
計算を行なうことができる。

【０１１７】

【数３５】

【０１１８】ドップラ方程式（２）により、そしてもし
この方程式でcos θ＝１（cos θ＝１は、半径の軸上に
投影された速度をうることを意味する。）であるなら
ば、下記の式をうることができる。

【０１１９】

【数３６】 φ₁ ＝ω_d1・Δｔ＝−ω₁ Δｔ ２ｖ／ｃ （３４） φ₂ ＝ω_d2・Δｔ＝−ω₂ Δｔ ２ｖ／ｃ （３５）

【０１２０】ここで、ω₁ ＝２πｆ₁ およびω₂ ＝２π
ｆ₂ である。

【０１２１】方程式（３５）から方程式（３４）を減算
することにより、下記の式をうる。

【０１２２】

【数３７】 Δφ＝φ₂ −φ₁ ＝−Δｔ（ω₂ −ω₁ ）２ｖ_a ／ｃ （３６）

【０１２３】このことから、２周波数法で決定すること
ができる速度は、下記の式でうることができる。

【０１２４】

【数３８】 ｖ_a ＝−Δφ ^. c／〔２Δｔ（ω₂ −ω₁ ）〕 （３７）

【０１２５】Δφの値は、−πに等しいかまたは−πよ
り大きいから、πに等しいまたはπより小さい領域内に
おいて、一義的に定まる。この領域はナイキストの領域
に対応する。Δφはφ₁ およびφ₂ よりはるかに小さい
から、ｖ_a の一義的測定領域は拡大される。

【０１２６】経験によれば、方程式（２９）および（３
２）において、もし４個の値ないし８個の平均が実行さ
れるならば、すなわち、Ｎが例えば４ないし８の値をも
つことができるならば、良好な結果のえられることがわ
かっている。

【０１２７】ドップラ信号の補間のためのアルゴリズム 補間は、ｖの一義的測定領域を拡大するために、ドップ
ラ信号を定められた周波数に変更する役割りを果たす。
図２および図３のところで説明したように、これらの測
定値の間でまた別の値が補間されることにより、このこ
とがえられる。１つの例として、周波数ω₁ における信
号Ｐ₁ （ｔ）の２つの測定された値の間での値の補間を
考える。測定された信号Ｐ₁ （ｔ）は下記の値から成
る。

【０１２８】

【数３９】 Ｐ₁ （０），Ｐ₁ （Δｔ），Ｐ₁ （２Δｔ），Ｐ₁ （３Δｔ），・・・(38)

【０１２９】補間により、Δｔ／２、３Δｔ／２、５Δ
ｔ／２、などのまた別の値を計算することができる。下
記において、Ｐ₁ （Δｔ／２）の計算を説明する。以前
には、Ｐ₁ （△ｔ）とＰ₁ （０）との間の位相差φ₁ は
方程式（２８）により計算される。この位相φ₁ は−π
から＋πの範囲内において一義的である。すなわち、下
記の式が成り立つ。

【０１３０】

【数４０】 −π≦φ₁ ≦＋π （３９）

【０１３１】補間のためのこの他の情報として、方程式
（３７）からわかるように、２周波数法によって計算さ
れた速度ｖの平均値ｖ_a が用いられる。ｖ_a に基づい
て、方程式（３４）により位相φ₁ を計算することがで
きる。

【０１３２】

【数４１】 φ₁ ′＝−ω₁ （２Ｖ_a ／ｃ）Δｔ （４０）

【０１３３】値φ₁ ′を用いて、値φ₁ が−πから＋π
までの領域内に実効的にとどまっているかどうかがわか
る。｜φ₁ ′−φ₁ ｜＞πである時、Ｐ₁ （Δｔ）とＰ
₁ （０）との間の位相差φ₁ に２πの倍数ｋが加えられ
る。したがって、ｋの値は下記の式によって決定され
る。

【０１３４】

【数４２】 ｋ＝round （〔φ₁ ′−φ₁ 〕／２π） （４１）

【０１３５】ここで、round （ ）は次の整数
に丸めることを意味する。したがって、値Ｐ₁ （Δｔ／
２）は下記のように補間される。

【０１３６】

【数４３】 arg Ｐ₁ （Δｔ／２） ＝arg Ｐ₁ （０）＋（１／２）（φ₁ ＋k ^. 2 π) （４２） ｜Ｐ₁ （Δｔ／２）｜ ＝（｜Ｐ₁ （０）｜＋｜Ｐ₁ （Δｔ）｜）／２ （４３）

【０１３７】したがって、信号Ｐ₁ （Δｔ／２）の偏角
と振幅が、周知の方法で、直角位相信号Ｉ₁ （Δｔ／
２）およびＱ₁ （Δｔ／２）に再び変換される。

【０１３８】もしＰ₁ （０）とＰ₁ （Δｔ）の間で（ｍ
−１）個の中間値を補間することが必要であるならば、
その時には間隔ｔはｍ個の長さに分割される。したがっ
て、ｎ番目の補間値は下記の式によってえられる。

【０１３９】

【数４４】

【０１４０】測定された値またはこれらの測定された値
の間で補間された値に対応する信号から成る信号のシー
ケンスが、コンピュータ装置２１４に送られる。コンピ
ュータ装置２１４では、高速フーリエ変換が計算され
る。

【０１４１】全体的な所見として、平均速度の値が比較
的少数個の値、典型的には４個ないし８個の値、から決
定されることが、再び繰り返されなければならない。こ
れとは違って、フーリエ変換は大幅に多数個の値から計
算される。したがって、平均周波数およびしたがって速
度の平均値が、フーリエ変換が計算される信号シーケン
スに沿って、変動することは確かにありうる。このこと
は、フーリエ変換の周波数領域がもとのナイキスト領域
より大きいことを意味する。これは前記トラッキング法
に対するまた別の相異点である。これらのことは、信号
帯域の幅がもとのナイキスト領域より小さい時にのみあ
てはまる。

【０１４２】送信パルスの生成装置の説明 図７は、超音波を反射する粒子１２を運んでいる流体１
１の流速を測定するドップラ流速計のブロック線図であ
る。この流体は、例えば、血管の中を流れている血液で
あることができる。この装置の原理は、ドイツ国特許出
願番号第ＤＥ−Ａ−２ ４０６６３０号、および対応す
る米国特許第ＵＳ−Ａ−３ ９１４９９０号に開示され
ている。このよく知られている原理により、例えば管１
３の中を流れている流体１１が、変換器１４から送信さ
れた少なくとも２個の引き続く超音波パルスで照射され
る。送信された超音波ビーム内にある流体中の反射体
（例えば粒子）によって対応するドップラ周波数だけ変
化した反射波が、同じ変換器１４によって受信される。
第１送信超音波パルスと第２送信超音波パルスと同じ通
過時間を有する反射波の間の位相差から、評価装置１８
の中で適切な処理を行なうことにより、その時間的ペー
スが検査される横断面内の流体の速度分布に対応する出
力信号がえられる。

【０１４３】この装置は、超音波変換器１４と、送信器
１５と、受信器１６と、中央制御装置１７と、評価装置
１８とを有する。

【０１４４】点線７で示された方向に、超音波変換器１
４が超音波パルスを放射して流体１１を照射する。これ
らの送信パルスは予め定められた繰り返し周波数で送り
出される。変換器１４は、また流体内の粒子によって反
射された反射波を受信し、そして対応する反射波信号を
送り出す。

【０１４５】変換器１４は例えば下記の技術データを有
する超音波変換器である。送信周波数は約３．０ＭＨｚ
でなければならない。 −セラミックの共振周波数，３．２ＭＨｚ −円形セラミックの直径，１３ｍｍ −セラミックの曲率半径，１２０ｍｍ −変換器と近視野および遠視野の限界との間の距離，８
８ｍｍ。

【０１４６】変換器１５は超音波変換器１４に接続さ
れ、そして送信パルスを生成する機能を有する。この送
信パルスによって超音波変換器が励振され、そして超音
波パルスが放射される。

【０１４７】受信器１６が変換器１４に同じように接続
され、そしてそこで反射信号の処理が行なわれる。これ
らの反射信号は、放射された超音波の第１パルスと第２
パルスに応答して流体内の粒子によって反射された、少
なくとも２群の反射信号に対応する。反射信号の隣接し
ているがしかし分離している周波数帯域が、受信器１６
の中で、分離した信号処理チャンネル２１または２２の
おのおのの中で処理される。

【０１４８】受信器１６の出力に接続された評価装置１
８により、流速に関する情報を有する少なくとも１つの
出力信号が、受信器でえられたドップラ情報からえられ
る。

【０１４９】中央制御装置１７は、このドップラ計のプ
ログラム可能な動作のために必要である、すべての制御
信号およびタイミング信号を生ずる。したがって、制御
装置１７は実質的にプログラム可能デイジタル信号発生
器である。中央制御装置１７の最も重要な機能は下記の
通りである。中央制御装置１７は送信器１５を制御し
て、送信パルスの時刻と持続時間を決定する。中央処理
装置１７は受信器１６に対してすべてのタイミング信号
を供給し、したがって、受信器の中での信号の流れを制
御する。中央処理装置１７は評価装置１８に対しタイミ
ング信号を供給し、したがって、評価装置１８の中の信
号の流れを制御する。

【０１５０】本発明により、送信パルスが周期的繰り返
されたシーケンスの周波数スペクトルが、図８に示され
た、２つの隣接しているがしかし分離した周波数帯域３
１および３２から成るように、送信器１５が設計され
る。送信パルスはパルスの繰り返し周波数で周期的に放
射されるから、送信信号のスペクトルは周波数線スペク
トルから成り、これらのスペクトル線の間の間隔はパル
ス繰り返し周波数に対応する。送信信号の周波数スペク
トルの形または構造により、送られるエネルギの大部分
は２つの周波数ｆ₁ またはｆ₂ の付近の、これらの周波
数帯域の中に集中するという結果になる。

【０１５１】周波数帯域３１または３２は相互に非常に
接近していて、これらの周波数帯域を同じ超音波変換器
で放射することができ、そしてこれらの帯域が同じ減衰
を受けることができる。周波数ｆ₀ ，ｆ₁ およびｆ
₂ は、例えば次の値、すなわち、ｆ₀ ＝３．２ＭＨｚ，
ｆ₁ ＝２．９２ＭＨｚおよびｆ₂ ＝３．４５ＭＨｚの値
を有する。

【０１５２】放射されるエネルギの大部分は２つの隣接
する周波数帯域３３または３４の中に集中することが好
ましい。搬送波周波数ｆ₀ はこれら２つの帯域の限界を
定め、そしてｆ₁ は周波数帯域３３の平均周波数であ
り、およびｆ₂ は周波数帯域３４の平均周波数である。

【０１５３】周波数ｆ₁ と周波数ｆ₂ との間の間隔は、
周波数間隔ｆ₀ として定められる。ｆ₀ に対してｆ₁ と
ｆ₂ は対称的に配置されているから、方程式ｆ₂ −ｆ₀
＝ｆ ₀ −ｆ₁ ＝Δｆ／２が成り立つ。

【０１５４】図９は、引き続く２個の送信パルス３５ま
たは３６の１つの例を示した図面である。これらの送信
パルスによって放射されるシーケンスは、図８による送
信スペクトルを有する。これらの送信パルス３５または
３６のおのおのは、主として周波数ｆ₀ の方形波信号で
構成され、この方形波信号の振動周期はｔ₀ ＝１／ｆ ₀
である。１つの送信パルスの持続時間はＴ_p であり、例
えばＴ_p ＝６マイクロ秒である。引き続く送信パルスの
間の間隔Ｔ_r は例えば２５０マイクロ秒である。この間
隔Ｔ_r はパルス繰り返し周波数ＰＲＦ＝４ＫＨｚに対応
する。

【０１５５】変換器に送られる送信パルスは、例えば、
５０ボルトと１００ボルトの間の電圧を有する。

【０１５６】図９に示されているように、送信パルスの
おのおのは周期的方形波信号からずれている。それは、
信号の中央の点３７のところで、位相が反転しているか
らである。この位相の反転は、送信パルスの周期的シー
ケンスが図８による周波数スペクトルを有するために必
要である。

【０１５７】図１５（ａ）および（ｂ）は送信パルスの
また別の２つの例を示している。図１５（ａ）による送
信パルスは、３つの要素送信パルスから成っている。こ
れらの要素送信パルスは、図９による送信パルスのペー
スを有する。図１５（ｂ）による送信パルスは、図９に
よる送信パルスのペースを有する２個の完全な送信パル
ス要素と、始めの部分と終りの部分とに１つのこのよう
な送信パルス要素の一部分とを有する。

【０１５８】図１５（ｂ）による送信パルスは完全な送
信パルスを３個以下有し、そしてその持続時間は図１５
による送信パルスの持続時間よりわずかに短い。図１５
（ｂ）に示されているように、送信パルスは送信パルス
要素の要求された端部と一致して開始することができ
る。

【０１５９】図１５（ａ）および（ｂ）による送信パル
スは、位相が反転する点１８１〜１８５を有する。ここ
で、位相反転点１８１、１８３、１８５のおのおのは、
図９の位相反転点３７に対応し、一方位相反転点１８
２、１８４は、信号パルスの２つの要素を接合すること
によって作られる。実際的な実施例では、送信パルスは
それ自身連続である信号の時間的区分わけによって作ら
れる。この送信パルスは、周期的に分布した位相反転点
を有する。

【０１６０】図１６は、図１５（ａ）または（ｂ）に示
された波形を有する送信パルスのシーケンスの周波数ス
ペクトルの概要図である。図面を簡単にするために、図
１６では、周波数スペクトル線の外郭線だけが示されて
いるが、実際は、図８のような周波数スペクトルを有す
る。そのおのおのが１つの送信パルス要素のみを有する
送信パルス・シーケンスのための図８による送信スペク
トルと、そのおのおのが種々の送信パルス要素を有する
送信パルス・シーケンスに対する図１６による周波数ス
ペクトルとを比較すると、図１５（ａ）および（ｂ）に
よるこれらの前記送信パルスを用いることにより、図８
のように間隔ｆが等しい周波数ｆ₁ および周波数ｆ₂ の
付近のより狭い周波数帯域の中に、放射される超音波エ
ネルギがより強く集中することがわかる。

【０１６１】図８におけるように、図１６において、ｆ
₁ とｆ₂ との間の間隔Δｆが示されている。

【０１６２】図１５（ａ）に示されているように、位相
反転点の周期は、すなわち、等しい位相の反転点の間の
時間間隔はＴ_u で示されている。ここで、Ｔ_u は２／Δ
ｆに等しい。このことはまた、図１５（ｂ）による送信
パルスに対してもいえる。

【０１６３】図１５（ａ）または（ｂ）による送信パル
スの中で、位相反転点は規則正しい時間間隔Ｔ_e ＝Ｔ_u
／２になければならない。

【０１６４】もし送信信号が持続時間Ｔ_p の送信パルス
のシーケンスで構成され、そして図１５（ａ）および
（ｂ）のようにＴ_p はＴ_u より非常に大きいならば、送
信パルスの開始の位置とその終了の位置は、その中に含
まれる位相の反転点の位置に関して必須のものではな
い。

【０１６５】送信パルスの持続時間Ｔ_p は位相反転点の
周期Ｔ_u より大きいことが好ましい。図８による送信信
号スペクトルを確実にうるために、Ｔ_p はＴ_u より大き
く選定されなければならない。Ｔ_u に対してＴ_p が大き
くなればなる程、送信信号の周波数スペクトルの中の平
均周波数ｆ₁ と平均周波数ｆ₂ の付近の帯域の中に、放
射されるエネルギがますます強く集中する。

【０１６６】図１０は、図７の送信器１５の第１実施例
を示す。この実施例は、実質上、デイジタル部品で構成
される。図１０による送信回路は、制御インタフェース
４１と、制御装置４２と、２個のスクロール・レジスタ
と、電力増幅器４５とを有する。電力増幅器４５の出力
は、図７の超音波変換器１４に接続される。図１０によ
る回路により、送信パルスが図９に従って発生される。
そこでは、送信パルスの送信された周期的シーケンス
は、図８に従う周波数スペクトルを有する。

【０１６７】必要な信号は、線路３９、５６、５１、ア
ドレス・バス４６およびデータ・バス４７によって、制
御インタフェース４１に送られる。インタフェース４１
は、線路４８または４９によって、レジスタ４３または
４４に接続される。レジスタ４３または４４の出力は、
おのおのが電力増幅器４５の入力に接続される。制御装
置４２は、線路５１、５２、５３によって、図７の中央
制御装置１７に接続される。制御装置４２はまた、線路
５６によってインタフェース４１に接続され、および線
路５４または５５によってレジスタ４３または４４に接
続される。

【０１６８】レジスタ４３およびレジスタ４４は同じ構
造を有し、かつ、対応した接続を有する。けれども、こ
れらのレジスタは異なるがしかし相補的であるバイトで
ロードされる。

【０１６９】図１１は、図１０のレジスタ４３の概要図
である。図１１に示されているように、このレジスタに
おいて、４個のバイト（ワード）６５、６６、６７、６
８が、バス６１、６２、６３、６４の対応する線路によ
って、レジスタ４３の中にロードされる。これらのバイ
トの内容が図１１のレジスタ４３のメモリ・セルの中に
示されている。

【０１７０】図１０からわかるように、レジスタ４３お
よびレジスタ４４のおのおのは、その出力からその入力
へのフィードバックを有する。下記で詳細に説明される
ように、図１０による送信回路が適切に動作するために
は、このことは重要である。レジスタ４３とレジスタ４
４の内容が相補的である、すなわち、例えばレジスタ４
３の中に１０１０１０‥が記憶されているならば、レジ
スタ４４の中には０１０１０１‥が記憶されることが、
同様に重要である。

【０１７１】図１０および図１１に従う送信回路の動作
を図１２の信号図を参照して説明する。

【０１７２】制御インタフェース４１は、制御装置４２
とスクロール・レジスタ４３および４４をプログラムす
る役割りを果たす。

【０１７３】この目的のために、インタフェースは、ア
ドレス・バスまたはデータ・バス４６、４７から、およ
び線路３９から、図７による原理装置の復号化マイクロ
プロセッサのバスからの図１２の付随するストロボ信号
７２を受け取る。信号７２は、レジスタ４３、４４に送
られるバイトのロードを実行するタイミング信号であ
る。レジスタ４３またはレジスタ４４に供給されるバイ
トのロードを制御するために、インタフェース４１は、
線路５１および線路５６からまた別の信号７１を受け取
る。供給されたバイトをロードすることは、信号７１が
図に示された状態にある時間間隔１７１においてのみ可
能である。レジスタ４３またはレジスタ４４に供給され
たバイトのロードを制御するために、信号７１に応答し
て制御装置４２の中に、また別の信号７３が発生され
る。この信号７３は、線路５４または線路５５により、
これらのレジスタに送られる。線路５２により、制御装
置４２はまた、図７の中央制御装置１７から、タイミン
グ信号を受け取る。このタイミング信号に応答して、制
御装置４２は対応するタイミング信号７４を発生する。
このタイミング信号７４は、線路５４または線路５５に
よって、レジスタ４３またはレジスタ４４に供給され
る。

【０１７４】図１２に示されているように、レジスタ４
３またはレジスタ４４のロードのためにえられた時間間
隔１７１は、時刻１７２と時刻１７３の間持続する。こ
の時間間隔において、信号７１および信号７３の状態の
変動がその度に検査される。タイミング信号７４のより
近い前端を有する時間間隔１７１の後、時刻１７４にお
いて、スクロール・サイクルがレジスタ４３またはレジ
スタ４４の中で始まる。タイミング信号７４によって実
行されるこのスクロール・サイクルにより、出力信号７
６または出力信号７７がレジスタ４３またはレジスタ４
４の出力のおのおのに生ずる。図１０の電力増幅器４５
により、出力信号７６と出力信号７７との差がえられ
る。この方式で送信パルス７８がえられる。この送信パ
ルスは図９による要求されたペースを有する。

【０１７５】図１２に示されているように、レジスタ４
３またはレジスタ４４の中のスクロール・サイクルは、
時間間隔１７６に拡大される。時間間隔１７６は、時刻
１７４から時刻１７５まで持続する。時間間隔１７６の
持続時間は、したがって送信パルス７８の持続時間は、
制御装置４２の中のカウンタによって決定される。制御
装置４２は、時刻１７５に、信号の変化を生じ、それが
送信パルス７８の終了を生ずる。制御装置４２の線路５
３により、信号７５が図７の中央制御装置１７に送信さ
れる。

【０１７６】レジスタ４３またはレジスタ４４は、レジ
スタ４３に対して図１１に示された構造を有する、スク
ロール・レジスタである。図１１に示されているよう
に、これらのレジスタのおのおのは４×８メモリ・セル
を有し、そのおのおのが１バイト用である。これらは、
おのおのが８ビットの４バイトでもって、インタフェー
ス４１によりロードされる。前記で説明したように、そ
の後、３２ビット長のサンプル（送信サンプル）が電力
増幅器４５の１つの入力へ移動される。レジスタの入力
へ出力をフィードバックすることにより、要求される長
さの間、スクロール・サイクルを実行することができ
る。

【０１７７】図１０に従っていま説明した送信回路は好
ましいものである。それは、比較的経済的なデイジタル
部品で主として構成されているからであり、そして操作
パラメータを比較的簡単に変えることにより、送信パル
スの周波数スペクトルを柔軟に構成することができるか
らである。けれども、送信パルスを発生するために、図
１０による送信回路の代わりに、図１３による送信回路
を用いることができる。図１３の送信回路は、実質的
に、振幅変調を行なう変調器と、電子スイッチと、電力
増幅器とから成る。

【０１７８】図１３による送信回路は、次の機能群を有
する。すなわち、周波数ｆ₀ の搬送波信号を発生する搬
送波周波数発振器８１と、周波数ｆ_m ＝Δｆ／２を有す
る変調信号を発生する変調周波数発振器８２を有する。
この実施例では、ｆ_m ＝２５０ＫＨｚである。図１０に
よる送信器は、乗算器８３と、案内可能な電子ゲート回
路８４と、電力増幅器８５とを有する。電力増幅器８５
の出力は、図７の超音波変換器に接続される。乗算器８
３により、発振器８１から送られてくる搬送波信号と発
振器８２から送られてくる変調信号とが乗算され、その
結果図１４に示された出力信号がえられる。図１４から
わかるように、この方法で作られた振幅変調信号は、従
来の振幅変調の結果である。この振幅変調信号が、線路
８６によって、ゲート回路８４の入力に送られる。図７
の中央制御装置１７から送られてきたパルス８８が、線
路８７入力に、ゲート回路の第２入力に供給され、その
結果、振幅変調信号が設定された時間間隔にゲート回路
８４を通過する。このようにして、送信パルスが作られ
る。この送信パルスが電力増幅器８５で増幅された後、
超音波変換器１４に送られる。図１３のパルス８８の持
続時間は、送信パルスの持続時間を決定する。

【０１７９】図１４からわかるように、ゲート回路８４
の入力における振幅変調された波形は、周期的に位相が
反転する点８９を有する。図１３のパルス８８と発振器
８１または発振器８２から供給される信号との間の適切
な時間関係により、ゲート回路８４の出力に生ずる送信
パルスは、１つのこのような位相の反転点８９を有す
る。このように（図９による送信パルスの場合のよう
に）、超音波変換器に供給される送信パルスは、図８に
従う周波数スペクトルを確実に有する。

【図面の簡単な説明】

【図１】本発明によるドップラ流速計の原理的構造を表
す非常に単純化されたブロック線図。

【図２】図１のコンピュータ装置２１３を示す図。

【図３】図１のコンピュータ装置２１３を示す図。

【図４】図１による本発明のドップラ流速計の第１実施
例の図。

【図５】図１による本発明のドップラ流速計の第２実施
例の図。

【図６】指標のグラフを示す図。

【図７】ドップラ流速計のブロック線図。

【図８】図９に従う波形を有する送信パルスのシーケン
スの周波数スペクトル図。

【図９】送信パルス図。

【図１０】図７の送信器１５の第１実施例の１つのブロ
ック線図。

【図１１】図１０のレジスタ４３のブロック線図。

【図１２】図１０に従う送信動作を示した信号図。

【図１３】図７の送信器１５の第２実施例のブロック線
図。

【図１４】図１３のゲート回路８４の入力信号の図。

【図１５】送信パルス図。

【図１６】図１５（ａ）または（ｂ）による波形を有す
る送信パルスのシーケンスの周波数スペクトルの概要
図。

【符号の説明】

１４ 超音波変換器 １５ 送信器 １６ 受信器 １８ 評価装置

Claims (10)

【特許請求の範囲】
1. 【請求項１】 （イ） 予め定められたパルス繰り返し
   周波数で供給される送信パルスに対応して応答する超音
   波パルスで流体を照射し、かつ、流体内の粒子によって
   反射された反射波を受け取り、かつ、対応する反射信号
   を発生する、超音波変換器と、 （ロ） 送信パルスを発生し、かつ、前記送信パルスに
   よって前記超音波変換器が励振されて前記超音波パルス
   が放射される、超音波変換器に接続された送信器と、 （ハ） 放射された第１超音波パルスおよび第２超音波
   パルスに応答して流体路内の設定点の中にある粒子によ
   って反射された少なくとも２つの異なる反射波に対応す
   る反射信号を受け取りかつ処理するために超音波変換器
   に接続され、かつ、反射信号の隣接しているがしかし分
   離した周波数帯域が別々の信号処理チャンネル内でそれ
   ぞれ処理される、受信器と、 （ニ） 前記受信器の出力に接続され、かつ、前記受信
   器でえられたドップラ情報に基づく流速に対応する信号
   出力がえられる、評価装置と、を有し、および （ホ） 前記受信器が（ホ１） 前記信号処理チャンネ
   ルのそれぞれにそなえられ、かつ、前記各周波数帯域中
   で前記超音波変換器から受け取った前記反射信号の直角
   位相復調を実行し、各処理チャンネルの出力に、測定さ
   れた複素値（Ｐ（０）、Ｐ（Δｔ）、Ｐ（２Δｔ），
   …）を定める各直角位相信号の対をその都度うることが
   できる装置と、（ホ２） ２つの前記信号チャンネルか
   らの直角位相信号を処理し、かつ、流体路の中の解析さ
   れる点における流速の平均値に対応する第１出力信号を
   生ずることができる装置と、（ホ３） 第１出力信号と
   前記信号処理チャンネルのうちの１つのチャンネルの出
   力の直角位相信号とを論理的に相関させ、かつ、この論
   理的相関により、流速によって決定された平均値および
   種々の測定された測定値と一致する計算された複素値
   （Ｐ（Δｔ／２）、Ｐ（３Δｔ／２）…）を定める信号
   を生ずることができる装置と、（ホ４） 測定された複
   素値に対応する信号と計算された複素値に対応する信号
   （Ｐ（０）、Ｐ（Δｔ／２）、Ｐ（Δｔ）、Ｐ（３Δｔ
   ／２）、Ｐ（２Δｔ）、…）から成る信号シーケンスを
   評価し、かつ、この評価によって流速の瞬間値に関する
   情報がえられる装置と、を有することを特徴とする、超
   音波を反射する粒子を含む流体の流速を測定するための
   ドップラ流速計。
2. 【請求項２】 請求項１において、前記送信パルスの周
   期的シーケンスの周波数スペクトルが２つの隣接してい
   るがしかし分離した周波数帯域から成ることを実現する
   ように前記送信器が作成されることを特徴とする、ドッ
   プラ流速計。
3. 【請求項３】 請求項１において、測定された複素値に
   対応する信号および計算された複素値に対応する信号
   （Ｐ（０）、Ｐ（Δｔ／２）、Ｐ（Δｔ）、Ｐ（３Δｔ
   ／２）、Ｐ（２Δｔ）、…）から成る信号シーケンスの
   解析を実行し、かつ、この解析によりその振幅が信号シ
   ーケンスの中の信号によって定められる指標の位相変化
   に対応し、したがって流速の瞬間値に対応する第２出力
   信号を生ずる装置を前記送信器が有することを特徴とす
   るドップラ流速計。
4. 【請求項４】 請求項１において、測定値と計算値に
   対応する信号から成る信号シーケンスの中の信号はこれ
   らの信号の間の時間間隔が測定された値に対応する信号
   に対するものよりも短いことを特徴とするドップラ流速
   計。
5. 【請求項５】 請求項１において、前記指標の位相変化
   の解析を実行する装置がフーリエ変換を実行する装置を
   有することを特徴とするドップラ流速計。
6. 【請求項６】 請求項１において、前記第１出力信号を
   生ずる装置が４個ないし８個の反射信号の中に含まれて
   いる情報に基づいて流速の平均値を計算する装置を有す
   ることを特徴とするドップラ流速計。
7. 【請求項７】 請求項５において、フーリエ変換の実行
   のさいに流速の平均値の計算のためにその情報が考慮さ
   れる反射信号の数よりは大きな一定数の反射信号の中に
   含まれている情報が考慮されることを特徴とするドップ
   ラ流速計。
8. 【請求項８】 請求項１において、信号シーケンスの音
   響的再生のための装置を有することを特徴とするドップ
   ラ流速計。
9. 【請求項９】 請求項１において、２つの続いて起こる
   値の間で１つの値を計算することができかつ前記計算値
   を測定された値の補間によって計算することができる装
   置を有することを特徴とするドップラ流速計。
10. 【請求項１０】 請求項１において、２つの続いて起こ
    る測定値の間で種々の計算値をうることができかつ前記
    計算値を測定値の補間によって計算することができる装
    置を有することを特徴とするドップラ流速計。