ウェーブレットベースの相関測定は、追加の計算オーバーヘッドの価値がありますか?


9

信号間の相関の尺度として、相関とコヒーレンスの両方を使用しました。時間-周波数アプローチが私にこれらの世界のベストを与えると思っていました。

私の質問は、この余分なデータが信号の全体像に十分に追加され、計算の一部としてウェーブレット変換を行うことに関連する計算コストの増加を正当化できるかどうかです。

参照:ArXiv論文:S.Klimenko、G.Mitselmakher、A.Sazonovによる「確率的重力波の検出のためのウェーブレットドメインでの相互相関手法」


追加の計算コストは​​いくらですか?FFTまたはFWTを使用して高速化できますか?
エンドリス

@endolithこれらのアルゴリズムをすでに組み込んでいると仮定すると、私は思います。
jonsca

1
まあ、コヒーレンスと相関の両方がFFTを使用する可能性がありますが、これはO(N log N)ですが、FWTはO(N)なので、ウェーブレット法は実際には高速かもしれません。:私は二度尋ねるにもかかわらず、しかし、これを明確に理解していないmath.stackexchange.com/questions/28581/... stackoverflow.com/questions/1787536/...
endolith

1
とにかく、あなたがしようとしていることに最も適切な方を使用する必要があります。これは、「どちらが良いですか?ドライバーかハンマーか?」と尋ねるようなものです。
endolith

1
@jonscaあなたの直感は実際に正しいです。どうやらDWT変換は時変であり、このプロパティは何らかの悪用につながる可能性があります。自分が取り組んでいるプロジェクトでも、まったく同じことをしています。目的は2つの信号間のTDOA(到達時間遅延)を推定することです。そのため、最初に(手書きの)DWTを使用してそれらを変換し、次にそれらを相互相関させます。これは、私の公開ドロップボックスから読むことができる論文へのリンクです。(dl.dropbox.com/u/4724281/waveletBasedTDOA.pdf
スペイシー

回答:


5

まず、ジョブに適したツールを使用する必要があります。 相関コヒーレンス対ウェーブレットベースの相関はすべて異なるものなので、この質問は「どちらが良いですか?ドライバーかハンマーか?」と尋ねるようなものです。それは、何をしようとしているのか、時間、周波数スペクトル、またはその両方の類似性を気にするかどうかによって異なります。

Oログ O

経験的にn個の実際の入力からn 個の出力を生成すると、nが約4096より大きい場合、PyWaveletのマルチレベルウェーブレット変換はNumPyのFFTよりも高速になります。

ここに画像の説明を入力してください

しかしながら

  1. これはPythonであり、2つの実装の効率は非常に異なる可能性があります。wavedec()高速ウェーブレット変換と見なされるかどうかさえわかりません。彼らは彼らのドキュメンテーションで略語DWTを使います。Haar DWTとFWTは同じものですか?
  2. 時間は使用するウェーブレットによって異なります。Meyerウェーブレットは、Daubechiesが同じ量のデータを生成するのに6倍の時間がかかります。
  3. FWTがどのように時間-周波数平面をタイリングする、またはn個の出力を生成することで、FFTを使用したn点の円形相互相関と同じ種類の類似性測定値を取得するのに十分かどうかは、まだわかりません。(技術的には、時間-周波数ではなく時間スケール平面ですが、複雑なモーレットウェーブレットでも同じだと思いますか?)FWTは平面の「クリティカルサンプリング」であり、FFTと同じ量のデータを生成します。したがって、それらを比較するのは公平なようです。

重要な点は、計算時間は少なくともどちらでもほぼ同じであるため、どちらを使用するかを決定するときに心配する必要はないと思います。


3

これは非常に遅いですが、とにかくそれは価値があります...

バツtバツΔstΔtΔsΔtバツtバツtΔteΔωtΔωバツt

ON

したがって、DWTを使用して時間スケール平面を調べても、それほど遠くまでは行きません。これは、DWTによって「訪問された」スケールが2の係数で分離され、FFTを使用して時間周波数平面で取得できるカバレッジよりもはるかに密度が低いため、特に当てはまります。他の多くの名前の中でも、並進不変のウェーブレット変換(非デシメートウェーブレット変換と呼ばれることもあります)を使用する必要があります。それでも、計算されたスケールサンプルのスパース性が残っています。

さらに、タイムスケール平面内の位置をエネルギー密度があると考えることが望ましい場合がよくあります。このアプローチは、前述の複雑なMorletウェーブレットなどの分析ウェーブレットを使用することで容易になります。変換不変性と分析性を計算時間に対してバランスさせる1つの方法は、複雑なデュアルツリーウェーブレット変換です。時間-周波数平面で同じことをするのはおそらくもっと簡単です。最初にFFTを実行し、すべての負の周波数をゼロにしてからIFFTを実行することにより、信号に対して近似ヒルベルト変換を実行します。

相関が時間の類似性を探し、コヒーレンスが周波数の類似性を探すという直感が正しい場合は、時間-周波数平面に固執する方がよいでしょう。計算は確かに簡単で、周波数軸に沿ってサンプリングを調整するのも簡単です。上記のどのアプローチも、スケール軸をより高密度にサンプリングすることを扱っていません。それを行うには、ほとんどの場合、継続的なウェーブレット変換に移動する必要がありますが、他にも気づいていないものがあるかもしれません。Matlabをお持ちの場合は、上記のリンクにアクセスしてください。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.