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)文章:
主站蜘蛛池模板: 午夜精品久久久久久久星辰影院 | 亚洲日本三级 | 日韩三区在线观看 | 成人黄色在线 | 久久久久久免费毛片精品 | 欧美 日韩 国产 一区 | 91在线精品视频 | 日本淫视频 | 亚洲人成人一区二区在线观看 | 白浆在线 | 99综合网 | 国产天堂 | 亚洲视频二区 | 久草中文在线 | 一级片视频免费观看 | 欧美一区2区三区4区公司二百 | 国产小视频在线 | 国产精品无码久久久久 | 成人免费视屏 | 国产乱码精品一品二品 | 免费av观看 | 91精品久久久久久久久久入口 | 欧美日韩亚洲系列 | 成人免费看片 | 日韩国产精品一区二区三区 | 伊人无码高清 | 久久国产精品色av免费观看 | 久久久久久高清 | 7777精品伊人久久精品影视 | 久久久久国产一区二区三区四区 | 国产精品久久久久久久久久 | 一区二区三区四区电影视频在线观看 | 精品一区视频 | 中文字幕亚洲区一区二 | 亚洲午夜电影 | 91爱啪啪 | 日韩免费毛片视频 | 免费午夜视频在线观看 | 国产精品高清一区二区三区 | 久操伊人 | 爱爱无遮挡 |