位置参数和功能出现问题

时间:2018-10-26 02:13:17

标签: python scipy

我正在使用scipy.odr分析三个数据集,而如何调试此特定问题我完全不知所措。当我运行它时,它告诉我函数*thin_lens1()*接受一个位置参数,但是给出了两个。我已经尝试了很多重新排列,但我仍然不知道我做错了什么,因为*thin_lens1()*的唯一论点是p1。有人会提供的任何帮助将不胜感激,因为我完全被困在这里。

import numpy as np
import matplotlib.pyplot as plt
import scipy
import scipy.odr, scipy.special, scipy.stats


#importing data & defining focus length constants

p1,q1,up1,uq1,p2,q2,up2,uq2,p3,q3,up3,uq3 = 
np.loadtxt('lenses_data_processed.txt', 
                                     comments='#', unpack=True, 
usecols = range(0,12))

f1 = 158.695


#defining thin lens equation

def thin_lens1(p1):
    return (p1-f1)/(p1*f1)

#parameter guesses
p_guess = [1, 2]

#applying ODR
fit = scipy.odr.Model(thin_lens1)
mydata= odr.RealData(p1, q1, sx=up1, sy=uq1)
myodr = odr.ODR(mydata, fit, beta0 = p_guess, maxit=500)
myoutput = myodr.run()
myoutput.pprint()

0 个答案:

没有答案