av一区二区在线观看_亚洲男人的天堂网站_日韩亚洲视频_在线成人免费_欧美日韩精品免费观看视频_久草视

您的位置:首頁技術(shù)文章
文章詳情頁

用Python實現(xiàn)Newton插值法

瀏覽:2日期:2022-06-22 10:19:10
1. n階差商實現(xiàn)

def diff(xi,yi,n): ''' param xi:插值節(jié)點xi param yi:插值節(jié)點yi param n: 求幾階差商 return: n階差商 ''' if len(xi) != len(yi): #xi和yi必須保證長度一致return else:diff_quot = [[] for i in range(n)]for j in range(1,n+1): if j == 1:for i in range(n+1-j): diff_quot[j-1].append((yi[i]-yi[i+1]) / (xi[i] - xi[i + 1])) else:for i in range(n+1-j): diff_quot[j-1].append((diff_quot[j-2][i]-diff_quot[j-2][i+1]) / (xi[i] - xi[i + j])) return diff_quot

測試一下:

xi = [1.615,1.634,1.702,1.828]yi = [2.41450,2.46259,2.65271,3.03035]n = 3print(diff(xi,yi,n))

返回的差商結(jié)果為:

[[2.53105263157897, 2.7958823529411716, 2.997142857142854], [3.0440197857724347, 1.0374252793901158], [-9.420631485362996]]

2. 牛頓插值實現(xiàn)

def Newton(x): f = yi[0] v = [] r = 1 for i in range(n):r *= (x - xi[i])v.append(r)f += diff_quot[i][0] * v[i] return f

測試一下:

x = 1.682print(Newton(x))

結(jié)果為:

2.5944760289639732

3.完整Python代碼

def Newton(xi,yi,n,x): ''' param xi:插值節(jié)點xi param yi:插值節(jié)點yi param n: 求幾階差商 param x: 代求近似值 return: n階差商 ''' if len(xi) != len(yi): #xi和yi必須保證長度一致return else:diff_quot = [[] for i in range(n)]for j in range(1,n+1): if j == 1:for i in range(n+1-j): diff_quot[j-1].append((yi[i]-yi[i+1]) / (xi[i] - xi[i + 1])) else:for i in range(n+1-j): diff_quot[j-1].append((diff_quot[j-2][i]-diff_quot[j-2][i+1]) / (xi[i] - xi[i + j])) print(diff_quot)f = yi[0] v = [] r = 1 for i in range(n):r *= (x - xi[i])v.append(r)f += diff_quot[i][0] * v[i] return f

到此這篇關(guān)于用Python實現(xiàn)牛頓插值法的文章就介紹到這了,更多相關(guān)python牛頓插值法內(nèi)容請搜索好吧啦網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持好吧啦網(wǎng)!

標(biāo)簽: Python 編程
相關(guān)文章:
主站蜘蛛池模板: 久久久涩 | 国产视频1 | 日韩国产专区 | 黄网站在线观看 | 欧美日韩精品中文字幕 | 九九九国产 | 在线免费观看视频你懂的 | 日日草夜夜草 | 国产在线a | 超碰免费观看 | 婷婷中文字幕 | 精品国产一区二区三区日日嗨 | 成人午夜av | 一区二区三区视频 | 一级做a毛片 | 成人免费视频久久 | 亚洲免费视频一区二区 | 综合九九 | 免费视频久久 | 中文字幕一区二区三区不卡 | a亚洲精品 | 成人性视频免费网站 | 国产精品一区二区三区99 | 午夜视频在线免费观看 | 国产精品99久久久久久宅男 | 免费av观看| 91极品尤物在线播放国产 | 高清国产午夜精品久久久久久 | 久草资源 | 国产精品福利在线观看 | 亚洲天堂精品一区 | 伊人看片 | 天久久 | 亚洲a视频| 欧美日在线 | 日本又色又爽又黄的大片 | 日韩国产精品一区二区三区 | 在线免费观看黄a | 天天澡天天狠天天天做 | 黄色三级免费网站 | 一区二区国产在线观看 |