在python/numpy/scipy中查找两个数组之间的插值交集

按照SER的回答:

import numpy as np
x = np.array([0,1])
y1 = np.array([0,2])
y2 = np.array([1,0])

def solve(f,x):
    s = np.sign(f)
    z = np.where(s == 0)[0]
    if z:
        return z
    else:
        s = s[0:-1] + s[1:]
        z = np.where(s == 0)[0]
        return z

def interp(f,x,z):
    m = (f[z+1] - f[z]) / (x[z+1] - x[z])
    return x[z] - f[z]/m

f = y1-y2
z = solve(f,x)
ans = interp(f,x,z)

print(ans)

这个问题可以通过假设你找到一个零,然后对两个级数的差执行函数来简化。首先,“solve”查找符号转换发生的位置(意味着在中间的某个位置出现零),然后“interp”执行线性插值以查找解决方案。

在…上面 Digitizing an analog signal ,我创建了一个名为 find_transition_times . 您可以通过传递 y_1 - y_2 对于 y 和 0 对于 threshold :

In [5]: xSupport = np.array([0,1])
   ...: y_1 = np.array([0,2])
   ...: y_2 = np.array([1,0])
   ...: 

In [6]: find_transition_times(xSupport, y_1 - y_2, 0)
Out[6]: array([ 0.33333333])