JP2004163889A - 入力信号の統計的特性を解析する装置、並びに入力信号を解析する集積回路および方法 - Google Patents
入力信号の統計的特性を解析する装置、並びに入力信号を解析する集積回路および方法 Download PDFInfo
- Publication number
- JP2004163889A JP2004163889A JP2003189764A JP2003189764A JP2004163889A JP 2004163889 A JP2004163889 A JP 2004163889A JP 2003189764 A JP2003189764 A JP 2003189764A JP 2003189764 A JP2003189764 A JP 2003189764A JP 2004163889 A JP2004163889 A JP 2004163889A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- crosslation
- function
- level
- event
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 41
- 230000004044 response Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 claims 2
- 230000006870 function Effects 0.000 description 133
- 230000000630 rising effect Effects 0.000 description 23
- 238000005070 sampling Methods 0.000 description 23
- 230000008569 process Effects 0.000 description 21
- 210000004027 cell Anatomy 0.000 description 19
- 239000000523 sample Substances 0.000 description 16
- 238000013139 quantization Methods 0.000 description 10
- 238000005311 autocorrelation function Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000012935 Averaging Methods 0.000 description 5
- 238000011835 investigation Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000002441 reversible effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 3
- 238000012806 monitoring device Methods 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000003750 conditioning effect Effects 0.000 description 2
- 239000013256 coordination polymer Substances 0.000 description 2
- 238000005314 correlation function Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000012905 input function Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 210000000352 storage cell Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06G—ANALOGUE COMPUTERS
- G06G7/00—Devices in which the computing operation is performed by varying electric or magnetic quantities
- G06G7/12—Arrangements for performing computing operations, e.g. operational amplifiers
- G06G7/14—Arrangements for performing computing operations, e.g. operational amplifiers for addition or subtraction
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Error Detection And Correction (AREA)
- Analogue/Digital Conversion (AREA)
Abstract
【解決手段】所定の傾斜をもって所定のレベルとの交差に対応するイベントを突きとめるために信号が評価される。信号の複数のバージョンであって、検出されたイベントの間隔に対応する量だけ相互に関してシフトされた複数のバージョンが結合される。これにより、得られる表示の形状により、その信号に関する統計的な情報が提供される。
【選択図】図1
Description
【0001】
本発明は、信号の統計的特性を決定する方法および装置に関し、特に、時間連続的にランダムな、無秩序な、または、不規則に振舞う信号の特性化に専ら適用できるものである。
【背景技術】
【0002】
例えば、信号の分類や信号の挙動のモニタリングや予測のために、その信号の統計的な特性を解析することが必要となる多くの状況がある。以下により詳細に説明するように、このような決定が有用である一例は、例えば暗号技術に使用される乱数生成における決定である。ランダムな、または、無秩序な雑音信号は、デジタイザに供給され、デジタイザは、所定のサンプリング間隔で信号をサンプリングし、乱数から成る信号のデジタル表示を出力することができる。効率のためには、サンプリング間隔は短くあるべきである。しかし、短いサンプリング間隔は、相互に統計的に独立でない乱数を生成する可能性がある。したがって、統計的に独立な乱数を生成するために必要な最小のサンプリング間隔の決定を可能にするために、雑音信号の統計的な特性を解析することが望ましい。
【0003】
他にも信号の統計データの決定が有用である多くの状況がある。信号が信号源の物理パラメータの変動を表す場合には、統計的解析は、その信号源を分類することに使用できる。例えば、各信号が画像中の変動を表しており、統計的な評価がその画像の対象を分類することに使用できる可能性がある。同様に、統計的解析が、話し言葉や音楽などの音の分類に使用できる可能性がある。
【0004】
既知の解析技法には、周波数領域(または、スペクトル)法や時間領域法が含まれる。時間領域法は、要求される情報を提供するためにしばしば必要であり、通常は信号の自己相関に基づいている。
【考案の開示】
【発明が解決しようとする課題】
【0005】
しかし、従来の相関技法は、関心の信号がガウシアンであり、時間の進行方向で見た場合の信号の統計的挙動が時間の逆進行方向で見た場合の信号の統計的挙動に対応するという暗黙の仮定に基づいている。すなわち、挙動のいかなる非対称性も、相関関数が時間の方向に影響されないという事実によって失われる。現実には、モニタされる信号の多くは、実際に非ガウシアンである。このような信号における非線形依存性は、標準的な相関技法によっては検出できない。
【0006】
したがって、従来技術よりも有用な結果を提供する信号の統計的挙動を解析するための方法および装置が提供されることが望ましい。
【課題を解決するための手段】
【0007】
本発明のいくつかの態様が、添付の特許請求の範囲において記載されている。
【0008】
他の態様によれば、複数のイベントを検出するために信号が調べられ、それぞれのイベントは、あるしきい値レベルと交わるときに所定の傾斜をとる信号に対応する。(好ましい実施形態では、信号は、例えば傾斜が負と異なって正である場合に、所定の傾斜を有すると見なされる。このように、信号がそのレベルを上昇させつつしきい値と交差する場合(すなわち、各「上昇交差」において)、あるいは、信号がそのレベルを下降させつつしきい値と交差する場合(すなわち、各「下降交差」)に各イベントが生じる。
【0009】
信号の複数のバージョン(好ましくは同一の複製)がその単一の信号から抽出され、それぞれのバージョンが他のバージョンにおけるそれぞれの異なるイベントと一致するイベントを包含するように、それぞれ相互にシフトしている。そして、複数のバージョンは、例えば平均化(ここでは「平均化」の語は本明細書において加算を包含することが意図されている)によって結合される。
【0010】
結果としての関数は、検出されたイベントの前または後の信号の平均的挙動の測度である。便宜上、この関数を本明細書中では「クロスレーション関数」と称するものとし、そのような関数を導出するために構成された装置を「クロスレータ」と称するものとする。関数は、関数が基づいているイベントが上昇交差の場合には、「前方クロスレーション」関数と称され、関数が基づいているイベントが下降交差の場合には、「後方クロスレーション」関数と称される。
【0011】
しきい値レベル、および、クロスレーション関数が基づいているイベントのタイプによって決まるであろう、信号のクロスレーション関数の形状は、入力信号に関する有用な情報を含んでいるであろう。(それぞれのイベントが結合される時点として定義される)原点に対する所与の時点において、関数の振幅は、各イベントに関する対応する時間における、特定の値に向かう入力信号の偏りを表す。
【0012】
更に、異なるクロスレーション(特に前方、および、後方クロスレーション)の形状間の関係は、他の有用な情報を含む。信号が時間的に反転した場合、下降交差は上昇交差に等価であることが理解されるだろう。したがって、時間反転可能な信号は、任意の所与のしきい値レベルに対して、対称な前方、および、後方クロスレーション関数を示すことになる。よって、これらの関数間の関係は、入力信号の時間反転性についての標識となる。
【0013】
更に、1つまたは複数のクロスレーション関数の形状の変化は、入力信号の性質に関する有用な情報をも含んでいる場合がある。
【0014】
したがって、本発明の装置は、1つまたは複数のクロスレーション関数の形状によって決まる1つまたは複数のパラメータを抽出して、入力信号の統計的特性を示す値、あるいは、一連の値を提供することが好ましい。
【0015】
例えば、以下で説明される実施形態では、前方、および、後方クロスレーション関数が調査されて、乱数生成のために、ランダムな入力信号のサンプリングに使用されるサンプリングパルス間の間隔に対応する時点におけるこれらの関数の振幅が決定される。この振幅が入力信号の平均値から有意な程度外れている場合には、この間隔でのサンプリングにより、連続するサンプルの値に、それらの独立性を低下させる偏りが生じることが示唆される。したがって、この解析装置の出力は、この好ましくない状況を告知し、または、修正することに使用できる。
【発明を実施するための最良の形態】
【0016】
図1を参照すると、図1は、本発明による信号解析装置を使用する乱数発生器を示す。
【0017】
この乱数発生器は、無秩序な出力信号x(t)を生成する物理的ランダム信号源(PRS)を有する。典型的な信号x(t)の波形が図2a)および図2b)のそれぞれに示されている。
【0018】
信号x(t)は、サンプリングパルス生成器(SPG)からのサンプリングパルスをも受信するアナログ−デジタル変換器(ADC)に送られる。無秩序な信号x(t)は、サンプリングパルス間の周期に対応する間隔でサンプラー(SMP)によりサンプリングされ、それぞれのアナログ出力は、振幅量子化器(QUA)に印加される。量子化器は、アナログ入力サンプルが比較されるJ個の異なる量子化レベルを生成する。出力(OP)においては、アナログサンプルのレベルに応じてデジタルの数値が生成される。
【0019】
したがって、乱数発生器は、サンプリングパルス間の周期に対応する間隔で、0からJの範囲で分布する乱数を生成する。
【0020】
ここまでに述べたシステムは既知である。図1の実施形態においては、モニタリング装置(MON)が設けられている。このモニタリング装置(MON)は、以下に詳述するように、無秩序な信号x(t)と量子化器(QUA)からの1からJの量子化レベルとを受信し、乱数が統計的に独立であることが期待できるか否かを示すモニタ出力(MOP)を生成する。
【0021】
モニタリング装置(MON)は図7に示されており、本発明による信号解析装置(本明細書ではクロスレータとも称する)(CRS)を備えている。このクロスレータ(CRS)は、信号x(t)を受信し、パラレル−シリアル変換器(PTS)を介して1からJのそれぞれの量子化レベル信号をも継続して受信する。クロスレータ(CRS)は、(以下説明するように)クロスレーション関数を出力(CFO)において時間シフト比較器(TSC)に出力する。この時間シフト比較器(TSC)は、統計的に独立なサンプルを得るために必要な最小サンプリング間隔を示す信号MSIを導出する。比較器(CMP)は、この値を現在のサンプリングパルス間隔を示す値SPIと比較する。比較器は、現在のサンプリングパルス間隔が計算された最小サンプリング間隔を超えているか否かを示すモニタ出力(MOP)を生成する。これは、正確な動作のためには超えているべきだからである。
【0022】
図2ないし図6を参照しつつ、クロスレータ(CRS)の動作原理を説明する。
【0023】
図2a)を参照すれば、時間的に連続したランダムな、無秩序な、あるいは、他の不規則なプロセスを示す信号x(t)と一定レベル(しきい値)の値Lとが示されている。そこには、信号x(t)が正の傾斜をもってレベルLに交差する複数の時間的瞬間が存在する。そのような時間的瞬間
t+ 1,t+ 2,・・・,t+ k-1,t+ k,t+ k+1,・・・
は、レベルLについての一組の上昇交差を形成し、これらの上昇交差は図2a)においてドットで示されている。
【0024】
これらの上昇交差の任意の一つ、例えばt+ kのものを選択し、時間的瞬間t+ kの前後の信号x(t)を考える。t+ kにおける上昇交差に関連付けられた信号軌跡x+ k(τ)は、τを相対時間として、
x+ k(τ)=x(t+ k+τ)
で定義される。したがって、図3に示される選択された軌跡x+ k(τ)は、調査対象である信号x(t)の単なる時間シフトされた複製である。基本信号x(t)の時間原点t=0に関係なく、相対時間τの関数である軌跡x+ k(τ)は、常に原点τ=0を含む。
【0025】
上記構成によれば、レベルLの各上昇交差は、基本信号x(t)の対応する時間シフトされた複製を規定する。図3は、図示された信号セグメントx(t)におけるレベルLの全ての上昇交差により生成された軌跡を連続的にかつ分離して示している。全ての上昇交差は、これらが相対時間τの同一原点τ=0を共同して規定するとともに共有している点で一致している。
【0026】
図3の軌跡は、図4において、相対時間τの関数として重ね合わせても示されている。
【0027】
{t+ i,i=1,2,・・・,k−1,k,k+1,・・・}における対応する上昇交差と関連付けられた軌跡{x+ i(τ),i=1,2,・・・,k−1,k,k+1,・・・}は、本明細書では前方クロスレーション(FC)関数と称される関数C+ L(τ)を導くために平均化することができる。説明上、図5に、図4に示される軌跡を平均することで得られる実験的前方クロスレーション関数C+ L(τ)を示す。この関数は、レベルLの上昇交差に調整された信号x(t)の平均的挙動を特徴付け、相対時間τによって決まる。特に、τ=0における値は、構成上、図4から推論されるように、単純にLに等しい。τの大きい値に対しては、上昇交差に対する依存性が消滅するため、C+ L(τ)は、基本一次プロセスx(t)の平均AVに近づく。
【0028】
同様に、信号x(t)がレベルLと負の傾斜をもって交差する時間的瞬間が決定される。図2b)に示される結果としての時間的瞬間
t-1,t- 2,・・・,t- m-1,t- m,t- m+1,・・・
は、レベルLの一組の下降交差を形成する。
【0029】
図3ないし図5を参照して説明したプロセスと類似のプロセスにより、上昇交差ではなく下降交差に基づく点を除いて前方クロスレーション関数C+ L(τ)に対応する、図6に示す関数C- L(τ)を導くことが可能である。したがって、この関数は、レベルLの下降交差に調整されたx(t)の平均的挙動を表す。
【0030】
信号x(t)によるレベルLの下降交差は、基本信号x(t)の時間反転されたレプリカx(−t)によるレベルLの上昇交差と一致することに留意すべきである。したがって、下降交差に基づくクロスレーション関数C- L(τ)は、後方クロスレーション(BC)関数と称されることになる。この場合も、C- L(0)=Lであり、C- L(|τ|)は大きい値のτに対して平均値AVに近づく。
【0031】
正の値のみをとる単極信号について前方および後方クロスレーション関数が決定される場合には、しきい値レベルLは常に正である。しかし、2極信号の場合は、複数の手法が可能である。
1.非負(または非正)のしきい値レベルのみが使用される。
2.正負(ゼロを含む)のしきい値レベルが信号処理に使用できる。
3.非負(または非正)のしきい値レベルのみが使用されるが、元の信号とその極性が反転されたレプリカの両方が処理される。
【0032】
前方クロスレーション(FC)関数および後方クロスレーション(BC)関数は、調査対象のプロセスの有用な特性化をもたらす。例えば、相対時間τの正の値に対しては、プロセスが正の傾斜をもって所定のレベルをいずれかの時間的瞬間に交差した場合には、前方クロスレーション(FC)関数は、そのプロセスの将来値の予測を容易にする。τの負の値に対しては、前方クロスレーション(FC)関数は、上昇交差の時間的瞬間の前のそのプロセスの平均的挙動を記述する。
【0033】
同様に、プロセスが負の傾斜をもって所定のレベルに交差した場合には、後方クロスレーション(BC)関数は、そのプロセスの将来値の予測を容易にする。相対時間τの負の値に対しては、後方クロスレーション(BC)関数は、下降交差の時間的瞬間の前のそのプロセスの平均的挙動を記述する。
【0034】
反転時間においてプロセスを調査した場合、前方クロスレーション(FC)関数と後方クロスレーション(BC)関数の役割は入れ替わる。したがって、時間反転可能なプロセスにおいては、前方クロスレーション(FC)関数と後方クロスレーション(BC)関数は、相互に鏡像である。このように、前方クロスレーション(FC)関数および後方クロスレーション(BC)関数は、対象のプロセスの時間反転性の評価に利用できる。
【0035】
本発明の一実施形態によれば、前方クロスレーション(FC)関数および/または後方クロスレーション(BC)関数は、図7に示すクロスレータ(CRS)を使用して導出することができる。図7のモニタ(MON)の部分を形成するクロスレータ(CRS)、および、以降説明される改良型クロスレータは、多くの異なるアプリケーションで使用できる、別個の集積回路上に形成することも考えられる、汎用装置として形成することができる点が注記される。これらのクロスレータにより提供される機能のいくつかは、ある種のアプリケーションでは必要でない場合があり、もちろん、以下で説明される全ての機能が図7のモニタ(MON)における使用に必要な訳ではない。
【0036】
クロスレータ(CRS)は、極性反転回路(PRC)、複数のタップを有するアナログ遅延線(TDL)、レベル交差検出器(LCD)、2つのパルス遅延回路(PDLおよびDEL)、パルスカウンタ(PCT)、複数のサンプル・ホールド回路(SHC)、複数のアキュムレータ(ACC)、および、ストレージレジスタ(SRG)を備える。ストレージレジスタ(SRG)は、適当な波形補間器を有しても良い。
【0037】
時間的に変動する入力信号x(t)の極性(正または負)は、極性反転回路(PRC)の2値(binary)極性選択入力(PS)において保持される適切な値により設定される。選択された極性の信号は、その後、遅延線(TDL)の入力(IP)に印加される。図示の構成では、遅延線(TDL)のM個のタップのそれぞれは、入力(IP)において生じる信号の時間遅延されたレプリカを供給する。任意の時間的瞬間において、遅延線(TDL)のM個のタップで観測される信号サンプルは、共同して、遅延線(TDL)に沿って伝播する信号の有限部分の離散時間表現を形成する。好ましくは、遅延線(TDL)の連続するタップ間の相対遅延は、一定値を有する。
【0038】
遅延線(TDL)のM個のタップのそれぞれは、サンプル・ホールド回路(SHC)のそれぞれに接続され、選択されたタップ(CT)、好ましくは中央のタップもまた、レベル交差検出器(LCD)に接続される。
【0039】
レベル交差検出器(LCD)は、2値選択器入力(UD)に保持される値に応じて、上昇交差か下降交差のいずれかを検出する。所望の交差レベルLは、レベル交差検出器(LCD)のしきい値入力(LV)に適切なしきい値を与えることにより設定される。前方クロスレーション(FC)関数を決定する場合には、レベル交差検出器(LCD)は上昇交差の検出器として動作する。同様に、後方クロスレーション(BC)関数を決定する場合には、レベル交差検出器(LCD)は下降交差を検出する。
【0040】
前方クロスレーション(FC)関数を決定する場合、レベル交差検出器(LCD)により中央タップ(CT)において指定のレベルLの上昇交差が検出される毎に、レベル交差検出器(LCD)の出力において短いトリガパルス(TP)が生成される。このトリガパルス(TP)は、共通のトリガパルス(TP)入力を介して、全てのサンプル・ホールド回路(SHC)の同時動作を開始させる。各サンプル・ホールド回路(SHC)は、その入力において生じる信号の瞬間値を捉えて、この値をそれぞれのアキュムレータ(ACC)に供給する。
【0041】
トリガパルス(TP)はまた、パルスカウンタ(PCT)の現在状態を1だけ増加させる。パルスカウンタ(PCT)の容量は、所定のレベル交差の数N(すなわち、処理される信号軌跡の数N)と等しい。トリガパルス(TP)はまた、遅延が好ましくはサンプル・ホールド回路(SHC)の安定化時間に等しい適当なパルス遅延回路(PDL)に与えられる。
【0042】
パルス遅延回路PDLから得られた遅延されたトリガパルスは、共通のアキュムレータ入力(DT)を介して、それぞれのサンプル・ホールド回路(SHC)により駆動される全てのアキュムレータ(ACC)の同時動作を開始させる。各アキュムレータ(ACC)の機能は、クロスレータ(CRS)の一つの全動作サイクルの間に、その入力に連続して生じる全てのN個のサンプルの加算または平均を実行することである。
【0043】
所定数Nのレベル交差がレベル交差検出器(LCD)により検出され、パルスカウンタ(PCT)により記録されたとき、サイクル終了(EC)パルスがパルスカウンタ(PCT)の出力に生成される。サイクル終了(EC)パルスは、パルスカウンタ(PCT)を、そのリセット入力(RT)を介してリセットし、これはまた、アキュムレータの内容のストレージレジスタ(SRG)への伝送を開始させる。パルス遅延回路(DEL)により適切に遅延されたそれぞれのサイクル終了(EC)パルスは、全てのアキュムレータ(ACC)を共通の入力リセット(RS)を介してその初期ゼロ状態に設定する。サイクル終了(EC)パルスの発生後すぐに、決定された前方クロスレーション(FC)関数を時間離散化したものがストレージレジスタ(SRG)の出力(CFO)において利用可能となる。
【0044】
ストレージレジスタ(SRG)において波形補間が行われない場合には、決定された前方クロスレーション(FC)関数はM個の値で表示される。しかし、ストレージレジスタ(SRG)において何らかの追加的な信号処理を実行して、アキュムレータ(ACC)により供給されるM個より多い一次値を有する前方クロスレーション(FC)関数の補間された(平滑化された)表示を生成することができる。
【0045】
図8は、σを処理信号の実効値として、3つの異なる上昇交差レベルLの値:L=σ、L=2σ、および、L=3σに対して実験的に定めた実験的前方クロスレーション(FC)関数の形状を示す。この場合では、クロスレータにより処理される信号x(t)は、物理的雑音源により生成されたものであった。
【0046】
後方クロスレーション(BC)関数を決定する場合、レベル交差検出器(LCD)によりタップ(CT)においてレベルLの下降交差が検出される毎に、レベル交差検出器(LCD)の出力において短いトリガパルス(TP)が生成される。残りの機能および動作は、前方クロスレーション(FC)関数を決定する場合にクロスレータが行う機能、動作と同一である。
【0047】
高速で変化する信号を処理する場合、レベル交差検出器(LCD)により導入される遅延は過大であるかもしれず、補償することが必要である。遅延補償は、例えば、次の2つの手法の1つを採用して達成することができる。
1.レベル交差検出器(LCD)が中央タップ(CT)の手前のタップにより駆動され、このようにして得られた前トリガパルスが補助回路によりレベル交差検出器(LCD)の出力で追加的に遅延され、これにより、(レベル交差検出器(LCD)およびこの回路により)導入される全遅延がこれら2つのタップ間の相対遅延と一致する。
2.遅延時間(TDL)、中央タップ(CT)の手前の前トリガタップ、および、レベル交差検出器(LCD)の遅延に一致する2タップ間の相対遅延により、専用の前トリガタップを設ける。
【0048】
以上では、入力信号x(t)は単極性であるという仮定の下で、動作の説明を行った。しかし、クロスレータ(CSR)は、2極性の信号を処理し、正負のしきい値交差双方に基づくそれぞれの関数を導出するために動作することも可能である。これを達成するために、負のしきい値に基づく関数が導出される毎に、レベル交差検出器(LCD)が要求される関数を導出するために対応する正の交差レベルを使用できるように、極性選択入力(PS)における信号により極性反転回路(PRC)が入力信号x(t)の極性を反転させる。
【0049】
ここで、図7のモニタ(MON)の動作を説明する。
【0050】
当初は、パラレル−シリアル変換器(PTS)は、レベル交差検出器(LCD)のしきい値入力(LV)に量子化レベルの値1を送るように構成される。レベル交差検出器の信号入力(UD)は、クロスレータがその出力(CFO)に前方クロスレーション関数を生成するように設定される。
【0051】
図5を参照して、クロスレーション関数の値と入力関数x(t)の平均値AV間の差の絶対値がしきい値THより大きい場合は、クロスレーション関数は有意値を持つものと仮定する。したがって、その値は、−τaから+τbまでの範囲で有意である。
【0052】
サンプリング間隔が|τb|よりも小さい場合は、連続するランダム値がそれに先行する値に依存する偏りを有する危険がある。これは、τの正の値に対する有意な前方クロスレーション関数のレベルはこの関数の前方予測性を表すからである。これに対応して、サンプリングレベルが|τa|よりも小さければ、その場合は、先行の乱数がそれに後続する値に関連する偏りを有する、すなわち、後方予測性、すなわち、後続の値から先行の値を決定できることの危険がある。乱数生成においては、乱数の「シード値」の予測を防ぐためにこれを避けることが重要である場合がある。
【0053】
したがって、最小サンプリング間隔が|τa|と|τb|の最大値よりも大きいことを確実にすることが望ましい。時間シフト比較器(TSC)は、クロスレーション関数を評価して、クロスレーション関数と入力信号x(t)の平均値AVとの間に有意な差が存在する|τ|の最大値を決定する。
【0054】
その後、入力(UD)は、クロスレータが後方クロスレーション関数をその出力において生成するように切り替えられ、時間シフト比較器は、再度、クロスレータの出力が有意となる最大値|τ|を見つけるよう動作する。
【0055】
そして、パラレル−シリアル変換器(PTS)は、レベル交差検出器(LCD)に第2の量子化レベルを伝送するよう動作し、前方および後方クロスレーション関数を得るためにクロスレータの動作が繰り返される。このシーケンスは、各量子化レベル1〜Jに対して実行される。
【0056】
したがって、時間シフト比較器(TSC)は、i=0(前方クロスレーションの場合)または1(後方クロスレーションの場合)で、J=1〜Jとして、前方および後方クロスレーション関数の双方および全ての量子化レベル1〜Jに対する複数の値τijを計算する。ここで、各値τijは、それぞれのクロスレータ関数が平均値AVから有意に異なっている最大値|τ|を表している。
【0057】
最小サンプリング間隔MSIは、ここで、
i=0,1、および、J=1〜Jに対して、MIS=τijの最大値
として計算される。
【0058】
この動作が、より詳細に図9のフローチャートに示されている。ステップ900で第1の量子化レベル(j=1)が選択され、ステップ902で前方クロスレーション(i=0)が選択される。ブロック904中で示される手順は、値τijを導出することを意図したものである。ステップ906において、(後方クロスレーションを選択するために)iが増分され、ステップ908において、iが1を超えてないかどうか確かめるためにチェックされる。超えていない場合は、904の手順が、後方クロスレーションに対する値τijを導出するために反復される。
【0059】
値iは、ステップ906において再度増分され、今回は、ステップ908がiが1を越えたことを検出し、そのため、プログラムはステップ910に進む。ここで、値jは次の量子化レベルを選択するために増分される。ステップ912において、プログラムは最後の量子化値Jを超過していないことを判断し、したがって、902から910のステップが反復される。このように、904の手順の間に、jの全ての値、および、前方および後方クロスレーション関数の双方に対して値τijが計算される。
【0060】
手順904は、ステップ914において、起こり得る最大のτの値、τmaxに等しい変数τHを初期的に設定することを含む。
【0061】
ステップ916において、プログラムは、この時点τHにおけるクロスレーション関数の値、すなわち、V(τH)から入力信号x(t)の平均値AVを引いた間の差を決定する。プログラムはその後、この差の絶対値が所定のしきい値THより大きいか否かを判断する。プログラムは、τの最高値、τmaxを調査することからスタートするため、クロスレーション関数は、平均レベルAVにほぼ等しくなり、したがって、プログラムはステップ918に進むことになる。この時点で、τHの値は(遅延線(DTL)の連続する段間の遅延を表す)ある増分量τiだけ減少させられる。ステップ916は反復される。
【0062】
このように、プログラムは、最高値τmaxからスタートして、クロスレーション関数がしきい値THから外れることをステップ916が検出するまで、クロスレーション関数を調査する。この時点で、プログラムはステップ920に進む。
【0063】
ステップ920において、プログラムは、他の変数τLを、起こり得る最小のτの値、tminに等しいように設定する。プログラムは、その後、ステップ922に進む。ここで、プログラムは、現在の値τLに対する相関関数と平均値AVとの差がしきい値THを超えているかどうかを判断する。超えていない場合は、プログラムはステップ924に進み、τLは増分値τiだけ増大させられる。そして、プログラムはステップ922に戻る。プログラムが増加する値τに対するクロスレーション関数を引き続きチェックしてこの値がしきい値領域外となるまで、これが継続する。その後、プログラムは、ステップ926に進む。
【0064】
ステップ926において、プログラムは、値τijをτHとτLの最大値に等しいように設定して、この値τijを後の使用の為に記録する。
【0065】
図9に示す手順の最後に、プログラムは、ステップ912からステップ928に進み、ここで、最小サンプリング間隔MISが、記録された全てのτij値の最大値に等しいように設定される。
【0066】
この値は、この値と実際のサンプリング間隔を表す値SPIとの比較を行う比較器(CMP)に送られる。実際のサンプリング間隔がMSIよりも大きい場合には、比較器の出力(MSP)は、連続する乱数が統計的に独立であると予測されることを示す。所望により、比較器の出力は、サンプリング間隔を制御するため、すなわち、現在のサンプリング間隔がMSIよりも小さいと判断された場合にはこれを増加させるために使用することができる。
【0067】
前方クロスレーション(FC)関数、および、後方クロスレーション(BC)関数は、調査対象のプロセスの有用な特性化を提供するものであるが、実際の応用では、前方クロスレーション(FC)関数と後方クロスレーション(BC)関数の合計や差などのある種の組み合わせがより有益であることが明らかとなる場合がある。
【0068】
前方クロスレーション(FC)関数C+ L(τ)と後方クロスレーション(BC)関数C- L(τ)の合計SL(τ)
SL(τ)=C+ L(τ)+C- L(τ)
は、クロスレーション合計(CS)関数と称され、典型例が図10に示されている。クロスレーション合計(CS)関数SL(τ)は、従来の自己相関関数により提供される情報と幾分類似した情報を提供する。特に、ガウシアンプロセスのクロスレーション合計関数は、このプロセスの自己相関関数に比例する。更に、任意の時間反転可能なプロセスのクロスレーション合計(CS)関数は、その引数である相対遅延τの偶関数である。
【0069】
前方クロスレーション(FC)関数C+ L(τ)と後方クロスレーション(BC)関数C- L(τ)の差DL(τ)
DL(τ)=C+ L(τ)−C- L(τ)
は、クロスレーション差(CD)関数と称される。典型例はやはり図10に示されている。クロスレーション差(CD)関数DL(τ)は、従来の自己相関関数の導関数により提供される情報に関連する情報を提供する。特にガウシアンプロセスのクロスレーション差(CD)関数は、当該プロセスの否定(negated)導関数に比例する。また、任意の時間反転可能なプロセスのクロスレーション差(CD)関数は、その引数である相対遅延τの奇関数である。
【0070】
クロスレーション合計(CS)関数、および、クロスレーション差(CD)関数は、図11に示す改良型クロスレータ(CRS)を使用して時間連続信号x(t)に対して決定できる。このシステムは、極性反転回路(PRC)、複数のタップを有するアナログ遅延線(TDL)、レベル交差プロセッサ(LCP)、2つのパルス遅延回路(PDL、および、DEL)、パルスカウンタ(PCT)、複数のサンプル・ホールド回路(SHC)、複数の加算/減算アキュムレータ(ASA)、および、ストレージレジスタ(SRG)を備える。ストレージレジスタ(SRG)はまた、適当な波形補間器を備えていてもよい。
【0071】
改良型クロスレータにより実行される動作は、図7の基本的なクロスレータ(CRS)により実行される動作と次の点で相違する。
【0072】
遅延線(TDL)の中央タップ(CT)においてレベル交差(上昇交差、または、下降交差)が検出される毎に、レベル交差プロセッサ(LCP)は、短いトリガパルス(TP)を生成する。所望の交差レベルLは、レベル交差プロセッサ(LCP)のしきい値入力(LV)に適切なしきい値を与えることにより設定される。クロスレーション合計関数、または、クロスレーション差関数を決定するために必要な動作モードは、レベル交差プロセッサ(LCP)の2値選択器入力(SD)に適切な値を与えることにより選択される。
【0073】
それぞれの加算/減算アキュムレータ(ASA)は、その制御入力(AS)に現れる「加算」または「減算」の指令に応じて、それぞれのサンプル・ホールド回路(SHC)により供給されるサンプル値を加算または減算する。
【0074】
クロスレーション合計(CS)関数が改良型クロスレータにより決定される場合、レベル交差プロセッサ(LCP)は、検出されたレベル交差のタイプ(上昇交差、または、下降交差)に関係なく、「加算」指令を共通の制御入力(AS)を介して全ての加算/減算アキュムレータ(ASA)に送る。しかし、クロスレーション差(CD)関数が決定される場合は、レベル交差プロセッサ(LCP)は,検出された各上昇交差に対しては「加算」指令を、検出された各下降交差に対しては、「減算」指令を送る。時間連続信号においては、(同一レベルの)上昇交差と下降交差が交番するため、加算および減算の動作も交差パターンに従って交番する。
【0075】
改良型クロスレータのシステムでは、パルスカウンタ(PCT)は全てのレベル交差をカウントするが、処理された上昇交差の数N+が処理された下降交差の数N-に正確に一致し、したがって、N+=N-=Nとなることを確実にするために、その容量は常に偶数2Nに設定される。
【0076】
図11のクロスレータ(CRS)は、例えば、それぞれの量子化レベルに対してクロスレーション合計のみを生成し、クロスレーション合計が入力信号x(t)の平均値から有意の差を示す最大遅延値|τ|を計算するために時間シフト比較器(TSC)を使用することにより、図7のモニタ(MON)において使用できる。
【0077】
図7の基本的なクロスレータ、または、図11の改良型クロスレータに採用される複数のタップを有するアナログ遅延線(TDL)は、アナログまたはデジタルのシリアル・イン−パラレル・アウト(SIPO)シフトレジスタで置き換えることができる。図12は、SIPOシフトレジスタ(SIPOSR)を組み込んだ図7の基本的なクロスレータのブロックダイアグラムである。このシステムもまた、信号調整回路(SCU)、クロック発生器(CKG)、レベル交差検出器(LCD)、2つのパルス遅延回路(PDL、および、DEL)、パルスカウンタ(PCT)、複数のサンプル・ホールド回路(SHC)、複数のアキュムレータ(ACC)、および、ストレージレジスタ(SRG)を備えている。このストレージレジスタ(SRG)も、適当な波形補間器を備えていても良い。
【0078】
アナログ時間連続信号x(t)は、信号調整回路(SCU)により適切な(アナログまたはデジタルの)形態に変換され、その後、SIPOSRのシリアル入力(IP)に印加される。
【0079】
SIPOシフトレジスタは、C1、C2、・・・、CMのM個の記憶セルからなっている。それぞれのセルは、入力端子、出力端子、および、クロック端子(CP)を有する。これらのセルは、最初のセル(C1)と最後のセル(CM)を除いて、各セルのその入力端子が先行のセルの出力端子に接続され、その出力端子が後続のセルの入力端子に接続されるように、直列に接続されている。セルC1の入力端子は、SIPOシフトレジスタのシリアル入力として使用される。全てのM個のセルの出力端子は、SIPOシフトレジスタのパラレル出力端子とみなされる。全てのセルのクロック端子(CP)は、一緒に接続されて、SIPOシフトレジスタのクロック端子を形成する。
【0080】
適切なクロックパルスのシーケンスがクロック発生器(CKG)から与えられる。時間的瞬間t0においてクロックパルスがSIPOシフトレジスタのクロック端子に加えられた場合、各セルに記憶された信号サンプルは、後続のセルに転送(シフト)されてこれに記憶され、セルC1は、入力信号x(t)の値x(t0)を記憶する。シフトレジスタは、デジタル素子、あるいは、例えば、「バケットブリゲード」型の電荷結合素子(CCD)の形態の時間離散型アナログ素子として実施することができる。
【0081】
SIPOシフトレジスタのパラレル出力は、それぞれのM個のサンプル・ホールド回路(SHC)に接続される。2つの選択された隣接するSIPOSR出力は、レベル交差検出器(LCD)の2つの入力にも接続されている。図12に示すシステムでは、その選択された出力は、セルCYとセルCZの出力である。
【0082】
SIPOSR出力の個数Mが奇数である場合には、好ましくは、上記2つの選択された出力の一方は、SIPOSRの中央の出力、すなわち、出力(M+1)/2である。しかしながら、SIPOSR出力の個数が偶数である場合には、上記2つの選択された出力は、好ましくは、出力M/2と出力M/2+1である。
【0083】
SIPOシフトレジスタは、クロック発生器(CKG)により供給されるクロックパルスで規定される離散時間において動作するため、信号サンプルによる所定レベルLの交差の検出は、幾分より複雑である。しかしながら、交差の検出は、次の判断ルールを適用することにより達成することができる。
A.CYの出力<Lであり、かつ、CZの出力>Lである場合、レベルの上昇交差が、セルCYとセルCZの間に位置する「仮想」セルVCにおいて生じている。
B.CYの出力>Lであり、かつ、CZの出力<Lである場合、レベルの下降交差が、セルCYとセルCZの間に位置するセルVCにおいて生じている。
C.それ以外であれば、セルVCにおいてレベルの交差は生じていない。
【0084】
統計的考察により、クロック発生器の周期が処理されている信号の時間的変動に比較して小さい場合には、仮想セルVCの「時間的」位置は、そのクロック周期に渡って均一に分散されていることになる。したがって、仮想セルVCは、セルCYとセルCZの中間に「位置」している。
【0085】
上述したクロスレータ(CRS)は、前方および後方クロスレーション関数(これらから、クロスレーション合計関数およびクロスレーション差関数が導出できる)を別個に生成し、または、クロスレーション合計関数およびクロスレーション差関数を直接に生成することを可能にする。これらの関数は、正負の双方であり得る、異なるそれぞれの交差レベルに対して生成することができる。特に便宜な構成においては、入力信号x(t)は、クロスレーション関数の処理の簡略化を可能とする、ゼロの平均値AVを有する。
【0086】
どのクロスレーション関数、あるいは、関数の組み合わせを使用すべきかの選択は、クロスレータの用途に依存する。前方および後方クロスレーション関数の両方を別個に生成することは、信号の予測性の判断に有用であることが予見される。しかし、信号の分類などの他の状況は、クロスレーション合計関数および/またはクロスレーション差関数の使用を正当化するかもしれない。いずれにせよ、これらの関数は、単一の交差レベル、あるいは、複数の交差レベルに対して導出できる。一般的には、非ガウシアン信号に対しては、信号x(t)の平均AVから有意に異なる1つまたは複数の交差レベルを使用することがより有益である。
【0087】
他のタイプのクロスレーション関数を導出することも可能である。上述した構成では、各関数は、それぞれの交差レベルに対応している。それぞれ異なる交差レベルに関連するクロスレーション関数の組み合わせ(例えば、それらの差)に関連する追加的な関数を導出することが可能である。例えば、平均値AVの交差レベルに基づくクロスレーション関数(すなわち、前方、または、後方クロスレーション関数のいずれか)を正のレベルLに対する対応するクロスレーション関数から差し引くことができる。これの結果関数は、ガウシアン信号に対しては、自己相関関数の拡縮されたレプリカである。結果物を別個に導出した自己相関関数と比較することにより、入力信号の特性が、ガウシアンの特性から外れている程度を判断することが可能である。更に、クロスレーションの技法をガウシアンの信号に対する自己相関関数の導出に用いることもまた、独立して有用であると見なされる。
【0088】
上述した構成では、入力信号x(t)の振幅ではなく、その傾斜の符号のみが考慮された。しかし、これは必須ではなく、代わりに、クロスレータは、正負の各方向における異なる振幅の傾斜同士を区別するよう構成することができる。すなわち、傾斜は、(正か負かの傾斜のいずれかを表示する)単一ビットではなく、2つまたはそれ以上のビットにより表現することができる。この状況では、別個のクロスレーション関数をそれぞれの量子化された傾斜レベルに対して導くことができる。あるいは、(例えば、最も急峻な傾斜など)ある種の量子化された傾斜レベルのみをクロスレーション関数の導出において考慮するような構成であってもよい。
【0089】
入力信号x(t)は、例えば、雑音、圧力、変位、速度、温度など、興味のある任意の物理量を表すことができる。したがって、本発明は、例えば、通信、電波天文学、遠隔検知、水中音響学、地球物理学、会話分析、バイオ医薬など、広範な適用分野を有する。上で提示した特定の例は、時間とともに変動する入力信号に言及したものだが、当該関数の引数は、例えば、相対時間、距離、空間的場所、角度位置など、任意の独立変数を表すことができる。
【0090】
上で示したように、クロスレータ(CRS)が別個の集積回路素子で形成されている場合は、入力信号x(t)のための入力端子、交差レベルを示す信号(LV)を受信するためのしきい値端子、および、パラレル、あるいは、シリアルの形態のいずれかにおける出力機能(CSO)を提供する少なくとも一つの出力端子を備えることが好ましい。
【0091】
導出されたクロスレーション関数は、分類の目的に使用することができ、よって、例えば、導出されたクロスレーション加算波形などのクロスレーション波形は、その信号を発生させた対象を最もよく表す特定のクラスを示すために使用される。この目的のために、クロスレーション波形の一組の代表的な「テンプレート」(各テンプレートがそれぞれのクラスに対応し、そのクラスのクロスレーション関数の形状を表す)を記録するために適当なメモリーを備えることができる。分類は、決定されたクロスレーション関数の適当な表示と記録されたテンプレートとの最適適合を見つけることにより実行できる。
【0092】
クロスレーション波形の形状は、いくつかの(「未知」を含む)信号を発信する物体のクラスの区別に使用する「指紋」サインと考えることができる。
【0093】
上記の本発明の好ましい実施形態の記述は、説明と記述の目的で提示されたものである。網羅的とすることも、本発明を開示されたそのままの形態に限定することも意図するものではない。上述の説明によれば、多くの変更、改良、変形により、当業者が意図された特定の用途に適合した多様な態様で本発明を利用することが可能となることは明らかである。
【図面の簡単な説明】
【0094】
【図1】本発明による信号解析装置を組み込んだ乱数発生器を示す図である。
【図2】図1の発生器に使用される無秩序な信号x(t)を示す図である。
【図3】無秩序な信号x(t)のセグメント、および、この信号セグメント中に観測されるレベルの全ての上昇交差に関連する複数の軌跡を示す図である。
【図4】重ね合わされた場合の図3の軌跡を示す図である。
【図5】図4の軌跡を平均することで得られる無秩序な信号x(t)の実験的前方クロスレーション関数C+ L(τ)を示す図である。
【図6】無秩序な信号x(t)の実験的後方クロスレーション関数C- L(τ)を示す図である。
【図7】信号解析装置を組み込まれた図1の発生器のモニタリングユニットのブロックダイアグラムである。
【図8】σを調査対象の信号の実効値とした場合の3つの異なる交差レベルL:(a)L=3σ、(b)L=2σ、(c)L=σに対して実験的に得られた実験的前方クロスレーション関数C+ L(τ)の形状を示す図である。
【図9】図7のユニットの時間シフト比較器の動作を示すフローチャートである。
【図10】クロスレーション合計関数SL(τ)、および、クロスレーション差関数DL(τ)の形状を示す図である。
【図11】図7の信号解析装置の改良版のブロックダイアグラムである。
【図12】信号解析装置の他の改良版を示す図である。
Claims (11)
- 入力信号の統計的特性を解析する装置であって、
前記信号を受信する信号入力と、
前記入力に接続され、前記信号のレベルが所定のレベルと所定の傾斜をもって交差するイベントを検出する手段と、
前記イベントの間隔に対応する量だけ互いにシフトされた前記信号の複数のバージョンを結合して、前記信号の表示を形成する手段と、
前記表示の形状に依存し、前記信号の統計的特性を示すパラメータを測定する手段と
を備える装置。 - 前記傾斜が所定の符号を有する場合に、信号が所定の傾斜を有すると見なされるよう構成された
請求項1に記載の装置。 - 検出された第1の所定の傾斜のイベントに応じて第1の表示を形成し、検出された第2の異なる所定の傾斜のイベントに応じて第2の表示を形成するよう構成された
請求項1または2に記載の装置。 - 前記パラメータは、結合された前記第1と第2の表示の形状に依存する
請求項3に記載の装置。 - 前記表示を形成するために、前記イベントを検出する手段は、第1と第2の異なるイベントのタイプの検出を実行でき、前記結合する手段は、所定の態様で前記第1のイベントのタイプに対応する量だけシフトされた前記信号のバージョンを、前記第2のイベントのタイプの間隔に対応する量だけ相互にシフトされた前記信号のバージョンに結合することを実行できる
請求項1〜4のいずれか一項に記載の装置。 - 前記結合の所定の態様の変更を実行できるモード切替手段を有する
請求項5に記載の装置。 - 前記所定のレベルは、前記信号の平均レベルと実質的に相違する
請求項1〜6のいずれか一項に記載の装置。 - 前記所定のレベルを規定する信号を受信する交差レベル入力手段を備える
請求項1〜7のいずれか一項に記載の装置。 - 請求項1〜8のいずれか一項に記載の装置と、
前記入力信号を受信する第1の入力端子と、
前記所定のレベルを表すしきい値信号を受信する第2の入力端子と、
前記表示を形成する出力信号を供給する少なくとも一つの出力端子と
を備える集積回路。 - 入力信号を解析する方法であって、
前記信号のレベルが所定の傾斜をもって所定のレベルに交差するイベントを検出することと、
前記イベントの間隔に対応する量だけ相互にシフトされた前記信号の複数のバージョンの結合の表示を形成することと
を含み、
更に、前記表示の形状に依存するパラメータを測定するステップ
を含む方法。 - 前記パラメータは、前記形状と記録された表示の形状との間の類似度の度合いを示す請求項10に記載の方法。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP02254612A EP1378854B1 (en) | 2002-07-01 | 2002-07-01 | Signal statistics determination |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004163889A true JP2004163889A (ja) | 2004-06-10 |
JP4612286B2 JP4612286B2 (ja) | 2011-01-12 |
Family
ID=29719768
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003189764A Expired - Fee Related JP4612286B2 (ja) | 2002-07-01 | 2003-07-01 | 入力信号の統計的特性を解析する装置、並びに入力信号を解析する集積回路および方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US7120555B2 (ja) |
EP (1) | EP1378854B1 (ja) |
JP (1) | JP4612286B2 (ja) |
CN (1) | CN1330943C (ja) |
DE (1) | DE60224263T2 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010204102A (ja) * | 2009-03-03 | 2010-09-16 | Mitsubishi Electric R&D Centre Europe Bv | スペクトル解析および同スペクトル解析を利用したfmcw自動車レーダ |
JP2010210620A (ja) * | 2009-03-03 | 2010-09-24 | Mitsubishi Electric R&D Centre Europe Bv | スペクトル解析 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1306805A1 (en) * | 2001-10-25 | 2003-05-02 | Mitsubishi Electric Information Technology Centre Europe B.V. | Image Analysis |
EP1528407A1 (en) | 2003-10-31 | 2005-05-04 | Mitsubishi Electric Information Technology Centre Europe B.V. | Decomposition of a wideband random signal |
EP1596219A1 (en) | 2004-05-13 | 2005-11-16 | Mitsubishi Electric Information Technology Centre Europe B.V. | Signal processing circuit for time delay determination |
US7453765B2 (en) * | 2006-05-16 | 2008-11-18 | Ikelle Luc T | Scattering diagrams in seismic imaging |
JP2008182425A (ja) * | 2007-01-24 | 2008-08-07 | Denso Corp | フィルタ回路 |
EP2221732A1 (en) | 2009-02-23 | 2010-08-25 | Mitsubishi Electric R&D Centre Europe B.V. | Signal processing |
EP2287748A1 (en) * | 2009-08-21 | 2011-02-23 | Mitsubishi Electric R&D Centre Europe B.V. | Determination of system characteristics |
KR20130116857A (ko) * | 2010-06-25 | 2013-10-24 | 가부시키가이샤 한도오따이 에네루기 켄큐쇼 | 액정 표시 장치 및 전자 기기 |
CN112422360A (zh) * | 2020-10-14 | 2021-02-26 | 锐捷网络股份有限公司 | 一种报文采样方法、装置、设备及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1185472A (ja) * | 1997-09-12 | 1999-03-30 | Toshiba Corp | 物理乱数発生装置及び方法並びに物理乱数記録媒体 |
JP2000066592A (ja) * | 1998-08-19 | 2000-03-03 | Syst Kogaku Kk | 乱数生成装置 |
JP2000089672A (ja) * | 1998-09-09 | 2000-03-31 | Ntt Data Corp | 暗号化回路及びランダムパスワード発生回路 |
JP2000276330A (ja) * | 1999-03-26 | 2000-10-06 | Syst Kogaku Kk | 故障判断機能を備えた乱数生成装置 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3737788A (en) * | 1965-06-11 | 1973-06-05 | North American Rockwell | Slope responsive signal identification means |
US3732405A (en) * | 1971-08-11 | 1973-05-08 | Nasa | Apparatus for statistical time-series analysis of electrical signals |
US3896395A (en) * | 1974-07-18 | 1975-07-22 | Bell Telephone Labor Inc | Linear amplification using quantized envelope components to phase reverse modulate quadrature reference signals |
US4122393A (en) * | 1977-05-11 | 1978-10-24 | Ncr Corporation | Spread spectrum detector |
US4387465A (en) * | 1981-04-13 | 1983-06-07 | Trw Inc. | Sequential threshold detector |
JPS60105972A (ja) * | 1983-11-15 | 1985-06-11 | Iwatsu Electric Co Ltd | アナログ信号測定装置 |
US4897855A (en) * | 1987-12-01 | 1990-01-30 | General Electric Company | DPCM system with adaptive quantizer having unchanging bin number ensemble |
US4843309A (en) * | 1988-03-21 | 1989-06-27 | Tektronix, Inc. | Waveform timing alignment system for digital oscilloscopes |
US5272725A (en) * | 1991-02-25 | 1993-12-21 | Alliedsignal Inc. | Digital video quantizer |
US6249556B1 (en) * | 1998-05-27 | 2001-06-19 | Intel Corporation | Dynamic thresholding for input receivers |
US6205168B1 (en) * | 1998-11-12 | 2001-03-20 | Sharp Laboratories Of America, Inc. | Sequential detection system and method with adaptive bias |
GB9828693D0 (en) * | 1998-12-24 | 1999-02-17 | Mitsubishi Electric Inf Tech | Time delay determination |
US6571186B1 (en) * | 1999-09-14 | 2003-05-27 | Textronix, Inc. | Method of waveform time stamping for minimizing digitization artifacts in time interval distribution measurements |
US20030035549A1 (en) * | 1999-11-29 | 2003-02-20 | Bizjak Karl M. | Signal processing system and method |
SE0102475D0 (sv) * | 2001-07-10 | 2001-07-10 | Totalfoersvarets Forskningsins | Method of exploiting nonlinear dynamics in a fluxgate magnetometer and a low-power fluxgate magnetometer |
JP4345331B2 (ja) * | 2003-03-17 | 2009-10-14 | セイコーエプソン株式会社 | 遮光手段を用いた露光装置及び露光方法 |
-
2002
- 2002-07-01 DE DE60224263T patent/DE60224263T2/de not_active Expired - Lifetime
- 2002-07-01 EP EP02254612A patent/EP1378854B1/en not_active Expired - Lifetime
-
2003
- 2003-07-01 JP JP2003189764A patent/JP4612286B2/ja not_active Expired - Fee Related
- 2003-07-01 CN CNB031545033A patent/CN1330943C/zh not_active Expired - Fee Related
- 2003-07-01 US US10/609,375 patent/US7120555B2/en not_active Expired - Lifetime
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1185472A (ja) * | 1997-09-12 | 1999-03-30 | Toshiba Corp | 物理乱数発生装置及び方法並びに物理乱数記録媒体 |
JP2000066592A (ja) * | 1998-08-19 | 2000-03-03 | Syst Kogaku Kk | 乱数生成装置 |
JP2000089672A (ja) * | 1998-09-09 | 2000-03-31 | Ntt Data Corp | 暗号化回路及びランダムパスワード発生回路 |
JP2000276330A (ja) * | 1999-03-26 | 2000-10-06 | Syst Kogaku Kk | 故障判断機能を備えた乱数生成装置 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010204102A (ja) * | 2009-03-03 | 2010-09-16 | Mitsubishi Electric R&D Centre Europe Bv | スペクトル解析および同スペクトル解析を利用したfmcw自動車レーダ |
JP2010210620A (ja) * | 2009-03-03 | 2010-09-24 | Mitsubishi Electric R&D Centre Europe Bv | スペクトル解析 |
Also Published As
Publication number | Publication date |
---|---|
JP4612286B2 (ja) | 2011-01-12 |
DE60224263T2 (de) | 2008-12-11 |
EP1378854B1 (en) | 2007-12-26 |
EP1378854A1 (en) | 2004-01-07 |
US7120555B2 (en) | 2006-10-10 |
DE60224263D1 (de) | 2008-02-07 |
CN1330943C (zh) | 2007-08-08 |
CN1495414A (zh) | 2004-05-12 |
US20040059517A1 (en) | 2004-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9506951B2 (en) | Method and apparatus for data acquisition with waveform trigger | |
JP4612286B2 (ja) | 入力信号の統計的特性を解析する装置、並びに入力信号を解析する集積回路および方法 | |
JP2002533732A (ja) | 時間遅延の決定および信号シフトの決定 | |
JPH07154378A (ja) | 光伝送特性測定装置 | |
TWI667486B (zh) | Partial discharge discriminating device and partial discharge discriminating method | |
JP2005354676A (ja) | 信号処理回路及び方法ないしこの方法を用いた時間遅延検出装置及び物体位置特定装置 | |
JP2014077791A (ja) | 複数信号の相互関係を示す方法及び試験測定装置 | |
JPH10260209A (ja) | パルス信号分類装置 | |
JP4507245B2 (ja) | 時間測定システム、物体検出システム、シフト測定方法 | |
CN107643434B (zh) | 一种基于分段切比雪夫距离的复杂波形触发方法 | |
AU2019229361B2 (en) | Input count rate estimation in radiation pulse detectors | |
TWI573401B (zh) | 估測抖動容忍度的時脈資料回復電路與方法 | |
JP5828617B2 (ja) | スペクトル解析 | |
TWI429225B (zh) | 目標訊號判斷方法及其裝置 | |
JP2001013180A (ja) | 信号分析装置 | |
RU2169378C1 (ru) | Способ обнаружения сигнала и устройство для его осуществления | |
KR101712893B1 (ko) | 펄스 신호의 펄스 정보 추정을 위한 신호 감지 장치 및 그 제어 방법 | |
US7031412B2 (en) | Digital signal processing apparatus | |
RU2663881C1 (ru) | Устройство определения времени прихода оптического сигнала | |
RU2292114C2 (ru) | Способ управления порогом в оптимальном обнаружителе | |
JP2603758B2 (ja) | キャビテーション雑音検出方法及びその装置 | |
JP3479154B2 (ja) | 信号弁別回路 | |
JPH08304556A (ja) | 磁探信号処理器 | |
RU2100822C1 (ru) | Ранговый адаптивный последовательный обнаружитель сигналов | |
RU2342781C2 (ru) | Устройство для компенсации амплитудно-частотных искажений в линейных трактах с использованием двухмерных метрических пространств |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060206 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090623 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20090924 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20090929 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20091023 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100608 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100901 |
|
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: 20100921 |
|
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: 20101015 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131022 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4612286 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
S631 | Written request for registration of reclamation of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313631 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131022 Year of fee payment: 3 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
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 |
|
LAPS | Cancellation because of no payment of annual fees |