Python实现牛顿插值法(差商表)-创新互联
def func(x,y,X,infor=True):
list2=[y[0]] # 差商表的对角线的第一个元素始终是y0
count=1
while(True):
if len(y)>1:
list=[] # 空列表用来保存,每次计算后差商表的行
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用来记录差商表的对角线元素,每计算一次,取行的第一个元素
count += 1
y = list
else:
break
if infor: # 判断是否要继续计算,结果
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)
运行结果:
另外有需要云服务器可以了解下创新互联scvps.cn,海内外云服务器15元起步,三天无理由+7*72小时售后在线,公司持有idc许可证,提供“云服务器、裸金属服务器、高防服务器、香港服务器、美国服务器、虚拟主机、免备案服务器”等云主机租用服务以及企业上云的综合解决方案,具有“安全稳定、简单易用、服务可用性高、性价比高”等特点与优势,专为企业上云打造定制,能够满足用户丰富、多元化的应用场景需求。
分享题目:Python实现牛顿插值法(差商表)-创新互联
文章出自:http://azwzsj.com/article/eesgd.html