Py学习  »  Python

python中线性插值的逻辑误差

BCKN • 5 年前 • 1352 次点击  

我的线性插值有一些逻辑错误,它适用于某些情况,但不能完全工作。

def interpolate(x, y, x_test):
    for i in range(len(x)):
        if x[i] > x_test:   #extrapolated condition: when the largest value of
            x_below = i - 1 #list x is greater than x_test 
            x_above = i 
            y_below = i - 1
            y_above = i
            break
        elif x[i] < x_test: #extrapolated condition: when the largest value of 
            x_below = i + 1 #list x is greater than x_test 
            x_above = i 
            y_below = i + 1
            y_above = i
            break                
        else:             #interpolated condition: when x_test lies between  
            return y[i]    #two sample points.

    #a = (yabove - ybelow) / (xabove - xbelow)         
    a = (y[y_above] - y[y_below]) / (x[x_above] - x[x_below])  
    #b = ybelow - a * xbelow
    b = y[y_below] - a * x[x_below]
    #y’ = a * x’ + b
    return a * x_test + b  

interpolate([1, 3, 5], [1, 9, 25], 5.0) 我预计产量是25,但实际产量是17.0。

Python社区是高质量的Python/Django开发社区
本文地址:http://www.python88.com/topic/56378
 
1352 次点击  
文章 [ 1 ]  |  最新文章 5 年前
R4444
Reply   •   1 楼
R4444    6 年前

我想你在找这样的东西:

def interpolate(x, y, x_test):
    for i in range(len(x)):
        if x[i] > x_test:   #extrapolated condition: when the largest value of
            x_below = i - 1 #list x is greater than x_test
            x_above = i
            y_below = i - 1
            y_above = i
            continue # <---- I changed break to continue
        elif x[i] < x_test: #extrapolated condition: when the largest value of
            x_below = i + 1 #list x is greater than x_test
            x_above = i
            y_below = i + 1
            y_above = i
            continue # <---- I changed break to continue
        else:             #interpolated condition: when x_test lies between
            return y[i]    #two sample points.

    #a = (yabove - ybelow) / (xabove - xbelow)
    a = (y[y_above] - y[y_below]) / (x[x_above] - x[x_below])
    #b = ybelow - a * xbelow
    b = y[y_below] - a * x[x_below]
    #y’ = a * x’ + b
    return (a * x_test + b)

print(interpolate([1, 3, 5], [1, 9, 25], 5.0))

输出:

25

breaks continues