ピーク検出は、pythonでオクターブまたはoct2pyを使用すると異なる結果を示します

Question

固定幅ガウスの数を1つの広いピークに合わせようとしています。私がこれを行うことができる唯一の方法はpeakfit.mです。サンプルデータはhereからダウンロードできます。

オクターブでは私が使用するコードは次のとおりです。

>>peakfit(data,90.5,3,3,11,0,0,0,0,[1,1,1],0,0) 
ans = 
    1.0000e+000 9.0012e+001 1.0185e+004 1.0000e+000 1.0749e+004 
    2.0000e+000 9.0500e+001 9.3075e+003 1.0000e+000 9.9035e+003 
    3.0000e+000 9.0988e+001 1.0186e+004 1.0000e+000 1.0749e+004 

これらは私が最初のデータセットを作成するために使用される3つのピークにかなり近いです。私は

import oct2py as op 
import numpy as np 

data=np.loadtxt(file) 

octave=op.Oct2Py() 
octave.peakfit(data,90.5,3,3,11,0,0,0,0,[1,1,1],0,0) 

#Out: 
array([[ 1, 91, 8873,  1, 14583], 
     [ 2, 88, 17314,  1, 400], 
     [ 3, 90, 11011,  1, 18459]]) 

をoct2py使用する場合

これらは、私は両方のケースで単一のピークに合わせた場合、結果は同じである結果です。

オクターブ

1.0000e+000 9.0500e+001 2.0576e+004 1.4670e+000 3.2130e+004 

Oct2py

array([[ 1.00000000e+00, 9.05000004e+01, 2.05763986e+04, 
      1.46695147e+00, 3.21304879e+04]]) 

違いから来るだろうか？

Answer 1

解決策が見つかりました。 oct2pyの引数を浮動小数点として渡す必要があります。

octave.peakfit(data,90.5,3.0,3.0,11.0,0.0,0.0,0.0,0.0,[1.0,1.0,1.0],0.0,0.0)