Python高等工程数学

Posted 2021-07-31 yangbocsu

牛顿插商表（Python）

#yangbocsu 2021.06.2 民主楼
 
import numpy as np
 
def Newton(x,y):     
    diagonal_Line = [y[0]]  #差商表对角线第一个元素
    gap, n, d=1, len(x),[]
    D = np.zeros((n,n+1))  #矩阵大小 n x n+1
 
    for k in range(n):#装入第一列和第二列的值即x y 
        D[k][0] = x[k]
        D[k][1] = y[k]
    
    #计算插商值
    while(True):
        if len(y)>1:  
            l=[]                 # 空列表用来保存,每次计算后差商表的行
            for i in range(len(y)-1):
                xn=x[i+gap]-x[i]
                yn=y[i+1]-y[i]
                l.append(yn/xn)
                d.append(yn/xn)
            diagonal_Line.append(l[0])
            gap += 1
            y = l # 每列的y值
        else:
            break
    #牛顿插商表
    h, n, index1=0, len(x),1
    for k in range(2,n+1):#列  2,3,...,n
        index = index1
        while True:
            D[index][k] = d[h] #diagonal_L