所以我试图计算一些非常简单的东西,但是我的代码给了我错误的答案。我有一组抛物线数据(变量“ chi”)(其卡方函数)。数据的顶点(最小的chi值)位于y = 13.671。我试图找出该抛物线将在哪里截取y = 14.671,但到目前为止,我得到的是假想的根,考虑到y = 14.671线在两个不同的点截取了数据,这是没有意义的。
我已经尝试过'numpy.polynomial import Polynomial'包。基本上,我已经要求python给我抛物线的方程式,然后告诉我减去y值的根在哪里。
fit = P.fit(betatest,chi,2)
yval = min(chi) + 1.
rts = (fit - yval).roots()
sigma = abs(beta_bg - rts[0])
sig = "%.3f" % sigma
##note: P is from 'numpy.polynomial import Polynomial', and
## betatest and chi are both 1D vectors of length 1000. Betatest is
## just a linspace. Again, the minimum for chi is 13.671, and it's a
## positive parabola (smiley face).
我没有收到任何错误消息,但是我的想法是虚构的。当我运行代码并要求rts时,我得到了:
在[24]中:rts Out [24]:数组([0.37239476-0.15967081j,0.37239476 + 0.15967081j])
我拍了一张我的数据的图片以及要在其中找到xvalue的yvalue。 enter image description here
我做错什么了吗?
谢谢!