我试图使用这里引用的fsolve:http://glowingpython.blogspot.gr/2011/05/hot-to-find-intersection-of-two.html
找到两条曲线的交集。两条曲线基本上都是两个浮动数组。
第一个是一维数组Pmech ( Pmech(x) )
,第二个是二维数组Pair ( Pair(x,y) )
x轴对于两个阵列都是通用的,所以我要做的是让每一个y轴都能看到Pair和Pmech相交的地方。
我知道fsolve()
采用参数函数,而不是数组,因此我编写了两个基本函数来实现此功能:
def Pmix(x):
return Pmech[x]
def Paera(x,y):
return Pair[x,y]
如上面的链接所示,我实现了findIntersection
函数:
def findIntersection(fun1,fun2,x0):
return fsolve(lambda x: (fun1(x) - fun2(x,y) for y in range(1,100)),x0)
但我得到了以下错误:
TypeError: float() argument must be a string or a number
Traceback (most recent call last):
File "batteries.py", line 261, in <module>
findIntersection(Pmix,Paera,0)
File "batteries.py", line 238, in findIntersection
fsolve(lambda x: (fun1(x) - fun2(x,y) for y in range(1,100) ),x0)
File "/usr/lib/python2.7/dist-packages/scipy/optimize/minpack.py", line 125, in fsolve
maxfev, ml, mu, epsfcn, factor, diag)
minpack.error: Result from function call is not a proper array of floats.
在您的示例中,如果不使用
fsolve
,您可能会更容易地解决这个问题:(fun1(x) - fun2(x,y) for y in range(1,100))
是发电机
[fun1(x) - fun2(x,y) for y in range(1,100)]
是一个列表。你需要后者。
但是,正如btel在另一个答案中提到的,对于数组中的交集,不能只重用用于查找函数交集的代码。
相关问题 更多 >
编程相关推荐