0

まず、私の数学は限られているので、この質問には簡単な答えがあるかもしれません。したがって、私は次の方程式を使用してグアシアン分布を作成しています。

height * np.exp( - ((x-mean)/width)**2 )

上記の式でグシアンを作成する場合、ピークの幅はどこに適用されますか？半値全幅ですか？次の値で次のガウス分布を作成しました。

height = 5
mean = 100
width = 10

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

次にFWHMを計算すると16.6510941453、ピーク幅をFWHMに適用できないようになります。どこに適用されますか？

FWHMが平均の10分の1になるように、FWHMを制限しようとしています。したがって、上記の例では、ガウス分布がピーク高さ5で平均100でFWHMが10になるようにしたいと思います。

4

1 に答える 1

2

あなたの方程式では、幅パラメータは実際sigmaには です。これはガウス分布の標準偏差であり、FWHM ではありません。以下は、これらの 2 つのプロパティ間で変換する関数です。

from numpy import sqrt, log

def sigma2Gamma(sigma):
    '''Function to convert standard deviation (sigma) to FWHM (Gamma)'''
    return sigma * sqrt(2 * log(2)) * 2 / sqrt(2)

Gamma = sigma2Gamma(10)
print Gamma
# prints 16.651092223153956, which is what you saw in your graph

def Gamma2sigma(Gamma):
    '''Function to convert FWHM (Gamma) to standard deviation (sigma)'''
    return Gamma * sqrt(2) / ( sqrt(2 * log(2)) * 2 )

sigma = Gamma2sigma(10)
print sigma
# prints 6.0056120439322491, which is the standard deviation that will
# give a FWHM of 10

方程式を次のように変更することをお勧めします

height * np.exp( - ((x-mean)/Gamma2sigma(width))**2 )

標準偏差ではなくFWHMを入力したい場合

于 2012-05-16T18:02:36.647 に答える