JPH0559383B2 - - Google Patents

Info

Publication number
JPH0559383B2
JPH0559383B2 JP1005092A JP509289A JPH0559383B2 JP H0559383 B2 JPH0559383 B2 JP H0559383B2 JP 1005092 A JP1005092 A JP 1005092A JP 509289 A JP509289 A JP 509289A JP H0559383 B2 JPH0559383 B2 JP H0559383B2
Authority
JP
Japan
Prior art keywords
signal
correlation
frequency
calculation means
pulse signal
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 - Lifetime
Application number
JP1005092A
Other languages
English (en)
Other versions
JPH02236171A (ja
Inventor
Takeshi Nakajima
Juji Ikeda
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.)
Ono Sokki Co Ltd
Original Assignee
Ono Sokki 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 Ono Sokki Co Ltd filed Critical Ono Sokki Co Ltd
Priority to JP1005092A priority Critical patent/JPH02236171A/ja
Priority to US07/378,098 priority patent/US5005144A/en
Publication of JPH02236171A publication Critical patent/JPH02236171A/ja
Publication of JPH0559383B2 publication Critical patent/JPH0559383B2/ja
Granted legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/26Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting optical wave

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optical Radar Systems And Details Thereof (AREA)
  • Measuring Volume Flow (AREA)

Description

