2017-08-26 12 views
1

NumPy配列に格納された数値的に計算された光パルスの時間的形状を持っています。私はガウスエンベロープ関数のパラメータ(特に、FWHM)を見つけたいと思いますが、最初はそれに合わせなければなりません。
エンベロープフィッティングのためのSciPy関数または関数セットが組み込まれているか、このコードを自分で記述する必要がありますか? (精度はあまり重要ではないので、私はおそらくちょうど地元maxesを取り、それらにガウスに合うと思います...)SciPy - ガウスエンベロープをフィッティングする

光パルスは、次のようになります。ちょうどあなたの場合は

The light pulse in question

答えて

1

最大値を見つけたいと思って、Scipyのargrelextremaを使うことをお勧めします。

しかし、封筒の解析フォームを見つけるために、ScipyのHilbert transform を使用すると思われます。次に、ガウス分布に合うようにエンベロープの配列を生成することができます。

ガウスフィットの場合、良い答えはhereです。

+0

ありがとうございます!これはヒルベルト変換との私の最初の遭遇なので、私は読んでいないだろう...それはかなりクールだ。 – Neinstein

関連する問題