1: ..."/>
溫馨提示×

溫馨提示×

您好,登錄后才能下訂單哦!

密碼登錄×
登錄注冊×
其他方式登錄
點(diǎn)擊 登錄注冊 即表示同意《億速云用戶服務(wù)條款》

Python實(shí)現(xiàn)牛頓插值法(差商表)

發(fā)布時間:2020-04-29 04:27:38 來源:網(wǎng)絡(luò) 閱讀:2822 作者:even_07 欄目:編程語言

def func(x,y,X,infor=True):
    list2=[y[0]]       #  差商表的對角線的第一個元素始終是y0
    count=1
    while(True):
        if len(y)>1:
            list=[]                 # 空列表用來保存,每次計(jì)算后差商表的行
            for i in range(len(y)-1):
                n=x[i+count]-x[i]
                m=y[i+1]-y[i]
                l=m/n
                list.append(l)
            list2.append(list[0])        # list2用來記錄差商表的對角線元素,每計(jì)算一次,取行的第一個元素
            count += 1
            y = list
        else:
            break
    if infor:                           # 判斷是否要繼續(xù)計(jì)算,結(jié)果
        W=0
        for i in range(len(list2)):
            if i==0:
                w=list2[i]
            else:
                w = list2[i]
                for j in range(i):
                    w*=(X-x[j])
            W+=w
        print('牛頓插值:', W)
    return '牛頓差商表對角線列:%s' %list2

ret=func([0.32, 0.34, 0.36],[0.314567, 0.333487, 0.352274],'',infor=False)
print(ret)
ret=func([0.32, 0.34, 0.36],[0.314567, 0.333487, 0.352274],0.3367)
print(ret)

運(yùn)行結(jié)果:
Python實(shí)現(xiàn)牛頓插值法(差商表)

向AI問一下細(xì)節(jié)

免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點(diǎn)不代表本網(wǎng)站立場,如果涉及侵權(quán)請聯(lián)系站長郵箱:is@yisu.com進(jìn)行舉報(bào),并提供相關(guān)證據(jù),一經(jīng)查實(shí),將立刻刪除涉嫌侵權(quán)內(nèi)容。

AI