【発明の詳細な説明】 (産業上の利用分野) 本発明は、例えば、エンジン内のガス流速やタ
ーボ機械の内部流速などをレーザビームの干渉を
利用して測定するレーザドツプラ流速計の流速測
定装置に関し、特に、ドツプラ・バースト信号の
信号処理系に係るものである。
(従来の技術) 一般に、この種のレーザドツプラ流速計は、レ
ーザ光源より分光器や収束器を介して2本のレー
ザビームを予め設定された測定点に照射収束さ
せ、該測定点に生じる干渉縞内を粒子(被測定
体)が通過すると散乱光が生じ、該散乱光を集光
器を介して光電変換器で受光し、該光電変換器が
出力するドツプラ・バースト信号(バースト状の
ドツプラ信号)を信号処理系で処理して上記粒子
の速度を測定するように構成されている。
この信号処理系は、従来、カウンタ方式或いは
トラツカ方式が利用されており、該カウンタ方式
は入力周波数の周期の数を内蔵したクロツクで計
測し、逆数演算を行つてドツプラ周波数を求めて
いる。一方、上記トラツカ方式は入力周波数と内
蔵した電圧制御発振器の周波数とが常に一定にな
るようにした周波数追尾回路で周波数負帰還をか
け、この追尾周波数よりドツプラ周波数が求めて
いる。そして、上記ドツプラ周波数が粒子の速度
に比例しているので、該ドツプラ周波数より流速
を算出するようにしている。
(発明が解決しようとする課題) 上述したレーザドツプラ流速計において、信号
処理系にカウンタ方式やトラツカ方式を採用して
いるが、上記ドツプラ・バースト信号にはシヨツ
トノイズやホワイトノイズなどの各種のノイズが
重畳しており、SN比の悪いドツプラ信号の場合、
上記両方式の何れにおいても流速を測定すること
ができず、測定範囲が狭いという問題があつた。
更に、上記カウンタ方式にあつては乱流の渦スケ
ールなどを求める際に必要となる時系列データを
得ることができないという欠点があり、一方、上
記トラツカ方式にあつてはアナログ処理であるの
で、粒子数が少なく乱れの大きい流れの場合に適
用することができないという欠点があつた。
本発明は、斯かる点に鑑みてなされたもので、
ドツプラ・バースト信号を2値化して相関関数演
算し、この相関信号の周波数を算出することによ
り、SN比の悪いドツプラ・バースト信号に対し
ても高精度で且つ実時間で流速の測定を可能にす
ることを目的とするものである。
(課題を解決するための手段) 上記目的を達成するために、請求項1に係る発
明が講じた手段は、第1図に示すように、先ず、
レーザトツプラ流速計の光電変換器が出力するド
ツプラ・バースト信号を受け、該ドツプラ・バー
スト信号を所定レベルでクリツピングして2値化
し、パルス信号を出力するクリツピング手段が設
けられている。そして、該クリツピング手段のパ
ルス信号を受け、予め設定されたクロツク信号に
より該パルス信号をサンプリングしてパルス信号
を出力するサンプリング手段が設けられている。
更に、該サンプリング手段のパルス信号を受
け、該パルス信号の相関関数を演算して相関信号
を出力するデジタル相関関数演算手段が設けられ
ている。加えて、該デジタル関関数演算手段の相
関信号を受け、該相関信号が予め設定された基準
レベルをクロス点を検出し、該クロス点に基づい
て相関信号の周波数を算出する周波数算出手段が
設けられると共に、該周波数算出手段の算出周波
数に基づいて被測定体の流速を算出する流速算出
手段が設けられた構成としている。
また、請求項2に係る発明が講じた手段は、請
求項1の発明に加えて、上記デジタル相関関数演
算手段が出力する相関信号を滑らかな整形相関信
号に整形して周波数算出手段に出力する波形整形
手段が設けられた構成としている。
また、請求項3に係る発明が講じた手段は、請
求項1又は2の発明における上記周波数算出手段
はデジタル相関関数演算手段が演算した相関関数
のサンプル数の1/2を基準レベルとした構成とし
ている。
(作用) 上記構成により、請求項1に係る発明では、例
えば、レーザ光源より予め設定された測定点に向
つて2本のレーザビームを照射収束させており、
この測定点で干渉縞が生じ、この干渉縞内を被測
定体、例えば、粒子が通過すると、該粒子によつ
て散乱光が生じ、該散乱光を集光して受光し、光
電変換してドツプラ・バースト信号を出力する。
そして、このドツプラ・バースト信号をクリツピ
ング手段が所定レベルでクリツピングして2値化
し、パルス信号を出力し、該パルス信号をサンプ
リング手段が予め設定されたクロツク信号により
サンプリングし、その結果をデジタル信号として
パルス信号を出力する。続いて、該パルス信号を
受けてデジタル相関関数演算手段が該パルス信号
と相関関数、例えば、自己相関関数を演算し、三
角波の相関信号を出力する。その際、請求項2に
係る発明では、上記相関信号がノイズで凹凸状に
なつているので、波形整形手段が、例えば、移動
平均法で滑らかな整形相関信号に整形して出力す
る。
その後、周波数算出手段が上記相関信号又は整
形相関信号を受け、該相関信号が予め設定された
基準レベル、特に、請求項3に係る発明では、演
算に用いた相関信号のτ=0のときのサンプル数
の1/2を基準レベルとし、該基準レベルをクロス
するクロス点を検出し、該クロス点に基づいて上
記相関信号の周波数を算出する。そして、この算
出周波数がドツプラ・バースト信号の周波数に対
応しているので、流速算出手段が上記算出周波数
より粒子の流速を算出する。
(発明の効果) 従つて、請求項1に係る発明によれば、トツプ
ラ・バースト信号を2値化して分割すると共に、
相関関数を用いて導出した相関信号と基準レベル
のクロス点とより該相関信号の周波数を算出する
ようにしたために、SN比の悪いドツプラ信号に
対しても正確に流速を測定することができ、従来
に比して測定精度を著しく向上させることができ
るので、広範囲な流速測定を行うことができる。
また、従来のカウンタ方式等に比し、流速を実
時間で測定することができると共に、被測定体が
少なく乱れの大きい流れに対しても確実に流速を
測定することができる。
また、請求項2に係る発明によれば、相関信号
を波形整形するので、正確に周波数を算出するこ
とができることから、より測定精度を向上させる
ことができる。
また、請求項3の発明によれば、演算に用いた
相関関数のτ=0のときのサンプル数の1/2を基
準レベルとしているので、ノイズが多くて相関関
数の振幅が小さい場合でもクロス点を確実に検出
することができ、測定精度の向上を図ることがで
きる。
(実施例) 以下、本発明の実施例を図面に基づいて詳細に
説明する。
第2図に示すように、1はレーザドツプラ流速
計であつて、レーザビームの干渉性を利用してエ
ンジン内のガス流速やターボ機械の内部流速を測
定するものである。
該レーザドツプラ流速計1は、レーザビームな
どを送受光する差動型の光学系2と、該光学系2
より出力されるドツプラ・バースト信号(バース
ト状のドツプラ信号)を処理する信号処理系3と
より構成されている。該光学系2は後方散乱形で
あつて、レーザビームの処理手段4と散乱光の受
光手段5とより構成されている。更に、該照射手
段4は、レーザビーム41を発するレーザ光源4
2と、該レーザビーム41の偏光特性を補正する
偏波面回転器43と、該偏波面回転器43で補正
された1本のレーザビーム41を2本の平行なビ
ーム41a,41aに分けるビームスプリツタ4
4と、該両ビーム41a,41aを予め設定され
た測定点0に照射収束させる収束レンズ45とよ
り構成されている。そして、該測定点0におい
て、上記両ビーム41a,41aの干渉により干
渉縞が生じ、該干渉縞内を被測定体である粒子が
通過すると、散乱光51が生じることになる。
上記受光手段5は該散乱光51を受光して光電
変換しており、該散乱光51を集光する集光レン
ズ52と、該集光レンズ52で集光した散乱光5
1の伝播方向を変える反射鏡53と、上記散乱光
51を電気信号に変換してドツプラ・バースト信
号を出力する光電変換器54とより構成されてい
る。そして、このドツプラ・バースト信号の周波
数、つまり、散乱光の強弱変化が上記粒子の速度
に比例している。
そこで、上記信号処理系3は、第4図の波形図
に示すように、上記ドツプラ・バースト信号の周
波数より粒子速度を算出するように構成されてお
り、先ず、上記光電変換器54が出力するドツプ
ラ・バースト信号をバンドパスフイルタ31が所
定周波数帯域の信号のみを取り出し、所定のノイ
ズを除去した後、第4図aに示すドツプラ・バー
スト信号をクリツピング手段であるクリツピング
回路32に出力するように構成されている。該ク
リツピング回路32は、ノイズの影響を除去する
ため、ドツプラ・バースト信号のペデスタル又は
エンベロープを予め設定された基準レベルL1
比較して信号処理の開始のトリガとしており、該
ドツプラ・バースト信号が所定レベルL1になる
と立上げ、ゼロクロスで2値化し、第4図bに示
すパルス信号をサンプリング手段であるサンプリ
ング回路33に出力するように成つている。
該サンプリング回路32は、第3図にも示すよ
うにシフトレジスタで構成されると共に、クロツ
ク回路34が接続され、該クロツク回路34が、
例えば、500MHzのクロツク信号を出力しており、
該サンプリング回路33は上記クロツク信号でパ
ルス信号を一定周期でサンプリングし、第4図c
に示すパルス信号をデジタル相関関数演算手段で
あるデジタル相関関数演算回路35に出力するよ
うに構成されている。
該デジタル相関関数演算回路35は、上記パル
ス信号の自己相関関数R(τ)を演算して、第4
図dに示す三角波の相関信号を波形整形手段であ
る波形整形回路36に出力するように構成されて
いる。該デジタル相関関数演算回路35は、第3
図にも示すように、2つの相関用第1シフトレジ
スタ35a及び第2シフトレジスタ35bと該両
シフトレジスタ35a,35bのビツト数NT
対応した複数個のエクスクルーシブノア回路35
c,35c,……とより構成され、該両シフトレ
ジスタ35a,35bは、例えば、64ビツトのシ
フトレジスタで、上記サンプリング回路33のシ
フトレジスタに直列に接続される一方、上記エク
スクルシーブノア回路35c,35c,……は両
シフトレジスタ35a,35bの各ビツトからの
出力信号を受けて動作し、上記波形整形回路36
に論理信号を出力するように成つている。
上記波形整形回路36は、SN比(SNR)の悪
いドツプラ・バースト信号の場合、相関信号が凹
凸状になるので、相関信号の波形を滑らかにして
ノイズを除去し、整形相関信号を周波数算出手段
である周波数算出回路37に出力するように構成
されている。ここで、上記SN比(SNR)は次に
示すように定義する。
SNR=20・log(0.47・A/α) …… A:ドツプラ信号の最大振幅 0.47・A:ドツプラ・バースト信号の実効値 α:ノイズの実効値 そして、上記波形整形回路36は移動平均法によ
り相関信号を整形しており、次式に示すように、 (t)={R(t−1)+R(t)+R(t+1)}/3
…… で演算し、求める相関関数値(t)を前後の相関関
数値R(t−1)、R(t+1)を用いて平均化処
理するように構成されている。
上記周波数算出回路37は、第4図dに示すよ
うに、上記整形相関信号を予め設定された基準レ
ベルL2とクロスするクロス点より該整形相関信
号の周波数を算出して流速算出手段である流速算
出回路38に出力するように構成されている。つ
まり、上記整形相関信号と基準レベルL2との2
つのクロス点Mi、Mi+1の平均値(Mi+Mi+
1)/2よりピーク位置Piを算出し、このピーク
位置Piをカウントして上記整形相関信号の周波数
を算出するようにしている。そして、上記周波数
算出回路37の基準レベルL2は自己相関関数R
(τ)の遅れ時間τが零のときの値R(0)(該値
R(0)は相関関数演算のサンプル数NTに等し
い)の1/2の値に設定されている。すなわち、上
記デジタル相関関数演算回路35の両シフトレジ
スタ35a,35bにおける各ビツトの出力状態
が0−0及び1−1になる確率と0−1及び1−
0になる確率とは共に1/2であり、自己相関関数
はR(0)/2(=NT/2)を中心に三角波を形
成するので、上記基準レベルL2をR(0)/2に
設定している。
更に、上記周波数算出回路37は、算出周波数
の有効性(バリデーシヨン:VALB)を判定す
るように構成されており、この有効性(VAL)
は次式に示すように、 VAL=|1−(S2/S1)|
…… S1:最初の三角波より求めた周波数 S2:最後の三角波より求めた周波数 であり、この有効性(VAL)が次式で充足する
か否かにより判定している。
VAL<β …… β:設定された許容誤差 そして、上記流速算出手段38は、上記整形相
関関数の周波数がドツプラ・バースト信号の周波
数fに対応しているので、該数波数fより粒子速
度Vを次式に基づいて導出するように構成されて
いる。
f=(V→/λ)・(k→1−k→2)…… λ:レーザビームの波長 k→1:第1レーザビームの進行方向の単位ベクト
ル k→2:第2レーザビームの進行方向の単位ベクト
ル 尚、本実施例においては、波形整形回路36を
設けたが、ノイズの少ない場合該波形整形回路3
6を設けることなくデジタル相関信号を相関関係
演算回路35より周波数算出回路37に直接入力
するようにしてもよい。
次に、このレーザ・トツプラ流速計1の流速測
定動について具体的な実験結果に基づいて説明す
る。
先ず、レーザ光源42より出射されたレーザビ
ーム41は偏波面回転器43で偏向特性が補正さ
れた後、ビームスプリツタ44で2本の平行なビ
ーム41a,41aに分割され、両ビーム41
a,41aは収束レンズ45を介して測定点0に
照射収束する。この測定点0を粒子が通過する
と、散乱光51が生じ、該散乱光51は収束レン
ズ45、集光レンズ52及び反射鏡53を介して
光電変換器54に入射し、該光電変換器54がド
ツプラ・バースト信号を出力する。
続いて、このドツプラ・バースト信号の処理動
作について第5図に示す具体的な波形図に基づき
説明する。
一般的に、上記ドツプラ・バースト信号はノイ
ズが重畳しており、ペデスタル成分が完全に除去
されている場合、該ドツプラ・バースト信号は次
式で表わされる。
B(t)=A・exp[−{(2√2・t)/W・Nf}2] ・cos{(2π・t/W)+ψ} …… B(t):ドツプラ・バースト信号 t:サンプリング周期で無次元化された時間 A:振幅を示す定数 W:ドツプラ周波数に対するクロツク信号による
サンプリング周波数の比(以下、便宜上入力周
波数比という。) Nf:フリンジ数 ψ:位相遅れ そこで、第5図aはホワイトノイズをノイズと
して第式に加えたドツプラ・バースト信号を示
しており、この第5図の各波形は、SN比が−
3dB、入力周波数比Wが20、フリンジ数Nfが12、
最大振幅Aが1.0、トリガレベルL1が0.3、デジタ
ル相関関数演算回路35のビツト数NTが128、
無次元最大遅延時間の設定地Nτが128とした場合
の実験結果である。そして、このドツプラ・バー
スト信号をバンドパスフイルタ31を通した後、
該ドツプラ・バースト信号のベデスタル成分をト
リガとして用い、クリツピング回路32が上記ド
ツプラ・バースト信号を所定レベルでクリツピン
グして2値化し、パルス信号を出力する(第4図
b参照)、続いて、該パルス信号をサンプリング
回路33がクロツク回路34のクロツク信号によ
り一定周期でサンプリングし、第5図bに示すよ
うに、所定の数のパルス信号を出力する。その
際、第5図bおけるZで示すように、ノズルによ
るパルス信号も出力されることになる。
その後、上記パルス信号を受けてデジタル相関
関数演算回路35がデジタル信号の自己相関関数
R(τ)を演算し、第5図cに示すように、相関
信号を出力する。この相関信号は三角波となると
共に、上記ノイズにより凹凸状となつているの
で、該相関信号を波形整形回路36が第式に基
づく移動平均法により波形整形し、第5図dに整
形相関信号を出力する。
次いで、上記整形相関信号を受けて周波数算出
回路37が、τ=0の自己相関関数R(0)の値
1/2の基準レベルL2とする基準線のクロス点Miを
検出し、相関信号のピーク位置Piを算出してカウ
ントし、相関信号の周波数を算出する。そして、
上記相関信号の有効性(VAL)を第式に基づ
いて判定し、該有効性(VAL)が第式を充足
する場合、周波数算出回路37の算出周波数に基
づいて流速算出回路38が第式により粒子速度
Vを算出する。
第6図a〜eは上記レーザドツプラ流速計1を
用いた数値実験結果を示しており、a〜eは入力
周波数比Wを5、10、20、30、40に各々設定した
ものであり、横軸にSN比、縦軸に有効性(%)
並びに精度を示している。尚、この精度は、入力
に対する出力の誤差で、1−(W′/W)で示し、
このW′は周波数算出回路37が算出した測定ド
ツプラ周波数に対するクロツク信号によるサンプ
リング周波数の比(以下、便宜上測定周波数比と
いう。)である。また、第6図の白丸は波形整形
しない場合で、黒丸は波形整形した場合を示して
いる。
この第6図より明らかなように、入力周波数比
Wが10〜20の範囲内で、SN比が−6dBまでは3
%以内の測定誤差と、80%以上の有効性とを得る
ことができ、図示しないが、従来、SN比が0dB
における測定誤差が16%であつたのに対し、著し
く精度が向上している。
また、第7図a,bは、SN比が−6dB、入力
周波数比Wが20、フリンジ数Nfが12において、
100個の実験デーのうち第式を充足した個数を
示すヒストグラムであり、第7図aは波形整形し
ない場内、第7図bは波形整形した場合であり、
第7図aにあつては平均測定周波数比W′が
19.31、準備偏差σが1.74であり、第7図bにあ
つては平均測定周波数比W′が19.91、標準偏差σ
が0.65である。この第7図より明らかなように、
波形整形回路36により波形整形すると、より精
度が向上することになる。
また、第8図は波形整形した場合において、各
SN比における測定誤差及び有効性の特性を示し
ており、この第8図から明らかなように、入力周
波数比Wとして5〜40までの値で用いると、SN
比が−6dBまでのドツプラ・バースト信号を65%
以上の有効性でもつて処理できる。換言すれば、
500MHzのクロツク信号を用いると、12.5MHz〜
100MHzのドツプラ・バースト信号を処理するこ
とができる。更に、SN比が−10dBのドツプラ・
バースト信号を入力周波数比Wが10〜20の範囲内
で45%以上の有効性で処理でき、換言すれば、
500MHzのクロツク信号で、25MHz〜50MHzのド
ツプラ・バースト信号を処理することができる。
また、第9図は従来例と比較したSN比に対す
る標準偏差σの特性図であり、白丸が本発明によ
るもの、黒丸がカウンタ方式による従来例を示
し、この第9図より明らかなように、本発明によ
るものは標準偏差で小さく、測定精度が著しく向
上している。
従つて、上記ドツプラ・バースト信号を2値化
してデジタル相関を演算し、その結果の相関関数
と基準レベルL2とのクロス点より該相関信号の
周波数を算出するようにしたために、SN比の悪
いドツプラ・バースト信号であつても正確に流速
を測定することができ、従来に比して測定精度を
著しく向上させることができると共に、ダイナミ
ツクレンジを拡大させることができるので、広範
囲な流速測定を行うことができる。
また、従来のカウンタ方式等に比し、流速を実
時間で測定することができると共に、粒子が少な
く乱れの大きい流れに対しても確実に流速を測定
することができる。
また、上記相関信号を波形整形するので、より
正確に周波数を算出することができ、より測定精
度を向上させることができる。
また、上記周波数算出手段37はτ=0の自己
相関関係R(0)の値の1/2を基準レベルL2とし
ているので、クロス点を確実に検出でき、測定精
度の向上を計ることができる。
尚、本実施例は、1つの光学系2を設け、自己
相関関数を利用して流速を求めるようにしたが、
2方向から散乱光51を受光し、相互相関関数を
利用してその位相を求めると、粒子の粒径と速度
とを同時に実時間で測定することができる。具体
的には、2つの散乱光を各々光電変換器54、バ
ンドパスフイルタ31、クリツピング回路32で
信号処理し、第3図一点鎖線で示すように、各シ
フトレジスタからなるサンプリリング回路33,
33を介してデジタル相関関数演算回路35の各
シフトレジスタ35a,35bにパルス信号を入
力させて、相互相関関数を演算するようにして、
エンジン内における燃料の粒径などを同時に実時
間で測定するようにしてもよい。
また、本実施例は一点の流速を求めるようにし
たが、乱流などの渦構造を把握するようにしても
よく、その際、流れ場における速度の時・空間相
関を求める必要がある。そのために、空気内の2
点もしくはそれ以上の点の速度を同時に測定し、
空間相関を求めなくてはならない。例えば、10点
の測定を同時に行う場合、従来の信号処理器では
10台必要になりシステムとしては膨大になる。し
かしながら、本発明によれば、それを、簡素化、
小型化することができる。つまり、第10図に示
すように、10台の光電変換器54からの信号をそ
れぞ二値化し、それを一旦高速メモリして、速度
とその時間情報をメモリに記憶させる。
その後に、相関関係演算回路35をマルチプレ
クサを用いて作動し、各々の情報をコンピユータ
39で処理する。この実施例によれば、ICは
17MHzで作動するので、256ビツトのデータを
88.5kHzのデータレートで処理できる。10台並列
にあると約1/10として、8KHzのデータレートに
対応することができる。実際の測定で得られるド
ツプラ・バースト信号のデータレートは1kHz以
下であるので、この例示した8kHz以下は実時間
ということになる。従つて、光電変換器54等を
10台備えれば良く、一台のの装置で10チヤンネル
に対応できることになり、コンパクト化を著しく
図ることができる。
尚また、上記光学系2は後方散乱形にしたが、
前方散乱形でもよく、また、参照光方式や単一ビ
ーム方式などでもよい。
また、サンプリング回路33及びデジタル相関
関数演算回路35はシフトレジスタに限られるも
のではなく、また、波形整形回路36は移動平均
法に限られるものではない。
【図面の簡単な説明】
第1図は本発明の構成を示すブロツク図であ
る。第2図〜第10図は本発明の実施例を示し、
第2図はレーザドツプラ流速計の概略構成図、第
3図はデジタル相関関数演算回路の概略構成図で
ある。第4図a,b,c,dはそれぞれドツプ
ラ・バースト信号、クリツピングされたパルス信
号、サンプリングされたパルス信号、相関信号を
示す波形図、第5図a,b,c,dはそれぞれ具
体的実験結果のドツプラ・バースト信号、サンプ
リングされたパルス信号、相関信号、整形相関信
号を示す波形図である。第6図a,d,c,d,
eはそれぞれ入力周波数比Wが異なる条件下にお
ける測定誤差及び有効性の特性図、第7図a,b
は波形整形しないものとしたものとにおける有効
性を示す柱状図、第8図は波形整形した場合の測
定誤差及び有効性の特性図、第9図は本発明と従
来例とを比較して示す標準偏差の特性図、第10
図は他の実施例を示すブロツク図である。 1……レーザドツプラ流速計、2……光学系、
3……信号処理系、4……照射手段、5……受光
手段、32……クリツピング回路、33……サン
プリング回路、34……クロツク回路、35……
デジタル相関関数演算回路、36……波形整形回
路、37……周波数算出回路、38……流速算出
回路、41……レーザビーム、42……レーザ光
源、51……散乱、54……光電変換器。

