如何为二次数据点集求根

时间:2019-05-31 08:38:10

标签: python polynomial-math quadratic poly

所以我试图计算一些非常简单的东西,但是我的代码给了我错误的答案。我有一组抛物线数据(变量“ 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

我做错什么了吗?

谢谢!

0 个答案:

没有答案