def func(x,y,X,infor=True):
list2=[y[0]] # 差商表的對(duì)角線的第一個(gè)元素始終是y0
count=1
while(True):
if len(y)>1:
list=[] # 空列表用來(lái)保存,每次計(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用來(lái)記錄差商表的對(duì)角線元素,每計(jì)算一次,取行的第一個(gè)元素
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 '牛頓差商表對(duì)角線列:%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é)果:
十余年創(chuàng)新互聯(lián)網(wǎng)站建設(shè),由一走到現(xiàn)在,當(dāng)中離不開團(tuán)隊(duì)頑強(qiáng)的創(chuàng)業(yè)精神,離不開伴隨我們同行的客戶與專業(yè)的合作伙伴,創(chuàng)力信息一直秉承以“見一個(gè)客戶,了解一個(gè)行業(yè),交一個(gè)朋友”的方式為經(jīng)營(yíng)理念,提出“讓每一個(gè)客戶成為我們的終身客戶”為目標(biāo),以為用戶提供精細(xì)化服務(wù),全面滿足用戶需求為宗旨,誠(chéng)信經(jīng)營(yíng),更大限度為用戶創(chuàng)造價(jià)值。期待邁向下一個(gè)更好的十余年。
當(dāng)前題目:Python實(shí)現(xiàn)牛頓插值法(差商表)
本文地址:http://www.2m8n56k.cn/article46/jcgceg.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供品牌網(wǎng)站設(shè)計(jì)、軟件開發(fā)、品牌網(wǎng)站制作、企業(yè)網(wǎng)站制作、網(wǎng)站排名、面包屑導(dǎo)航
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請(qǐng)盡快告知,我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng),如需處理請(qǐng)聯(lián)系客服。電話:028-86922220;郵箱:[email protected]。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時(shí)需注明來(lái)源: 創(chuàng)新互聯(lián)