私は実験から、そして理論式から2組の周波数データを持っています。私はscipyの最小限の機能を使いたい。 ここに私のコードスニペットがあります。 ここで、gは私が探したいカップリングです。 Ad indはx軸のプロットのインダクタンスです。エラー機能で複数のパラメータを使用してscipy.minimizeを使用するにはどうすればよいですか?
from scipy.optimize import minimize
def eigenfreq1_func(ind,w_q,w_r,g):
return (w_q+w_r)+np.sqrt((w_q+w_r)**2.0-4*(w_q+w_r-g**2.0))/2
def eigenfreq2_func(ind,w_q,w_r,g):
return (w_q+w_r)-np.sqrt((w_q+w_r)**2.0-4*(w_q+w_r-g**2))/2.0
def err_func(y1,y1_fit,y2,y2_fit):
return np.sqrt((y1-y1_fit)**2+(y2-y2_fit)**2)
g_init=80e6
res1=eigenfreq1_func(ind,qubit_freq,readout_freq,g_init)
print res1
res2=eigenfreq2_func(ind,qubit_freq,readout_freq,g_init)
print res2
fit=minimize(err_func,args=[qubit_freq,res1,readout_freq,res2])
しかし、それは次のエラーを示しています:
:
初期パラメータ推定値の配列を忘れてしまった – Uvar