Claims (1)

  1. 【特許請求の範囲】 1 レーザドツプラ流速計の光電変換器が出力す
    るドツプラ・バースト信号を受け、該ドツプラ・
    バースト信号を所定レベルでクリツピングして2
    値化し、パルス信号を出力するクリツピング手段
    と、 該クリツピング手段のパルス信号を受け、予め
    設定されたクロツク信号により該パルス信号をサ
    ンプリングしてパルス信号を出力するサンプリン
    グ手段と、 該サンプリング手段のパルス信号を受け、該パ
    ルス信号の相関関数を演算して相関信号を出力す
    るデジタル相関関数演算手段と、 該デジタル相関関数演算手段の相関信号を受
    け、該相関信号が予め設定された基準レベルをク
    ロスするクロス点を検出し、該クロス点に基づい
    て相関信号の周波数を算出する周波数算出手段
    と、 該周波数算出手段の算出周波数に基づいて被測
    定体の流速を算出する流速算出手段とを備えてい
    ることを特徴とするレーザドツプラ流速計の流速
    測定装置。 2 レーザドツプラ流速計の光電変換器が出力す
    るドツプラ・バースト信号を受け、該ドツプラ・
    バースト信号を所定レベルでクリツピングして2
    値化し、パルス信号を出力するクリツピング手段
    と、 該クリツピング手段のパルス信号を受け、予め
    設定されたクロツク信号により該パルス信号をサ
    ンプリングしてパルス信号を出力するサンプリン
    グ手段と、 該サンプリング手段のパルス信号を受け、該パ
    ルス信号の相関関数を演算して相関信号を出力す
    るデジタル相関関数演算手段と、 該デジタル相関関数演算手段の相関信号を受
    け、該相関信号を滑らかな整形相関信号に整形し
    て出力する波形整形手段と、 該波形整形手段の整形相関信号を受け、該整形
    相関信号が予め設定された基準レベルをクロスす
    るクロス点を検出し、該クロス点に基づいて整形
    相関信号の周波数を算出する周波数算出手段と、 該周波数算出手段の算出周波数に基づいて被測
    定体の流速を算出する流速算出手段とを備えてい
    ることを特徴とするレーザドツプラ流速計の流速
    測定装置。 3 周波数算出手段は、デジタル相関関数演算手
    段が演算した相関関数のサンプル数の1/2を基準
    レベルとしていることを特徴とする請求項1又は
    2記載のレーザドツプラ流速計の流速測定装置。
