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
Links
- 238000004364 calculation method Methods 0.000 claims description 53
- 238000005314 correlation function Methods 0.000 claims description 36
- 238000005259 measurement Methods 0.000 claims description 27
- 238000007493 shaping process Methods 0.000 claims description 23
- 238000005070 sampling Methods 0.000 claims description 22
- 238000000034 method Methods 0.000 description 21
- 239000002245 particle Substances 0.000 description 17
- 238000010586 diagram Methods 0.000 description 13
- 238000012545 processing Methods 0.000 description 12
- 238000005311 autocorrelation function Methods 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 6
- 230000010287 polarization Effects 0.000 description 5
- 238000007796 conventional method Methods 0.000 description 2
- NJPPVKZQTLUDBO-UHFFFAOYSA-N novaluron Chemical compound C1=C(Cl)C(OC(F)(F)C(OC(F)(F)F)F)=CC=C1NC(=O)NC(=O)C1=C(F)C=CC=C1F NJPPVKZQTLUDBO-UHFFFAOYSA-N 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/26—Measuring 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図に示すように、先ず、
レーザトツプラ流速計の光電変換器が出力するド
ツプラ・バースト信号を受け、該ドツプラ・バー
スト信号を所定レベルでクリツピングして2値化
し、パルス信号を出力するクリツピング手段が設
けられている。そして、該クリツピング手段のパ
ルス信号を受け、予め設定されたクロツク信号に
より該パルス信号をサンプリングしてパルス信号
を出力するサンプリング手段が設けられている。
け、該パルス信号の相関関数を演算して相関信号
を出力するデジタル相関関数演算手段が設けられ
ている。加えて、該デジタル関関数演算手段の相
関信号を受け、該相関信号が予め設定された基準
レベルをクロス点を検出し、該クロス点に基づい
て相関信号の周波数を算出する周波数算出手段が
設けられると共に、該周波数算出手段の算出周波
数に基づいて被測定体の流速を算出する流速算出
手段が設けられた構成としている。
求項1の発明に加えて、上記デジタル相関関数演
算手段が出力する相関信号を滑らかな整形相関信
号に整形して周波数算出手段に出力する波形整形
手段が設けられた構成としている。
求項1又は2の発明における上記周波数算出手段
はデジタル相関関数演算手段が演算した相関関数
のサンプル数の1/2を基準レベルとした構成とし
ている。
えば、レーザ光源より予め設定された測定点に向
つて2本のレーザビームを照射収束させており、
この測定点で干渉縞が生じ、この干渉縞内を被測
定体、例えば、粒子が通過すると、該粒子によつ
て散乱光が生じ、該散乱光を集光して受光し、光
電変換してドツプラ・バースト信号を出力する。
そして、このドツプラ・バースト信号をクリツピ
ング手段が所定レベルでクリツピングして2値化
し、パルス信号を出力し、該パルス信号をサンプ
リング手段が予め設定されたクロツク信号により
サンプリングし、その結果をデジタル信号として
パルス信号を出力する。続いて、該パルス信号を
受けてデジタル相関関数演算手段が該パルス信号
と相関関数、例えば、自己相関関数を演算し、三
角波の相関信号を出力する。その際、請求項2に
係る発明では、上記相関信号がノイズで凹凸状に
なつているので、波形整形手段が、例えば、移動
平均法で滑らかな整形相関信号に整形して出力す
る。
形相関信号を受け、該相関信号が予め設定された
基準レベル、特に、請求項3に係る発明では、演
算に用いた相関信号のτ=0のときのサンプル数
の1/2を基準レベルとし、該基準レベルをクロス
するクロス点を検出し、該クロス点に基づいて上
記相関信号の周波数を算出する。そして、この算
出周波数がドツプラ・バースト信号の周波数に対
応しているので、流速算出手段が上記算出周波数
より粒子の流速を算出する。
ラ・バースト信号を2値化して分割すると共に、
相関関数を用いて導出した相関信号と基準レベル
のクロス点とより該相関信号の周波数を算出する
ようにしたために、SN比の悪いドツプラ信号に
対しても正確に流速を測定することができ、従来
に比して測定精度を著しく向上させることができ
るので、広範囲な流速測定を行うことができる。
時間で測定することができると共に、被測定体が
少なく乱れの大きい流れに対しても確実に流速を
測定することができる。
を波形整形するので、正確に周波数を算出するこ
とができることから、より測定精度を向上させる
ことができる。
相関関数のτ=0のときのサンプル数の1/2を基
準レベルとしているので、ノイズが多くて相関関
数の振幅が小さい場合でもクロス点を確実に検出
することができ、測定精度の向上を図ることがで
きる。
説明する。
計であつて、レーザビームの干渉性を利用してエ
ンジン内のガス流速やターボ機械の内部流速を測
定するものである。
どを送受光する差動型の光学系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が生じることになる。
変換しており、該散乱光51を集光する集光レン
ズ52と、該集光レンズ52で集光した散乱光5
1の伝播方向を変える反射鏡53と、上記散乱光
51を電気信号に変換してドツプラ・バースト信
号を出力する光電変換器54とより構成されてい
る。そして、このドツプラ・バースト信号の周波
数、つまり、散乱光の強弱変化が上記粒子の速度
に比例している。
に示すように、上記ドツプラ・バースト信号の周
波数より粒子速度を算出するように構成されてお
り、先ず、上記光電変換器54が出力するドツプ
ラ・バースト信号をバンドパスフイルタ31が所
定周波数帯域の信号のみを取り出し、所定のノイ
ズを除去した後、第4図aに示すドツプラ・バー
スト信号をクリツピング手段であるクリツピング
回路32に出力するように構成されている。該ク
リツピング回路32は、ノイズの影響を除去する
ため、ドツプラ・バースト信号のペデスタル又は
エンベロープを予め設定された基準レベルL1と
比較して信号処理の開始のトリガとしており、該
ドツプラ・バースト信号が所定レベルL1になる
と立上げ、ゼロクロスで2値化し、第4図bに示
すパルス信号をサンプリング手段であるサンプリ
ング回路33に出力するように成つている。
うにシフトレジスタで構成されると共に、クロツ
ク回路34が接続され、該クロツク回路34が、
例えば、500MHzのクロツク信号を出力しており、
該サンプリング回路33は上記クロツク信号でパ
ルス信号を一定周期でサンプリングし、第4図c
に示すパルス信号をデジタル相関関数演算手段で
あるデジタル相関関数演算回路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
に論理信号を出力するように成つている。
いドツプラ・バースト信号の場合、相関信号が凹
凸状になるので、相関信号の波形を滑らかにして
ノイズを除去し、整形相関信号を周波数算出手段
である周波数算出回路37に出力するように構成
されている。ここで、上記SN比(SNR)は次に
示すように定義する。
り相関信号を整形しており、次式に示すように、 (t)={R(t−1)+R(t)+R(t+1)}/3
…… で演算し、求める相関関数値(t)を前後の相関関
数値R(t−1)、R(t+1)を用いて平均化処
理するように構成されている。
うに、上記整形相関信号を予め設定された基準レ
ベル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に
設定している。
の有効性(バリデーシヨン:VALB)を判定す
るように構成されており、この有効性(VAL)
は次式に示すように、 VAL=|1−(S2/S1)|
…… S1:最初の三角波より求めた周波数 S2:最後の三角波より求めた周波数 であり、この有効性(VAL)が次式で充足する
か否かにより判定している。
関関数の周波数がドツプラ・バースト信号の周波
数fに対応しているので、該数波数fより粒子速
度Vを次式に基づいて導出するように構成されて
いる。
ル k→2:第2レーザビームの進行方向の単位ベクト
ル 尚、本実施例においては、波形整形回路36を
設けたが、ノイズの少ない場合該波形整形回路3
6を設けることなくデジタル相関信号を相関関係
演算回路35より周波数算出回路37に直接入力
するようにしてもよい。
定動について具体的な実験結果に基づいて説明す
る。
ーム41は偏波面回転器43で偏向特性が補正さ
れた後、ビームスプリツタ44で2本の平行なビ
ーム41a,41aに分割され、両ビーム41
a,41aは収束レンズ45を介して測定点0に
照射収束する。この測定点0を粒子が通過する
と、散乱光51が生じ、該散乱光51は収束レン
ズ45、集光レンズ52及び反射鏡53を介して
光電変換器54に入射し、該光電変換器54がド
ツプラ・バースト信号を出力する。
作について第5図に示す具体的な波形図に基づき
説明する。
ズが重畳しており、ペデスタル成分が完全に除去
されている場合、該ドツプラ・バースト信号は次
式で表わされる。
サンプリング周波数の比(以下、便宜上入力周
波数比という。) 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を算出する。
用いた数値実験結果を示しており、a〜eは入力
周波数比Wを5、10、20、30、40に各々設定した
ものであり、横軸にSN比、縦軸に有効性(%)
並びに精度を示している。尚、この精度は、入力
に対する出力の誤差で、1−(W′/W)で示し、
このW′は周波数算出回路37が算出した測定ド
ツプラ周波数に対するクロツク信号によるサンプ
リング周波数の比(以下、便宜上測定周波数比と
いう。)である。また、第6図の白丸は波形整形
しない場合で、黒丸は波形整形した場合を示して
いる。
Wが10〜20の範囲内で、SN比が−6dBまでは3
%以内の測定誤差と、80%以上の有効性とを得る
ことができ、図示しないが、従来、SN比が0dB
における測定誤差が16%であつたのに対し、著し
く精度が向上している。
周波数比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により波形整形すると、より精
度が向上することになる。
SN比における測定誤差及び有効性の特性を示し
ており、この第8図から明らかなように、入力周
波数比Wとして5〜40までの値で用いると、SN
比が−6dBまでのドツプラ・バースト信号を65%
以上の有効性でもつて処理できる。換言すれば、
500MHzのクロツク信号を用いると、12.5MHz〜
100MHzのドツプラ・バースト信号を処理するこ
とができる。更に、SN比が−10dBのドツプラ・
バースト信号を入力周波数比Wが10〜20の範囲内
で45%以上の有効性で処理でき、換言すれば、
500MHzのクロツク信号で、25MHz〜50MHzのド
ツプラ・バースト信号を処理することができる。
る標準偏差σの特性図であり、白丸が本発明によ
るもの、黒丸がカウンタ方式による従来例を示
し、この第9図より明らかなように、本発明によ
るものは標準偏差で小さく、測定精度が著しく向
上している。
してデジタル相関を演算し、その結果の相関関数
と基準レベルL2とのクロス点より該相関信号の
周波数を算出するようにしたために、SN比の悪
いドツプラ・バースト信号であつても正確に流速
を測定することができ、従来に比して測定精度を
著しく向上させることができると共に、ダイナミ
ツクレンジを拡大させることができるので、広範
囲な流速測定を行うことができる。
時間で測定することができると共に、粒子が少な
く乱れの大きい流れに対しても確実に流速を測定
することができる。
正確に周波数を算出することができ、より測定精
度を向上させることができる。
相関関係R(0)の値の1/2を基準レベルL2とし
ているので、クロス点を確実に検出でき、測定精
度の向上を計ることができる。
相関関数を利用して流速を求めるようにしたが、
2方向から散乱光51を受光し、相互相関関数を
利用してその位相を求めると、粒子の粒径と速度
とを同時に実時間で測定することができる。具体
的には、2つの散乱光を各々光電変換器54、バ
ンドパスフイルタ31、クリツピング回路32で
信号処理し、第3図一点鎖線で示すように、各シ
フトレジスタからなるサンプリリング回路33,
33を介してデジタル相関関数演算回路35の各
シフトレジスタ35a,35bにパルス信号を入
力させて、相互相関関数を演算するようにして、
エンジン内における燃料の粒径などを同時に実時
間で測定するようにしてもよい。
たが、乱流などの渦構造を把握するようにしても
よく、その際、流れ場における速度の時・空間相
関を求める必要がある。そのために、空気内の2
点もしくはそれ以上の点の速度を同時に測定し、
空間相関を求めなくてはならない。例えば、10点
の測定を同時に行う場合、従来の信号処理器では
10台必要になりシステムとしては膨大になる。し
かしながら、本発明によれば、それを、簡素化、
小型化することができる。つまり、第10図に示
すように、10台の光電変換器54からの信号をそ
れぞ二値化し、それを一旦高速メモリして、速度
とその時間情報をメモリに記憶させる。
クサを用いて作動し、各々の情報をコンピユータ
39で処理する。この実施例によれば、ICは
17MHzで作動するので、256ビツトのデータを
88.5kHzのデータレートで処理できる。10台並列
にあると約1/10として、8KHzのデータレートに
対応することができる。実際の測定で得られるド
ツプラ・バースト信号のデータレートは1kHz以
下であるので、この例示した8kHz以下は実時間
ということになる。従つて、光電変換器54等を
10台備えれば良く、一台のの装置で10チヤンネル
に対応できることになり、コンパクト化を著しく
図ることができる。
前方散乱形でもよく、また、参照光方式や単一ビ
ーム方式などでもよい。
関数演算回路35はシフトレジスタに限られるも
のではなく、また、波形整形回路36は移動平均
法に限られるものではない。
る。第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 レーザドツプラ流速計の光電変換器が出力す
るドツプラ・バースト信号を受け、該ドツプラ・
バースト信号を所定レベルでクリツピングして2
値化し、パルス信号を出力するクリツピング手段
と、 該クリツピング手段のパルス信号を受け、予め
設定されたクロツク信号により該パルス信号をサ
ンプリングしてパルス信号を出力するサンプリン
グ手段と、 該サンプリング手段のパルス信号を受け、該パ
ルス信号の相関関数を演算して相関信号を出力す
るデジタル相関関数演算手段と、 該デジタル相関関数演算手段の相関信号を受
け、該相関信号が予め設定された基準レベルをク
ロスするクロス点を検出し、該クロス点に基づい
て相関信号の周波数を算出する周波数算出手段
と、 該周波数算出手段の算出周波数に基づいて被測
定体の流速を算出する流速算出手段とを備えてい
ることを特徴とするレーザドツプラ流速計の流速
測定装置。 2 レーザドツプラ流速計の光電変換器が出力す
るドツプラ・バースト信号を受け、該ドツプラ・
バースト信号を所定レベルでクリツピングして2
値化し、パルス信号を出力するクリツピング手段
と、 該クリツピング手段のパルス信号を受け、予め
設定されたクロツク信号により該パルス信号をサ
ンプリングしてパルス信号を出力するサンプリン
グ手段と、 該サンプリング手段のパルス信号を受け、該パ
ルス信号の相関関数を演算して相関信号を出力す
るデジタル相関関数演算手段と、 該デジタル相関関数演算手段の相関信号を受
け、該相関信号を滑らかな整形相関信号に整形し
て出力する波形整形手段と、 該波形整形手段の整形相関信号を受け、該整形
相関信号が予め設定された基準レベルをクロスす
るクロス点を検出し、該クロス点に基づいて整形
相関信号の周波数を算出する周波数算出手段と、 該周波数算出手段の算出周波数に基づいて被測
定体の流速を算出する流速算出手段とを備えてい
ることを特徴とするレーザドツプラ流速計の流速
測定装置。 3 周波数算出手段は、デジタル相関関数演算手
段が演算した相関関数のサンプル数の1/2を基準
レベルとしていることを特徴とする請求項1又は
2記載のレーザドツプラ流速計の流速測定装置。
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)
| 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)
| 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 |
-
1989
- 1989-01-11 JP JP1005092A patent/JPH02236171A/ja active Granted
- 1989-07-10 US US07/378,098 patent/US5005144A/en not_active Expired - Lifetime
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 |