JP1005092A 1989-01-11 1989-01-11 レーザドップラ流速計の流速測定装置 Granted JPH02236171A (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP1005092A JPH02236171A (ja) 1989-01-11 1989-01-11 レーザドップラ流速計の流速測定装置
US07/378,098 US5005144A (en) 1989-01-11 1989-07-10 Signal processing system for flow velocity in laser doppler velocimeter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP1005092A JPH02236171A (ja) 1989-01-11 1989-01-11 レーザドップラ流速計の流速測定装置

Publications (2)

Publication Number Publication Date
JPH02236171A JPH02236171A (ja) 1990-09-19
JPH0559383B2 true JPH0559383B2 (ja) 1993-08-30

Family

ID=11601749

Family Applications (1)

Application Number Title Priority Date Filing Date
JP1005092A Granted JPH02236171A (ja) 1989-01-11 1989-01-11 レーザドップラ流速計の流速測定装置

Country Status (2)

Country Link
US (1) US5005144A (ja)
JP (1) JPH02236171A (ja)

Families Citing this family (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5467320A (en) * 1993-01-08 1995-11-14 Halliburton Company Acoustic measuring method for borehole formation testing
DE4311182A1 (de) * 1993-04-06 1994-10-13 Bosch Gmbh Robert Verfahren zur berührungslosen Drehzahlmessung
DE4323553A1 (de) * 1993-07-14 1995-01-19 Deutsch Franz Forsch Inst Verfahren und Vorrichtung zur Signalverarbeitung
US5453834A (en) * 1994-04-11 1995-09-26 Evenstad; James O. Digital burst frequency translator
FR2749942B1 (fr) * 1996-06-18 1998-09-11 Sextant Avionique Sonde velocimetrique optique
FR2749943B1 (fr) * 1996-06-18 1998-09-11 Sextant Avionique Sonde velocimetrique optique
EP1717984B1 (en) * 1999-02-24 2010-05-19 Sony Deutschland Gmbh Receiving apparatus and synchronising method for a digital telecommunication system
US8861622B2 (en) * 1999-02-24 2014-10-14 Sony Deutschland Gmbh Transmitting apparatus and method for a digital telecommunication system
US6553518B1 (en) 1999-03-08 2003-04-22 International Business Machines Corporation Severe error detectors, methods and computer program products that use constellation specific error event thresholds to detect severe error events during demodulation of a signal comprising symbols from a plurality of symbol constellations
US6661847B1 (en) 1999-05-20 2003-12-09 International Business Machines Corporation Systems methods and computer program products for generating and optimizing signal constellations
US6839382B1 (en) 1999-10-29 2005-01-04 International Business Machines Corporation System, methods and computer program products for identifying digital impairments in modem signals using signature analysis and signal level comparison analysis
US6505222B1 (en) 1999-10-29 2003-01-07 International Business Machines Corporation Systems methods and computer program products for controlling undesirable bias in an equalizer
US6967995B1 (en) 1999-10-29 2005-11-22 International Business Machines Corporation Methods, systems and computer program products for carrier drop detection using a variable threshold
US6792004B1 (en) 1999-10-29 2004-09-14 International Business Machines Corporation Systems, methods and computer program products for averaging learned levels in the presence of robbed-bit signaling based on proximity
US6611563B1 (en) 1999-10-29 2003-08-26 International Business Machines Corporation Systems, methods and computer program products for data mode refinement of modem constellation points
US6816545B1 (en) 1999-10-29 2004-11-09 International Business Machines Corporation Systems, methods and computer program products for identifying digital impairments in modems based on clusters and/or skips in pulse code modulation signal levels
US6823017B1 (en) 1999-10-29 2004-11-23 International Business Machines Corporation Systems, methods and computer program products for filtering glitches from measured values in a sequence of code points
US6792040B1 (en) 1999-10-29 2004-09-14 International Business Machines Corporation Modems having a dual power mode capability and methods of operating same
US6650657B1 (en) 1999-10-29 2003-11-18 International Business Machines Corporation Systems, methods and computer program products for identifying digital impairments in modem signals
US6826157B1 (en) 1999-10-29 2004-11-30 International Business Machines Corporation Systems, methods, and computer program products for controlling data rate reductions in a communication device by using a plurality of filters to detect short-term bursts of errors and long-term sustainable errors
US6662322B1 (en) 1999-10-29 2003-12-09 International Business Machines Corporation Systems, methods, and computer program products for controlling the error rate in a communication device by adjusting the distance between signal constellation points
US6765955B1 (en) 1999-10-29 2004-07-20 International Business Machines Corporation Methods, systems and computer program products establishing a communication configuration for a modem connection to compensate for echo noise
US6823004B1 (en) 1999-10-29 2004-11-23 International Business Machines Corporation Methods, systems and computer program products for monitoring performance of a modem during a connection
US6754258B1 (en) 1999-10-29 2004-06-22 International Business Machines Corporation Systems, methods and computer program products for averaging learned levels in the presence of digital impairments based on patterns
CA2449551A1 (en) * 2003-11-17 2005-05-17 Photon Control Inc. Optical device and method for sensing multiphase flow
JP4364777B2 (ja) * 2004-12-02 2009-11-18 本田技研工業株式会社 内燃機関の空燃比制御装置
US7450222B1 (en) 2007-07-02 2008-11-11 The United States Of America As Represented By The United States Department Of Energy Correlated-intensity velocimeter for arbitrary reflector
FR2951275B1 (fr) * 2009-10-09 2012-11-02 Epsiline Dispositif de mesure de la vitesse du vent
JP5518110B2 (ja) * 2012-02-16 2014-06-11 株式会社堀場製作所 粒径分布測定装置
US20140073877A1 (en) * 2012-09-11 2014-03-13 Nellcor Puritan Bennett Llc Methods and systems for determining physiological information using autocorrelation and an envelope technique
JP5963087B2 (ja) * 2012-12-19 2016-08-03 イマジニアリング株式会社 計測器具、計測システム、及び計測方法
JP2019076641A (ja) 2017-10-27 2019-05-23 ソニー株式会社 計測装置、及び計測方法
CN114487467B (zh) * 2022-01-30 2023-04-18 中国工程物理研究院流体物理研究所 一种具有复测功能的激光干涉测速装置及方法
CN116298363B (zh) * 2023-01-13 2025-11-07 北京理工大学 一种激光多普勒冲击测试信号解调方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3842252A (en) * 1969-03-26 1974-10-15 Nat Res Dev Optical signal processing
US3777133A (en) * 1971-01-26 1973-12-04 C Wormald Cross correlator
US4019038A (en) * 1971-06-10 1977-04-19 Kent Instruments Limited Correlators
BE790341A (ja) * 1971-10-29 1973-04-20 Siderurgie Fse Inst Rech
US4244026A (en) * 1978-11-06 1981-01-06 General Electric Company Velocity measuring correlation sonar
US4596254A (en) * 1984-12-18 1986-06-24 Tsi Research Associates Limited Partnership Laser Doppler flow monitor
US4843564A (en) * 1987-04-23 1989-06-27 Tsi Incorporated Apparatus and method for measuring frequency of coherent component of a composite signal

Also Published As

Publication number Publication date
US5005144A (en) 1991-04-02
JPH02236171A (ja) 1990-09-19

Similar Documents

Publication Publication Date Title
JPH0559383B2 (ja)
US4540283A (en) Apparatus and method for determining the size and velocity of particles, droplets, bubbles or the like using laser light scattering
TWI651525B (zh) 粒子計數器
JPH05264434A (ja) 光学式サブミクロン粒子検出器
JPH02134540A (ja) 微粒子測定装置
US20230168079A1 (en) Method for obtaining the profile of a surface moving in relation to the system
US4397550A (en) Laser doppler velocimeter
CN103211588B (zh) 基于横向流速模型的血流流速成像方法
US4843564A (en) Apparatus and method for measuring frequency of coherent component of a composite signal
JP3235368B2 (ja) レーザドップラー速度測定装置
US5453834A (en) Digital burst frequency translator
US6311559B1 (en) Vibration measurement method and apparatus
JPH0290037A (ja) レーザ多焦点法による粒子の速度,径,屈折率の同時測定システム
CN108549086A (zh) 激光多普勒信号滤波带自适应选择及测试方法
JP2825480B2 (ja) 光学速度計プローブ
JPS59128449A (ja) 物体の移動速度検出方法
CN116148219B (zh) 一种基于涡旋光干涉的微小折射率变化率测量方法及系统
JP3282746B2 (ja) 半導体レーザデジタル振動変位計測装置
Lading Principles of laser anemometry
Lehmann et al. Features of a combined FFT and Hilbert transform for phase Doppler signal processing
JPH0933556A (ja) レーザ2焦点流速計
JPS6047965A (ja) レ−ザ・ドツプラ速度計
JP2735020B2 (ja) 集光ビームスポット径測定装置
Evenstad et al. A device for phase shift measurement in an advanced phase Doppler velocimeter
JPH0815436A (ja) エラー信号検出装置及びそれを用いたドップラー速度 計

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

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

Free format text: PAYMENT UNTIL: 20080830

Year of fee payment: 15

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

Free format text: PAYMENT UNTIL: 20090830

Year of fee payment: 16

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090830

Year of fee payment: 16